用c语言求定积分

#include
#include
double integral(double(*fun)(double x),double a,double b,int,n){
double s,h,y;
int i;
s=(fun(a)+fun(b))/2;
h=(b-a)/n; /*积分步长*/
for(i=1;i s=s+fun(a+i*h);
y=s*h;
return y;/*返回积分值*/
}
double f(double x){
return(x*sinx) /*修改此处可以改变被积函数*/
}
int main(){
double y;
y=integral(f,1.0,2.0,150);/*修改此处可以改变积分上下限和步数,步长=(上限-下限)/步数*/
printf("y=%f\n",y);
return 0;
}
int main()

#include
#include
int main()
{
double a,b,n,k,i,x,y;
double EPS=1e-7;
a=0;//a为积分下限
b=1.57079632679459661;//b为积分上限
n=1e-7;//n控制精度
k=x=y=0;//k为辅助变量
i=0;//i为积分的值

while(k<=1/n)
{
k++;
x=k*n*b-k*n*a;
y=sin(x);//定义被积函数
i=i+b*y*n-a*y*n;
}
printf("%lf\n",i);
getchar();
}