Here we study an ensemble of \(M\)  realizations of disorder each with \(P\) 
eigenvalues.
For the SVD analysis we construct
a matrix \(X\)  of size \(M \times P\)  where \(X_{mp}\)  is the \(p\)  level of
the \(m\) -th realization. After carrying out SVD on \(X\) ,
the singular values squared \(\lambda _k\)  are ranked from the largest
to the smallest. This approach has been applied to the
spectrum of disordered systems in several studies
[1]}, [2]}, [3]}, [4]}.
As is usual in the SVD analysis in these studies
the first few \(\lambda _k\) 
(\(k \le O(1)\) ) correspond to global features of the spectra.
Larger SV (\(\lambda _k\) ) show a power law behavior \(k^{-\alpha }\)  with \(\alpha =2\) 
at the Poisson regime and \(\alpha =1\)  for the Wigner regime.
