OFFSET
1,3
COMMENTS
The problem of calculating this average was apparently first published by Hans Jörg Oser in 1975, at the suggestion of the American mathematician Charles Royal Johnson. - Amiram Eldar, Mar 26 2026
LINKS
Vangalur S. Alagar, The Distribution of the Distance between Random Points, Journal of Applied Probability, Vol. 13, No. 3 (1976), pp. 558-566. See p. 565, eq. (18).
Allan J. Gottlieb (ed.), Puzzle Corner, Technology Review, MIT , March/April 1976, p. 69.
H. J. Oser, Problem 75-12, An Average Distance, SIAM Review, Vol. 17, No. 3 (1975), p. 566; Solution, by J. D. Murchland and D. J. Daley, ibid., Vol. 18, No. 3 (1976), pp. 497-501.
Kenneth J. Travers and Kenneth G. Gray, The Monte Carlo Method: A Fresh Approach to Teaching Probabilistic Concepts, The Mathematics Teacher, Vol. 74, No. 5 (1981), pp. 327-334. See p. 332.
FORMULA
Equals (116 - 8*sqrt(2) - 20*sqrt(5) + 140*log(1/2 + sqrt(5)/2) - 40*log(1 + sqrt(2)) + 80*log(2 + sqrt(5)) + log(32) + 10*log(-1 + sqrt(5)) - 15*log(123 + 55*sqrt(5))) / 120.
EXAMPLE
1.0881382498612540260133800398822187615937033186594483235371...
MATHEMATICA
d[x1_, y1_, x2_, y2_] := Sqrt[ (x1-x2)^2+(y1-y2)^2 ]
Integrate[ d[x1, y1, x2, y2], {x1, -1, 0}, {x2, 0, 1}, {y1, 0, 1}, {y2, 0, 1} ]
(* Alternative: *)
RealDigits[(116 - 8 Sqrt[2] - 20 Sqrt[5] + 140 ArcCsch[2] - 40 ArcSinh[1] + 80 ArcSinh[2] + Log[32] + 10 Log[-1 + Sqrt[5]] - 15 Log[123 + 55 Sqrt[5]]) / 120, 10, 120][[1]]
PROG
(PARI) (116 - 8*sqrt(2) - 20*sqrt(5) + 140*log(1/2 + sqrt(5)/2) - 40*log(1 + sqrt(2)) + 80*log(2 + sqrt(5)) + log(32) + 10*log(-1 + sqrt(5)) - 15*log(123 + 55*sqrt(5))) / 120 \\ Amiram Eldar, Mar 26 2026
CROSSREFS
KEYWORD
nonn,cons
AUTHOR
Richard C. Schroeppel and Michael Kleber, Mar 05 2008
STATUS
approved
