cp's OEIS Frontend

A310000 Decimal expansion of AGM(1, phi/2), where phi is the golden ratio (A001622).

%I A310000 #28 Dec 14 2024 04:00:27
%S A310000 9,0,1,9,7,9,3,3,8,1,1,4,3,4,3,1,2,3,3,9,7,2,7,1,5,3,6,5,8,7,7,9,8,6,
%T A310000 2,7,5,5,1,6,2,3,7,4,6,7,3,6,9,9,0,1,4,0,7,9,8,4,7,7,9,4,2,9,1,1,9,4,
%U A310000 1,4,2,6,2,6,2,0,5,7,7,2,7,5,4,1,8
%N A310000 Decimal expansion of AGM(1, phi/2), where phi is the golden ratio (A001622).
%C A310000 Related to the pendulum acceleration relation at 72 degrees. 2*Pi*sqrt(l/g)/AGM(1, phi/2) gives the period T of a mathematical pendulum with a maximum deflection angle of 72 degrees from the downward vertical. The length of the pendulum is l and g is the gravitational acceleration.
%F A310000 Equals AGM(1, cos(Pi/5)).
%e A310000 0.9019793381143431233972715365...
%t A310000 RealDigits[ArithmeticGeometricMean[1, GoldenRatio/2], 10, 100][[1]] (* _Amiram Eldar_, Aug 26 2019 *)
%o A310000 (Python)
%o A310000 import decimal
%o A310000 iters = int(input('Precision: '))
%o A310000 decimal.getcontext().prec = iters
%o A310000 D = decimal.Decimal
%o A310000 def agm(a, b):
%o A310000     for x in range(iters):
%o A310000         a, b = (a + b) / 2,(a * b).sqrt()
%o A310000     return a
%o A310000 print(agm(1, (D(5).sqrt()+1)/4))
%o A310000 (PARI) agm(1, cos(Pi/5)) \\ _Michel Marcus_, Apr 05 2020
%Y A310000 Cf. A001622, A309893, A053004, A014549, A068521.
%K A310000 nonn,cons
%O A310000 0,1
%A A310000 _Daniel Hoyt_, Aug 26 2019