第八章习题13-写一个用矩形法求定积分的通用函数,分别求积分区间为[0,1]sinx,cosx,e的x方的定积分

在这里插入图片描述

 🌏个人博客:尹蓝锐的博客

在这里插入图片描述

希望文章能够给到初学的你一些启发~ 如果觉得文章对你有帮助的话,点赞 + 关注+ 收藏支持一下笔者吧~

1、题目要求:

这个积分是一个比较特殊的积分,它没有原函数,需要通过数值积分法来求解。一种比较常用的数值积分法是梯形法则,其基本思路是将积分区间划分成若干个小区间,然后在每个小区间上使用梯形面积来近似计算积分值。 具体来说,在区间[0,1]上均匀取n个点,将区间划分成n-1个小区间,每个小区间的宽度为h=1/n。然后,我们可以将这些小区间上的积分近似地表示为梯形面积之和,即: ∫[0,1] sinx/x dx ≈ h/2 * (y0 + yn + 2(y1+y2+...+yn-1)) 其中,y0 = sin(0)/0 = 1,yn = sin(1)/1 ≈ 0.8415,yi = sin(xi)/xi,xi = i*h。 将上式代入,得到: ∫[0,1] sinx/x dx ≈ h/2 * (1 + 0.8415 + 2(y1+y2+...+yn-1)) 因此,只需要计算出y1+y2+...+yn-1的值即可求得积分的近似值。这可以通过数值计算来实现,例如使用循环来累加每个小区间上的积分值即可。 代码实现如下: ``` #include <stdio.h> #include <math.h> int main() { int n = 100; // 区间划分数 double h = 1.0/n; // 小区间宽度 double sum = 0.0; // 积分值的累加器 // 计算每个小区间上的积分值并累加 for (int i = 1; i < n; i++) { double x = i*h; double y = sin(x)/x; sum += y; } // 计算积分的近似值 double integral = h/2 * (1 + sin(1)/1 + 2*sum); printf("∫[0,1] sinx/x dx ≈ %.6f\n", integral); return 0; } ``` 运行结果为:∫[0,1] sinx/x dx ≈ 0.946083。 注意:由于这是一个数值积分法,所以计算结果只是一个近似值,并不是精确值。积分的精确值可能需要使用其他更加高级的数值积分法来求解。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

尹蓝锐

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值