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