python计算方程式根的方法
作者：songguo 发布时间：2023-08-09 09:54:47
标签：python,计算
本文实例讲述了python计算方程式根的方法。分享给大家供大家参考。具体实现方法如下：
''' roots = polyRoots(a).
Uses Laguerre's method to compute all the roots of
a[0] + a[1]*x + a[2]*x^2 +...+ a[n]*x^n = 0.
The roots are returned in the array 'roots',
'''
from evalPoly import *
from numpy import zeros,complex
from cmath import sqrt
from random import random
def polyRoots(a,tol=1.0e-12):
def laguerre(a,tol):
x = random()
# Starting value (random number)
n = len(a) - 1
for i in range(30):
p,dp,ddp = evalPoly(a,x)
if abs(p) < tol: return x
g = dp/p
h = g*g - ddp/p
f = sqrt((n - 1)*(n*h - g*g))
if abs(g + f) > abs(g - f): dx = n/(g + f)
else: dx = n/(g - f)
x = x - dx
if abs(dx) < tol: return x
print 'Too many iterations'
def deflPoly(a,root): # Deflates a polynomial
n = len(a)-1
b = [(0.0 + 0.0j)]*n
b[n-1] = a[n]
for i in range(n-2,-1,-1):
b[i] = a[i+1] + root*b[i+1]
return b
n = len(a) - 1
roots = zeros((n),dtype=complex)
for i in range(n):
x = laguerre(a,tol)
if abs(x.imag) < tol: x = x.real
roots[i] = x
a = deflPoly(a,x)
return roots
raw_input("\nPress return to exit")
希望本文所述对大家的Python程序设计有所帮助。
0
投稿猜你喜欢
5月3日晚，央视在《新闻联播》前播放了B站青年宣言片《后浪》，这是B站首次登陆央视黄金时段，今天在朋友圈陆续看到相关的视频。最早用B站的同学- 1. 匿名函数匿名函数（anonymous function）是指未与任何标识符绑定的函数，多用在functional programmin
- 计算机一般来说是需要定期的清理，系统的内存不能无限延伸，同时有一些不需要的文件也可以得以清除掉。有些人会使用os.remove来进行文件的清
- yaml 文件的应用场景与格式介绍yaml 文件的应用场景yaml其实也类似于 json、txt ，它们都属于一种文本格式。在我们的实际工作
- 这带来了很多好处，但是也要求学习一些新的概念、命令、使用和管理任务。所以在投入你的生产系统之前，看看它解决什么、管理什么，以及对它的正反面评
- JavaScript： <script type="text/javascript"> var level1
- 为了吸引更多的用户，设计好一个分享海报还是很有必要的。而小程序要生成一个海报还是有点坑的，下面分享下我们打卡小程序的一些经验。分享海报的效果
- 今天在使用MySQL时却不知如何处理，插入记录后不知怎样获得刚刚插入的id，查过文档后发现了select last_insert_id()，
- 1.手动协程操作：# pip install geventfrom greenlet import greenletdef test():
- 本文实例讲述了mysql外键基本功能与用法。分享给大家供大家参考，具体如下：本文内容：什么是外键外键的增加外键的修改和删除外键的约束模式首发
- 本文实例讲述了PHP在线打包下载功能实现方法。分享给大家供大家参考，具体如下：昨天晚上，为了弄这个打包下载的事，弄的事焦头烂额。有几个问题，
- 腾讯“月捐计划”倡导爱心人士，通过每月小额捐款的形式，长期关注和支持公益项目。并和亿万爱心网友一起，每人每月1份爱，点滴付出，汇成爱海，形成
- 一.引言先看获取到的效果拍摄时间：2021:12:18 16:22:13照片拍摄地址:('内蒙古自治区包头市昆都仑区', &
- 在python的类中会使用到“self”，本文主要是谈谈类中的self和变量的关系。先看示例代码I。
- Windows环境下一、开启 Imagick 扩展1、安装PHP扩展：Imagick，下载地址 https://pecl.php.net/p
- 如下所示：import arcpy... from arcpy import env... env.workspace="C:\\
- 本文实例讲述了Flask-Mail用法。分享给大家供大家参考，具体如下：很多类型的应用程序都需要在特定事件发生时提醒用户，而常用的通信方法是
- go和golang之间有区别吗？go和golang之间没有区别，go就是golang，全名是golang，简称go语言。Go（又称 Gola
- 由于在Python2 中的默认编码为ASCII，但是在Python3中的默认编码为UTF-8。问题：所以在使用np.load(det.npy
- 什么多态：同一事物有多种形态为何要有多态=》多态会带来什么样的特性，多态性多态性指的是可以在不考虑对象具体类型的情况下而直接使用对象多态指的
