B(t) = \left( \begin{array}{cc}0 &1\\ - w^2(t) & 0\end{array} \right) où w^2(t) = \omega_0^2(1+\varepsilon \cos \omega t ) (qui est positive pour \varepsilon assez petit)

Z = \left( \begin{array}{c} \theta \\ \dot{\theta} \end{array} \right)