08-Oct-2025 12:16:52

percolation_simulation_test():
  MATLAB/Octave version 6.4.0
  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.6
plist =

   2   2

plist =

   3   2

plist =

   3   2
   2   1
   3   3

plist =

   3   2
   2   1
   3   3
   3   5

plist =

   3   2
   2   1
   3   3
   4   5

plist =

   3   2
   2   1
   4   3

plist =

   3   2
   2   1
   4   3
   4   2

plist =

   3   2
   2   1
   4   3
   4   1
   5   3

plist =

   3   2
   2   1
   4   3
   4   1
   6   3

plist =

   3   2
   2   1
   4   3
   4   1
   7   3

plist =

   3   2
   2   1
   4   3
   4   1
   7   3
   7   4

plist =

   3   2
   2   1
   4   3
   4   1
   7   3
   7   4
   7   5

plist =

   3   2
   2   1
   4   3
   4   1
   7   3
   7   4
   7   5
   5   6
   7   6

plist =

   3   2
   2   1
   4   3
   4   1
   7   3
   7   4
   7   5
   5   6
   7   6
   5   8
   7   9

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    7   10

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    7   10
    5   12
    7   12

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    7   10
    3   10
    5   10

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    8    9

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
    9    7

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   10    6

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   11    5

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    7

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   13    7

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   14    3

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   15    3

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   16    1

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   17    2

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    3

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   19    4

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   20    4

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   20    4
   20    6

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   20    4
   20    6
   18    7
   20    7

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   20    4
   20    6
   18    7
   20    7
   18    8
   20    8

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   20    4
   20    6
   18    7
   20    7
   18    8
   20    8
   18    9
   20    9

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   20    4
   20    6
   18    7
   20    7
   18    8
   20    8
   18    9
   20    9
   18   10
   20   10

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   20    4
   20    6
   18    7
   20    7
   18    8
   20    8
   18    9
   20    9
   18   10
   21    9

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   20    4
   20    6
   18    7
   20    7
   18    8
   20    8
   18    9
   20    9
   18   10
   21    9
   21    8

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   20    4
   20    6
   18    7
   20    7
   18    8
   20    8
   18    9
   20    9
   18   10
   21    9
   21    8
   21    7

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   20    4
   20    6
   18    7
   20    7
   18    8
   20    8
   18    9
   20    9
   18   10
   21    9
   21    8
   22    7

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   20    4
   20    6
   18    7
   20    7
   18    8
   20    8
   18    9
   20    9
   18   10
   21    9
   21    8
   22    7
   22    8

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   20    4
   20    6
   18    7
   20    7
   18    8
   20    8
   18    9
   20    9
   18   10
   21    9
   21    8
   22    7
   22    8
   22    9

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   20    4
   20    6
   18    7
   20    7
   18    8
   20    8
   18    9
   20    9
   18   10
   21    9
   21    8
   22    7
   22    8
   23    9

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   20    4
   20    6
   18    7
   20    7
   18    8
   20    8
   18    9
   20    9
   18   10
   21    9
   21    8
   22    7
   22    8
   23    9
   22    8
   21   11
   23   11

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   20    4
   20    6
   18    7
   20    7
   18    8
   20    8
   18    9
   20    9
   18   10
   21    9
   21    8
   22    7
   22    8
   23    9
   22    8
   21   11
   24   11

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   20    4
   20    6
   18    7
   20    7
   18    8
   20    8
   18    9
   20    9
   18   10
   21    9
   21    8
   22    7
   22    8
   23    9
   22    8
   21   11
   24   11
   24   12

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   20    4
   20    6
   18    7
   20    7
   18    8
   20    8
   18    9
   20    9
   18   10
   21    9
   21    8
   22    7
   22    8
   23    9
   22    8
   21   11
   24   11
   25   12

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   20    4
   20    6
   18    7
   20    7
   18    8
   20    8
   18    9
   20    9
   18   10
   21    9
   21    8
   22    7
   22    8
   23    9
   23    7

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   20    4
   20    6
   18    7
   20    7
   18    8
   20    8
   18    9
   20    9
   18   10
   21    9
   21    8
   22    7
   22    8
   24    9

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   21    4

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   22    4

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   22    4
   22    2

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   22    4
   23    2

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   23    4

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   24    4

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   25    4

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   11    4
   12    5
   14    5
   16    2
   17    1
   18    2
   17    1
   18    2
   25    4
   25    3

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   10    3
   10    2
   12    2

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   10    3
   10    2
   12    2
    9    1
    9    2
    8    3
    9    2
    7    4
    7    3
    7    2
    9    2

