Ex = 1.09
p = 0.00109

r = Ex*p
# print(r)    # r=0.0011881

# 由于r需为正整数，取r=1
r = 1
res = 1

# P(X>=10) = 1-P(X=0)-P(X=1)-...-P(X=9)
p0_9 = pnbinom(9,r,p)

print(res-p0_9)