Abstract. Long-time stability and convergence properties of two time-discretization methods for an integro-differential equation of parabolic type are studied. The methods are based on the standard backward Euler and second-order backward differencing methods. The memory term is approximated by a quadrature rule, with emphasis on such rules with reduced computational memory requirements. Discretization of the spatial partial differential operators by the finite element method is also considered.