i1 : R = QQ[x, y, z] / ideal(x * y - z^2) o1 = R o1 : QuotientRing
i2 : D = zeroDivisor( R ) o2 = 0, the zero divisor of R o2 : WDiv