firstOrderDeformation -- compute the first order very flat deformations of the N2 variety X

i1 : loadPackage("RandomCanonicalCurves",Reload=>true);

i2 : kk=ZZ/101

o2 = kk

o2 : QuotientRing

i3 : g=7

o3 = 7

i4 : S=kk[x_0..x_(g-1)]

o4 = S

o4 : PolynomialRing

i5 : Ican=(random canonicalCurve)(g,S);

o5 : Ideal of S

i6 : betti(fX=res(Ican,LengthLimit=>3))

            0  1  2  3
o6 = total: 1 10 16 16
         0: 1  .  .  .
         1: . 10 16  .
         2: .  .  . 16

o6 : BettiTally

i7 : T = firstOrderDeformation fX;

             10       18
o7 : Matrix S   <--- S

i8 : rank source T == 3*g-3

o8 = true