Procedure from library poly.lib (see poly_lib).
poly.lib
rad_con(g,I); g polynomial, I ideal
1 (TRUE) (type int) if g is contained in the radical of I 0 (FALSE) (type int) otherwise
Example:
LIB "poly.lib"; ring R=0,(x,y,z),dp; ideal I=x2+y2,z2; poly f=x4+y4; rad_con(f,I); → 0 ideal J=x2+y2,z2,x4+y4; poly g=z; rad_con(g,I); → 1