Stiff systems of ordinary differential equations (ODEs) arise frequently in science, engineering, and applied mathematics, and their numerical solutions require methods with strong stability properties to avoid severe step-size restrictions. This motivates the development of efficient and stable numerical schemes tailored for such problems. In this paper, a fully implicit hybrid block method, termed the Two-Point Superclass Block Backward Differentiation Formula with Off-step Points (2SBBDFO), is proposed for solving stiff ODE systems. The method extends the classical two-point block BDF by incorporating a non-zero coefficient $\beta_{k-1}$, enabling the simultaneous computation of two grid-... More >

Stiff systems of ordinary differential equations (ODEs) arise frequently in several areas of applied mathematics, science, and engineering, and their numerical integration requires methods with strong stability characteristics. In this paper, a block extended backward differentiation formula (BEBDF) with a stability control parameter is developed for the numerical solution of stiff systems of ODEs. The proposed method is constructed within the framework of block multistep methods, allowing the simultaneous computation of solution values at several grid points. A free parameter $\rho$ is incorporated into the formulation of the method to regulate and improve its stability behavior. The presen... More >

This paper presents an adaptive block time-stepping algorithm designed for the numerical integration of first-order stiff oscillatory problems, which typically involve multiple time scales and stringent stability conditions. The proposed method integrates a block formulation enabling the simultaneous evaluation of several solution points with an adaptive step size mechanism to effectively regulate local truncation errors. A predictor-corrector strategy is employed, where an explicit predictor generates initial estimates, and an implicit block corrector enhances stability for stiff systems. The resulting nonlinear equations are resolved using Newton's iterative method. The scheme is proven to... More >

This paper develops a third-order fully implicit adaptive variable step-size superclass block backward differentiation formula (VSBBDF3) for the efficient numerical simulation of nonlinear stiff dynamical systems, including oscillatory, chaotic, and reaction-kinetics systems governed by ordinary differential equations. The method extends the classical block BDF framework through a structured superclass coefficient formulation while preserving full implicitness. By computing multiple solution approximations simultaneously within each block, the scheme enhances both stability and computational efficiency. An adaptive step-size strategy controls local truncation errors, enabling dynamic respons... More >

In this paper, we present a modified three-point superclass of Block Backward Differentiation Formula (BBDF) for the efficient numerical solution of stiff systems of ordinary differential equations (ODEs). The principal enhancement of this work is a structural modification of the classical BBDF that forms a new parameterized superclass of methods, leading to improved stability and reduced error constants compared with the standard three-point BBDFs. The proposed scheme is formulated as a fully implicit block method capable of simultaneously producing three solution approximations within each integration step. A detailed theoretical analysis is conducted to establish the order of accuracy, co... More >