C» who

   fem

C» pp=postinterp(fem,'p',[1.85;-0.37;0.73],'solnum',1:200)
C» pp2=[0;pp]
C» h=ifft(pp2,400,'symmetric')
C» h2=h/max(h)
C» wavwrite(h2,800,'pippo.wav')
C» dir