%I #19 Oct 27 2023 10:41:48
%S 4,2,1,2,7,9,5,4,3,9,8,3,9,0,3,4,3,2,7,6,8,8,2,1,7,6,0,6,5,0,2,9,8,0,
%T 9,1,6,1,0,3,6,7,2,1,4,0,7,2,6,1,2,2,3,2,1,6,5,4,3,7,5,4,5,4,0,6,5,1,
%U 7,2,9,3,9,2,2,4,3,7,7,9,1,5,3,6,3,2,9,0,6,8,8,4,7,1,9,2,4,6,2,4,3,9
%N Decimal expansion of the greatest minimal separation between ten points in a unit square.
%C The corresponding values for two to nine points have simple expressions:
%C N ... d_min
%C 2 ... sqrt(2) (A002193)
%C 3 ... sqrt(6) - sqrt(2) (A120683)
%C 4 ... 1 (A000007)
%C 5 ... sqrt(2) / 2 (A010503)
%C 6 ... sqrt(13) / 6
%C 7 ... 4 - 2*sqrt(3)
%C 8 ... sqrt(2 - sqrt(3)) (A101263)
%C 9 ... 1 / 2 (A020761)
%C In contrast, the value for ten points has a minimal polynomial of degree 18.
%C The smallest square ten unit circles will fit into has side length s = 2 + 2/d = 6.74744152... and the maximum radius of ten non-overlapping circles in the unit square is 1 / s = 0.14820432...
%H C. de Groot, R. Peikert, D. Würtz, <a href="https://www.researchgate.net/publication/228327432_The_Optimal_Packing_of_Ten_Equal_Circles_in_a_Square">The Optimal Packing of Ten Equal Circles in a Square</a>, IPS Research Report, ETH Zürich, No. 90-12, August 1990.
%H Eckard Specht, <a href="http://hydra.nat.uni-magdeburg.de/packing/csq/">The best known packings of equal circles in a square</a>
%H Jeremy Tan, <a href="https://gist.github.com/Parclytaxel/fe51678ea07cc448c56c3927afc44ac1">Sympy (Python) program</a>
%H <a href="/index/Al#algebraic_18">Index entries for algebraic numbers, degree 18</a>
%F d is the smallest real root of 1180129*d^18 - 11436428*d^17 + 98015844*d^16 - 462103584*d^15 + 1145811528*d^14 - 1398966480*d^13 + 227573920*d^12 + 1526909568*d^11 - 1038261808*d^10 - 2960321792*d^9 + 7803109440*d^8 - 9722063488*d^7 + 7918461504*d^6 - 4564076288*d^5 + 1899131648*d^4 - 563649536*d^3 + 114038784*d^2 - 14172160*d + 819200.
%e 0.421279543983903432768821760650298...
%o (PARI) p = Pol([1180129, -11436428, 98015844, -462103584, 1145811528, -1398966480, 227573920, 1526909568, -1038261808, -2960321792, 7803109440, -9722063488, 7918461504, -4564076288, 1899131648, -563649536, 114038784, -14172160, 819200]); polrootsreal(p)[1]
%Y Cf. A281115 (10 points in unit circle).
%K nonn,cons
%O 0,1
%A _Jeremy Tan_, Jan 14 2017
|