188
11
•
Problèmes
-13.9541
55.4278
-198.1569
✝
✆
3. Construire une fonction u0.m qui permet l’initialisation.T est u0((0 :0.2 :1)’.
✞
☎
>> u0((0:0.2:1)’)
ans =
0
0.6198
1.0471
1.0951
0.7158
0.0000
✝
✆
4. Compléterlepremier programme pour calculer la solution approchée V en utilisant l’instruction ode23 puis dessiner les solutions des problèmes exact et approché en dimension 3. Solution
exacte : u ( x , t ) = e
− p
2 t
sin p x + x
2
(1 − x ). On tracera aussi les lignes de niveaux. Utiliser
meshgrid,s urf, contour.Sauvegarder ELF2.m.T est pour N = 4, afficher V (3, :), puis graphique avec N = 20.
N = 4
✞
☎
>> V(3,:)
ans =
00 .5256
0.8882
0.9334
0.6073
0
✝
✆
N = 20
11
•
Problèmes
-13.9541
55.4278
-198.1569
✝
✆
3. Construire une fonction u0.m qui permet l’initialisation.T est u0((0 :0.2 :1)’.
✞
☎
>> u0((0:0.2:1)’)
ans =
0
0.6198
1.0471
1.0951
0.7158
0.0000
✝
✆
4. Compléterlepremier programme pour calculer la solution approchée V en utilisant l’instruction ode23 puis dessiner les solutions des problèmes exact et approché en dimension 3. Solution
exacte : u ( x , t ) = e
− p
2 t
sin p x + x
2
(1 − x ). On tracera aussi les lignes de niveaux. Utiliser
meshgrid,s urf, contour.Sauvegarder ELF2.m.T est pour N = 4, afficher V (3, :), puis graphique avec N = 20.
N = 4
✞
☎
>> V(3,:)
ans =
00 .5256
0.8882
0.9334
0.6073
0
✝
✆
N = 20