python计算方程式根的方法
更新时间:2015年05月07日 11:19:44 作者:songguo
这篇文章主要介绍了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-typing: 类型标注与支持 Any类型详解
这篇文章主要介绍了Python-typing: 类型标注与支持 Any类型详解,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧2021-05-05python GUI框架pyqt5 对图片进行流式布局的方法(瀑布流flowlayout)
这篇文章主要介绍了python GUI框架pyqt5 对图片进行流式布局的方法(瀑布流flowlayout),本文通过实例代码给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友可以参考下2020-03-03Linux下将Python的Django项目部署到Apache服务器
这篇文章主要介绍了Python的Django项目部署到Apache服务器上的要点总结,文中针对的是wsgi连接方式,需要的朋友可以参考下2015-12-12
最新评论