Edited the definition of the limb darkening laws so that
the intensity at the disc center was normalised to 1 before
performing the least squares fit. There still remains a
discrepancy between the output from this code and the IDL
code of David Sing. This seems to be due to the different
results obtained with the numpy trapz() and the
IDL int_tabulated() routines. This issue is currently unresolved.