[pslpsl]

n=50;
a=[2.0 3.2 3.5 3.8];
x(1)=0.5;
[p,k]=size(a);
for j=1:k
for i=1:n
x(i+1)=a(j)*(x(i)-x(i)^2);
end
end
plot(1:n+1,x)