|  | LIB "gmspoly.lib";
ring R=0,(x,y,z),dp;
poly f=x+y+z+x2y2z2;
def Rs=goodBasis(f);
setring(Rs);
b;
==> b[1]=1
==> b[2]=s2x-3sx2+x3
==> b[3]=5/2x
==> b[4]=10s2x2-25/2sx3+5/2x4
==> b[5]=-25/4sx+25/4x2
print(jet(A,0));
==> 0,0,0,-25/8,0,    
==> 0,0,0,0,    125/8,
==> 1,0,0,0,    0,    
==> 0,1,0,0,    0,    
==> 0,0,1,0,    0     
print(jet(A/var(1),0));
==> 1/2,0,0,  0,0, 
==> 0,  1,0,  0,0, 
==> 0,  0,3/2,0,0, 
==> 0,  0,0,  2,0, 
==> 0,  0,0,  0,5/2
 |