The numerical solution of boundary value problems for stiff differential equations

Abstract

The numerical solution of boundary value problems for certain stiff ordinary differential equations is studied. The methods developed use singular perturbation theory to construct approximate numerical solutions which are valid asymptotically; hence, they have the desirable feature of becoming more accurate as the equations become stiffer. Several numerical examples are presented which demonstrate the effectiveness of these methods.