Numerical integration of differential viscoelastic models