python求方程的根_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程序设计有所帮助。

python求方程的根_python计算方程式根的方法相关推荐

python牛顿法求方程的根_python计算方程式根的方法如何用python计算三元方程
已知一元二次方程的3个参数a,b,c,编写python程序分享import math a = float(raw_input('Enter coefficient a')) b = float(raw ...
python编程求1!+2!+…+n!_python计算阶乘和的方法(1!+2!+3!+...+n!)
python计算阶乘和的方法(1!+2!+3!+...+n!) 更新时间:2019年02月01日 08:50:15 作者:一米阳光里的晴天娃娃今天小编就为大家分享一篇python计算阶乘和的方法(1 ...
python求平方根的代码_Python求解平方根的方法
本文实例讲述了Python求解平方根的方法.分享给大家供大家参考.具体如下: 主要通过SICP的内容改写而来.基于newton method求解平方根.代码如下: #!/usr/bin/python ...
写一个判断方程根的C语言程序,写出求方程ax2+bx+c=0的根的算法,画出相应的程序框图,并要求输出它的实根....
写出求方程ax2+bx+c=0的根的算法,画出相应的程序框图,并要求输出它的实根. 答案思路分析:主要考查含参问题的讨论方法及条件嵌套结构的应用.输入三个实数a,b,c,先判断a是否为零,当a≠0时 ...
python求数组平均值_用python求一个数组的和与平均值的实现方法
用python求一个数组的和与平均值的实现方法如下所示: # coding = GBK a =[1,2,3,4,5] sum=0 b = len(a) print("这个数组的长度为:&q ...
python求n的阶乘_python求n的阶乘
阶乘是基斯顿·卡曼(Christian Kramp,1760-1826)于1808年发明的运算符号,是数学术语.一个正整数的阶乘(factorial)是所有小于及等于该数的正整数的积,并且0的阶乘为1 ...
python求50的阶乘_python中求阶乘
广告关闭腾讯云11.11云上盛惠 ,精选热门产品助力上云,云服务器首年88元起,买的越多返的越多,最高返5000元! 我如何去计算python中的一个整数的阶乘?... 写一个猜数字的游戏,预先设定 ...
python简单小案例列表_python计算列表内各元素的个数实例
python计算列表内各元素的个数实例如下所示: list = [1,2,3,4,5,6,7,5,4,3,2,12] set = set(list) dict = {} for item in se ...
python 求平面两点距离_Python求平面内点到直线距离的实现
近期遇到个问题,需要计算平面内点到直线的距离,发现数学知识都还给老师了,度娘后找到计算方法,特此记录. 点到直线的计算公式: 通过公式推导,得到信息: A:直线斜率 B:固定值-1 C:直线截距b 转 ...
python求两数最大公因数_Python求两个数的最大公约数
Python求两个数的最大公约数一.求最大公约数算法: 1. 整数A对整数B进行取整, 余数用整数C来表示举例: C = A % B 2. 如果C等于0,则C就是整数A和整数B的最大公约数 3. ...

python求方程的根_python计算方程式根的方法

python求方程的根_python计算方程式根的方法相关推荐

最新文章

热门文章