équation de convection-diffusion stationnaire

Bonjour, j'ai essayé de résoudre une équation de convection-diffusion par la méthode des différences finis, je me suis trouvée en fin à résoudre un système d'équation linéaire A c = B (c: vecteur concentration que je dois trouvée, et B vecteur qui présente les conditions aux limites), j'ai essayé de programmer ça avec MATLAB (c'est mon premier pas avec MATLAB alors je m'érite comeme des encouragements ) mais j'ai un problème qelq part, je ne sais pas ou
voila ce que j'ai fais:

% résolution d'une équation de convection-diffusion stationnaire
L = 1 ; %longeur du domaine
N = 5 ; %nombre d'itiration
deltax = L / N; %pas d'espace
n = N+L ; %nombre de noeds
D = 0.01; %diffusion (m2/s)
C0 = 3.5e-3; %concentration initiale (g/l)
C1 = 3.e-3; %concentration finale (g/l)
u = 0.91; %vitesse (m/s)
alfa = (D + u/(deltax/2));
beta = (D - u/(deltax/2));
% remplissage de la matrice A
A = [n,n];
for i = 1:(n);
    A(i,i) = -2 * D;
end
for i = 1:(n-1);
    A(i,i+1) = beta;
    A(i+1,i) = alfa,
end
%vecteur B
B = [1,n];
for i = 1;
    B(i) = (- alfa * C0);
end
    for i = n;
    B(i) = (- beta * C1);
    end
    for i = 2 : (n-1);
    B(i) = 0;
    end
%résolution
vecterc = [1,n];
for i = 1 : n;
    vecterc = A\B;
end
c = vecterc;
plot(c)
    end