“…This phenomenon determines the evolution of parameter-uniform numerical methods i.e., the methods in which the error constant is independent of ε and of the mesh parameter. Various ε-uniform numerical schemes such as the variational method, the finite difference methods (FDMs), the rational spectral collocation methods, the finite element methods (FEMs), the adaptive mesh methods, and the layer-adapted mesh methods have been developed in the literature for singularly perturbed boundary value problems (SPBVPs) (readers are referred to [1,[11][12][13][18][19][20]32] and the references therein). Although the Shishkin mesh is one of the simplest non-uniform meshes; it has a drawback, that is before one attempts to solve the differential equation, significant information about the exact solution must be known.…”