function S=solequadiff(n)

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

// coefficients pour le systeme lineaire AX=Y
b=2+k/(n*n);
a=(-1)*ones(n-1,1);
c=a;
y=f/(n*n); y(1)=y(1)+alpha; y(n-1)=y(n-1)+beta;

// resolution du systeme lineaire
getf('solLU.sci')
u=solLU(n-1,a,b,c,y);
v=(1-x).*exp(x);
S=[x,u,v];