|  |  D.4.22.12 radical Procedure from libraryprimdec.lib(see  primdec_lib).
 
Example:Usage:
radical(I[, l]); I ideal, l list (optional)
Optional parameters in list l (can be entered in any order):
 0, "fullRad" -> full radical is computed (default)
 1, "equiRad" -> equiRadical is computed
 "KL" -> Krick/Logar algorithm is used
 "SL" -> modifications by Laplagne are used (default)
 "facstd" -> uses facstd to first decompose the ideal (default for non homogeneous ideals)
 "noFacstd" -> does not use facstd (default for homogeneous ideals)
 
Return:
ideal, the radical of I (or the equiradical if required in the input parameters)
Note:
A combination of the algorithms of Krick/Logar (with modifications by Laplagne) and Kemper is used.
Works also in positive characteristic (Kempers algorithm).
 |  | 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+z3+2z2+2
 | 
 
 |