Analysis of some numerical methods on layer adapted meshes for singularly perturbed quasilinear systems

We consider a coupled system of first-order singularly perturbed quasilinear differential equations with given initial conditions. The leading term of each equation is multiplied by a distinct small positive parameter, which induces overlapping layers. The quasilinear system is discretized by using first and second order accurate finite difference schemes for which we derive general error estimates in the discrete maximum norm. As consequences of these error estimates we establish nodal convergence of O((N ^?1 lnN) ^p ),p=1,2, on the Shishkin mesh and O(N ^?p ),p=1,2, on the Bakhvalov mesh, where N is the number of mesh intervals and the convergence is robust in all of the parameters. Numerical computations are included which confirm the theoretical results.