Added into density the ability of calculating the total charge of the
system via intergration of the density. (trap. rule).

At the moment the code seems to work for positive Z0 values
but not for negative values due to a bad calc. of the density.