Solving Directly Higher Order Ordinary Differential Equations by Using Variable Order Block Backward Differentiation Formulae

Variable order block backward differentiation formulae (VOHOBBDF) method is employed for treating numerically higher order Ordinary Differential Equations (ODEs). In this respect, the purpose of this research is to treat initial value problem (IVP) of higher order stiff ODEs directly. BBDF method is symmetrical to BDF method but it has the advantage of producing more than one solutions simultaneously. Order three, four, and five of VOHOBBDF are developed and implemented as a single code by applying adaptive order approach to enhance the computational efficiency. This approach enables the selection of the least computed LTE among the three orders of VOHOBBDF and switch the code to the method that produces the least LTE for the next step. A few numerical experiments on the focused problem were performed to investigate the numerical efficiency of implementing VOHOBBDF methods in a single code. The analysis of the experimental results reveals the numerical efficiency of this approach as it yielded better performances with less computational effort when compared with built-in stiff Matlab codes. The superior performances demonstrated by the application of adaptive orders selection in a single code thus indicate its reliability as a direct solver for higher order stiff ODEs.

[1]  Direct fifth order block backward differentiation formulas for solving second order ordinary differential equations , 2016 .

[2]  The Numerical Computation of Three Step Hybrid Block Method for Directly Solving Third Order Ordinary Differential Equations , 2017 .

[3]  Zarina Bibi Ibrahim,et al.  Fixed Coeficient A(α) Stable Block Backward Differentiation Formulas for Stiff Ordinary Differential Equations , 2019, Symmetry.

[4]  B. Sumithra,et al.  Numerical solution of stiff system by backward Euler method , 2015 .

[5]  M. Suleiman,et al.  Solution of Higher-Order ODEs Using Backward Difference Method , 2011 .

[7]  C. W. Gear.  The numerical integration of ordinary differential equations , 1967 .

[8]  F. Ismail,et al.  Fourth-Order Improved Runge–Kutta Method for Directly Solving Special Third-Order Ordinary Differential Equations , 2015 .

[9]  Fudziah Bt. Ismail,et al.  A new variable step size block backward differentiation formula for solving stiff initial value problems , 2013, Int. J. Comput. Math..

[10]  HIGH ORDER BLOCK IMPLICIT MULTI-STEP (HOBIM) METHODS FOR THE SOLUTION OF STIFF ORDINARY DIFFERENTIAL EQUATIONS , 2014 .

[11]  Z. Ibrahim,et al.  Adaptive order of block backward differentiation formulas for stiff ODEs , 2017 .

[12]  R. Ponalagusamy,et al.  A Parallel Fourth Order Rosenbrock Method: Construction, Analysis and Numerical Comparison , 2015 .

[13]  Ali Abdi,et al.  Construction of high-order quadratically stable second-derivative general linear methods for the numerical integration of stiff ODEs , 2016, J. Comput. Appl. Math..

[14]  Asma Izzati Asnor,et al.  Algorithm of modified variable step block backward differentiation formulae for solving first order stiff ODEs , 2018 .

[15]  C. W. Gear,et al.  The automatic integration of stiff ordinary differential equations , 1968, IFIP Congress.