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)技巧,需要的朋友可以參考下
本文實例講述了Python查找函數(shù)f（x）=0根的解決方法。分享給大家供大家參考。具體實現(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è)計有所幫助。
您可能感興趣的文章:
相關(guān)文章
Boosting是機(jī)器學(xué)習(xí)的三大框架之一。Boost也被稱為增強(qiáng)學(xué)習(xí)或提升法，其中典型的代表算法是AdaBoost算法。本文介紹了AdaBoost算法及python實現(xiàn)，感興趣的可以學(xué)習(xí)一下2022-10-10
淺談哪個Python庫才最適合做數(shù)據(jù)可視化
數(shù)據(jù)可視化是任何探索性數(shù)據(jù)分析或報告的關(guān)鍵步驟，目前有許多非常好的商業(yè)智能工具，比如Tableau、googledatastudio和PowerBI等，本文就詳細(xì)的進(jìn)行對比，感興趣的可以了解一下2021-06-06
這篇文章主要介紹了Python?高效辦公第三方庫?openpyxl?的安裝，本文介紹了什么是庫及庫的分類，本文結(jié)合示例代碼給大家介紹的非常詳細(xì)，對大家的學(xué)習(xí)或工作具有一定的參考借鑒價值,需要的朋友可以參考下2022-12-12
這篇文章主要為大家詳細(xì)介紹了python3使用matplotlib繪制條形圖，文中示例代碼介紹的非常詳細(xì)，具有一定的參考價值，感興趣的小伙伴們可以參考一下2019-03-03
最近在做一個項目，這個項目需要使用網(wǎng)絡(luò)爬蟲從特定網(wǎng)站上爬取數(shù)據(jù)，于是乎，我打算寫一個爬蟲系列的文章，與大家分享如何編寫一個爬蟲。下面這篇文章給大家介紹了python爬蟲基本知識，感興趣的朋友一起看看吧2018-03-03
Python腳本,標(biāo)識符,變量使用,腳本語句,注釋,模塊引用詳解
這篇文章主要為大家詳細(xì)介紹了Python腳本,標(biāo)識符,變量使用,腳本語句,注釋,模塊引用，文中示例代碼介紹的非常詳細(xì)，具有一定的參考價值，感興趣的小伙伴們可以參考一下2022-02-02
這篇文章主要介紹了Python-基礎(chǔ)-簡介入門的相關(guān)資料,需要的朋友可以參考下2014-08-08
這篇文章主要為大家詳細(xì)介紹了python實現(xiàn)圖片識別汽車功能，具有一定的參考價值，感興趣的小伙伴們可以參考一下2018-11-11
今天小編就為大家分享一篇關(guān)于Python中extend和append的區(qū)別講解，小編覺得內(nèi)容挺不錯的，現(xiàn)在分享給大家，具有很好的參考價值，需要的朋友一起跟隨小編來看看吧2019-01-01
這篇文章主要介紹了python雙向鏈表原理與實現(xiàn)方法,結(jié)合實例形式分析了Python雙向鏈表的定義、以及節(jié)點(diǎn)的判斷、遍歷、添加、刪除等相關(guān)操作技巧,需要的朋友可以參考下2019-12-12