Matlab - Urgente problemas con integrales

clear

syms x

a=0; b=1; h=1/16;

n=(b-a)/h;

sp=0;  si=0;

f(x)=exp(x);

for i=1:n-1

    if rem(i,2)==0

      sp=sp+2*f(a+i*h);

    else

       si=si+4*f(a+i*h);

    end

end

s=f(a)+f(b)+sp+si;

Qh=double(h*s/3)

VE=double(f(b)-f(a));

Er=abs((VE-Qh)/VE)

>> simpson_compuesto

Qh =

    1.7183

Er =

   8.4732e-08