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