function [x,g,g2,g3,g4]=bosn(alpha);

x=zeros(201,1);
g=x; 
g2=x; 
g3=x; 
g4=x;

for i=1:201;
x(i)=(i-1)*0.005;
g(i)=alpha*x(i)*(1-x(i));
g2(i)=alpha*g(i)*(1-g(i));
g3(i)=alpha*g2(i)*(1-g2(i));
g4(i)=alpha*g3(i)*(1-g3(i));
g5(i)=alpha*g4(i)*(1-g4(i));
g6(i)=alpha*g5(i)*(1-g5(i));
g7(i)=alpha*g6(i)*(1-g6(i));
g8(i)=alpha*g7(i)*(1-g7(i));
g9(i)=alpha*g8(i)*(1-g8(i));
g10(i)=alpha*g9(i)*(1-g9(i));
end

plot(x,x,x,g,x,g2,x,g3,x,g4);
title('x, g(x), g^2(x), g^3(x), g^4(x),')

%"Complex and Chaotic Nonlinear Dynamics. 
% Advances in Economics and Finance, 
% Mathematics and Statistics" 
% T.Vialar, Springer 2009.
% Copyright(c).