We present an analytic proof of the existence of phase transition in the large N limit of certain random noncommutaitve geometries. These geometries can be expressed as ensembles of Dirac operators. When they reduce to single matrix ensembles, one can apply the Coulomb gas method to find the empirical spectral distribution. We elaborate on the nature of the large N spectral distribution of the Dirac operator itself. Furthermore, we show that these models exhibit both a single and double cut region for certain values of the order parameter and find the exact value where the transition occurs.