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