Wyznaczanie wartości własnych - metoda potęgowa

clear;
clc;
close all;

A = [2 2 2     %Macierz z polecenia zadania
     2 5 7
     2 7 9];
X = ones(3,1); % Tworzymy wektor początkowy [1 1 1]
k = 20; %Ilosc krokow, im wiecej tym blizej do prawdziwej wartosci

wartoscWlasna = zeros(1,k); %przygotowanie miejsca pod wektor wartości własnych
wektorWlasny = zeros(3,k);  %przygotowanie miejsca pod wektor wektorów własnych

for p=1:k
V=A*X;
m=max(abs(V)); %określanie wartości własnej jako maksimum z modułu wektora
%Macierz * wektor własny
X=V./m;  % Zapisywanie nowego wektora
wartoscWlasna(1,p)= m; %wektor wartości własnych
wektorWlasny(:,p)=X;   %wektor wektorów własnych
end