function manualfourier()
 one = wavread('one.wav');
  fig=figure; 
 title('manual reconstruction');
 subplot(3,1,1); plot(one);
 [x,y,but]=ginput(1);
 
while (but==1)
   sound(one); 
   [x,y,but]=ginput(1);
end

 ff = fft(one);
 reconstructedone = zeros(size(one));
 subplot(3,1,2); plot(reconstructedone);
 for t=1:floor((size(ff,1)-1)/2)
   partialfft = bandpassfilter(ff,t,t);
   addtoreconstruction = real(ifft(partialfft));
   
   subplot(3,1,3); plot(addtoreconstruction); 
   subplot(3,1,2); hold off; plot(reconstructedone,'r');
   reconstructedone = reconstructedone + addtoreconstruction;
   subplot(3,1,2); hold on; plot(reconstructedone,'b');
   ax=axis;subplot(3,1,3); axis(ax);
   [x,y,but]=ginput(1);
 end
