10 Think Outside the Box

We found that a very similar problem could be easily solved and that it agrees with the given problem to within 10^(-21). Since the answer is about 10^(-7), this is enough. The variation is to consider the region with no right boundary and ask for the probability that a particle hits the left before hitting the top or bottom. This turns out to have the simple and exact answer: 4/πarctan(^(-5 π)). Doubling this, and doing some very rough estimation on the discrepancy between our variation and the given problem, leads to the slightly better formula 4/πarctan (^(-5 π)) (1 + 1/(1 + 4/πarctan (^(-10π)))), which gives the answer to the problem posed to 14 digits. The answer to 100 digits is:
0.00000003837587979251226103407133186204839100793005594072509569030022799173436606852743276500842845647269910.

So, given the mathematical work we did, one can easily get 14 digits as follows.

NumberForm[N[4/π ArcTan[^(-5 π)] (1 + 1/(1 + 4/π ArcTan[^(-10 π)])) ], 15]

3.83758797925126 × 10^-7

Or 13 digits from:

NumberForm[N[8/π ArcTan[^(-5 π)] ], 13]

3.837587979251 × 10^-7

Because the argument is small, one gets the same digits even if the arctangent is ignored. We were lucky, since our method would fail if the rectangle was closer to a square!

Bornemann showed that the exact answer to the given question is 4/πUnderoverscript[∑, k = 0, arg3] (-1)^k /(2k + 1) sech[5π (2k + 1)], which is very similar to the steady-state series in #8. One term of this series gives 14 digits. Moreover, Bornemann observed that this series equals 8/πUnderoverscript[∑, k = 0, arg3] (-1)^k arctan (e^(-(2k + 1)   5  π)) (the first term of which is our arctan approximation) and it is well known (Borwein & Borwein, Pi and the AGM, §3.2, exer. 12) that these series equal 2/πarccosλ ( /10)^(1/2), where λ is the modular λ elliptic function. This formula is fast, taking only a small fraction of a second to get 500 digits.

N[2/π ArcCos[ModularLambda[/10]^(1/2)], 20]

3.8375879792512261034*10^-7


Created by Mathematica  (June 27, 2004)