我终于写出求一元四次方程解的程序了

革命年

我终于写出求一元四次方程解的程序了
这是一次到四次的代码

1. def ycfc(a,b):
   
2.     #一次方程
   
3.     return -b/a
   
4. def ecfc(a,b,c):
   
5.     #二次方程
   
6.     return (-b+((b**2)-4*a*c)**0.5)/a/2,(-b-((b**2)-4*a*c)**0.5)/a/2
   
7. def scfc(a,b,c,d):
   
8.     #三次方程
   
9.     d=b*c/(6*(a**2))-(b**3)/(27*(a**3))-d/2/a
   
10.     e=(d+(c/3/a-(b**2)/9/(a**2))**3)**0.5
    
11.     q=-b/3/a
    
12.     u=0+1j
    
13.     v=(-1+(3**0.5)*u)/2
    
14.     v1=(-1-(3**0.5)*u)/2
    
15.     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))
    
16. def hcfc(a,b,c,d,e):
    
17.     #四次方程
    
18.     d1=(c**2)-3*b*d+12*a*e
    
19.     d2=2*(c**3)-9*b*c*d+27*a*(d**2)+27*(b**2)*e-72*a*e
    
20.     d3=((2**(1/3))*d1)/(3*a*(d2+(-4*d1**3+d2*2)**0.5)*(1/3))
    
21.     d4=((d2+(-4*(d1**3)+d2**2)**0.5)*1/3)/(3*(2**1/3)*a)
    
22.     d=d3+d4
    
23.     q=-b/4/a
    
24.     w=(b**2)/4*(a**2)-2*c/3/a
    
25.     r=-b*b*b/a/a/a+4*b*c/a/a-8*d/a
    
26.     t=(w+d)**0.25
    
27.     y=r/t
    
28.     u=0.5*(w+d)**0.5
    
29.     return q-u-0.5*(w-d-y)**0.5,q-u+0.5*(w-d-y)**0.5,q+u-0.5*(w-d+y)**0.5,q+u+0.5*(w-d+y)**0.5

复制代码

PiotrPoczupow

牛批

yaoyang1001

66666666