Bak_Sneppen model

Bak_Sneppen.py

def bak_sneppen(data, steps=10):
    N = len(data)

    for i in range(steps*N):
        pos = np.argmin(data) # Find the position of the minimum value
        
        # Replace the minimum valuea and both it's neighbours
        data[(pos-1) % N] = np.random.random()
        data[pos] = np.random.random()
        data[(pos+1) % N] = np.random.random()

    return data