-- -*- M2-comint -*- {* hash: -539445018 *} i1 : R = ZZ/32003[a..d]; i2 : I = monomialCurveIdeal(R,{1,3,4}) 3 2 2 2 3 2 o2 = ideal (b*c - a*d, c - b*d , a*c - b d, b - a c) o2 : Ideal of R i3 : J = ideal(a^3,b^3,c^3-d^3) 3 3 3 3 o3 = ideal (a , b , c - d ) o3 : Ideal of R i4 : I = intersect(I,J) 4 3 3 3 4 3 3 4 6 3 2 o4 = ideal (b - a d, a*b - a c, b*c - a*c d - b*c*d + a*d , c - b*c d - ------------------------------------------------------------------------ 3 3 5 5 2 3 2 3 2 4 2 4 3 3 3 3 2 3 c d + b*d , a*c - b c d - a*c d + b d , a c - a d + b d - a c*d , ------------------------------------------------------------------------ 3 3 3 3 2 3 3 2 3 2 2 3 2 3 3 2 3 2 b c - a d , a*b c - a c*d + b c*d - a*b d , a b*c - a c d + b c d - ------------------------------------------------------------------------ 2 3 3 3 3 2 4 2 3 2 a b*d , a c - a b*d , a c - a b d) o4 : Ideal of R i5 : removeLowestDimension I 3 2 2 2 3 2 o5 = ideal (b*c - a*d, c - b*d , a*c - b d, b - a c) o5 : Ideal of R i6 : topComponents I 3 2 2 2 3 2 o6 = ideal (b*c - a*d, c - b*d , a*c - b d, b - a c) o6 : Ideal of R i7 : radical I 2 2 3 2 6 3 3 2 4 5 o7 = ideal (b*c - a*d, a*c - b d, b - a c, c - c d - b d + b*d ) o7 : Ideal of R i8 : minimalPrimes I 3 2 2 2 3 2 o8 = {ideal (- b*c + a*d, - c + b*d , a*c - b d, - b + a c), ideal (- c + ------------------------------------------------------------------------ 2 2 d, b, a), ideal (c + c*d + d , b, a)} o8 : List i9 :