Previous numerical methods for solving systems of higher order ordinary differential equations (ODEs) directly require calculating the integration coefficients at every step. This research provides a block multi step method for solving orbital problems with periodic solutions in the form of higher order ODEs directly. The advantage of the proposed method is, it requires calculating the integration coefficients only once at the beginning of the integration is presented. The derived formulae is then validated by running simulations with known higher order orbital equations. To provide further efficiency, a relationship between integration coefficients of various order is obtained.

Real life phenomena found in various fields such as engineering, physics, biology and communication theory can be modeled as nonlinear higher order ordinary differential equations, particularly the Duffing oscillator. Analytical solutions for these differential equations can be time consuming whereas, conventional numerical solutions may lack accuracy. This research propose a block multistep method integrated with a variable order step size (VOS) algorithm for solving these Duffing oscillators directly. The proposed VOS Block method provides an alternative numerical solution by reducing computational cost (time) but without loss of accuracy. Numerical simulations are compared with known exact solutions for proof of accuracy and against current numerical methods for proof of efficiency (steps taken).