f=@(x)1/5*x^2;
a=2;
b=7;
n=32;
dx=(b-a)/n;
equ=0;
sum = 0;
for i=0:n
equ = a + i*dx;
if i==0 || i==n 
  sum=sum+f(equ);
elseif mod(i,2)==0 
  sum=sum+4*f(equ);
else
  sum=sum+2*f(equ);
endif
end
Res=sum*(dx/3);
Res