code: def refractory(times, refract=0.002):
  	
  	times_refract = []
  	times_refract.append(times[0])
  	for i in range(1,len(times)):
  		if times_refract[-1]+refract <= times[i]:
  			times_refract.append(times[i])        
  	return times_refract
docs: Removes spikes in times list that do not satisfy refractor period
  
  	:param times: list(float) of spike times in seconds
  	:type times: list(float)
  	:param refract: Refractory period in seconds
  	:type refract: float
  	:returns: list(float) of spike times in seconds
  
  	For every interspike interval < refract, 
  	removes the second spike time in list and returns the result