function y = intgnd2(x,t) global x0 y = exp(-3*x.^2).*exp(-i*t*(x-x0).^2);