function x=direktno(L,b)
%DIREKTNO resuje spodnji trikotni sistem
%x=DIREKTNO(L,b)
%L je spodnja trikotna matrika
%b je desna stran
%x je resitev sistema

n=size(L,1);
x=zeros(n,1);
for i=1:n
  x(i)=(b(i)-L(i,1:i-1)*x(1:i-1))/L(i,i);
end