format long

n=64;
h=1/64;
x=0:h:1;
fx=exp(-x);
w=2*ones(1,n+1);
w(1)=1;
w(n+1)=1;
s=w*fx'*h/2;
fprintf("s=%.15f",s);