Analysis and implementation of a computational technique for a coupled system of two singularly perturbed parabolic semilinear reaction-diffusion equations having discontinuous source terms.