Home License -- for personal use only. Not for government, academic, research, commercial, or other organizational use. 13-May-2025 16:30:27 percolation_simulation_test(): MATLAB/Octave version 9.11.0.2358333 (R2021b) Update 7 percolation_simulation() simulates percolation. A rectangular region is formed by a grid of M x N squares. Each square may be porous or solid. We are interested in paths of porous squares connecting the top and bottom, or the left and right boundaries. Grid uses 32 rows x 32 columns Prescribed probability of occupation = 0.57 plist = 2 2 plist = 2 2 2 4 plist = 2 2 3 4 plist = 2 2 4 4 plist = 2 2 5 4 plist = 2 2 5 4 5 3 plist = 2 2 5 4 6 4 plist = 2 2 5 4 7 4 plist = 2 2 5 4 8 4 plist = 2 2 5 4 8 4 7 3 8 6 plist = 2 2 5 4 8 4 7 3 9 6 plist = 2 2 5 4 8 4 7 3 9 6 9 7 plist = 2 2 5 4 8 4 7 3 9 6 10 7 plist = 2 2 5 4 8 4 7 3 9 6 10 7 9 6 8 9 10 9 plist = 2 2 5 4 8 4 7 3 9 6 10 7 10 6 plist = 2 2 5 4 8 4 7 3 9 6 10 7 10 6 10 5 plist = 2 2 5 4 8 4 7 3 9 6 10 7 10 6 10 4 11 6 plist = 2 2 5 4 8 4 7 3 9 6 10 7 10 6 10 4 11 6 11 7 plist = 2 2 5 4 8 4 7 3 9 6 10 7 10 6 10 4 11 6 12 6 plist = 2 2 5 4 8 4 7 3 9 6 10 7 10 6 10 4 11 6 13 6 plist = 2 2 5 4 8 4 7 3 9 6 10 7 10 6 10 4 11 6 13 6 13 4 plist = 2 2 5 4 8 4 7 3 9 6 10 7 10 6 10 4 11 6 13 6 13 4 11 3 13 3 plist = 2 2 5 4 8 4 7 3 9 6 10 7 10 6 10 4 11 6 13 6 13 4 11 3 14 3 plist = 2 2 5 4 8 4 7 3 9 6 10 7 10 6 10 4 11 6 13 6 13 4 11 3 14 3 14 4 plist = 2 2 5 4 8 4 7 3 9 6 10 7 10 6 10 4 11 6 13 6 13 4 11 3 14 3 15 4 plist = 2 2 5 4 8 4 7 3 9 6 10 7 10 6 10 4 11 6 13 6 13 4 11 3 14 3 15 4 14 3 15 5 plist = 2 2 5 4 8 4 7 3 9 6 10 7 10 6 10 4 11 6 13 6 13 4 11 3 14 3 15 4 14 3 15 4 16 6 plist = 2 2 5 4 8 4 7 3 9 6 10 7 10 6 10 4 11 6 13 6 13 4 11 3 14 3 15 4 14 3 15 4 16 6 14 7 16 7 plist = 2 2 5 4 8 4 7 3 9 6 10 7 10 6 10 4 11 6 13 6 13 4 11 3 14 3 15 4 14 3 15 4 16 6 14 7 16 6 15 9 17 10 plist = 2 2 5 4 8 4 7 3 9 6 10 7 10 6 10 4 11 6 13 6 13 4 11 3 14 3 15 4 14 3 15 4 16 6 14 7 16 6 12 11 14 12 plist = 2 2 5 4 8 4 8 3 plist = 2 2 5 4 8 4 8 3 8 1 plist = 2 2 5 4 8 4 9 3 plist = 3 2 plist = 3 2 3 1 plist = 3 2 4 1 plist = 2 9 4 9 plist = 2 9 4 8 5 10 plist = 2 9 4 8 6 10 plist = 2 9 4 8 6 10 6 11 plist = 2 9 4 8 6 10 6 11 4 12 6 12 plist = 2 9 4 8 6 10 6 11 4 12 6 12 6 13 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 7 15 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 6 24 8 24 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 6 24 8 24 8 25 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 6 24 8 24 8 25 8 26 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 6 24 8 24 8 25 9 26 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 6 24 8 24 8 25 9 26 9 24 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 6 24 8 24 8 25 10 28 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 6 24 8 24 8 25 11 28 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 6 24 8 24 8 25 12 28 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 6 24 8 24 8 25 12 28 11 27 12 29 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 6 24 8 24 8 25 12 28 11 27 13 29 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 6 24 8 24 8 25 12 28 11 27 13 29 12 28 13 30 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 6 24 8 24 8 25 12 28 11 27 13 29 12 28 13 30 11 31 13 31 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 6 24 8 24 8 25 12 28 11 27 13 29 12 28 13 30 11 31 12 28 14 28 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 6 24 8 24 8 25 12 28 11 27 13 29 12 28 13 30 11 31 12 28 15 28 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 6 24 8 24 8 25 12 28 11 27 13 29 12 28 13 30 11 31 12 28 16 31 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 5 23 4 22 5 21 3 23 2 22 3 21 2 23 2 24 2 25 4 25 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 5 23 4 22 5 21 3 23 2 22 3 21 2 23 2 24 2 25 4 25 2 26 4 27 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 5 23 4 22 5 21 3 23 2 22 3 21 2 23 2 24 2 25 4 25 2 26 4 27 2 28 4 28 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 5 23 4 22 5 21 3 23 2 22 3 21 2 23 2 24 2 25 4 25 2 26 4 27 2 28 4 28 2 29 4 29 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 5 23 4 22 5 21 3 23 2 22 3 21 2 23 2 24 2 25 4 25 2 26 4 27 2 28 4 28 2 29 4 29 4 30 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 5 23 4 22 5 21 3 23 2 22 3 21 2 23 2 24 2 25 4 25 2 26 4 27 2 28 4 28 2 29 4 29 5 30 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 5 23 4 22 5 21 3 23 2 22 3 21 2 23 2 24 2 25 4 25 2 26 4 27 2 28 4 28 2 29 4 29 5 30 5 29 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 5 23 4 22 5 21 3 23 2 22 3 21 2 23 2 24 2 25 4 25 2 26 4 27 2 28 4 28 2 29 4 29 5 30 5 29 5 27 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 5 23 4 22 5 21 3 23 2 22 3 21 2 23 2 24 2 25 4 25 2 26 4 27 2 28 4 28 2 29 4 29 5 30 5 29 6 27 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 5 23 4 22 5 21 3 23 2 22 3 21 2 23 2 24 2 25 4 25 2 26 4 27 2 28 4 28 2 29 4 29 5 30 5 29 7 28 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 5 23 4 22 5 21 3 23 2 22 3 21 2 23 2 24 2 25 4 25 2 26 4 27 2 28 4 28 2 29 4 29 5 30 5 29 7 28 5 29 7 29 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 5 23 4 22 5 21 3 23 2 22 3 21 2 23 2 24 2 25 4 25 2 26 4 27 2 28 4 28 2 29 4 29 5 30 5 29 7 28 5 29 8 29 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 5 23 4 22 5 21 3 23 2 22 3 21 2 23 2 24 2 25 4 25 2 26 4 27 2 28 4 28 2 29 4 29 5 30 5 29 7 28 6 31 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 5 23 4 22 5 21 3 23 2 22 3 21 2 23 2 24 2 25 5 25 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 5 23 4 22 5 21 3 23 2 22 4 19 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 5 23 4 22 5 21 3 23 2 22 4 19 2 18 2 17 4 17 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 8 22 6 23 5 23 4 22 6 20 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 9 19 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 9 19 9 18 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 9 19 10 18 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 9 19 10 18 10 17 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 9 19 10 18 10 16 11 18 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 9 19 10 18 10 16 12 18 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 9 19 10 18 10 16 12 18 12 19 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 9 19 10 18 10 16 12 18 13 18 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 7 14 6 18 8 18 6 19 8 19 6 20 8 20 8 21 6 22 9 19 10 18 11 16 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 8 14 plist = 2 9 4 8 6 10 6 11 4 12 6 12 7 13 6 12 7 14 5 15 8 15 9 14 plist = 2 9 5 8 plist = 2 9 6 8 plist = 13 1 plist = 14 16 plist = 15 16 plist = 15 16 15 17 plist = 15 16 16 16 plist = 15 16 16 16 16 15 plist = 15 16 16 16 17 16 plist = 15 16 16 16 18 17 plist = 13 24 15 24 plist = 13 24 16 24 plist = 13 24 16 24 16 23 plist = 13 24 16 24 16 23 16 22 plist = 13 24 16 24 16 23 17 22 plist = 13 24 16 24 16 23 17 22 17 25 plist = 13 24 16 24 16 23 17 22 17 25 17 26 plist = 13 24 16 24 16 23 17 22 17 25 18 26 plist = 13 24 16 24 16 23 17 22 17 25 18 26 17 25 18 29 plist = 13 24 16 24 16 23 17 22 17 25 18 26 17 25 18 29 18 30 plist = 13 24 16 24 16 23 17 22 17 25 18 26 17 25 18 29 19 30 plist = 13 24 16 24 16 23 17 22 17 25 18 26 17 25 18 29 19 30 19 29 plist = 13 24 16 24 16 23 17 22 17 25 18 26 17 25 18 29 19 30 20 29 plist = 13 24 16 24 16 23 17 22 17 25 18 26 17 25 18 29 19 30 20 29 19 28 20 31 plist = 13 24 16 24 16 23 17 22 17 25 18 26 17 25 18 29 19 30 20 29 19 28 20 31 18 32 20 32 plist = 13 24 16 24 16 23 18 21 plist = 13 24 16 24 16 23 18 21 16 20 18 20 plist = 13 24 16 24 16 23 18 21 16 20 19 20 plist = 13 24 16 24 16 23 18 21 16 20 20 20 plist = 13 24 16 24 16 23 18 21 16 20 20 19 21 21 plist = 13 24 16 24 16 23 18 21 16 20 20 19 21 21 19 22 21 22 plist = 13 24 16 24 16 23 18 21 16 20 20 19 21 21 19 22 22 22 plist = 13 24 16 24 16 23 18 21 16 20 20 19 21 21 19 22 22 22 21 21 22 23 plist = 13 24 16 24 16 23 18 21 16 20 20 19 21 21 19 22 22 22 21 21 22 23 19 24 19 25 21 26 plist = 13 24 16 24 16 23 18 21 16 20 21 18 plist = 13 24 16 24 16 23 18 21 16 20 22 18 plist = 13 24 16 24 16 23 18 21 16 20 23 18 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 19 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 24 20 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 25 20 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 26 22 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 27 23 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 28 23 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 28 23 27 22 26 25 28 25 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 28 23 27 22 26 25 29 25 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 28 23 27 22 26 25 30 25 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 28 23 27 22 26 25 30 25 30 26 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 28 23 27 22 26 25 30 25 30 25 31 27 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 28 23 27 22 26 25 30 25 30 25 31 27 31 28 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 28 23 27 22 26 25 30 25 30 25 31 27 31 28 31 29 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 28 23 27 22 26 25 30 25 30 25 31 27 31 28 31 29 31 30 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 28 23 27 22 26 25 30 25 30 25 31 27 31 28 31 29 31 30 29 31 31 31 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 28 23 27 22 26 25 30 25 30 25 31 27 31 28 31 29 31 30 29 31 32 31 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 28 23 27 22 26 25 30 25 30 25 31 27 31 28 31 29 31 30 29 31 32 31 32 30 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 28 23 27 22 26 25 30 25 30 25 31 27 31 28 31 29 31 30 29 31 32 31 32 30 32 27 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 28 23 27 22 26 25 30 25 30 25 31 27 31 28 31 29 31 30 26 30 28 30 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 28 23 27 22 26 25 30 25 30 25 31 27 31 28 31 29 31 30 26 30 26 27 25 29 21 30 23 28 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 28 23 27 22 26 25 30 25 30 25 31 27 31 28 31 29 31 30 26 30 26 27 25 29 21 30 23 28 23 27 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 28 23 27 22 26 25 30 25 30 25 31 27 31 28 31 29 31 30 26 30 26 27 25 29 21 30 23 28 24 27 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 28 23 27 22 26 25 30 25 30 25 31 27 31 28 31 29 31 30 26 30 26 27 25 29 21 30 23 28 25 27 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 28 23 27 22 26 25 30 25 30 25 31 27 31 28 31 29 31 30 26 30 26 27 25 29 21 30 23 28 26 27 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 28 23 27 22 26 25 30 25 30 25 31 27 31 28 31 29 31 30 26 30 26 27 25 29 21 30 23 28 26 27 26 26 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 28 23 27 22 26 25 30 25 30 25 31 27 31 28 31 29 31 30 26 30 26 27 25 29 21 30 23 28 26 27 26 26 24 25 26 25 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 29 23 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 30 23 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 31 23 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 32 23 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 23 18 26 20 26 21 24 22 27 22 26 21 32 23 32 24 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 24 18 plist = 13 24 16 24 16 23 18 21 16 20 23 18 22 17 25 18 plist = 18 3 plist = 18 3 18 4 plist = 18 3 19 4 plist = 18 3 19 4 18 3 19 5 plist = 18 3 19 4 18 3 19 5 19 6 plist = 18 3 19 4 18 3 19 5 20 6 plist = 18 3 19 4 18 3 19 5 20 6 19 5 20 7 plist = 18 3 19 4 18 3 19 5 20 6 19 5 21 6 plist = 18 3 19 4 18 3 19 5 20 6 19 5 21 6 19 5 21 5 plist = 18 3 19 4 18 3 19 5 20 6 19 5 21 6 19 5 22 5 plist = 18 3 19 4 18 3 19 5 20 6 19 5 21 6 19 5 22 5 22 3 plist = 18 3 19 4 18 3 19 5 20 6 19 5 21 6 19 5 22 5 22 3 22 2 plist = 18 3 19 4 18 3 19 5 20 6 19 5 21 6 19 5 22 5 22 3 22 2 20 1 22 1 plist = 18 3 19 4 18 3 19 5 20 6 19 5 21 6 19 5 22 5 22 3 22 2 20 1 23 1 plist = 18 3 19 4 18 3 19 5 20 6 19 5 21 6 19 5 22 5 22 3 22 2 20 1 23 1 23 2 plist = 18 3 19 4 18 3 19 5 20 6 19 5 21 6 19 5 22 5 22 3 22 2 20 1 23 1 23 2 23 3 plist = 18 3 19 4 18 3 19 5 20 6 19 5 21 6 19 5 22 5 22 3 22 2 20 1 23 1 23 2 24 3 plist = 18 3 19 4 18 3 19 5 20 6 19 5 21 6 19 5 22 5 22 3 22 2 20 1 23 1 23 2 24 3 24 2 plist = 18 3 19 4 18 3 19 5 20 6 19 5 21 6 19 5 22 5 22 3 22 2 20 1 23 1 23 2 24 3 24 2 24 1 plist = 18 3 19 4 18 3 19 5 20 6 19 5 21 6 19 5 22 5 22 3 22 2 20 1 23 1 23 2 24 3 24 2 25 1 plist = 19 11 plist = 19 11 17 12 17 14 19 14 plist = 20 16 plist = 21 10 plist = 22 10 plist = 22 10 21 9 22 11 plist = 22 10 21 9 23 11 plist = 22 10 21 9 23 11 22 10 23 12 plist = 22 10 21 9 23 11 22 10 24 12 plist = 22 10 21 9 23 11 22 10 24 12 24 11 plist = 22 10 21 9 23 11 22 10 24 12 24 11 22 10 22 9 24 9 plist = 22 10 21 9 23 11 22 10 24 12 24 11 22 10 22 9 25 9 plist = 22 10 21 9 23 11 22 10 24 12 24 11 22 10 22 9 25 9 25 8 plist = 22 10 21 9 23 11 22 10 24 12 24 11 22 10 22 9 25 9 26 8 plist = 22 10 21 9 23 11 22 10 24 12 24 11 22 10 22 9 25 9 27 7 plist = 22 10 21 9 23 11 22 10 24 12 24 11 22 10 22 9 25 9 28 7 plist = 22 10 21 9 23 11 22 10 24 12 24 11 22 10 22 9 25 9 29 7 plist = 22 10 21 9 23 11 22 10 24 12 24 11 22 10 22 9 25 9 30 6 plist = 22 10 21 9 23 11 22 10 24 12 24 11 22 10 22 9 25 9 31 6 plist = 22 10 21 9 23 11 22 10 24 12 24 11 22 10 22 9 25 9 31 6 31 5 plist = 22 10 21 9 23 11 22 10 24 12 24 11 22 10 22 9 25 9 31 6 31 5 29 4 31 4 plist = 22 10 21 9 23 11 22 10 24 12 24 11 22 10 22 9 25 9 31 6 31 5 29 4 32 4 plist = 22 10 21 9 23 11 22 10 24 12 24 11 22 10 22 9 25 9 31 6 31 5 29 4 32 4 32 6 plist = 22 10 21 9 23 11 22 10 24 12 25 11 plist = 22 10 21 9 23 11 22 10 24 12 25 11 25 12 plist = 22 10 21 9 23 11 22 10 24 12 25 11 26 12 plist = 22 10 21 9 23 11 22 10 24 12 25 11 26 12 26 11 plist = 26 4 plist = 25 16 27 16 plist = 25 16 28 17 plist = 28 19 plist = 29 12 plist = 29 12 29 13 plist = 29 12 29 13 29 14 plist = 29 12 29 13 29 14 29 15 plist = 29 12 29 13 29 14 30 15 plist = 30 20 plist = 31 20 plist = 31 20 31 19 plist = 31 20 32 18 plist = 30 9 plist = 30 9 30 10 plist = 30 9 30 9 31 11 plist = 30 9 30 9 32 11 plist = 32 13 Total number of sites = 1024 Number of occupied sites = 595 Observed probability of occupation = 0.581055 Number of components = 34 Mean component size = 17.3529 component size distribution Size Number Probability 1 14 0.424242 2 6 0.181818 3 1 0.030303 4 1 0.030303 7 1 0.030303 8 1 0.030303 9 1 0.030303 10 1 0.030303 11 2 0.0606061 36 1 0.030303 55 1 0.030303 89 1 0.030303 157 1 0.030303 164 1 0.030303 Number of horizontal spanning components = 0 Number of vertical spanning components = 0 Graphics saved as "occupation.png" k = 0 k = 1 k = 1 k = 1 k = 0 k = 0 k = 1 k = 1 k = 0 k = 1 k = 0 k = 12 k = 12 k = 0 k = 0 k = 18 k = 0 k = 0 k = 20 k = 20 k = 20 k = 20 k = 20 k = 20 k = 20 k = 0 k = 0 k = 0 k = 0 k = 0 k = 23 k = 23 k = 1 k = 1 k = 1 k = 0 k = 0 k = 0 k = 1 k = 0 k = 0 k = 1 k = 1 k = 0 k = 0 k = 1 k = 0 k = 0 k = 0 k = 0 k = 0 k = 0 k = 20 k = 20 k = 20 k = 20 k = 0 k = 0 k = 23 k = 23 k = 23 k = 0 k = 0 k = 23 k = 1 k = 0 k = 0 k = 1 k = 1 k = 0 k = 1 k = 1 k = 1 k = 0 k = 1 k = 1 k = 1 k = 1 k = 1 k = 0 k = 20 k = 20 k = 20 k = 0 k = 20 k = 20 k = 20 k = 20 k = 0 k = 25 k = 0 k = 0 k = 23 k = 23 k = 0 k = 23 k = 1 k = 1 k = 1 k = 1 k = 1 k = 1 k = 1 k = 1 k = 0 k = 1 k = 0 k = 1 k = 1 k = 1 k = 1 k = 0 k = 20 k = 20 k = 20 k = 0 k = 20 k = 0 k = 0 k = 0 k = 25 k = 25 k = 0 k = 0 k = 23 k = 23 k = 23 k = 23 k = 0 k = 0 k = 1 k = 0 k = 0 k = 0 k = 1 k = 0 k = 1 k = 1 k = 0 k = 1 k = 0 k = 1 k = 1 k = 0 k = 0 k = 20 k = 20 k = 20 k = 20 k = 20 k = 0 k = 0 k = 0 k = 25 k = 0 k = 0 k = 0 k = 23 k = 23 k = 0 k = 0 k = 0 k = 0 k = 0 k = 0 k = 0 k = 1 k = 1 k = 1 k = 1 k = 1 k = 1 k = 1 k = 0 k = 1 k = 1 k = 0 k = 20 k = 20 k = 20 k = 20 k = 0 k = 23 k = 0 k = 0 k = 0 k = 23 k = 0 k = 23 k = 23 k = 23 k = 23 k = 7 k = 7 k = 0 k = 0 k = 2 k = 0 k = 0 k = 1 k = 1 k = 1 k = 1 k = 0 k = 0 k = 1 k = 1 k = 1 k = 0 k = 0 k = 20 k = 20 k = 0 k = 0 k = 23 k = 23 k = 0 k = 23 k = 23 k = 23 k = 23 k = 0 k = 0 k = 0 k = 0 k = 0 k = 0 k = 2 k = 2 k = 2 k = 0 k = 0 k = 1 k = 0 k = 1 k = 0 k = 13 k = 0 k = 0 k = 1 k = 0 k = 0 k = 20 k = 0 k = 23 k = 23 k = 0 k = 23 k = 23 k = 23 k = 0 k = 23 k = 0 k = 31 k = 0 k = 0 k = 2 k = 2 k = 2 k = 2 k = 0 k = 0 k = 1 k = 1 k = 1 k = 1 k = 0 k = 0 k = 0 k = 1 k = 1 k = 1 k = 0 k = 20 k = 20 k = 0 k = 23 k = 23 k = 23 k = 23 k = 23 k = 0 k = 0 k = 0 k = 31 k = 31 k = 0 k = 0 k = 0 k = 0 k = 0 k = 2 k = 2 k = 2 k = 0 k = 0 k = 1 k = 0 k = 1 k = 0 k = 1 k = 1 k = 0 k = 1 k = 1 k = 0 k = 0 k = 23 k = 23 k = 23 k = 23 k = 0 k = 0 k = 0 k = 0 k = 0 k = 31 k = 31 k = 0 k = 31 k = 2 k = 2 k = 2 k = 0 k = 2 k = 2 k = 0 k = 1 k = 0 k = 0 k = 1 k = 1 k = 1 k = 0 k = 17 k = 0 k = 0 k = 21 k = 21 k = 0 k = 23 k = 23 k = 23 k = 23 k = 23 k = 23 k = 0 k = 29 k = 0 k = 31 k = 31 k = 31 k = 0 k = 0 k = 2 k = 2 k = 2 k = 2 k = 0 k = 1 k = 1 k = 1 k = 1 k = 0 k = 1 k = 1 k = 0 k = 0 k = 21 k = 21 k = 0 k = 0 k = 0 k = 23 k = 23 k = 23 k = 23 k = 23 k = 0 k = 29 k = 29 k = 0 k = 0 k = 0 k = 0 k = 0 k = 0 k = 0 k = 2 k = 2 k = 2 k = 0 k = 0 k = 1 k = 0 k = 1 k = 1 k = 0 k = 0 k = 19 k = 0 k = 21 k = 0 k = 0 k = 0 k = 23 k = 0 k = 23 k = 0 k = 0 k = 0 k = 29 k = 29 k = 0 k = 32 k = 32 k = 0 k = 0 k = 0 k = 8 k = 0 k = 2 k = 2 k = 2 k = 2 k = 0 k = 0 k = 0 k = 0 k = 0 k = 0 k = 0 k = 21 k = 21 k = 21 k = 0 k = 24 k = 0 k = 16 k = 0 k = 26 k = 0 k = 0 k = 29 k = 29 k = 0 k = 32 k = 0 k = 0 k = 0 k = 0 k = 0 k = 2 k = 2 k = 2 k = 2 k = 0 k = 2 k = 0 k = 0 k = 0 k = 0 k = 14 k = 14 k = 0 k = 0 k = 0 k = 0 k = 24 k = 0 k = 16 k = 0 k = 26 k = 0 k = 0 k = 29 k = 29 k = 29 k = 0 k = 33 k = 3 k = 0 k = 0 k = 0 k = 2 k = 0 k = 2 k = 0 k = 0 k = 2 k = 2 k = 0 k = 14 k = 14 k = 14 k = 14 k = 14 k = 0 k = 22 k = 22 k = 0 k = 0 k = 16 k = 0 k = 26 k = 26 k = 26 k = 0 k = 29 k = 0 k = 0 k = 0 k = 0 k = 2 k = 2 k = 2 k = 0 k = 0 k = 2 k = 0 k = 2 k = 2 k = 0 k = 2 k = 0 k = 14 k = 14 k = 0 k = 14 k = 14 k = 0 k = 0 k = 0 k = 16 k = 16 k = 0 k = 0 k = 0 k = 26 k = 26 k = 0 k = 0 k = 0 k = 30 k = 2 k = 2 k = 2 k = 0 k = 2 k = 2 k = 2 k = 2 k = 2 k = 2 k = 2 k = 2 k = 2 k = 0 k = 0 k = 0 k = 0 k = 0 k = 0 k = 16 k = 16 k = 16 k = 16 k = 16 k = 16 k = 0 k = 0 k = 0 k = 0 k = 0 k = 30 k = 30 k = 0 k = 0 k = 2 k = 2 k = 0 k = 2 k = 2 k = 2 k = 2 k = 0 k = 2 k = 2 k = 0 k = 15 k = 15 k = 0 k = 16 k = 0 k = 0 k = 16 k = 0 k = 16 k = 16 k = 0 k = 0 k = 0 k = 28 k = 28 k = 0 k = 30 k = 30 k = 0 k = 4 k = 0 k = 2 k = 0 k = 2 k = 2 k = 2 k = 2 k = 0 k = 0 k = 0 k = 2 k = 0 k = 0 k = 0 k = 16 k = 16 k = 16 k = 16 k = 16 k = 0 k = 0 k = 16 k = 16 k = 16 k = 16 k = 0 k = 0 k = 30 k = 30 k = 30 k = 0 k = 0 k = 0 k = 2 k = 0 k = 2 k = 0 k = 2 k = 2 k = 0 k = 0 k = 10 k = 0 k = 0 k = 0 k = 16 k = 0 k = 16 k = 16 k = 0 k = 16 k = 16 k = 0 k = 16 k = 0 k = 16 k = 16 k = 0 k = 0 k = 30 k = 0 k = 30 k = 0 k = 2 k = 2 k = 2 k = 2 k = 2 k = 2 k = 2 k = 2 k = 0 k = 0 k = 0 k = 16 k = 16 k = 0 k = 16 k = 16 k = 16 k = 0 k = 16 k = 16 k = 16 k = 16 k = 0 k = 16 k = 16 k = 16 k = 16 k = 0 k = 0 k = 16 k = 0 k = 16 k = 2 k = 2 k = 2 k = 2 k = 2 k = 2 k = 2 k = 0 k = 2 k = 0 k = 11 k = 0 k = 16 k = 0 k = 16 k = 16 k = 0 k = 16 k = 16 k = 0 k = 16 k = 16 k = 0 k = 16 k = 0 k = 16 k = 16 k = 16 k = 16 k = 16 k = 16 k = 16 k = 0 k = 2 k = 2 k = 0 k = 0 k = 2 k = 2 k = 2 k = 2 k = 0 k = 0 k = 0 k = 16 k = 16 k = 16 k = 16 k = 0 k = 16 k = 16 k = 16 k = 16 k = 0 k = 16 k = 16 k = 0 k = 0 k = 16 k = 0 k = 0 k = 0 k = 16 k = 16 k = 0 k = 2 k = 2 k = 2 k = 2 k = 0 k = 2 k = 2 k = 0 k = 0 k = 2 k = 0 k = 16 k = 0 k = 0 k = 16 k = 16 k = 0 k = 16 k = 16 k = 0 k = 0 k = 0 k = 16 k = 16 k = 16 k = 16 k = 16 k = 16 k = 16 k = 0 k = 0 k = 2 k = 2 k = 2 k = 0 k = 0 k = 0 k = 2 k = 2 k = 2 k = 0 k = 2 k = 2 k = 0 k = 2 k = 0 k = 16 k = 16 k = 16 k = 0 k = 16 k = 16 k = 0 k = 16 k = 0 k = 16 k = 16 k = 0 k = 0 k = 16 k = 16 k = 0 k = 16 k = 0 k = 0 k = 2 k = 2 k = 2 k = 2 k = 0 k = 0 k = 2 k = 0 k = 2 k = 2 k = 0 k = 2 k = 0 k = 0 k = 16 k = 0 k = 16 k = 0 k = 0 k = 16 k = 16 k = 16 k = 16 k = 16 k = 0 k = 16 k = 0 k = 16 k = 16 k = 16 k = 0 k = 2 k = 2 k = 2 k = 0 k = 2 k = 2 k = 0 k = 2 k = 2 k = 2 k = 2 k = 2 k = 2 k = 2 k = 0 k = 16 k = 0 k = 16 k = 0 k = 0 k = 16 k = 16 k = 0 k = 0 k = 16 k = 16 k = 16 k = 0 k = 16 k = 16 k = 0 k = 0 k = 2 k = 2 k = 2 k = 2 k = 2 k = 2 k = 2 k = 0 k = 0 k = 2 k = 2 k = 2 k = 0 k = 2 k = 0 k = 16 k = 16 k = 16 k = 16 k = 0 k = 16 k = 0 k = 0 k = 16 k = 16 k = 0 k = 16 k = 0 k = 16 k = 16 k = 0 k = 5 k = 0 k = 2 k = 2 k = 2 k = 0 k = 2 k = 0 k = 0 k = 0 k = 0 k = 2 k = 2 k = 0 k = 2 k = 0 k = 16 k = 16 k = 16 k = 0 k = 16 k = 16 k = 16 k = 16 k = 16 k = 16 k = 16 k = 16 k = 0 k = 16 k = 16 k = 16 k = 0 k = 2 k = 2 k = 0 k = 2 k = 2 k = 0 k = 0 k = 0 k = 2 k = 2 k = 2 k = 2 k = 0 k = 2 k = 2 k = 0 k = 0 k = 16 k = 16 k = 0 k = 0 k = 16 k = 0 k = 0 k = 0 k = 16 k = 0 k = 16 k = 16 k = 16 k = 16 k = 6 k = 0 k = 0 k = 0 k = 0 k = 0 k = 0 k = 9 k = 0 k = 2 k = 2 k = 2 k = 0 k = 2 k = 2 k = 0 k = 0 k = 16 k = 16 k = 16 k = 0 k = 0 k = 16 k = 0 k = 27 k = 0 k = 16 k = 16 k = 16 k = 0 k = 0 k = 0 Graphics saved as "components.png" percolation_simulation_test(): Normal end of execution. 13-May-2025 16:30:40