i1 : R = QQ[x,y,z]/ideal(x^2-y*z) o1 = R o1 : QuotientRing
i2 : m = ideal(x,y,z) o2 = ideal (x, y, z) o2 : Ideal of R
i3 : isReflexive(m) o3 = false
i4 : I = ideal(x,y) o4 = ideal (x, y) o4 : Ideal of R
i5 : isReflexive(I) o5 = true