Version imprimable
Voici deux posts avec deux implémentations (C++ et Java) de la méthode de Runge-Kutta :
C++ :
http://www.developpez.net/forums/sho...77&postcount=2
Java :
http://www.developpez.net/forums/sho....php?p=2528269
Ces posts datent, je viens seulement de penser à les poster ici.
A noter que j'ai proposé une solution pour résoudre une équation de la forme :
dy/dx = f(x,y) (je vous passe les détails sur les hypothèses à faire sur les fonctions)
cool... :king:
On va pouvoir poster des filtres a base de PDE...
Pour tous ceux qui cherchent des algorithmes performants pour le calcul numérique, je ne peux que recommander "Numerical Recipes" déjà souvent cité sur ce forum, et plus encore le site www.netlib.org .
En ce qui concerne l'intégration des équations différentielles ordinaires, il n'y a pas que Runge-Kutta, et je vous recommande de faire un tour sur www.unige.ch/~hairer/software.html
Bonne chance à tous.
Jean-Marc Blanc
Bonjour,
je viens d'essayer le lien java qui ne marche pas :(
Grrrrrrrrr :(Citation:
Envoyé par ToTo13
Je sais pourquoi... J'ai une sauvegarde quelque part, je remettrais quand je pourrais.
Si quelqu'un est intéressé à l'utiliser en Fortran, j'ai retrouvé ça dans mon répertoire à vieilleries:
Subroutine D006(SM,N,T,DT,Y)
C
C***********************************************************************
C
C Bibliothèque JMBPACK
C
C Série D: Equations différentielles
C
C***********************************************************************
C
C Méthode de Runge-Kutta d'ordre 4 (1/6 + 2/6 + 2/6 + 1/6)
C
C***********************************************************************
C
C Version 3.0: Jean-Marc Blanc, juillet 2004
C
C **********************************************************************
C
C SM Nom d'un sous-programme fourni par l'utilisateur pour le
C calcul des dérivées à partir des valeurs de la variable
C indépendante et des variables d'état. Ce sous-programme
C doit être de la forme
C Subroutine SM(N,T,Y,DY)
C DY(1)= ...
C DY(2)= ...
C .
C .
C DY(N)= ...
C Return
C End
C et avoir été déclaré dans une instruction External.
C
C N Ordre du système différentiel.
C
C T Variable indépendante.
C
C DT Pas d'intégration.
C
C Y Vecteur des variables d'état.
C
C **********************************************************************
C
Implicit None
C
Integer N
Real*8 DT,T,Y(N)
C
Integer I
Real*8 DY(N),K1(N),K2(N),K3(N),K4(N),Z(N)
C
Call SM(N,T,Y,DY)
Do I=1,N
K1(I)=DT*DY(I)
Z(I)=Y(I)+K1(I)/2.d0
End Do
C
Call SM(N,T+DT/2.d0,Z,DY)
Do I=1,N
K2(I)=DT*DY(I)
Z(I)=Y(I)+K2(I)/2.d0
End Do
C
Call SM(N,T+DT/2.d0,Z,DY)
Do I=1,N
K3(I)=DT*DY(I)
Z(I)=Y(I)+K3(I)
End Do
C
Call SM(N,T+DT,Z,DY)
Do I=1,N
K4(I)=DT*DY(I)
Y(I)=Y(I)+(K1(I)+2.d0*K2(I)+2.d0*K3(I)+K4(I))/6.d0
End Do
C
T=T+DT
Return
C
End
Salut à tous
Jean-Marc Blanc
On a réussi à recupérer les données :Citation:
Envoyé par ToTo13
http://www.developpez.net/forums/sho....php?p=2528269