function[iter,iRk,wektorn]=iteracjiodwrotnej
%west-wektor startowy, mian-mianownik w wektorze znormalizowanym,
%wenor-wektor znormalizowany, wenor1- wektor znormalizownay nastepny, wesn-
%wektor startowy nastepny iter -liczba iteracji, iRk-iloczyn Redleya
%(wartość wałasna maksymalny co do wartosci bezwglednej) ,wektorn- wektor znormalizownany na koncu
% epsL i epsR - szcowane błedy
A=input('Podaj macierz A: ')
west=input('Podaj wektor startowy: ')
blad=input('Podaj wartość błedu n: ')
lamk=0.0;
mian=(west'*west)^0.5;
wenor=west/mian;
epsL=1;
epsR=1;
iter=1;
while (epsL>blad)&(epsR>blad)
wesn=inv(A)*wenor;
lamk1=wenor'*wesn;
k=epsR;
mian=(wesn'*wesn)^0.5;
wenor1=wesn/mian;
epsL=abs((lamk1-lamk)/lamk1);
epsR=((wenor1-wenor)'*(wenor1-wenor))^0.5;
iter=iter+1;
if (iter>200)
printf('rozbiezny\n');
return
end %if
lamk=lamk1;
wenor=wenor1;
end %while
iRk=lamk^(-1);
wektorn=wenor;
end %function




function[iter,iRk,wektorn]=potegowa
%west-wektor startowy, mian-mianownik w wektorze znormalizowanym,
%wenor-wektor znormalizowany, wenor1- wektor znormalizownay nastepny, wesn-
%wektor startowy nastepny iter -liczba iteracji, iRk-iloczyn Redleya
%(wartość wałasna maksymalny co do wartosci bezwglednej) ,wektorn- wektor znormalizownany na koncu
% epsL i epsR - szcowane błedy
A=input('Podaj macierz A: ')
west=input('Podaj wektor startowy: ')
blad=input('Podaj wartość błedu n: ')
lamk=0.0;
mian=(west'*west)^0.5;
wenor=west/mian;
epsL=1;
epsR=1;
iter=1;
while (epsL>blad)&(epsR>blad)
wesn=A*wenor;
lamk1=wenor'*wesn;
k=epsR;
mian=(wesn'*wesn)^0.5;
wenor1=wesn/mian;
epsL=abs((lamk1-lamk)/lamk1);
epsR=((wenor1-wenor)'*(wenor1-wenor))^0.5;
iter=iter+1;
if (iter>200)
printf('rozbiezny\n');
return
end %if
lamk=lamk1;
wenor=wenor1;
end %while
iRk=lamk^(-1);
wektorn=wenor;
end %function