Solving Universal law of Gravitation as differential equation

  • Thread starter Treborman
  • Start date
  • #1
6
0
Hi, I'm trying to solve Newton Universal law of Gravitation as a differential equation in one dimension with two objects. One mass is fixed at distance r = 0, whilst the other is some initial distance away. I'm using distance rather than displacement and also ignored the proportionality factor to make the formula simpler. So in essence, this is the differential equation I want to solve:

[tex] a = \frac{d^2r}{dt^2} = -\frac{1}{r^2} [/tex]

However, I've hit a snag when trying to integrate for the 2nd time (see the image file). I have used two substitutions (a and b), and the last line in the file suggests that I will obtain the inverse function that will contain hyperbolic and polynomial terms, i.e. I won't be able to determine the original function of r with respect to t.

Can anyone help?
 

Attachments

Last edited:

Answers and Replies

  • #2
6
0
I've just discovered the Lagrange inversion theorem - would anybody be able to give me resources for this and limits if possible (preferably free, i.e. the internet).
 
  • #3
150
0
Mult. both sides by dr/dt: Integrate.
 
  • #4
6
0
I don't get it - do you mean integrate with respect to dr/dt? Would you be able to get me started - I was (hoping) to expand this investigation fairly quickly, not get stuck at the first hurdle :frown:
 
  • #5
150
0
wrt. both sides. Use subst.
 
  • #6
HallsofIvy
Science Advisor
Homework Helper
41,833
962
gammamcc, I think you are missing the point. The problem is not to integrate [itex]d^2r/dt^2= 1/r^2[/itex]. He's alread done that. That results in
[tex]v= \frac{dr}{dt}= \sqrt{c_1- \frac{2}{r}}[/tex]
The question is how to integrate that.
 
  • #7
150
0
ditto

Use separation of variables. (You wind up solving t explicitly, instead of r)
:bugeye:
 
  • #8
HallsofIvy
Science Advisor
Homework Helper
41,833
962
Yes, of course. The whole question is about integrating
[tex]\int \frac{dr}{\sqrt{c_1- \frac{2}{r}}}dr[/tex]
 
  • #9
150
0
Your attachment was not immediately available.......

This is all one expects typically. If you want more fun, put in some angular momentum.
 
Last edited:
  • #10
Ok, so start with

[tex]r^{\prime \prime} = - \frac{1}{r^2} [/tex]

times by r' and integrate to get

[tex]\frac{1}{2}r^{\prime 2} + C = \frac{1}{r} [/tex]

square it to get

[tex]\frac{1}{r^2} = (\frac{r^{\prime 2}}{2} + C)^2 [/tex]

but we can sub r'' for the [tex]1/r^2[/tex], so

[tex]-r^{\prime \prime} = (\frac{r^{\prime 2}}{2} + C)^2 [/tex]

Now make the sub

[tex]r^{\prime} = \sqrt{2C} \tan{v} [/tex]

to give us

[tex]\sqrt{2C} v^{\prime} \sec^2{v} + C^2 \sec^4{v} = 0[/tex]

or in other words

[tex]v^{\prime} + \sqrt{\frac{C^3}{2}} \sec^2{v} = 0[/tex]

now multiply by cos^2 and integrate, and you get

[tex]v + \frac{1}{2} \sin{2v} = \kappa - \sqrt{2 C^3} t [/tex]

(kappa is a constant). So, what you can do is create a function, [tex]F[/tex], such that it's inverse
is defined by

[tex]F^{-1}(p) = p + \frac{1}{2} \sin{2p} [/tex]

and you can express v as

[tex]v(t) = F(\kappa - \sqrt{2C^3}t) [/tex]

and then

[tex]r(t) = \sqrt{2C} \int{ \tan{F(\kappa - \sqrt{2C^3}t)} dt} [/tex]

or more elegently written as

[tex]r(t) = (2 \varphi)^{\frac{1}{3}} \int{ \tan{F(\kappa + \varphi t)} dt} [/tex]

where [tex]\varphi = -\sqrt{2 C^3} [/tex].


...as far as I can tell... although it's late and I'm rather tired right now...
 
Last edited:
  • #11
150
0
not trivial

This does not seem to make an explicit (or even convenient) solution for r any easier. Know this is not a trivial problem cf. Goldstein's "Classical Mechanics" pp. 98-102 (2nd ed.).
 
  • #12
I've got a feeling it won't be possible to get a completely explicit solution for r(t). It's more a question of getting the most convenient implicit solution. One thing: the function

[tex]F^{-1}(p) = p + \sin{p} [/tex]

is, as far as I can tell, one-to-one: i.e., the gradient is never negative, and there shouldn't be any ambiguity when inverting it.

I think it makes a good candidate for a named function. :wink:
 
Last edited:
  • #13
HallsofIvy
Science Advisor
Homework Helper
41,833
962
What it reduces to is an "elliptic" integral that cannt be done in terms of elementary functions. I remember seeing ("B.C."- "before computers) a 6-foot shelf of "Values of the elliptic integrals of the first and second kinds" in a university library.
 
  • #14
arildno
Science Advisor
Homework Helper
Gold Member
Dearly Missed
9,970
133
I will do tihs for you:
[tex]\ddot{r}=-\frac{K}{r^{2}}\to\frac{1}{2}\dot{r}^{2}=\frac{K}{r}+C[/tex]
Or, not bothering about the sign right now, or the name on the constants, we have:
[tex]\dot{r}=\sqrt{\frac{2K+Cr}{r}}\to\int\sqrt{\frac{r}{2K+Cr}}dr=t+A[/tex]
Now, make the following substitution:
[tex]u=\sqrt{\frac{r}{2K+Cr}}\to{r}=\frac{2Ku^{2}}{1-Cu^{2}}[/tex]
Therefore, we have:
[tex]\frac{dr}{du}=\frac{4Ku(1-Cu^{2})+4KCu^{3}}{(1-Cu^{2})^{2}}=\frac{4Ku}{(1-CKu^{2})^{2}}[/tex]
Therefore, you get the integral in the equation transformed to:
[tex]\int\frac{4Ku^{2}du}{(1-KCu^{2})^{2}}=A+t[/tex]
In order to perform the integration, we note that:
[tex]\int\frac{4Ku^{2}du}{(1-KCu^{2})^{2}}=\frac{2u}{C(1-KCu^{2})}-\frac{2}{C}\int\frac{du}{1-KCu^{2}}[/tex]

You will therefore get an implicit equation for u as a function of t, and hence, r will also be only implicitly determined as a function of time.
The indefinite integral remaining will either be an arctan- or artanh-expression, depending on the sign of the product KC.
 
Last edited:
  • #15
441
0
hmm maple can produce an explicit formula for r, however the formula is quite long and involves a number of "Rootof" operations. So you most likely will onlly have an eact solution for a couple of select points.
 
  • #16
Hi, I'm trying to solve Newton Universal law of Gravitation as a differential equation in one dimension with two objects. One mass is fixed at distance r = 0, whilst the other is some initial distance away. I'm using distance rather than displacement and also ignored the proportionality factor to make the formula simpler. So in essence, this is the differential equation I want to solve:

[tex] a = \frac{d^2r}{dt^2} = -\frac{1}{r^2} [/tex]

However, I've hit a snag when trying to integrate for the 2nd time (see the image file). I have used two substitutions (a and b), and the last line in the file suggests that I will obtain the inverse function that will contain hyperbolic and polynomial terms, i.e. I won't be able to determine the original function of r with respect to t.

Can anyone help?
solutidn a=1/r^2 is r=b t^2/3
 
Top