SciPy

scipy.special.bdtr

scipy.special.bdtr(k, n, p) returns the sum of the terms 0 through k of the Binomial probability density: sum(nCj p**j (1-p)**(n-j), j=0..k) = <ufunc 'bdtr'>