我终于写求一元三次方程的程序了

革命年

我终于写求一元三次方程的程序了。
这是源代码

1. def scfc(a,b,c,d):
   
2.     #三次方程
   
3.     d=b*c/(6*(a**2))-(b**3)/(27*(a**3))-d/2/a
   
4.     e=(d+(c/3/a-(b**2)/9/(a**2))**3)**0.5
   
5.     q=-b/3/a
   
6.     u=0+1j
   
7.     v=(-1+(3**0.5)*u)/2
   
8.     v1=(-1-(3**0.5)*u)/2
   
9.     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))

复制代码