The linear, homogeneous, parabolic equation is solved by applying finite element discretizations in space and A0 —stable, linear multistep, multiderivative (L.M.S.D.) methods in time. Such schemes are unconditionally stable. An error analysis establishes an optimal bound in the L2 —norm. Methods typifying the class of L.M.S.D. schemes are derived and their implementation examined