如果v0和[s]是标量,上面的公式可写为f(Km, Vmax)=0的形式。假定km在区间[a,b]内方程有解,可在该区间内均匀取n点,计算一元方程得到n个Vmax的解。下面的代码,演示了非线性方程的解法。
from scipy import optimize
def f(x):
return 3/(x*5) + 1/x - 1/7
result = optimize.root_scalar(f, bracket=[-2,2])
result.root
-1.2566741077549764e-12
不知道你这个问题是否已经解决, 如果还没有解决的话: