How to generate zero-truncated negative binomial distribution?

we know that pk(truncated)=pk/(1-p0) so Fk(truncated)=p0(truncated)+p1(truncated)+p2(truncated)+...+pk(truncated) = 0 + [p1/(1-p0)] + [p2/(1-p0)] + ... + [pk/(1-p0)] = 0 + (p1+p2+...+pk)/(1-p0) Fk(truncated)= 0 + (p1+p2+...+pk+p0-p0)/(1-p0) = 0 + (Fk-p0)/(1-p0) we get [(1-p0)*Fk(truncated)]+p0 = Fk y=rand(100) so can we use nbininv([(1-p0)*y+p0],R,P)??