LIB "stdmodule.lib";
// Example 1:
ring r=QQ,(x,y,z),Dp;
ideal A=x2-yz,xy+x,xz;
module G=[x2yz+x2z-2,xz+xyz2-x3z],[x3y-xy2z+x3-xyz,-1],[xy+x,x2z2-x2+yz];
MSPOLY(G,A);
==> _[1]=-x7y*gen(1)+3x5y2z*gen(1)-x5z3*gen(2)+x4yz3*gen(1)-3x3y3z2*gen(1)+x3\
   yz4*gen(2)+xy4z3*gen(1)-x7*gen(1)+3x5yz*gen(1)+x4z3*gen(1)-3x3y2z2*gen(1)\
   +xy3z3*gen(1)-x6*gen(2)+x5y*gen(1)+3x4yz*gen(2)+x4z2*gen(2)-2x3y2z*gen(1)\
   +x3z3*gen(2)-3x2y2z2*gen(2)-x2yz3*gen(2)+xy3z2*gen(1)+y3z3*gen(2)+x5*gen(\
   1)-2x3yz*gen(1)+xy2z2*gen(1)+x4*gen(2)-2x3z*gen(1)-2x2yz*gen(2)+x2z2*gen(\
   2)-2x2z2*gen(1)+2xyz2*gen(1)+y2z2*gen(2)
// Example 2:
ring r2=QQ,(x,y,z),dp;
ideal A=x2-y,xz+z,y2-z;
module G=[x2yz+x2z-2, xyz2+xz-x3z],[x3y+x3-xy2z-xyz,xz-xz-1],[xy+x,x2z2-x2+yz];
MSPOLY(G,A);
==> _[1]=-x3z*gen(2)+x2yz*gen(1)+xyz2*gen(2)+x2z*gen(1)+xz*gen(2)-2*gen(1)
==> _[2]=x3y*gen(1)-xy2z*gen(1)+x3*gen(1)-xyz*gen(1)-gen(2)