function y=quadfun(x)
y=exp(-x.^2);