for 1000 simulations overnight

no v_depoly regulation
nucrates = 0:0.2:1.8;
nucrates = repmat(nucrates,1,100);