def est_pair(n): return n%2==0 def mélange(n): if est_pair(n): return n//2 else: return 3*n+1 def Collatz(u0,n): u = [u0] for k in range(n): u.append(mélange(u[-1])) return u