def Integrand( gamma, x ,A, eCrit, eta, index, gammaTh): try: val = EDist(A,gamma,eta,gammaTh,index) * synchrotron_1(x/(eCrit*gamma*gamma))[0] except pygsl.errors.gsl_Error, err: #print err val = 0.
def fastInt(gamma, x, A, gamma_pl, eStar, index): try: val = fastEDist(A,gamma,gamma_pl,index) * synchrotron_1(x/(eStar*gamma*gamma))[0] except pygsl.errors.gsl_Error, err: #print err val = 0.
def Integrand( gamma, x ,A, eCrit, eta, index, gammaTh): try: val = EDist(A,gamma,eta,gammaTh,index) * synchrotron_1(x/(eCrit*gamma*gamma))[0] except pygsl.errors.gsl_Error, err: # print err val = 0.
def fastInt(gamma, x, A, gamma_pl, eStar, index): try: val = fastEDist(A, gamma, gamma_pl, index) * synchrotron_1( x / (eStar * gamma * gamma))[0] except pygsl.errors.gsl_Error, err: print err val = 0.