Make a copy of ``poisson2.m'' called ``poissonGe.m''. In ``poissonGe.m'', remove all references to the second Gaussian at (diag(R)'/2+[x,0,0]), keeping all terms pertaining to the first Gaussian in all mathematical expressions. Be sure to keep the printouts checking the normalization of the Gaussian and the total charge in the cell, and to keep accounting for the self-energy part for the remaining Gaussian. Finally, change the norm on the remaining Gaussian at (diag(R)'/2) so that it represents a net charge of four instead of one.