Quote by Anthony
I think the OP might be making life difficult for himself. With any form of variational problem, you should try to use Noether's theorem to integrate your equations up, rather than try to go from the EulerLagrange equations themselves. This will invariably make your life easier. Your Lagrangian density reads:
[tex] \mathcal{L} = \left( 1 \frac{1}{r}\right) \left( \frac{\mathrm{d} t}{\mathrm{d} \lambda}\right)^2  \left( 1 \frac{1}{r}\right)^{1} \left( \frac{\mathrm{d} r}{\mathrm{d} \lambda}\right)^2  r^2 \left( \frac{\mathrm{d} \theta}{\mathrm{d} \lambda}\right)^2  r^2 \sin^2\theta \left( \frac{\mathrm{d} \phi}{\mathrm{d} \lambda}\right)^2 [/tex]
using natural units. It is clear that the vector fields:
[tex] \frac{\partial}{\partial t} \quad \textrm{and}\quad \frac{\partial}{\partial \phi}[/tex]
are Killing, so Noether's theorem integrates up two of the EulerLagrange equations for you and gives you two constants of motion:
[tex] \left(1\frac{1}{r}\right) \frac{\mathrm{d} t}{\mathrm{d}\lambda} = \mathrm{const} \,(=E) \quad \textrm{and} \quad r^2 \sin^2\theta \frac{\mathrm{d} \phi}{\mathrm{d} \lambda} = \mathrm{const} \,(=h) \qquad (*) [/tex]
i.e. on a given geodesic, these quantities remain unchanged. Similarly, since [tex]\partial_\lambda \mathcal{L}=0[/tex], we know [tex] \mathcal{L}[/tex] remains constant, and we set it to {+1,1,0} depending on whether you're interested in timelike, spacelike or nulll geodesics. Call this constant k. Note that all our ODEs are now 1st order. Setting [tex]\theta = \pi/2[/tex] (validity can be deduced from the [tex]\theta[/tex] EL equation) and using (*) in [tex]\mathcal{L}=k[/tex] gives the ODE:
[tex] \left( \frac{\mathrm{d} r}{\mathrm{d}\lambda}\right)^2 = E^2  \left( 1\frac{1}{r}\right) \left(k + \frac{h^2}{r^2}\right)[/tex]
If you'd prefer to parameterise your geodesics using [tex]\phi[/tex], use the second of the constraints in (*) again and you get:
[tex] \frac{h^2}{r^4} \left( \frac{\mathrm{d} r}{\mathrm{d} \phi}\right)^2 = E^2  \left( 1\frac{1}{r}\right) \left( k + \frac{h^2}{r^2}\right) [/tex]
If you'd prefer to do all this using tensors, just apply Noether's theorem in the form: if [tex]L_V g=0[/tex] (i.e. [tex]V[/tex] is a Killing vector) then [tex]V^\mu \dot{x}^\nu g_{\mu\nu} = \mathrm{const}[/tex].

Exactly! This simply is based on the proposition given in post #389 and by this perfect explanation we are done here.
AB