Gravitational force with 1/r^3

[nearc]

apologies if this has been covered elsewhere in the forums, but i was unable to find it. I've seen the usual gravitational force equation [i.e. F=gMm/r^2] with essentially r^3 in the bottom. an example of this can be seen in http://www.mathworks.com/moler/exm/chapters/orbits.pdf on pages 8 and 9.

i'm looking for any info [paper/webpages] that might explain how to use this approach, thanks

 

[rbj]

it's a vector equation and the author (who happens to be the guy who created the first MATLAB) is not being conventional with his math symbols.

Newtons gravitational law can first be stated:

$$F = G \frac{m_1 m_2}{r^2}$$

where $F$ is the magnitude of the force vector that runs along the line that connects $m_1$ and $m_2$. and $r$ is the Euclidian distance between them. this is a scaler equation and is the simple inverse-square gravitational equation you see most everywhere.

now suppose that the body with mass $m_1$ is located at point $\vec{p_1}$ and the body with mass $m_2$ is located at point $\vec{p_2}$. the coordinates of the points are

$$\vec{p_1} = ( \ x_1, \ y_1, \ z_1 \ )$$

$$\vec{p_2} = ( \ x_2, \ y_2, \ z_2 \ )$$

the distance between the two points $\vec{p_2}$ and $\vec{p_1}$ is

$$r = | \vec{p_1} - \vec{p_2} | = \sqrt{(x_1-x_2)^2 +(y_1-y_2)^2 +(z_1-z_2)^2 }$$

now, the difference vector that connects $\vec{p_2}$ to $\vec{p_1}$ is $\vec{p_1}-\vec{p_2}$, so $r$ is just the magnitude of that difference vector.

lastly to get a force vector $\vec{F}$ that has magnitude $F$, we need to create a unit vector that runs along the line that connects $\vec{p_2}$ to $\vec{p_1}$. that unit vector is the difference vector $\vec{p_1}-\vec{p_2}$ normalized by dividing by its length $r=|\vec{p_1}-\vec{p_2}|$. this is where the extra power in $r^3$ comes from.

so the force acting on the body with mass $m_1$ located at point $\vec{p_1}$ by the body with mass $m_2$ located at point $\vec{p_2}$ is

$$\vec{F} = -G \frac{m_1 m_2}{r^2} \frac{\vec{p_1}-\vec{p_2}}{|\vec{p_1}-\vec{p_2}|}$$

or

$$\vec{F} = -G \frac{m_1 m_2}{|\vec{p_1}-\vec{p_2}|^3} (\vec{p_1}-\vec{p_2})$$

that is a vector equation, but it is still the inverse-square law. the extra power of $r$ in the denominator is needed to normalize the vector connecting point $\vec{p_2}$ to point $\vec{p_1}$