“…This can be made possible by the application of suitable numerical methods. Researchers like [6,3,1,13,10,11,14,20] , etc have proposed different numerical methods for solving IVPs in ODEs but some of these methods are not self-starting methods and required single step methods to generate other starting values which is computationally cumbersome and prone to enormous error since the starting method and the method are not of the same order. Our interest in this work is to formulate a family of block methods with optimized region of stability suitable for stiff and non-stiff Initial Value problems in Ordinary Differential Equations of the form π¦ β² = π(π₯), π¦(π₯ 0 ) = π¦ 0, π₯ β [π, π]…”