Numerical treatment of a quasilinear initial value problem with boundary layer

The paper deals with the singularly perturbed quasilinear initial value problem exhibiting initial layer. First the nature of solution of differential problem before presenting method for its numerical solution is discussed. The numerical solution of the problem is performed with the use of a finite-fitted difference scheme on an appropriate piecewise uniform mesh (Shishkin-type mesh). An error analysis shows that the method is first-order convergent except for a logarithmic factor, in the discrete maximum norm, independently of the perturbation parameter. Finally, numerical results supporting the theory are presented. © 2015 Taylor & Francis.