The science
G(N,p): each of the N(N−1)/2 pairs is linked independently with probability p. Mean degree ⟨c⟩ = p(N−1).
Degrees are Binomial(N−1,p), which tends to Poisson: P(k) = e−c ck/k!.
Giant component: its fraction S obeys S = 1 − e−cS. S = 0 for c ≤ 1; just above threshold S ≈ 2(c−1). At c = 1 the largest cluster scales as N2/3.
Below threshold the mean size of the cluster containing a random node is ⟨s⟩ = 1/(1−c); it diverges at c = 1 (finite N rounds it off).
Full connectivity needs c ≈ ln N: the expected number of isolated nodes is N e−c.
Edges use a fixed random number per pair (edge present if Uij < p), so raising ⟨c⟩ only adds edges — a true monotone percolation process. Node positions come from a Fruchterman–Reingold force layout and carry no meaning. Simplifications: the theory is for N→∞ while the sim has finite N, so deviations near c = 1 are expected.