In this paper, we derive an algorithm to solve the linear quadratic (LQ) optimal regulator problems. The new approach is based on efficient Legendre and Chebyshev formulae at the Chebyshev-Gauss-Lobatto points. The technique enjoys advantages of both the Legendre and Chebyshev approximations near the end points. To show the numerical behavior of the proposed method, the simulation results of an example are presented. [ABSTRACT FROM AUTHOR]