Here we have the new file for solving the Dirac equation, using functions.

The function finding the eigenvectors/values appears to work reasonably,
producing a constant charge density when all negative energy states are
included. The density of states function appears to be flawed.