#!/usr/bin/python3

import numpy as np
import matplotlib.pyplot as plt
import math as math

#Q4a
def bracketing(f,a,b,n):
    brackets= []
    for i in range(0,n):
        if f(a+(b-a)*i/n)*f(a+(b-a)*(i+1)/n)<0:
            brackets.append([a+(b-a)*i/n,a+(b-a)*(i+1)/n])
    return brackets


def f(x):
    return x**4-x-1

print(bracketing(f,-2,2,1000))

#Q4b
def secant(f,x,y,tol):
    while abs(x-y) > tol:
        temp = x
        x = (y*f(x)-x*f(y))/(f(x)-f(y))
        y = temp
    return x