The initial-boundary value problem for a pseudo-parabolic equation exhibiting initial layer is considered. For solving this problem numerically independence of the perturbation parameter, we propose a difference scheme which consists of the implicit-Euler method for the time derivative and a central difference method for the spatial derivative on uniform mesh. The time domain is discretized with a nonuniform grid generated by equidistributing a positive monitor function. The performance of the numerical scheme is tested which confirms the expected behavior of the method. The existing method is compared with other methods available in the recent literature. © 2020 J. Mohapatra and D. Shakti, published by De Gruyter 2020.