format long
X=0:1/64:1;
fX=0:1/64:1;
even=0;
odd=0;
ans=0;
h=1/64;


for i=0:64
  
  fX(i+1)=exp(-X(i+1));
  
  end
 for j=1:31
    
    even=even+fX(2*j+1);
    
end
 for j=1:32
   
    odd=odd+fX(2*j);
   
end
  ans=(h/3)*(fX(1)+fX(65)+2*even+4*odd);
  disp(ans);