Necesito un programa que realice el método de falsa posición. La función es (euler^3x) - 4 y el primer intervalo es [0 1]. Iterando hasta conseguir un porcentaje de error del 1% o menos.

def falsa_posicion(funcion, intervalo_inicial, intervalo_final, iteraciones):
    for i in range(iteraciones):
        x = random.randint(intervalo_inicial, intervalo_final)
        y = funcion(x)
        if y < 0:
            return x
    return -1

# Funcion de prueba
def funcion(x):
    return x**2