tp3 PORC MULD

from pylab import *

def mydicho(f, a, b, N):
    alpha = a
    beta = b
    x = zeros(N)

    for n in range(0,N):
        val = (alpha + beta)/2.
        prod = f(alpha) * f(val)

        if prod<0:
            beta = val
        else:
            alpha = val
        x[n] = val
    return x

def f(x):
    return exp(x)-2


def mydicho2(f, a, b, N, eps):
    alpha = a
    beta = b
    x = []
    val = (alpha + beta)/2.
    while abs(f(val))>eps:
        val = (alpha + beta)/2.
        prod = f(alpha) * f(val)

        if prod<0:
            beta = val
        else:
            alpha = val
        x.append(val)
    return x

def mypointfixe(Xinit, N, f):
    x = zeros(N)
    x[0] = Xinit

    for n in range(1,N):
        x[n] = f(x[n-1])

    return x