1

我如何评估 xe^x/(e^x-1) 与零附近的数值稳定性以及当 x 非常积极或消极？我可以访问 和 中的所有常用数学numpy函数scipy。

4

1 回答 1

2 
def f(x):
  if abs(x) > 0.1: return x*exp(x)/(exp(x)-1)
  else: return 1/(1.-x/2.+x**2/6.-x**3/24.)

如果需要更高的精度，可以以明显的方式扩展最后一行的扩展，并且可以通过重新措辞来加快速度。就目前而言，它的误差高达 1e-6。

于 2016-09-03T05:03:40.340 回答