def conv_mp21(a0, a1, ne):
    a2 = ccl_term(a0, a1)
    eps = 1/10**ne
    # while abs(a2-a1)/2.0 > eps/2.0:
    while abs(a2-a1)/2.0 > eps:
        a0 = a1
        a1 = a2
        a2 = ccl_term(a0, a1)
    return round(a2, ne)

def ccl_term(a0, a1):
    return (2*a1+a0) / 3.0