源代码如下:
#include#includefloat f1(float x)
{
return(1.0+x)
}
float f2(float x)
{
return(2.0*x+3.0)
}
float f3(float x)
{
return(exp(x)+1)
}
float f4(float x)
{
return(pow(1+x,2))
}
float f5(float x)
{
return(pow(x,3))
}
float fsimp(float a,float b,float (*p)(float))
{
float c,s
c=(a+b)/2
s=(b-a)/6*(p(a)+4*p(c)+p(b))
return s
}
int main()
{
float a,b
scanf("%f",&a)
printf("请输入积分上限b的值:")
scanf("%f",&b)
printf("%f\n",fsimp(a,b,f1))
}
扩展资料
1、对应于一个积分式要有一段程序,可以改变程序的一小部分来改变所要求的积分式。
2、除数不能位0。
3、两个整数相除,结果仍是整数。
4、若被除数其中有一个为浮点数或者两个都为浮点数,则结果为浮点类型。操作数必须为整数,不能是浮点数。
#include #include float f1(float x) { return(1.0+x)} float f2(float x) { return(2.0*x+3.0)} float f3(float x) { return(exp(x)+1)} float f4(float x) { return(pow(1+x,2))} float f5(float x) { return(pow(x,3))} float fsimp(float a,float b,float (*p)(float)) { float c,sc=(a+b)/2s=(b-a)/6*(p(a)+4*p(c)+p(b))return s} int main() { float a,bprintf("请输入积分下限a的值:")scanf("%f",&a)printf("请输入积分上限b的值:")scanf("%f",&b)printf("%f\n",fsimp(a,b,f1))printf("%f\n",fsimp(a,b,f2))printf("%f\n",fsimp(a,b,f3))printf("%f\n",fsimp(a,b,f4))printf("%f\n",fsimp(a,b,f5))}基本是这样的,用梯形发求定积分,对应于一个积分式就要有一段程序,不过你可以改变程序的一小部分来改变你所要求的积分式。以c为例:求f(x)=xsinx从1到2的积分
#include <math.h>
float integral(float(*fun)(float x),float a,float b,int,n)
{float s,h,y
int i
s=(fun(a)+fun(b))/2
h=(b-a)/n/*积分步长*/
for(i=1i<ni++)
s=s+fun(a+i*h)
y=s*h
return y/*返回积分值*/
}
float f(float x)
{return(x*sinx) /*修改此处可以改变被积函数*/
}
main()
{float y
y=integral(f,1.0,2.0,150)/*修改此处可以改变积分上下限和步长*/
printf("y=%f\n",y)
}