IsPrimary(I: IDEAL): BOOL
/**/ Use P ::= QQ[x,y,z]; /**/ IsPrimary(ideal(x^2, y^2, z)); true /**/ IsPrimary(ideal(x*(x-1), y^2, z)); false