Abstract: We consider an approximation scheme for integro-partial differential equations which arise in the theory of linear viscoelasticity. This scheme is based on a modification (to account for the singular kernel) of certain averaging type approximation methods for delay equations. We use this scheme to investigate the effects of a history parameter (the delay length) on the behavior of the eigenvalues, and to consider the numerical solution of an optimal control problem.