This numerical study presents the diagonal block method of order four for solving the second-order boundary value problems (BVPs) with Robin boundary conditions at two-point concurrently using constant step size. The solution is obtained directly without reducing to a system of first-order differential equations using a combination of predictor-corrector mode via shooting technique. The shooting method was adapted with the Newton divided difference interpolation approach as the strategy of seeking for the new init…