A<-matrix(c(1:9)/10, ncol=1)
A_derivs<-list2derivs(list(A, A^0, A^2, A^3, A^4), deriv_order=4)
B_derivs<-transform(A, type="exp", par=0, deriv_order=4)
ind2joint (list(A_derivs,B_derivs),
list(trind_generator(1),trind_generator(1)),
list(trind_generator(1),trind_generator(1+1)), 4)
Run the code above in your browser using DataLab