A ghost perturbation scheme to solve ordinary differential equations

arXiv2206.02445MaRDI QIDQ6401229FDOQ6401229

Authors: P. L. Garrido 

Publication date: 6 June 2022

Abstract: We propose an algebraic method that finds a sequence of functions that exponentially approach the solution of any second-order ordinary differential equation (ODE) with any boundary conditions. We define an extended ODE (eODE) composed of a linear generic differential operator that depends on free parameters, $p$, plus an $epsilon$ perturbation formed by the original ODE minus the same linear term. After the eODE's formal $epsilon$ expansion of the solution, we can solve order by order a hierarchy of linear ODEs, and we get a sequence of functions $y_n(x;epsilon,p)$ where $n$ indicates the number of terms that we keep in the $epsilon$-expansion. We fix the parameters to the optimal values $p^{\ast }(n)$ by minimizing a distance function of $y_n$ to the ODE's solution, $y$, over a given $x$-interval. We see that the eODE's perturbative solution converges exponentially fast in $n$ to the ODE solution when $epsilon=1$: $vert{y}_n(x;epsilon=1,p^{\ast }(n))-y(x)vert<Cdelt{a}^{n+1}$ with $delta<1$. The method permits knowing the number of solutions for Boundary Value Problems just by looking at the number of minima of the distance function at each order in $n$, $p^{\ast ,alpha}(n)$, where each $alpha$ defines a sequence of functions $y_n$ that converges to one of the ODE's solutions. We present the method by its application to several cases where we discuss its properties, benefits and shortcomings, and some practical algorithmic improvements.

This page was built for publication: A ghost perturbation scheme to solve ordinary differential equations

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6401229)