“…No-flux boundary conditions were used at both outer ends of the domain: We generated large numbers of variable morphogen gradients by numerically solving Eq. 9 with kinetic parameters p, d and D , and cell areas A independently drawn from log-normal distributions for each cell in the domain, as described before [2, 19]. The diffusion coefficient, D , and the degradation rate, d , set the steady-state patterning length scale, , and we report positional quantities relative to the average λ or the average cell diameter, which in turn is chosen to be a fixed multiple of the average λ .…”