ABSTRACT
In this paper, a new explicit fourth-order Runge–Kutta method adapted to oscillatory problems with two frequencies (μ and ν) is developed. The new method can integrate exactly the different harmonic oscillators and
simultaneously. Numerical stability and phase properties of the new method are analysed. Numerical experiments are conducted to illustrate the high accuracy and competence of our new method in comparison with some highly effective methods from the recent literature.
Disclosure statement
No potential conflict of interest was reported by the authors.
ORCID
Yonglei Fang http://orcid.org/0000-0002-5615-5971
Yanping Yang http://orcid.org/0000-0002-8245-334X
Xiong You http://orcid.org/0000-0002-0343-3878