On utilise un code de calcul avec intégration en temps. Le schéma retenu est: Simul. (les données ne sont pas exactement celles du Rafale)
\[ \boxed {\frac{\alpha ^{n+1}-2\alpha ^ n+\alpha ^{n-1}}{\Delta t^2}+c(\alpha ^ n-\alpha _0)=\frac{\varrho SL \vert V_ a\vert ^2}{2}c_{m0}(\alpha _ a).} \]
Dans un second temps nous avons utilisé l’algorithme des formes normales ce qui a permis de trouver un terme résonnant à l’ordre 3 qui donne naissance à un cycle limite (TSVP).