Final debugging

Repeat the above calculations in a larger cell, R=diag([12 12 12]), with a corresponding larger number of sampling points, S=[48; 48; 48]. You should find significantly improved agreement for both the total energy and the s-p energy difference.
Hint: Don't forget to rerun ``poissonGe.m'' for the larger cell as well!