plist =

    3    2
    2    1
    4    3
    4    1
    7    3
    7    4
    7    5
    5    6
    7    6
    5    8
    7    9
    5   10
    8   10
    9    9
    9    8
   11    6
   12    5
   10    3
   10    2
   12    2
    9    1
    9    2
    8    3
    9    2
    7    4
    7    3
    7    2
    9    2
    7    1
    9    1

plist =

    2   12

plist =

    2   12
    2   14

plist =

    1   18
    3   18

plist =

    1   18
    4   18

plist =

    1   18
    5   18

plist =

    1   18
    5   18
    5   16

plist =

    1   18
    5   18
    5   16
    3   15
    5   15

plist =

    1   18
    5   18
    5   16
    3   15
    6   15

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    6   16

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16
    6   15
    7   17

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16
    6   15
    7   16
    8   18

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16
    6   15
    7   16
    8   18
    6   19
    8   19

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16
    6   15
    7   16
    8   18
    6   19
    8   18
    9   21

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16
    6   15
    7   16
    8   18
    6   19
    8   18
    9   21
    7   22
    9   22

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16
    6   15
    7   16
    8   18
    6   19
    8   18
    9   21
    7   22
    9   21
   10   23

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16
    6   15
    7   16
    8   18
    6   19
    8   18
    9   21
    4   22
    4   21
    6   19

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16
    6   15
    7   16
    8   18
    6   19
    8   18
    9   21
    4   22
    3   22
    2   23
    3   22
    4   25

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16
    6   15
    7   16
    8   18
    6   19
    8   18
    9   21
    4   22
    3   22
    2   23
    1   22
    2   21
    2   21

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16
    6   15
    6   15
    8   15

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16
    6   15
    6   15
    9   15

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16
    6   15
    6   15
   10   15

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16
    6   15
    6   15
   11   15

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16
    6   15
    6   15
   12   15

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16
    6   15
    6   15
   13   15

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16
    6   15
    6   15
   13   14
   14   16

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16
    6   15
    6   15
   13   14
   15   16

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16
    6   15
    6   15
   13   14
   16   16

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16
    6   15
    6   15
   13   14
   16   16
   16   15

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16
    6   15
    6   15
   13   14
   16   16
   17   15

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16
    6   15
    6   15
   13   14
   16   16
   17   15
   16   14
   17   17

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16
    6   15
    6   15
   13   14
   16   16
   18   15

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16
    6   15
    6   15
   13   14
   16   16
   19   14

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16
    6   15
    6   15
   14   14

plist =

    1   18
    5   18
    5   16
    3   15
    6   15
    5   14
    7   16
    6   15
    6   15
   14   14
   12   13
   14   13

plist =

    2   27

plist =

    1   30
    3   30

plist =

    1   30
    3   29
    2   32
    4   32

plist =

    1   30
    3   29
    2   32
    5   32

plist =

    1   30
    3   29
    2   32
    6   32

plist =

    1   30
    3   29
    2   32
    6   32
    6   31

plist =

    1   30
    3   29
    2   32
    6   32
    7   31

plist =

    1   30
    4   29

plist =

    1   30
    4   29
    4   28

plist =

    1   30
    4   29
    5   28

plist =

    1   30
    4   29
    5   28
    5   27

plist =

    1   30
    4   29
    5   28
    6   27

plist =

    1   30
    4   29
    5   28
    6   27
    5   26
    6   28

plist =

    1   30
    4   29
    5   28
    6   27
    5   26
    7   28

plist =

    1   30
    4   29
    5   28
    6   27
    5   26
    7   28
    6   27
    7   29

plist =

    1   30
    4   29
    5   28
    6   27
    5   26
    7   28
    6   27
    8   28

plist =

    1   30
    4   29
    5   28
    6   27
    5   26
    7   28
    5   26
    7   26

plist =

    1   30
    4   29
    5   28
    6   27
    5   26
    7   28
    5   26
    7   26
    7   25

plist =

    1   30
    4   29
    5   28
    6   27
    5   26
    7   28
    5   26
    7   26
    8   26

plist =

    9   30

plist =

   10   30

plist =

   10   30
    9   29
   10   31

plist =

   10   30
    9   29
   10   31
    8   32
   10   32

