function [x,y,z]=gdttest(a,t0,si)
%wavelet transform Mexscro Hat
t=(1:1:si);    %給出時間序列
x=(1-((t-t0)/a).^2).*exp(-((t-t0)/a).^2/2)/a;   %二階導數
y=exp(-((t-t0)/a).^2/2)/a;     %高斯函數
z=(t-t0).*exp(-((t-t0)/a).^2/2)/a;   %一階導數