  We demonstrate a method to solve a general class of random matrix ensembles
numerically. The method is suitable for solving log-gas models with
biorthogonal type two-body interactions and arbitrary potentials. We reproduce
standard results for a variety of well-known ensembles and show some new
results for the Muttalib-Borodin ensembles and recently introduced
$\gamma$-ensemble for which analytic results are not yet available.
