function[R,s]=choleskyredd(B)
A = [60 30 20; 30 20 15; 20 15 12]

condicion=true;
for i=1:3;
    for j=1:3;
    if (A(i,j) ~= A(j,i) )
        condicion=false;
        i=3;
        break
        end
    end
end
if (condicion)
disp ("la matriz es simétrica")
else
disp ("la matriz NO es simétrica")
end


e = eig(A)
condicion=true;
for i=1:3;
    for j=1:1;
    if (e(i,j)<0)
        condicion=false;
        i=3;
        break
         end
    end
end
if (condicion)
    disp("La matriz es definida positiva")
    disp ("Admite factorizacion de Cholesky")
    B=A;
else
    disp("la matriz no es definida positiva")

end
s=0;
R=B;
[m,n]=size(B);
for k=1:m
    s=s+1;
    for j=k+1:m
        R(j,j:m)=R(j,j:m)-(R(k,k:m)*R(k,j)/R(k,k));
    end
    R(k,k:m)=R(k,k:m)/sqrt(R(k,k));
end
R=tri(R);
end