怎么用C语言表示积分

怎么用C语言表示积分,第1张

源代码如下:

#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

    由上面的比较可知,梯形法的精度要高于矩形法。


欢迎分享,转载请注明来源:内存溢出

原文地址: http://outofmemory.cn/yw/11656438.html

(0)
打赏 微信扫一扫 微信扫一扫 支付宝扫一扫 支付宝扫一扫
上一篇 2023-05-17
下一篇 2023-05-17

发表评论

登录后才能评论

评论列表(0条)

保存