In this paper, we discuss the numerical simulation for a class of constrained optimal control problems governed by integral equations. The Galerkin method is used for the approximation of the problem. A priori error estimates and a superconvergence analysis for the approximation scheme are presented. Based on the results of the superconvergence analysis, a recovery type a posteriori error estimator is provided, which can be used for adaptive mesh refinement.