#include<stdio.h>

#include<math.h>

#define F(x) exp(x)

 

main(void)

{

 float a,b,h,x,s;

 int n,i;

 

 a=0.;

 b=1.;

 n=1;

 

 do{

    s=F(a)+F(b);

    h=(b-a)/n;

    for(i=1;i<=(n-1);++i){

       x=a+i*h;

       s+=2.*F(x);}

    s*=h/2.;

    printf("n = %4d:  S = %7.5e\n",n,s);

    n*=2;

   }while(n<=pow(2,20));

}