Python查找函数f（x）=0根的解决方法
更新时间：2015年05月07日 11:28:51 作者：songguo
这篇文章主要介绍了Python查找函数f（x）=0根的解决方法,涉及Python数学运算函数求解的相关技巧,需要的朋友可以参考下
本文实例讲述了Python查找函数f（x）=0根的解决方法。分享给大家供大家参考。具体实现方法如下：
''' 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程序设计有所帮助。
您可能感兴趣的文章:
相关文章
Boosting是机器学习的三大框架之一。Boost也被称为增强学习或提升法，其中典型的代表算法是AdaBoost算法。本文介绍了AdaBoost算法及python实现，感兴趣的可以学习一下2022-10-10
数据可视化是任何探索性数据分析或报告的关键步骤，目前有许多非常好的商业智能工具，比如Tableau、googledatastudio和PowerBI等，本文就详细的进行对比，感兴趣的可以了解一下2021-06-06
这篇文章主要介绍了Python 高效办公第三方库 openpyxl 的安装，本文介绍了什么是库及库的分类，本文结合示例代码给大家介绍的非常详细，对大家的学习或工作具有一定的参考借鉴价值,需要的朋友可以参考下2022-12-12
这篇文章主要为大家详细介绍了python3使用matplotlib绘制条形图，文中示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下2019-03-03
最近在做一个项目，这个项目需要使用网络爬虫从特定网站上爬取数据，于是乎，我打算写一个爬虫系列的文章，与大家分享如何编写一个爬虫。下面这篇文章给大家介绍了python爬虫基本知识，感兴趣的朋友一起看看吧2018-03-03
Python脚本,标识符,变量使用,脚本语句,注释,模块引用详解
这篇文章主要为大家详细介绍了Python脚本,标识符,变量使用,脚本语句,注释,模块引用，文中示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下2022-02-02
这篇文章主要介绍了Python-基础-简介入门的相关资料,需要的朋友可以参考下2014-08-08
这篇文章主要为大家详细介绍了python实现图片识别汽车功能，具有一定的参考价值，感兴趣的小伙伴们可以参考一下2018-11-11
今天小编就为大家分享一篇关于Python中extend和append的区别讲解，小编觉得内容挺不错的，现在分享给大家，具有很好的参考价值，需要的朋友一起跟随小编来看看吧2019-01-01
这篇文章主要介绍了python双向链表原理与实现方法,结合实例形式分析了Python双向链表的定义、以及节点的判断、遍历、添加、删除等相关操作技巧,需要的朋友可以参考下2019-12-12
最新评论