Code restructuring, second grid.
Minor adjustments, B in file name
added rho via solution; added ode solving
fix in solve_ode
Added fitting routines
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.
rhoplotter and cutoff new python files

rhoplotter simply loads and quickly plots charge densities saved by
diracsolver.
cutoff evaluates a charge density for a number of energy windows
from the LDOS.
The last train
	arraytest.dat.npy
Merge branch 'summertime'