Diracsolver now with complete change in charge density.

Diracsolver finds drho including Fermi sea contribution and Thomas-
Fermi quadratic term. "Correction" factor included in raw band drho
data. High m-channel contribution included in total drho.

drhotest takes data of total drho for a number of betas, comparing
the symmetry of the data.