“…In the present years, several numerical methods for approximating the solution of delay differential equations have been established. Many authors have tried to find the solutions of delay differential equations using various techniques such as a collocation method based on Bernoulli operational matrix [10], Taylor polynomials [11,12,13,14], Euler bases together with operational matrices [15], perturbation-iteration algorithms [16], Laguerre series [17], Walsh stretch matrices and functional differential equation [18], Bernstein polynomials [19], Fourier operational matrices of differentiation and transmission [21] polynomial interpolation [20], Spline functions approximation [22], Adomian decomposition method [23,24], Hermite interpolation [25], collocation method [26], Chebyshev polynomials [27], Legendre polynomial approximation [28], differential transform method [29], block-pulse functions and Bernstein polynomials [30], Variational iteration method(VIM) [31,32], Jacobi rational-Gauss collocation (JRC) [33], successive interpolations [34], an efficient transferred Legendre pseudospectral method [35], Muntz-Legendre basis and operational matrices of fractional derivatives [36] etc.. Methods based on the wavelets are more attractive and considerable. Some of wavelets techniques are applied in order to solve the equation (1) namely, Chebyshev wavelets [37,38,39], Hermite wavelets [40,41], Legendre wavelets method [42,43], Haar wavelets method [44,…”