plist =

   10   30
    9   29
   10   31
    8   32
   11   32

plist =

   10   30
    9   29
   10   31
    8   32
   11   32
   11   31

plist =

   10   30
    9   29
   10   31
    8   32
   11   32
   11   31
   11   30

plist =

   10   30
    9   29
   10   31
    8   32
   11   32
   11   31
   11   30
    9   29
   11   29

plist =

   10   30
    9   29
   10   31
    8   32
   11   32
   11   31
   11   30
    9   29
   12   29

plist =

   10   30
    9   29
   10   31
    8   32
   11   32
   11   31
   11   30
    9   29
   12   29
   12   30

plist =

   10   30
    9   29
   10   31
    8   32
   11   32
   11   31
   11   30
    9   29
   12   29
   12   30
   12   31

plist =

   10   30
    9   29
   10   31
    8   32
   11   32
   11   31
   11   30
    9   29
   12   29
   12   30
   13   31

plist =

   10   30
    9   29
   10   31
    8   32
   11   32
   11   31
   11   30
    9   29
   12   29
   12   30
   13   31
   13   30

plist =

   10   17

plist =

   11   19

plist =

   11   18
   12   20

plist =

   11   18
   12   20
   12   21

plist =

   11   18
   12   20
   13   21

plist =

   11   18
   12   20
   13   21
   13   20

plist =

   11   18
   12   20
   13   21
   14   20

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   14   21

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   15   20

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   16   20

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21
   17   23

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21
   18   23

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21
   18   23
   18   24

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21
   18   23
   18   24
   16   25
   18   25

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21
   18   23
   18   24
   16   25
   19   23

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21
   18   23
   18   24
   16   25
   19   23
   19   22

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21
   18   23
   18   24
   16   25
   19   23
   20   22

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21
   18   23
   18   24
   16   25
   19   23
   20   22
   19   21
   20   23

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21
   18   23
   18   24
   16   25
   19   23
   20   22
   19   21
   21   23

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21
   18   23
   18   24
   16   25
   19   23
   20   22
   19   21
   21   23
   20   22
   21   24

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21
   18   23
   18   24
   16   25
   19   23
   20   22
   19   21
   21   23
   20   22
   21   24
   21   25

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21
   18   23
   18   24
   16   25
   19   23
   20   22
   19   21
   21   23
   20   22
   21   24
   21   25
   21   26

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21
   18   23
   18   24
   16   25
   19   23
   20   22
   19   21
   21   23
   20   22
   21   24
   21   25
   22   26

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21
   18   23
   18   24
   16   25
   19   23
   20   22
   19   21
   21   23
   20   22
   21   24
   21   25
   22   26
   21   25
   22   27

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21
   18   23
   18   24
   16   25
   19   23
   20   22
   19   21
   21   23
   20   22
   21   24
   21   25
   22   26
   22   24

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21
   18   23
   18   24
   16   25
   19   23
   20   22
   19   21
   21   23
   20   22
   21   24
   21   25
   22   26
   22   24
   22   23

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21
   18   23
   18   24
   16   25
   19   23
   20   22
   19   21
   21   23
   20   22
   21   24
   21   25
   22   26
   22   24
   23   24

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21
   18   23
   18   24
   16   25
   19   23
   20   22
   19   21
   21   23
   20   22
   21   24
   21   25
   22   26
   22   24
   24   24

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21
   18   23
   18   24
   16   25
   19   23
   20   22
   19   21
   21   23
   20   22
   21   24
   21   25
   22   26
   22   24
   25   27

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21
   18   23
   18   24
   16   25
   19   23
   20   22
   18   20
   20   20

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21
   18   23
   18   24
   16   25
   19   23
   20   22
   18   20
   21   20

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21
   18   23
   18   24
   16   25
   19   23
   20   22
   18   20
   22   21

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21
   18   23
   18   24
   15   25
   17   27

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21
   18   23
   18   24
   15   25
   17   27
   15   28
   17   28

plist =

   11   18
   12   20
   13   21
   14   20
   13   19
   17   20
   17   21
   18   23
   18   24
   15   25
   17   27
   15   28
   18   27

plist =

   13   10

plist =

   14    9

plist =

   17   32

plist =

   20   29

plist =

   21   29

plist =

   22   31

plist =

   23   31

plist =

   24   14

plist =

   24   14
   24   15

plist =

   24   14
   24   14
   25   17

