Systèmes différentiels : Méthode de Runge Kutta ordre 4 - équation de Lorenz / Navier-Stockes & attracteur de Lorenz

[Chaka]

Bonjour,
Dans le cadre d'un projet d'étude (3ème année de Licence), je voudrais mettre en place, sur Python, la méthode de RK4 pour résoudre les équations de Lorentz numériquement.
J'ai réussi à le faire sans trop de difficulté avec les méthodes de Newton, mais pour RK4, je bloque !! Le calcul des intermédiaires m'est compliqué.
Je vous mets ci-joint mes scripts Python.
Merci par avance pour vos retours.

Chaka 

Projet P.py

[Chaka]

Il m'est compliqué notamment car ma fonction lorenz est de R^3 dans R^3 mais chaque composante nécessite les autres pour son calcul et ma fonction RK4 prend des valeurs dans R si je ne dis pas d'idioties.

[Denis CAMUS]

Voir https://www.e-bahut.com/topic/57697-systèmes-différentiels-méthode-de-runge-kutta-ordre-4-équation-de-lorenz-navier-stockes-attracteur-de-lorenz/