function y = kernel(x,t)

     y = exp(i*(x.^2/t - pi)/4 )/sqrt(4*pi*abs(t) );
