PRO fourier, a, o, psname
psopen, psname,/landscape,/bold,/isolatin1,/color

!p.multi=[0,1,2]
r=read_wav(a)
q=fft(r)
n=n_elements(q)
z=shift(x,n/2.)
zip=n_elements(z)
c=float(zip)
abs=abs(z)
p=(abs)^2

freq=44100

nmin=-freq/2.
nmax=freq/2.
bin=freq/c


tum=((nmax-nmin)/bin)+1
num=findgen(tum)
kum=num*bin
sum=kum+nmin


xaxis=sum
plot, xaxis, p, $
title='Power Spectrum Of A Tuning Fork', xtitle='Frequency [Hz]',$
 ytitle=textoidl('decibels^2',font=0), thick=2, charthick=3,font=0, xrange=[-5000, 5000], yrange=[0,15000]


fr=read_wav(o)
xr=fft(fr)
nr=n_elements(xr)
zr=shift(xr,nr/2.)
zipr=n_elements(zr)
cr=float(zipr)
absr=abs(zr)
pr=(abs)^2


plot, xaxis, pr, $
title='Power Spectrum Of A Tuning Fork and Harmoinca', xtitle='Frequency [Hz]',$
ytitle=textoidl('decibels^2',font=0), thick=2, charthick=3,font=0, xrange=[-5000, 5000], yrange=[0, 15000]
psclose
end