“…In [18,7], the computation of approximate solution for nonlinear singular initial value problems was considered using Hermite wavelets operational matrix of integration and Chebyshev wavelets operational matrix of integration respectively while in [9], the optimal control problems was solved with the aid of Legendre mother wavelets operational matrix of integration. A general expression for the operational matrix of integration of Bessel functions was derived in [16] to solve some examples in optimal control, for more works, see [8,14,15,16,17,18,19,22,24,25,26,27,28,29,30,31,32,33,34,35,36,37].…”