这不是多项式,因为d升至分数级。充满浮点常数的方程式不太可能受益于像SymPy这样的符号包。使用SciPy可以很容易地获得数字实值解决方案:
from scipy.optimize import fsolve
from math import pi
fsolve(lambda d: 0.125567481 * (d+0.04) - d**2.25 * (pi*2.11+0.10053 + 2.513274*d), 1)
退货 array([ 0.06086999])
请注意,我清除了分母(将两边都乘以),因为除法使数值解更难。大致而言:表达式看起来像多项式越多越好。然后我按照fsolve语法的要求重新排列了表达式,使其等于0 。
数字1是解决方案的猜测;它不一定是好人;这只需要花点时间:我们需要一个正数根,所以让我们从一个正数开始。
本文收集自互联网,转载请注明来源。
如有侵权,请联系 [email protected] 删除。
我来说两句