“…In the last two decades, many published papers have focused on solving (1.1) by analytical (e.g., Adomian decomposition method, homotopy analysis method, variational iteration methods, Laplace transform decomposition method or differential transformation method) and numerical (e.g., B-spline method, the finite difference method, weighted residual method, the shooting method, multigrid-based methods, the Sinc-Galerkin method, collocation methods based on B-spline basis functions, Bessel collocation method) methods, see [13], [14], [15], [51], [52], [7], [8], [9], [16], [17], [29], [32], [35],... An extension of Bratu's problem is the following boundary value problem, socalled Liouville-Bratu-Gelfand problem (see [26], [21], [27], [28], [30], [35], [41]) :…”