i1 : R = QQ[x,y,z] o1 = R o1 : PolynomialRing |
i2 : I = ideal (x*y-z, x*z-y^2) 2 o2 = ideal (x*y - z, - y + x*z) o2 : Ideal of R |
i3 : bpd = binomialPrimaryDecomposition I 2 2 o3 = {ideal (y - x*z, x*y - z, x - y), ideal (z, y)} o3 : List |
i4 : intersect bpd == I o4 = true |