;  Given a number of points and a spacing, return an
;  array of the wavelet fourier series and the
;  array of frequencies (2 pi/wavelength) for the
;  mother wavelet as defined in Sato.

function sato_mother_wL_g, w
; Function for use by refr_mother_wL
; Only defined for 2 pi/3 < w < 8 pi/3

  c = 2.*!pi/3.
  h1 = exp(-(2.*c - w)^(-2.))
  h2 = exp(-(w-c)^(-2.))
  return, h1/(h1+h2)
end


function sato_mother_wL, n, d, wL

  psihatL = sato_wL(n, d, 1, 0, wL)
  
  return, psihatL

end

pro sato_mother_wL_test

  n = 256
  d = 20./n
  psihatL = sato_mother_wL(n, d, wL)
  window, 0, xs=1000,ys=500
  !p.multi=[0,2,0]
  plot, wL, abs(psihatL), xr=[-10,10], yr=[0,1.2], /ys

end