i1 : R = QQ[x,y,z]/ideal(x^2-y*z) o1 = R o1 : QuotientRing
i2 : isDomain(R) o2 = true
i3 : S = ZZ/5[x,y]/ideal(x^2*y^3) o3 = S o3 : QuotientRing
i4 : isDomain(S) o4 = false