Systeme d'equadiffs

[lajoie]

salut

svp, je fait un programme sur la résolution d'un système d'équations différentielles non linéaire à deux entrées, où j'utilise la méthode de Runge-Kutta d'ordre 4, mon système est de la forme
du/dt= f(u,v)
dv/dt= g(u,v)

dans l'algorithme, je calcul, les kf1, kf2, kf3, kf4 mais indépendamment du temps càd,

kf1= f( u(i), v(i) ) puis,
kf2=f( u(i)+h*kf1/2 ,v(i)+h*kf1/2 ), ...etc et je fais de même pour la fonction g ,

après je calcule : les solutions :u(i+1), v(i+1) et le temps:t(i)=i*h, ensuite je passe a l'itération suivante .

est ce que c'est juste cet algorithme. aidez moi svp, car j'en ai vraiment besoin

j'attends vos réponses, et je vous remercie infiniment .