def _int(f, a, b, args=()): #return quad(f, a, b, args=args, limit=1000, epsrel = 1e-7, # full_output=(0 if _feedback else 1))[0] #return Integrate_Patterson(f, a, b, eps = 1e-7, abs=1e-10, args=args)[0] return integrate.chebyshev(f, a, b, epsrel=1e-7, epsabs=1e-9, args=args)
def _int(f, a, b, args=()): #return quad(f, a, b, args=args, limit=1000, epsrel = 1e-7, # full_output=(0 if _feedback else 1))[0] #return Integrate_Patterson(f, a, b, eps = 1e-7, abs=1e-10, args=args)[0] return integrate.chebyshev(f, a, b, epsrel = 1e-7, epsabs=1e-9, args=args)
def _integrator(f, a, b): return integrate.chebyshev(f, a, b, epsrel = 1e-8, epsabs=1e-10)
def _integrator(f, a, b): return integrate.chebyshev(f, a, b, epsrel=1e-8, epsabs=1e-10)