如何为我的Fourier系列避免RuntimeWarning

Question

对于我的偏微分方程的特定解析解，它可以表示为傅里叶级数，我将其编码为：

def wire_anal(x, y, L, N):
series = 0
for n in range(1, N):
    series += (4000/((2*n)-1))*np.sin(((2*n)-1)*np.pi*x/L)*((np.sinh(((2*n)-1)*np.pi*y/L))/(np.sinh(np.pi*((2*n)-1))))
return series

当我将其运行为特定值时，我想说：

print(wire_anal(0.25, 0, 0.5, 150))

我得到

的错误

RuntimeWarning: overflow encountered in sinh

更多术语，请说：

print(wire_anal(0.25, 0.5, 0.5, 1500))

我得到一个NaN值。有没有办法解决这个问题？