Power Law Distribution Fitting in python (and fortran and cython)