def METHOD_NAME(state):
    """
    Return the set of states reachable from the given state
    by epsilon moves.
    """
    # Cache the result
    result = state.METHOD_NAME
    if result is None:
        result = {}
        state.METHOD_NAME = result
        add_to_epsilon_closure(result, state)
    return result