例题:编写一个用矩形求定积分的通用库函数。
p.s.这个只是近似值。。。
代码如下
#include<math.h>
#include<stdio.h>
double integral(double(*fun)(double),double a,double b)
{double h,y=0;
int n,i;
n=100;
h=(b-a)/n;
for(i=0;i<n;i++)
{y+=(*fun)(a+i*h)*h;
}
return (y);
}
int main()
{double y1,y2,y3;
y1=integral(sin,0.0,1.0);
y2=integral(cos,-1.0,1.0);
y3=integral(exp,0.0,2.0);
printf("y1=%lf\n",y1);
printf("y2=%lf\n",y2);
printf("y3=%lf\n",y3);
return 0;
}
运行结果
下一篇 指针15(指针数组)