EN
It is shown that mean values of electron operators in variational Gutzwiller state are equal to mean values of corresponding classical quantities calculated by means of a hermitian matrix model. In cases with small number of electrons in the system this property enables exact calculation of the mean values. In case of large number of electrons a simple and effective Monte Carlo method is formulated (within matrix model).