function s = logstep(r,x,n)

a=zeros(n);
b=zeros(n);

for i=1:n
     x = r * x * (1-x);
    a(i)=i;
    b(i)=x;
end

z=[a,b];
plot(z);
s=x;