Numerically integrate whole area

Added two lines to integrate full area on the request of Oksana and Romain, might be changed back after some more testing, ergo I did not delete the previous code and only change the final result regardless of what the code produces. I find this least invasive, but a bit messy.

pyzebra/fit2.py

@@ -191,6 +191,8 @@ def fitccl(
         else:
             pass
               
+        numfit_min = 0
+        numfit_max = len(x)
         count_errors = create_uncertanities(y, y_err)
         # create error vector for numerical integration propagation
         num_int_area = simps(count_errors[numfit_min:numfit_max], x[numfit_min:numfit_max])