//SINGULAR ExampleA.4.7




ring R = 0,(x,y,z),dp;
poly f = x2z+y3;
ideal I = maxideal(1)*f;    //<xf,yf,zf>
LIB"primdec.lib";
radical(I);

LIB"surf.lib";
plot(f);