%diary hatvanyModszer_demo.diary

clc
clear all
disp(' *** Véletlen egész mátrixokkal teszteljük a hatványmódszert:');
n = 5
maxAmax = 10
A = randi([-maxAmax, maxAmax],n)

disp(' *** Kiszámítjuk az eig-gel a sajátértékeket:');
Lambda = eig(A)

disp(' *** Milyen az A mátrix spektruma?');
disp(' *** A föntiek alapján várható-e  a hatványmódszer konvergenciája?');
input('Folytatás ENTER-rel');

maxk = 10^4
tol = 10^(-15)
[lambda1_koz, x1_koz, k] = hatvanyModszer(A, maxk, tol)

disp(' *** Elfogadható-e a hatványmódszerrel kapott eredmény?');
disp(' *** Ellenõrizzük:');
hibav = A*x1_koz - lambda1_koz*x1_koz

input('Folytatás ENTER-rel');
disp(' *** A  hatványmódszer Raleigh-hányadossal számított közelítése:');
[lambda1_kozR, x1_kozR, kR] = hatvanyModszerRaleigh(A, maxk, tol)

disp(' *** Mit tapasztalunk? Elfogadható-e a kapott eredmény?');
disp(' *** Ellenõrizzük:');
hibav = A*x1_kozR - lambda1_kozR*x1_kozR
diary off
