{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import mpbn\n",
    "from colomoto_jupyter import tabulate"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "mbn = mpbn.MPBooleanNetwork({\n",
    "    \"a\": \"!b\",\n",
    "    \"b\": \"!a\",\n",
    "    \"c\": \"!a & b\"\n",
    "})"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "c0 = dict(a=0, b=0, c=0)\n",
    "c1 = dict(a=1, b=1, c=1)\n",
    "c2 = dict(a=0, b=1, c=0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>a</th>\n",
       "      <th>b</th>\n",
       "      <th>c</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   a  b  c\n",
       "0  0  1  1\n",
       "1  1  0  0"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "a = list(mbn.attractors())\n",
    "tabulate(a)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "True"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "mbn.reachability(c0, c1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "True"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "mbn.reachability(c0, a[0])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "False"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "mbn.reachability(a[1], a[0])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>a</th>\n",
       "      <th>b</th>\n",
       "      <th>c</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   a  b  c\n",
       "0  0  1  1"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ra = mbn.attractors(reachable_from=c2)\n",
    "tabulate(ra)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>a</th>\n",
       "      <th>b</th>\n",
       "      <th>c</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   a  b  c\n",
       "0  0  1  0\n",
       "1  0  1  1"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "tabulate(mbn.reachable_from(c2), reindex=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>a</th>\n",
       "      <th>b</th>\n",
       "      <th>c</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>0</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   a  b  c\n",
       "0  0  0  0\n",
       "1  0  0  1\n",
       "2  0  1  0\n",
       "3  1  1  0\n",
       "4  1  1  1"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "tabulate(mbn.reachable_from(c2, reversed=True), reindex=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "# computing graph layout...\n"
     ]
    },
    {
     "data": {
      "image/svg+xml": [
       "<?xml version=\"1.0\" encoding=\"UTF-8\" standalone=\"no\"?>\n",
       "<!DOCTYPE svg PUBLIC \"-//W3C//DTD SVG 1.1//EN\"\n",
       " \"http://www.w3.org/Graphics/SVG/1.1/DTD/svg11.dtd\">\n",
       "<!-- Generated by graphviz version 2.44.1 (20200629.0846)\n",
       " -->\n",
       "<!-- Pages: 1 -->\n",
       "<svg width=\"516pt\" height=\"404pt\"\n",
       " viewBox=\"0.00 0.00 516.00 404.00\" xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\">\n",
       "<g id=\"graph0\" class=\"graph\" transform=\"scale(1 1) rotate(0) translate(4 400)\">\n",
       "<polygon fill=\"white\" stroke=\"transparent\" points=\"-4,4 -4,-400 512,-400 512,4 -4,4\"/>\n",
       "<!-- 000 -->\n",
       "<g id=\"node1\" class=\"node\">\n",
       "<title>000</title>\n",
       "<ellipse fill=\"none\" stroke=\"black\" cx=\"265\" cy=\"-378\" rx=\"27\" ry=\"18\"/>\n",
       "<text text-anchor=\"middle\" x=\"265\" y=\"-374.3\" font-family=\"Times,serif\" font-size=\"14.00\">000</text>\n",
       "</g>\n",
       "<!-- 000&#45;&gt;000 -->\n",
       "<g id=\"edge2\" class=\"edge\">\n",
       "<title>000&#45;&gt;000</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M284.9,-390.43C297.69,-393.68 310,-389.53 310,-378 310,-369.62 303.5,-365.14 295.04,-364.56\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"294.5,-361.1 284.9,-365.57 295.19,-368.07 294.5,-361.1\"/>\n",
       "</g>\n",
       "<!-- 100 -->\n",
       "<g id=\"node2\" class=\"node\">\n",
       "<title>100</title>\n",
       "<ellipse fill=\"none\" stroke=\"black\" cx=\"152\" cy=\"-18\" rx=\"27\" ry=\"18\"/>\n",
       "<text text-anchor=\"middle\" x=\"152\" y=\"-14.3\" font-family=\"Times,serif\" font-size=\"14.00\">100</text>\n",
       "</g>\n",
       "<!-- 000&#45;&gt;100 -->\n",
       "<g id=\"edge1\" class=\"edge\">\n",
       "<title>000&#45;&gt;100</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M238.13,-375.12C170.76,-368.82 0,-342.26 0,-235 0,-235 0,-235 0,-161 0,-117.07 11.98,-103.1 43,-72 63.8,-51.14 94.45,-37.25 117.96,-28.94\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"119.18,-32.22 127.54,-25.71 116.95,-25.58 119.18,-32.22\"/>\n",
       "</g>\n",
       "<!-- 110 -->\n",
       "<g id=\"node3\" class=\"node\">\n",
       "<title>110</title>\n",
       "<ellipse fill=\"none\" stroke=\"black\" cx=\"265\" cy=\"-306\" rx=\"27\" ry=\"18\"/>\n",
       "<text text-anchor=\"middle\" x=\"265\" y=\"-302.3\" font-family=\"Times,serif\" font-size=\"14.00\">110</text>\n",
       "</g>\n",
       "<!-- 000&#45;&gt;110 -->\n",
       "<g id=\"edge3\" class=\"edge\">\n",
       "<title>000&#45;&gt;110</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M259.16,-360.41C258.3,-352.51 258.05,-342.85 258.41,-333.94\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"261.9,-334.18 259.12,-323.96 254.92,-333.68 261.9,-334.18\"/>\n",
       "</g>\n",
       "<!-- 010 -->\n",
       "<g id=\"node4\" class=\"node\">\n",
       "<title>010</title>\n",
       "<ellipse fill=\"none\" stroke=\"black\" cx=\"389\" cy=\"-90\" rx=\"27\" ry=\"18\"/>\n",
       "<text text-anchor=\"middle\" x=\"389\" y=\"-86.3\" font-family=\"Times,serif\" font-size=\"14.00\">010</text>\n",
       "</g>\n",
       "<!-- 000&#45;&gt;010 -->\n",
       "<g id=\"edge4\" class=\"edge\">\n",
       "<title>000&#45;&gt;010</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M287.83,-368.27C306.71,-359.83 332.66,-345.2 347,-324 390.16,-260.17 392.63,-164.87 390.86,-118.43\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"394.35,-118.06 390.37,-108.24 387.36,-118.4 394.35,-118.06\"/>\n",
       "</g>\n",
       "<!-- 101 -->\n",
       "<g id=\"node5\" class=\"node\">\n",
       "<title>101</title>\n",
       "<ellipse fill=\"none\" stroke=\"black\" cx=\"79\" cy=\"-90\" rx=\"27\" ry=\"18\"/>\n",
       "<text text-anchor=\"middle\" x=\"79\" y=\"-86.3\" font-family=\"Times,serif\" font-size=\"14.00\">101</text>\n",
       "</g>\n",
       "<!-- 000&#45;&gt;101 -->\n",
       "<g id=\"edge5\" class=\"edge\">\n",
       "<title>000&#45;&gt;101</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M244.94,-365.82C212.24,-346.62 147.33,-304.36 110,-252 81.52,-212.05 77.34,-152.92 77.64,-118.54\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"81.15,-118.38 77.88,-108.31 74.15,-118.22 81.15,-118.38\"/>\n",
       "</g>\n",
       "<!-- 111 -->\n",
       "<g id=\"node6\" class=\"node\">\n",
       "<title>111</title>\n",
       "<ellipse fill=\"none\" stroke=\"black\" cx=\"209\" cy=\"-234\" rx=\"27\" ry=\"18\"/>\n",
       "<text text-anchor=\"middle\" x=\"209\" y=\"-230.3\" font-family=\"Times,serif\" font-size=\"14.00\">111</text>\n",
       "</g>\n",
       "<!-- 000&#45;&gt;111 -->\n",
       "<g id=\"edge6\" class=\"edge\">\n",
       "<title>000&#45;&gt;111</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M247.02,-364.31C237.26,-354.04 226.46,-338.9 220,-324 211.41,-304.17 206.29,-280.2 204.75,-262.13\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"208.24,-261.81 204.26,-251.99 201.25,-262.15 208.24,-261.81\"/>\n",
       "</g>\n",
       "<!-- 001 -->\n",
       "<g id=\"node7\" class=\"node\">\n",
       "<title>001</title>\n",
       "<ellipse fill=\"none\" stroke=\"black\" cx=\"321\" cy=\"-162\" rx=\"27\" ry=\"18\"/>\n",
       "<text text-anchor=\"middle\" x=\"321\" y=\"-158.3\" font-family=\"Times,serif\" font-size=\"14.00\">001</text>\n",
       "</g>\n",
       "<!-- 000&#45;&gt;001 -->\n",
       "<g id=\"edge7\" class=\"edge\">\n",
       "<title>000&#45;&gt;001</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M279.17,-362.1C289.9,-352.71 303.4,-339.46 310,-324 328.64,-280.35 323.43,-223.77 320.44,-190.62\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"323.9,-189.99 319.55,-180.33 316.92,-190.59 323.9,-189.99\"/>\n",
       "</g>\n",
       "<!-- 011 -->\n",
       "<g id=\"node8\" class=\"node\">\n",
       "<title>011</title>\n",
       "<ellipse fill=\"none\" stroke=\"black\" cx=\"389\" cy=\"-18\" rx=\"27\" ry=\"18\"/>\n",
       "<text text-anchor=\"middle\" x=\"389\" y=\"-14.3\" font-family=\"Times,serif\" font-size=\"14.00\">011</text>\n",
       "</g>\n",
       "<!-- 000&#45;&gt;011 -->\n",
       "<g id=\"edge8\" class=\"edge\">\n",
       "<title>000&#45;&gt;011</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M291.47,-374.08C354.45,-365.51 508,-334.44 508,-235 508,-235 508,-235 508,-161 508,-104.81 452.58,-59.18 417.39,-35.87\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"419.22,-32.89 408.92,-30.43 415.44,-38.78 419.22,-32.89\"/>\n",
       "</g>\n",
       "<!-- 100&#45;&gt;100 -->\n",
       "<g id=\"edge9\" class=\"edge\">\n",
       "<title>100&#45;&gt;100</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M171.9,-30.43C184.69,-33.68 197,-29.53 197,-18 197,-9.62 190.5,-5.14 182.04,-4.56\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"181.5,-1.1 171.9,-5.57 182.19,-8.07 181.5,-1.1\"/>\n",
       "</g>\n",
       "<!-- 110&#45;&gt;000 -->\n",
       "<g id=\"edge12\" class=\"edge\">\n",
       "<title>110&#45;&gt;000</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M270.88,-323.96C271.71,-331.83 271.95,-341.37 271.58,-350.19\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"268.07,-350.18 270.84,-360.41 275.06,-350.69 268.07,-350.18\"/>\n",
       "</g>\n",
       "<!-- 110&#45;&gt;100 -->\n",
       "<g id=\"edge10\" class=\"edge\">\n",
       "<title>110&#45;&gt;100</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M240.53,-298C218.87,-290.42 188.39,-276.07 173,-252 131.87,-187.7 139.81,-92.24 146.94,-46.02\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"150.44,-46.33 148.62,-35.89 143.53,-45.19 150.44,-46.33\"/>\n",
       "</g>\n",
       "<!-- 110&#45;&gt;110 -->\n",
       "<g id=\"edge11\" class=\"edge\">\n",
       "<title>110&#45;&gt;110</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M284.9,-318.43C297.69,-321.68 310,-317.53 310,-306 310,-297.62 303.5,-293.14 295.04,-292.56\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"294.5,-289.1 284.9,-293.57 295.19,-296.07 294.5,-289.1\"/>\n",
       "</g>\n",
       "<!-- 110&#45;&gt;010 -->\n",
       "<g id=\"edge13\" class=\"edge\">\n",
       "<title>110&#45;&gt;010</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M281.87,-291.74C306.97,-270.94 353.49,-228.01 375,-180 383.73,-160.52 387.13,-136.52 388.4,-118.36\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"391.91,-118.32 388.94,-108.15 384.92,-117.95 391.91,-118.32\"/>\n",
       "</g>\n",
       "<!-- 110&#45;&gt;101 -->\n",
       "<g id=\"edge14\" class=\"edge\">\n",
       "<title>110&#45;&gt;101</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M239.11,-300.58C210.68,-294.43 165.42,-280.66 138,-252 102.19,-214.57 87.8,-153.65 82.26,-118.45\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"85.7,-117.81 80.81,-108.41 78.77,-118.81 85.7,-117.81\"/>\n",
       "</g>\n",
       "<!-- 110&#45;&gt;111 -->\n",
       "<g id=\"edge16\" class=\"edge\">\n",
       "<title>110&#45;&gt;111</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M248.74,-291.5C240.27,-282.47 230.21,-270.33 222.21,-259.61\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"224.99,-257.48 216.31,-251.41 219.31,-261.57 224.99,-257.48\"/>\n",
       "</g>\n",
       "<!-- 110&#45;&gt;001 -->\n",
       "<g id=\"edge15\" class=\"edge\">\n",
       "<title>110&#45;&gt;001</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M267.17,-287.87C274.36,-262.95 292.43,-216.56 305.91,-187.87\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"309.06,-189.4 310.29,-178.88 302.77,-186.34 309.06,-189.4\"/>\n",
       "</g>\n",
       "<!-- 110&#45;&gt;011 -->\n",
       "<g id=\"edge17\" class=\"edge\">\n",
       "<title>110&#45;&gt;011</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M287.2,-295.28C306.27,-286.13 333.67,-271.01 353,-252 406.81,-199.08 422.19,-180.55 443,-108 447.41,-92.62 449.19,-86.75 443,-72 437.29,-58.39 426.19,-46.52 415.53,-37.48\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"417.57,-34.63 407.56,-31.14 413.21,-40.1 417.57,-34.63\"/>\n",
       "</g>\n",
       "<!-- 010&#45;&gt;010 -->\n",
       "<g id=\"edge19\" class=\"edge\">\n",
       "<title>010&#45;&gt;010</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M408.9,-102.43C421.69,-105.68 434,-101.53 434,-90 434,-81.62 427.5,-77.14 419.04,-76.56\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"418.5,-73.1 408.9,-77.57 419.19,-80.07 418.5,-73.1\"/>\n",
       "</g>\n",
       "<!-- 010&#45;&gt;011 -->\n",
       "<g id=\"edge18\" class=\"edge\">\n",
       "<title>010&#45;&gt;011</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M389,-71.7C389,-63.98 389,-54.71 389,-46.11\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"392.5,-46.1 389,-36.1 385.5,-46.1 392.5,-46.1\"/>\n",
       "</g>\n",
       "<!-- 101&#45;&gt;100 -->\n",
       "<g id=\"edge20\" class=\"edge\">\n",
       "<title>101&#45;&gt;100</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M93.77,-74.83C104.16,-64.87 118.24,-51.37 129.96,-40.14\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"132.45,-42.59 137.25,-33.14 127.61,-37.54 132.45,-42.59\"/>\n",
       "</g>\n",
       "<!-- 101&#45;&gt;101 -->\n",
       "<g id=\"edge21\" class=\"edge\">\n",
       "<title>101&#45;&gt;101</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M98.9,-102.43C111.69,-105.68 124,-101.53 124,-90 124,-81.62 117.5,-77.14 109.04,-76.56\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"108.5,-73.1 98.9,-77.57 109.19,-80.07 108.5,-73.1\"/>\n",
       "</g>\n",
       "<!-- 111&#45;&gt;000 -->\n",
       "<g id=\"edge24\" class=\"edge\">\n",
       "<title>111&#45;&gt;000</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M217.35,-251.32C222.24,-269.62 227.55,-299.88 238,-324 242.16,-333.59 248.11,-343.29 253.55,-351.72\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"250.74,-353.82 259.12,-360.3 256.61,-350.01 250.74,-353.82\"/>\n",
       "</g>\n",
       "<!-- 111&#45;&gt;100 -->\n",
       "<g id=\"edge22\" class=\"edge\">\n",
       "<title>111&#45;&gt;100</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M204.53,-216.21C194.69,-179.26 171.09,-90.67 159.16,-45.89\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"162.49,-44.77 156.53,-36.01 155.72,-46.57 162.49,-44.77\"/>\n",
       "</g>\n",
       "<!-- 111&#45;&gt;110 -->\n",
       "<g id=\"edge23\" class=\"edge\">\n",
       "<title>111&#45;&gt;110</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M225.4,-248.65C233.82,-257.63 243.76,-269.64 251.7,-280.27\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"248.88,-282.35 257.57,-288.41 254.56,-278.25 248.88,-282.35\"/>\n",
       "</g>\n",
       "<!-- 111&#45;&gt;010 -->\n",
       "<g id=\"edge25\" class=\"edge\">\n",
       "<title>111&#45;&gt;010</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M219.55,-217.2C232.97,-198 257.72,-165.51 285,-144 306.8,-126.81 335.06,-112.84 356.62,-103.56\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"358.13,-106.73 366,-99.64 355.42,-100.27 358.13,-106.73\"/>\n",
       "</g>\n",
       "<!-- 111&#45;&gt;101 -->\n",
       "<g id=\"edge26\" class=\"edge\">\n",
       "<title>111&#45;&gt;101</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M195.58,-218.34C172.81,-193.47 126.47,-142.85 99.58,-113.48\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"101.88,-110.81 92.54,-105.79 96.72,-115.53 101.88,-110.81\"/>\n",
       "</g>\n",
       "<!-- 111&#45;&gt;111 -->\n",
       "<g id=\"edge28\" class=\"edge\">\n",
       "<title>111&#45;&gt;111</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M228.9,-246.43C241.69,-249.68 254,-245.53 254,-234 254,-225.62 247.5,-221.14 239.04,-220.56\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"238.5,-217.1 228.9,-221.57 239.19,-224.07 238.5,-217.1\"/>\n",
       "</g>\n",
       "<!-- 111&#45;&gt;001 -->\n",
       "<g id=\"edge27\" class=\"edge\">\n",
       "<title>111&#45;&gt;001</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M225.14,-219.5C242.33,-207.74 269.54,-190.68 290.68,-178.38\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"292.68,-181.28 299.62,-173.27 289.2,-175.2 292.68,-181.28\"/>\n",
       "</g>\n",
       "<!-- 111&#45;&gt;011 -->\n",
       "<g id=\"edge29\" class=\"edge\">\n",
       "<title>111&#45;&gt;011</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M213.02,-216.07C218.07,-197.23 228.05,-166.48 244,-144 276.52,-98.16 329.34,-58.34 361.46,-36.59\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"363.55,-39.4 369.94,-30.94 359.67,-33.57 363.55,-39.4\"/>\n",
       "</g>\n",
       "<!-- 001&#45;&gt;000 -->\n",
       "<g id=\"edge31\" class=\"edge\">\n",
       "<title>001&#45;&gt;000</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M330,-179.37C339.76,-209.44 348.97,-274.9 328,-324 321.55,-339.12 308.5,-352.12 295.8,-361.46\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"293.46,-358.82 287.14,-367.32 297.39,-364.61 293.46,-358.82\"/>\n",
       "</g>\n",
       "<!-- 001&#45;&gt;100 -->\n",
       "<g id=\"edge30\" class=\"edge\">\n",
       "<title>001&#45;&gt;100</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M304.8,-147.38C274.89,-122.26 210.9,-68.49 175.96,-39.14\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"178.06,-36.33 168.16,-32.57 173.56,-41.69 178.06,-36.33\"/>\n",
       "</g>\n",
       "<!-- 001&#45;&gt;110 -->\n",
       "<g id=\"edge32\" class=\"edge\">\n",
       "<title>001&#45;&gt;110</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M318.83,-180.13C311.65,-204.99 293.66,-251.2 280.19,-279.9\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"276.87,-278.74 275.64,-289.26 283.16,-281.8 276.87,-278.74\"/>\n",
       "</g>\n",
       "<!-- 001&#45;&gt;010 -->\n",
       "<g id=\"edge33\" class=\"edge\">\n",
       "<title>001&#45;&gt;010</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M335.09,-146.5C344.57,-136.73 357.24,-123.69 367.92,-112.7\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"370.47,-115.1 374.93,-105.49 365.45,-110.22 370.47,-115.1\"/>\n",
       "</g>\n",
       "<!-- 001&#45;&gt;101 -->\n",
       "<g id=\"edge34\" class=\"edge\">\n",
       "<title>001&#45;&gt;101</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M296.46,-153.9C253.05,-141.34 162.62,-115.19 113.14,-100.88\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"113.82,-97.43 103.24,-98.01 111.87,-104.15 113.82,-97.43\"/>\n",
       "</g>\n",
       "<!-- 001&#45;&gt;111 -->\n",
       "<g id=\"edge35\" class=\"edge\">\n",
       "<title>001&#45;&gt;111</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M304.78,-176.56C287.54,-188.35 260.27,-205.44 239.14,-217.72\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"237.16,-214.82 230.21,-222.82 240.63,-220.9 237.16,-214.82\"/>\n",
       "</g>\n",
       "<!-- 001&#45;&gt;001 -->\n",
       "<g id=\"edge36\" class=\"edge\">\n",
       "<title>001&#45;&gt;001</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M340.9,-174.43C353.69,-177.68 366,-173.53 366,-162 366,-153.62 359.5,-149.14 351.04,-148.56\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"350.5,-145.1 340.9,-149.57 351.19,-152.07 350.5,-145.1\"/>\n",
       "</g>\n",
       "<!-- 001&#45;&gt;011 -->\n",
       "<g id=\"edge37\" class=\"edge\">\n",
       "<title>001&#45;&gt;011</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M325.69,-144.08C331.09,-125.77 340.78,-96 353,-72 358.07,-62.04 364.71,-51.73 370.9,-42.88\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"373.87,-44.75 376.87,-34.59 368.19,-40.66 373.87,-44.75\"/>\n",
       "</g>\n",
       "<!-- 011&#45;&gt;011 -->\n",
       "<g id=\"edge38\" class=\"edge\">\n",
       "<title>011&#45;&gt;011</title>\n",
       "<path fill=\"none\" stroke=\"black\" d=\"M408.9,-30.43C421.69,-33.68 434,-29.53 434,-18 434,-9.62 427.5,-5.14 419.04,-4.56\"/>\n",
       "<polygon fill=\"black\" stroke=\"black\" points=\"418.5,-1.1 408.9,-5.57 419.19,-8.07 418.5,-1.1\"/>\n",
       "</g>\n",
       "</g>\n",
       "</svg>\n"
      ],
      "text/plain": [
       "<networkx.classes.digraph.DiGraph at 0x7fdba5628520>"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "mbn.dynamics()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}