Taik-cours1_AN3

A. Taik Cours AN3 LST-MI
h=h’;
while(j < nt + 2)
for i=1:nx-1
w(i,j)=h(i);
end
h=M*h+N’;
j=j+1;
end
for i=nx:-1:2
for j=nt+1:-1:1
w(i,j)=w(i-1,j);
end
end
for j=1:nt+1
w(1,j)=0;
w(nx+1,j)=0;
end
mesh(t,x,w);
*********************************
Voici les courbes obtenues après compilation pour les deux cas:
Figure 2.2: Courbe de la solution numérique,
méthode explicite pour r < 1 2
Figure 2.3: Courbe de la solution numérique,
méthode explicite pour r > 1 2
Constatations:
Nous avons rémarqué, d’après les deux graphes corresondantes à r < 1 et r > 1 , que la méthode
2 2
présente deux régions de stabilité. Si r < 1 la méthode est stable et converge vers la solution
2
Département de Mathématiques
FST-Mohammedia, (2008)
26