function sato_jk_xL, n, d, j, k, xL

;  Given a number of points and a spacing, return an
;  array of the wavelet values and the
;  array of spatial positions for the
;  j, k wavelet as defined by Sato.

  dn = d * n
  xL=[findgen(n/2)-n/2,findgen(n/2)] * d
  
  psihatL=sato_jk_wL(n, d, j, k, wL)
  psiL=shift(float(fft( shift(psihatL,n/2)/dn, /inverse)),n/2)

  return, psiL
  
end

pro sato_jk_xL_test

  n = 256
  d = 20./n
  j=1
  k=1
  psihatL = sato_mother_wL(n, d, wL)
  psiL = sato_mother_xL(n, d, xL)
  psijkhatL = sato_jk_wL(n, d, j, k, wL)
  psijkL = sato_jk_xL(n, d, j, k, xL)
  window, 0, xs=1000,ys=500
  !p.multi=[0,2,0]
  plot, xL, psiL, xr=[-5,5], yr=4*[-1.5, 1.5], /xs,/ys
  oplot, xL, psijkL, li=2

  plot, wL, abs(psihatL), xr=[-20,20], yr=[0,1.2], /xs,/ys
  oplot, wL, abs(psijkhatL), /li

end

pro sato_jk_xL_test2

  n = 1024
  d = 20./n
  k=0
  for j=0, 4 do begin
    psijkL = sato_jk_xL(n, d, j, k, xL)
    print, j, max(psijkL), max(psijkL) * 2.^(-j/2.)
  end

end