Я пытаюсь решить два длинных нелинейных уравнения одновременно, используя 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?