We present a block algorithm for the general solution of y f  x, y, y, y . The numerical algorithmis developed by the methods of interpolation of the power series approximate solution and collocation of thedifferential system of the approximant at selected grid points to generate a continuous method. Block method isadopted to simultaneously generate all the…