cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 26 results. Next

A048056 Duplicate of A038534.

Original entry on oeis.org

1, 1, 9, 25, 1225, 3969, 53361, 184041, 41409225, 147744025, 2133423721
Offset: 0

Views

Author

Keywords

A056982 a(n) = 4^A005187(n). The denominators of the Landau constants.

Original entry on oeis.org

1, 4, 64, 256, 16384, 65536, 1048576, 4194304, 1073741824, 4294967296, 68719476736, 274877906944, 17592186044416, 70368744177664, 1125899906842624, 4503599627370496, 4611686018427387904, 18446744073709551616, 295147905179352825856, 1180591620717411303424
Offset: 0

Views

Author

Keywords

Comments

Also equal to A046161(n)^2.
Let W(n) = Product_{k=1..n} (1- 1/(4*k^2)), the partial Wallis product with lim n -> infinity W(n) = 2/Pi; a(n) = denominator(W(n)). The numerators are in A069955.
Equivalently, denominators in partial products of the following approximation to Pi: Pi = Product_{n >= 1} 4*n^2/(4*n^2-1). Numerators are in A069955.
Denominator of h^(2n) in the Kummer-Gauss series for the perimeter of an ellipse.
Denominators of coefficients in hypergeometric([1/2,-1/2],[1],x). The numerators are given in A038535. hypergeom([1/2,-1/2],[1],e^2) = L/(2*Pi*a) with the perimeter L of an ellipse with major axis a and numerical eccentricity e (Maclaurin 1742). - Wolfdieter Lang, Nov 08 2010
Also denominators of coefficients in hypergeometric([1/2,1/2],[1],x). The numerators are given in A038534. - Wolfdieter Lang, May 29 2016
Also denominators of A277233. - Wolfdieter Lang, Nov 16 2016
A277233(n)/a(n) are the Landau constants. These constants are defined as G(n) = Sum_{j=0..n} g(j)^2, where g(n) = (2*n)!/(2^n*n!)^2 = A001790(n)/A046161(n). - Peter Luschny, Sep 27 2019

References

  • J.-P. Delahaye, Pi - die Story (German translation), Birkhäuser, 1999 Basel, p. 84. French original: Le fascinant nombre Pi, Pour la Science, Paris, 1997.
  • O. J. Farrell and B. Ross, Solved Problems in Analysis, Dover, NY, 1971; p. 77.

Crossrefs

Apart from offset, identical to A110258.
Equals (1/2)*A038533(n), A038534, A277233.

Programs

  • Maple
    A056982 := n -> denom(binomial(1/2, n))^2:
    seq(A056982(n), n=0..19); # Peter Luschny, Apr 08 2016
    # Alternatively:
    G := proc(x) hypergeom([1/2,1/2], [1], x)/(1-x) end: ser := series(G(x), x, 20):
    [seq(coeff(ser,x,n), n=0..19)]: denom(%); # Peter Luschny, Sep 28 2019
  • Mathematica
    Table[Power[4, 2 n - DigitCount[2 n, 2, 1]], {n, 0, 19}] (* Michael De Vlieger, May 30 2016, after Harvey P. Dale at A005187 *)
    G[x_] := (2 EllipticK[x])/(Pi (1 - x));
    CoefficientList[Series[G[x], {x, 0, 19}], x] // Denominator (* Peter Luschny, Sep 28 2019 *)
  • PARI
    a(n)=my(s=n); while(n>>=1, s+=n); 4^s \\ Charles R Greathouse IV, Apr 07 2012

Formula

a(n) = (denominator(binomial(1/2, n)))^2. - Peter Luschny, Sep 27 2019

Extensions

Edited by N. J. A. Sloane, Feb 18 2004, Jun 05 2007

A038533 Denominator of coefficients of both EllipticK/Pi and EllipticE/Pi.

Original entry on oeis.org

2, 8, 128, 512, 32768, 131072, 2097152, 8388608, 2147483648, 8589934592, 137438953472, 549755813888, 35184372088832, 140737488355328, 2251799813685248, 9007199254740992, 9223372036854775808, 36893488147419103232, 590295810358705651712, 2361183241434822606848
Offset: 0

Views

Author

Wouter Meeussen, revised Jan 03 2001

Keywords

Comments

