fl=1;
ini=0;
fin=1.42;
wid=fin-ini;
winw=512;
winh=384;
frames=20000;
xx=zeros(frames,1);
yy=zeros(frames,1);
stp=1;
b=0.3;
for i=(1:winw)
	r=ini+i*wid/winw;
	x=1/10;
	y=0.1;
	for j=(1:100)
		t=1+y-r*x^2;
		y=b*x;
		x=t;
	end
	win=zeros(winh);
	for j=(1:100)
		t=1+y-r*x^2;
		y=b*x;
		x=t;
		win(ceil(winh*(x+1.5)/3))=1;
	end
	for j=(1:winh);
		if win(j)==1 
			xx(stp)=i;
			yy(stp)=j;
			stp=stp+1;
		end
	end
end
%plot(xx,yy,'b.');
plot(xx,yy,'ok','MarkerSize',1);
samp=2000;
while (fl==1)
[x,y]=ginput(1);
par=ini+x*wid/winw;
a=zeros(samp,1);
x=1/10;
y=0.1;
	
r=par;
for i=(1:100)
		t=1+y-r*x^2;
		y=b*x;
		x=t;
end
for i=(1:samp)
		t=1+y-r*x^2;
		y=b*x;
		x=t;
	a(i)=x;
end
sound(a,samp);
end;