We introduce a new matrix model, which counts exactly (not assymptoticaly) plane partitions and random lozenge tilings of an arbitrary domain. This provides a simple explanation of why asymptotic properties of those models are always some random matrix limit laws. We will also briefly mention how to compute the matrix integral, and make the link with algebraic geometry