Hallo
Die 4 Nullstellen eines normierten Polynoms 4. Grades x4 +ax3 +bx2 +cx +d = 0 lassen sich wie folgt finden: man ersetzt zuerst 4x +a = z, also
(4x)4 +4a (4x)3 +16b (4x)2 +64c (4x) +256d = 0
(4x +a)4 +(16b -6a2) (4x +a)2 +(64c -4a3 -32ab +12a3) (4x +a) +256d -a4 -16ba2 +6a4 -64ac +4a4 +32a2b -12a4 = 0
z4 +(16b -6a2) z2 +(8a3 -32ab +64c) z -3a4 +16a2b --64ac +256d = 0
z4 +k z2 +m z +n = 0
Dieses Polynom in z, das keinen Koeffizienten 3. Grades mehr besitzt, kann durch folgenden Ansatz von Tartaglia geloest werden:
(z2 +P)2 = (Q z +R)2 => z2 +P = +/- (Q z +R)
==> z1,2 = +Q/2 +/- (Q2/4 -P +R)1/2
∨ z3,4 = -Q/2 +/- (Q2/4 -P -R)1/2
z4 +(2P -Q2) z2 -2QR z +P2 -R2 = 0
Der Koeffizienten-Vergleich ergibt:
(1) 2 P -Q2 = k ==> Q2 = 2 P -k
(2) P2 -R2 = n ==> R2 = P2 -n
(3) -2 Q R = m
==> 4 Q2 R2 = m2 ==> 4 (2 P -k) * (P2 -n) = m2 => 4 (2 P3 -k P2 -2n P +kn) = m2
P ist also eine Nullstelle eines normierten Polynoms 3. Grades: P3 -k/2 P2 -n P +k*n/2 -m2/8 = 0
==> (2P)3 -k (2P)2 -4n (2P) +4 kn -m2 = 0
==> (6P)3 -3k (6P)2 -36n (6P) +108 kn -27 m2 = 0
==> (6P -k)3 +(-36n -3k2) (6P -k) +108kn -27m2 +k3 -36nk -3k3 = 0
=> (6P -k)3 +(-36n -3k2) (6P -k) -2k3 +72kn -27m2 = 0
Es gilt also folgenden Fälle zu unterscheiden
(a) -3k2 -36n = 0
==> 6P -k = (+2k3 -72kn +27m2) 1/3
(b) -3k2 -36n > 0
==> 6P -k = 2*(-k2 -12n)1/2*sinh 1/3 arsinh (2 k3 -72 kn +27 m2)/2/(-k2 -12 n)3/2
(c) -3k2 -36n < -1
==> 6P -k = 2*(+k2 +12n)1/2
* |2 k3 -72 kn +27 m2| / (2k3 -72 kn +27 m2)
*cosh 1/3 arcosh |2 k3 -72 kn +27 m2|/2/(-k2 -12 n)3/2
Falls der Zähler < 0 ist, ist auch 6 P -k < 0 und der arcosh wird vom Betrag des Quotienten gerechnet
(d) -1 < -3k2 -36n < 0
==> 6P -k = 2*(+k2 +12n)1/2*cos 1/3 arccos (2 k3 -72 kn +27 m2)/2/(-k2 -12 n)3/2
Die Lösungen der Gleichung 4. Grades sind
z1,2 = +Q/2 +/- (Q2/4 -P +R)1/2
∨ z3,4 = -Q/2 +/- (Q2/4 -P -R)1/2
dabei gilt
(1) 2 P -Q2 = k ==> Q2 = 2 P -k ==> Q = +/- (2 P -k)1/2 ==> Q/2 = +/- (P/2 -k/4)
(2) P2 -R2 = n ==> R2 = P2 -n ==> R = +/- (P2 -n)1/2
(3) - 2 Q R = m ==> Q = -m/R/2
Wird R wird mit positivem Vorzeichen gewaehlt, erhält Q das negierte Vorzeichen des Koeffizienten m erhält.
Wird R mit negativem Vorzeichen gewaehlt, erhält Q das nicht-negierte Vorzeichen des Koeffizienten m.
z4 +k z2 +m k +n = 0
z1 = -Q/2 + (Q2/4 -P +R)1/2
z2 = -Q/2 - (Q2/4 -P +R)1/2
z1 = +Q/2 + (Q2/4 -P -R)1/2
z2 = +Q/2 - (Q2/4 -P -R)1/2
Aus z erhält man x : 4x +a = z -> x = (z -a)/4
Beispiel:
(A) x4 +4 x3 +3 x2 +2 x +1 = 0
(4x)4 +4*4 (4x)3 +16*3 (4x)2 +64*2 (4x) +256 = 0
z4 +(16*3 -6*16) z2 +(8*64 -32*4*3 +64*2) z -3*256 +16*16*3 --64*4*2 +256*1 = 0
z4 -48 z2 +256 z -256 = 0
(6P)3 -3*(-48) (6P)2 -36*(-256) (6P) +108 *(-48)*(-256) -27 * 65536 = 0
(6P -(-48))3 +(-36*(-256) -3*2304) (6P -k) -2*(-110592) +72*(-48)*(-256) -27*65536 = 0
(6P +48)3 +2304 (6P +48) -663552 = 0
2304 > 0 ==> 6P1 +48 = 2 *(768)1/2 * sinh 1/3 arsinh +331776/(768)3/2 = 78,44953406
P = 5,074922343
z1,2 = -(P/2 -k/4)1/2 +/- ((P2 -n)1/2 -k/4 -P/2)1/2 = -(P/2 +12)1/2 +/- ( (P2 +256)1/2 +12 -P/2)1/2
= -3.812802273 +/- 5.123289298
∨ z3,4 = +(P/2 -k/4)1/2 +/- (-(P2 -n)1/2 -k/4 -P/2)1/2 = +(P/2 +12)1/2 +/- ( -(P2 +256)1/2 +12 -P/2)1/2 = +3.812802273 +/- i*2.706107088
==> z1 = +1.310487025 ∨ z2 = -8.936091571 ∨ z3,4 = +3.812802273 +/- i*2.706107088
==> x1 = 1/4(z1 -a) = 1/4(z1 -4) = -0.672378243 ∨ x2 = 1/4(z2 -a) = 1/4(z2 -4) = -3.234022893
∨ x3,4 = 1/4(z3,4 -a) = 1/4*z3,4 -1 = +3.812802273/4 -1 +/- i*2.706107088/4 = -0.046799431 +/- i*0.676526772
==> x1 = -0.672378243 ∨ x2 = -3.234022893
∨ x3,4 = -0.046799431 +/- i*0.676526772 = -0.678143539*ei*arccos 0.069011099 = -0.678143539*ei*pi/180*86.04281 = rr*eiphi
4phi = 4*86.04281° = 344.17124°
3phi = 3*86.04281° = 258.12843°
2phi = 2*86.04281° = 172.08562°
Probe: x14 +4 x13 +3 x12 +2 x1 +1 = 0 (w) ∨ x24 +4 x23 +3 x22 +2 x2 +1 = 0 (w)
∨ x34 +4 x33 +3 x32 +2 x3 +1 = (rr4*cos 4phi +4*rr3*cos 3phi +3*rr2 *cos 2phi +2*rr*cos 1phi +1)
+i(rr4*sin 4phi +4*rr3 *sin 3phi +3*rr2 *sin 2phi +2*rr*sin 1phi) = 0 +i*0 (w)
∨ x44 +4 x43 +3 x42 +2 x4 +1 = (rr4*cos (-4phi) +4*rr3*cos (-3phi) +3*rr2 *cos (-2phi) +2*rr*cos (-1phi) +1)
+i(rr4*sin (-4phi) +4*rr3 *sin (-3phi) +3*rr2 *sin (-2phi) +2*rr*sin (-1phi)) = 0 +i*0 (w)