Routhian reduction

Mixed Lagrangian/Hamiltonian for cyclic coordinates: R = Σ p_c q̇_c - L. Eliminates ignorable coords while keeping others in Lagrangian form.