
sigma=1;x0=0;
a=0;
b=2.33*sigma;nn=100;

x1=linspace(a,b,nn);
for i=1:nn
   ff(i)= (1/((2*pi)^.5*sigma))*exp((x1(i)-x0)^2/(2*sigma^2));
end
quad('fff',0,2.33)