Numerical Solution of Stiff and Oscillatory Differential Equations Using a Block Integrator

Abstract

This paper presents the derivation and implementation of a block integrator for the solution of stiff and oscillatory first-order initial value problems of Ordinary Differential Equations (ODEs). The integrator was derived by collocation and interpolation of the combination of power series and exponential function to generate a continuous implicit Linear Multistep Method (LMM). The basic properties of the derived integrator were investigated and the integrator was implemented on some sampled stiff and oscillatory problems. From the results obtained, it is obvious that the block integrator gives better approximation than some existing ones.