Les équations de la dynamique de l'électron et du noyau dans le référentiel barycentrique s'écrivent:

((m_n*m_e)/(m_e+m_n))*(dtemps(vector(r);2)=-e^2/4*pi*epsilon_0*r^3)*vector(r)

-((m_n*m_e)/(m_e+m_n))*(dtemps(vector(r);2)=e^2/4*pi*epsilon_0*r^3)*vector(r)

Il s'agit de deux fois la même équation de la dynamique d'une particule fictive de masse

mu=(m_e*m_n)/(m_e+m_n)et de position vector(r)