function [lambdamax , z] = Mises (A , nmax)

if nargin < 2
    maxit = 100 ;
else
    maxit = nmax;       % Maximale Anzahl der Iterationen
end
[n , n] = size (A) ;

eps = 1e-8 ;            % Genauigkeitsforderung
z   = ones (n,1) ;      % Startvektor
lam = 0 ;

for i=1:maxit
    zny = z/sqrt(z'*z) ;    % Normierung
    z   = A*zny ;           % Iterierter Vektor
    lambdamax = (zny'*z) ;  % Rayleigh-Quotient
    if abs((lambdamax-lam)/(lambdamax+lam)) < eps break; end
    lam = lambdamax ;
end

z = z/sqrt(z'*z) ;          % Normierung des Ergebnisvektors