“…The state and/or control involved in the equation are approximated by finite terms of orthogonal series and by using the operational matrix of integration the integral operations are eliminated. The form of the operational matrix of integration depends on the particular choice of the orthogonal functions like Walsh functions [4], Block-pulse functions [8], Laguerre series [9], Jacobi series [10], Fourier series [11], Bessel series [12], Taylor series [13], Shifted Legendre [14], Chebyshev polynomials [15] and Hermite polynomials [16]. In this study, we use wavelet functions to approximate both the control and state functions.…”