calculate the geodesic equations for a symmetric linear connection on the tangent bundle - Maple Programming Help

Online Help

All Products    Maple    MapleSim

Home : Support : Online Help : Mathematics : DifferentialGeometry : Tensor : DifferentialGeometry/Tensor/GeodesicEquations

Tensor[GeodesicEquations] - calculate the geodesic equations for a symmetric linear connection on the tangent bundle

Calling Sequences

     GeodesicEquations (C, Gamma, t)


   C       - a list of functions of a single variable, defining the components of a curve on a manifold M with respect to a given system of coordinates

   Gamma   - a connection on the tangent bundle to a manifold M

   t       - the curve parameter




See Also



Let M be a manifold and let 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  is the system of second order ODEs defined by TT=0.


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


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.


with(DifferentialGeometry): with(Tensor):


Example 1.

First create a 2-dimensional manifold M and define a connection on the tangent space of M.

DGsetup([x, y], M);

frame name: M

M > 

Gamma := Connection(a*(D_x &t dx &t dy + D_x &t dy &t dx) + b*y*D_x &t dy &t dy);




To determine the geodesic equations for this connection we first define a curve on M by specifying a list of functions of a single variable t.

M > 

C := [x(t), y(t)];




The program GeodesicEquations returns a vector whose components are the components of the geodesic equations.

M > 

V := GeodesicEquations(C, Gamma, t);




To solve these geodesic equations use DGinfo to obtain the coefficients of V as a list. Pass the result to dsolve to solve this system of 2 second order ODEs. See also DGsolve.

M > 

DE := Tools:-DGinfo(V, "CoefficientSet");


M > 

simplify(dsolve(DE, explicit));



See Also

DifferentialGeometry, Tensor, Christoffel, Connection, CovariantDerivative, DGinfo, DirectionalCovariantDerivative, ParallelTransportEquations