inout.lib
LIB "inout.lib"; ring r=0,(x,y,z),ls; ideal i= x,z+3y,x+y,z; matrix m[3][3]=i^2; pmat(m); ==> x2, xz+3xy, xy+x2, ==> xz, z2+6yz+9y2, yz+3y2+xz+3xy, ==> z2+3yz, y2+2xy+x2, yz+xz pmat(m,5); ==> x2, xz+.., xy+x2, ==> xz, z2+.., yz+.., ==> z2+.., y2+.., yz+xz