numerDecom.lib
LIB "numerDecom.lib"; ring r=0,(x,y,z),dp; poly f1=z^2; poly f2=z*(x^2+y); ideal I=f1,f2; def D=defl(I,1); setring D; DI;