plist =

   24   14
   24   14
   26   17

plist =

   24   14
   24   14
   27   17

plist =

   24   14
   24   14
   27   17
   27   16

plist =

   24   14
   24   14
   27   17
   27   15
   28   17

plist =

   24   14
   24   14
   27   17
   27   15
   29   17

plist =

   24   14
   24   14
   27   17
   27   15
   29   17
   29   18

plist =

   24   14
   24   14
   27   17
   27   15
   29   17
   29   18
   27   19
   29   19

plist =

   24   14
   24   14
   27   17
   27   15
   29   17
   29   18
   27   19
   29   19
   27   20
   29   20

plist =

   24   14
   24   14
   27   17
   27   15
   29   17
   29   18
   27   19
   29   19
   27   20
   29   20
   29   21

plist =

   24   14
   24   14
   27   17
   27   15
   29   17
   29   18
   27   19
   29   19
   27   20
   29   20
   30   21

plist =

   24   14
   24   14
   27   17
   27   15
   29   17
   29   18
   27   19
   29   19
   27   20
   29   20
   30   21
   29   20
   30   22

plist =

   24   14
   24   14
   27   17
   27   15
   29   17
   29   18
   27   19
   29   19
   27   20
   29   20
   30   21
   29   20
   30   22
   28   23
   27   27
   27   30
   29   31

plist =

   24   14
   24   14
   27   17
   27   15
   29   17
   29   18
   27   19
   29   19
   27   20
   29   20
   30   21
   29   20
   30   22
   28   23
   27   27
   27   30
   29   31
   27   32
   29   32

plist =

   24   14
   24   14
   27   17
   27   15
   29   17
   29   18
   27   19
   29   19
   27   20
   29   20
   30   21
   29   20
   31   22

plist =

   24   14
   24   14
   27   17
   27   15
   29   17
   29   18
   27   19
   29   19
   27   20
   29   20
   30   21
   29   20
   31   22
   29   20
   31   20

plist =

   24   14
   24   14
   27   17
   27   15
   29   17
   29   18
   27   19
   29   19
   27   20
   29   20
   30   21
   29   20
   31   22
   29   20
   31   20
   29   19
   31   19

plist =

   24   14
   24   14
   27   17
   27   15
   29   17
   29   18
   27   19
   29   19
   27   20
   29   20
   30   21
   29   20
   31   22
   29   20
   31   20
   29   19
   31   18
   32   20

plist =

   24   14
   24   14
   27   17
   27   15
   29   17
   29   18
   27   19
   29   19
   27   20
   29   20
   30   21
   29   20
   31   22
   29   20
   31   20
   29   19
   32   18

plist =

   24   14
   24   14
   27   17
   27   15
   29   17
   29   18
   27   19
   29   19
   27   20
   29   20
   30   21
   29   20
   31   22
   29   20
   31   20
   29   19
   32   18
   30   17
   30   16
   32   16

plist =

   24   14
   24   14
   27   17
   27   15
   29   17
   29   18
   27   19
   29   19
   27   20
   29   20
   30   21
   29   20
   32   23

plist =

   24   14
   24   14
   27   17
   27   15
   29   17
   29   18
   27   19
   29   19
   27   20
   29   20
   30   21
   29   20
   32   23
   32   24

plist =

   24   14
   24   14
   27   17
   27   15
   29   17
   29   18
   27   19
   29   19
   27   20
   29   20
   30   21
   29   20
   32   23
   32   24
   32   25

plist =

   24   14
   24   14
   27   17
   28   15

plist =

   24   14
   24   14
   27   17
   29   14

plist =

   24   14
   24   14
   27   17
   29   14
   27   13
   29   13

plist =

   27    1

plist =

   28    1

plist =

   29    3

plist =

   30    3

plist =

   30    3
   30    2

plist =

   30    3
   31    4

plist =

   30    3
   32    5

plist =

   31   28

plist =

   31   28
   31   29


  Total number of sites = 1024
  Number of occupied sites = 584
  Observed probability of occupation = 0.570312
  Number of components = 46
  Mean component size = 12.587

  component size distribution
    Size  Number  Probability

        1  24  0.533333
        2   4  0.0888889
        3   2  0.0444444
        4   3  0.0666667
        5   3  0.0666667
        6   1  0.0222222
        7   1  0.0222222
       20   1  0.0222222
       24   1  0.0222222
       33   1  0.0222222
       78   1  0.0222222
       88   1  0.0222222
       94   1  0.0222222
       164   1  0.0222222

  Number of horizontal spanning components = 0
  Number of vertical spanning components   = 0
  Graphics saved as "occupation.png"
