// modelo poblacional

a = input('alfa =  ')
x=[];
x(1) = 0.1;
for k = 1:99;
  x(k+1) = a*x(k).*(1-x(k));
end

clf()
plot2d(x)
xlabel('k')
ylabel('x(k)')
