function x=tridiagsolve(a,b,c,y);

% Resi sistem s tridiagonalno matriko z glavno diagonalo a, poddiagonalo b
% in naddidagolano c ter desno stranjo y
% Za resevanje se uporablja LU razcep brez pivotiranja

% Bor Plestenjak 29.2.2004

n=length(a);
u=zeros(n,1);
l=zeros(n-1,1);
x=zeros(n,1);

%% LU razcep
u(1)=a(1);
for j=1:n-1
    l(j)=c(j)/u(j);
    u(j+1)=a(j+1)-l(j)*b(j);
end

%% resimo sp. trikotni sistem
z(1)=y(1);
for j=1:n-1
    z(j+1)=y(j+1)-l(j)*z(j);
end

%% resimo zg. trikotni sistem
x(n)=z(n)/u(n);
for j=n-1:-1:1
    x(j)=(z(j)-b(j)*x(j+1))/u(j);
end




    