Sympy提取幂级数第k项系数：复杂表达式结果异常问题

测试场景1（正常工作） #

当目标表达式为exp(y)，提取y的第2项系数时，以下代码可正确返回结果1/2：

from sympy import symbols, series, exp
x, y, z = symbols('x y z')
expr = exp(y)
k=2
series_exp = series(expr, y, n=10)  # 避免与sympy.series函数重名
kth_coeff = series_exp.coeff(y, k) 
print(kth_coeff)

测试场景2（存在异常） #

当目标表达式为-exp(3*x/(z*y))*exp(z*y/x + 1/x)，提取y的第2项系数时，使用类似逻辑得到的结果为-z**2*exp(3*x/(z*y) + 1/x)/(2*x**2)，该结果明显不符合预期的第2项系数。

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