clear
N=400;
x(1)=0.4;
k=2.1
for n=2:N,
    x(n)=k*x(n-1)*(1-x(n-1));
end
plot(x,'b')
