In this paper, a numerical method for solving the optimal control （OC） problems is presented. The method is enlightened by the Chebyshev-Legendre （CL） method for solving the partial differential equations （PDEs）. The Legendre expansions are used to approximate both the control and the state functions. The constraints are discretized over the Chebyshev-Gauss-Lobatto （CGL） collocation points. A Legendre technique is used to approximate the integral involved in the performance index. The OC problem is changed into an equivalent nonlinear programming problem which is directly solved. The fast Legendre transform is employed to reduce the computation time. Several further illustrative examples demonstrate the efficiency of the proposed method.
