【python求一元三次方程】

python求一元三次方程

# 方程求解
import sympy as sy
import sympy
import sympy as sp # 导入sympy包
import math

a=math.pi
#-----------------方法一----------------------
x = sp.Symbol('x') # 定义符号变量
f = 0.8*x**3+47.6*x**2 + 303.8*x + 514.5-1286.25*a # 定义要求解的一元三次方程
x = sp.solve(f) # 调用solve函数求解方程
print("x=",x)
V=0.1
r=3.5
l=105
L=4*l
N=(V*L**3)/(l*a*r**2)
print("N=",N)
#----------方法二-----------
dd=sy.symbols('dd')
fx=V*4*(l+2*dd)*((r+dd)**2)-l*a*r**2 
answer=sy.solve([fx],[dd]) 
print('dd=',answer) #将answer打印出来即可。
#---------方法三---------------------
d = sympy.Symbol('d')
result = sympy.solve(V*4*(l+2*d)*((r+d)**2)-l*a*r**2,d)
print("result=",result)