transfer_money.ipynb

import numpy as np
import seaborn as sns
import pandas as pd
import matplotlib.pyplot as plt
from IPython.core.display import display

%matplotlib inline

N = 10
M = 10

def transfer_money(moneys):
    to_moneys = np.random.randint(0, N, N)
    for i in range(N):
        if moneys[i] > 0:
            moneys[i] -= 1
            moneys[to_moneys[i]] += 1

moneys = np.ones((N, )) * M

for i in range(10001):
    if i % 1000 == 0:
        ax = sns.distplot(moneys, kde=False, bins=20)
        ax.set_title("i=%d" % i)
        ax.set_xlabel("money")
        ax.set_ylabel("person count")
        plt.show()
    transfer_money(moneys)