This function should probably be in the Core of Macaulay2.
i1 : S = ZZ/101[a..d] o1 = S o1 : PolynomialRing |
i2 : I = ideal(a^2, a*b^3, c*d)
2 3
o2 = ideal (a , a*b , c*d)
o2 : Ideal of S
|
i3 : f = random(3, I)
3 2 2 2 2 2
o3 = 24a - 36a b - 30a c - 29a d + 19a*c*d + 19b*c*d - 10c d - 29c*d
o3 : S
|
i4 : f % I == 0 -- so f is in the ideal I o4 = true |
i5 : degree f == {3}
o5 = true
|