  We consider a wide range of matrix models and study them using the Monte
Carlo technique in the large $N$ limit. The results we obtain agree with exact
analytic expressions and recent numerical bootstrap methods for models with one
and two matrices. We then present new results for several unsolved multi-matrix
models where no other tool is yet available. In order to encourage an exchange
of ideas between different numerical approaches to matrix models, we provide
programs in Python that can be easily modified to study potentials other than
the ones discussed. These programs were tested on a laptop and took between a
few minutes to several hours to finish depending on the model, $N$, and the
required precision.
