import math
import numpy as np
import matplotlib.pyplot as plt
def func2 (x):
    return((2*x**2)-5*x-12)

def false_position(fx, xl, xu, iter_num):
    ea_xr = []
    for i in range(iter_num):
        xr = (xl*fx(xu) - xu*fx(xl)/(fx(u)- fx(xl))
        ea_xr.append(xr)
        if (fx(xl)*fx(xr) < 0):
              xu=xr
        else:
              xl = xr
        return ea_xr
false_position(func2,0,7,20)
              
def bisect(fc, xl, xu, iternum):
              estimates_xr = []
              for i in range(iternum):
              xr = (xl + xu)/2 
              estimates_xr.append(xr)
              if (fx(xl)*fx(xr)
                 xu = xr
              else:
                  xl = xr
            return estimates_xr
    bisect(func2, 0, 7, 20)
                  
def new_raph1(func, dfunc, xinit, ea, imax=20):
                  iter=0
                  xold=xinit
                  roots=[]
                  errors=[]
                  while (iter < imax):
                  xnew= xold- (func(xold)/dfunc(xold))
                  roots.appends(xnew)
                  
                  err = (np.abs((xnew - xold)/xnew)) * 100
                  xold = xnew
                  errors.append(err)
                  if (err<ea):
                iter += 1
                  return roots,errors
      x,y = new_raph(func2, der2, 0, 0)
      print("roots = ", np.round(x,7), ,"\n", ,"errors = ", np.round(y,7))