为什么这个积分很难近似?



我编写了一个用高斯正交法逼近积分的程序。我在作业中遇到的问题是为什么∫∫x,从0到1的积分很难近似?近似值在156.08左右。我认为它与范围有关,但无论我将n设置为25还是90,最后一次迭代的输出总是在156左右。如果这个问题不属于这里,请原谅。

可能是由于基本的浮点运算。在某种程度上,如果你把n设置得足够高,你最终会得到156.0796666……数值稳定性的一个很好的替代方法是在对数空间中执行计算,如果/只要可能。