二分法求解

import numpy as np #导入NumPy工具包
#定义几个不同的函数用于后续的二分法求解验证
f0=lambda x:2**x-4
f1=lambda x:np.log10(x)
f2=lambda x:x**0.3-1
f3=lambda x:x**2-x*2-1
# Python代码中的log函数如果没有设置底数的话，默认为e，即自然对数。
#文中其他非代码部分的log如果没有特别说明，一般也是指自然对数
#定义二分算法函数，其中a,b用于定义函数区间[a,b]，sec=(a+b)/2为区间中点
def bisection(a,b,fun):
    sec=(a+b)/2
    while True:
        if fun(sec)*fun(a)<0:
            b=sec
        elif fun(sec)*fun(a)>0:
            a=sec
        else:
            return sec
        least_sec=sec
        sec=(a+b)/2
        try:
            #给定一个精度ε，用于函数求解的近似计算
            if abs((sec-least_sec)/sec)<0.000005:
                return sec
        except ZeroDivisionError as e: #除数为0出错处理
            pass
    pass#输出上述几个函数基于给定精度的近似根
print(bisection(-5,5,f0))
print(bisection(0.5,5,f1))
print(bisection(0,15.0,f2))
print(bisection(-10.0,12.0,f3))