################################################################################
################################################################################
#   PERCOLATION
################################################################################
################################################################################

from ClassePile import *
import numpy as np, numpy.random as rd, matplotlib.pyplot as plt

###############################################################################
#   
################################################################################


p=.3
n=10
milieu=[[int(rd.rand()<p)*2 for j in range(n)] for i in range(n)]

plt.pcolormesh(np.array(milieu),cmap='Greys',edgecolors='k', linewidth=2)
ax = plt.gca();ax.invert_yaxis();ax.set_aspect('equal');plt.show()