FrbPrimaryDecomposition(I: IDEAL): LIST
/**/ Use R ::= QQ[x,y,z]; /**/ I := ideal(x^2, x*y^2, z^2); /**/ FrbPrimaryDecomposition(I); [ideal(x^2, y^2, z^2), ideal(x, z^2)]