We consider a hyperbolic type integro-differential equation, with weakly singular and smooth kernels. Well-posedness of the problem is studied in the framework of linear semigroup theory and Galerkin approximation method. Then numerical approximation of such problems, mainly based on finite element method, is studied and optimal order error estimates are proved