Denominators are powers of 2 since EllipticK(x) = Pi * Sum_{n>=0} 2^(-4*n-1) * binomial(2*n,n)^2 * x^n and EllipticE(x) = Pi * Sum_{n>=0} 2^(-4*n-1) (-1)^(2*n) * binomial(2*n,n)^2 /(-2*n+1) * x^n.

Crossrefs

Equals 2*A056982(n).

Programs

  • Mathematica
    a[n_] := 2^(4*n - 2*DigitCount[n, 2, 1] + 1); Array[a, 20, 0] (* Amiram Eldar, Aug 03 2023 *)
  • PARI
    a(n)=my(s=n); while(n>>=1, s+=n); 2<<(2*s) \\ Charles R Greathouse IV, Apr 07 2012

Formula

a(n) = 2^(1+4*n-2*w(n)) with w(n) = A000120(n) = number of 1's in binary expansion of n.

A273506 T(n,m), numerators of coefficients in a power/Fourier series expansion of the plane pendulum's exact phase space trajectory.

Original entry on oeis.org

1, -1, 7, 1, -1, 11, -1, 319, -143, 715, 1, -26, 559, -221, 4199, -2, 139, -323, 6137, -2261, 52003, 1, -10897, 135983, -4199, 527459, -52003, 37145, -1, 15409, -317281, 21586489, -52877, 7429, -88711, 1964315, 1, -76, 269123, -100901, 274873, -8671, 227447, -227447, 39803225, -2, 466003, -213739, 522629, -59074189, 226061641, -10690009, 25701511, -42077695, 547010035
Offset: 1

Views

Author

Bradley Klee, May 23 2016

Keywords

Comments

Triangle read by rows ( see examples ). The phase space trajectory of a simple pendulum can be written as (q,p) = (R(Q)cos(Q),R(Q)sin(Q)), with scaled, canonical coordinates q and p. The present triangle and A273507 determine a power / Fourier series of R(Q): R(Q) = sqrt(4 *k) * (1 + sum k^n * (A273506(n,m)/A273507(n,m)) * cos(Q)^(2(n+m)) ); where the sum runs over n = 1,2,3 ... and m = 1,2,3...n. The period of an oscillator can be computed by T(k) = dA/dE, where A is the phase area enclosed by the phase space trajectory of conserved, total energy E. As we choose expansion parameter "k" proportional to E, the series expansion of the complete elliptic integral of the first kind follows from T(k) with very little technical difficulty ( see examples and Mathematica function R2ToEllK ). For more details read "Plane Pendulum and Beyond by Phase Space Geometry" (Klee, 2016).
For some remarks on this pendulum problem and an alternative way to compute a(n,m) / A273507(n,m) using Lagrange inversion see the two W. Lang links. - Wolfdieter Lang, Jun 11 2016

Examples

			n/m  1    2     3     4
------------------------------
1  |  1
2  | -1,  7
3  |  1, -1,    11
4  | -1,  319, -143, 715
------------------------------
R2(Q) = sqrt(4 k) (1 + (1/6) cos(Q)^4 k +  (-(1/45) cos(Q)^6 + (7/72) cos(Q)^8) k^2)
R2(Q)^2 = 4 k + (4/3) cos(Q)^4 k^2 + ( -(8/45) cos(Q)^6 + (8/9) cos(Q)^8)k^3 + ...
I2 = (1/(2 Pi)) Int dQ (1/2)R2(Q)^2 = 2 k + (1/4) k^2 + (3/32) k^3 + ...
(2/Pi) K(k) ~ (1/2)d/dk(I2) = 1 + (1/4) k + (9/64) k^2 + ...
From _Wolfdieter Lang_, Jun 11 2016 (Start):
The rational triangle r(n,m) = a(n, m) / A273507(n,m) begins:
n\m   1          2          3         4   ...
1:   1/6
2: -1/45        7/72
3:  1/630      -1/30      11/144
4: -1/14175   319/56700 -143/3240  715/10368
... ,
row n = 5: 1/467775 -26/42525 559/45360 -221/3888 4199/62208,
row 6: -2/42567525 139/2910600 -323/145800 6137/272160 -2261/31104 52003/746496,
row 7: 1/1277025750 -10897/3831077250 135983/471517200 -4199/729000 527459/13996800 -52003/559872 37145/497664,
row 8:
-1/97692469875 15409/114932317500 -317281/10945935000 21586489/20207880000 -52877/4199040 7429/124416 -88711/746496 1964315/23887872.
... (End)
		

Crossrefs

