f=@(x)exp(x.^2);
quadv(f,-1,1)