def NNSD2D():
a = 1
b = np.sqrt(np.pi / 3)
parameters = Box2D.Parameters(a, b)
title = "2D box\n" + str(parameters)
spectrum = Box2D.CalculateSpectrum(parameters, 1000000)
unfolded = StretchSpectrum(spectrum)
spacings = LevelSpacing(unfolded)
PlotLevelDensity(spectrum,
lambda x: Box2D.LevelDensity(parameters, x),
extraTitle=title)
PlotCummulativeLevelDensity(
spectrum,
lambda x: Box2D.CummulativeLevelDensity(parameters, x),
extraTitle=title)
PlotNNSD(spacings, Poisson, extraTitle=title)
