function quadraticMap(P0,lambda,T)
%
%  P0 = rand; lambda = 0.9; T = 100;
%  figure(1); quadraticMap(P0,lambda,T)
%
%  P0 = rand; lambda = 2; T = 100;
%  figure(1); quadraticMap(P0,lambda,T)
%
%  P0 = rand; lambda = 3.2; T = 100;
%  figure(1); quadraticMap(P0,lambda,T)
%
%  P0 = rand; lambda = 3.5; T = 100;
%  figure(1); quadraticMap(P0,lambda,T)
%
%  P0 = rand; lambda = 4; T = 10000;
%  figure(1); quadraticMap(P0,lambda,T)
%
P(1) =P0;
gP = [];
xPlot = linspace(0,1);
yPlot = lambda*xPlot.*(1-xPlot);
for t = 1:T
    P(t+1) = lambda*P(t)*(1-P(t));
    gP = [gP P(t) P(t)];
    plot(gP(1:t),gP(2:t+1),xPlot,yPlot,xPlot,xPlot,gP(t),gP(t+1),'ro');
    pause
end