function Lambda = GersgorinLambdaMax(A)

n = size(A,1);
for i=1:n
    Lambda(i) = A(i,i) + sum(abs(A(i,:))) ...
        -abs(A(i,i));
end
Lambda = max(Lambda);