def integrand(x):
return -pow(ec, ap1) * gamma_inc(ap1, old_div(x, ec))
def _integral(a, b, index, ec):
ap1 = index + 1
integrand = lambda x: -pow(ec, ap1) * gamma_inc(ap1, x / ec)
return integrand(b) - integrand(a)
#f15
gf15 = lambda x, y: sf.poch(x, y)
rf15 = lambda x, y: gamma(fadd(x, y, exact=True)) / gamma(x)
#f16
gf16 = lambda x, y: sf.lnpoch(x, y)
rf16 = lambda x, y: fsub(
loggamma(fadd(x, y, exact=True)), loggamma(x), exact=True)
#f17
gf17 = lambda x, y: sf.pochrel(x, y)
rf17 = lambda x, y: fsub(
fdiv(gamma(fadd(x, y, exact=True)), gamma(x)), 1, exact=True) / y
#f18
gf18 = lambda x, y: sf.gamma_inc(x, y)
rf18 = lambda x, y: gammainc(x, y)
# print gf18(-338.122,0.122)
# print rf18(-338.122,0.122)
#f19
gf19 = lambda x, y: sf.gamma_inc_Q(x, y)
rf19 = lambda x, y: gammainc(x, y, regularized=True)
#f20
gf20 = lambda x, y: sf.gamma_inc_P(x, y)
rf20 = lambda x, y: gammainc(x, 0, y, regularized=True)
#f21
gf21 = lambda x, y: sf.beta(x, y)
rf21 = lambda x, y: beta(x, y)
def _integral(a, b, index, ec):
ap1 = index + 1
integrand = lambda x: -pow(ec, ap1) * gamma_inc(ap1, x / ec)
return integrand(b) - integrand(a)
def integrand(x): return -pow(ec, ap1) * \
gamma_inc(ap1, old_div(x, ec))
return integrand(b) - integrand(a)