I have added the file odedirac to the wet_polsi branch so that it may be used in
conjucture with 'density'. No changes have been made to the code.

I have edited density by adding a false switch for the old version of calc
the band density and added the cal from ode. I have also changed the fudge factor
to 1.0 and Ef to -1e-4