function out = gyaku1(A,alpha)
% Compute Eigenvalue by Inverse Iteration Method called by loop
%               H. Hasegawa;  July 2, 2006
s=size(A); I=eye(s(1)); xx=rand(s(1),1); k=0; res=1;
while ( res > 1.0e-12)
   k = k+1;
   x = xx/norm(xx,2);
   xx = (A-alpha*I) \ x;
   mu = x'*xx;
   res = norm(xx-mu*x,2);
   lambda = alpha+1/mu;
   out = [k lambda res]; % No output
end
 

   