x=linspace(0,1,50);
y=3.48^2*x.*(1-x).*(1-3.48*x+3.48*x.*x);
z=3.48^2*y.*(1-y).*(1-3.48*y+3.48*y.*y);
subplot(1,2,1);
plot(x,x,x,z);
title('b=3.48');
a=zeros(1,1001);
a(1)=0.5;
n=1000;
for i=1:10000
a(i+1)=3.48*a(i)*(1-a(i));
end
subplot(1,2,2);
plot(a,'.');

for i=260:400
k=i/100;
a=zeros(1,151);
a(1)=0.2;
for j=1:150
a(j+1)=k*a(j)*(1-a(j));
b=a(51:150);
c=k*ones(size(b));
plot(c,b,'.');
hold on
end;
end;
hold off