clear all

n_vs_r_matrix = [];

rlist=linspace(0,1,100);

figure;


for i=1:length(rlist)

    r=rlist(i);
    xlist(1)=rand;
    
    for t=1:1000
    xlist(t+1)=r*sin(pi*xlist(t)) ;  
    end

plot(r,xlist(901:1000),'k.','MarkerSize',16) ; hold on ; axis([0 1 0 1]);  
ylabel('x','FontSize',20)
xlabel('r','FontSize',20)
set(gca,'FontSize',20)
drawnow

end