数值积分:C中的梯形法则

时间:2015-03-12 14:46:17

标签: c integration numerical

我试图使用梯形法则将函数1 /((1 + x)(x ^ 0.5))在0和无穷大之间进行积分。我知道这是一个不正确的整体,梯形方法不是理想的,但它是我被告知使用的。我没有尝试直接在0和无穷大之间进行计算,而是试图将积分分成不同的极限并总结它们。由于某种原因,while循环的第一次迭代(限制,0和1)起作用并产生正确的答案(pi / 2),但后来的答案却没有。

有人可以告诉我我的代码在哪里

1 个答案:

答案 0 :(得分:2)

sum未初始化。你在第一次传球时很幸运,它是0.但是下一次通过外环时,它仍然具有第一次的值。在进入内循环并开始重新添加之前,需要将其设置为0.