# Integrating the inverse square law

1. Apr 11, 2009

### haaj86

Hi, I am trying to solve the following problem (but note you can skip the details and look at the last equation and help me with the integral): a small planet (mass m) is found at rest relative to the Sun (mass M) with a separation distance $$x_{0}$$. Consider now the small planet move in response to the gravitational force of the Sun (the sun is assumed so massive such that it doesnt move) what is the position function of the small planet as a function of time t? So,

$$m\frac{dv}{dt}=-\frac{GMm}{x^{2}}$$

$$\int^{v}_{0}{mvdv}=-\int^{x}_{x_{0}}{\frac{GMm}{x^{2}}}$$

$$\frac{1}{2}mv^{2}=GMm(\frac{1}{x}-\frac{1}{x_{0}})$$

Therefore,
$$\frac{dx}{dt}= \pm[2GM(\frac{1}{x}-\frac{1}{x_{0}})]^{\frac{1}{2}}$$

We take the minus sign because we know that the velocity is in the negative direction (the planet moves toward the sun), hence:
$$t=-\frac{1}{\sqrt{2GM}}\int^{x}_{x_{0}}{(\frac{1}{x'}-\frac{1}{x_{0}})^{-\frac{1}{2}}}dx'$$
So far I haven't found this integral solved in any text book, please if the calculation is long and you can't type it then tell me where to find it.

2. Apr 11, 2009

### Cantab Morgan

I'm not sure I followed your very first step. You may call me a stickler, but... what is the differential on the right hand side integral?

3. Apr 11, 2009

### Staff: Mentor

It's much easier to use conservation of energy. The potential energy in this case is

$$U(x) = - \frac {GMm}{x}$$

(where we define U = 0 at x = infinity). The total of potential plus kinetic energy is constant.

4. Apr 12, 2009

### haaj86

It is simply F=ma, where $$F=-\frac{GMm}{x^{2}}$$ and $$ma=m\frac{dv}{dt}$$

You are definitely right, I was going to get to equation three by using the fact that the total energy is constant, but I decided to start by integrating the e.o.m to show exactly what I mean by integrating the inverse square law. I am really stuck with the last integral and I am shocked to not being able to find it solved. The closest solution I found was to assume that $$x_{0}$$ is so large such that we can ignore the $$1/x_{0}$$ term and approximate the integral. But I really want to find the exact solution as x(t)= ....

5. Apr 12, 2009

### Hootenanny

Staff Emeritus
I believe he was referring to the right hand side of
which does not appear to have a variable of integration. However, I suspect that missing off the dx was just a typo.

6. Apr 12, 2009

### Cyosis

Defining a=1/x0 the integral can be written as:

$$\int \frac{1}{\sqrt{\frac{1}{x}-a}}\, dx =\int \frac{\sqrt{x}}{\sqrt{1-ax}}\, dx[/itex] Now use the substitution $u=\sqrt{1-ax}$. Ignoring constants you should arrive at an integral of the form [tex]\int \sqrt{1-u^2}\,du$$ which is easily solved by using another substitution, u=sin(z).

Edit: After some work I got to the following expression. Barring any mistakes it doesn't look like you can solve this for x.
$$t=-\frac{1}{\sqrt{2GM}}\int^{x}_{x_{0}}{(\frac{1}{x}-\frac{1}{x_{0}})^{-\frac{1}{2}}}dx= \frac{x_0 \sqrt{x_0}}{\sqrt{2GM}} \left(\sqrt{\frac{x}{x_0}(1-\frac{x}{x_0})}+\arcsin(\sqrt{1-\frac{x}{x_0}})\right)$$

Last edited: Apr 12, 2009
7. Apr 12, 2009

### haaj86

You are absolutely right, sorry about that.

Thank you so much, I worked through it following your steps and I got a very similar answer which is

$$t=-\frac{1}{\sqrt{2GM}}\int^{x}_{x_{0}}{(\frac{1}{x}-\frac{1}{x_{0}})^{-\frac{1}{2}}}dx= \frac{x_0 \sqrt{x_0}}{\sqrt{2GM}} \left(\frac{1}{2} sin(2arcsin(\sqrt{1-\frac{x}{x_{0}}})+\arcsin(\sqrt{1-\frac{x}{x_0}})\right)$$

I am not sure if this is what you basically got and you managed to reduce the first term to the form you have in your solution or that I have basically made a mistake somewhere. If you are willing to help me further to get it to the form of your solution then please let me know so that I can type my full solution and you can help me to spot the mistake I made. But for now, regardless whether I got the first term right or not the fact will remain that I will not be able solve it for x, so you have answered my main question, thank you for your help.

8. Apr 12, 2009

### Cyosis

Your answer is exactly the same. Putting $$u=\arcsin(\sqrt{1-\frac{x}{x_0}})$$ for convenience. $$\frac{1}{2} \sin(2u)=\cos(u) \sin(u)=\sqrt{1-\sin^2 u} \sin u$$. If you fill in u now you will see that the arcsin cancels and you will get the same answer as me.

9. Apr 12, 2009

### haaj86

Yes it works perfectly, thank you for all your help.

haaj86

10. Apr 12, 2009

### Cyosis

You're welcome.

11. Apr 12, 2009

### Cantab Morgan

Please look again, it's not a typo. Can you see what I was gently driving at? The differential is not dx. It is actually dt. This attempt to solve the differential equation fails, and produces such an intractable integral, because of this error in the very first step.

12. Apr 13, 2009

### Cyosis

You missed the chain rule $$\frac{dv}{dt}=\frac{dv}{dx} \frac{dx}{dt}=v \frac{dv}{dx}$$. So the differential is dx not dt.

13. Apr 13, 2009

### Hootenanny

Staff Emeritus
If it is not a typo, then it is just plain wrong for the reason that Cyosis points out:

14. Apr 14, 2009

### Cantab Morgan

Aha. I am wrong. I was going to suggest multiplying both sides by v, but I now see more carefully that you have already done that.