{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<!-- -*- coding: utf-8 -*- -->\n",
    "## Tutorial\n",
    "\n",
    "Start here to begin working with NetworkX.\n",
    "\n",
    "### Creating a graph\n",
    "\n",
    "Create an empty graph with no nodes and no edges."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import networkx as nx\n",
    "G = nx.Graph()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "By definition, a `Graph` is a collection of nodes (vertices) along with\n",
    "identified pairs of nodes (called edges, links, etc).  In NetworkX, nodes can\n",
    "be any hashable object e.g., a text string, an image, an XML object, another\n",
    "Graph, a customized node object, etc.\n",
    "\n",
    "# Nodes\n",
    "\n",
    "The graph `G` can be grown in several ways.  NetworkX includes many graph\n",
    "generator functions and facilities to read and write graphs in many formats.\n",
    "To get started though we'll look at simple manipulations.  You can add one node\n",
    "at a time,"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "G.add_node(1)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "add a list of nodes,"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "G.add_nodes_from([2, 3])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "or add any nbunch of nodes.  An *nbunch* is any iterable container of\n",
    "nodes that is not itself a node in the graph (e.g., a `list`, `set`, `graph`,\n",
    "`file`, etc.)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "H = nx.path_graph(10)\n",
    "G.add_nodes_from(H)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Note that `G` now contains the nodes of `H` as nodes of `G`.\n",
    "In contrast, you could use the graph `H` as a node in `G`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "G.add_node(H)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The graph `G` now contains `H` as a node.  This flexibility is very powerful as\n",
    "it allows graphs of graphs, graphs of files, graphs of functions and much more.\n",
    "It is worth thinking about how to structure your application so that the nodes\n",
    "are useful entities.  Of course you can always use a unique identifier in `G`\n",
    "and have a separate dictionary keyed by identifier to the node information if\n",
    "you prefer.\n",
    "\n",
    "# Edges\n",
    "\n",
    "`G` can also be grown by adding one edge at a time,"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "G.add_edge(1, 2)\n",
    "e = (2, 3)\n",
    "G.add_edge(*e)  # unpack edge tuple*"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "by adding a list of edges,"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "G.add_edges_from([(1, 2), (1, 3)])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "or by adding any ebunch of edges.  An *ebunch* is any iterable\n",
    "container of edge-tuples.  An edge-tuple can be a 2-tuple of nodes or a 3-tuple\n",
    "with 2 nodes followed by an edge attribute dictionary, e.g.,\n",
    "`(2, 3, {'weight': 3.1415})`.  Edge attributes are discussed further below"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "G.add_edges_from(H.edges())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "One can demolish the graph in a similar fashion; using\n",
    "`Graph.remove_node()`,\n",
    "`Graph.remove_nodes_from()`,\n",
    "`Graph.remove_edge()`\n",
    "and\n",
    "`Graph.remove_edges_from()`, e.g."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "G.remove_node(H)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "There are no complaints when adding existing nodes or edges. For example,\n",
    "after removing all nodes and edges,"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "G.clear()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "we add new nodes/edges and NetworkX quietly ignores any that are\n",
    "already present."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "G.add_edges_from([(1, 2), (1, 3)])\n",
    "G.add_node(1)\n",
    "G.add_edge(1, 2)\n",
    "G.add_node(\"spam\")        # adds node \"spam\"\n",
    "G.add_nodes_from(\"spam\")  # adds 4 nodes: 's', 'p', 'a', 'm'"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "At this stage the graph `G` consists of 8 nodes and 2 edges, as can be seen by:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "2"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "G.number_of_nodes()\n",
    "G.number_of_edges()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can examine the nodes and edges. The methods return iterators of nodes,\n",
    "edges, neighbors, etc. This is typically more memory efficient, but it does\n",
    "mean we need to specify what type of container to put the objects in. Here\n",
    "we use lists, though sets, dicts, tuples and other containers may be better\n",
    "in other contexts."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[2, 3]"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "list(G.nodes())\n",
    "list(G.edges())\n",
    "list(G.neighbors(1))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Removing nodes or edges has similar syntax to adding:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [],
   "source": [
    "G.remove_nodes_from(\"spam\")\n",
    "list(G.nodes())\n",
    "G.remove_edge(1, 3)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "When creating a graph structure by instantiating one of the graph\n",
    "classes you can specify data in several formats."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "H = nx.DiGraph(G)   # create a DiGraph using the connections from G\n",
    "list(H.edges())\n",
    "edgelist = [(0, 1), (1, 2), (2, 3)]\n",
    "H = nx.Graph(edgelist)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# What to use as nodes and edges\n",
    "\n",
    "You might notice that nodes and edges are not specified as NetworkX\n",
    "objects.  This leaves you free to use meaningful items as nodes and\n",
    "edges. The most common choices are numbers or strings, but a node can\n",
    "be any hashable object (except `None`), and an edge can be associated\n",
    "with any object `x` using `G.add_edge(n1, n2, object=x)`.\n",
    "\n",
    "As an example, `n1` and `n2` could be protein objects from the RCSB Protein\n",
    "Data Bank, and `x` could refer to an XML record of publications detailing\n",
    "experimental observations of their interaction.\n",
    "\n",
    "We have found this power quite useful, but its abuse\n",
    "can lead to unexpected surprises unless one is familiar with Python.\n",
    "If in doubt, consider using `convert_node_labels_to_integers()` to obtain\n",
    "a more traditional graph with integer labels.\n",
    "\n",
    "# Accessing edges\n",
    "\n",
    "In addition to the methods\n",
    "`Graph.nodes()`,\n",
    "`Graph.edges()`, and\n",
    "`Graph.neighbors()`,\n",
    "fast direct access to the graph data structure is also possible\n",
    "using subscript notation."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{}"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "G[1]  # Warning: do not change the resulting dict\n",
    "G[1][2]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "You can safely set the attributes of an edge using subscript notation\n",
    "if the edge already exists."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [],
   "source": [
    "G.add_edge(1, 3)\n",
    "G[1][3]['color'] = \"blue\""
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Fast examination of all edges is achieved using adjacency iterators.\n",
    "Note that for undirected graphs this actually looks at each edge twice."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(1, 2, 0.125)\n",
      "(2, 1, 0.125)\n",
      "(3, 4, 0.375)\n",
      "(4, 3, 0.375)\n"
     ]
    }
   ],
   "source": [
    "FG = nx.Graph()\n",
    "FG.add_weighted_edges_from([(1, 2, 0.125), (1, 3, 0.75), (2, 4, 1.2), (3, 4, 0.375)])\n",
    "for n, nbrs in FG.adjacency():\n",
    "   for nbr, eattr in nbrs.items():\n",
    "       data = eattr['weight']\n",
    "       if data < 0.5: print('(%d, %d, %.3f)' % (n, nbr, data))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Convenient access to all edges is achieved with the edges method."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(1, 2, 0.125)\n",
      "(3, 4, 0.375)\n"
     ]
    }
   ],
   "source": [
    "for (u, v, d) in FG.edges(data='weight'):\n",
    "    if d < 0.5: print('(%d, %d, %.3f)' % (u, v, d))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Adding attributes to graphs, nodes, and edges\n",
    "\n",
    "Attributes such as weights, labels, colors, or whatever Python object you like,\n",
    "can be attached to graphs, nodes, or edges.\n",
    "\n",
    "Each graph, node, and edge can hold key/value attribute pairs in an associated\n",
    "attribute dictionary (the keys must be hashable).  By default these are empty,\n",
    "but attributes can be added or changed using `add_edge`, `add_node` or direct\n",
    "manipulation of the attribute dictionaries named `G.graph`, `G.node`, and\n",
    "`G.edge` for a graph `G`.\n",
    "\n",
    "## Graph attributes\n",
    "\n",
    "Assign graph attributes when creating a new graph"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'day': 'Friday'}"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "G = nx.Graph(day=\"Friday\")\n",
    "G.graph"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Or you can modify attributes later"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'day': 'Monday'}"
      ]
     },
     "execution_count": 21,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "G.graph['day'] = \"Monday\"\n",
    "G.graph"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Node attributes\n",
    "\n",
    "Add node attributes using `add_node()`, `add_nodes_from()`, or `G.node`"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[(1, {'room': 714, 'time': '5pm'}), (3, {'time': '2pm'})]"
      ]
     },
     "execution_count": 22,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "G.add_node(1, time='5pm')\n",
    "G.add_nodes_from([3], time='2pm')\n",
    "G.nodes[1]\n",
    "G.nodes[1]['room'] = 714\n",
    "list(G.nodes(data=True))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Note that adding a node to `G.node` does not add it to the graph, use\n",
    "`G.add_node()` to add new nodes.\n",
    "\n",
    "# Edge Attributes\n",
    "\n",
    "Add edge attributes using `add_edge()`, `add_edges_from()`, or subscript notation."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [],
   "source": [
    "G.add_edge(1, 2, weight=4.7 )\n",
    "G.add_edges_from([(3, 4), (4, 5)], color='red')\n",
    "G.add_edges_from([(1, 2, {'color': 'blue'}), (2, 3, {'weight': 8})])\n",
    "G[1][2]['weight'] = 4.7"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The special attribute `weight` should be numeric and holds values used by\n",
    "algorithms requiring weighted edges.\n",
    "\n",
    " Directed graphs\n",
    "\n",
    "The `DiGraph` class provides additional methods specific to directed\n",
    "edges, e.g.,\n",
    "`DiGraph.out_edges()`,\n",
    "`DiGraph.in_degree()`,\n",
    "`DiGraph.predecessors()`,\n",
    "`DiGraph.successors()` etc.\n",
    "To allow algorithms to work with both classes easily, the directed versions of\n",
    "`neighbors()` and `degree()` are equivalent to `successors()` and the sum of\n",
    "`in_degree()` and `out_degree()` respectively even though that may feel\n",
    "inconsistent at times."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[2]"
      ]
     },
     "execution_count": 24,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "DG = nx.DiGraph()\n",
    "DG.add_weighted_edges_from([(1, 2, 0.5), (3, 1, 0.75)])\n",
    "DG.out_degree(1, weight='weight')\n",
    "DG.degree(1, weight='weight')\n",
    "list(DG.successors(1))\n",
    "list(DG.neighbors(1))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Some algorithms work only for directed graphs and others are not well\n",
    "defined for directed graphs.  Indeed the tendency to lump directed\n",
    "and undirected graphs together is dangerous.  If you want to treat\n",
    "a directed graph as undirected for some measurement you should probably\n",
    "convert it using `Graph.to_undirected()` or with"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [],
   "source": [
    "H = nx.Graph(G)  # convert G to undirected graph"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Multigraphs\n",
    "\n",
    "NetworkX provides classes for graphs which allow multiple edges\n",
    "between any pair of nodes.  The `MultiGraph` and\n",
    "`MultiDiGraph`\n",
    "classes allow you to add the same edge twice, possibly with different\n",
    "edge data.  This can be powerful for some applications, but many\n",
    "algorithms are not well defined on such graphs.\n",
    "Where results are well defined,\n",
    "e.g., `MultiGraph.degree()` we provide the function.  Otherwise you\n",
    "should convert to a standard graph in a way that makes the measurement\n",
    "well defined."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[1, 2, 3]"
      ]
     },
     "execution_count": 26,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "MG = nx.MultiGraph()\n",
    "MG.add_weighted_edges_from([(1, 2, 0.5), (1, 2, 0.75), (2, 3, 0.5)])\n",
    "dict(MG.degree(weight='weight'))\n",
    "GG = nx.Graph()\n",
    "for n, nbrs in MG.adjacency():\n",
    "   for nbr, edict in nbrs.items():\n",
    "       minvalue = min([d['weight'] for d in edict.values()])\n",
    "       GG.add_edge(n, nbr, weight = minvalue)\n",
    "\n",
    "nx.shortest_path(GG, 1, 3)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Graph generators and graph operations\n",
    "\n",
    "In addition to constructing graphs node-by-node or edge-by-edge, they\n",
    "can also be generated by\n",
    "\n",
    "1. Applying classic graph operations, such as:\n",
    "\n",
    "   ```\n",
    "   subgraph(G, nbunch)      - induce subgraph of G on nodes in nbunch\n",
    "   union(G1,G2)             - graph union\n",
    "   disjoint_union(G1,G2)    - graph union assuming all nodes are different\n",
    "   cartesian_product(G1,G2) - return Cartesian product graph\n",
    "   compose(G1,G2)           - combine graphs identifying nodes common to both\n",
    "   complement(G)            - graph complement\n",
    "   create_empty_copy(G)     - return an empty copy of the same graph class\n",
    "   convert_to_undirected(G) - return an undirected representation of G\n",
    "   convert_to_directed(G)   - return a directed representation of G\n",
    "   ```\n",
    "\n",
    "1. Using a call to one of the classic small graphs, e.g.,"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [],
   "source": [
    "petersen = nx.petersen_graph()\n",
    "tutte = nx.tutte_graph()\n",
    "maze = nx.sedgewick_maze_graph()\n",
    "tet = nx.tetrahedral_graph()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "1. Using a (constructive) generator for a classic graph, e.g.,"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [],
   "source": [
    "K_5 = nx.complete_graph(5)\n",
    "K_3_5 = nx.complete_bipartite_graph(3, 5)\n",
    "barbell = nx.barbell_graph(10, 10)\n",
    "lollipop = nx.lollipop_graph(10, 20)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "1. Using a stochastic graph generator, e.g.,"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [],
   "source": [
    "er = nx.erdos_renyi_graph(100, 0.15)\n",
    "ws = nx.watts_strogatz_graph(30, 3, 0.1)\n",
    "ba = nx.barabasi_albert_graph(100, 5)\n",
    "red = nx.random_lobster(100, 0.9, 0.9)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "1. Reading a graph stored in a file using common graph formats,\n",
    "   such as edge lists, adjacency lists, GML, GraphML, pickle, LEDA and others."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [],
   "source": [
    "nx.write_gml(red, \"path.to.file\")\n",
    "mygraph = nx.read_gml(\"path.to.file\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "For details on graph formats see Reading and writing graphs\n",
    "and for graph generator functions see Graph generators\n",
    "\n",
    "# Analyzing graphs\n",
    "\n",
    "The structure of `G` can be analyzed using various graph-theoretic\n",
    "functions such as:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{1: 0, 2: 0, 3: 0, 'spam': 0}"
      ]
     },
     "execution_count": 31,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "G = nx.Graph()\n",
    "G.add_edges_from([(1, 2), (1, 3)])\n",
    "G.add_node(\"spam\")       # adds node \"spam\"\n",
    "list(nx.connected_components(G))\n",
    "sorted(d for n, d in G.degree())\n",
    "nx.clustering(G)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Functions that return node properties return iterators over node, value\n",
    "2-tuples. These are easily stored in a [dict](https://docs.python.org/2/library/stdtypes.html#dict) structure if you desire."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{1: 2, 2: 1, 3: 1, 'spam': 0}"
      ]
     },
     "execution_count": 32,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dict(nx.degree(G))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "For values of specific nodes, you can provide a single node or an nbunch of\n",
    "nodes as argument.  If a single node is specified, then a single value is\n",
    "returned.  If an nbunch is specified, then the function will return a\n",
    "dictionary."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[0, 1, 1, 2]"
      ]
     },
     "execution_count": 33,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "nx.degree(G, 1)\n",
    "G.degree(1)\n",
    "dict(G.degree([1, 2]))\n",
    "sorted(d for n, d in G.degree([1, 2]))\n",
    "sorted(d for n, d in G.degree())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "See Algorithms for details on graph algorithms\n",
    "supported.\n",
    "\n",
    "# Drawing graphs\n",
    "\n",
    "NetworkX is not primarily a graph drawing package but basic drawing with\n",
    "Matplotlib as well as an interface to use the open source Graphviz software\n",
    "package are included.  These are part of the `networkx.drawing` module and will\n",
    "be imported if possible.\n",
    "\n",
    "First import Matplotlib's plot interface (pylab works too)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "You may find it useful to interactively test code using `ipython -pylab`,\n",
    "which combines the power of ipython and matplotlib and provides a convenient\n",
    "interactive mode.\n",
    "\n",
    "To test if the import of `networkx.drawing` was successful draw `G` using one of"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {},
   "outputs": [],
   "source": [
    "G = nx.petersen_graph()\n",
    "plt.subplot(121)\n",
    "nx.draw(G, with_labels=True, font_weight='bold')\n",
    "plt.subplot(122)\n",
    "nx.draw_shell(G, nlist=[range(5, 10), range(5)], with_labels=True, font_weight='bold')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "when drawing to an interactive display.  Note that you may need to issue a\n",
    "Matplotlib"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAD8CAYAAABzTgP2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXdYFFcXh39ro2iUjqgIGMDeS+w1sWKvwRh7LBGIjYAl\nihpQ7AYL1mg0ElEhavSLRrFXYhcrRVAjIAJSlrZ7vj8GlLLLtpmdXZj3efZhd+bec88Oc+fsvffc\nc0RERBAQEBAQEMinAt8KCAgICAjoFoJhEBAQEBAogmAYBAQEBASKIBgGAQEBAYEiCIZBQEBAQKAI\ngmEQEBAQECiCYBgEBAQEBIogGAYBAQEBgSIIhkFAQEBAoAiCYRAQEBAQKIJgGAQEBAQEiiAYBgEB\nAQGBIgiGQUBAQECgCIJhEBAQEBAogmAYBAQEBASKIBgGAQEBAYEiCIZBQEBAQKAIgmEQEBAQECiC\nYBgEBAQEBIpQiW8FBOSQkAD8+itw/z6QmgrUqAE0awZMnAhYWvKtnYCA/iH0KaURERHxrYRAIW7d\nAvz8gFOnmM9ZWZ/OGRkBREC/foC3N9C2LT86CgjoE0KfUhnBMOgSW7cC8+YBYjFzs8pDJGJu6DVr\ngBkztKefgIC+IfQptRAMg65QcANnZipfx9hYuJEFBOQh9Cm1EQyDLnDrFtC9e5Eb+C6AHwGEAxAD\nsAcwC8DM4nWNjYELF4A2bbSiqoCAXiCjTwHABgC7ATwCIAWwBMDS4nWFPiV4JekEfn7MULcQQwCc\nBvA5gOEAngD4HkBY8bpiMVNfQEDgEzL6FAD8C8AMgG1pdYU+JRgG3klIYBbFCg3ccgHE5b/fDeA3\nAK3yP8cUr08EnDwJJCZyrKiAgJ4go08V8BuA8wBalFZf6FOCYeCdX38tcagyAI/895MAjANwG0Bz\nAENlyRCJZMoRECiXsNEXynmfEgwD39y/X9R9Lp8hYNYVbgHYD2bDyRAAn8mSIRYDDx5wp6OAgD4h\np0+pRDnvU8IGN75JTS1xKAlAPwCZAC4BaAygDwAfAFaQsQANAMnJ7OolbAYS0Fdk9Cm1YLtP6RGC\nYeCbGjVKHIoGYxQqA2gLwABAQzCjh8fy5JiasqNPaZuBjh4FliwRNgMJ6DYy+pRasNWn9BBhKolv\nmjUDDA2LHGoIxnMiF0AvAOMBHMw/11mWDCMjoGlTzXXZupVx8QsNZQxC8eG4WMwcCw1lym3dqnmb\nAgJs06wZqFifKmAngAlg1uwAIDT/c2jxgmz1KT1FMAx8M2FCiUNVAZwE8CWYEUIwAEcwPtijZckg\nkilHJQpvBlK0tYWIKTdvnmAcBHSKzMxM7BWJkJOdLfP8ZQB78cnr717+57vFC7LRp/QYwTDwjZUV\nMzUjEhU5/AWAM2DWGzIBROCTp1IRRCKgf3/N5v1v3SqxQzQCwEAAFmAWvIcAeFm8XoFxCA9Xv20B\nARZ48+YNFi5cCDs7Oxy5dAkfOnUCFetTAPArAJLxWlq4EBt9Ss8RDIMu4O3NDF3VwciIqa8JxTYD\npQD4CsAJAB0A9AbwJ4ABYHaLFkHYDCTAI7dv38a4cePQuHFjpKam4urVqzh27Bgs162DiM8+pecI\nhkEXaNsWL2bMgAoRXRgK4rposnVfxmagKwDegHGXPQ7gCJg9FI8AhBSvL2wGEtAyEokEoaGh6Nat\nGwYPHoymTZsiKioKAQEBcHJyYgq1bcv0DWNj1YSz0afKAIJh0AFiYmLQ5cABRM+axdyYMobAhZGK\nROwF+5Kxiadg2S4JQBSAV2AMBcDMyZagnG8GEtAOaWlp2LRpE5ydneHr64sZM2YgKioKnp6eMJXl\nQTRjxifjoKBPSQCQEEDvEyTAKx8+fKAmTZrQhg0bmAO3bhENG0ZkaEhkZETE/CZnXkZGJDEwoGOV\nK1PiqVPsKDB2bNE2AMoFqJPsqViaWqzsx9e4cezoIyBQjJiYGJo7dy6ZmZnR8OHD6fLlyySVSpUX\nkN+npAYGlFn8vjUyIqmhIYWZmdGJpUu5+xJ6hjBi4BGJRAJXV1d07NgR7u7uzME2bYAjR4DYWMDH\nBxg3Dnfq1MHz9u0BHx9UiIvD5dmzMf+PP9hRQsZmoEpggvXtB7AAwDYArvnn5C7HlePNQALccO3a\nNYwaNQotW7aEVCpFeHg4Dh8+jE6dOkGkYARQhPw+dXTDBhxq0gQYNw5wcWH++vhAFBuLSn/+ie/3\n7EG2HG+mcgfflqk8M2/ePOrevTvl5OSUWm7p0qW0cOHCj59TU1PJxsaGbt68qbkSMkYMBFB2ofcJ\nAFnkjxjOCSMGAQ7Jzc2loKAg+uKLL8jBwYHWr19PqamprMieNm0arV+/Xu55FxeXUs+XJ4QRA0/s\n3r0boaGhOHz4MCpXrlxqWQcHB0RHR3/8XL16daxYsQIeHh4gTdNpyNkM1BtMuO/JYBae34HxSuoh\nS0Y53wwkoDkpKSlYvXo16tWrh4CAAHh6euL58+f44YcfUL16dVbauHTpErp06SL3vK+vL/z8/PDh\nwwdW2tNnBMPAAxcvXoSXlxeOHz8Oc3NzheWLGwYAmDBhAnJycnDw4EE5tRSTk5ODgwYGMjcDNQOz\nGWgfmKmlHwEclieonG8GElCfFy9ewM3NDQ4ODrh79y6OHj2KS5cuYdiwYahYsSJr7SQlJeHVq1do\n3ry53DJNmzZFnz59sGbNGtba1Vv4HrKUNyIjI8na2ppOnz6tdJ24uDiqWbNmieOXL1+mOnXqUHp6\nuko6pKen04YNG8jW1pa+/PJLiu/UiaQikewpIkUvkYhZLBcQUBKpVEphYWE0aNAgMjc3Jy8vL4qL\ni+O0zdDQUOrdu7fCctHR0WRmZkZv377lVB9dRxgxaJHU1FQMHDgQixcvxldffaV0vVq1aiE5ORni\nYhmpOnXqhM6dO2PVqlVKyUlOTsby5cvh4OCAixcv4siRIzhz5gys1q8XNgMJcE5OTg727duH1q1b\nY9q0aejbty9evnwJPz8/1KlTh9O2L168iK5duyosZ29vj3HjxmH58uWc6qPz8G2Zygu5ubnUt29f\nmjlzplr1nZycKCIiosTxly9fkpmZGcXExMit++bNG5o3bx6ZmprShAkTZMqhLVuIjI1VGy0YGzP1\nBARKITExkZYvX042NjbUq1cvOnHiBEkkEq3q0LZtW7pw4YJSZRMSEsjMzIwiIyM51kp3EUYMWmL+\n/PnIy8vDhg0b1Kova50BAOrWrQt3d3d4enqWOBcZGYlp06ahUaNGyM7Oxt27d7Fnzx40bNiwZAMz\nZuDhhAnIFIlkxpgpjBRgb4OdQJklIiIC3333HZycnBAVFYX//e9/+OeffzBgwABUqKC9R096ejoi\nIiLQrl07pcpbWlrCw8MDixcv5lgzHYZvy1QeCAwMJGdnZ3r//r3aMr777jvavHmzzHMZGRlUt27d\nj7+I7t27R2PGjCFzc3NauHAhJSQkKJT/4cMHqlu3Ll0PCCh1g53UwICOV6lCkX/8ofZ3ESi7SKVS\n+t///kd9+vQhKysrWrJkCe/z9WfOnKFOnTqpVCctLY1q1qxJd+7c4Ugr3UZI1MMxYWFhWLx4MS5f\nvix7276SyBsxAICxsTFWrVqFyZMnw9nZGXfu3MEPP/yAwMBApV39Fi1ahB49euCL778Hvv+eiX30\n669MesPkZCZpSdOmEE2YgGe//Ybt+/fj2KhRan8fgbKFWCzGgQMHsH79elSoUAGzZ89GaGgoDOXk\nRdAmitxUZVGtWjUsXLgQ3t7eOFWQtKo8wbdlKss8e/aMrKys6OzZsxrLCgoKouHDh5c4LpVK6eTJ\nk9SlSxcyMDCgsWPHklgsVkn2tWvXqGbNmvTu3TulymdlZZGdnZ3Sc7YCZZf//vuPFi9eTFZWVtS/\nf386c+aMauEqtED37t3p5MmTKtfLzs4mBwcHCgsLY18pHUcwDByRnJxM9evXp23btrEi78aNG9Sq\nVauPn/Py8igoKIhatGhBTZo0oQMHDtD169fJ2tqaUlJSlJabnZ1NjRs3pqCgIJX02b9/P33xxRc6\n9xAQ0A53796l8ePHk4mJCU2fPp0eP37Mt0oyyc7OpqpVq6rUJwqzf/9+ateuXbm7zwXDwAG5ubn0\n1Vdfkbu7O2syExISyNTUlLKysmj79u3k6OhIHTt2pOPHjxfx8Jg8eTLNnTtXabnLli0jFxcXlW98\niURCLVu2pODgYJXqCegvEomEjh07Rj169KBatWqRr6+v0qNMvrh69Sq1aNFC7foSiYSaN29OR44c\nYVEr3UcwDBwwa9Ys6tOnD+Xm5rIm88OHD1SlShWysbGhPn360IULF2Q+zN++fUvm5ub09OlThTIj\nIiLIwsKCYmNj1dLpzJkz5OjoqDDWk4B+k5aWRgEBAeTk5EStWrWi/fv3U3Z2Nt9qKcWqVavIzc1N\nIxknT56k+vXrs9qfdR3BMLDM5s2bqWHDhmoPXYvz7t07WrJkCVlYWFD16tXp999/V1jH39+fXFxc\nSi0jkUioY8eOFBAQoJF+vXv31liGgG4SGxtLnp6eZG5uTkOGDJH7Y0SXcXFx0XhUK5VKqVu3brRz\n506WtNJ9BMPAImfOnCFra2t68eKFxrLi4uJo9uzZZGpqSpMnT6anT5+Si4sLhYSEKKyblZVFjo6O\ndKqUnA2bN2+mDh06aLzR6M6dO2Rtbc1aBEwB/rlx4waNGTOGTE1Nyd3dnZX7mQ8kEgmZmJiw4i57\n7do1qlOnDmVmZrKgme4jGAaWePLkCVlaWtL58+c1kvP06VOaPHkymZqa0uzZs4vEkHFzc6N169Yp\nJefPP/+kBg0ayJzmiY2NJXNzc3r06JFGuhYwbtw4Wrx4MSuyBPghNzeXgoODqWPHjlS3bl1as2YN\na6Nevrh37x45OTmxJm/o0KHk7+/PmjxdRjAMLJCUlEROTk4aDTVv375NI0eOJAsLC1qyZInMRb11\n69YpPV8qlUrpq6++oo0bN5Y47uLiQktZzFYVExNDZmZm9ObNG9ZkCiggPp5o1Somn4aLC/N31Soi\nJTYzFiYlJYXWrl1LdnZ21KFDBzp06FCZmUsPCAigSZMmsSavYE0uOTmZNZm6imAYNCQnJ4d69uxJ\nc+bMUbmuVCqlCxcuUJ8+fahWrVq0du1aSktLk1s+JCRE4dpBYR4+fEgWFhaUmJj48dgff/xBjRo1\noqysLJX1LY158+bRd999x6pMARncvEk0dCizM93QsMTOdDI0ZM4rSOIUGRlJHh4eZGpqSqNHj6br\n169r6Qtoj9GjR9OePXtYlTlp0iTy8vJiVaYuIhgGDZBKpTR9+nTq378/5eXlqVTv2LFj1LFjR3J0\ndKTt27cr9aC+e/cuNW7cWCUd3dzcaMaMGUTEjGxq1qxJV69eVUmGMrx//54sLS1lB+gTYIeCQIeK\nQqSLRDIDHEqlUrp48SINHTqUzMzMaP78+fTy5Uuevgy3SKVSsrGxYT0QXmxsLJmZmdHr169Zlatr\nCIZBA3755Rdq3Lix0guvubm5dODAAWrSpAk1b96cgoKCVDIoKSkpZGxsrJJnSFJSEllaWtK9e/do\n4sSJGrvulcbq1atp8ODBnMkv1xSKfjsZoAYAVQXIDKB+AD0oJfptTk4O7d+/n9q0aUOOjo70yy+/\nlDoyLQu8ePGCatWqxYkXVXkYHQuGQU3+/vtvqlmzJkVFRSksKxaLaevWrVSvXj3q0qULnTx5Uu0b\n1tTUlOLj41WqExAQQM2bN6e6devShw8f1GpXGcRiMdnZ2dGlS5c4a6NccvNmkZDoAKg9QFMAss//\nXBsgsQzjkFO5MvWztKTu3bvTn3/+qdIPEX1mz549NHr0aE5kJyUlkbm5OT158oQT+bqAYBiKo8Si\nXkREBFlaWip8AKamptKqVavIxsaGBgwYQJcvX9ZYvVatWtGNGzdUqpOamkqVK1emBQsWaNy+Ivbt\n20ft27fXO393nWbo0CLTR+GFHvzR+YYBAP0rwzBIAEru2ZPvb6B1Jk2axOn+Gl9fXxo5ciRn8vlG\nMAwFKLmol3z6NH3++eelLmolJCTQwoULydzcnL7++mu6d+8ea2oOHz5c5bhG8+bNox49epCDg4PK\nAfZUpbyGEOCM+PiS92Oh19N8o1ARoDfy1hwMDVX2VtJ3nJycWO13xcnIyKBatWrRrVu3OGuDTwTD\nQKT0op5UJCJxhQp05KuvZIp5+fIlubm5kampKU2bNo2TjUHz5s0jPz8/pcuHh4eTlZUVxcfH0+DB\ng1Wqqy5///03OTk5CaEy2GDVKrmGIQ2gjvmGYX5pi9FGRkTlxP+eiIn4ampqynmWuG3btlGvXr04\nbYMvhAxuW7cC8+YBmZlMNyoFEREMpVIMvXKFqZfPkydPMHHiRLRo0QIGBgZ4+PAhtm3bhs8//5x1\ndUvLy1Cc3NxcTJkyBatXr4aVlRXWrFmDNWvW4M2bN6zrVZjevXvDzs4OO3fu5LSdcsH9+0BWVonD\niQB6ALgKYCqAUrN+i8VMXo1ywqVLl9CpUyfOs8RNmjQJsbGxOHPmDKft8EH5Ngy3bn0yCsUIAdAW\ngBGAGgA6A0jOPyfKzATmzUPEvn0YPnw4unbtinr16iEyMhKrV69GrVq1OFNZFcOwbt06WFpaYty4\ncQAAR0dHTJkyBQsWLOBMvwL8/f2xbNkypKWlcd5WmSY1tcShl2Dux3AA3gC2Ayg9GSuYZEvlBHUS\n86hD5cqVsWLFCnh5eUEqlXLenlbhe8jCK8UW9Qpev+cPzw0AGp3vHtgUoFfFFvX+MjKiDRs2UHp6\nutZUfvz4MTk6Oios9+zZMzI3Ny/hNfXhwweysbFReQFbHcaOHUtLlizhvJ0yzdixJe7PWvn3Z12A\nPAq9bpQ2nTRuHN/fhDuKOYycMDWl6JkztbKuIpFIqHXr1iqv++k65dcwyFnUkwJkm9/xwkrraABJ\neVjUy8zMJAMDg1LnT6VSKXXv3p3Wrl0r8/zu3bu14jkUHR1NZmZm9N9//3HaTplGxhoD5Lz2lLc1\nhlIcRqQq7ALXlLIYfr78GgY5i3oFXh5GYDYOVQXoc4ACdKjD1axZs0hwveLs3LmT2rZtK9dnveBX\nzv79+7lS8SNz5syh6dOnc95OmUWBV5JSr7LolaThLnC2+fLLL2kLx21ok/JrGGQM0QmgK4V+gdUH\naCJAhvmfQ3RkiN6hQwe6ePGizHNv3rwhS0tLunv3bqkyrly5QrVr1+Z8B2xSUhJZWFiU6c1AnCNn\nylOpl0hENGwY39+AXQrtAlf6xbFxCA8PJxsbG61OK3NJ+V18lrGoBwCWhd7/BmA3gEn5n4/JqsDD\nol5pC9Bubm6YOnUqmjdvXqqMjh07omvXrli1qlR/Fo0xMzPD/Pnz4e3tzWk7ZRpvb8DISL26RkZM\n/bKCHIeR0pxFADDl580DwsM5Uat169bo0qULNm7cyIl8bVN+DUONGjIP2wGoLqdKNVkHTU1ZUkh5\n5BmGkJAQPHjwAIsXL1ZKzqpVq7BlyxbExMSwrGFR3NzcEB4ejitXrnDaTpmlbVtgzRrA2Filapki\nES4MGgS0acORYjzg58e43xbiIIBhAB4AGAxgJIAPAEr4GorFTH2OWLFiBdatW4ekpCTO2tAW5dcw\nNGsGGBqWOFwFwA/5778FM1rYDaAigLHFCxsZAU2bcqikbGQZhpSUFLi5uWH79u0wlPG9ZGFrawsP\nDw94enpyoeZHjIyMsHz5cnh6eoIU7BURkMOMGZ+Mg0iBc6pIBBgbI33pUnx79SoCAgK0oyPXJCQA\np06h8H4jAvBj/vv/AQgCsBPAfQC1i9cnAk6eBBITOVHPyckJI0aMgB+Hxkdr8D2XxRulLOrlAuQF\nUE2AjAFqA9BfOrSod/bsWeratWuRY9OmTVMr4mNGRgbVrVtX48xzisjLy6NmzZoplZpUoBRu3WLW\nDAwNGeeH4s4QhobM+fxQDVFRUWRvb08bNmzgWXEWkOEwopKziBYcRl6/fk2mpqZ6H868/BoGIr1d\n1IuMjKS6det+/HzhwgWqVauW2pmlgoKCqEWLFpxH3jx16hTVr1+/zGQI45WEBOYBN24c/VngBOHv\nL/OHSkxMDDk4OMh1X9YbZDiMqOwsogWHEW9vb5o4cSKnbXBN+TYMxcIZq+zlwFMArZycHKpSpQrl\n5OSQWCwmZ2dnOnr0qNrypFIpde7cmbZv386ilrLb6dWrF23bto3Tdsobygz8X758SZ9//rl+5yx2\ncSnRD58VMgw384/NzP88UV7fVSELojokJyeTpaUlaznV+aD8rjEAai/qwdiYqcfTol7lypVhY2OD\n2NhYrFixAk2aNMHQoUPVlicSibBx40YsXrwYKSkpLGpasp1Vq1bBx8cH6enpnLUjUJK6devi/Pnz\n2L59u/7OgctwGFHZWQTg3GHExMQEnp6eWLhwIaftcArflkkn0LHNMsrQvXt3CgwMJAsLC9bSDE6Z\nMkWt3NWq4urqSj4+Ppy3U15QpRu/fv2anJ2dafny5RxqxBFyNqX+lD9CaFBoKqkiQNd5WGMoQCwW\nk62tLSdpdLWBYBgKKGVRTypjUY9vJkyYQPb29rRjxw7WZL59+1YrmamioqLIzMyM3r59y2k75QVV\nf9+9efOGGjRoQEuXLuVII46Q4zCitLOIlh1Gdu3aRV27dtXLpFUiIsF/sAiJicCvvzJhipOTcezS\nJbSdNAk23t6ApaXC6tqiX79+iIiIQExMDESK3BdVYM2aNTh//jxOnDjBmkxZzJkzB9nZ2di8eTOn\n7ZQHRCKRym7Ab9++Ra9evTB8+HD4+Piweg9xyrBhQGgooM5jSyQChg4FjhxhXy8Z5OXloVmzZliz\nZg369++vlTZZg2fDpPMMGzaMDh06xLcaRYiKiqJq1arRgAEDWJednZ1NTk5OdPLkSdZlF+bdu3dk\nYWFBT58+5bSd8oC63Tg+Pp6aNGlCCxYs0J9ftXrmMHL06FFq1qwZ50mD2KZ8Lz4rQYMGDfDkyRO+\n1fgIEWH69OkYO3YsJzssq1SpgrVr12LOnDnIzc1lXX4B5ubmmDt3rlZyQwjIxsrKCufOncOJEyfg\n7e2tH5sP9cxhZMiQITA2Nsbvv/+u1XY1hm/LpOvs3buXXF1d+VbjI/v27aPmzZtTdHQ0WVtbc9KG\nVCql3r17c74pKjMzk+rUqaO3C3S6gqbdODExkVq0aEFz587Vm5FDXkAAZYpEJNEDh5Hz58+Tvb09\nZWVl8aaDqgiGQQE3btygVq1a8a0GETFDfysrK7p16xZJJBIyNDSkjIwMTtp69OgRWVhYUALHC3W7\nd++mzp07680DSRdh4/ddUlIStWrVin744QfZ/4tiyXBo7FjmM0/hvAMDA2la69YkLcjHoMQucD7p\n168fbdy4kW81lEYwDApISUmhatWq6cSDy9XVlebOnfvxs7OzM6ebaNzd3TnPpZCXl0dNmjSh0NBQ\nTtspy7A18H///j21adOG3N3dP93vpSTDIS0mwylMRkYG1apV61MWQhV2gfPF3bt3ydramj58+MC3\nKkohGAYlsLGxKTUxjjb466+/yMHBoUi89z59+tCJEyc4azMpKYmsrKzo3r17nLVBxHy3Bg0aCKEy\n1ITNGeHk5GRq164dzZw5kySbN+vk/p6VK1fS8OHDZZ7T5dnxsWPH6o2LsO5eRR2iR48edObMGd7a\n//DhA9WtW7eEDtOnT6dffvmF07Y3b95MPXr04HTEJJVKqUePHhQYGMhZG2UZth+GKSkptMrBgbIq\nVVLd64dj4/D+/ftSEz/psmGIjIwkMzMzio+P51sVhQheSUrAt2fSwoUL0bNnT3z55ZdFjtvb28tN\n2MMW3333HRITExESEsJZGyKRCP7+/vDx8UFGRgZn7QgoR41nzzD/7VsY5OUBAM4DEMl5/Vq4IsfJ\ncABg5cqVGDJkCOrXr89ZG1xRr149uLq64ueff+ZbFYUIhkEJ+DQM165dQ3BwMNasWVPinIODA+dJ\ndipVqoQNGzZg3rx5yMrK4qydNm3aoGvXrli3bh1nbQgoiZ8fRIX+13UAeBR6TS5U1LF4XQ6T4bx6\n9Qo7d+7E0qVLOZGvDRYtWoT9+/dz/oNOY/gesugDf//9N/Xq1Uvr7WZnZ1OjRo0oKChI5vmbN29S\ny5YttaLLkCFDyNfXl9M2IiMjydzcXC+G2roEq924lDwlBa9N+bGJWmk57MSUKVPI09Oz1DL68Ej7\n6aef6JtvvuFbjVLR/auoA8TExFDt2rW13q6Pjw+5uLjInd9PTEwkExMTrejy4sULMjc3Zy1gnzw8\nPDxo1qxZnLZR1mD1YSgnUN3HuGEAOeYbhn3yDAMHgeoeP35MFhYW9P79+1LL6YNhSE1N1YpThybo\n/lXUASQSCRkbG2vV1ezRo0dkbm5OsbGxcstIpVKqVq2a2gl6VMXLy4u+/fZbTttITEwkc3Nzevbs\nGaftlCVYfRjKSIZT+HUs3yjYAJRd2kI0y8lwhg8fTitXrlRYTh8MAxHRhg0bqH///nyrIRdhjUEJ\nKlSoAGdnZzx9+lQr7UmlUkydOhU+Pj6wtbWVW04kEsnM/8wVCxYswJkzZ3Djxg3O2rCwsMDcuXP1\nO5a9PpOaWurpDfl/Z4DJjy6X5GSWFAJu3ryJ69evw83NjTWZfDN9+nRERETg4sWLfKsiE8EwKEmD\nBg20Zhi2bdsGIsKMGTMUltWmYfjss8/g6+sLDw8PSKVSztrx8PDA1atXOTVAAnKQkQyngAcAzgEw\nBDBdkRyWkuEQEby8vLBkyRIYqxofSYcxMDDAsmXL4OXlpZMxqgTDoCTa8kyKi4vDTz/9hJ07d6JC\nBcX/Hnt7e849kwrz7bffQiKRcBoUzNjYGMuWLYOnp6dOdpoyTbNmgKGhzFMFo4WxAEoNQG9kBDRt\nyoo6p0+fxuvXrzFx4kRW5OkSrq6uSE9Px7Fjx/hWpQSCYVASbRgGIsLMmTPh5uaGRo0aKVVHmyMG\ngJlW27RpE7y8vJj0nAkJgL8/8M03wMCBzF9/fyavhQaMHz8eSUlJnOeFECjGhAkyD78DUPBT4AcF\nIvLy8pBZD6qKAAAgAElEQVQ1ZozGqkilUnh5eeHnn39GpUqVNJana1SsWBG+vr5YsGABJBIJ3+oU\nhdcVDj3i7t271KRJE07bCAoKokaNGqkUhTE0NJSTvAyKWNSnDz2qX5/TGDonTpygRo0aCaEyFMB6\nNx46VHEYDDkviUhEl6ysyNrampYsWaJRlr7ff/+d2rZtq9Kue317pEmlUurcuTPt3r2bb1WKoF9X\nkUcyMjLI0NCQ8vLyOJH/7t07qlmzpsohqO/du0eNGjXiRCe5bNlCEiMjylP0oNAwho5UKqVu3bqx\nmr60LML6w5CFZDiPHj2iqVOnkomJCU2cOJHu37+vkgrZ2dlUr149Onv2rEr19M0wEBFduXKFbG1t\nSSwW863KR/TvKvKInZ0dvXjxghPZEyZMIDc3N5XrffjwgYyNjbUX/XXLFtUfGhoYhxs3blDt2rWL\nBA8UKAonD0OW/s+JiYm0fPlysrGxoV69etGJEyeUymYWEBBAvXv3VlltfTQMRESDBg2iNWvW8K3G\nR/TzKvIEV9FMz5w5Q3Xr1lV7n4S5ublGQ3alkfNL8luAagFUBSBzgPoAdFvOL0l1GDVqFK1YsYLl\nL1N24OxhWGAcFEwrSZUYGWZlZdHevXupRYsWVL9+fdqyZYtcY5+WlkY1a9ak27dvq6yyvhqGhw8f\nkqWlJaWkpPCtChEJhkElPDw8WLfqGRkZVK9ePY0MTps2bejatWssaiUHOXPP3QD6GqDpADnnb4Cq\nK2taadgwtZp9/vw5mZubc540SF/h9GF46xbzf5OTDCe3UiW6ZGWltNGXSqUUFhZGgwYNIgsLC/Ly\n8qJXr14VKbN8+XIaM2aMWurqq2EgYmYNFi5cyLcaRCQYBpXYunUrTZkyhVWZ8+bNo6+//lojGSNG\njKCDBw+ypJEclIihQwD9m28YKgCUU/y8BjF03Nzc1JpqKw9o5WFYKBkOubh8TIaT/eoV1atXj86d\nO6eyyGfPntGsWbPIxMSExo4dS+Hh4R93vj9//lwtNfXZMLx8+ZLMzMzozZs3fKsiGAZVCAsLoy5d\nuqhWqZSUiOHh4WRlZaVx0Lj58+dzHuBOUQydXwCaUWjEMF9WOQ1i6CQkJGj0wCjL8P0wPHDgALVr\n107tda7379+Tv78/2draUu3atalPnz5qO3nwfS00Zc6cOTRjxgy+1RAMgyr8999/ZGlpqVxhBSkR\npYaG9E+NGnRiyRKN9dqyZQtNnTpVYzmloiCGTrd8gwCA6gAUKq+sBjF0VqxYQaNGjWLxS5UN+H4Y\nSiQSatGiBR05ckQjOc+fP6dq1apRy5YtqV69erRhwwaV1934vhaaoumIiS30+ypqGalUSjVq1KB3\n796VXlDJRTsJQFIWsl6dOnWKvvzyS41kFJCXl0exsbF04cIF2rt3L/n4+NCECRPoipmZwmkkcb5B\nqABQJYBiZJVzcVFbt/T0dKpVqxbd1GJ+YX1AFx6Gp06dovr162u052T8+PEf59ivXr1KI0eOJDMz\nM5ozZw5FR0crJUMXroWmLF++nEaPHs2rDmVvOyGHiESijzGTOnbsKLvQ1q1MFqvMTIXyKgCfsl4B\ngBKxkWShyu5nqVSKt2/fIiYmBtHR0UX+xsTEIC4uDhYWFrC3t4e9vT0cHBzQqVMnOMTFAWfPlpAn\nBhNMrSKYGDp9AVQD8AFAFAC74hU0iKFTtWpV+Pj4YP78+QgLC4NIJFJblgC79OnTBzVr1sSvv/6K\nKVOmqFz/4cOHOHnyJJ4/fw4A6NChAzp06ICXL1/il19+QevWrdGzZ0/Mnj0bHTp0KNP/+9mzZ8PJ\nyQn//vsvWrduzYsOIiIiXlrWUyZMmICuXbti0qRJJU/eugV07y7TKBwE4Jr/3gOf4s58xNgYuHAB\naNNGZZ2ysrJQo0YNZGZmokKFCkhISCjxwC94Hxsbixo1ahR58Bd+X7duXRjKipXj7w8sWQIUy+J2\nPv97dQVgCuASgEdgYuk8B1AkJJuREeDjA8yfr/J3LCAvLw/NmjXD6tWrMWDAALXllCVEIhF0oRtf\nv34dI0aMwPPnz2FkZKRS3cGDB6Nbt26YM2eOzPNpaWnYs2cPNm7cCAsLC8yePRvDhw9H5cqVi5TT\nlWuhKVu2bEFoaChOnz7NS/uCYVARPz8/JCcnw9/fv+TJYcOA0PzZ9UK8AtAUQDqAPMgxDCIRMHQo\ncORIqe0TEZKSkko8+Hft2gVbW1u8fv0aVatWlfvgt7OzUy9KZUICYGdXwjA8A/AdmMibaWAMQicA\nPwFoUlyGoSEQGwtYlhqCTSHHjx+Ht7c37t27h4oVK2okqyygSw/DYcOGoX379vD09FS6zpUrV+Dq\n6oqnT5/K/lFSCIlEguPHj2P9+vWIioqCm5sbpk6dCtP8kaguXQtNyM3NRcOGDREYGIhevXppXwHe\nJrH0lKNHj9LAgQNLnpDjzikFqCdAjQAanb846yFvnj7fnfP9+/d0+/ZtOnLkCK1du5bc3NzIxcWF\nmjRpQtWqVSNTU1Nq2bIlDR06lObMmUObNm2iBg0a0J49e7hNJqRBDB1N9jEURyqVUpcuXWjXrl2s\nyNN3dKkbR0REKJVprYCCWEF79uxRua3w8HD65ptvyNTUlL7//nt69uyZTl0LTTl48KDKsaLYouxc\nRS0RERFBTk5OJU/IcedcB2ZH8B2AxiswDJkiES0yMKDPPvuMmjVrRoMGDSIPDw9av349hYSE0N27\nd+Vmaxs7diz9+uuv3H55FmLosMX169epTp06lJGRwZpMfUXXHoaTJ08mLy8vpcoeP36cGjVqpFEM\nstevX9OCBQvIwsKCANC5c+d4eZiyjUQioZYtW1JwcLDW29atO0oPyM7OJgMDA8rOzi56QoY75wOA\nDABamf9ZkWEggLJGjVLrpl60aBEtYcH1VSFajpVUGiNGjOB+/4YeoGuGIS4ujszMzBTmB8/Ly6Mm\nTZpQaGgoK+1mZGQQAGrQoAG1aNGCfv31V5UiFesi//vf/8jZ2VnrEYaFfAwqUqVKFdjZ2SEyMrLo\nCRkpEY8AyAFwAYALgAKfnmMAvOXIN8jMVMvjQmt5GWbMANasYRbLFekpEjHl1qxR2+OqNHx9fbF2\n7Vq8e/eOddkC6lOnTh1MnjwZPj4+pZY7cOAAqlevjkGDBrHSbsHa2aNHj+Dr64sDBw7A3t4ey5cv\nR6Im+UE4yjmiDL1790bt2rWxe/du7eqhVTNURhg0aBAdPXq06EEZI4Yl+LTpq/irm7xf12puADt3\n7pzqu7I1QUEMHTI0ZM6zOH0ki++//548PDw4bUPX0cVunJSURObm5vTkyROZ57OyssjOzo4uXrzI\narvFr8X9+/dp0qRJZGJiQlOmTKGHDx8qL0zBJlU2co4ow8M9e+ikoSFJDQy0pofu3VF6gKenZ8kp\nDAUhI5SaStIgZER0dDTVqVOHhW+nIoVi6EgGDKDfKlSgXF9ftWMiqUp8fDyZm5tTZGSkVtrTRXTR\nMBAR+fr60ogRI2SeW79+PScJpuRdi/j4eFq6dClZW1tT79696dSpU6VP2Sq5SVXTnCMKyddDooxz\nB4t66OYdpePs3r2bvv3226IHlQgyp9AwaBBkLjc3lypXrsz7nKqTkxM9evRIq20uW7ZM7WicZQFd\nNQwZGRkyd6qnpqaSlZWVysl7lEHRtRCLxbR7925q2rQpNWzYkAIDAykzM7NoIV1ZR+NRD928o3Sc\nK1euULt27Uqe4Nmd097envcYK/3792dtMVFZCkJl3OJ42kpX0VXDQES0bds26tmzZ5Ff5z/99BON\n0yBmVmkoey2kUin9888/NGDAALK0tKSFCxcyUU1leN51kzEV3FjeQ5mte1COB6AYoFkAWQJkCFBH\ngK5zoIfu3lE6zLt376h69eolh6I8u3P26NGDTp8+rZEMTfHw8KDVq1drvd3AwEDq0aMH8z8pJaJt\nWUSXDUNOTg45OTl9vC/fvn1LZmZmSsc+UhV1rsWTJ09oxowZZGJiQrdsbZnEQzIMg0eh10pZfZjF\nvTryfmROy9elCUBjABIB9BlAiSzrobt3lI5jaWkpO246j8O/SZMmUWBgoMZyNGHz5s3cR3qVQW5u\nLo2ws6P/OnTgfbFQ2+iyYSAiOnToELVq1YokEgnNmjWL3N3dOWtLk2vx/skTyq1UqUT/LDAMSvVl\nDaaDPyJnWjoeoMpgglTG5x/7Jl+3JSzrIbirqkmDBg3w5MmTkicKuXNKFQlh2Z1Tay6rpeDk5PQx\nEJo2qbRjBw6+fQvLa9eYsB3FQndALGaOhYYy8ay2btW6juWVESNGoEKFCggICMDvv/+OhQsX8q2S\nTEz//BOVKsmPK2oKwARALwC35BUSiYBff9VMETn1HwHIBVAXgFX+sYLIandZ1kMwDGoi1zAAwIwZ\nyPr7b/xVuTKkBgZM8LjCGBkxcYOGDmUC57Hk419uDUN+RNtK2dlQGDmJ6FNEW8E4aAWRSISVK1di\n4cKF+P7772FlZaW4Eh/cv1/yBwWAz8DsQxoNJlrwOQB9ALyVJUMsBh484ESP+Py/1Qodq5r/t4Qu\nGuohhN1Wk4Lw2/I4FBWFoC+/xMC9exnL/eABkJzMhJ1u2hSYMEHjYHLF0QXDYGtri6SkJGRkZKBq\n1aqKK2jKrVtFwpyfB9BDTtE9ACYUfCgwDm3bqhXRVkA1zM3NkZOT8zHYnU4iY5MqwGxILdjKmQPA\nGcBLAGEAvpZV/rffMPi339RW4xiAgTKOW+f/TS90rOB9TVmCkpPV1kEYMahJqSMGANu2bcO0adOY\nh//8+cC+fcDx48zf+fNZNwoAYG9vz7thqFixIurVq4cXL15op0E/P+bXUT51wESvLXhNLlTUsXhd\nsZipL8A5CxYsgLu7O1avXo2MjAy+1ZFNjRolDmUC+E9OcXkPz0HjxoGY9Vu1XgPHjpUptxGAygBi\n8Wn0UDCl1VxWBQ2MsGAY1KQ0w3Dv3j3ExsZqPV9AzZo1kZaWxnvH09p0UkICcOoUMz2UjyOYkOYF\nr4IO0wpA5+L1iYCTJ7US2qA8c+HCBTx58gQ///wzunTpgg0bSgSd1wmkTZogr0qVIscSADgA6Adg\nOoC2YEYL1mDWGkpgZMTMCGhCs2bMVHMxrMGMeKX5bY8Bk+elGoBZLOshGAY1sbOzQ3x8PDJlJOUJ\nDAzE1KlTS13I4oIKFSrAzs4OMTExWm23OM7Oznj27Bn3DSlYXCMAm/Lf/yCvEBuLhQJyISL8+OOP\nWL58OapUqYIVK1Zg/fr1SEpK4lu1j6SlpWHTpk1ov20bJLm5Rc6ZA/gWTN6RvWB+qQ8BE/fMQpYw\nImaaWBNKqb8RwMx8PUIBtAdwGkweFDb1EAyDmlSsWBGOjo4lHoDp6ekICgpSK70hG+jCOoPWRgxy\nFukKOAHgBQAbMAuHMmFjsVBALqGhoRCLxfj6a2Y23snJCSNHjoSfDkzhvXz5EvPmzYO9vT0uXryI\n9QcOoMrgwUWCQ34GYAeASDBpbN8CCAHQWJZAkQjo31/zaWIrK6BfP5lBKo0AbAaQCCALwFUAHTjQ\nQzAMGiBrOungwYPo2rUrateuzYtO5cowyFksLKBgwmIGmLzUctFgkU5APnl5eViwYAH8/PxQocKn\nR81PP/2EPXv2IDY2lhe9rl27hlGjRqFly5aQSqUIDw/H4cOH0alTJ4gWLCjpRagsRkaAt7y4ySri\n7c2rHoJh0IDihoGIsHXrVkyfPp03nXTBMGhtKknGYmEBD8C4FRqCmRsuFV32lNFj9u7dCysrK/Tr\n16/IcRsbG0yfPh1Lly7Vmi55eXn4448/0L59e4wdOxYdO3ZETEwM1q1bBwcHh08F27b9FFZeFQr2\nI7Hl4ca3HupvzxPYv38/jR49+uPnmzdvkr29PUkkEt50Cg4OpiFDhvDWPhETh6Zq1apys81pSlRU\nFPn7+9PGOnUoU84O1En5O0InK9qpqkFEW11BF7txZmYm1alTh65duybzfHJyMllaWrIecLH4tUhO\nTiZ/f3+ytbWlzp0705EjR5TLFqdD0VUlhoaUp+g+FqKr6g7h4eHUokWLj58nTZpEfn5+PGpEdOvW\nrSI68UXz5s1LRNXUhMjISFq1ahW1adOGLCwsaOrUqRT2xx9MjPpinSQRTIAxgMmix3kIA57RRcPg\n7++v8AfK6tWrafDgway2W3Atnj9/TrNmzSITExNydXVVL8CijuQc8f7qK4po1EireujeHaVHpKWl\nkbGxMUkkEkpOTiYTExN6+/Ytrzq9e/eOatSowasOREQjR46kAwcOaCTjxYsX5OfnR61atSIrKyua\nNm0a/fPPP0XTHPIc0VYX0DXDkJycTBYWFhQREVFqObFYTLa2tnTlyhVW2pVKpQSABg8eTObm5uTl\n5UVxcXGaCy6Uc4RcXJi//v5a+UERHh5ONjY2lJ6erlU9dOuO0kPq1KlDMTEx9Msvv9CoUaP4Voek\nUil99tln9P79e171WLBgAS1dulTles+ePSNfX19q2bIlWVtb04wZM+jcuXPyc97yHNFWF9A1w+Dt\n7U2TJk1SquyuXbuoS5cuauU5LyA7O5v27dtHLVu2JAC0ZcsW5kFaBujduzdt3rxZ6+3q1h2lh3z5\n5Zd06tQpaty4MZ07d45vdYiIqFmzZvTvv//yqsOePXvI1dVVqbJPnz6lFStWUPPmzcna2ppmzpxJ\nYWFhys0FE+lOYhWe0CXD8Pr1azIzM6PY2Filyufm5lLDhg3pr7/+UrmtxMREWrFiBdnY2FCvXr3o\nxIkTOnUtNOXs2bP0+eefU05OjtbbLjtXkSe+//57cnd3J2dnZ41+9bDJ4MGD6fDhw7zqcOXKFWrb\ntq3c848fP6Zly5ZR06ZNycbGhmbNmkUXLlxQ3hgUR1cWC3lAlx6G06ZNo7lz56pUJyQkhJo1a6a0\n00ZERAR99913ZGJiQhMnTqR79+59PKdL10ITpFIptW3blg4ePMhL+2XjKvLIL7/8Qk5OTrRu3Tq+\nVfnIDz/8wEuynMIkJCSQiYlJEWP56NEj8vHxoSZNmlCtWrXIzc2NLl68yJ4X161bJB06lMQASeTl\nY9DCYqG20ZWH4dOnT8nc3JzevXunUj2pVErt27en3377rdQyf//9N/Xt25esrKxoyZIlMtfzdOVa\naEpwcDC1bNmSNw9HIbqqhtjY2CAqKgrjx4/nW5WP2Nvba2cfQSlYWDABAy5fvoyzZ88iODgYqamp\nGDFiBLZt24YOHToU2fTECm3a4MnPP2PsrVu47e6ulYi2Ap9YtGgR5syZA3Nzc5XqFYTlnjBhAkaO\nHAkDA4OP58RiMQ4cOIANGzZAJBJh9uzZCAkJgaGMWEJlhby8PCxcuBCbNm1iv48oiWAYNOTu3buo\nUqUKzMzM+FblIw4ODjh9+jQvbRMRHj58iODgYGRlZWHkyJFwdXXFjh070L59e85v9LCwMLT46ism\ngq2A1ggPD8fly5exZ88etep369YNDRs2RGBgINzd3fH27Vts2bIFgYGBaNOmDTZs2IBevXpBJCNM\nRFlj9+7dqF27Nnr37s2bDoJh0AAiQlBQEAAgNTUVNUrZiatNtL37mYhw//59HD58GMHBwRCLxRgx\nYgS6du2KMWPGYOLEiVrT5dy5cxgyZIjW2hNg8Pb2xk8//aRRDg4/Pz/06tULN27cwMmTJzFmzBhc\nuHABDRo0YFFT3SYzMxM+Pj4ICQnh1QgKhkEDwsLCYGRkhEaNGuHp06do164d3yoBAOpVq4Yhz5+D\nvvkGotRUJnREs2bAxImsTaUQEe7du4fg4GAEBwcjJycHI0eOxL59+9C2bVuIRCL4+PhoLy8DAKlU\nivPnz2Pjxo1aa1MA+Oeff/Dy5UtMnjxZcWEZSKVS/PXXX1i/fj0yMzPx9u1bvHjxQuUpqbLApk2b\n0KFDB/6fJbysbJQRRo4cSQEBATR27Fjau3cv3+owPv1DhxIZGpYMFVGw+Dp0KFNODaRSKf3777/k\n7e1Njo6O5ODgQJ6ennTz5k2ZHlkHDhygkSNHavqtlObOnTvk7OystfZ0BT67sUQiodatW9Mff/yh\nct309HQKCAggJycnatWqFe3fv5+ePHlCZmZmFB8fr5Y++vxIe//+PVlYWNCTJ0/4VkVYfFaXt2/f\n4syZM9ixYweSk5NLzeamFfLzHkMsBohQIi5jQZaz0FDg77+ZQFtK5JomIty+fRvBwcE4fPgwiAgj\nR45EUFAQWrVqVepwV9v5n8+dO4eePXtqrb1yRUICk7fi/n0mqm3+KPS4mRmICCNGjFBaVFxcHAIC\nArBr1y506dIFO3fuRJcuXT7eS2PHjsXPP/9c7kZ+K1euxNChQ1G/fn2+VdFj88ozP//8M02ZMoWI\niA4dOkRDhw7lTxmWN3hJpVK6efMmzZ8/nxwcHMjR0ZG8vb3p9u3bKu3VSE5OpqpVq2ptf4eLiwsd\nOnRIK23pEpx240KjUCrmAiw1MiKxSERvO3VSahR648YNGjNmDJmampK7uzu9ePFCZrn4+HgyMzOj\nqKgoldXV10daXFwcmZmZ0atXr/hWhYiEfQxqkZeXR/b29hQeHk5ERPfv36eGDRvyo4yMkBDd8gPI\nFX41VhASQiqV0vXr12nevHlkb29PTk5OtGDBArpz545GD3ZLS0t6/fo1W99WLrm5uVSjRg1K0POA\neOrA2cOQhU2Dubm5FBwcTB07dqS6devSmjVrKCUlRWHTS5YsoW+++UZllfXVMEyZMoU8PT35VuMj\n+nkVeebkyZPUunXrj5/FYjEZGBjwsnVdVhC5AsPgUei1UkaHlopE9K57d5ozZw7VrVuX6tevT4sW\nLaJ79+6x9iu/Y8eOdP78eVZklcb169epadOmnLeji3DyMNRwFJqSkkJr164le3t76tChAx06dEh+\nvCsZpKamkpWVVZFdzcqgj4bh8ePHZGFhwXt8s8Lo31XUAQYNGkQ7duwocqxevXr09OlT7SoSH19i\neF/YMCjTmbNEIlo5dy7dv3+fkymfCRMm0Pbt21mXWxxfX1/y8PDgvB1dhPWHoYxR6B2AegNkBpAR\nQA0B2izjfpIYGdGa/Omi0aNH0/Xr19VWY8OGDdS/f3+V6uijYRg2bBitXLmSbzWKIGRwU5G4uDhc\nunQJY8aMKXJcVppPzlGQxN4UgAmAXgBuySljYGiIH62t0bRpU078prW1AB0WFiYsPLOFn98nZ4V8\nhoBJOv85gOEAngD4HkBYsaokFqP3v//i7t27CAoKwhdffKG2GtOnT0dERAQuXryotgxd58aNG7hx\n4wbc3Nz4VqUIgmFQkZ07d8LV1RXVqlUrcpwXw3D/PpCVVeLwZwBcAIwGYAcmxWUfMInMSyAWM6Ej\nOEIbhiE7OxvXrl1D165dOW2nXJCQAJw6xfz+zycXQFz++90AfgPQKv9zTLHqFQE0jYtDXXXzFRfC\nwMAAy5Ytw48//ggqpE9ZgYjg5eWFJUuWwFjVFJ4cIxgGFcjLy8OuXbswbdq0EucaNGiAp0+faleh\n1FSZh48BOA5gG5iRgh2AZJT8dfeR5GQOlGPQRv7nGzduoEGDBjAxMeG0nXKBjFFoZQAe+e8nARgH\n4DaA5gCGypIhEikczSqLq6srMjIycOzYMVbk6RKnT5/Gf//9p9XIAMoiGAYVOHHiBOzs7NC0adMS\n53gZMcgIwZEJ4D85xeX+s01NWVKoJI6OjoiKioJUKuWsDWH/AovIGYUOAWAP5ofGfjAhE4aAGZ2W\ngMVRaMWKFeHn54cFCxZAIpGwIlMXkEql8PLyws8//4xKlXRvO5lgGFRg27ZtmD59usxz9evXx+PH\nj7U75G3WDCgWZTIBgAOAfgCmA2gL4CUAazBrDSUwMmIij3JE1apVYWZmhri4OMWF1URYX2ARGaPQ\nJDD3UwyASwDeA2gBwAdAoDw5LI5C+/fvD3Nzc+zbt481mXzzxx9/oHLlyhg2bBjfqshEMAxKEhUV\nhX///VfuDk/L/BhE7969055SEyaUOGQO4FsAzwDsBRAP5pfdWQAWsmQQyZTDJlxOJ2VmZuLff/9F\n586dOZFf7pAxCo0GMxKtDOaHhimAhvnnHsuTw+IotCAs95IlS5AlYzSjb+Tk5GDRokVYuXKlzkaL\nFQyDkuzYsQPjxo2DkZxFNZFIpP3pJCsrZHTtisID7M8A7AAQCUAMZsE5BEBjWfVFIqB/f85zFHC5\nAH3lyhW0aNFCo6ieAoWQMQptCMAMzCJ0LwDjARzMPyfTHHMwCu3YsSNatmyJzZs3syqXD3bs2AFH\nR0edHuUKhkEJcnJysHv3bpmLzoXRtmH4888/MfTmTUgrV1ZPgJER4O3NrlIy4NIwCOsLLCNj9FgV\nwEkAX4IZIQQDcASwAYznWwk4GoX6+vpi1apVSJXjdKEPpKenY8WKFVi5ciXfqpSKYBiUICQkBI0b\nN1YY3EpbhiEnJwdz5syBh4cHlp06hcobNwKqursZGzOB9Nq04UbJQnA5lSSsL7CMlRXQrx+o2BTH\nFwDOgFlvyAQQgU+eSkXgcBTauHFjDBgwAP7+/qzL1hbr169H9+7d0bJlS75VKR2eN9jpBd27d1cq\nrPCff/6p8k5NVYmOjqZ27drRwIEDKSkp6dMJFuLacMWjR4/IycmJdbkpKSlUrVo1EovFrMvWJ9ju\nxrlXr1JWxYqqhcOQEX+LC16+fElmZmb05s0bmed1+ZGWmJhI5ubm9Pz5c75VUYjuXkUd4fHjx2Rt\nbU3Z2dkKyz59+pTq1avHmS4hISFkaWlJa9eulR2+4tYtJtm9oSGTf6Fwhy3IxzBsGKcdVxZisZiq\nVKnCeiyp48ePU8+ePVmVqY+w+TDMycmhESNG0KZGjUjKYsReNpkzZw7NmDFD5jldNgyzZ8+mmTNn\n8q2GUoiIyuCWQhaZM2cODAwM4Ofnp7Bsbm4uqlevjuTkZFaTlefk5ODHH39ESEgIgoKC0L59+9Ir\nJD7YFOsAABnASURBVCYyG4wePGDcBk1NmcXACRM4X2iWR0EeaicnJ9ZkFiSeX7hwIWsy9RGRSMSK\nm3ROTg7GjBmD7OxsHDlyBIZ79hTJ8VGKAsx6lZI5PjQlKSkJ9evXx/Xr1+Ho6FhMFXauBdu8fPkS\nrVq1wqNHj1CzZk2+1VEMv3ZJt8nMzCRzc3OKjIxUuk7Dhg3pwYMHrOkQHR1Nbdu2LTl1pGd89dVX\n9Ndff7Eqs0WLFnT16lVWZeojbHTj7OxsGjx4MA0cOJCysrI+ndDRUeiKFSto9OjRzIf4eKJVq4jG\njqVjANHYscxnHQrBPn78eFq0aBHfaiiNYBhKYe/evdS3b1+V6gwZMoSCg4NZaV/h1JEeMXPmTNqw\nYQNr8t69e0fVq1fnJ9S5jqGpYcjKyiIXFxcaMmSI/CnThAQif3+iceOIXFyYv/7+vD1809PTqa+5\nOb3v3l1mEiE2UtmyxYMHD8jS0lKpPBS6gu7txdYhtm3bBk9PT5XqsOGZVHjq6NixY4qnjvQAtj2T\nzp8/j06dOqGyuq66AgCArKwsDB8+HIaGhggKCpJ/PS0tgfnztatcKVTdtw/H0tJQ8fx52QXUTGXL\nBQsWLICXlxdqyNg8qKsI7qpyuH//PmJjY+Hi4qJSPU0NQ3R0NDp37oyoqCjcvn27TBgFgP29DML+\nBc0Ri8UYMmQIqlatWrpR0DXy85tXzslR/AAjAjIzmbWSrVu1oV0RLl++jHv37mHmzJlab1sTBMMg\nh8DAQEyZMkXlAFeaGIbQ0FB88cUX+PrrrxEaGgozMzO15OgiTk5OrI4YhP0LmpGZmYnBgwfD1NQU\nv//+u/4YhVu3mId8ZubHQxEABoIJ+fIZmBAwL4vXKzAO4eHa0vRjWG0fHx9WnVG0At9zWbpIWloa\nmZqaUlxcnMp1k5OTqVq1aiqtCWRnZ5OHhwfZ2dlplPFKl8nJyaEqVaqwsufgzZs3ZGpqSnl5eSxo\npv+o2o0zMjKoZ8+e5OrqqlK6TZ2gWCrbZIBq5WcsdAFoGD7lOJfI2sMzbJjWVD1+/Dg1btxYL+9T\nYcQgg4MHD6JLly6oU6eOynVNTExQrVo1vH79WqnyBVNH0dHRuH37tkYZr3SZypUrw87ODlFRURrL\nCgsLQ7du3VCxYkUWNCtfZGRkYMCAAahduzb27dunkyGf5SIjidAVAG/AhAQ/DuAImDwRj8DECCsC\nEXDyJOPOzTESiQTe3t7w9fXVy/tUMAwyCAwMlBteWxmUnU4KCQkps1NHsmBrAVpYX1CP9PR09OvX\nD/b29tizZ4/+PbBkJP8pmKBJAhAF4BUYQwEA92TJYDGJUGkcOHAA1atXx8CBAzlviwsEw1CM8PBw\nJCUloXfv3mrLUJTNLScnBz/88APmzJmD48ePY/bs2TobfpdN2FqAFtYXVCctLQ19+/aFs7Mzdu3a\npX9GAZCZRKgbgE4A0sDko7YFUDAe4COVLcCkmv3pp590Oqy2IvRoHKkdtm3bhu+++06jjlO/fn25\nI4bo6GiMHj0atWrVwu3bt2HKYfY0XcPJyQl3797VSEZMTAzS09PRqFEjlrQq+6SmpqJfv35o2rQp\ntm7digoV9PT3oIyoqpXApKw9BGYRui6AiwB+ByB3jz+HqWwBYOvWrWjSpAm6dOnCaTtcoqd3CDek\npqbi8OHDGudglTeVVDB15OrqipCQkHJlFAB2ppLCwsLQo0cPvf0lpm1SU1PRp08ftGjRQr+NAiAz\niRDArDaPBfAzgGEATucf/1KeHA773YcPH+Dn56dUCB1dRhgxFGL//v3o3bu3xrFMihuG7OxseHp6\n4tixYzh+/HiZXWBWBBtTScL6Qj4JCcxc+f37OAYA33zDJNmZOPFjPKyUlBT07t0bX3zxBTZt2qT/\nxrRZM+DIkRLTSb3BZC40AXAKwDsAAwD0kCWD41S2a9asQZ8+fWTmhdcr+HaL0hWkUik1adKEzp49\nq7EsiURCRkZG9OHDB4qMjKQ2bdrQ4MGD6f379yxoqr9IJBIyNDSktLQ0tepLpVKqXbu2XoQt5oyb\nNxmXTQVhIFLOnKHWrVuTh4eH3odT+Uh8fMnvDJAbQFYAVQLIFqAfARLLiwBraMhZGI+3b9+SmZkZ\nRUdHcyJfmwiGIZ/Lly+Ts7Mza52oefPm5O/vT5aWlrR+/fqy0zk1pHHjxnTnzh216j59+pRsbW3L\n77VUMueGVCSiTJGIgnv1KnvXqtg+BpVeHO9jmDVrFnl4eHAmX5vo8YQjuxQsOrMx3M7OzkZmZib8\n/f1x4sQJ/PDDD/o/jGcJTaaTzp07V37XF/LDQCAzEyDCeAC1ARiA2fHbF8Cd/KIiIhgRYfi1axBt\n28aXxtzg7c1MB6kDh6lso6KicPDgwTITAl4wDGDiux8/fhzjx4/XWFZUVBQ6d+6MChUqYPz48WjX\nrh0LGpYdNDUM5XJ9QUYYiJdgXDUngZlf/xtMKIjCiHgIA8E5bdsyAfF0LJXt4sWL4e7uDkue8p2w\njWAYAOzduxcDBw6EhYWFRnKOHj2K9u3b45tvvsGSJUvw8mWJiC3lHnU9k6RSKc6fP48ePWQuKZZt\n/Pw+RQvN5zwYl8ytAA7mH3sFILd4XbGYqV+WmDHjk3FQNHoUiT4ZBY6iq969exdnz57F7NmzOZHP\nC3zPZfGNVColJycnunz5stoysrKyyM3Njezt7enGjRtERHT79m1q2rQpW2qWGc6fP08dO3ZUud79\n+/fJ0dGRA410HDkLrgTQLwDNAMg5Pz7QfB4WXHmllCRCEgMDrSUR6tu3L23atInTNrRNuXdXDQsL\nQ5UqVdCxY0e16kdFRWHUqFGwtbUtsmHN2dkZz58/h0Qi0c9dphyh7lRSwfpCuaOU8A2HAVzIf18H\nzA5gmRSEgdChfAqs0KYN475aLJXttSdPIG3cGJ127OA8le358+fx9OlT/Pnnn5y2o23K/VTStm3b\nMH36dLUWNI8cOYL27dtj3LhxOHr0aJENa1WrVoWVlZUwnVQMGxsbZGZmIiUlRaV65XZ9QUYYiALO\nAxADCAUTH2gEZISbBrQSBoJXCpII7dsHHD+OZ4sWYZOBAedGgYjw448/Yvny5ahSpQqnbWmbcm0Y\n4uPjcebMGYwbN06letnZ2XB3d8f8+fPx119/wcPDQ6ZhYSObW1lDJBKpPGqQSCS4ePFi+RwxyAgD\nIQYgyX9vCMYjqRqAPDCB5GTCcRgIXaJr1664ePEiqFAUVi4IDQ1FdnY2vv76a07b4YNybRh2796N\n4cOHq5RyLyoqCp06dcKrV69w+/ZttG3bVm5ZwTDIRlXDcOfOHdSuXRvW1tYcaqWjyLg3b4AJFjcG\nwAwArQF8ABMbqJU8OeUo/Iq9vT0qVqyIyMhIztrIy8vDggUL4Ofnp99hRuRQ9r6RkkgkEmzfvh3T\npk1Tuk7B1NG3336LI0eOwMTEpNTygmGQjaqeSeV2fQFgwkAUy/5VC4AzgDMAdgFIBjASwDkAsn7i\nkKEhp2EgdA2RSIQuXbrg0qVLnLWxd+9eWFtbo2/fvpy1wSfl1jCcPn0aZmZmaKOEX3N2djbc3Nw+\nTh25u7srtSYhGAbZqDpiKLfrCwAwYUKJQ85g1heSAOQAeA0mumgTOSKys7Iw/9GjUkPBlzW4NAxi\nsRhLly7V67Daiii3hqEgGY+if2xkZCQ6deqE169fK5w6Ko5gGGSjimHIycnB1atX0a1bN4610lGs\nrIB+/RT768tDJAL16wdDW1t06dIFAwYMwD///MP5/DvfcGkYAgIC0KZNG7Rv354T+ToBz+6yvBAX\nF0empqYKg7kFBweTpaUlbdy4Ua2YM1KplKpXr05JSUnqqlomSUxMpBo1aih1TS9fvkwtW7bUglY6\nzM2bTIwkdeIDGRt/9OPPzMyk7du3U6NGjahp06a0a9cuVnJw6yISiYRMTU3pzZs3rMpNTk4mCwsL\nioiIYFWurlEuRww7d+7E119/jWrVqsk8XzB15OnpqdLUUXFEIhHq169frobwymBubg6RSIR3794p\nLFuup5EKYCkMhJGREaZOnYqHDx9i9erVOHToEOzt7bF06VLEx8dzoDh/VKhQAZ07d2Z91LBq1SoM\nGjQIDRs2ZFWurlG2DUNCAuDvz8SqHzgQ+OYbSFauxNHAQLmLzpGRkf9v785jojzzOIB/BwpyaKyi\n21Ur4sFgdz1oZV0UaEy96hFSWUmTEiLWVeuxttlKA9K1xAtLunUtu+JVj1h2W5HUA7SReKBcRWps\nlbRaa4FaohBjUpBT5tk/XkZ5mXdgZnhfhnn5fhIzwzvzPjyD7zO/571+P0yfPh1VVVV2HzpSwsNJ\nluy5ZJWBoY2KaSAMBgPmzp2Lr776CufPn0dVVRXGjx+PN998E9d1dL+D2oeTqqqqsHfvXiQnJ6vW\nZq/l7F0WTXSSs77F01M0GgzS6yUlstXMh44++eQT1dIVb926Vbz33nuqtKUnMTEx4uDBg52+p76+\nXvj6+orffvutZzrlCjpJA/GkHoMDaSBqamrE5s2bxbBhw8TMmTNFdna2aG1t1ehD9IyioiIxefJk\n1dpbuXKlWL9+vWrt9Wb6Cww25qwXBoP0vl27RGNjo1izZo0YPXq0KOkQLLorKytLREZGqtqmHiQn\nJ4sNGzZ0+p5z586J0NDQHuqRi6muFiI1VYjYWCEWLpQeU1O7nROpsbFRHD58WAQHB4ugoCCxa9cu\nUVdXp1Kne1ZTU5Po37+/ePjwYbfbunnzphgyZEifOV+or8BgDgp2nJxr9fISW0eOFFFRUapsQB2V\nlZUJo9GoeruuLiMjQ0RHR3f6nqSkpC6DB2nDZDKJCxcuiMjISDFkyBCRkJAg7t696+xu2e2VV14R\n2dnZ3W4nOjpabNu2TYUeuQb9nGNQyFnfCOBvAH4HwBtSkrGvO6zm1tiI+Hv3cCwhocsb1hwxduxY\nVFRUoLm5WfW2XZktN7nx/ILzGAwGzJgxAydOnEBhYSHq6uowYcIExMTEoNSF6juocZ6htLQU+fn5\nWLdunUq96v30ExgUcta/A+DfAJ6DVMSkCMBsSMXC2/N4/BiG7ds16Va/fv3g7++v6e35rigwMBC3\nb9+2ej19bW0tvvvuO4ez3pJ6AgMDkZaWhjt37iA4OBhRUVGIiIjAl19+idbW1q4bcCI1AkNCQgI2\nbtwIX19flXrV++kjMFRXA2fOSAeHzIsAHID0Ac9BKmYSA6AWUrCQEQI4fVpK36sBXplkaeDAgfDx\n8UFVVZXi6/n5+QgJCYG3o2UcSXWDBg1CfHw8fvrpJ6xduxYffvghjEYjdu7cidraWmd3T1FoaCiu\nXbuGhg6TRlvl5uaisrISy5YtU7lnvZs+AoNCzvoySNWs/CEdSgIAc/KLa0ptmHPWayAoKIiBQYG5\nZoUSHkbqvTw8PPD666+juLgYn332GQoKChAQEIB3330X5eXlzu6ejK+vLyZOnIiSkhK71zWZTEhI\nSMCWLVvg4eGhQe96L30EBoWc9ebbddrfwmbeEbyn1IaGOevHjx/Pm9wUdHYvAwODa5g2bRqOHj2K\nq1evwmAwYMqUKYiOjkZhYWGvSbsRERGBS5cu2b1eZmYmDAYDFi9erEGvejd9BAaFnPXmBM117ZaZ\nn//eWjsa5aznoSRlgYGBiiegHz58iB9//BFTp051Qq/IEaNGjcJHH32E8vJyREREIDY2FqGhofj8\n88/R0mJRibpHOXKeoaWlBe+//z62b9+uy7TaXdHHJ1bIWf8HAB4AKvF07+FK2+Nka+1olLPeHBh6\nywyqt7B2KCkvLw/Tpk3TXVWsvmDAgAFYt24dbt26hcTERKSnp2PMmDFITU3FQycVCwoLC0NxcTEe\nP35s8zqffvopAgICMGvWLA171nvpIzAo5Kx/DkAcABOAmZCKmvwP0qGltUpteHtrlrPez88PHh4e\nustH013WDiXxMJLrc3d3x2uvvYa8vDwcP34c169fx9ixY7F27VqHan53h5+fH/z9/XHtmuLZRQuP\nHj3Cpk2bkJKSonHPei99BAaFnPUAsBPAakh7DMcBhAI4C6nSVUfCZLLajhp4OMnSuHHjcOfOHYtL\nHhkY9GXKlCk4cuQIbty4gYEDB2L69OmIjIzEhQsXemwv2lzu0xY7d+5EeHi4TbVa9EofgcFKznpv\nAP8BUAPpZrdCANMUVjcBOGMwYFdmJpqamjTpIgODJV9fX/j5+eGXX355suz+/fv49ddf8eKLLzqx\nZ6SF4cOHY+vWraioqMCCBQuwevVqvPTSSzh8+LBm487M1vMMDx48wMcff4wtW7Zo2p/eTh+BAQAS\nE6XDQQ5w8/HBqN27kZOTg8DAQKSnp6u+oTIwKOt4OOnixYuIiIjAM88848RekZZ8fHywcuVKlJWV\nYdu2bcjIyEBAQAA2b96MGo3uJYqIiEB+fn6XeygpKSlYvHgxjEajJv1wFfoJDN3MWf/HJUuQk5OD\nY8eOITs7G4GBgdi9e7dqAYKBQVnH1Bg8jNR3uLm5Yd68eTh79ixyc3NRUVEBo9GI5cuXo6ysTNXf\n9fzzz6N///6djsHKykocOHAAGzduVPV3uyL9BAZAlZz1U6dORU5ODjIzM3Hy5EkYjUbs2bOn27mO\neJObgupqvHH3LqampT2plzEmKwuzg4Od3TPqYRMmTMD+/ftx8+ZNjBw5ErNmzXpSM0Kt8xAvv/xy\np4eTkpOT8dZbb2H48OGq/D6X5swMfppRMWd9UVGRePXVV4W/v7/YvXu3aGpqcqhLLS0twsvLSzx6\n9Mih9XWlXb2MFk9P2f9PPSBMXl6K9TKo72hsbBQHDx4UkyZNEi+88ILYs2ePqK+v71ab+/btEzEx\nMYqvlZWViaFDh2qSYdkVGYTQ8cX1NTVSmovr16Wb1wYNki5JjYsDhipdm2RdcXExkpOT8f333yMp\nKQlxcXH2XWdfXY1/TpyIpSEhGOzmJt17MWkSsHSp3X1xaenpUhbchgZZbisLBoN0zshKFTLqG4QQ\nOH/+PHbs2IGSkhKsWLECa9aswbBhw+xu63ZhIY7On48NCxdKN8W2G4OLVqxAWFgY1q9fr8GncEFO\nDkwup7CwUMydO1eMGjVK7N27t+s9iHaz40Y3N+W9l74yO3agXoa5mBLRDz/8IFatWiWeffZZERsb\nK65evWrbim1j0OTlJeo7bl/e3qLV01Oc9vISDZcuafsBXAgDg4MKCgrEnDlzngSI5uZmyzc5UE1O\nt0pKLILCMkCMB4QvIAYDYh4grlsLDnaWqiT9evDggUhJSREjRowQM2bMECdOnLBehtTGMdjaF8ag\nHRgYuqmgoEDMnj1bBAQEiH379j0NEJwdyy1aZDE4AYhQQPwVEAFtP48ARINS4IyKcvYnoF6mublZ\nZGRkiJCQEDFu3DiRlpYmamtrn76BY9BhDAwqyc/PfxIgjiclCVOHDXIHICYCwq3tC/CDzjZMvc2O\n79+XDpl1+Kyl7Z7/3PZ3ASC+Ufq7eHl1u54x6ZPJZBKXL18WUVFRws/PT8THx4t7p051GhT+2257\ne7svjEE76etyVScKCwvD2bNnceTIEfjt3QtTuxKjAPANgMEARnbVUEODVI1OT6zUuZjS7rn5YmB3\nAIqnFTWsl0GuzWAwIDw8HFlZWbhy5QpaWlpwJSrKYgya3YWUKsfqLZR6HIN2YmBQWbjRiPDaWrh3\nWH4EwEUAXV6hr3E1OadQqJfRXh2ApW3P/w4rgUHDehmkH6NHj8aOxEQscHNT/HITAJYAGA7gL9Ya\n0eMYtBPzDqhNjVmteXYcH2/T200mE1paWtDc3Kz4aOsyR1/r6v1p5eWwdi9zDYD5AEoBLAfwYWcf\n1Elpm8nFHDoEg5UbXP8FIB/A123PrbJzDOoNA4Paupgd26ShAWdSU/GPL76w6QvYZDLB09MTHh4e\nTx7bP+9smS2v9evXDwMGDHC4rZGJicCpUxYfswLAHAC3ACQC2NbV30WjehmkM1bG4A1I29km2LDn\n3sf3UBkY1KZQTc4Rfw4KQvqOHTZ90bu7u1udIfUK4eFAbq7FYJ0OoApSXe56AO+0LX8DgEXtNg3r\nZZDOWBmDWZDOZeUBuAzg27blJyFlYrY4q9CH91AZGNSmUE3OEYPHjMHgP/1JlbacLi4O+OADi8VV\nbY+VkGpnmAVDITAIoWm9DNIRK2PQfBnSmQ7LfwZQpLRCH95D5clntSlUkwOA/ZAqyl1t+/l428/H\nldrQ2+zYSr0MYeVfXMf1DQZg/vy+lTqEHGdlDCZDvp0taVv+NqQLQ2T0NgbtxMCgNiuz2nwAhwGY\nS9J82/azYrFBPc6Ou1EvA97e0vpEtlBj7OhxDNqBgUFtVmbHh6A8O07uuL5eZ8fdrJeBPlxmkexk\nZQx2dAjSGLS4OkmvY9AO+s6u6ixXrgAzZgBWbrDplI8PkJen3y9CZlelnsAx2C3cY9ACZ8fWrVol\nDbpFi6TjwB0PL3l7S8sXLZLex6BAjuAY7BbuMWiJs+POqVgvg0gRx6BDGBi0Vloq5V05fVra+Boa\nnr7m7S1trPPnSydX+/gshUgTHIN2Y2DoKZwdEzkXx6DNGBiIiEiGJ5+JiEiGgYGIiGQYGIiISIaB\ngYiIZBgYiIhIhoGBiIhkGBiIiEiGgYGIiGQYGIiISIaBgYiIZBgYiIhIhoGBiIhkGBiIiEiGgYGI\niGQYGIiISIaBgYiIZBgYiIhIhoGBiIhkGBiIiEiGgYGIiGQYGIiISIaBgYiIZBgYiIhIhoGBiIhk\nGBiIiEiGgYGIiGQYGIiISIaBgYiIZBgYiIhIhoGBiIhkGBiIiEiGgYGIiGT+D4xFoNfHLEG9AAAA\nAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f3bf769a450>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "command if you are not using matplotlib in interactive mode (see\n",
    "[Matplotlib FAQ](http://matplotlib.org/faq/installing_faq.html#matplotlib-compiled-fine-but-nothing-shows-up-when-i-use-it)\n",
    ")."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {},
   "outputs": [],
   "source": [
    "options = {\n",
    "    'node_color': 'black',\n",
    "    'node_size': 100,\n",
    "    'width': 3,\n",
    "}\n",
    "plt.subplot(221)\n",
    "nx.draw_random(G, **options)\n",
    "plt.subplot(222)\n",
    "nx.draw_circular(G, **options)\n",
    "plt.subplot(223)\n",
    "nx.draw_spectral(G, **options)\n",
    "plt.subplot(224)\n",
    "nx.draw_shell(G, nlist=[range(5,10), range(5)], **options)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "You can find additional options via `draw_networkx()` and layouts\n",
    "via `layout()`.\n",
    "You can use multiple shells with `draw_shell()`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {},
   "outputs": [],
   "source": [
    "G = nx.dodecahedral_graph()\n",
    "shells = [[2, 3, 4, 5, 6], [8, 1, 0, 19, 18, 17, 16, 15, 14, 7], [9, 10, 11, 12, 13]]\n",
    "nx.draw_shell(G, nlist=shells, **options)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "To save drawings to a file, use, for example"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {},
   "outputs": [],
   "source": [
    "nx.draw(G)\n",
    "plt.savefig(\"path.png\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "writes to the file `path.png` in the local directory. If Graphviz and\n",
    "PyGraphviz or pydot, are available on your system, you can also use\n",
    "`nx_agraph.graphviz_layout(G)` or `nx_pydot.graphviz_layout(G)` to get the\n",
    "node positions, or write the graph in dot format for further processing."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "metadata": {},
   "outputs": [],
   "source": [
    "from networkx.drawing.nx_pydot import write_dot\n",
    "pos = nx.nx_agraph.graphviz_layout(G)\n",
    "nx.draw(G, pos=pos)\n",
    "write_dot(G, 'file.dot')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "See Drawing for additional details."
   ]
  }
 ],
 "metadata": {},
 "nbformat": 4,
 "nbformat_minor": 2
}