This manuscript presents a step by step guide on derivation and analysis of a new numerical method to solve initial value problem of fourth order ordinary differential equations. The method adopted hybrid techniques using power series as the basic function. Collocation of the fourth derivatives was done at both grid and off-grid points. The interpolation of the approximate function is also taken at the first four points. The complete derivation of the new technique is introduced and shown here, as well as the full analysis of the method. The discrete schemes and its first, second, and third derivatives were combined together and solved simultaneously to obtain the required 32 family of block integrators. The block integrators are then applied to solve problem. The method was tested on a linear system of equations of fourth order ordinary differential equation in order to check the practicability and reliability of the proposed method. The results are displaced in tables; it converges faster and uses smaller time for its computations. The basic properties of the method were examined, the method has order of accuracy p=10, the method is zero stable, consistence, convergence and absolutely stable. In future study, we will investigate the feasibility, convergence, and accuracy of the method by on some standard complex boundary value problems of fourth order ordinary differential equations. The extension of this new numerical method will be illustrated and comparison will also be made with some existing methods.