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.

Systems of linear ordinary differential equations with bounded coefficients may have very oscillating solutions

Non-oscillation of linear differential equations with coefficients containing powers of natural logarithm

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

Numerical solution of highly oscillatory ordinary differential equations

Boundedness of Solutions of Nonlinear Differential Equations

Bounded and exponentially decaying solutions of almost linear dynamic systems on time scales

On the Boundedness of Solutions of a Rational System with a Variable Coefficient

Boundedness of Solutions of A Class of Nonlinear Systems

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

Infinitely Many Periodic Solutions of Ordinary Differential Equations

Non-oscillation of modified Euler type linear and half-linear differential equations

Oscillatory criteria of noncanonical even-order differential equations with a superlinear neutral term

Oscillation for systems of nonlinear neutral type parabolic partial functional differential equations

Oscillation for Systems of a Class of Second Order Partial Eifferential Equations with Deviating Arguments

Oscillation of Solutions of Linear Differential Equations

Zero distribution of solutions of some linear differential equations

Oriented Rotatability Exponents of Solutions to Homogeneous Autonomous Linear Differential Systems

A Note on the Stability of Exponential Dichotomy of Linear Differential Equations

On a biharmonic coupled system with non-standard nonlinearity: Existence, blow up and numerics

Oscillation Result for Some Non-homogeneous Linear Differential Equations and Its Application

Sharp estimates for spectrum of Cauchy operator and periods of solutions of Lipschitz differential equations