%Aufgabe 1 a
%
%Eingabe A = positiv definite Matrix
%Ausgabe R = Cholesky-Zerlegung (äquivalent chol in MATLAB)

function R = cholesky(A)

[m, n] = size(A);

if m == n;
R = zeros(n,n);
for i = 1:n;
R(i,i) = ( A(i,i) - sum( R(:,i) .* R(:,i) ) )^(1/2);
for j = i+1:n;
R(i,j) = ( A(i,j) - sum( R(:,i) .* R(:,j) ) ) / R(i,i);
end;
end;
else
fprintf('Matrix muss quadratisch sein.');
end;