se=zeros(1, 100);
av=linspace(0.1, 1, 100)%1/4
for i=1:100
%a=av(i);
a=0.8;%0.5;%0.52;%0.45;%1/4;
%0.5: +- -ban 1 hullám szõrösödik ez 2 határ marad meg legtovább, szimetria 
b=1-2*a;
meret=20;

PFmx=(diag(b*ones(meret,1),0)+diag(a*ones(meret-1,1),1)+diag(a*ones(meret-1,1),-1))^i; %Perron-Flobenius mátrix elõállítás
[svektorok, sajatertekek]=eig(PFmx); %sajátérték számítás
sajatertekek=eig(PFmx)
se(i)=min(sajatertekek);

end

%{
se
plot(av, se);
%svektorok(:,17:20)
%plot(svektorok(:,(meret-4):meret))
legend({});
%xlim([1,50]);
%}