“…Most of these equations do not have analytical solutions, but numerical techniques may be used to obtain approximate solutions. 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, integ...…”