Abstract:When the eigenvalues of the coefficient matrix for a linear scalar ordinary differential equation are of large magnitude, its solutions exhibit complicated behaviour, such as high-frequency oscillations, rapid growth or rapid decay. The cost of representing such solutions using standard techniques grows with the magnitudes of the eigenvalues. As a consequence, the running times of most solvers for ordinary differential equations also grow with these eigenvalues. However, a large class of scalar ordinary differential equations with slowly-varying coefficients admit slowly-varying phase functions that can be represented at a cost which is bounded independent of the magnitudes of the eigenvalues of the corresponding coefficient matrix. Here, we introduce a numerical algorithm for constructing slowly-varying phase functions which represent the solutions of a linear scalar ordinary differential equation. Our method's running time depends on the complexity of the equation's coefficients, but is bounded independent of the magnitudes of the equation's eigenvalues. Once the phase functions have been constructed, essentially any reasonable initial or boundary value problem for the scalar equation can be easily solved. We present the results of numerical experiments showing that, despite its greater generality, our algorithm is competitive with state-of-the-art methods for solving highly-oscillatory second order differential equations. We also compare our method with Magnus-type exponential integrators and find that our approach is orders of magnitude faster in the high-frequency regime.

Solving linear differential equations with mixed arguments

Solving nonlinear neutral delay integro-differential equations via general linear methods

Regularity, synthesis, rigidity and analytic classification for linear ordinary differential equations of second order

Least-squares Solutions of Linear Differential Equations

On the computation of accurate initial conditions for linear higher-index differential-algebraic equations and its application in initial value solvers

General Solution to a Nonlocal Linear Differential Equation of First-Order

A solver for linear scalar ordinary differential equations whose running time is bounded independent of frequency

Differential Operator Method of Finding A Particular Solution to An Ordinary Nonhomogeneous Linear Differential Equation with Constant Coefficients

The Solvability of the Discrete Boundary Value Problem on the Half-Line

Initial-boundary value problems for discrete evolution equations: discrete linear Schrodinger and integrable discrete nonlinear Schrodinger equations

Analytical solution of linear ordinary differential equations by differential transfer matrix method

Rough differential equations with path-dependent coefficients

Evolution Equations with Liouville Derivative on R without Initial Conditions

On Solutions of Systems of Differential Equations on Half-Line with Summable Coefficients

Discontinuous functional differential equations with delayed or advanced arguments

Infinite first order differential systems with nonlocal initial conditions

Solving first order differential equations presenting elementary functions

Novel criterion for the existence of solutions with positive coordinates to a system of linear delayed differential equations with multiple delays

On the numerical solution of non-linear first order ordinary differential equation systems

Linear determining equations, differential constraints and invariant solutions

On the solvability of boundary value problems for linear differential-algebraic equations with constant coefficients