def expand_region_with_neighbors(region, neighbors):
    """Expands a given region with its neighboring nodes according to the graph.
    
    Args:
         region: The set of points currently in the region.
         neighbors: The set of border points that are within distance eps from at least one point currently in the region.
                  
    Returns:
         The expanded region.
    """
    
    region_ext = set(region.copy())
    # update region
    for n in neighbors:
        region_ext.add(n)
    
    return region_ext