A family of improved Falkner-type methods for oscillatory systems

For solving general second-order initial value problems u′′(t)=f(t,u(t),u′(t)),u′′(t)=f(t,u(t),u′(t)), the classical Falkner methods can date back to Falkner’s work in 1936. In this paper, we propose and study a family of improved Falkner-type methods for the oscillatory system u′′(t)+Mu(t)=g(t,u(t)),u′′(t)+Mu(t)=g(t,u(t)), where g: R × Rd → Rd, in which the first derivative does not appear explicitly, and M ∈ Rd × d is a symmetric positive semi-definite matrix. The new methods take into account the oscillatory structures of the problem and exactly integrate the unperturbed problem u′′(t)+Mu(t)=0u′′(t)+Mu(t)=0. The global error bounds of the new methods are presented. Numerical experiments are performed to show that the new methods are more efficient than other effective methods appeared in the scientific literature.