A fifth-order interpolant for the Dormand and Prince Runge-Kutta method (Q1263257): Difference between revisions

On considère le problème avec la condition initial \[ y'(t)=f(t,y(t)),\quad t\geq t_ 0,\quad y\in {\mathbb{R}}^ N,\quad y(t_ 0)=y_ 0. \] Dans ce travail les AA. présentent une famille d'interpolants de l'ordre 5 pour la solution d'ordre 5, offerte par la paire Runge-Kutta RK5(4)7M de \textit{J. R. Dormand} et \textit{P. J. Prince} [ibid. 6, 19-26 (1980; Zbl 0448.65045)]. L'application dans les calculs de cette famille d'interpolants nécessite deux évaluations additionnelles, à chaque pas, relativement à la fonction f. On apprécie que la paire RK5(4)7M ayant les ordres 4 et 5, qui a été établie par Dormand et Prince, paraît être supérieure à la paire classique de Fehlberg. Les AA. montrent que, si l'on ajoute une étape à DOPRI5(4), on peut obtenir une approximation de l'ordre 5 dans tout point de l'intervalle \([t_ n,t_{n+1}]\). Puisque dans ce processus sont 4 degrés de liberté, on obtient une famille d'interpolants RK pour DOPRI5(4), qui dépend de 4 paramètres. Puis on détermine une méthode (un interpolant) optimale dans cette famille par la choix des paramètres libres pour minimiser, dans une certaine norme, les coefficients dominants de l'erreur de troncation locale de la solution continue. En particulier, avec cette mésure d'erreur, l'interpolant proposé dans ce travail est plus exact par rapport à un interpolant antérieur (DPS), proposé par \textit{L. F. Shampine} [Math. Comput. 173, 135-150 (1986; Zbl 0594.65046)]. Quelques-uns résultats numériques, présentés, dans la section 3 du travail, montrent l'exactité de la solution continue dans tous les points de l'intervalle d'intégration.