This paper presents a spectral method of solving the controlled Duffing oscillator. The method is based upon constructing the Mth degree interpolation polynomials, using Chebyshevs nodes, to approximate the state and the control vectors. The differential and integral expressions that arise from the system dynamics and the performance index are converted into some algebraic equations. The optimum condition is obtained by applying the method of constrained extremum. (C) 1997 by Elsevier Science Inc.