clc
clear all

f=inline('exp(-x)*cos(4*x)')
a=1;
b=3;
n=12;
h=(b-a)/n;
x0=a;
%para la suma 1
sum1=f(x0);
%para la suma 2
sum2=0;
x0=a+h;
for i=1:3:n-2
sum2=sum2+f(x0);
x0=x0+3*h;
end
%para la suma 3
sum3=0;
x0=a+2*h;
for i=2:3:n-1
sum3=sum3+f(x0);
x0=x0+3*h;
end
%para la suma 4
sum4=0;
x0=a+3*h;
for i=3:3:n-3
sum4=sum4+f(x0);
x0=x0+3*h;
end
%para la suma 5
sum5=0;
x0=a+n*h;
sum5=f(x0);
%para la suma total
sumtot=(3*h/8)*(sum1+3*sum2+3*sum3+2*sum4+sum5);
sumtot