#include<iostream.h>
#include<math.h>
#include<stdio.h>
 


double f(double s)
 {      double t;
        t=1/(1+s*s);
        return t;
 }


void main()
{   int n1,k;
    double h,a,b,I,m,n2=0,n3=0;

	cout<<"劃分的區間數：";
	cin>>n1;

	cout<<"請輸入上限值：";
	cin>>b;

	cout<<"請輸入下限值：";
	cin>>a;

    h=(b-a)/n1;
    

  m=f(a)+f(b);
	
 for( k=1;k<n1;k++)
 {
	 n2=n2+4*f(a+k*h);
	 k++;
 }
 for( k=2;k<n1;k++)
 {
     n3=n3+2*f(a+k*h);
	 k++;
 }

     I=(h/3)*(m+n2+n3);

	 cout<<endl;
	 cout<<"積分結果是："<<I;


}