def test_incomplete_gamma_sum(self):
"""Test if the sum of lower and upper incomplete gamma functions correctly produce the
gamma function"""
rand = Random(372924)
for _ in range(10):
a, x = rand.random() * 10, rand.random() * 10
l = incomplete_gamma_lower(a, x)
r = incomplete_gamma_upper(a, x)
g = gamma(a)
self.assertAlmostEqual(l + r, g, 10)
评论列表
文章目录