Haha, thanks for that was having a little trouble with that code.
Hmm, I think I might have been really stupid here.
Previously I calculated the eigenvalues and eigenvectors for the three spin operators ##S_x , S_y , S_z##, which I presumed was unrelated, but thinking about it they probably...