Make symmetrization work

The generation of matrices B and H will now include proper
symmetrization. The matrices will be symmetric and be able to run
through the eigenValue algorithm unlike previous symmetrization
Fixed a off by one error in coulombpotential and fouled around a little