  Hubbard model is an important model in theory of strongly correlated electron
systems. In this contribution we introduce this model along with numerically
exact method of diagonalization of the model.
