primdec.lib
LIB "primdec.lib"; ring r = 0,(x,y,z),dp; poly p = z2+1; poly q = z3+2; ideal i = p*q^2,y-z2; list pr = prepareAss(i); pr; ==> [1]: ==> _[1]=z2-y ==> _[2]=y2z+z3+2z2+2