Solution de la série n°05

Exercice n°01 :
= +
0 =1

1. Calcul de la solution approximative de l'équation différentielle en x = 1 à l'aide de la méthode

d'Euler

L'intervalle de travail est [0,1].

Le pas h est ℎ = = 0.1

De plus on a , = +

On peut donc utiliser la méthode d’Euler :

≅ +ℎ ,

Et obtenir successivement des approximations de :

0.1 , 0.2 , 0.3 , 0.4 , 0.5 , 0.6 , 0.7 , 0.8 , 0.9 , 1.0 .

La première itération produit :

0.1 ≅ 0 + ℎ 0,1 = 1 + 0.1 1 + 0 = 1.1

De manière similaire, la deuxième itération donne :

0.2 ≅ 0.1 + ℎ 0.1,1.1 = 1.1 + 0.1 1.1 + 0.1 = 1.22

Et la troisième itération donne :

0.3 ≅ 0.2 + ℎ 0.2,1.22 = 1.22 + 0.1 1.22 + 0.2 = 1.362

Le tableau ci-dessous rassemble les valeurs des dix premières itérations :

0.0 0.1 0.2 0.3 0.4 0.5 0.6

1.0000 1.1000 1.2200 1.3620 1.5282 1.72102 1.943122

0.7 0.8 0.9 1.0

2.1974342 2.48717762 2.815895382 3.1874849202

Donc, la solution approximative de l’équation différentielle en = 1:

 1 ≅ 3.1874849202

2. Comparaison du résultat obtenu avec la solution exacte

La solution exacte est: !"# = −1 − + 2% . Ainsi, la solution exacte en & = 1 donne :

!"# = −1 − + 2% = 3.4365636569

Donc l'erreur commise est :

' =| !"# 1 − 1 | = |3.4365636569 − 3.1874849202| = 0.2490787367

Exercice 2 :
= + %)
0 =2

1. Calcul de la première itération de la méthode d'Euler modifiée

Le pas est h = +. ,.
De plus, on a f x, y = y + e)2
On peut donc utiliser la méthode d'Euler modifiée :
ℎ
= + 3 , + , 4
2
Avec :
= +ℎ ,
Et obtenir l’approximation de 0.1 . La première itération produit :
= 0 + ℎ 0,2 = 2 + 0.1 2 + % )∙ = 2.3

qui est le résultat obtenu à l'aide de la méthode d'Euler. La deuxième étape donne :

ℎ 0.1
0.1 ≅ 0 + 6 0,2 + 0.1,2.3 7 = 2 + 62 + % )∙ + 2.3 + % )∙ . 7 = 2.326070138
2 2

Erreur commise
la solution exacte !"# = % + %) .

La solution exacte en x = 0.1 donne : !"# 0.1 = % .

+ % )∙ .
= 2.326573676, donc l'erreur commise
est :

' " 89 . 9|:;<=>?; . @ . | = |2.326573676 − 2.326070138| = 0.000503538

2 . Calcul des deux premières itérations de la méthode d'Euler modifiée

Le pas est h = 0.05
De plus, on a &, A = A + % )
On peut donc utiliser la méthode d'Euler modifiée et obtenir successivement les approximations de y(0.05) et
y(0.1). La première itération produit :
 = 0 +ℎ 0,2 = 2 + 0.05 2 + % )∙ = 2.15

qui est le résultat obtenu à l'aide de la méthode d'Euler. La deuxième étape donne :

8 . B
0.05 ≅ 0 + 6 0,05 + 0.05,2.15 7 = 2 + 62 + % )∙ + 2.15 + % )∙ . B7
= 2.156379273
) )

De manière similaire, la deuxième itération donne :

= 0.05 + ℎ 0.05,2.156379273 = 2.156379273 + 0.05 2.156379273 + % )∙ . B

= 2.3194567

La correction conduit à son tour à :

8
0.1 ≅ 0.05 + 6 0.05,2.156379273 + 0.1,2.319456782 7 = 2.156379273 +
)
. B
62.156379273 + % )∙ . B
+ 2.319456782 + % )∙ . 7 = 2.326439516
)

Erreur commise

L’erreur commise est :

' " 89 . B =| !"# 0.1 − C 0.1 | = |2.326573676 − 2.326439516| = 0.00013416

3. Rapport des erreurs :

Le rapport des erreurs est :

E;> FGH.I
D= = 3.75 ≅ 2)
E;> FGH.HJ

On voie nettement que l'écart entre la solution exacte et la solution approximative diminue d'un facteur de
3.75 ≅ 2) lorsque le pas h est divisé par 2, ce qui confirme que la méthode d'Euler modifiée est d'ordre 2.

Exercice n°03 :
=− + +1
0 =1

1. Calcul de l'approximation de 0,2 :

L'intervalle de travail est [0,0.2],

Le pas h e s t : h = 0.1.
De plus on a : , =− + + 1.
a. Méthode d'Euler
On peut donc utiliser la méthode d'Euler :
≅ +ℎ ,
et obtenir successivement des approximations de y 0.1 et y 0.2 . La première itération produit :
y 0.1 ≅ y 0 + h 0.1 = 1 + 0.1 −1 + 0 + 1 = 1.

De manière similaire, la deuxième itération donne :

