format long;
h=1/64;
f=@(x)exp(-x);
sum1=0;
sum2=0;
for(i=h:2*h:63*h)
sum1=sum1+f(i);
end
for(j=2*h:2*h:62*h)
sum2=sum2+f(j);
end
ans=h/3*(f(0)+4*sum1+2*sum2+f(1))