I got tripped up by this recently when doing eigenvalue calculations in python. I wanted to evaluate the incomplete gamma function . After using the SciPy function I was scratching my head as to why I was seeing a discrepancy between my numerical calculations for the eigenvalues and my theoretical calculation. Then I came across this post by John D Cook that helped clarify things. The SciPy function actually calculates .