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

 

Logo

Please make a donation to keep the OEIS running. We are now in our 56th year. In the past year we added 10000 new sequences and reached almost 9000 citations (which often say "discovered thanks to the OEIS").
Other ways to donate

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A335828 Numerators of coefficients in a power series expansion of the distance between two bodies falling freely towards each other along a straight line under the influence of their mutual gravitational attraction. 2
1, 1, 11, 73, 887, 136883, 7680089, 26838347, 14893630313, 1908777537383, 2422889987331397, 233104477447558811, 2430782624763507659, 14420190617640617313953, 4515429325405165295004389, 812454316441781379614873497, 166481868581561511154267399013 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,3

COMMENTS

Consider two point objects with masses m_1 and m_2 that are starting to fall towards each other from rest at time t = 0 and initial distance r_0. Foong (2008) gave the solution for the distance as a function of time, r(t) = r_0 * f(t/t_0), where t_0 = sqrt(r_0^3/(G*(m1+m2)), G is the gravitational constant (A070058), and f(x) = 1 - Sum_{n>=1} c(n) * x^(2*n) is a dimensionless function. c(n) are the rational coefficients whose numerators are given in this sequence. The denominators are given in A335829. The collision occurs when f(x) = 0, at x = Pi/(2*sqrt(2)) (A093954), which corresponds to the time t = (Pi/(2*sqrt(2))) * t_0.

A similar expansion was given by Ernst Meissel in his study of the three-body problem in 1882. In Meissel's expansion the coefficients are c(n)/2^n.

REFERENCES

Sudhir Ranjan Jain, Mechanics, Waves and Thermodynamics: An Example-based Approach, Cambridge University Press, 2016. See page 97.

Ernst Meissel, Über Reihen, denen man bei der numerischen Lösung des Problems der Dreikörperproblems begegnet, wenn die Anfangsgeschwindigkeiten Null sind, in: Jahresbericht über die Realschule in Kiel: Während des Schuljahres 1881/82, A. F. Jensen, Kiel, 1882, pp. 1-11.

LINKS

Amiram Eldar, Table of n, a(n) for n = 1..243

S. K. Foong, From Moon-fall to motions under inverse square laws, European journal of physics, Vol. 29, No. 5 (2008), pp. 987-1003, alternative link.

Jaak Peetre, Ernst Meissel and the Pythagorean problem - the Drei-Körper-Problem in the Nachlass Meissel, draft, 1997.

FORMULA

a(n) = numerator(c(n)), c(1) = 1/2, c(n) = (2 * Sum_{k=1..n-1} (n-k)*(2*n-2*k-1)*c(n-k)*c(k) - Sum_{m=2..n-1} (n-m)*(2*n-2*m-1)*c(n-m) * Sum_{k=1..m-1} c(m-k)*c(k))/(n*(2*n - 1)).

c(n) ~ c_0 * n^(-5/3) * (Pi/(2*sqrt(2))^(-2*n), where c_0 = (3*Pi)^(2/3) / (18*Gamma(4/3)) = 0.277587...

EXAMPLE

The series begins with f(x) = 1 - (1/2)*x^2 - (1/12)*x^4 - (11/360)*x^6 - ...

MATHEMATICA

c[1] = 1/2; c[n_] := c[n] = (2*Sum[(n - k)*(2*n - 2*k - 1)*c[n - k]*c[k], {k, 1, n - 1}] - Sum[(n - m)*(2*n - 2*m - 1)*c[n - m]*c[m - k]*c[k], {m, 2, n - 1}, {k, 1, m - 1}])/(n*(2*n - 1)); Numerator @ Array[c, 17]

(* or *)

Quiet[-Numerator @ CoefficientList[AsymptoticDSolveValue[{y[x]*y'[x]^2 == 2*(1-y[x]), y[0] == 1}, y[x], {x, 0, 25}], x][[3;; -1;; 2]]] (* requires Mathematica 11.3+ *)

CROSSREFS

Cf. A070058, A093954, A202623, A335829 (denominators).

Sequence in context: A163775 A092244 A155634 * A003367 A121784 A213165

Adjacent sequences:  A335825 A335826 A335827 * A335829 A335830 A335831

KEYWORD

nonn,frac

AUTHOR

Amiram Eldar, Jun 25 2020

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.

License Agreements, Terms of Use, Privacy Policy. .

Last modified December 5 04:13 EST 2020. Contains 338943 sequences. (Running on oeis4.)