version 1.4 add subroutine drawbandfull(drawband1.f90). change delta_norm to delta_a delta_b, change the particle number compute method, using qusiparticle. I am not sure this version is right.
If it is right, why the previous one doesn't work