在使用Python的工程中的数值方法,第2版,作者:Jaan Kiusalaas,我写了第146页的相同模块,它使用二分法计算根,f(x)= 0:
from math import log,ceil
import error
def bisection2(f,x1,x2,switch=0,tol=1.0e-9):
f1 = f(x1)
if f1 == 0.0: return x1
f2 = f(x2)
if f2 == 0.0: return x2
if f1*f2 > 0.0: error.err('Root is not bracketed')
n = ceil(log(abs(x2 - x1)/epsilon)/log(2.0))
for i in range(n):
x3 = 0.5*(x1 + x2); f3 = f(x3)
if (switch == 1) and (abs(f3) > abs(f1)) \
(abs(f3) > abs(f2)):
return None
if f3 == 0.0:return x3
if f2*f3 < 0.0:
x1 = x3; f1 = f3
else:
x2 = x3; f2 = f3
return (x1 + x2)/2.0
我遇到以下错误:
importError Traceback(最近一次调用最后一次)execfile中的/usr/lib/python2.7/dist-packages/IPython/utils/py3compat.pyc(fname,* where)173 else:174 filename = fname - > 175 builtin.execfile (文件名,*在哪里)
/home/uwhpsc/Desktop/bisection2/bisection2.py in()1来自数学导入日志,ceil ----> 2导入错误3 4 def bisection2(f,x1,x2,switch = 0,tol = 1.0e- 9):5 f1 = f(x1)
ImportError:没有名为error的模块
请有人可以告诉我如何解决这个问题?
要使此代码有效,您可以删除说:
import error
在线上说:
if f1*f2 > 0.0: error.err('Root is not bracketed')
改为:
if f1*f2 > 0.0: quit('Root is not bracketed')