- #1

My problem involves the precession of the perihelion of Mecury

F sub g = - (GMm)/r^2 * (1 + a/r) where a << r

I proved in previous parts d^2r/dt^2 - r*(dtheta/dt)^2 =

-GM/r^2 * (1 + a/r) [eqn 1] and r* d^2theta/dt^2 + 2*dr/dt*dtheta/dt = 0

I also used u(theta) = 1/r(t) to turn eqn 1 to d^2u/dtheta^2 +

u(1-GMa/l^2) = GM/l^2 where l = L/m.

Where I'm stuck is showing the solution is u(theta) = u(sub 0) *

(1+ ecos(n(theta - theta (sub 0))) where e and theta (sub 0) are constants of integration and u(sub 0) and n are in terms of a, G, M and l. I've tried many times but can't get it to work out.

F sub g = - (GMm)/r^2 * (1 + a/r) where a << r

I proved in previous parts d^2r/dt^2 - r*(dtheta/dt)^2 =

-GM/r^2 * (1 + a/r) [eqn 1] and r* d^2theta/dt^2 + 2*dr/dt*dtheta/dt = 0

I also used u(theta) = 1/r(t) to turn eqn 1 to d^2u/dtheta^2 +

u(1-GMa/l^2) = GM/l^2 where l = L/m.

Where I'm stuck is showing the solution is u(theta) = u(sub 0) *

(1+ ecos(n(theta - theta (sub 0))) where e and theta (sub 0) are constants of integration and u(sub 0) and n are in terms of a, G, M and l. I've tried many times but can't get it to work out.

Last edited by a moderator: