Hermite and Laguerre β-ensembles: asymptotic corrections to the eigenvalue density

We consider Hermite and Laguerre $\beta$-ensembles of large $N\times N$ random matrices. For all $\beta$ even, corrections to the limiting global density are obtained, and the limiting density at the soft edge is evaluated. We use the saddle point method on multidimensional integral representations of the density which are based on special realizations of the generalized (multivariate) classical orthogonal polynomials. The corrections to the bulk density are oscillatory terms that depends on $\beta$. At the edges, the density can be expressed as a multiple integral of the Konstevich type which constitutes a $\beta$-deformation of the Airy function. This allows us to obtain the main contribution to the soft edge density when the spectral parameter tends to $\pm\infty$.