Ti = np.arange(0.0,N+dt,dt)
tlist = []
for t in Ti:
	t0 = -hlife*log(1-random())/log(2)
	tlist.append(t0)

time = np.sort(tlist)