The function geodesic_eqns(coord, Tau, Cf2) generates (but does not solve) the Euler-Lagrange equations of the geodesics for a metric with Christoffel symbols of the second kind Cf2 and coordinate variables coord.  The equations are written in terms of the coordinate variable names as functions of the given parameter Tau. They are returned in the format of a list of equations.

Cf2 should be indexed using the cf2 indexing function provided by the tensor package.  It can be computed using the Christoffel2 routine.