Répondre à la discussion
Affichage des résultats 1 à 2 sur 2

Runge Kutta (et Fortran)



  1. #1
    maki00

    Runge Kutta (et Fortran)

    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

    -----


  2. #2
    maki00

    Talking Re : Runge Kutta (et Fortran)

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

Sur le même thème :

Discussions similaires

  1. Runge Kutta
    Par ydethe dans le forum Mathématiques du supérieur
    Réponses: 2
    Dernier message: 11/10/2011, 15h14
  2. runge kutta et equa diff du second ordre
    Par Nertius dans le forum Mathématiques du supérieur
    Réponses: 4
    Dernier message: 05/01/2008, 16h44
  3. Fortran .....
    Par epsilox dans le forum Logiciel - Software - Open Source
    Réponses: 0
    Dernier message: 07/12/2007, 21h04
  4. Runge kutta pour la RFD
    Par machinefr dans le forum Mathématiques du supérieur
    Réponses: 4
    Dernier message: 28/10/2007, 21h58
  5. Matlab/fortran
    Par briscouille dans le forum Logiciel - Software - Open Source
    Réponses: 0
    Dernier message: 26/09/2007, 09h19