A nonlinear singularly perturbed boundary value problem depending on a parameter is considered. Two numerical methods are applied to solve this problem. First, we solve the problem using backward Euler finite difference scheme on layer adapted meshes. Then, Richardson extrapolation technique is applied to improve the accuracy of the computed solution. Numerical experiment is carried out to validate the theoretical estimates. © 2015 The Authors.