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 $\beta$ is simply a parameter in the boundary-value problem, any $\beta> 0$ can be used, in principle. The limiting distribution of the $n$th largest eigenvalue can also be computed. Our methods are available in the Julia package TracyWidomBeta.jl.