“…Its non-linearity trait makes it very difficult to find an exact solution, and the discrete point findings of numerical approaches do not provide consistency. To address this problem, a number of researchers have provided techniques, such as the sine Gordon expansion method [16][17][18][19], the Homotopy analysis method [20][21], the Bernoulli sub-equation function method [22], the generalized exponential rational function method [23], the He's variational iteration method [24][25][26], the Homotopy perturbation technique [27][28][29], the Chebyshev spectral collocation method [30], the rational sine-cosine and rational sinh-cosh methods [31], the conformable derivative [32], the timefractional Caputo derivative [33], the Atangana-Baleanu derivatives [34], the fractional residual power series method [35], the Residual Power Series Method [36], semi-analytic technique to deal with nonlinear fractional differential equations [37]. For more details we can refer to [38].…”