#Métodos Numéricos - 2023. IMERL, Fing
#Función para el ejercicio 8 del práctico 2
function x=tridiagonal(A,b)
  [n, m] = size(A);
  L = eye(n);
  U = A;

  for i=2:n
     L(i,i-1)=U(i,i-1)/U(i-1,i-1);
     U(i,i-1)=0;
     U(i,i)=U(i,i)-L(i,i-1)*U(i-1,i);
     if (i<n)
        U(i,i+1)=U(i,i+1)-L(i,i-1)*U(i-1,i+1);
     endif
  endfor

  #Resuelvo usando las funciones que ya implementamos.
  y=adelante(L,b);
  x=atras(U,y);

  #Pregunta: cómo se podría resolver de manera eficiente?
  #Respuesta: ver thomas.m
endfunction