This paper contributes a scheme for solving linear optimal control problem in systems run by Volterra integral equations. At first, the original optimal control problem of Volterra integral equations is transformed into an integral problem via the Pontryagins maximum principle, and then, we applied homotopy analysis method (HAM) to solve the integral problem. The proposed method is compared with other methods. Numerical results show that the HAM method is more powerful and faster.