为何基于Chebyshev-Gauss近似的erf函数正常，erfc函数始终返回0？

问题原因 #

• 你在erfc函数里将积分上限设为999999，这个值远超出实际计算需要。Chebyshev-Gauss积分的节点xi会分布在[x, 999999]区间内，当右端点极大时，绝大多数节点都会落在数值极大的区域，此时exp(-xi²)的值会趋近于机器精度下的0（比如xi>20时，exp(-400)已经远小于浮点数最小可表示非零值）。
• 所有积分项wi * f(xi)都近似为0，最终积分结果自然为0。

解决办法 #

方法一：利用erf与erfc的数学关系（推荐） #

根据定义，erfc(x) = 1 - erf(x)，直接调用已验证正确的erf函数即可，无需重新积分：

def erfc(x):
    return 1 - erf(x)

这种方式既避免了无穷区间积分的问题，还能保证计算效率与精度。

方法二：转换无穷积分到有限区间（若必须通过积分实现） #

对于∫(x到∞) exp(-t²) dt，可通过变量替换将无穷区间转为有限区间，比如令t = x + s/(1-s)（s∈[0,1)），转换后积分区间变为[0,1]，再用Chebyshev-Gauss积分计算：

def erfc(x):
    if x < 0:
        return 2 - erfc(-x)
    # 变量替换：将x到∞的积分转为0到1的积分
    def integrand(s):
        t = x + s/(1 - s)
        dt_ds = 1/(1 - s)**2
        return exp(-t**2) * dt_ds
    return 2 * integral(integrand, 0, 1) / sqrt(pi)

这种方式能避免因积分区间过大导致的数值下溢问题。

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