geothmetic meandian (https://xkcd.com/2435/)

geothmetic_meandian.py

import statistics


def has_converged(xs):
  return all(x == xs[0] for x in xs)


def f(xs):
  return [
      statistics.mean(xs),
      statistics.median(xs),
      statistics.geometric_mean(xs)
  ]


def gmdn(xs):
  while not has_converged(xs):
    print('.', end='')
    xs = f(xs)
  print('\n')
  return xs


print(gmdn([1, 1, 2, 3, 5]))