We generally study the density of eigenvalues in unitary ensembles of random matrices from the recurrence coefficients with regularly varying conditions for the orthogonal polynomials. First we calculate directly the moments of the density. Then, by studying some deformation of the moments, we get a family of differential equations of first order which the densities satisfy