我终于写求一元三次方程的程序了
我终于写求一元三次方程的程序了。这是源代码
def scfc(a,b,c,d):
#三次方程
d=b*c/(6*(a**2))-(b**3)/(27*(a**3))-d/2/a
e=(d+(c/3/a-(b**2)/9/(a**2))**3)**0.5
q=-b/3/a
u=0+1j
v=(-1+(3**0.5)*u)/2
v1=(-1-(3**0.5)*u)/2
return q+(d+e)**(1/3)+(d-1)**(1/3),q+v*((d+e)**(1/3))+v1*((d-e)**(1/3)),q+v1*((d+e)**(1/3))+v*((d-e)**(1/3))
页:
[1]