C++辛普森规则积分算法异常：部分函数与区间计算失效求助

1. 区间数n未保证为偶数 #

辛普森规则要求将积分区间划分为偶数个子区间（对应奇数个采样点），如果代码中n是奇数，会导致最后一段退化为梯形近似，直接破坏计算精度。比如测试sqrt(1+x)在[1,3]时，若n取奇数，结果必然偏差。

• 检查代码中n的设定：无论是固定值还是传入参数，必须确保n % 2 == 0。如果是动态传入，建议加断言或强制修正为最近的偶数。

2. 采样点权重的奇偶判断逻辑错误 #

辛普森规则的权重规则是：

• 首尾两个端点（第0个和第n个采样点）权重为1
• 下标为奇数的中间点（第1、3、...、n-1个）权重为4
• 下标为偶数的中间点（第2、4、...、n-2个）权重为2
  如果代码中把奇偶下标搞反（比如把偶数点设为权重4），会导致积分结果系统性偏差，尤其是多项式这类对权重敏感的函数，直接会算出完全错误的结果。
• 核对循环中权重赋值的逻辑：比如用i % 2判断时，确认i的起始计数是从0还是1开始，避免奇偶判断错位。

3. 自定义多项式的函数实现错误 #

你提到自定义多项式在所有区间都错误，大概率不是辛普森规则的问题，而是多项式本身的求值代码写错了。比如：

• 多项式项的系数、指数写错
• 运算优先级错误（比如漏加括号）
• 整数溢出或类型不匹配（比如用int计算高次项，导致结果截断）
• 建议单独测试多项式函数：代入几个已知x值，手动计算结果和代码输出对比，确认函数本身没问题。

4. 步长计算或最终缩放错误 #

虽然你排查过步长和h/3的乘法，但再仔细核对：

• 步长h = (b - a) / n，这里的n是子区间数，不是采样点数。如果误把采样点数当成n，会导致步长计算错误。
• 最终结果必须是(h/3) * 求和项，如果写成h * 求和项 / 2或者其他系数，结果必然错误。

5. 浮点数精度的隐性问题 #

对于sqrt(1+x)这类函数，在某些区间如果采样点过少，辛普森规则的误差会被放大，但你说仅[0,3]正确，更可能是前面的奇偶区间数问题。另外，如果代码中用float而非double，精度不足也会导致偏差，建议统一用double类型计算。

内容的提问来源于stack exchange，提问作者Filippo Pavarino