function G=Chols_1(A)
% Factorización de Cholesky
%
 n=size(A,1);
 for i=1:n
    j=i+1:n; A(i,i)=sqrt(A(i,i));
             A(i,j)=A(i,j)/A(i,i);
             A(j,j)=A(j,j)-A(i,j)'*A(i,j);
 end
 G=triu(A);
end
