Bonjour, j'ai un souci je n'arrive pas à simuler l'équation adimensionnée
Voila mon programme :
clear all; close all; format compact;
Le = 1;
Z = 5; %Nombre de Zel'dovich
nx=101;
nt=101;
h=1/(nx-1);
k=1/(nt-1);
T=zeros(nx,nt);
omega=zeros(nx,nt);
for n=1:nt
T(1,n)=0;
T(nx,n)=1;
end %Conditions aux limites
for n=1:nt
for i=2:nx-1
omega(i,n) = Z^2/2/Le*exp(Z*(T(i,n)-1));
T(i,n+1)= T(i,n) + k/h^2*(T(i-1,n)-2*T(i,n)+T(i+1,n))+k*omega(i,n );
end
end
T
Je dois à la fin obtenir un graphe de T fonction de x mais je n'arrive pas à obtenir T.
Pouvez-vous me donner un coup de pouce?
Merci
-----