The direct product of integrable Hamiltonian systems is integrable. Pull each family's commuting integrals of motion back to the product phase space. Integrals from different factors have zero Poisson bracket because they use disjoint canonical coordinates, and their differentials are independent wherever the differentials within each factor are independent.
In canonical coordinates, Hamilton's equations are
A function is an integral of motion if it is constant along every Hamiltonian flow, equivalently
for a time-independent . The system is an integrable Hamiltonian system when it has first integrals whose differentials are linearly independent almost everywhere and which are in mutual Poisson involution: for all .
Differentiating the spherical coordinate system
gives the orthogonal decomposition
where
Thus the scale factors of orthogonal coordinates are
Since the total differential satisfies
comparison of coefficients gives the gradient