Fix #100

Use the algorithm found in the R source to calculate the centralF pdf,
which avoids large exponentiation and hence doesn't suffer from NaN
overflow errors.