Procedure from library primdec.lib (see primdec_lib).
primdec.lib
prepareAss(i); i ideal
list, the radicals of the maximal dimensional components of i.
Uses algorithm of Eisenbud/Huneke/Vasconcelos.
Example:
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+yz+2y+2