如何通过Symphy计算python中对数方程的解

Question

我试着用这个公式，用sympy求Moody图（f）的摩擦系数

from sympy import *
f = Symbol ('f')
Re = 3000
relative_roughness=0.1
solve(f**(-0.5)+2*log(relative_roughness/3.7+2.51/(Re*(f)**(0.5)), 10), domain=S.Reals)

但程序结束时出现错误NotImplementedError:multiple generators[log（1/37+83667/（10000000000000000sqrt（f）），sqrt（f）] 方程2log（1/37+83667/（10000000000000000*sqrt（f）））/log（10）+1/sqrt（f）未实现任何算法

当相对粗糙度=0时，它工作，但我需要大于0。请帮忙。如何修复它

Answer 1

使用sympy的nsolve函数。这里有一个0.1的初始猜测：

In [29]: eq = f**(-0.5)+2*log(relative_roughness/3.7+2.51/(Re*(f)**(0.5)), 10)                                                    

In [30]: nsolve(eq, f, 0.1)                                                                                                       
Out[30]: 0.106947153535321