get_initial_n function

Calculate initial population abundances