在Python Fsolve中求解具有动态未知数的非线性方程系统方程

Question

我想解决一个非线性方程组，它由数量众多的未知数组成。我已经开始使用“静态”方法，并且效果很好。基本上，代码大致如下：

from scipy.optimize import fsolve

def fixEqSyst(z):
v1 = z[0]
v2 = z[1]
v3 = z[2]
v4 = z[3]
f=np.zeros(4)
f[0] = 2*v1-3*v2+7*v3**2
f[1] = v1+3*v2**2-9*v3
f[2] = -3v1**2+12*v2+7*v3
f[3] = 4*v1+5*V2*v3
return f

z = fsolve(fixEqSyst, [0, 0, 0, 0])

如何转换函数以同时管理动态数量的变量和函数？我曾尝试应用for-loop策略，但没有成功，您还有其他解决方案吗？

此外，让我们考虑一下，作为输入，我将永远知道#of变量和#of函数。