All published worksheets from http://sagenb.org
n = 160 N = 212480000 r = 7000 d = 52 p = (1-(falling_factorial(N-r, n) / falling_factorial(N, n))^d) * 1.0 p
1/p - 1