sing.lib
LIB "sing.lib"; ring r = 32003,(x,y,z),ds; ideal i = xy,xz,yz; matrix T = deform(i); print(T); ==> x,0,0, ==> 0,0,z, ==> 0,y,0 print(deform(x3+y5+z2)); ==> xy3,y3,xy2,y2,xy,y,x,1