function demopd
% demonstrace metody pevneho bodu

format compact
A=input('Zadej hodnotu parametru A: ');
x0=input('Zadej pocatecni iteraci: ');
x=0:0.01:1;
f=A*x.*(1-x);
f1=figure;
plot(x,f,x,x)
hold on
plot([x0,x0],[0,x0],'r*--');
axis([0,1,0,1]);
disp('Neco stiskni')
pause
xn=x0;
ni=0;
konec=0;
klavesa=0;
while ni<100 & konec~=1;
 ni=ni+1;
 xo=xn; xn=A*xo*(1-xo)
 plot([xo,xo,xn],[xo,xn,xn],'r*--');
 pause(2)
 klavesa=get(f1,'CurrentCharacter');
 if isempty(klavesa), klavesa=0; end
 if klavesa==27, konec=1; end
% pause
end
disp('Neco stiskni')
pause
hold off
close(f1)
