Chiang Mai Journal of Science

A code based on predictor-corrector 2-point block of Backward Differntiation Formulas using variable step size is developed for solving second order Ordinary Differntial Equations (ODEs). The new method solved the second order differential equations directly without reduction to its first order differential equations. Further, the new method allows the differentiation coefficients to be stored and thus avoiding the coefficients to be calculated at each step but robust enough to allow for step size variation. Computational advantages are presented comparing the results obtained by the new method with solver ode 15s and ode 23s in MATLAB. Numerical results prove that the new developed method is reliable and efficient.