コード例 #1
0
 def integrand(x):
     return -pow(ec, ap1) * gamma_inc(ap1, old_div(x, ec))
コード例 #2
0
        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)
コード例 #3
0
#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)
コード例 #4
0
ファイル: functions.py プロジェクト: giacomov/astromodels 
        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)
コード例 #5
0
            def integrand(x):                return -pow(ec, ap1) * \
  gamma_inc(ap1, old_div(x, ec))

            return integrand(b) - integrand(a)