From d3b26617d5d2bb1f64919dd465af73fb76e4f1a1 Mon Sep 17 00:00:00 2001
From: "David A. Madore" <david@procyon.(none)>
Date: Wed, 8 Jun 2011 18:33:23 +0200
Subject: =?UTF-8?q?[divers/sageries]=20Calcul=20de=20la=20r=C3=A9solvante?=
 =?UTF-8?q?=20pour=20C4=20dans=20D4,=20pour=20m=C3=A9moire.?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 divers/sageries/resolvante-c4-dans-d4 | 15 +++++++++++++++
 1 file changed, 15 insertions(+)
 create mode 100644 divers/sageries/resolvante-c4-dans-d4

(limited to 'divers')

diff --git a/divers/sageries/resolvante-c4-dans-d4 b/divers/sageries/resolvante-c4-dans-d4
new file mode 100644
index 0000000..5e7cfa4
--- /dev/null
+++ b/divers/sageries/resolvante-c4-dans-d4
@@ -0,0 +1,15 @@
+R.<z1,z2,z3,z4,p,a1,a2,a3,a4> = PolynomialRing(QQ,9,order='lex(4),lex(1),deglex(4)')
+asym1 = -(z1+z2+z3+z4)
+asym2 = +(z1*z2+z1*z3+z1*z4+z2*z3+z2*z4+z3*z4)
+asym3 = -(z1*z2*z3+z1*z2*z4+z1*z3*z4+z2*z3*z4)
+asym4 = +(z1*z2*z3*z4)
+Isym = R.ideal([a1-asym1, a2-asym2, a3-asym3, a4-asym4])
+pexpr = z1*z3+z2*z4
+Ip = R.ideal([p-pexpr])
+fexpr = z1*z2^2+z2*z3^2+z3*z4^2+z4*z1^2
+fexprConj = fexpr.subs({z1:z3,z3:z1})
+coef0 = fexpr*fexprConj
+coef1 = -(fexpr+fexprConj)
+B = (Isym+Ip).groebner_basis()
+coef1.reduce(B)
+coef0.reduce(B)
-- 
cgit v1.2.1