International Journal of
Physical Sciences

  • Abbreviation: Int. J. Phys. Sci.
  • Language: English
  • ISSN: 1992-1950
  • DOI: 10.5897/IJPS
  • Start Year: 2006
  • Published Articles: 2572

Full Length Research Paper

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

Aashikpelokhai, U. S. U1 and Momodu, I. B. A2*
1Department of Mathematics, Ambrose Alli University, Ekpoma, Nigeria. 2Department of Computer Science, Ambrose Alli University, Ekpoma, Nigeria.
Email: [email protected]

  •  Accepted: 01 April 2008
  •  Published: 30 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= (q1, q2………qk) are obtained from the simultaneous linear equation (SLE) Sq=b where

           

           

 

 

 

P0 = yn       

 

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.