y=nbdtrc(k,n,p) returns the sum of the terms k+1 to infinity of the negative binomial distribution.
scipy.special.nbdtr
scipy.special.nbdtri