Denominators: A273507. Time Dependence: A274076, A274078, A274130, A274131. Elliptic K: A038534, A056982. Cf. A000984, A001790, A038533, A046161, A273496.

Programs

  • Mathematica
    R[n_] := Sqrt[4 k] Plus[1, Total[k^# R[#, Q] & /@ Range[n]]]
    Vq[n_] :=  Total[(-1)^(# - 1) (r Cos[Q] )^(2 #)/((2 #)!) & /@ Range[2, n]]
    RRules[n_] :=  With[{H = ReplaceAll[1/2 r^2 + (Vq[n + 1]), {r -> R[n]}]},
    Function[{rules}, Nest[Rule[#[[1]], ReplaceAll[#[[2]], rules]] & /@ # &, rules, n]][
       Flatten[R[#, Q] ->  Expand[(-1/4) ReplaceAll[ Coefficient[H, k^(# + 1)], {R[#, Q] -> 0}]] & /@ Range[n]]]]
    RCoefficients[n_] :=  With[{Rn = ReplaceAll[R[n], RRules[n]]}, Function[{a},
        Coefficient[Coefficient[Rn/2/Sqrt[k], k^a],
           Cos[Q]^(2 (a + #))] & /@ Range[a]] /@ Range[n]]
    R2ToEllK[NMax_] := D[Expand[(2)^(-2) ReplaceAll[R[NMax], RRules[NMax]]^2] /. {Cos[Q]^n_ :> Divide[Binomial[n, n/2], (2^(n))], k^n_ /; n > NMax -> 0},k]
    Flatten[Numerator@RCoefficients[10]]
    R2ToEllK[10]

A273496 Triangle read by rows: coefficients in the expansion cos(x)^n = (1/2)^n * Sum_{k=0..n} T(n,k) * cos(k*x).

Original entry on oeis.org

1, 0, 2, 2, 0, 2, 0, 6, 0, 2, 6, 0, 8, 0, 2, 0, 20, 0, 10, 0, 2, 20, 0, 30, 0, 12, 0, 2, 0, 70, 0, 42, 0, 14, 0, 2, 70, 0, 112, 0, 56, 0, 16, 0, 2, 0, 252, 0, 168, 0, 72, 0, 18, 0, 2, 252, 0, 420, 0, 240, 0, 90, 0, 20, 0, 2
Offset: 0

Views

Author

Bradley Klee, May 23 2016

Keywords

Comments

These coefficients are especially useful when integrating powers of cosine x (see examples).
Nonzero, even elements of the first column are given by A000984; T(2n,0) = binomial(2n,n).
For the rational triangles for even and odd powers of cos(x) see A273167/A273168 and A244420/A244421, respectively. - Wolfdieter Lang, Jun 13 2016
Mathematica needs no TrigReduce to integrate Cos[x]^k. See link. - Zak Seidov, Jun 13 2016

Examples

			n/k|  0   1   2   3   4   5   6
-------------------------------
0  |  1
1  |  0   2
2  |  2   0   2
3  |  0   6   0   2
4  |  6   0   8   0   2
5  |  0   20  0   10  0   2
6  |  20  0   30  0   12  0   2
-------------------------------
cos(x)^4 = (1/2)^4 (6 + 8 cos(2x) + 2 cos(4x)).
I4 = Int dx cos(x)^4 = (1/2)^4 Int dx ( 6 + 8 cos(2x) + 2 cos(4x) ) = C + 3/8 x + 1/4 sin(2x) + 1/32 sin(4x).
Over range [0,2Pi], I4 = (3/4) Pi.
		

Crossrefs

Programs

  • Mathematica
    T[MaxN_] := Function[{n}, With[
           {exp = Expand[Times[ 2^n, TrigReduce[Cos[x]^n]]]},
           Prepend[Coefficient[exp, Cos[# x]] & /@ Range[1, n],
            exp /. {Cos[_] -> 0}]]][#] & /@ Range[0, MaxN];Flatten@T[10]
    (* alternate program *)
    T2[MaxN_] := Function[{n}, With[{exp = Expand[(Exp[I x] + Exp[-I x])^n]}, Prepend[2 Coefficient[exp, Exp[I # x]] & /@ Range[1, n], exp /. {Exp[] -> 0}]]][#] & /@ Range[0, MaxN]; T2[10] // ColumnForm (* _Bradley Klee, Jun 13 2016 *)

Formula

From Robert Israel, May 24 2016: (Start)
T(n,k) = 0 if n-k is odd.
T(n,0) = binomial(n,n/2) if n is even.
T(n,k) = 2*binomial(n,(n-k)/2) otherwise. (End)

A273507 T(n, m), denominators of coefficients in a power/Fourier series expansion of the plane pendulum's exact phase space trajectory.

Original entry on oeis.org

6, 45, 72, 630, 30, 144, 14175, 56700, 3240, 10368, 467775, 42525, 45360, 3888, 62208, 42567525, 2910600, 145800, 272160, 31104, 746496, 1277025750, 3831077250, 471517200, 729000, 13996800, 559872, 497664, 97692469875, 114932317500, 10945935000, 20207880000, 4199040, 124416, 746496, 23887872
Offset: 1

Views

Author

Bradley Klee, May 23 2016

Keywords

Comments

Triangle read by rows ( see example ). The numerator triangle is A274076.
Comments of A273506 give a definition of the fraction triangle, which determines an arbitrary-precision solution to the simple pendulum equations of motion. For more details see "Plane Pendulum and Beyond by Phase Space Geometry" (Klee, 2016).

Examples

			n/m  1      2      3     4
------------------------------
1  | 6
2  | 45,    72
3  | 630,   30,    144
4  | 14175, 56700, 3240, 10368
------------------------------
		

Crossrefs

Numerators: A273506. Time Dependence: A274076, A274078, A274130, A274131. Elliptic K: A038534, A056982. Cf. A000984, A001790, A038533, A046161, A273496.

Programs

  • Mathematica
    R[n_] := Sqrt[4 k] Plus[1, Total[k^# R[#, Q] & /@ Range[n]]]
    Vq[n_] :=  Total[(-1)^(# - 1) (r Cos[Q] )^(2 #)/((2 #)!) & /@ Range[2, n]]
    RRules[n_] :=  With[{H = ReplaceAll[1/2 r^2 + (Vq[n + 1]), {r -> R[n]}]},
    Function[{rules}, Nest[Rule[#[[1]], ReplaceAll[#[[2]], rules]] & /@ # &, rules, n]][
       Flatten[R[#, Q] ->  Expand[(-1/4) ReplaceAll[
              Coefficient[H, k^(# + 1)], {R[#, Q] -> 0}]] & /@ Range[n]]]]
    RCoefficients[n_] :=  With[{Rn = ReplaceAll[R[n], RRules[n]]}, Function[{a},
        Coefficient[Coefficient[Rn/2/Sqrt[k], k^a],
           Cos[Q]^(2 (a + #))] & /@ Range[a]] /@ Range[n]]
    Flatten[Denominator@RCoefficients[10]]

A274076 T(n, m), numerators of coefficients in a power/Fourier series expansion of the plane pendulum's exact differential time dependence.

Original entry on oeis.org

-2, 2, -2, -4, 8, -20, 2, -58, 14, -70, -4, 16, -344, 112, -28, 4, -556, 1064, -152, 308, -308, -8, 10256, -3368, 4576, -6248, 2288, -1144, 2, -1622, 33398, -98794, 34606, -4862, 2002, -1430, -4, 6688, -187216, 140384, -1242904, 59488, -25168, 77792, -48620
Offset: 1

Views

Author

Bradley Klee, Jun 09 2016

Keywords

Comments

Triangle read by rows ( see examples ). The denominators are given in A274078.
The rational triangle A273506 / A273507 gives the coefficients for an exact solution of the plane pendulum's phase space trajectory. Differential time dependence for this solution also adheres to the simple form of a triangular summation: dt = dQ(-1+ sum k^n * (T(n, m)/A274078(n, m)) * cos(Q)^(2(n+m)) ); where the sum runs over n = 1,2,3 ... and m = 1,2,3...n. Expanding powers of cosine ( Cf. A273496 ) it is relatively easy to integrate dt ( cf. A274130 ). One period of motion takes Q through the range [ 0 , -2 pi]. Integrating dt over this domain gives another (Cf. A273506) calculation of the series expansion for Elliptic K ( see examples and Mathematica function dtToEllK ). For more details read "Plane Pendulum and Beyond by Phase Space Geometry" (Klee, 2016).

Examples

			The triangle T(n, m) begins:
n/m  1    2     3     4
------------------------------
1  | -2
2  |  2, -2
3  | -4,  8,  -20
4  |  2, -58,  14,  -70
------------------------------
The rational triangle T(n, m) / A274078(n, m) begins:
n/m    1        2         3       4
------------------------------------------
1  | -2/3
2  |  2/15,   -2/3
3  | -4/315,   8/27,   -20/27
4  |  2/2835, -58/945,  14/27,  -70/81
------------------------------------------
dt2(Q) = dQ(-1 - (2/3) cos(Q)^4 k +  ((2/15) cos(Q)^6  - (2/3) cos(Q)^8) k^2 ) + ...
dt2(Q) = dQ(-1 - (1/4) k - (9/64) k^2 + cosine series ) + ...
(2/Pi) K(k) ~ I2 = (1/(2 Pi)) Int dt2(Q) =  1 + (1/4) k + (9/64) k^2+ ...
		

Crossrefs

Denominators: A274078. Phase Space Trajectory: A273506, A273507. Time Dependence: A274130, A274131. Elliptic K: A038534, A056982. Cf. A000984, A001790, A038533, A046161, A273496.

Programs

  • Mathematica
    R[n_] := Sqrt[4 k] Plus[1, Total[k^# R[#, Q] & /@ Range[n]]]
    Vq[n_] :=  Total[(-1)^(# - 1) (r Cos[Q] )^(2 #)/((2 #)!) & /@ Range[2, n]]
    RRules[n_] :=  With[{H = ReplaceAll[1/2 r^2 + (Vq[n + 1]), {r -> R[n]}]},
    Function[{rules}, Nest[Rule[#[[1]], ReplaceAll[#[[2]], rules]] & /@ # &, rules, n]][
       Flatten[R[#, Q] ->  Expand[(-1/4) ReplaceAll[ Coefficient[H, k^(# + 1)], {R[#, Q] -> 0}]] & /@ Range[n]]]]
    dt[n_] := With[{rules = RRules[n]}, Expand[Subtract[ Times[Expand[D[R[n] /. rules, Q]], Normal@Series[1/R[n], {k, 0, n}] /. rules, Cot[Q] ], 1]]]
    dtCoefficients[n_] :=  With[{dtn = dt[n]}, Function[{a}, Coefficient[ Coefficient[dtn, k^a], Cos[Q]^(2 (a + #))] & /@ Range[a]] /@ Range[n]]
    dtToEllK[NMax_] := ReplaceAll[-dt[NMax], {Cos[Q]^n_ :> Divide[Binomial[n, n/2], (2^(n))], k^n_ /; n > NMax -> 0} ]
    Flatten[Numerator[dtCoefficients[10]]]
    dtToEllK[5]

A274078 T(n,m), denominators of coefficients in a power/Fourier series expansion of the plane pendulum's exact differential time dependence.

Original entry on oeis.org

3, 15, 3, 315, 27, 27, 2835, 945, 27, 81, 155925, 2025, 2025, 135, 27, 6081075, 779625, 30375, 405, 243, 243, 638512875, 212837625, 654885, 42525, 8505, 1215, 729, 10854718875, 638512875, 58046625, 4465125, 127575, 3645, 729, 729
Offset: 1

Views

Author

Bradley Klee, Jun 09 2016

Keywords

Comments

Triangle read by rows (see example). Comments of A274076 give a definition of the fraction triangle, which determines to arbitrary precision the differential time dependence for the time-independent solution (cf. A273506, A273507) of the plane pendulum's equations of motion. For more details see "Plane Pendulum and Beyond by Phase Space Geometry" (Klee, 2016).

Examples

			n\m|    1    2    3    4
---+---------------------
1  |    3;
2  |   15,   3;
3  |  315,  27,  27;
4  | 2835, 945,  27,  81;
		

Crossrefs

Numerators: A274076. Phase Space Trajectory: A273506, A273507. Time Dependence: A274130, A274131. Elliptic K: A038534, A056982. Cf. A000984, A001790, A038533, A046161, A273496.

Programs

  • Mathematica
    R[n_] := Sqrt[4 k] Plus[1, Total[k^# R[#, Q] & /@ Range[n]]]
    Vq[n_] :=  Total[(-1)^(# - 1) (r Cos[Q] )^(2 #)/((2 #)!) & /@ Range[2, n]]
    RRules[n_] :=  With[{H = ReplaceAll[1/2 r^2 + (Vq[n + 1]), {r -> R[n]}]},
    Function[{rules}, Nest[Rule[#[[1]], ReplaceAll[#[[2]], rules]] & /@ # &, rules, n]][
       Flatten[R[#, Q] ->  Expand[(-1/4) ReplaceAll[ Coefficient[H, k^(# + 1)], {R[#, Q] -> 0}]] & /@ Range[n]]]]
    dt[n_] := With[{rules = RRules[n]}, Expand[Subtract[ Times[Expand[D[R[n] /. rules, Q]], Normal@Series[1/R[n], {k, 0, n}] /. rules, Cot[Q] ], 1]]]
    dtCoefficients[n_] :=  With[{dtn = dt[n]}, Function[{a}, Coefficient[ Coefficient[dtn, k^a], Cos[Q]^(2 (a + #))] & /@ Range[a]] /@ Range[n]]
    Flatten[Denominator[dtCoefficients[10]]]

A274130 Irregular triangle T(n,m), numerators of coefficients in a power/Fourier series expansion of the plane pendulum's exact time dependence.

Original entry on oeis.org

1, 1, 11, 29, 1, 1, 491, 863, 6571, 4399, 13, 5, 1568551, 28783, 45187, 312643, 4351, 1117, 17, 35, 25935757, 81123251, 2226193, 2440117, 16025, 34246631, 18161, 35443, 49, 7, 5301974777, 22870237, 1603483793, 23507881213, 122574691, 122330761339, 903325919, 1976751869, 956873, 18551, 35, 77
Offset: 1

Views

Author

Bradley Klee, Jun 10 2016

Keywords

Comments

Irregular triangle read by rows ( see examples ). The row length sequence is 2*n = A005843(n), n >= 1.The denominators are given in A274131.
The triangles A274076 and A274078 give the coefficients for the exact, differential time dependence of the plane pendulum's equations of motion. Integrating, we obtain time dependence as a Fourier sine series: t = -( (2/pi)K(k) Q + sum k^n * (T(n,m)/A274131(n,m)) * sin(2 m Q) ); where the sum runs over n = 1,2,3 ... and m = 1,2,3,...,2 n. Combining the phase space trajectory and time dependence, it is possible to express Jacobian elliptic functions {cn,dn} in parametric form. For more details read "Plane Pendulum and Beyond by Phase Space Geometry" (Klee, 2016).

Examples

			n\m  1     2     3      4    5   6 ...
-----------------------------------------
1  | 1    1
2  | 11   29    1      1
3  | 491  863   6571   4399  13  5
row n=4: 1568551, 28783, 45187, 312643, 4351, 1117, 17, 35,
row n=5: 25935757, 81123251, 2226193, 2440117, 16025, 34246631, 18161, 35443, 49, 7.
-----------------------------------------
The rational irregular triangle T(n, m) / A274131(n, m) begins:
n\m    1          2           3             4            5         6
-----------------------------------------------------------------------------
1  |  1/6,      1/48
2  |  11/96,    29/960,    1/160,          1/1536
3  |  491/5760, 863/30720, 6571/725760, 4399/1935360, 13/34560, 5/165888
row n=4: 1568551/23224320, 28783/1161216, 45187/4644864, 312643/92897280, 4351/4644864, 1117/5806080, 17/663552, 35/21233664,
row n=5: 25935757/464486400, 81123251/3715891200, 2226193/232243200, 2440117/619315200, 16025/11354112, 34246631/81749606400, 18161/185794560, 35443/2123366400, 49/26542080, 7/70778880.
-----------------------------------------------------------------------------
t1(Q) =-Q -(1/4)*k*Q -k*((1/6)*Sin[2*Q]+(1/48)*Sin[4*Q])+...
(2/Pi) K(k) ~ (1/(2 Pi)) t1(-2*Pi) =  1+(1/4)*k+...
		

Crossrefs

Denominators: A274131. Phase Space Trajectory: A273506, A273507. Time Dependence: A274076, A274078. Elliptic K: A038534, A056982. Cf. A000984, A001790, A038533, A046161, A273496.

Programs

  • Mathematica
    R[n_] := Sqrt[4 k] Plus[1, Total[k^# R[#, Q] & /@ Range[n]]]
    Vq[n_] :=  Total[(-1)^(# - 1) (r Cos[Q] )^(2 #)/((2 #)!) & /@ Range[2, n]]
    RRules[n_] :=  With[{H = ReplaceAll[1/2 r^2 + (Vq[n + 1]), {r -> R[n]}]},
    Function[{rules}, Nest[Rule[#[[1]], ReplaceAll[#[[2]], rules]] & /@ # &, rules, n]][
       Flatten[R[#, Q] ->  Expand[(-1/4) ReplaceAll[ Coefficient[H, k^(# + 1)], {R[#, Q] -> 0}]] & /@ Range[n]]]]
    dt[n_] := With[{rules = RRules[n]}, Expand[Subtract[ Times[Expand[D[R[n] /. rules, Q]], Normal@Series[1/R[n], {k, 0, n}] /. rules, Cot[Q] ], 1]]]
    t[n_] := Expand[ReplaceAll[Q TrigReduce[dt[n]], Cos[x_ Q] :> (1/x/Q) Sin[x Q]]]
    tCoefficients[n_] := With[{tn = t[n]},Function[{a}, Coefficient[Coefficient[tn, k^a], Sin[2 # Q] ] & /@ Range[2 a]] /@ Range[n]]
    tToEllK[NMax_]:= Expand[((t[NMax] /. Q -> -2 Pi)/2/Pi) /. k^n_ /; n > NMax -> 0]
    Flatten[Numerator[-tCoefficients[10]]]
    tToEllK[5]

A274131 Irregular triangle T(n,m), denominators of coefficients in a power/Fourier series expansion of the plane pendulum's exact time dependence.

Original entry on oeis.org

6, 48, 96, 960, 160, 1536, 5760, 30720, 725760, 1935360, 34560, 165888, 23224320, 1161216, 4644864, 92897280, 4644864, 5806080, 663552, 21233664, 464486400, 3715891200, 232243200, 619315200, 11354112, 81749606400, 185794560, 2123366400, 26542080, 70778880
Offset: 1

Views

Author

Bradley Klee, Jun 10 2016

Keywords

Comments

Irregular triangle read by rows (see example). The row length sequence is 2*n = A005843(n), n >= 1.
The numerator triangle is A274130.
Comments of A274130 give a definition of the fraction triangle, which determines to arbitrary precision the time dependence for the time-independent solution (cf. A273506, A273507) of the plane pendulum's equations of motion. For more details see "Plane Pendulum and Beyond by Phase Space Geometry" (Klee, 2016).

Examples

			n\m  1      2      3          4       5       6
------------------------------------------------------
1  | 6     48
2  | 96    960    160      1536
3  | 5760  30720  725760   1935360  34560   165888
------------------------------------------------------
row 4: 23224320, 1161216, 4644864, 92897280, 4644864, 5806080, 663552, 21233664,
row 5: 464486400, 3715891200, 232243200, 619315200, 11354112, 81749606400, 185794560, 2123366400, 26542080, 70778880.
		

Crossrefs

Numerators: A274130. Phase Space Trajectory: A273506, A273507. Time Dependence: A274076, A274078. Elliptic K: A038534, A056982. Cf. A000984, A001790, A038533, A046161, A273496.

Programs

  • Mathematica
    R[n_] := Sqrt[4 k] Plus[1, Total[k^# R[#, Q] & /@ Range[n]]]
    Vq[n_] :=  Total[(-1)^(# - 1) (r Cos[Q] )^(2 #)/((2 #)!) & /@ Range[2, n]]
    RRules[n_] :=  With[{H = ReplaceAll[1/2 r^2 + (Vq[n + 1]), {r -> R[n]}]},
    Function[{rules}, Nest[Rule[#[[1]], ReplaceAll[#[[2]], rules]] & /@ # &, rules, n]][
       Flatten[R[#, Q] ->  Expand[(-1/4) ReplaceAll[ Coefficient[H, k^(# + 1)], {R[#, Q] -> 0}]] & /@ Range[n]]]]
    dt[n_] := With[{rules = RRules[n]}, Expand[Subtract[ Times[Expand[D[R[n] /. rules, Q]], Normal@Series[1/R[n], {k, 0, n}] /. rules, Cot[Q] ], 1]]]
    t[n_] := Expand[ReplaceAll[Q TrigReduce[dt[n]], Cos[x_ Q] :> (1/x/Q) Sin[x Q]]]
    tCoefficients[n_] := With[{tn = t[n]},Function[{a}, Coefficient[Coefficient[tn, k^a], Sin[2 # Q] ] & /@ Range[2 a]] /@ Range[n]]
    Flatten[Denominator[-tCoefficients[10]]]
Showing 1-10 of 26 results. Next