用nsolve解两个非线性方程

Question

我正在尝试使用nsolve同时求解两个长非线性方程。我得到一个ZeroDivisionError，我认为这是因为我没有为任何一个变量设置初始猜测。在方程式中，Rci，Rfo，T_ci和qp_in都是常数。

from sympy import Symbol, nsolve
x=Symbol('x')
y=Symbol('y')
nsolve([x-T_ci-qp_in/(np.pi*(Rci+Rfo)*y), y-(1.58e-3* \\
((x+T_ci)/2.0)**.79)/delta+(sigma/(1/2.0))*(x**4-T_ci**4)/(x-T_ci)], \\ 
[x, y], [0, 0])

如何使用nsolve方法对x和y进行初始猜测，这样我就不会得到ZeroDivisionError？除了nsolve之外，还有解决这些方程的更好方法吗？