samp=2^16;
%ini=3.95;
fin=4;
ini=2.5;
%fin=3.5699456;

fstp=(fin-ini)/samp;
a=zeros(samp,1);
x=1/10;
r=ini;
b=rand(size(a));
for i=(2:samp)
	x=r*x*(1-x);
	r=r+fstp;
	b(i)=x;
end
%sound(b,2000);
%player=audioplayer(b,2000);
%play(player);
wavwrite(b,2000,'logistic1');
c=fft(b);
plot(imag(c));
hold on
x=linspace(1,length(b),length(b));
scatter(x,max(c)/10*b/max(b),'r','.');
plot(imag(c));
hold off
%mymat=b';
%save cwtmat mymat