clc
clear all
f=inline('(1/sqrt(2*pi))*exp((-x^2)/2)')
a=-1;
b=1;
n=4;
h=(b-a)/n;
x0=a;
%para la suma 1
sum1=f(x0);
%para la suma 2
sum2=0;
x0=a+h;
for i=2:2:n
sum2=sum2+f(x0);
x0=x0+2*h;
end
%para la suma 3
sum3=0;
x0=a+2*h;
for i=2:2:n-2
sum3=sum3+f(x0);
x0=x0+2*h;
end
%para la suma 4
sum4=0;
x0=a+n*h;
sum4=f(x0);
%para la suma total
sumtot=(h/3)*(sum1+4*sum2+2*sum3+sum4);
sumtot