ring r=0,(x,y,z),ds;
ideal i=x^4-y*z^2,x*y-z^3,y^2-x^3*z;  // the space curve singularity
qhweight(i);

ring r1 = 32003,(x,y,z),ds;
int a,b,c,t=11,5,3,0;
poly f = x^a+y^b+z^(3*c)+x^(c+2)*y^(c-1)+x^(c-1)*y^(c-1)*z3+
       x^(c-2)*y^c*(y^2+t*x)^2;
option(noprot);
timer=1;
ring r2 = 32003,(x,y,z),dp;
poly f=imap(r1,f);
ideal j=jacob(f);
vdim(std(j));

ring s = 0, (x1,x2,x3,x4,x5), wp(42,54,58,78,231);
ideal i = imap(r,i);
i;

ring r=0,(x,y,z),ds;
ideal j=(x-z,x^3-y^2);
ideal i=(x-y,x^3-z^2);
ideal k=intersect(i,j);

minbase(k);