a=0;
b=1;
h=1/64;
n=(b-a)/h;
x=a:h:b;
f=e.^(-x);
w=ones(1,n+1);
w(2:2:n)=4;
w(3:2:n-1)=2;
s=w*f'*h/3