global tref wref
tref= -500:1:500;
wref= wavelet_ft(tref);
plot(tref,wref,'b-')
save wavelet.mat tref wref