  A method of eliminating the memory from the equations of motion of linear
viscoelasticity is presented. Replacing the unbounded memory by a quadrature
over a finite or semi-finite interval leads to considerable reduction of
computational effort and storage. The method applies to viscoelastic media with
separable completely monotonic relaxation moduli with an explicitly known
retardation spectrum. In the seismological Strick-Mainardi model the quadrature
is a Gauss-Jacobi quaddrature. The relation to fractional-order viscoelasticity
is shown
