Ich bin auf der Suche nach reelen Ergebnissen.
Also dann wollen wir mal.
Zuerst bringst du die Formel in die Form einer kubischen Gleichung. Glücklicherweise ergibt sich dabei direkt sie sogenannte reduzierte Form (ohne quadratischen Term)
R^3+pR+q=0
mit
p=-4\frac{C_4C_9(M-O-C_7)(M-2O)^3}{C_8C_{10}}
und
q=4\frac{C_4C_9(LC_7+2NO)(M-2O)^3}{C_8C_{10}}
Diese Gleichung lässt sich mit den Formel von Cardano lösen.
Dazu rechnest du zuerst die Diskriminante aus
D=\left(\frac{q}{2}\right)^2+\left(\frac{p}{3}\right)^3
Jetzt kommst du ohne Fallunterscheidung nicht mehr weiter.
Der erste und einfachste Fall ist, dass D>0. Dann gibt es genau eine relle Lösung, und zwar
R=\sqrt[3]{-\frac{q}{2}+\sqrt{D}}+\sqrt[3]{-\frac{q}{2}-\sqrt{D}}
Der zweite Fall ist, dass D=0. Dann gibt es eine einfache und eine doppelte Lösung
R_1=\sqrt[3]{-4q}
R_{2,3}=\sqrt[3]{\frac{q}{2}}
Wenn q=0 ist, dann fallen diese Lösungen zur dreifachen Lösung R=0 zusammen.
Der dritte und komplizierteste Fall ist, dass DR_1=\sqrt{-\frac{4}{3}p}\cos\left(\frac{1}{3}\arccos\left(-\frac{q}{2}\sqrt{-\frac{27}{p^3}}\right)\right)
R_2=-\sqrt{-\frac{4}{3}p}\cos\left(\frac{1}{3}\left(\arccos\left(-\frac{q}{2}\sqrt{-\frac{27}{p^3}}\right)+\pi\right)\right)
R_2=-\sqrt{-\frac{4}{3}p}\cos\left(\frac{1}{3}\left(\arccos\left(-\frac{q}{2}\sqrt{-\frac{27}{p^3}}\right)-\pi\right)\right)
Als Algorithmus sieht das ganze dann ungefähr so aus.
Berechne
s=4\frac{C_4C_9(M-2O)^3}{C_8C_{10}}
p=-s(M-O-C_7)
q=s(LC_7+2NO)
D=\left(\frac{q}{2}\right)^2+\left(\frac{p}{3}\right)^3
Falls D\>0
R=\sqrt[3]{-\frac{q}{2}+\sqrt{D}}+\sqrt[3]{-\frac{q}{2}-\sqrt{D}}
Falls D=0
Falls q=0
R=0
Sonst
R_1=\sqrt[3]{-4q}
R_{2,3}=\sqrt[3]{\frac{q}{2}}
Falls D
R\_1=\sqrt{-\frac{4}{3}p}\cos\left(\frac{1}{3}\arccos\left(-\frac{q}{2}\sqrt{-\frac{27}{p^3}}\right)\right)
R\_2=-\sqrt{-\frac{4}{3}p}\cos\left(\frac{1}{3}\left(\arccos\left(-\frac{q}{2}\sqrt{-\frac{27}{p^3}}\right)+\pi\right)\right)
R\_2=-\sqrt{-\frac{4}{3}p}\cos\left(\frac{1}{3}\left(\arccos\left(-\frac{q}{2}\sqrt{-\frac{27}{p^3}}\right)-\pi\right)\right)
Und jetzt viel Spaß beim Programmieren !
hendrik