python解多项式方程_python – SymPy不能求解四阶多项式方程
我有一个四阶多项式方程,我需要找到所有的根.
简单的例子:
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不能求解四阶多项式方程相关推荐
- 如何利用python解方程_Python 解方程的三种方法
# 首发于我的博客 The North. 新年第一篇,搞起. 这回写一个好久之前想做,一直搁着没做的东西-- Python 解方程(其实是放假回家,趁着家里电脑重装 LOL 的时间过来写一篇). 咱这 ...
- python计算矩阵方程_python/sympy求解矩阵方程的方法
sympy版本:1.2 假设求解矩阵方程 AX=A+2X 其中 求解之前对矩阵方程化简为 (A−2E)X=A 令 B=(A−2E) 使用qtconsole输入下面程序进行求解 In [26]: fro ...
- python解不定积分_python快速求解不定积分和定积分
欢迎点击「算法与编程之美」↑关注我们! 本文首发于微信公众号:"算法与编程之美",欢迎关注,及时了解更多此系列博客. 基本概念 定积分的定义如下: 不定积分定义如下: 如果想了解更 ...
- python解常微分方程龙格库_求解二阶常微分方程的RungeKutta四阶方法
我试着做一个简谐振子的例子,它将用龙格-库塔四阶法求解.要求解的二阶常微分方程(ODE)和初始条件为: y''+y=0 y(0)=0和y'(0)=1/pi 范围在0到1之间,共有100步.我用u作为辅 ...
- python 解方程_python如何解方程的三种方法
python求解方程组的三种方法: Numpy求解方程组x + 2y = 3 4x + 5y = 6 当然我们可以手动写出解析解,然后写一个函数来求解,这实际上只是用 Python 来单纯做" ...
- python解非线性方程_python解非线性方程的问题
废话不多说,上代码 方法1: from sympy import * t=symbols('t') eq=eq(10**20*exp(-12000/t)/(0.01+10**14*exp(-12000 ...
- 如何用python解方程组_python如何解方程组
解方程的基本思想是: aa = solve([f],[x]) f是方程,x是变量,这个代码解的是关于x的方程f=0 既然可以解方程(组),就一定可以解微分方程(组). from sympy impor ...
- python解压_python解压缩
解压缩 如果我们给出一个列表,我们需要一次性取出多个值,我们是不是可以用下面的方式实现呢? name_list = ['chen', 'python', 'jason'] x = name_list[ ...
- python图像数独_Python图像识别+KNN求解数独的实现
Python-opencv+KNN求解数独 最近一直在玩数独,突发奇想实现图像识别求解数独,输入到输出平均需要0.5s. 整体思路大概就是识别出图中数字生成list,然后求解. 输入输出demo 数独 ...
最新文章
- maven 插件之 AutoConfig 工具使用笔记
- iptables中文man文档
- hihoCoder 1227 2015 北京网络赛 A题
- Tomcat工作原理及简单模拟实现
- C#LeetCode刷题之#342-4的幂(Power of Four)
- MFC开发IM--MFC非模态对话框的创建及显示
- 【vue开发问题-解决方法】(九)使用element upload自定义接口上传文件,input多文件上传
- [二维压缩] | 条件图像游程编码
- 软件项目管理实验一补充
- 您的计算机无法加入域,Win7系统电脑无法加入域提示找不到网络路径的应对方案...
- rac多scan-ip配置
- python基础(八):封装、继承、多态
- 计算机网络应用层题库
- Android app 后台被杀恢复
- 【工大SCIR】对话中的情感分析与生成简述
- 相似度计算之Jaccard系数
- vue.js 回显textarea入的空格和换行
- Android错误:unexpected text found in layout file
- 南都周刊 别了 老兵乔布斯
- Android studio 改app图标,名字及一键拨号
热门文章
- 组合数函数-快速提取所有可能的组合数
- ML基石_9_LinearRegression
- 便携式不锈钢管道焊接机器人_304不锈钢管居然可以发黑!?
- python 函数 默认参数
- 回归模型-线性回归算法
- Nvidia TX2 安装中文输入法
- 语义分割--End-to-End Instance Segmentation with Recurrent Attention
- 目标定位--Deep Self-Taught Learning for Weakly Supervised Object Localization
- [librosa]0.8版本后不再支持write_wav
- Linux: what's the difference between /proc/self and /proc/$$?