def get_subnetwork (matrix, cell):
	"""
	subset the matrix to only include specified cell and its neightbors 
	"""
	neighbors = []
	from_cell = matrix[cell]
	to_cell = matrix[:, cell]
	neighbors = [idx for idx, val in enumerate(from_cell) if val != 0] + [idx for idx, val in enumerate(to_cell) if val != 0]
	return sorted(list(set(neighbors)) + [cell])