Bonjour,
J'ai un système d'équations différentielle non linéaire que je souhaite résoudre :
x''(t) = b.z - a.y (1)
x.y'(t) - c.y.x'(t) = - d.y (2)
x.z'(t) - c.z.x'(t) = e (3)
Avec x, y et z sont les fonctions et a, b, c, d et e des constantes.

Pour la résoudre j'utilise une méthode approchée, j'utilise l'approximation de la dérivée d'une fonction :
g'(t)=[g(t+h)-g(t)]/h
g''(t)=[g(t+2.h)-2.g(t+h)+g(t)]/h² avec h=0.001

J'injecte donc ces expressions dans mes équations, je réalise une résolution numérique par itération.
Mais je voudrais savoir à quel point je peux faire confiance à mes résultats, est-ce que l'erreur ne seras pas trop grande par hasard ??

Merci d'avance.