“…For example, various nonlinear differential equations have been solved using the Taylor matrix method [9], the closed-form method [10], Chebyshev polynomial approximations [11], the variational iteration method [12], the subdomain finite element method [13], the differential quadrature method [14,15], the variational iteration method [16,17], He's variational iteration method [18], the cubic B-spline scaling functions and Chebyshev cardinal functions [19], the homotopy perturbation method [20,21], the variation of parameters method [22], the Adomian decomposition method [23,24], the quintic B-spline differential quadrature method [25], the variational iteration method, the power series method [26], the Adomian decomposition method [27], the Pade series method [28], the Legendre polynomial function approximation [29], the Taylor method [30], the Chebyshev series method [31] and the modified variational iteration method [32]. Recently, the Taylor, Chebyshev, Legendre, Bernstein and Bessel matrixcollocation methods have been used [33][34][35][36][37][38][39][40][41][42][43][44][45][46][47] to solve some types of differential, integral, integro-differential-difference equations. Additionally, Sheikholeslami et al…”