Deux oscillateurs parallèles couplés par un ressort

Les deux solides sont identiques et de masse m, la raideur du ressort est k, l₀ est la longueur du ressort à vide et l à un instant quelconque. La position des mobiles sont repérées respectivement, par (x₁ , d) et (x₂ ,−d). L’application de la deuxième de Newton conduit à un système d’équations différentielles couplées :

La résolution numérique de ce système est faite avec le package pst-ode. En fonction des conditions initiales, on obtient différents types de comportements pour le système.
Les 2 mobiles oscillent en opposition :

La danse des mobiles : si on modifie les conditions initiales, en particulier les vitesses on a un comportement un peu moins monotone des mobiles : globalement les 2 mobiles progressent vers l’avant mais avec des petits retours en arrière, comme dans une espèce de danse(il manque la musique).

Tous les fichiers sont dans l'archive :
Les-ressorts

Ce document est la suite de :
http://pstricks.blogspot.fr/2016/06/un-oscillateur-ressort-non-sinusoidal.html