diff options
author | David A. Madore <david+git@madore.org> | 2011-11-27 08:21:49 +0100 |
---|---|---|
committer | David A. Madore <david+git@madore.org> | 2011-11-27 08:21:49 +0100 |
commit | 3488be4c25d9a8bc59952e84e2fb1179798a84bf (patch) | |
tree | 471e79faa4029e3c4e02f5faee18867bec99750f /divers | |
parent | e7b595293338f69f9ee7b7fe7391a7306a5ca551 (diff) | |
download | galois-3488be4c25d9a8bc59952e84e2fb1179798a84bf.tar.gz galois-3488be4c25d9a8bc59952e84e2fb1179798a84bf.tar.bz2 galois-3488be4c25d9a8bc59952e84e2fb1179798a84bf.zip |
sageries: une correction + je continue le calcul un peu plus loin.
Diffstat (limited to 'divers')
-rw-r--r-- | divers/sageries/resolvante-d5-dans-m20 | 10 |
1 files changed, 8 insertions, 2 deletions
diff --git a/divers/sageries/resolvante-d5-dans-m20 b/divers/sageries/resolvante-d5-dans-m20 index 2277bea..be91ce2 100644 --- a/divers/sageries/resolvante-d5-dans-m20 +++ b/divers/sageries/resolvante-d5-dans-m20 @@ -61,7 +61,7 @@ Isym = R.ideal([a1-asym1, a2-asym2, a3-asym3, a4-asym4, a5-asym5]) B = Isym.groebner_basis() pexpr = z1^2*(z2*z5 + z3*z4) + z2^2*(z1*z3 + z4*z5) + z3^2*(z1*z5 + z2*z4) + z4^2*(z1*z2 + z3*z5) + z5^2*(z1*z4 + z2*z3) qexpr = z1*z2 + z2*z3 + z3*z4 + z4*z5 + z1*z5 -dexpr = (z1-z2)*(z1-z3)*(z1-z4)*(z2-z3)*(z2-z4)*(z3-z4) +dexpr = (z1-z2)*(z1-z3)*(z1-z4)*(z1-z5)*(z2-z3)*(z2-z4)*(z2-z5)*(z3-z4)*(z3-z5)*(z4-z5) def debase(poly): r = poly.reduce(B) return [r.coefficient({z1:0,z2:i2,z3:i3,z4:i4,z5:i5}) for i2 in range(2) for i3 in range(3) for i4 in range(4) for i5 in range(5)] @@ -79,4 +79,10 @@ results = [[None for j in range(12)] for i in range(12)] for i in range(12): for j in range(12): print (i,j) - results[i][j] = qsubmat[complements[i],complements[j]].determinant() + results[i][j] = (-1)^(i+j)*(qsubmat[complements[j],complements[i]].determinant()) + +qsubadj = matrix(R, 12, 12, results) +tmp = [sum([pdsubmat[1][k]*qsubadj[k][j] for k in range(12)]) for j in range(12)] +subfactor = -tmp[11]/1280 +qsubdetfact = R(qsubdet/subfactor) +donq = [R(pol/subfactor)/qsubdetfact for pol in tmp] |