An Advanced Numerical Approach for Solving Stiff Initial Value Problems Using a Self-Starting Two-Stage Composite Block Scheme

Document Type : Research Paper

Authors

Faculty of Computer and Mathematical Sciences, Universiti Teknologi MARA, 40450 Shah Alam, Selangor, Malaysia.

Abstract

In this study, a two-point composite block method based on the backward differentiation formula (CBBDF) is introduced to solve stiff ordinary differential equations. The CBBDF method incorporates an additional intermediate point among the interpolating points, developed in two stages: the first stage employs the Euler’s method as a fundamental building block, while the second stage utilizes CBBDF of order three. A key distinction of the method with the classical block method is the introduction of an independent parameter γ, which eliminates the need for an external startup calculation, while maintaining the accuracy and stability of numerical solutions. The theoretical analysis verifies that the proposed method is convergent and A-stable. It fulfills the essential properties of consistency and zero-stability, and it lies within the A-stability region. To demonstrate the effectiveness of the proposed approach, several stiff initial value problems of linear and non-linear are solved. For validation, the results are compared with existing literature. While approximating the solution at multiple points simultaneously, the composite block method offers the ability to use larger step sizes for solution approximation. The CBBDF method shows promising results, achieving a reliable degree of accuracy as indicated by its maximum error and average error measurements.

Keywords

Main Subjects



Articles in Press, Accepted Manuscript
Available Online from 02 May 2025
  • Receive Date: 08 January 2025
  • Revise Date: 26 March 2025
  • Accept Date: 23 April 2025