Translating Szabo and Ostlund HF code to Python
Incomplete, need to experiment with different diagonalisation methods.
