“…It is a fact that the analytical solution of a classical fractional differential equation is difficult to obtain, thus many numerical methods have been developed such as finite difference method, 13 finite element method, 14 spectral method, 15 nonpolynomial spline method, 16‐21 and other methods involving fractional‐order Lagrange polynomials, 22 Legendre‐Laguerre polynomials, 23 Genocchi hybrid functions, 24 Bernoulli wavelets, 25 to name a few. As one would expect, it is even more difficult to obtain analytical solutions of equations involving generalized fractional derivatives, and numerical treatment is more appropriate for such problems.…”