复化梯形公式.cpp

来自「大家请多指教」· C++ 代码 · 共 59 行

#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,n=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++)
  n=n+2*f(a+k*h);

     I=(h/2)*(m+n);

	 cout<<endl;
	 cout<<"积分结果是："<<I;


}