L'équation de Duffing est décrite par cette équation différentielle non-linéaire :
Cet oscillateur forcé peut-être écrit comme un système d'équations d'état à trois dimensions :
Te=0.02 Tfin=200 g=0.3 a=0.15 w=1 To=2*%pi/w ci1=0.1 ci2=0.1 ci3=0