Error estimates in L∞(0,T;L 2(Ω)), L∞(0,T;L 2(Ω)2),L∞(0,T;L ∞(Ω)), L∞(0,T;L∞ (Ω)2), Ω in  ${\mathbb R}^2$  , are derived for a mixed finite element method for the initial-boundary value problem for integro-differential equation  $$u_t={\rm div}\{a\bigtriangledown u+\int^t_0b_1\bigtriangledown u{\rm d}\tau+\int^t_0{\bf c}u{\rm d}\tau\}+f$$
 , are derived for a mixed finite element method for the initial-boundary value problem for integro-differential equation  $$u_t={\rm div}\{a\bigtriangledown u+\int^t_0b_1\bigtriangledown u{\rm d}\tau+\int^t_0{\bf c}u{\rm d}\tau\}+f$$  based on the Raviart-Thomas space V h  x Wh  ⊂ H(div;Ω) x L2(Ω). Optimal order estimates are obtained for theapproximation of u,ut  in L∞(0,T;L 2(Ω)) and the associated velocity p in L∞(0,T;L 2(Ω)2), divp in L∞(0,T;L 2(Ω)). Quasi-optimal order estimates are obtained for the approximation of u in L∞  (0,T;L ∞(Ω)) and p in L∞(0,T;L ∞(Ω)2.
 based on the Raviart-Thomas space V h  x Wh  ⊂ H(div;Ω) x L2(Ω). Optimal order estimates are obtained for theapproximation of u,ut  in L∞(0,T;L 2(Ω)) and the associated velocity p in L∞(0,T;L 2(Ω)2), divp in L∞(0,T;L 2(Ω)). Quasi-optimal order estimates are obtained for the approximation of u in L∞  (0,T;L ∞(Ω)) and p in L∞(0,T;L ∞(Ω)2.