start to add a script that does spectral decomposition

for small T.  This can be used to precondition the matrix
the spectral decomposition now compute eigevectors

it doesn't agree yet with the directo computation.
Merge branch 'master' of github.com:js850/kmc_rates