Python查找函數(shù)f（x）=0根的解決方法
更新時間：2015年05月07日 11:28:51 作者：songguo
這篇文章主要介紹了Python查找函數(shù)f（x）=0根的解決方法,涉及Python數(shù)學(xué)運(yùn)算函數(shù)求解的相關(guān)技巧,需要的朋友可以參考下
本文實(shí)例講述了Python查找函數(shù)f（x）=0根的解決方法。分享給大家供大家參考。具體實(shí)現(xiàn)方法如下：
''' root = ridder(f,a,b,tol=1.0e-9).
Finds a root of f(x) = 0 with Ridder's method.
The root must be bracketed in (a,b).
'''
import error
from math import sqrt
def ridder(f,a,b,tol=1.0e-9):
fa = f(a)
if fa == 0.0: return a
fb = f(b)
if fb == 0.0: return b
if fa*fb > 0.0: error.err('Root is not bracketed')
for i in range(30):
# Compute the improved root x from Ridder's formula
c = 0.5*(a + b); fc = f(c)
s = sqrt(fc**2 - fa*fb)
if s == 0.0: return None
dx = (c - a)*fc/s
if (fa - fb) < 0.0: dx = -dx
x = c + dx; fx = f(x)
# Test for convergence
if i > 0:
if abs(x - xOld) < tol*max(abs(x),1.0): return x
xOld = x
# Re-bracket the root as tightly as possible
if fc*fx > 0.0:
if fa*fx < 0.0: b = x; fb = fx
else: a = x; fa = fx
else:
a = c; b = x; fa = fc; fb = fx
return None
print 'Too many iterations'
希望本文所述對大家的Python程序設(shè)計(jì)有所幫助。
您可能感興趣的文章:
相關(guān)文章
django 鏈接多個數(shù)據(jù)庫 并使用原生sql實(shí)現(xiàn)
這篇文章主要介紹了django 鏈接多個數(shù)據(jù)庫 并使用原生sql實(shí)現(xiàn)，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧2020-03-03
程序在運(yùn)行時,如果python解釋器遇到一個錯誤,會停止程序的執(zhí)行,并且提示一些錯誤信息,這就是異常,下面這篇文章主要給大家介紹了關(guān)于Python基礎(chǔ)教程之錯誤和異常的處理方法,需要的朋友可以參考下2022-05-05
pyspark給dataframe增加新的一列的實(shí)現(xiàn)示例
這篇文章主要介紹了pyspark給dataframe增加新的一列的實(shí)現(xiàn)示例,文中通過示例代碼介紹的非常詳細(xì)，對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值，需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2020-04-04
這篇文章是一個極簡的Python入門指引、原文在網(wǎng)絡(luò)上人氣相當(dāng)高，將一些基本知識用實(shí)例代碼說明，需要的朋友可以參考下2015-04-04
Python字典添加，刪除，查詢等相關(guān)操作方法詳解
這篇文章主要介紹了Python字典添加，刪除，查詢等相關(guān)操作方法詳解,需要的朋友可以參考下2020-02-02
Python 3.7中一個令人興奮的新特性是 data classes 。這篇文章主要介紹了Python3.7 新特性之dataclass裝飾器,需要的朋友可以參考下2019-05-05
Python元類編程實(shí)現(xiàn)一個簡單的ORM
本文主要介紹了Python元類編程實(shí)現(xiàn)一個簡單的ORM，文中通過示例代碼介紹的非常詳細(xì)，對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值，需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2023-03-03
這篇文章主要介紹了Python之web模板應(yīng)用，web模板可以更加靈活和方便的控制HTML的顯示，非常具有實(shí)用價值,有需要的小伙伴可以參考下2017-12-12
這篇文章主要介紹了python中asyncio異步編程學(xué)習(xí)，內(nèi)部就是基于協(xié)程實(shí)現(xiàn)的異步編程，如果想研究異步編程的同學(xué)，要仔細(xì)看哦2021-04-04
Python必備技巧之Pandas數(shù)據(jù)合并函數(shù)
Pandas中一共有五個數(shù)據(jù)合并函數(shù)，分別為：concat、append、merge、join、combine，本文詳細(xì)講解這五個函數(shù)的使用方法，需要的可以參考一下2022-03-03