IBRAHIM ABDILLAHI choukri

DEVOIR MAISON

1°) y''+ay=0 , y(0)=1 , y'(0)=0

On cherche à calculer numériquement l’équation différentiel au dessus.Pour cela, on applique des

quatre programmes par les différentes méthodes qui sont ( par la méthode d' Euler, par schéma de
krank-Nicholason, par la méthode d’Euler modifiée, et par la méthode de Runge et Kutta d'ordre).

On pose v=y' donc v'=y''=-ay et l’équation a résoudre s’écrit :

y'= v
avec y(0)=1 , v(0)=0
v'=-ay

Méthode krank-nicholason :
Méthode d'Euler modifie :

Méthode d’Euler
 méthode de Runge_Kutta d'ordre 4 :

n 100 1000 10000 100000

1.10-4
h 10-1 10-2 10-3
Méthode d’Euler 81,83 123,36 142,53 146 ,68
Méthode d’Euler 4,2.105 12,1 5 3,6
modifie
Méthode 6,285.10-6 3,5812.10-6 1,871.10-7 1,3347.10-7
krank_nicholson
Méthode de 6,1082.10-1 2,6092.10-2 8,3049.10-2 2,6347.10-5
Runge_Kutta
d'ordre 4
a 100 1000 10000 100000
La courbe de l’erreur en fonction de h en échelle logarithmique :
→ Pour la méthode d’Euler

Ln(h) Ln (Erreur d’Euler)

-2,30258509 4,4046439245
-4,6051 701 9 4,81 51 069098
-6,90775528 4,9595525039
-9,21 034037 4,9882533432

5,1

4,9

4,8

4,7
ln(err)

4,6
Ln (Erreur d’Euler)
4,5

4,4

4,3

4,2

4,1
-10 -9 -8 -7 -6 -5 -4 -3 -2 -1
ln(h)

→ Pour la Méthode d’Euler modifie

Ln(h) Ln(Erreur d’Euler modifie)

-2,302585093 1 2,9480099903
-4,6051 701 86 2,4932054526
-6,907755279 1 ,60943791 24
-9,21 0340372 1 ,2809338455
 14

12

10

8
ln(err)

6 Ln(Erreur d’Euler modifie)

0
-10 -9 -8 -7 -6 -5 -4 -3 -2 -1
ln(h)

→ Pour la Méthode krank-nicholson

Ln(h) Ln(Erreur krank_nicholson)

-2,302585093 -1 1 ,97734471 59
-4,6051 701 86 -1 2,53981 261 82
-6,907755279 -1 5,491 6226037
-9,21 0340372 -1 5,8293891 035

0
-10 -9 -8 -7 -6 -5 -4 -3 -2 -1
-2

-4

-6

-8
ln(err)

Ln(Erreur krank_nicholson)
-10

-12

-14

-16

-18
ln(h)

→ pour la Méthode de Runge_Kutta d'ordre 4

 Ln(h) Ln(Erreur de Runge_Kutta d'ordre 4)
-2,302585093 -0,4929529622
-4,6051 701 86 -3,6461 26525
-6,907755279 -2,4883244839
-9,21 0340372 -1 0,5441 561 41 5

0
-10 -9 -8 -7 -6 -5 -4 -3 -2 -1
-2

-4
ln(err)

-6 Ln(Erreur de Runge_Kutta
d'ordre 4)

-8

-10

-12
ln(h)

=> Une graphe qui rassemble tous les 4 courbes :

15

10

Ln (Erreur d’Euler)
0
Ln(Erreur d’Euler modifie)
ln(err)

-10 -9 -8 -7 -6 -5 -4 -3 -2 -1
Ln(Erreur krank_nicholson)
-5 Ln(Erreur de Runge_Kutta d'ordre
4)
-10

-15

-20
ln(h)