Modifications to find the error. It isn't there, the numerical
calculations are correct.

A few modifications to compare my data with Juns.

smp_rspheres now allows for the centres to be specified. A plain number
can still be specified though.