i1 : Q = ZZ/101[x,y,z] o1 = Q o1 : PolynomialRing
i2 : isPrimary ideal(y^6) o2 = true
i3 : isPrimary(ideal(y^6), ideal(y)) o3 = true
i4 : isPrimary ideal(x^4, y^7) o4 = true
i5 : isPrimary ideal(x*y, y^2) o5 = false