Numerical integration of ordinary differential equations. (Q1451369)

Das geschilderte Verfahren zur numerischen Integration einer Differentialgleichung \(y'= f(x, y)\) entsteht durch geringe Abänderung aus dem von \textit{J. C. Adams} angegebenen Verfahren. Ist die Integration so weit vorgeschritten, daß die zu den äquidistanten Abszissen \(\ldots, x_{n-2},x_{n-1}, x_n\) zugehörigen Ordinaten \(\ldots, y_{n-2}, y_{n-1}, y_n\) bereits bekannt sind, so wird durch die Punkte \(x_\nu, f_\nu\equiv f(x_\nu, y_\nu)\) einmal für \(\nu = n\), \(n-1, n - 2,n- 3\), und einmal für \(\nu= n+1, n, n- 1, n- 2\) die Interpolationsparabel dritten Grades gelegt und über das Intervall \(\langle x_{n-3},x_{n+1} \rangle\) bzw. \(\langle x_{n-1},x_{n+1} \rangle\) integriert. Es entstehen die Formeln \[ y_{n+1} = y_{n-3} + \frac 43 h(2f_n-f_{n-1} + 2f_{n-2} + 0\cdot f_{n-3}) \] und \[ y_{n+1} = y_{n-1} + \frac 13 h(f_{n+1}+4f_n + f_{n-1} + 0\cdot f_{n-2}), \] von denen die erste \(y_{n+1}\) zu berechnen und die zweite \(y_{n +1}\) durch Iteration zu verbessern gestattet. Die auftretenden Zahlenkoeffizienten, von denen zwei verschwinden, sind bequem für die praktische Rechnung. Es folgen einige Bemerkungen über die Berechnung der Ausgangswerte, über die zu wählende Stellenzahl, den begangenen Fehler, sowie ein Beispiel.