In this talk we explain the time integration of parabolic equations with block implicit methods (BIM). Depending on the size of the block, high-order BIM with A-stability are designed without the need of multiple initial guesses. Similar to Runge-Kutta methods, a BIM can be defined by a tableau including two matrices and two vectors. We discuss a special scheme defined by a positive definite matrix and a positive diagonal matrix; both matrix properties are desirable but not available in Runge-Kutta methods. Moreover, we show that the traditional finite element theory for parabolic problems discretized by the backward Euler or Crank-Nicolson schemes can also be extended to BIM.