Stability Analysis of Singly Diagonally Implicit Block Backward Differentiation Formulas for Stiff Ordinary Differential Equations

Abstract

In this research, a singly diagonally implicit block backward differentiation formulas (SDIBBDF) for solving stiff ordinary differential equations (ODEs) is proposed. The formula reduced a fully implicit method to lower triangular matrix with equal diagonal elements which will results in only one evaluation of the Jacobian and one LU decomposition for each time step. For the SDIBBDF method to have practical significance in solving stiff problems, its stability region must at least cover almost the whole of the negative half plane. Step size restriction of the proposed method have to be considered in order to ensure stability of the method in computing numerical results. Efficiency of the SDIBBDF method in solving stiff ODEs is justified when it managed to outperform the existing methods for both accuracy and computational time.