In this talk I will give a brief outline of the Grassmann integration technique (which is also called the supersymmetry approach) and its application to the rigorous study of main spectral characteristics of random matrices. I am going to start from a simplest case of GUE matrices, and then to show the relation between characteristic polynomials of random band matrices and some questions about the classical Heisenberg model.