在Python中如何计算数值稳定性？

Question

如何在数值稳定性为零的情况下，以及当x为正或负时，如何计算xe^x/（e^x-1）？我可以使用numpy和scipy中所有常用的数学函数。在

Answer 1

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。在