From 3488be4c25d9a8bc59952e84e2fb1179798a84bf Mon Sep 17 00:00:00 2001 From: "David A. Madore" Date: Sun, 27 Nov 2011 08:21:49 +0100 Subject: sageries: une correction + je continue le calcul un peu plus loin. --- divers/sageries/resolvante-d5-dans-m20 | 10 ++++++++-- 1 file changed, 8 insertions(+), 2 deletions(-) (limited to 'divers/sageries') 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] -- cgit v1.2.1