function S=solequadiff(n)

% vecteurs pour le problème aux limites
x=[0:1/n:1]; x=x(2:n)'; % valeurs utiles de l'abscisse
k=x;                    % valeurs utiles de la fonction k
f=(1+2*x-x.*x).*exp(x); % valeurs utiles de la fonction f
alpha=1; beta=0;

% coefficients pour le système linéaire AX=Y
b=2+k/(n*n);
a=-ones(n-2,1);
c=a;
y=f/(n*n);
y(1)=y(1)+alpha;
y(n-1)=y(n-1)+beta;

% résolution du système linéaire
u=solLU(n-1,a,b,c,y);
v=(1-x).*exp(x);
S=[x,u,v];