我有一个四阶多项式方程,我需要找到所有的根.

简单的例子:

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)

如果a,b,c,d,e是纯实的,那么

它工作得很好.

但就我而言,所有这些都是复杂的数字.然后我接到了电话:

PolynomialError: 'cannot return general quartic solution'

但它并没有真正帮助.存在某种奇怪的问题,因为现在调用(在修复中已更改):

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.

但是要确定上面的表达式是否非零是最简单的事情,所以不知道问题可能在哪里.

python解多项式方程_python – SymPy不能求解四阶多项式方程相关推荐

  1. 如何利用python解方程_Python 解方程的三种方法

    # 首发于我的博客 The North. 新年第一篇,搞起. 这回写一个好久之前想做,一直搁着没做的东西-- Python 解方程(其实是放假回家,趁着家里电脑重装 LOL 的时间过来写一篇). 咱这 ...

  2. python计算矩阵方程_python/sympy求解矩阵方程的方法

    sympy版本:1.2 假设求解矩阵方程 AX=A+2X 其中 求解之前对矩阵方程化简为 (A−2E)X=A 令 B=(A−2E) 使用qtconsole输入下面程序进行求解 In [26]: fro ...

  3. python解不定积分_python快速求解不定积分和定积分

    欢迎点击「算法与编程之美」↑关注我们! 本文首发于微信公众号:"算法与编程之美",欢迎关注,及时了解更多此系列博客. 基本概念 定积分的定义如下: 不定积分定义如下: 如果想了解更 ...

  4. python解常微分方程龙格库_求解二阶常微分方程的RungeKutta四阶方法

    我试着做一个简谐振子的例子,它将用龙格-库塔四阶法求解.要求解的二阶常微分方程(ODE)和初始条件为: y''+y=0 y(0)=0和y'(0)=1/pi 范围在0到1之间,共有100步.我用u作为辅 ...

  5. python 解方程_python如何解方程的三种方法

    python求解方程组的三种方法: Numpy求解方程组x + 2y = 3 4x + 5y = 6 当然我们可以手动写出解析解,然后写一个函数来求解,这实际上只是用 Python 来单纯做" ...

  6. python解非线性方程_python解非线性方程的问题

    废话不多说,上代码 方法1: from sympy import * t=symbols('t') eq=eq(10**20*exp(-12000/t)/(0.01+10**14*exp(-12000 ...

  7. 如何用python解方程组_python如何解方程组

    解方程的基本思想是: aa = solve([f],[x]) f是方程,x是变量,这个代码解的是关于x的方程f=0 既然可以解方程(组),就一定可以解微分方程(组). from sympy impor ...

  8. python解压_python解压缩

    解压缩 如果我们给出一个列表,我们需要一次性取出多个值,我们是不是可以用下面的方式实现呢? name_list = ['chen', 'python', 'jason'] x = name_list[ ...

  9. python图像数独_Python图像识别+KNN求解数独的实现

    Python-opencv+KNN求解数独 最近一直在玩数独,突发奇想实现图像识别求解数独,输入到输出平均需要0.5s. 整体思路大概就是识别出图中数字生成list,然后求解. 输入输出demo 数独 ...

最新文章

  1. maven 插件之 AutoConfig 工具使用笔记
  2. iptables中文man文档
  3. hihoCoder 1227 2015 北京网络赛 A题
  4. Tomcat工作原理及简单模拟实现
  5. C#LeetCode刷题之#342-4的幂(Power of Four)
  6. MFC开发IM--MFC非模态对话框的创建及显示
  7. 【vue开发问题-解决方法】(九)使用element upload自定义接口上传文件,input多文件上传
  8. [二维压缩] | 条件图像游程编码
  9. 软件项目管理实验一补充
  10. 您的计算机无法加入域,Win7系统电脑无法加入域提示找不到网络路径的应对方案...
  11. rac多scan-ip配置
  12. python基础(八):封装、继承、多态
  13. 计算机网络应用层题库
  14. Android app 后台被杀恢复
  15. 【工大SCIR】对话中的情感分析与生成简述
  16. 相似度计算之Jaccard系数
  17. vue.js 回显textarea入的空格和换行
  18. Android错误:unexpected text found in layout file
  19. 南都周刊 别了 老兵乔布斯
  20. Android studio 改app图标,名字及一键拨号

热门文章

  1. 组合数函数-快速提取所有可能的组合数
  2. ML基石_9_LinearRegression
  3. 便携式不锈钢管道焊接机器人_304不锈钢管居然可以发黑!?
  4. python 函数 默认参数
  5. 回归模型-线性回归算法
  6. Nvidia TX2 安装中文输入法
  7. 语义分割--End-to-End Instance Segmentation with Recurrent Attention
  8. 目标定位--Deep Self-Taught Learning for Weakly Supervised Object Localization
  9. [librosa]0.8版本后不再支持write_wav
  10. Linux: what's the difference between /proc/self and /proc/$$?