Procedure from library primdec.lib (see primdec_lib).
primdec.lib
radical(i); i ideal.
ideal, the radical of i.
A combination of the algorithms of Krick/Logar and Kemper is used. Works also in positive characteristic (Kempers algorithm).
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; ideal pr= radical(i); pr; → pr[1]=z2-y → pr[2]=y2z+yz+2y+2