On an n-dimensional Riemannian manifold ${\displaystyle M}$ , the geodesic equation written in a coordinate chart with coordinates ${\displaystyle x^{a}}$  is:

${\displaystyle {\frac {d^{2}x^{a}}{ds^{2}}}+\Gamma _{bc}^{a}{\frac {dx^{b}}{ds}}{\frac {dx^{c}}{ds}}=0}$ 

where the coordinates x^a(s) are regarded as the coordinates of a curve γ(s) in ${\displaystyle M}$  and ${\displaystyle \Gamma _{bc}^{a}}$  are the Christoffel symbols. The Christoffel symbols are functions of the metric and are given by:

${\displaystyle \Gamma _{bc}^{a}={\frac {1}{2}}g^{ad}\left(g_{cd,b}+g_{bd,c}-g_{bc,d}\right)}$ 

where the comma indicates a partial derivative with respect to the coordinates:

${\displaystyle g_{ab,c}={\frac {\partial {g_{ab}}}{\partial {x^{c}}}}}$ 

As the manifold has dimension ${\displaystyle n}$ , the geodesic equations are a system of ${\displaystyle n}$  ordinary differential equations for the ${\displaystyle n}$  coordinate variables. Thus, allied with initial conditions, the system can, according to the Picard–Lindelöf theorem, be solved. One can also use a Lagrangian approach to the problem: defining

${\displaystyle L={\sqrt {g_{\mu \nu }{\frac {dx^{\mu }}{ds}}{\frac {dx^{\nu }}{ds}}}}}$ 

and applying the Euler–Lagrange equation.

As the laws of physics can be written in any coordinate system, it is convenient to choose one that simplifies the geodesic equations. Mathematically, this means a coordinate chart is chosen in which the geodesic equations have a particularly tractable form.

When the geodesic equations can be separated into terms containing only an undifferentiated variable and terms containing only its derivative, the former may be consolidated into an effective potential dependent only on position. In this case, many of the heuristic methods of analysing energy diagrams apply, in particular the location of turning points.

Solving the geodesic equations means obtaining an exact solution, possibly even the general solution, of the geodesic equations. Most attacks secretly employ the point symmetry group of the system of geodesic equations. This often yields a result giving a family of solutions implicitly, but in many examples does yield the general solution in explicit form.

In general relativity, to obtain timelike geodesics it is often simplest to start from the spacetime metric, after dividing by ${\displaystyle ds^{2}}$  to obtain the form

${\displaystyle -1=g_{\mu \nu }{\dot {x}}^{\mu }{\dot {x}}^{\nu }}$ 

where the dot represents differentiation with respect to ${\displaystyle s}$ . Because timelike geodesics are maximal, one may apply the Euler–Lagrange equation directly, and thus obtain a set of equations equivalent to the geodesic equations. This method has the advantage of bypassing a tedious calculation of Christoffel symbols.

• Geodesics of the Schwarzschild vacuum
• Mathematics of general relativity
• Transition from special relativity to general relativity