function y = fun1(x)

y = 2*exp(-x.*x)/sqrt(pi);
