clear('x')
r=0.9;
x(1)=0.1;
for n=1:29
    x(n+1)=r*x(n)*(1-x(n))
end
plot([1:1:30],x,'-o')
