源代码如下:
#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
printf("请输入积分下限a的值:")
scanf("%f",&a)
printf("请输入积分上限b的值:")
scanf("%f",&b)
printf("%f\n",fsimp(a,b,f1))
}
扩展资料
1、对应于一个积分式要有一段程序,可以改变程序的一小部分来改变所要求的积分式。
2、除数不能位0。
3、两个整数相除,结果仍是整数。
4、若被除数其中有一个为浮点数或者两个都为浮点数,则结果为浮点类型。 *** 作数必须为整数,不能是浮点数。
Matlab计算定积分的具体步骤如下:1、以f(x)=e^2x+sin(x+π/3) ,积分下限:a=0,积分上限:b=π/4 为例。首先,建立被积函数M文件;点击New,选择Function;
2、在编辑器中输入指令:function f=f(x);f=exp(2*x).*sin(x+pi/3)(注意,用 .* )
3、保存函数文件到自己的工作路径,比如G:\matlab\work。
输入指令:cd G:\matlab\work,是该路径成为当前路径。
4、matlab中计算积分的两种指令:
1.F=quad('fname',a,b,tol,trace) Simpson数值积分法
2.F=quad8('fname',a,b,tol,trace) Newton-Cotes数值积分法
其中: fname是被积函数表达式或函数名,a,b分别是上下限,tol可以控制积分精度,省略则取0.001;trace=1则用图形表示积分过程,trace=0,没有图形。
两者 比较,quad8精度更高。
5、调用积分函数squad进行计算。输入指令: F=quad('f(x)',0,pi/4),如图,回车后可得到计算结果1.8612。其他函数,只要修改函数文件中的表达式即可。
实际问题描述:
求定积分近似值
程序代码如下:
#include
#include
void main()
{
int i,n=1000
float a,b,h,t1,t2,s1,s2,x
printf("请输入积分限a,b:")
scanf("%f,%f",&a,&b)
h=(b-a)/n
for(s1=0,s2=0,i=1i<=ni++)
{
x=a+(i-1)*h
t1=(float)exp(-x*x/2)t2(float)=exp(-(x+h)*(x+h)/2)
s1=s1+t1*h /*矩形面积累加*/
s2=s2+(t1+t2)*h/2 /*梯形面积累加*/
}
printf("矩形法算得积分值:%f.\n",s1)
printf("梯形法算得积分值:%f.\n",s2)
}
程序运行结果如下:
矩形法算得积分值:0.855821
梯形法算得积分值:0.855624
由上面的比较可知,梯形法的精度要高于矩形法。
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)