def EdNdE_gamma(x): sigma=crs.sigma_ic(x,Eph,Ee) return const.C*const.m2cm*np.dot(dN_ph,np.dot(sigma,dNe))
def sigma_fcn(x): sigma=crs.sigma_ic(x,source.Eph,Ee) #print 'sigma: ',np.size(sigma) return const.C*const.m2cm*np.dot(ph,np.dot(sigma,np.ones_like(Ee)))