def estabilitzacio(y0,k,epsilon):
    y1 = y0
    y2 = y1 + k*y1*(1 - y1)
    n = 1
    while abs(y2 - y1) >= epsilon:
        y1 = y2
        y2 = y1 + k*y1*(1  - y1)
        n = n + 1
    return n