我有 sympy 的问题，我需要求解两个方程组：
B1=2.51/(Re(f^(0.5)))
f=(1/(-2*log10((epsilon/D/3.7)+B1)))^2
math.log10我尝试使用 sympy、numpy 和 fsolve，但是使用返回浮点数的 时出现问题：
B1=2.51/(Re(f^(0.5)))
f=(1/(-2*log10((epsilon/D/3.7)+B1)))^2
返回：
TypeError: can't convert expression to float