# PlanetPhysics/Generalized Coordinates for Free Motion

The differential equations for the motion of a particle under any forces when we use rectangular coordinates are known from Newston's laws of motion
$m{\ddot {x}}=F_{x}$ $m{\ddot {y}}=F_{y}$ $m{\ddot {z}}=F_{z}$
where $F_{x},F_{y},F_{z}$ are the components of the actual forces on the particle resolved parallel to each of the fixed rectangular axes, or rather their equivalents $m{\ddot {x}},m{\ddot {y}},m{\ddot {z}}$, are called the effective forces on the particle. They are of course a set of forces mechanically equivalent to the actual forces acting on the particle.
The equations of motion of the particle in terms of any other system of coordinates are easily obtained.
Let $q_{1},q_{2},q_{3}$, be the coordinates in question. The appropriate formulas for transformation of coordinates express $x,y,z$ in terms of $q_{1},q_{2},q_{3}$.
$x=f_{1}(q_{1},q_{2},q_{3}),\,\,\,\,y=f_{2}(q_{1},q_{2},q_{3}),\,\,\,\,z=f_{3}(q_{1},q_{2},q_{3})$
For the component velocity ${\dot {x}}$ we have
${\dot {x}}={\frac {\partial x}{\partial q_{1}}}{\dot {q_{1}}}+{\frac {\partial x}{\partial q_{2}}}{\dot {q_{2}}}+{\frac {\partial x}{\partial q_{3}}}{\dot {q_{3}}}$
and ${\dot {x}},{\dot {y}},{\dot {z}}$ are explicit functions of $q_{1},q_{2},q_{3},{\dot {q_{1}}},{\dot {q_{2}}},{\dot {q_{3}}}$ linear and homogeneous in terms of ${\dot {q_{1}}},{\dot {q_{2}}},{\dot {q_{3}}}$.
We may note in passing that it follows from this fact that ${\dot {x}}^{2},{\dot {y}}^{2},{\dot {z}}^{2}$ are homogeneous quadratic functions of ${\dot {q_{1}}},{\dot {q_{2}}},{\dot {q_{3}}}$.
Obviously
${\frac {\partial {\dot {x}}}{\partial {\dot {q_{1}}}}}={\frac {\partial x}{\partial q_{1}}}$
and since
${\frac {d}{dt}}{\frac {\partial x}{\partial q_{1}}}={\frac {\partial ^{2}x}{\partial q_{1}^{2}}}{\dot {q_{1}}}+{\frac {\partial ^{2}x}{\partial q_{2}\partial q_{1}}}{\dot {q_{2}}}+{\frac {\partial ^{2}x}{\partial q_{3}\partial q_{1}}}{\dot {q_{3}}}$
and
${\frac {\partial {\dot {x}}}{\partial q_{1}}}={\frac {\partial ^{2}x}{\partial q_{1}^{2}}}{\dot {q_{1}}}+{\frac {\partial ^{2}x}{\partial q_{1}\partial q_{2}}}{\dot {q_{2}}}+{\frac {\partial ^{2}x}{\partial q_{1}\partial q_{3}}}{\dot {q_{3}}}$
${\frac {d}{dt}}{\frac {\partial x}{\partial q_{1}}}={\frac {\partial {\dot {x}}}{\partial q_{1}}}$
Let us now find an expression for the work $\delta _{q_{1}}W$ done by the effective forces when the coordinate $q_{1}$ is changed by an infinitesimal amount $\delta q_{1}$ without changing $q_{2}$ or $q_{3}$. If $\delta x,\delta y,\delta z$ are changes thus produced in $x,y,z$, obviously from the definition of work
$\delta _{q_{1}}W=m\left[{\ddot {x}}\delta x+{\ddot {y}}\delta y+{\ddot {z}}\delta z\right]$
if expressed in rectangular coordinates. We need, however, to express $\delta _{q_{1}}W$ in terms of our coordinates $q_{1},q_{2},q_{3}$.
$\delta _{q_{1}}W=m\left[{\ddot {x}}{\frac {\partial x}{\partial q_{1}}}+{\ddot {y}}{\frac {\partial y}{\partial q_{1}}}+{\ddot {z}}{\frac {\partial z}{\partial q_{1}}}\right]\delta q_{1}$
Now
${\ddot {x}}{\frac {\partial x}{\partial q_{1}}}={\frac {d}{dt}}\left({\dot {x}}{\frac {\partial x}{\partial q_{1}}}\right)-{\dot {x}}{\frac {d}{dt}}{\frac {\partial x}{\partial q_{1}}}$
but from earlier definitions
${\frac {\partial x}{\partial q_{1}}}={\frac {\partial {\dot {x}}}{\partial {\dot {q_{1}}}}}\,\,\,\,\,and\,\,\,\,\,{\frac {d}{dt}}{\frac {\partial x}{\partial q_{1}}}={\frac {\partial {\dot {x}}}{\partial q_{1}}}$
Hence
${\ddot {x}}{\frac {\partial x}{\partial q_{1}}}={\frac {d}{dt}}\left({\dot {x}}{\frac {\partial {\dot {x}}}{\partial {\dot {q_{1}}}}}\right)-{\dot {x}}{\frac {\partial {\dot {x}}}{\partial q_{1}}}={\frac {d}{dt}}{\frac {\partial }{\partial {\dot {q_{1}}}}}\left({\frac {{\dot {x}}^{2}}{2}}\right)-{\frac {\partial }{\partial q_{1}}}\left({\frac {{\dot {x}}^{2}}{2}}\right)$
and therefore
   $\delta _{q_{1}}W=\left[{\frac {d}{dt}}{\frac {\partial T}{\partial {\dot {q_{1}}}}}-{\frac {\partial T}{\partial q_{1}}}\right]\delta q_{1}$
where
$T={\frac {m}{2}}\left[{\dot {x}}^{2}+{\dot {y}}^{2}+{\dot {z}}^{2}\right]$
and is the kinetic energy of the particle.
To get our differential equation we have only to write the second member of (1) equal to the work done by the actual forces when $q_{1}$ is changed by $\delta q_{1}$.
If we represent the work in question by $Q_{1}\delta q_{1}$, our equation is
   ${\frac {d}{dt}}{\frac {\partial T}{\partial {\dot {q_{1}}}}}-{\frac {\partial T}{\partial q_{1}}}=Q_{1}$
and of course we get such an equation for every coordinate. Even though we derived this differential equation for a single particle in free motion, it is the same for a systems of particles, except the kinetic energy is for all the particles in the system, which brings us to Lagrange's equations
   $Q_{i}={\frac {d}{dt}}\left({\frac {\partial T}{\partial {\dot {q_{i}}}}}\right)-{\frac {\partial T}{\partial q_{i}}}$
In any concrete problem, $T$ must be expressed in terms of $q_{1},q_{2},q_{3}$, and their time derivatives before we can form the expression for the work done by the effective forces. $Q_{1}\delta q_{1},Q_{2}\delta q_{2},Q_{3}\delta q_{3}$, the work done by the actual forces, must be obtained from direct examination of the problem.