Méthode d'Euler, équation différentielle

[totoflingueur]

Bonjour, je suis nouveau sur ce forum et je tiens à dire qu'il est très intéressant.

Alors voici mon problème : je suis en train de créer un logiciel de simulation du botafumeiro (pendule dont la longueur varie). Il m'a donc fallu résoudre une équation différentielle du 2° ordre non linéaire à coefficients variables (ça fait mal :marteau: )
Jai donc utilisé la methode d' Euler :
Pour simplifier, voici l'équation : f ''(t) + a(t)*f '(t) +b(t)*sin(f(t)) = 0
j'ai des conditions initiales, puis je fais la récurrence :
f(t+h) = f(t)+h*f '(t)
f '(t+h) = f '(t)+h*f ''(t)
f ''(t+h) = -a(t+h)*f '(t) - b(t+h)*sin(f(t+h))

Mon programme me donne des résultats totalement délirant, j'ai donc tester simplement avec un tableur. meme resultat, par contre si je fixe les coefficients en constante, la méthode fonctionne, j'obtiens une courbe pseudo périodique(c'est le pendule simple)
Du coup j'aimerais bien savoir si la méthode d'Euler fonctionne aussi avec des coefficients variables, et si oui si vous ne voyez pas une erreur dans mon application de la méthode.

Si vous pouvez m'indiquer des axes de recherche ce serait sympa.

merci d'avance