The ultimate goal of this study is to develop a numerically effective approximation technique acquire numerical solutions the integer and fractional-order Bratu singular Lane–Emden-type problems especially with exponential nonlinearity. Both initial boundary conditions were considered fractional derivative being in Liouville–Caputo sense. In direct approach, generalized Bessel matrix method bas...