Simulation to explore https://twitter.com/juliagalef/status/887414654401630208

neqlt.py

import random 

def hist(x):
    h = [0 for i in range(11)]
    band = x[N-1]/10.0

    for i in range(N):
        dec = (int)(x[i]/band)
        h[dec] = h[dec] + 1
        
    return h
    
N = 100 # population size
W = 100 # initial wealth
I = 1000 # number of iterations

x = [W for i in range(N)]

for i in range(I):
    x.sort()
    print i
    print x
    print hist(x)
    print "-" * 60
    for j in range(N):
        if x[j] > 0:
            k = random.randint(0, N-1)
            x[k] = x[k] + 1
            x[j] = x[j] - 1