Application of backward differentiation methods to the finite element solution of time-dependent problems

The application of finite element methods to parabolic partial differential equations leads to large linear systems of first-order ordinary differential equations. Very often these systems are stiff and difficulties arise in their numerical solution. We attempt to analyse the problem of how to select numerical methods for the solution of such linear systems.