%
% Cholesky-Zerlegung einer SPD-Matrix
%
% Vorlesung "Einfuerung in die Numerik", Uni Basel, FS 11
%

function L = CholeskyDecomp(A)

n = size(A,1);
L = zeros(n);

eigenvalues = ... ; % Eigenwerte von A

if ... 
    % A ist SPD-Matrix, die Cholesky-Zerlegung existiert
    % Algorithmus 7 aus dem Skript
    for k = 1 : n
        L(k, k) = ... ;
        L(k+1:n, k) = ... ;
    end
else
    % A ist nicht positiv definit
    ... ;
end

end