R语言integrate无穷区间积分报non-finite function value错误

integrate()在无穷区间做数值积分时，会自动在全区间采样计算被积函数值，当采样到足够大的正x值时，被积函数中的指数项会发生数值溢出，触发非有限值错误：

• R使用双精度浮点数存储数值，能表示的最大正数值约为1.8e+308，对应exp(709.78)，一旦指数运算的输入超过709.78，结果就会返回Inf
• 你的被积函数中存在exp(0.5*x)项，当x>1420时，0.5*x>710，exp(0.5*x)直接计算为Inf，后续运算会产生NaN或其他非有限值，最终报错
• x取绝对值很大的负值时，exp(0.5*x)会快速趋近于0，不会触发溢出问题，不需要额外处理
• 另外你写的被积函数存在冗余计算：0.5*exp(0.5*x)*2可以直接化简为exp(0.5*x)，多余的乘法不影响结果但属于无效计算。

通过分段计算被积函数规避数值溢出即可，逻辑非常直接：

• 当x ≤ 1400时（留足够安全余量避免溢出），按原公式正常计算
• 当x > 1400时，exp(0.5*x)已经是极大值，exp(-exp(0.5*x))无限趋近于0，因此1-exp(-exp(0.5*x))项可以直接取值为1，不需要再计算指数项，从根源避免溢出

修正后的可运行代码如下：

# 修正后带溢出保护的被积函数
int_fixed <- function(x) {
  base_term <- x * exp(0.5*x) * dnorm(x)
  adjust_term <- ifelse(
    x > 1400,
    1,
    1 - exp(-exp(0.5*x))
  )
  base_term * adjust_term
}

# 执行无穷区间积分
integrate(int_fixed, lower = -Inf, upper = Inf)

运行上述代码即可正常得到积分结果，不会触发非有限值报错。如果需要更高的计算精度，也可以根据需求调整分段阈值，只要保证阈值小于1419（即计算exp(0.5*x)时不会溢出）即可。

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