Implementeer de halveringsmethode in Python, d.w.z., definieer de volgende Python functie:

def bisection_solve(f, a, b, tol=0.001, maxiter=100):
    """
    Find the zero of the function f between a and b using the
    Bisection Method with tolerance tol (default: 0.001) and 
    maximum number of iterations equal to maxiter (default: 100)
    """

Pas de functie bisection_solve toe op de veelterm \(x^3+2x-1\) op het interval \([0,1]\) een keer toe met de standaardwaarde van de tolerantie en met een tolerantie van \(10^{-5}\).

Het exacte en benaderde nulpunt binnen het interval is overigens \[\frac{1}{6}\,\sqrt[3]{108+12\sqrt {177}}-4{\frac {1}{\sqrt[3]{108+12\sqrt {177}}}}\approx 0.4533976515\]