“…18 In most cases, it is difficult or impossible to obtain an analytical solution for the integral equation (1). Therefore, in recent decades, there has been much attention to constructing some numerical techniques to achieve the approximate solutions with high accuracy for the VIDEs such as the homotopy perturbation method, 19 variational iteration method (VIM), 20,21 operational matrix method based on orthogonal polynomials 1,[22][23][24][25][26][27] operational tau method, 28,29 collocation approach, [30][31][32] block-pulse functions method, 33,34 radial basis function method, 3,35 wavelet method, 36,37 hat functions method, 38 least squares method, 10,39,40 hybrid functions method, [41][42][43][44] fixed point method, 45 successive approximations method, [46][47][48] Euler polynomials, 49,50 delta functions, 51 Taylor method, 52 and so on.…”