“…where ε ∈ (0, 1] is a small parameter, and b u (x, u) has a positive lower bound, which implies that the solution u(x) generally has an exponential boundary layer at x = 0 [18,31]. As a result, the accuracy of classical numerical methods depends not only on the number of the meshnodes N but also on the parameter ε (e.g., for (1.1) the error depends on N/ε) [22,26].…”