import numpy as np
import matplotlib.pyplot as plt

def bisection(a,b,f,err):
    c = (a+b)/2
    while c-a > err:
        if f(c) == 0:
            return c
        elif f(a)*f(c) > 0:
            a = c
        else:
            b = c
        c = (a+b)/2
    return c

x = np.linspace(.5,8,1000)
y =  np.log(x) - np.sin(x) #Don't think we covered this in class about how to write In(x) and sin(x)
plt.plot(x,y)

x = np.linspace(2,4,1000)
y =  np.log(x) - np.sin(x) 
plt.plot(x,y)

def f(x):
    return np.log(x) - np.sin(x)

bisection(2,4,f,0.001)