Esempio:

> with(linalg):

> T3d:=(x0,y0,z0)->matrix(4,4,[[1,0,0,0],[0,1,0,0],[0,0,1,0],[x0,y0,z0,1]]):

> trasforma:=(f,m)->delcols(evalm(augment(f,vector(rowdim(f),1))&*m),rowdim(m)..rowdim(m)):

> trasforma([[x,y,z]],T3d(h,k,l));

>