Parabolic systems: finite element method