C» who

   fem

C» pp=postinterp(fem,'p',[0.8;0.2;0.4],'solnum',1:200)
C» pp2=[0;pp]
C» pp2(201)=real(pp(200))
C» h=fftshift(ifft(pp2,400,'symmetric'))
C» h2=h/max(abs(h))
C» wavwrite(h2,4000,32,'h2.wav')
C» dir