k = 0
k = 1
k = 0
k = 1
k = 0
k = 1
k = 1
k = 1
k = 1
k = 1
k = 1
k = 0
k = 18
k = 0
k = 1
k = 1
k = 1
k = 1
k = 1
k = 1
k = 1
k = 0
k = 0
k = 0
k = 0
k = 34
k = 34
k = 34
k = 0
k = 0
k = 43
k = 0
k = 1
k = 1
k = 1
k = 1
k = 0
k = 0
k = 1
k = 1
k = 1
k = 1
k = 1
k = 1
k = 0
k = 0
k = 1
k = 1
k = 1
k = 1
k = 1
k = 0
k = 1
k = 1
k = 1
k = 0
k = 1
k = 0
k = 0
k = 0
k = 35
k = 35
k = 0
k = 44
k = 0
k = 1
k = 1
k = 1
k = 1
k = 1
k = 1
k = 1
k = 1
k = 1
k = 1
k = 0
k = 1
k = 1
k = 1
k = 0
k = 1
k = 1
k = 0
k = 0
k = 1
k = 0
k = 0
k = 1
k = 1
k = 0
k = 0
k = 35
k = 35
k = 35
k = 0
k = 0
k = 0
k = 1
k = 0
k = 0
k = 0
k = 1
k = 1
k = 1
k = 1
k = 0
k = 1
k = 0
k = 1
k = 0
k = 1
k = 0
k = 0
k = 1
k = 1
k = 1
k = 1
k = 1
k = 1
k = 1
k = 1
k = 0
k = 38
k = 0
k = 0
k = 35
k = 35
k = 0
k = 0
k = 1
k = 1
k = 1
k = 0
k = 1
k = 1
k = 0
k = 0
k = 1
k = 1
k = 1
k = 1
k = 1
k = 0
k = 0
k = 21
k = 0
k = 1
k = 0
k = 0
k = 0
k = 0
k = 1
k = 0
k = 0
k = 0
k = 41
k = 0
k = 0
k = 35
k = 35
k = 0
k = 0
k = 0
k = 0
k = 1
k = 1
k = 1
k = 0
k = 1
k = 1
k = 1
k = 0
k = 1
k = 0
k = 0
k = 19
k = 0
k = 0
k = 1
k = 1
k = 0
k = 0
k = 1
k = 0
k = 0
k = 0
k = 0
k = 41
k = 0
k = 0
k = 0
k = 35
k = 2
k = 0
k = 0
k = 0
k = 0
k = 1
k = 0
k = 1
k = 1
k = 0
k = 1
k = 1
k = 1
k = 0
k = 0
k = 0
k = 1
k = 1
k = 1
k = 1
k = 1
k = 1
k = 1
k = 0
k = 0
k = 0
k = 0
k = 0
k = 35
k = 35
k = 35
k = 35
k = 2
k = 2
k = 2
k = 0
k = 1
k = 1
k = 0
k = 1
k = 1
k = 0
k = 0
k = 0
k = 0
k = 17
k = 0
k = 0
k = 1
k = 1
k = 1
k = 1
k = 1
k = 1
k = 0
k = 1
k = 0
k = 35
k = 35
k = 35
k = 35
k = 0
k = 35
k = 35
k = 0
k = 0
k = 0
k = 0
k = 0
k = 1
k = 1
k = 1
k = 1
k = 0
k = 16
k = 0
k = 17
k = 17
k = 0
k = 0
k = 1
k = 1
k = 1
k = 1
k = 1
k = 1
k = 1
k = 1
k = 0
k = 35
k = 35
k = 0
k = 0
k = 0
k = 0
k = 0
k = 3
k = 0
k = 1
k = 1
k = 1
k = 1
k = 1
k = 1
k = 0
k = 15
k = 0
k = 17
k = 17
k = 0
k = 0
k = 0
k = 1
k = 1
k = 1
k = 1
k = 0
k = 1
k = 0
k = 0
k = 0
k = 35
k = 0
k = 0
k = 42
k = 0
k = 0
k = 31
k = 3
k = 0
k = 0
k = 1
k = 0
k = 1
k = 0
k = 1
k = 0
k = 0
k = 0
k = 0
k = 17
k = 0
k = 1
k = 1
k = 1
k = 0
k = 0
k = 0
k = 1
k = 1
k = 1
k = 1
k = 0
k = 35
k = 0
k = 0
k = 0
k = 0
k = 0
k = 31
k = 3
k = 3
k = 0
k = 1
k = 1
k = 1
k = 1
k = 0
k = 12
k = 0
k = 4
k = 4
k = 0
k = 0
k = 1
k = 0
k = 1
k = 0
k = 23
k = 0
k = 1
k = 0
k = 1
k = 1
k = 1
k = 0
k = 0
k = 31
k = 0
k = 31
k = 31
k = 31
k = 3
k = 0
k = 8
k = 0
k = 0
k = 0
k = 0
k = 0
k = 0
k = 0
k = 4
k = 4
k = 4
k = 4
k = 0
k = 0
k = 1
k = 0
k = 0
k = 1
k = 1
k = 0
k = 0
k = 0
k = 1
k = 0
k = 31
k = 31
k = 31
k = 0
k = 0
k = 31
k = 3
k = 3
k = 0
k = 0
k = 4
k = 4
k = 0
k = 4
k = 0
k = 0
k = 0
k = 0
k = 4
k = 4
k = 0
k = 4
k = 0
k = 4
k = 4
k = 0
k = 1
k = 0
k = 31
k = 31
k = 0
k = 0
k = 0
k = 31
k = 31
k = 0
k = 31
k = 31
k = 0
k = 0
k = 4
k = 4
k = 4
k = 4
k = 4
k = 4
k = 4
k = 4
k = 4
k = 4
k = 4
k = 0
k = 4
k = 4
k = 4
k = 4
k = 0
k = 1
k = 1
k = 0
k = 31
k = 31
k = 0
k = 0
k = 31
k = 31
k = 0
k = 0
k = 0
k = 31
k = 0
k = 0
k = 0
k = 4
k = 4
k = 4
k = 4
k = 0
k = 0
k = 0
k = 0
k = 0
k = 4
k = 4
k = 4
k = 4
k = 0
k = 0
k = 0
k = 0
k = 0
k = 29
k = 0
k = 31
k = 0
k = 31
k = 31
k = 0
k = 0
k = 31
k = 31
k = 31
k = 0
k = 0
k = 0
k = 4
k = 0
k = 4
k = 4
k = 0
k = 13
k = 13
k = 0
k = 0
k = 0
k = 4
k = 0
k = 4
k = 4
k = 0
k = 0
k = 26
k = 0
k = 0
k = 0
k = 31
k = 31
k = 31
k = 31
k = 31
k = 31
k = 31
k = 31
k = 0
k = 4
k = 4
k = 4
k = 4
k = 4
k = 0
k = 4
k = 4
k = 0
k = 0
k = 14
k = 0
k = 0
k = 0
k = 0
k = 4
k = 0
k = 14
k = 14
k = 0
k = 28
k = 0
k = 0
k = 0
k = 31
k = 0
k = 0
k = 31
k = 31
k = 0
k = 31
k = 31
k = 4
k = 0
k = 0
k = 0
k = 4
k = 4
k = 4
k = 4
k = 0
k = 14
k = 14
k = 0
k = 14
k = 0
k = 0
k = 0
k = 0
k = 14
k = 14
k = 0
k = 0
k = 30
k = 0
k = 0
k = 0
k = 31
k = 31
k = 31
k = 31
k = 31
k = 31
k = 0
k = 4
k = 0
k = 0
k = 0
k = 4
k = 0
k = 0
k = 4
k = 0
k = 0
k = 14
k = 14
k = 14
k = 14
k = 14
k = 14
k = 14
k = 14
k = 14
k = 14
k = 14
k = 0
k = 32
k = 32
k = 0
k = 31
k = 31
k = 31
k = 31
k = 31
k = 31
k = 31
k = 4
k = 4
k = 0
k = 4
k = 4
k = 0
k = 0
k = 4
k = 4
k = 0
k = 14
k = 14
k = 14
k = 14
k = 0
k = 14
k = 14
k = 0
k = 14
k = 0
k = 14
k = 14
k = 0
k = 0
k = 31
k = 31
k = 0
k = 31
k = 31
k = 31
k = 0
k = 0
k = 4
k = 4
k = 4
k = 4
k = 4
k = 4
k = 4
k = 4
k = 4
k = 0
k = 0
k = 0
k = 0
k = 14
k = 0
k = 14
k = 0
k = 14
k = 14
k = 14
k = 0
k = 0
k = 0
k = 0
k = 31
k = 0
k = 39
k = 0
k = 31
k = 31
k = 31
k = 0
k = 4
k = 4
k = 4
k = 4
k = 0
k = 0
k = 4
k = 0
k = 4
k = 4
k = 0
k = 0
k = 0
k = 14
k = 0
k = 14
k = 14
k = 14
k = 14
k = 14
k = 14
k = 14
k = 0
k = 0
k = 0
k = 0
k = 0
k = 31
k = 31
k = 0
k = 31
k = 31
k = 0
k = 0
k = 4
k = 0
k = 9
k = 0
k = 0
k = 10
k = 0
k = 4
k = 0
k = 0
k = 0
k = 0
k = 0
k = 0
k = 14
k = 14
k = 0
k = 14
k = 14
k = 14
k = 14
k = 14
k = 0
k = 0
k = 0
k = 0
k = 31
k = 0
k = 31
k = 31
k = 5
k = 0
k = 4
k = 4
k = 0
k = 7
k = 7
k = 0
k = 4
k = 4
k = 0
k = 0
k = 0
k = 14
k = 14
k = 14
k = 14
k = 14
k = 0
k = 14
k = 14
k = 0
k = 0
k = 14
k = 0
k = 40
k = 40
k = 0
k = 31
k = 0
k = 31
k = 31
k = 0
k = 6
k = 0
k = 0
k = 7
k = 7
k = 7
k = 7
k = 0
k = 4
k = 0
k = 0
k = 0
k = 14
k = 0
k = 14
k = 0
k = 0
k = 0
k = 14
k = 14
k = 14
k = 0
k = 14
k = 0
k = 40
k = 0
k = 0
k = 31
k = 0
k = 0
k = 31
k = 6
k = 6
k = 0
k = 7
k = 7
k = 7
k = 0
k = 0
k = 0
k = 0
k = 0
k = 0
k = 14
k = 14
k = 0
k = 14
k = 14
k = 14
k = 0
k = 0
k = 14
k = 14
k = 0
k = 14
k = 14
k = 0
k = 31
k = 31
k = 31
k = 0
k = 0
k = 0
k = 6
k = 0
k = 7
k = 7
k = 7
k = 7
k = 7
k = 7
k = 0
k = 0
k = 0
k = 0
k = 14
k = 0
k = 14
k = 14
k = 14
k = 0
k = 24
k = 0
k = 0
k = 0
k = 33
k = 0
k = 0
k = 31
k = 0
k = 31
k = 0
k = 36
k = 36
k = 0
k = 0
k = 0
k = 7
k = 7
k = 0
k = 7
k = 7
k = 0
k = 11
k = 11
k = 11
k = 11
k = 0
k = 14
k = 14
k = 0
k = 14
k = 0
k = 24
k = 24
k = 24
k = 0
k = 0
k = 0
k = 31
k = 31
k = 0
k = 31
k = 0
k = 36
k = 36
k = 0
k = 7
k = 7
k = 7
k = 0
k = 0
k = 0
k = 0
k = 11
k = 11
k = 11
k = 11
k = 11
k = 11
k = 0
k = 0
k = 20
k = 0
k = 22
k = 0
k = 24
k = 0
k = 25
k = 0
k = 31
k = 31
k = 31
k = 31
k = 31
k = 0
k = 36
k = 0
k = 45
k = 7
k = 0
k = 7
k = 0
k = 7
k = 7
k = 7
k = 0
k = 11
k = 11
k = 11
k = 11
k = 11
k = 0
k = 0
k = 20
k = 0
k = 0
k = 0
k = 0
k = 25
k = 25
k = 25
k = 0
k = 0
k = 0
k = 0
k = 31
k = 31
k = 0
k = 0
k = 45
k = 7
k = 7
k = 7
k = 7
k = 7
k = 7
k = 0
k = 11
k = 11
k = 11
k = 11
k = 0
k = 11
k = 0
k = 0
k = 20
k = 20
k = 0
k = 0
k = 27
k = 0
k = 25
k = 0
k = 0
k = 31
k = 31
k = 31
k = 31
k = 31
k = 0
k = 37
k = 0
  Graphics saved as "components.png"

percolation_simulation_test()
  Normal end of execution.

08-Oct-2025 12:17:57