The method proposed in this paper uses direct approach to the control system optimization. In order to reduce variational problem to ordinary extremum problem all state and control variables are approximated by a particular from of power polynomials in which the polynomial coefficients are given by linear combinations of all coordinate of the output state signal xN(t).