实际问题描述:
求定积分近似值
程序代码如下:
#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
由上面的比较可知,梯形法的精度要高于矩形法。
#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=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
}
int main()