  The binding energy spectra of the heavy quarkonia are calculated by solving
the Schr\"odinger equations with Coulomb plus confining potentials. Statistical
properties of the obtained spectra are studied by plotting nearest level
spacing distribution histograms.
