function [l,u,z] = my_tridiag_lu_31_03(a,b,c)
%calcola la fattorizzazione LU di una matrice tridiagonale
%definita da tre vettori a,b,c  i cui elementi rappresentano
% la diagonale principale,la sottodiagonale e la sopradiagonale
%rispettivamente ponendo il risultato in l, u, z .....
n=length(a);
n1=n-1;
l=zeros(n1,1);z=zeros(n1,1);
u=zeros(n,1);
u(1)=a(1);
for k=2:n
    z(k-1)=c(k-1);
    l(k-1)=b(k-1)/u(k-1);
    u(k)=a(k)-l(k-1)*z(k-1);
end
end