-- -*- M2-comint -*- hash: 4202898279937703506

i1 : setRandomSeed 0
 -- setting random seed to 0

o1 = 0

i2 : kk=ZZ/101; d=1;

i4 : n=2*d

o4 = 2

i5 : R=kk[a_0..a_(binomial(n+2,2)-1)]

o5 = R

o5 : PolynomialRing

i6 : S=kk[x_0..x_(n-1),a_0..a_(binomial(n+2,2)-1)]

o6 = S

o6 : PolynomialRing

i7 : M=genericSymmetricMatrix(S,a_0,n)

o7 = | a_0 a_1 |
     | a_1 a_2 |

             2      2
o7 : Matrix S  <-- S

i8 : X=(vars S)_{0..n-1}

o8 = | x_0 x_1 |

             1      2
o8 : Matrix S  <-- S

i9 : Y=X*M

o9 = | x_0a_0+x_1a_1 x_0a_1+x_1a_2 |

             1      2
o9 : Matrix S  <-- S

i10 : (M1,M2)=matrixFactorizationK(X,Y);

i11 : (eOdd,eEv)=cliffordOperators(M1,M2,R);

i12 : sM = symMatrix(eOdd,eEv)

o12 = | a_0 a_1 |
      | a_1 a_2 |

              2      2
o12 : Matrix R  <-- R

i13 : f = det sM

         2
o13 = - a  + a a
         1    0 2

o13 : R

i14 : f == (cliffordModule(eOdd,eEv)).hyperellipticBranchEquation

o14 = true

i15 : (c0,c1)=centers(eOdd,eEv)

o15 = ({-3} | 0  -a_1^2+a_0a_2 |, {-1} | -a_1 -a_2 |)
       {-2} | -1 0             |  {-1} | a_0  a_1  |

o15 : Sequence

i16 : assert(c0^2-(-1)^d*f*id_(source c0)==0)

i17 : assert(c1^2-(-1)^d*f*id_(source c1)==0)

i18 : 
