A computational method for solving singularly perturbed turning point problems exhibiting twin boundary layers

摘要

Singularly perturbed turning point problems (TPPs) for second order ordinary differential equations (DEs) exhibiting twin boundary layers are considered. In order to obtain numerical solution of these problems a computational method is suggested in which exponentially fitted difference schemes are combined with classical numerical methods. In this method, the given internal (domain of definition of differential equation) is divided into four subintervals and we solve the differential equation in each interval and combine the solution. The implementation of the method on parallel architectures is mentioned. Error estimates are derived and numerical examples are given.