function []=cholesky(A,b)
clc;
clear all;
	A=matrizA;
	b=vectorb;

	Aumentada=[A,b]
	[n,m]=size(A);

for k=1:n
    Acumulador = 0;
    for p=1:k-1
        Acumulador = Acumulador + L(k,p)*U(p,k);
    end
    L(k,k)= sqrt(A(k,k)- Acumulador);
    U(k,k) = L(k,k);
    
    for i=k:n
        Acumulador2 = 0;
        for p=1:k-1
            Acumulador2= Acumulador2 + L(i,p)*U(p,k);
        end
        L(i,k) = (A(i,k)- Acumulador2)/L(k,k);
    end
    for j= k:n
        Acumulador3 = 0;
        for p = 1:k-1
            Acumulador3 = Acumulador3 + L(k,p)*U(p,j);
        end
        U(k,j)= (A(k,j)- Acumulador3)/L(k,k);
    end
    
    
end
disp ('L=')
disp(L)
disp('U=')
disp (U)
[Lb]= Aumentada(L,b);
for i=1:n
    Sumatoria=0;
    for p=1:i-1
        Sumatoria=Sumatoria+Lb(i,p)*z(p,1);
    end
    z(i,1)=(Lb(i,n+1)-Sumatoria)/Lb(i,i);
end

[Uz]=AB(U,z);

function [X]=SustitucionRegresiva(UB,n)
for i=n:-1:1
    Sumatoria=0;
    for p=i+1:n
        Sumatoria=Sumatoria+UB(i,p)*X(p,1);
    end
    X(i,1)=(UB(i,n+1)-Sumatoria)/UB(i,i);
end

disp('X=')
disp(X)
end
