On Szebehely's problem for holonomic systems involving generalized potential functions

We consider the problem of finding the generalized potential function V = U _i(q ¹, q ²,..., q ⁿ)q ⁱ + U(q ¹, q ²,...;q ⁿ) compatible with prescribed dynamical trajectories of a holonomic system. We obtain conditions necessary for the existence of solutions to the problem: these can be cast into a system of n – 1 first order nonlinear partial differential equations in the unknown functions U ₁, U ₂,...;, U _n, U. In particular we study dynamical systems with two degrees of freedom. Using adapted coordinates on the configuration manifold M ₂ we obtain, for potential function U(q ¹, q ²), a classic first kind of Abel ordinary differential equation. Moreover, we show that, in special cases of dynamical interest, such an equation can be solved by quadrature. In particular we establish, for ordinary potential functions, a classical formula obtained in different way by Joukowsky for a particle moving on a surface.Work performed with the support of the Gruppo Nazionale di Fisica Matematica (G.N.F.M.) of the Italian National Research Council.