Runge Kutta (et Fortran)

[invite642d375a]

Ma question concerne la méthode de runge kutta ordre 4 pour résoudre un système d'équations différentielles ordinaires (8 au total)
Lorsque j'appelle la sous-routine (en Fortran) contenant mon système, faut-il le modifier ou bien puis-je le rentrer tel quel

Pour donner un ordre d'idée pour l'instant j'ai ça (et pour info le programme renvoie des aberrations )

SUBROUTINE deriv(k,tau,t,y,dy)

USE vars

IMPLICIT NONE
REAL::t, tau, y(8), dy(8), k(2)

dy(1)=-k(1)*y(1)*(y(5)*Ctot)**(1./2.)-k3*y(5)+y0_DCE/tau-Qs*y(1)/V
dy(2)=-k(2)*(y(5)*Ctot)**(1./2.)*y(2)+k(1)*(y(5)*Ctot)**(1./2.)*y(1)+k3*y(5)-Qs*y(2)/V
dy(3)=k(1)*(y(5)*Ctot)**(1./2.)*y(1)+k(2)*(y(5)*Ctot)**(1./2.)*y(2)-Qs*y(3)/V
dy(4)=k(2)*(y(5)*Ctot)**(1./2.)*y(2)-Qs*y(4)/V
dy(5)=-k3*y(5)-Qs*y(5)/V+y0_CCl4/tau
dy(6)=k3*y(5)-Qs*y(6)/V
dy(7)=k3*y(5)-Qs*y(7)/V
dy(8)=y0_N2/tau-Qs*y(8)/V

END SUBROUTINE deriv
-----

[invite642d375a]

oubliez la question , j'ai trouvé ma réponse finalement!