A two-stage fourth-order “almost” P-stable method for oscillatory problems

We propose a two-stage fourth-order “almost” P-stable method for y″ = f(t,y) based on a pair of discretizations in which the displacement is computed using the trapezoidal rule with end point correction and the velocity is computed using suitably modified Simpsonʹs rule.