Moin.
Ich habe mir jetzt die halbe Nacht um die Ohren geschlagen, und eine funktionstüchtige numerische Approximation auf die Beine gestellt.
Der Code kann unter anderem:
- symbolische Differentiation von Produkten nach Leibniz-Regel
- numerische Differentiation bis zur 4. Ordnung nach 5-point stencil
- numerische Integration nach adaptiver Simpson-Regel
Das Ergebnis entspräche einer Approximation der Ordnung 4 meines nicht funktionstüchtigen numerischen Ansatz.
Löblicherweise stimmen die Ergebnisse meiner funktionierenden Skripte im Rahmen des zu erwartenden Fehlers miteinander überein.
Schöne Arbeit! Glückwunsch!
Drr Simpson war ein Geniestreich von Dir; ich hab grade mal g[0] durch die Näherungsformel für s(phi) aus dem Skript (l, m, n...; Sinusterme) und das nähert den Hochwert um ca. 0.1% SCHLECHTER als der Simpson.
Nochmals: Glückwunsch!
Grüsse
Solkar