源代码如下:
#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、若被除数其中有一个为浮点数或者两个都为浮点数,则结果为浮点类型。操作数必须为整数,不能是浮点数。
积分分为两种,数值积分,公式积分。
公式积分:部分函数可以直接用公式求得其不定积分函数。C语言中可以直接用积分公式写出其积分函数。
数值积分:按照积分的定义,设置积分范围的步长,用梯形面积累加求得其积分。
以【f(x)=x*sin(x) 从1到2的积分】为例:
#include <math.h>#include <stdio.h>
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=1i<ni++)
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
}