Abstract: This paper is in order to solve the numerical computation of the periodic orbits of nonlinear autonomous differential equations. We design an optimal control method of the numerical computation. The solution method is basically a continue shooting method,under the guardianship of nonlinear least squaves that takes advantage of an optimal Central formulation, in which attached to the contral techniques of the penalty function.In this method,we use the Crank-Nicolson discretization schemes and Conjugate Gradient methods.This method has been applied to the solution of nonlinear differential equations modelling chemical reactions,the corresponding numerical results are presented here.