The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A038534 Numerators of coefficients of EllipticK/Pi. 25
 1, 1, 9, 25, 1225, 3969, 53361, 184041, 41409225, 147744025, 2133423721, 7775536041, 457028729521, 1690195005625, 25145962430625, 93990019574025, 90324408810638025, 340357374376418025, 5147380044581630625, 19520119892056100625, 1187604094232693162025 (list; graph; refs; listen; history; text; internal format)
 OFFSET 0,3 COMMENTS The denominators are given in A038533. Also numerators in expansion of the hypergeometric series 2F1(1/2,1/2; 1; x). This means numerators of the expansion coefficients of 2*K(k)/Pi = 2F1(1/2,1/2; 1; k^2) in powers of k^2, with K(k) the complete elliptic integral of the first kind. The denominators are given in A056982. The period T of the plane pendulum (mass m, length L, Earth's gravity g, energy E) is 4*sqrt(L/g)*K(sin(phi_0/2)) with cos(phi_0) = -E/(m*g*L) (maximal phi value). See the Landau - Lifschitz reference, p. 30. - Wolfdieter Lang, May 29 2016 It is easy and inexpensive to make a satisfactory precision measurement of a(1)/4, a(2)/64, and a(3)/256 using a pendulum rigged from a computer mouse. In "Digital Pendulum Data Analysis" (see links) amplitude vs. time data is transformed to period vs. sin(phi_0/2)^2 data, thus allowing extraction of expansion coefficients as fit parameters.- Bradley Klee, Dec 25 2016 REFERENCES B. C. Berndt, Ramanujan's Notebooks Part III, Springer-Verlag, see p. 91, Eq. 2.1. L. D. Landau und E. M. Lifschitz, Mechanik, Akademie Verlag, Berlin, 1967, p. 30 (Exercise 1 in chapter III, paragraph 11.) LINKS Vincenzo Librandi, Table of n, a(n) for n = 0..200 B. Klee, Digital Pendulum Data Analysis: Output, Github, 2016. G. N. Watson, A Note on Gamma Functions.Edinburgh Mathematical Notes, 42, 1959, pp 7-9. FORMULA a(n) = 2^(-2*w(n))*binomial(2*n,n)^2 with w(n) = A000120(n), the number of 1's in binary expansion of n. a(n) = A001790(n)^2. a(n) = (A056040(2*n)/A060632(2*n))^2. - Peter Luschny, Aug 06 2014 a(n) = (-1)^n*A056982(n)*C(-1/2,n)*C(n-1/2,n). - Peter Luschny, Apr 08 2016 MAPLE swing := proc(n) option remember; if n = 0 then 1 elif n mod 2 = 1 then swing(n-1)*n else 4*swing(n-1)/n fi end: sigma := n -> 2^(add(i, i = convert(iquo(n, 2), base, 2))): a := n -> (swing(2*n)/sigma(2*n))^2; seq(a(n), n=0..20); # Peter Luschny, Aug 06 2014 MATHEMATICA Numerator@ CoefficientList[ Series[ EllipticK@x, {x, 0, 19}]/Pi, x] (* Robert G. Wilson v, Jul 19 2007 *) CROSSREFS Cf. A001790, A038533, A038535, A056040, A056982, A060632. Sequence in context: A201807 A228465 A048056 * A117755 A144114 A209969 Adjacent sequences:  A038531 A038532 A038533 * A038535 A038536 A038537 KEYWORD frac,nonn AUTHOR Wouter Meeussen, revised Jan 03 2001 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified July 25 16:03 EDT 2021. Contains 346291 sequences. (Running on oeis4.)