function y = iterate(r);

xi = 0.5;
for n = 1:3000,
    xi = r*xi*(1-xi);
    x(n) = xi;
end
y = x(1001:3000)';