Better algorithm for estimating second derivative

Makes a little bit of difference. Extracted from

http://mathformeremortals.wordpress.com/2013/01/12/a-numerical-second-derivative-from-three-points/
Maximum radius changed
Useless variables removed
Bug in the potential calculation fixed