s=1-(1-p*x)/((r+1-r*x)*q) * 1/(x^(n+1)), x=1 + q*p^r + (r+1)(q*p^r)^2, q=0.9, p=0.1, r=5, n=3600