Debugging

For basic debugging, delete the statement ``gbl_f=2;'' from ``Geatoms.m'' and add the statement

gbl_f=2*ones(Ns,1); %# Constant fillings as a vector

just after the definition of Ns in ``Geatoms.m''. Rerun ``Geatoms.m'' and verify that your output (with perhaps changes in the fourteenth or fifteenth digits) is identical to before, at least through the finite-difference test.