Rm=35 t1=10 x = symbols('x') aah=400 Ki=500 r = (x / (Rm * t1) ** 0.5) / 2 Ftw = 1 + 0.07244 * r + 0.64856 * r ** 2 - 0.237 * r ** 3 + 0.038154 * r ** 4 - 0.002348 * r ** 5 # r为拉姆达,无参量因子 y=[Ftw*aah*(pi*(x/2)/1000)**0.5-Ki] print(y) v=solve(Eq(0.0223606797749979*pi**0.5*x**0.5*(-6.40336118377931e-9*x**5 + 3.89326530612245e-6*x**4 - 0.000904870204760636*x**3 + 0.0926514285714286*x**2 + 0.38720808728272*x + 200) - 500,0),x) y = [N(solution) for solution in v] print(y)
你输出你的v看看数据内容。每一步计算的结果验证过了吗?
您好,我是有问必答小助手,您的问题已经有小伙伴解答了,您看下是否解决,可以追评进行沟通哦~
如果有您比较满意的答案 / 帮您提供解决思路的答案,可以点击【采纳】按钮,给回答的小伙伴一些鼓励哦~~
ps:问答VIP仅需29元,即可享受5次/月 有问必答服务,了解详情>>>https://vip.csdn.net/askvip?utm_source=1146287632