Truncation Bounds for Differentially Finite Series
Bornes pour la troncature des séries différentiellement finies
Résumé
We describe a flexible symbolic-numeric algorithm for computing bounds on the tails of series solutions of linear differential equations with polynomial coefficients. Such bounds are useful in rigorous numerics, in particular in rigorous versions of the Taylor method of numerical integration of ODEs and related algorithms. The focus of this work is on obtaining tight bounds in practice at an acceptable computational cost, even for equations of high order with coefficients of large degree. Our algorithm fully covers the case of generalized series expansions at regular singular points. We provide a complete implementation in SageMath and use it to validate the method in practice.
We describe a flexible symbolic-numeric algorithm for computing bounds on the tails of series solutions of linear differential equations with polynomial coefficients. Such bounds are useful in rigorous numerics, in particular in rigorous versions of the Taylor method of numerical integration of ODEs and related algorithms. The focus of this work is on obtaining tight bounds in practice at an acceptable computational cost, even for equations of high order with coefficients of large degree. Our algorithm fully covers the case of generalized series expansions at regular singular points. We provide a complete implementation in SageMath and use it to validate the method in practice. Résumé.-Nous décrivons un algorithme symbolique-numérique souple pour le calcul de bornes sur les restes de solutions séries d'équations différentielles linéaires à coefficients polyno-miaux. Ces bornes sont destinées au calcul numérique rigoureux, et utiles notamment dans des versions rigoureuses de la méthode de Taylor d'intégration des EDO ou d'autres algorithmes apparentés. L'objectif principal de ce travail est d'obtenir des bornes fines en pratique pour un coût de calcul acceptable, y compris dans le cas d'équations d'ordre élevé à coefficients de grand degré. Notre algorithme couvre entièrement le cas des développement en séries généralisées au voisinage de points singuliers réguliers. Nous présentons une implémentation complète de la méthode en SageMath, et nous l'utilisons pour valider son bon comportement en pratique.
Origine : Accord explicite pour ce dépôt
Loading...
