function Mar25()
% The eigenvalues of A = inv(B - mu*I) are
%   1/(lambda(1)-mu),...,1/(lambda(n)-mu))
% Thus, if mu is closest to lambda(k), then when we apply the power method
% to A, it will try to compute 1/(lambda(k)-mu).