clc;
clear all;
close all;

len=2^19;
fs=44100;

tlen=len/fs
sig(1:len)=0;
sig(len/2)=1;


nprimes=60;
primss=list_primes(nprimes*20);
prims=primss(4:20:nprimes*20);

for ii=1:nprimes
	ii;
	ord=prims(ii);
	clear aa;
	clear bb;
	aa(1:ord+1)=0;
	bb(1:ord+1)=0;
	fac=.5;
	bb(1)=1;
	aa(ord+1)=1;
	bb(ord+1)=fac;
	aa(1)=fac;
	sig=filter(aa,bb,sig);
	sig=sig(len:-1:1);

end

sig=sig';
t=max(abs(sig));
sig=sig/t*.9;
wavwrite(sig,44100,16,'swisha.wav');
subplot(2,1,1);
plot(sig);
subplot(2,1,2);
sigt=fft(sig,2*len);
sigt=abs(sigt(1:len+1));
tmp=max(sigt);
sigt=sigt/tmp;
sigt=log10(sigt)*20;
plot(sigt);