Dear Prof. Hiary, Please see the attached for my thesis, data files, and code. The data files are pk(0,t)@[0,4] for N0=3,4,5,6. (nearest neighboring spacing distribution) pk(1,t)@[0,5] for N0=4,5,6,7. (2nd nearest neighboring spacing distribution) pk(2,t)@[1,7] for N0=5,6,7,8. (3rd) pk(3,t)@[1,7] for N0=5,6,7,8. (4th) In particular, @[a,b] means the data is for t =a+0.05, a+0.1, ..., b. More details could be referred to Chapter 6 of my thesis. The file Distribution_p(0,t)_archive.nb is for calculating p(0,t) and comparing with the empirical data. But it can be easily modified to calculate p(k,t) for general k. The file Spheroidal6.m is Noah Graham's spheroidal function package. Please let me know if you have any questions. Best, Lin