y 0.2 ≅ y 0.1 + h 0.1.1 = 1 + 0.1 −1 + 0.1 + 1 = 1.01
Donc l'approximation de y(0.2) en utilisant la méthode d'Euler esty 0.2 ≅ 1.01.

b. Méthode de Runge-Kutta d'ordre 4

On peut donc utiliser la méthode de Runge-Kutta d'ordre 4:
k = hf ,
ℎ k
k ) = hf L + , + M
2 2
ℎ k)
k N = hf L + , + M
2 2
k O = hf + ℎ, + kN

1
≅ + k + 2k ) + 2k N + k O
6

et obtenir successivement des approximations de 0.1 et 0.2 . La première itération produit :

k = hf 0,1 = 0.1 −1 + 0 + 1 = 0
0.1 0
k ) = hf L0 + , 1 + M = 0.1 0.05,1 = 0.1 −1 + 0.05 + 1 = 0.005
2 2
0.1 0.005
k N = hf L0 + ,1 + M = 0.1 0.05,1.0025 = 0.1 −1.0025 + 0.05 + 1 = 0.00475
2 2
k O = hf 0 + 0.1,1 + 0.00475 = 0.1 0.1,1.00475 = 0.1 −1.00475 + 0.1 + 1 = 0.009525

Ce qui entraine que :

0.1 ≅ 0 + P k + 2k ) + 2k N + k O = 1 + P 0 + 2 ∙ 0.005 + 2 ∙ 0.00475 + 0.009525 = 1.00483

De manière similaire, la deuxième itération donne :

k = hf 0.1,1.0048375 = 0.1 −1.0048375 + 0.1 + 1 = 0.00951625

. . RB P)B
k ) = hf Q0.1 + ) , 1.0048375 + S = 0.1 0.15,1.00959525 = 0.1 −1.00959525 + 0.15 +
)
1 = 0.014040437
. . O O ONT
k N = hf Q0.1 + ) , 1.0048375 + )
S = 0.1 0.15,1.011857718 = 0.1 −1.011857718 +
0.15 + 1 = 0.013814228

k O = hf 0.1 + 0.1,1.0048375 + 0.013814228 = 0.1 0.2,1.018651728 = 0.1 −1.018651728 +

0.2 + 1 = 0.018134827

Ce qui entraine que :

0.2 ≅ 0.1 + P k + 2k ) + 2k N + k O = 1.0048375 + P 0.00951625 + 2 ∙ 0.014040437 + 2 ∙

0.013814228 + 0.018134827 = 1.018730901

Donc l'approximation de y(0.2) en utilisant la méthode de Runge-Kutta d'ordre 4 est

0.2 ≅ 1.018730901
 2. Erreur commise

La solution exacte est : !"# 0.2 ≅ 1.018730780

a. Erreur commise par la méthode d'Euler :

EV = | !"# 0.2 − 0.2 | = |1.018730780 − 1.01| = 0.008730780 = 0.873078 ∙ 10 )

b. Erreur commise par la méthode de Runge-Kutta d'ordre 4 :

EWXO = | !"# 0.2 − 0.2 | = |1.01873078 − 1.018730901| = 0.000000121 = 0.121 ∙ 10 P

Il est clair que la méthode d'Euler est très imprécise par rapport à celle de Runge-Kutta d'ordre 4. Il est donc
préférable d'utiliser des méthodes d'ordre aussi élevé que possible.

Exercice n04 :
Étant donné,

Y /Y = )
+ )
, y(1) = 1.2

Donc, , = )
+ )

= 1 et = 1,2

De même, h = 0,05

Calculons les valeurs de [ , [) , [N et [O .

)7
[ =ℎ , = 0.056 )
+ = 0.056 1 )
+ 1.2 ) 7 = 0.0561 + 1.447 = 0.0562.447 = 0.122

[) = ℎ + ℎ⁄2 , + [ ⁄2 = 0.05 1 + 0.025,1.2 + 0.061

{puisque h/2 = 0,05/2 = 0,025 et [ /2 = 0,122/2 = 0,061}

= 0.056 1.025,1.261 7 = 0.056 1.025 )

+ 1.261 ) 7 = 0.0561.051 + 1.5907 = 0.0562.6417 = 0.1320

[N = ℎ + ℎ⁄2 , + [) ⁄2 = 0.05 1 + 0.025,1.2 + 0.066

{puisque h/2 = 0,05/2 = 0,025 et [) /2 = 0,132/2 = 0,066}

= 0.056 1.025,1.266 7 = 0.056 1.025 ) + 1.266 ) 7 = 0.0561.051 + 1.6027 = 0.0562.6537 = 0.1326

[O = ℎ + ℎ, + [N = 0.05 1 + 0.05,1.2 + 0.1326 = 0.05 1.05,1.3326 = 0.056 1.05 )

+
)7
1.3366 = 0.0561.1025 + 1.77587 = 0.0562.87837 = 0.1439

Par la méthode RK4, on a ;

= + P [ + 2[) + 2[N + [O
= 1.05 = + P [ + 2[) + 2[N + [O

En substituant les valeurs de , [ , [) , [N et [O , on obtient ;

1.05 = 1.2 + P 0.122 + 2 0.1320 + 2 0.1326 + 0.1439 = 1.2 + P 0.122 + 0.264 + 0.2652 +
0.1439 = 1.2 + P 0.7951 = 1.2 + 0.1325 = 1.3325