SymPy, 4. sıradaki polinom denklemini çözemez

Question

4. dereceden bir polinom denklemine sahibim ve tüm kökleri bulmam gerekiyor. Basit bir örnek:

from sympy import (Symbol,solve,I) 

a=4+5*I; b=3+7*I; c=12-56*I; d=33+56*I; e=345-67*I; x=Symbol('x') 
eq=a*x**4 + b*x**3 + c*x**2 + d*x +e 
solve(eq,x) 

bir olursa, b, c, d, e, sonra gayet işe, saf gerçek. Ama benim durumumda hepsi karmaşık sayılar.

PolynomialError: 'cannot return general quartic solution' 

ben tür benzer bir sorun bulmak ve düzeltmeyi uygulamak: O zaman diyoruz aldın Description of the issue. Fix of the issue

ama gerçekten yardımcı olmuyor. Garip tür sorun vardır artık çağrı (düzeltme değişti gibi) geçerli:

PolynomialError: Cannot determine if `-((12 - 56*I)/(4 + 5*I) - 3*(3 + 7*I)**2/(8*(4 + 5*I)**2))**2/12 + (3 + 7*I)*((33 + 56*I)/(4*(4 + 5*I)) + (3 + 7*I)*(3*(3 + 7*I)**2/(256*(4 + 5*I)**2) - (12 - 56*I)/(16*(4 + 5*I)))/(4 + 5*I))/(4 + 5*I) - (345 - 67*I)/(4 + 5*I)` is nonzero. 

Ama ifadesi yukarıda sıfırdan farklı en basit şey olup olmadığını belirlemek için, yani nereye sorunu bilmiyorum olabilirdi. keyfi quartic çözümleri destekler SymPy en son sürümüne

Answer 1

Yükseltme. Eğer daha esnek bir çözüm istiyorsanız, gibi bir şey ile ikili aramayı kullanarak x için

Answer 2

çözebiliriz aşağıdadır:

def maybeRightX(maybeX, polys): 
    sum = 0 
    for i in range(len(polys)): 
     sum += polys[i]*(maybeX ** i) 
    return sum 

def solve(y, polys): 
    lo = 0 
    hi = y 
    while lo <= hi: 
     mid = (lo + hi)//2 
     if (maybeRightX(mid, polys)) < y: 
      lo = mid + 1 
     else: 
      hi = mid - 1 
    return (hi + 1)