Solutions of the geodesic equations for bound test particle motion in a Schwarzschild field are expressed using Jacobian elliptic functions. Keplerian orbital elements are identified and related to a set of canonical constants. Relativistic Lagrange planetary perturbation equations are derived.