In this paper, we consider the integration of systems of second-order linear inhomogeneous initial value problems with constant coefficients. Hybrid Numerov methods are used that are constructed in the sense of Runge-Kutta ones. Thus, the Taylor expansions at the internal points are matched properly in the final expression. We present the order conditions taking advantage of the special structure of the problem at hand. These equations are solved using differential evolution technique, and we present a method with algebraic order eighth at a cost of only 5 function evaluations per step. Numerical results over some linear problems, especially arising from the semidiscretization of the wave equation indicate the superiority of the new method.