%Factorizacion A=LU sin pivoteo para una matriz tridiagonal usando 
%solamente la informacion de las diagonales.
%d es la diagonal principal
%l la diagonal inferior
%u la diagonal superior

%Las salidas son tambien l,d,u. 
%l representa la diagonal inferior de L
%d la diagonal principal de U y u la diagonal superior de U.

function [d,l,u]=factlutridiag(d,l,u)
% A matriz de NxN tri diagonal

N=length(d);

for i=2:N-1
    d(i)=d(i)-(l(i-1)/d(i-1))*u(i-1);
    l(i-1)=(l(i-1)/d(i-1));
end




%Aqui abajo el codigo si uno decide trabaja con matrices (no recomendado)

%function [L,U]=factlu3d(A)

%[N,N]=size(A);
%a=diag(A);
%b=diag(A,1);
%c=diag(A,-1);

%for i=2:N-1
%    a(i)=a(i)-(c(i-1)/a(i-1))*b(i-1);
%    c(i-1)=(c(i-1)/a(i-1));
%end
%L=diag(a)+diag(b,1);
%U=L