Problème diffusion matlab

Bonjour,
Je souhaiterais l'équation de diffusion de particules générale en une dimension et en fonction du temps. Physiquement, à l'instant t=0, l'échantillon est sous une concentration C0 à gauche, à alpha*C0 dans l'échantillon et 0 à droite.
Voila le programme que j'ai écrit :

clear all

%Entrée valeurs opérateurs%
C0=input('Concentration solution amont (en Bq/m^3) = ');
L0=input('Epaisseur échantillon (en mm) = ');
S=input('Surface (en m²) = ');
Eps=input('Taux de porosité accessible = ');
De=input('Coefficient de diffusion effectif (en m²/s) = ');
T=input('Horizon temporel (en jours)= ');
Al=input('Valeur de alpha = ');

%Definition des constantes%
L=L0/1000;
N=100;
h=L/N;
x=(0:h:L)';

n=100*T; To=T/n;
t=(0:To:T)';

K=(De*To)/(Eps*h*h);

%Conditions initiale et aux limites et Matrice%
C(1,1)=C0; C(2:N,1)=Al*C0; C(N+1,1)=0;

A=zeros(N+1,N+1); A(1,1)=1; A(2,1)=K; A(N+1,N+1)=1;
A(2:N,2:N)=diag(-(2*K+1)*ones(N-1,1))+diag(K*ones(N-2,1),1)+diag(K*ones(N-2,1),-1);


%Algorithme%
for i=1:n
C(:,i+1)=A\C(:,i);
end


%Graphes%
surf(t,x,C)

Mais celui-ci ne semble pas fonctionner (sous matlab, le programme tourne depuis plus d'une heure sans m'afficher le moindre résultat)