Computing the Tracy-Widom Distribution for Arbitrary 𝛽 > 0

Abstract

We compute the Tracy-Widom distribution describing the asymptotic distribution of the largest eigenvalue of a large random matrix by solving a boundary-value problem posed by Bloemendal in his Ph.D. Thesis (2011). The distribution is computed in two ways. The first method is a second-order finite-difference method, and the second is a highly accurate Fourier spectral method. Since 𝛽 is simply a parameter in the boundary-value problem, any 𝛽 > 0 can be used, in principle. The limiting distribution of the 𝑛th largest eigenvalue can also be computed. Our methods are available in the Julia package TracyWidomBeta.jl.