; 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, n, d, wL diln = 1. ; dilation posn = 0. ; position shift psihatL = sato_wL(n, d, diln, posn, 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