Modified the mardling function

The mardling function now returns the critical ratio.  For the critical
eccentricity I wrote a new function called mardling_ecc which numerically
solves the mardling function for the critical eccentricity.