A223068 A sequence related to the period T of a simple gravity pendulum for arbitrary amplitudes.
1, 16, 3072, 737280, 1321205760, 951268147200, 2009078326886400, 265928913086054400, 44931349155019751424000, 109991942731488351485952000, 668751011807449177034588160000, 2471703739640332158319837839360000
Offset: 0
Examples
T = 2*Pi*sqrt(L/g) * (1 + (1/16)*phi^2 + (11/3072)*phi^4 + (173/737280)*phi^6 + ... ).
Programs
-
Maple
nmax:=11: f := series(1/((Pi/4)*(1+cos(phi/2))/EllipticK((1-cos(phi/2))/(1+cos(phi/2)))), phi, 2*nmax+1): for n from 0 to nmax do a(n):= denom(coeff(f, phi, 2*n)) od: seq(a(n), n=0..nmax); # Johannes W. Meijer, Jan 05 2017
-
Mathematica
s = Series[EllipticK[Sin[t/2]^2 ], {t, 0, 50}]; CoefficientList[2*s, t^2] // Denominator (* Jean-François Alcover, Oct 07 2014 *)
Comments