A trigonometrically fitted block Numerov type method (TBNM), is proposed for solving y′′ = f(x, y, y′) directly without reducing it to an equivalent first order system. This is achieved by constructing a continuous representation of the trigonometrically fitted Numerov method (CTNM) and using it to generate the well known trigonometrically fitted Numerov method (TNUM) and three new additional methods, which are combined and applied in block form as simultaneous numerical integrators. The stability property of the TBNM is discussed and the performance of the method is demonstrated on some numerical examples to show accuracy and efficiency advantages.
