A parameter robust higher order numerical method for singularly perturbed two parameter problems with non-smooth data

A singularly perturbed second order ordinary differential equation having two parameters with a discontinuous source term is presented for numerical analysis. Theoretical bounds on the derivatives, regular and singular components of the solution are derived. A hybrid monotone difference scheme with the method of averaging at the discontinuous point is constructed on Shishkin mesh. Parameter-uniform error bounds for the numerical approximation are established. Numerical results are presented which support the theoretical results.