asdasdasdqwdas

    R = R(1:n, 1:n);

    % LZNK
    d = y;
    # x = inv(R) * d, gdzie d = Q^T * y, ale Q^T = Q = H_1 * H_2 * ...
    for k = fliplr(n:1)
        v = B(:,k);
        % uwaga: v juz sa znormalizowane
        H = eye(n,n) - (2 .* v * v');
        d = H * d;
    endfor
    x = inv(R) * d(1:n,:);