(* ::Package:: *)
% run 1000 simulations of mutation without cell growth. we start with a
% fixed number of cells, each of which has some probability to mutate. This
% corresponds to the situation in which cells plated on selective media are
% in some fashion "induced" to undergo mutations allowing survival.
mut_rate = 1e-5; % probability that a particular cells will mutate
n_cells = 128000;
n_simulations = 1000;
outcomes = zeros(n_simulations,1);
for i=1:n_simulations
outcomes(i) = sum(rand(n_cells,1)