def newton_côtes(f, domain: tuple[float, float], resolution: int = 1000) -> float:
    """
    >>> from math import pi, sin
    >>> newton_côtes(sin, (0, pi/2))
    """
    step = (domain[1]-domain[0])/resolution
    return step * sum(
                f(domain[0] + k * step) for k in range(resolution)
            )

def precision(methode, n):
    """
    >>> 
    """