I save the hamiltonian matrix to a file.
A routine to calculate the mean phonons.