Computing the Tracy-Widom distribution for arbitrary β>0
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. 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 nth largest eigenvalue can also be computed. Our methods are available in the Julia package TracyWidomBeta.jl.
READ FULL TEXT