scipy.special.bdtrc

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

y=bdtrc(k,n,p) returns the sum of the terms k+1 through n of the Binomial probability density: sum(nCj p**j (1-p)**(n-j), j=k+1..n)