//On va résoudre un système tridiagonal
function x=solLU(n,a,b,c,y)

//On décale les coordonnées de a
a=[0;a(:)];

//Les notations "a" et "b" vont désigner "l" et "d"
for i=2:n
   a(i)=a(i)/b(i-1);
   b(i)=b(i)-a(i)*c(i-1);
end

//La notation "y" va désigner "x'"
for i=2:n
   y(i)=y(i)-a(i)*y(i-1);
end

//On initialise x et le calcule
x=zeros(n,1);
x(n)=y(n)/b(n);
for i=n-1:-1:1
   x(i)=(y(i)-c(i)*x(i+1))/b(i);
end

endfunction