% Startvektor x % Eigenvektor = [1 0 1]' % lambda = 3 x = [ 1 ; 0 ; 1; 3 ] % zu iterieren v = x(1:3) lambda = x(4) df = [ A-lambda*eye(3) -v ; 2*v' 0 ] x = x - inv(df) * [ A*v-lambda*v ; v'*v-1 ] % ende der Schleife