用渐近法计算python对数方程的解
用渐近法计算python对数方程的解
提问于 2020-12-04 11:19:10
我试图用这个公式求出穆迪图(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结束:多个生成器日志(1/37+ 836666666666667/(1000000000000000000*sqrt(f)))/log(10) sqrt(f) )没有实现求解方程2_log(1/37 +836666666666667/(1000000000000000000*sqrt(f)))/log(10)+ 1/sqrt(f) )的算法。
它适用于relative_roughness=0，但我需要超过0。拜托救救我。怎么修呢？
回答 1
Stack Overflow用户
发布于 2020-12-04 11:47:46
使用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页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/65142641
复制相关文章
相似问题
腾讯云开发者
