Mathematical modeling of epithelial shape changes through apical constriction. Cell division probability was set as p=1.2×10−5. The degree of apical constriction was set as ka=5 initially and increased by 0.0012 at each time step to ka=293. The cell numbers changed from 40 (t=0) to 272 (t=24×104).