In this research, a six-order, fully implicit Block Backward Differentiation Formula with two off-step points (BBDFO(6)), for the integration of first-order ordinary differential equations (ODEs) that exhibit stiffness, is proposed. The order, consistency and stability properties of the method are discussed, and the method is found to be zero stable and consistent. Hence, the method is convergent. The numerical comparisons with the existing methods of a similar type are given to demonstrate the accuracy of the derived method.

Recently, block backward differentiation formulas (BBDFs) are used successfully for solving stiff differential equations. In this article, a class of hybrid block backward differentiation formulas (HBBDFs) methods that possessed $A$ –stability are constructed by reformulating the BBDFs for the numerical solution of stiff ordinary differential equations (ODEs). The stability and convergence of the new method are investigated. The methods are found to be zero-stable and consistent, hence the method is convergent. Comparisons between the proposed method with exact solutions and existing methods of similar type show that the new extension of the BBDFs improved the stability with acceptable degree of accuracy.