The Hamiltonian is now correct although we should test it for 3 bodies
not just 2. Also left in is the code to do traces of matrices. This
will be deleted but we will have a record of it just in case it isn't
accepted into repa.
Remove the code to calculate the trace of a matrix. We are almost
certainly *not* going to need it.