The authors suggest a numerical method to solve an optimal control problem for the Duffing oscillator. The base of the approach is a cell averaging method producing Chebyshev interpolating polynomials for the approximation of the state and control vectors. The system dynamics and the performance index are converted into some algebraic equations and their approximation is obtained. Tables illustrating numerical results are presented.