r=input('r-value: 0<r<4: ');
x=0.01*(0:100);
f=r.*x.*(1.0-x);
figure('Position',[50,50,600,600])
axis([-0.1,1.0,-0.1,1.0]);
plot([-0.1,1.0],[0.0,0.0],'-k');
hold on;
plot([0.0,0.0],[-0.1,1.0],'-k');
plot(x,f,'-r','LineWidth',1.5);
x1=[0.0,1.0];
y1=x1;
plot(x1,y1,'-g','LineWidth',1.5);
if (r < 1.0)
  xb=0.0;
else 
  xb=1.0-1.0/r;
end
xb1=[xb,xb];
yb=r*xb*(1.0-xb);
yb1=[yb,yb];
plot(xb1,yb1,'ob','MarkerSize',10);
xlabel('x(n)','FontSize',12)
ylabel('x(n+1)','FontSize',12)
title('Fixed Point','FontSize',12)
hold off;