An efficient non-uniform Haar wavelet method is proposed for the numerical solution of system of first order linear partially singularly perturbed initial value problem on piecewise uniform Shishkin mesh and ?-mesh. Further, ...