partie3_ex78.py

Created on April 22, 2022

289 Bytes

from math import exp

def f(x):
    return (-4 * x**2 - 10 * x + 8) * exp(-0.5 * x)

def sol(e):
    a = 0
    b = 1
    while b - a > e:
        m = (a + b) / 2
        p = f(a) * f(m)
        if p > 0:
            a = m
        else:
            b = m
    return a, b