Full Length Research Paper

A quadratic based integration scheme for the solution of singulo-stiff differential equations

Aashikpelokhai, U. S. U¹ and Momodu, I. B. A²*

¹Department of Mathematics, Ambrose Alli University, Ekpoma, Nigeria.

²Department of Computer Science, Ambrose Alli University, Ekpoma, Nigeria.

*Corresponding author. E-mail: bayomomoduphd@yahoo.com.

Accepted 1 April, 2008

   Abstract

In this term paper, we designed a quadratic based integration scheme for the solution of initial value problem (IVPs) in ordinary differential equation (ODEs). This was achieved by considering the rational interpolating operator

   satisfying       

A class of rational integrator formula given by

in Aashikpelokhai (1991) with K=14 was also implemented and the results compared.

The vector q= (q₁, q2………q_k) are obtained from the simultaneous linear equation (SLE) Sq=b where

P₀ = y_n,       

The results as analyzed with the computer show that the integrator copes favourably well with singular problems, stiff problem and singulo-stiff problems.

Keywords: Rational Integrators, region of absolute stability, Integration Scheme, Consistency, convergence, stiff and singular problems.