#include
#include
using namespace std;
//梯形法求函数f(x)的定积分
double integral(double a,double b, double (*f)(double),int n){
double h,x,sum;
int i;
h=(b-a)/n;
x=a;
sum=(f(a)+f(b))/2;
for(i=1; i
sum += f(x);
}
sum *= h;
return sum;
}
void main(){
double a,b,sum;
int i,n;
cout <<"Input a b and n: ";
cin >>a >>b >>n;
sum=integral(a,b,sin,n);
cout <