A variable order one-step scheme for numerical solution of ordinary differential equations

AbstractThe author proposes some stable and convergent two-point integration formulae which are particularly well suited to systems of ordinary differential equations with oscillating solutions. The numerical integration algorithms are based on the representation of the theoretical solution by the perturbation of a polynomial interpolating function with a trigonometric function. For non-oscillatory systems, the proposed schemes reduce to the normal Taylor series