The current paper proposes a technique for the numerical solution of linear control systems. The method is based on Galerkin method, which uses the interpolating scaling functions. For a highly accurate connection between functions and their derivatives, an operational matrix for the derivatives is established to reduce the problem to a set of algebraic equations. Several test problems are given, and the numerical results are reported to show the accuracy and eciency of this method.
