NumPy, Matplotlib, метод бисекций

from math import cos,exp
 
def dicho(f,a,b,eps):
    fa=f(a)
    fb=f(b)
    if fa*fb < 0:
        while (abs(a-b) > eps):
            c=0.5*(a+b)
            fc=f(c)
            if abs(fc)<eps:
                return c
            if fc*fa < 0:
                b=c
                fb=fc
            else:
                a=c
                fa=fc
        return 0.5*(a+b)
    else:
        return None
        
print(dicho(lambda x: exp(x)-2*cos(x),0,3,1.0e-6))

import numpy as np 
import matplotlib.pyplot as plt
 
def half_func(a, b, f):
    x = (a + b) / 2
    while np.abs(f(x)) >= eps:
        x = (a + b) / 2
        a, b = (a, x) if f(a) * f(x) < 0 else (x, b)
    return (a + b) / 2
 
def fun(x):
  return np.exp(x) - 2*np.cos(x)
 
 
eps = 1e-4
x = np.linspace(0, 1, 100)
y1 = np.exp(x)
y2 = 2*np.cos(x)
 
plt.xlabel('x')
plt.ylabel('y1, y2')
plt.grid()
plt.plot(x, y1, x, y2)
plt.show()
 
xo = half_func(0.4, 0.6, fun)
print(f'{xo:.5f}')
print(f'{fun(xo):.6f}')

def half_func(a, b, f):
    x = (a + b) / 2
    fa = fx = f(a)
    while np.abs(f(x)) >= eps:
        x = (a + b) / 2
        fx = f(x)
        if fa*fx < 0:
          b = x 
        else:
          a = x
          fa = fx
    return (a + b) / 2