[英]Cardano's Formula in Python
我需要有人檢查我在 Python 中使用卡爾達諾公式的嘗試。 我正在嘗試求解三次方程的根,我想知道到目前為止我所做的是否正確(或錯誤)。 TYIA
def solve(a,b,c,d):
Q = (3*a*c - (b**2)) / (9*(a**2))
R = (9*a*b*c - 27*(a**2)*d - 2*(b**3)) / (54*(a**3))
D = (Q**3) + (R**2)
S = (R + (D**(1/2)))**(1/3)
T = (R - (D**(1/2)))**(1/3)
x1 = S + T - (b/(3*a))
x2 = -((S + T)/2) - (b/(3*a)) + 0.5j * (3**(1/2)) * (S - T)
x3 = -((S + T)/2) - (b/(3*a)) - 0.5j * (3**(1/2)) * (S - T)
return (x1,x2,x3)
您可以將根代回方程式中,看看是否為零。
def verify(a,b,c,d,x):
return (a*(x**3) + b*(x**2) + c*x + d)
roots = solve(a,b,c,d)
for x in roots:
print(verify(a,b,c,d,x))
你寫錯了三次方程的解。 一般來說,沒有方程組是不可能的,所以我可以提供以下解決方案
from cmath import *
solutions = set()
def cbrt(polynomial):
solution = set()
root1 = polynomial ** (1 / 3)
root2 = (polynomial ** (1 / 3)) * (-1 / 2 + (sqrt(3) * 1j) / 2)
root3 = (polynomial ** (1 / 3)) * (-1 / 2 - (sqrt(3) * 1j) / 2)
solution.update({root1, root2, root3})
return solution
print('ax^3+bx^2+cx+d=0')
a, b, c, d = map(complex, input().split())
if a != 0:
p = (3 * a * c - b ** 2) / (3 * a ** 2)
q = (2 * b ** 3 - 9 * a * b * c + 27 * a ** 2 * d) / (27 * a ** 3)
alpha = cbrt(-q / 2 + sqrt((q / 2) ** 2 + (p / 3) ** 3))
beta = cbrt(-q / 2 - sqrt((q / 2) ** 2 + (p / 3) ** 3))
for i in alpha:
for j in beta:
if abs((i * j) + p / 3) <= 0.00001:
x = i + j - b / (3 * a)
solutions.add(x)
elif b != 0:
D = c ** 2 - 4 * b * d
solutions.add((-c - sqrt(D)) / 2 * b)
solutions.add((-c + sqrt(D)) / 2 * b)
elif c != 0:
solutions.add(-d / c)
elif d == 0:
solutions.add("True")
else:
solutions.add("False")
print(solutions)
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.