function L=choleski(A)

% FATTORIZZAZIONE DI CHOLESKI.

[n,m]=size(A);

L=zeros(size(A)); L=A;

for index=1:n,
    Ldiag=L(index,index); 
    if Ldiag <= 0
        fprintf('\n \t [WARNING][LA MATRICE NON E'' DEFINITA POSITIVA')
        return
    end 
    Ldiag=sqrt(L(index,index)); 
    L(index,index)=Ldiag;
    L(index,index+1:n)=L(index,index+1:n)./Ldiag; 
    L(index+1:n,index+1:n)=L(index+1:n,index+1:n)-L(index,index+1:n)'*L(index,index+1:n); 
end    
L=triu(L); 
