“…In the literature, there are many definitions of fractional derivative, the most popular definitions are of Riemann-Liouville, Liouville-Caputo, Caputo-Fabrizio, Atangana-Baleanu, Riesz, Hilfer, among others [8][9][10]. Recently, several numerical methods have been proposed to obtain approximate solutions of fractional ordinary differential equations and fractional partial differential equations, such as the fractional sub-equation method [11,12], the Adomian decomposition method [13][14][15], the Homotopy perturbation method [16,17], the variational iteration method [18][19][20][21], homotopy perturbation transform method [22,23], and so on.…”