for k = 1:n
    y = A(k,k) - sum(L(k, 1:k-1).*L(k, 1:k-1));
    if(y <= 0)
        disp('Blad, macierz nie jest dodatnio okreslona');
        return;
    end
    L(k, k) = sqrt(y);
    for j = (k+1):n
        L(j, k) = (A(j, k) - ( sum(L(j, 1:k-1) .* L(k, 1:k-1)) ))/L(k,k);
    end
end