function [l,u,z] = my_tridiag_lu_03_04(a,b,c)
%calcolo della fattorizzazione LU di matrici tridiagonali
%con l'algoritmo di Thomas
%tridiag(b,a,c)=L(l,1)*U(u,z)
n=length(a);
n1=n-1;
l=zeros(n1,1);
u=zeros(n,1);
z=zeros(n1,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