scipy.special.bdtr

scipy.special.bdtr(x1, x2, x3[, out])

y=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)