Let M be a manifold and let nabla be a symmetric linear connection on the tangent bundle of M.  If C is a curve in M with tangent vector T, then the geodesic equations for C with respect to the connection  nabla is the system of second order ODEs defined by nabla_T(T) = 0.

The procedure GeodesicEquations(C, Gamma, t) returns the vector nabla_T(T).

This command is part of the DifferentialGeometry:-Tensor package, and so can be used in the form GeodesicEquations(...) only after executing the command with(DifferentialGeometry) and with(Tensor) in that order.  It can always be used in the long form DifferentialGeometry:-Tensor:-GeodesicEquations.