import numpy as np
import math as mt
from sympy import *

x = Symbol('x')

formula = x**4 - 625

a = formula.diff(x1,2)
value = a.evalf(subs={x:2})

epsilon = 0.00001
a = 0
b = 3
def function(x):
    result = x**4 - 625
    return result

x1 = a - function(a)/(function(b) - function(a))*(b-a)