def directed_square():
    k = [K(x) for x in 'abcd']
    for i in range(-1,len(k)-1):
        k[i].x = k[i+1]
    seed_node,nsyms = k[0], len('abcd')