“…The main feature of using orthogonal functions is that they reduce dynamical system problems to those of solving a system of algebraic equations using differentiation or integration operational matrices. Special attention has been given to applications of the block pulse functions [38,39], spline wavelets [40], Legendre wavelets [41,42], Chebyshev polynomials [43,44], Walsh functions [45], and wavelet collocation method [46].…”