The synthesis of an optimal control function for deterministic systems described by integrodifferential equations is investigated. By using the elegant operational properties of shifted Chebyshev polynomials, a direct computational algorithm for evaluating the optimal control and trajectory of deterministic systems is developed. An example is given to illustrate the utility of this method, and other orthogonal functions are also shown for comparison.