#include <stdio.h>
#include <math.h>
#define f(x) (1/x)
//float f(float a)
//{return 1/a;}
main()
{float r1,r2,s1,s2,t1,t2,c1,c2,h,s,a,b,ep;
 int i,n=1,k=0;
 printf("請輸入積分區間：");
 scanf("%f,%f",&a,&b);
 printf("\n");
 printf("請輸入精度：ep=");
 scanf("%f",&ep);
 printf("\n");
 h=(b-a)/2;
 t2=(f(a)+f(b))*h;
 s2=0;c2=0;r2=0;
 do
 {
	 r1=r2;t1=t2;s1=s2;c1=c2;
	 h=h/2.0;
	 s=0;
	 for(i=1;i<=n;i++)
		 s=s+f((a+(2*i-1)*2*h));
	 t2=t1/2+s*2*h;
	 s2=(4*t2-t1)/3;
	 c2=(16*s2-s1)/15;
	 r2=(64*c2-c1)/63;
	 n=n*2;k++;
 }while(fabs(r2-r1)>ep);
 printf("積分結果為：%f",r2);
 printf("\n");
}