#Métodos Numéricos - 2023. IMERL, Fing
#Código de ejemplo del ejercicio 2 del práctico 2

function x=atras(U,b)
  %resolver Ux=b para U triangular superior
  [n,n]=size(U);
  x=zeros(n,1);
  x(n)=b(n)/U(n,n);
  for k=n-1:-1:1
    j=k+1:n;
    x(k)=(b(k)-U(k,j)*x(j))/U(k,k);
  endfor

  %comparar con U\b