function X = rh7(A,B)
    [m,n] = size(A);
    X = A\(eye(n)+B)/A;
end
