function x=obratno(U,b)
%OBRATNO resuje spodnji trikotni sistem
%x=OBRATNO(U,b)
%U je zgornja trikotna matrika
%b je desna stran
%x je resitev sistema

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