The bounded energy optimal control problem, for a class of heat conduction systems, is posed. The response of the system is approximated by using the method of moments. This approximation serves for constructing a non-linear integral equation, the solutions of which tire the optimal controls. A method for solving this equation is described and on example illustrates the computational scheme.