Centre Universitaire de
Tipaza Institut des Sciences
Licence LMD L2 (Semestre 4)
Filières : Génie des Procédés
Module : Méthodes Numériques
TP 4 : Résolution des équations différentielles ordinaires
Durée du TP : 2 séances de 1h30
But du TP :
Durant ce TP, nous allons mettre en œuvre les algorithmes des méthodes de résolution
des équations différentielles ordinaires étudiées pendant le cours : la méthode d’Euler
et la méthode de Range Kutta. Il sera ensuite question de comparer les deux méthodes
avec la solution exacte.
Rappel sur les différentes méthodes :
Méthode d’Euler :
$ % &' ( &')*
+1 % +0 , $- .&0 , +0 /
+2 % +1 , $- .&1 , +1 /
+'01 % +' , $-.&' , +' /
Méthode de Range Kutta d’ordre 4 :
11 % - .&0 , +0 /
$ $
12 % - 2&0 , , +0 , 3 11 4 %
2 2
$ $
13 % - 2&0 , , +0 , 3 12 4 %
2 2
14 % - .&0 , $, +0 , $13 / %
$
+1 % +0 , 3 .11 , 212 , 213 , 15 / %
6
&1 % &1 , $ %
11 % -.&1 , +1 / %
$ $
12 % - 2&1 , , +1 , 3 11 4 %
2 2
$ $
13 % - 2&1 , , +1 , 3 12 4 %
2 2
14 % -.&1 , $, +1 , $13 / %
$
+2 % +1 , 3 .11 , 212 , 213 , 15 / %
6
Travail demandé :
Soit l’équation différentielle
dy
= yt 3 − 1.5y
dt
avec +.0/ % 1 78 t 9 :0 2;
Résoudre cette équation avec un pas d’intégration $ % 0,2 en utilisant :
‐ La méthode d’Euler
‐ La méthode de Range Kutta d’ordre 4 (RK4)
Calculer la solution exacte de l’équation et comparer avec les approximations
précédentes, quel est votre commentaire.
TP Méthodes Numériques