A simple method for the numerical solution of differential equations. (Q566924)

Zur genäherten Integration eine Gleichung \(y'=f(x,y)\) wird die Formel \[ y_{n+1}=y_{n-1}+ \frac 13 h (y'_{n-1}+4 y'_n+y'_{n+1}) \] angewendet. Mit geschätztem \(y'_{n+1}\) wird \(y_{n+1}\) nach dieser Formel berechnet und damit ein verbessert Wert von \(y'_{n+1}\) aus der Gleichung gewonnen. Dieser Prozeßwird fortgesetzt, bis innerhalb der verlangten Genauigkeit keine Änderung mehr eintritt. Vor Beginn der Rechnung wird mittels der Formel für den Fehler der \textit{Simpson}schen Regel die Intervallbreite abgeschätzt, und durch Reihenentwicklung oder nach dem \textit{Piccard}schen Näherungsverfahren werden die Ausgangswerte berechnet. Das Verfahren läßt sich auf Gleichungen höherer Ordnung übertragen, wie am Beispiel der \textit{Bessel}schen Differentialgleichung gezeigt wird. Am Schlußwird kurz gezeigt, wie sich ein Rechenfehler geltend macht.