Discrete probability distribution calculation in Matlab...

Code returns minimal m s.t P(X>m)=1/n^2 for binomial experiment...

How can I efficiently calculate the negative binomial cumulative distribution function?...