| 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
|