From 59dae0bb2f96ebfd0493640d8898a0473dab5600 Mon Sep 17 00:00:00 2001
From: nbgitpuller <nbgitpuller@example.com>
Date: Wed, 7 Jul 2021 04:12:57 +0000
Subject: [PATCH] WIP

---
 .../4_1_PO_multilinear_unloading.ipynb        |  85 +++++++-----
 tour4_plastic_bond/4_2_BS_EP_SH_I_A.ipynb     | 128 +++++++++---------
 2 files changed, 117 insertions(+), 96 deletions(-)

diff --git a/tour4_plastic_bond/4_1_PO_multilinear_unloading.ipynb b/tour4_plastic_bond/4_1_PO_multilinear_unloading.ipynb
index ccbedd5..89c5ddc 100644
--- a/tour4_plastic_bond/4_1_PO_multilinear_unloading.ipynb
+++ b/tour4_plastic_bond/4_1_PO_multilinear_unloading.ipynb
@@ -42,8 +42,7 @@
    "cell_type": "markdown",
    "metadata": {},
    "source": [
-    "\n",
-    "To motivate a more general development of the model for the bond-slip behavior let us consider the case of non-monotonic loading. What happens in the material structure of the bond if the load is reduced and than it grows again?"
+    "To motivate a more general development of the model for the bond-slip behavior let us consider the case of non-monotonic loading. What happens in the material structure of the bond if the load is reduced and then it grows again?"
    ]
   },
   {
@@ -126,7 +125,7 @@
     {
      "data": {
       "application/vnd.jupyter.widget-view+json": {
-       "model_id": "be41afc544224f66bbbd0edbd988f4af",
+       "model_id": "600f2ec31e534e6d91b396a93fdbecd6",
        "version_major": 2,
        "version_minor": 0
       },
@@ -153,7 +152,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 3,
+   "execution_count": 9,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -161,7 +160,7 @@
     "                              unloading_ratio=0.0,\n",
     "                               numbe_of_increments=200,\n",
     "                              amplitude_type='constant',\n",
-    "                              loading_range='non-symmetric');"
+    "                              loading_range='non-symmetric');\n"
    ]
   },
   {
@@ -173,22 +172,33 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 4,
+   "execution_count": 8,
    "metadata": {},
    "outputs": [
     {
-     "data": {
-      "application/vnd.jupyter.widget-view+json": {
-       "model_id": "4183aa8efd0e4cf4824f5872c9a9b0d2",
-       "version_major": 2,
-       "version_minor": 0
-      },
-      "text/plain": [
-       "VBox(children=(HBox(children=(VBox(children=(Tree(layout=Layout(align_items='stretch', border='solid 1px black…"
-      ]
-     },
-     "metadata": {},
-     "output_type": "display_data"
+     "ename": "TypeError",
+     "evalue": "get_corr_pred() got an unexpected keyword argument 'kappa_n'",
+     "output_type": "error",
+     "traceback": [
+      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
+      "\u001b[0;31mTypeError\u001b[0m                                 Traceback (most recent call last)",
+      "\u001b[0;32m<ipython-input-8-30f7b3be9a30>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m\u001b[0m\n\u001b[1;32m      1\u001b[0m \u001b[0mpo_trc\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mrun\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 2\u001b[0;31m \u001b[0mpo_trc\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0minteract\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/bmcs_utils/model.py\u001b[0m in \u001b[0;36minteract\u001b[0;34m(self, **kw)\u001b[0m\n\u001b[1;32m     43\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     44\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0minteract\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m**\u001b[0m\u001b[0mkw\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 45\u001b[0;31m         \u001b[0;32mreturn\u001b[0m \u001b[0mAppWindow\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m**\u001b[0m\u001b[0mkw\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0minteract\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m     46\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     47\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0mapp\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m**\u001b[0m\u001b[0mkw\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/bmcs_utils/app_window.py\u001b[0m in \u001b[0;36minteract\u001b[0;34m(self)\u001b[0m\n\u001b[1;32m    133\u001b[0m                              layout=ipw.Layout(align_items=\"stretch\",\n\u001b[1;32m    134\u001b[0m                                                width=\"100%\"))\n\u001b[0;32m--> 135\u001b[0;31m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmodel_tree\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mselected\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;32mTrue\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    136\u001b[0m         \u001b[0mdisplay\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mapp_print\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    137\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/traitlets/traitlets.py\u001b[0m in \u001b[0;36m__set__\u001b[0;34m(self, obj, value)\u001b[0m\n\u001b[1;32m    602\u001b[0m             \u001b[0;32mraise\u001b[0m \u001b[0mTraitError\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m'The \"%s\" trait is read-only.'\u001b[0m \u001b[0;34m%\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mname\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    603\u001b[0m         \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 604\u001b[0;31m             \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mset\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mobj\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mvalue\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    605\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    606\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0m_validate\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mobj\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mvalue\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/traitlets/traitlets.py\u001b[0m in \u001b[0;36mset\u001b[0;34m(self, obj, value)\u001b[0m\n\u001b[1;32m    591\u001b[0m             \u001b[0;31m# we explicitly compare silent to True just in case the equality\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    592\u001b[0m             \u001b[0;31m# comparison above returns something other than True/False\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 593\u001b[0;31m             \u001b[0mobj\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_notify_trait\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mname\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mold_value\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mnew_value\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    594\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    595\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0m__set__\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mobj\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mvalue\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/traitlets/traitlets.py\u001b[0m in \u001b[0;36m_notify_trait\u001b[0;34m(self, name, old_value, new_value)\u001b[0m\n\u001b[1;32m   1215\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1216\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0m_notify_trait\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mname\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mold_value\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mnew_value\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1217\u001b[0;31m         self.notify_change(Bunch(\n\u001b[0m\u001b[1;32m   1218\u001b[0m             \u001b[0mname\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mname\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1219\u001b[0m             \u001b[0mold\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mold_value\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/ipywidgets/widgets/widget.py\u001b[0m in \u001b[0;36mnotify_change\u001b[0;34m(self, change)\u001b[0m\n\u001b[1;32m    604\u001b[0m                 \u001b[0;31m# Send new state to front-end\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    605\u001b[0m                 \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0msend_state\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mkey\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mname\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 606\u001b[0;31m         \u001b[0msuper\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mWidget\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mnotify_change\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mchange\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    607\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    608\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0m__repr__\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/traitlets/traitlets.py\u001b[0m in \u001b[0;36mnotify_change\u001b[0;34m(self, change)\u001b[0m\n\u001b[1;32m   1225\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0mnotify_change\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mchange\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1226\u001b[0m         \u001b[0;34m\"\"\"Notify observers of a change event\"\"\"\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1227\u001b[0;31m         \u001b[0;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_notify_observers\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mchange\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m   1228\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1229\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0m_notify_observers\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mevent\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/traitlets/traitlets.py\u001b[0m in \u001b[0;36m_notify_observers\u001b[0;34m(self, event)\u001b[0m\n\u001b[1;32m   1262\u001b[0m                 \u001b[0mc\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mgetattr\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mc\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mname\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1263\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1264\u001b[0;31m             \u001b[0mc\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mevent\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m   1265\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1266\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0m_add_notifiers\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mhandler\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mname\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mtype\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/bmcs_utils/app_window.py\u001b[0m in \u001b[0;36mselect_node\u001b[0;34m(self, event)\u001b[0m\n\u001b[1;32m    227\u001b[0m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mset_plot_backend\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mbackend\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    228\u001b[0m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0msetup_plot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mcontroller\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmodel\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 229\u001b[0;31m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mupdate_plot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mcontroller\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmodel\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    230\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    231\u001b[0m     \u001b[0mcurrent_plot_backend\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mtr\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mStr\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/bmcs_utils/app_window.py\u001b[0m in \u001b[0;36mupdate_plot\u001b[0;34m(self, model)\u001b[0m\n\u001b[1;32m    250\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0mupdate_plot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mmodel\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    251\u001b[0m         \u001b[0mpb\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot_backend_table\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mcurrent_plot_backend\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 252\u001b[0;31m         \u001b[0mpb\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mupdate_plot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mmodel\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    253\u001b[0m         \u001b[0mpb\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mshow_fig\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    254\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/bmcs_utils/app_window.py\u001b[0m in \u001b[0;36mupdate_plot\u001b[0;34m(self, model)\u001b[0m\n\u001b[1;32m     44\u001b[0m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot_fig\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mclf\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     45\u001b[0m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0maxes\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mmodel\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0msubplots\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot_fig\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 46\u001b[0;31m         \u001b[0mmodel\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mupdate_plot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0maxes\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m     47\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     48\u001b[0m \u001b[0;32mclass\u001b[0m \u001b[0mK3DBackend\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mPlotBackend\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/bmcs_cross_section/pullout/pullout_sim.py\u001b[0m in \u001b[0;36mupdate_plot\u001b[0;34m(self, axes)\u001b[0m\n\u001b[1;32m    934\u001b[0m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mhistory\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot_geo\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0max_geo\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    935\u001b[0m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mhistory\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot_Pw\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0max_Pw\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 936\u001b[0;31m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mhistory\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot_G_t\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0max_energy\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    937\u001b[0m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mhistory\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot_dG_t\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0max_dG_t\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/bmcs_cross_section/pullout/pullout_sim.py\u001b[0m in \u001b[0;36mplot_G_t\u001b[0;34m(self, ax, label_U, label_W, color_U, color_W)\u001b[0m\n\u001b[1;32m    259\u001b[0m         \u001b[0max\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mt\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mW_t\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mcolor\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mcolor_W\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mlabel\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mlabel_W\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    260\u001b[0m         \u001b[0max\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mt\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mU_bar_t\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mcolor\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mcolor_U\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mlabel\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mlabel_U\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 261\u001b[0;31m         \u001b[0max\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mt\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mU_bar_t\u001b[0m\u001b[0;34m+\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mG_omega_t\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mcolor\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'black'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mlinestyle\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'dashed'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mlabel\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'G_t'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    262\u001b[0m         ax.fill_between(t, W_t, U_bar_t, facecolor='gray', alpha=0.3,\n\u001b[1;32m    263\u001b[0m                         label='G(t)')\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/bmcs_cross_section/pullout/pullout_sim.py\u001b[0m in \u001b[0;36m_get_G_omega_t\u001b[0;34m(self)\u001b[0m\n\u001b[1;32m    178\u001b[0m         \u001b[0mA\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mxmodel\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mA\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    179\u001b[0m         \u001b[0meps_tEms\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget_eps_tEms\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mslice\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;32mNone\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 180\u001b[0;31m         _, D0_rs = mats.get_corr_pred(np.array([[0,0,0]]), 0,\n\u001b[0m\u001b[1;32m    181\u001b[0m                            kappa_n=np.array([0]), omega_n=np.array([0]))\n\u001b[1;32m    182\u001b[0m         \u001b[0mw_ip\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mfets\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mip_weights\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;31mTypeError\u001b[0m: get_corr_pred() got an unexpected keyword argument 'kappa_n'"
+     ]
     }
    ],
    "source": [
@@ -214,22 +224,33 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 5,
+   "execution_count": 10,
    "metadata": {},
    "outputs": [
     {
-     "data": {
-      "application/vnd.jupyter.widget-view+json": {
-       "model_id": "ab690fc65a954e10ba5fb6ae24ee19ff",
-       "version_major": 2,
-       "version_minor": 0
-      },
-      "text/plain": [
-       "VBox(children=(HBox(children=(VBox(children=(Tree(layout=Layout(align_items='stretch', border='solid 1px black…"
-      ]
-     },
-     "metadata": {},
-     "output_type": "display_data"
+     "ename": "TypeError",
+     "evalue": "get_corr_pred() got an unexpected keyword argument 'kappa_n'",
+     "output_type": "error",
+     "traceback": [
+      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
+      "\u001b[0;31mTypeError\u001b[0m                                 Traceback (most recent call last)",
+      "\u001b[0;32m<ipython-input-10-8fbd37e44f78>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m\u001b[0m\n\u001b[1;32m     20\u001b[0m                               loading_range='non-symmetric')\n\u001b[1;32m     21\u001b[0m \u001b[0mpm\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mrun\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 22\u001b[0;31m \u001b[0mpm\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0minteract\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/bmcs_utils/model.py\u001b[0m in \u001b[0;36minteract\u001b[0;34m(self, **kw)\u001b[0m\n\u001b[1;32m     43\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     44\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0minteract\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m**\u001b[0m\u001b[0mkw\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 45\u001b[0;31m         \u001b[0;32mreturn\u001b[0m \u001b[0mAppWindow\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m**\u001b[0m\u001b[0mkw\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0minteract\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m     46\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     47\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0mapp\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m**\u001b[0m\u001b[0mkw\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/bmcs_utils/app_window.py\u001b[0m in \u001b[0;36minteract\u001b[0;34m(self)\u001b[0m\n\u001b[1;32m    133\u001b[0m                              layout=ipw.Layout(align_items=\"stretch\",\n\u001b[1;32m    134\u001b[0m                                                width=\"100%\"))\n\u001b[0;32m--> 135\u001b[0;31m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmodel_tree\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mselected\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;32mTrue\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    136\u001b[0m         \u001b[0mdisplay\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mapp_print\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    137\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/traitlets/traitlets.py\u001b[0m in \u001b[0;36m__set__\u001b[0;34m(self, obj, value)\u001b[0m\n\u001b[1;32m    602\u001b[0m             \u001b[0;32mraise\u001b[0m \u001b[0mTraitError\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m'The \"%s\" trait is read-only.'\u001b[0m \u001b[0;34m%\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mname\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    603\u001b[0m         \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 604\u001b[0;31m             \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mset\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mobj\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mvalue\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    605\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    606\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0m_validate\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mobj\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mvalue\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/traitlets/traitlets.py\u001b[0m in \u001b[0;36mset\u001b[0;34m(self, obj, value)\u001b[0m\n\u001b[1;32m    591\u001b[0m             \u001b[0;31m# we explicitly compare silent to True just in case the equality\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    592\u001b[0m             \u001b[0;31m# comparison above returns something other than True/False\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 593\u001b[0;31m             \u001b[0mobj\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_notify_trait\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mname\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mold_value\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mnew_value\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    594\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    595\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0m__set__\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mobj\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mvalue\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/traitlets/traitlets.py\u001b[0m in \u001b[0;36m_notify_trait\u001b[0;34m(self, name, old_value, new_value)\u001b[0m\n\u001b[1;32m   1215\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1216\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0m_notify_trait\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mname\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mold_value\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mnew_value\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1217\u001b[0;31m         self.notify_change(Bunch(\n\u001b[0m\u001b[1;32m   1218\u001b[0m             \u001b[0mname\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mname\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1219\u001b[0m             \u001b[0mold\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mold_value\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/ipywidgets/widgets/widget.py\u001b[0m in \u001b[0;36mnotify_change\u001b[0;34m(self, change)\u001b[0m\n\u001b[1;32m    604\u001b[0m                 \u001b[0;31m# Send new state to front-end\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    605\u001b[0m                 \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0msend_state\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mkey\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mname\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 606\u001b[0;31m         \u001b[0msuper\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mWidget\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mnotify_change\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mchange\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    607\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    608\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0m__repr__\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/traitlets/traitlets.py\u001b[0m in \u001b[0;36mnotify_change\u001b[0;34m(self, change)\u001b[0m\n\u001b[1;32m   1225\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0mnotify_change\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mchange\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1226\u001b[0m         \u001b[0;34m\"\"\"Notify observers of a change event\"\"\"\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1227\u001b[0;31m         \u001b[0;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_notify_observers\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mchange\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m   1228\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1229\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0m_notify_observers\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mevent\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/traitlets/traitlets.py\u001b[0m in \u001b[0;36m_notify_observers\u001b[0;34m(self, event)\u001b[0m\n\u001b[1;32m   1262\u001b[0m                 \u001b[0mc\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mgetattr\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mc\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mname\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1263\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1264\u001b[0;31m             \u001b[0mc\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mevent\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m   1265\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1266\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0m_add_notifiers\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mhandler\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mname\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mtype\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/bmcs_utils/app_window.py\u001b[0m in \u001b[0;36mselect_node\u001b[0;34m(self, event)\u001b[0m\n\u001b[1;32m    227\u001b[0m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mset_plot_backend\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mbackend\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    228\u001b[0m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0msetup_plot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mcontroller\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmodel\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 229\u001b[0;31m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mupdate_plot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mcontroller\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmodel\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    230\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    231\u001b[0m     \u001b[0mcurrent_plot_backend\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mtr\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mStr\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/bmcs_utils/app_window.py\u001b[0m in \u001b[0;36mupdate_plot\u001b[0;34m(self, model)\u001b[0m\n\u001b[1;32m    250\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0mupdate_plot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mmodel\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    251\u001b[0m         \u001b[0mpb\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot_backend_table\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mcurrent_plot_backend\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 252\u001b[0;31m         \u001b[0mpb\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mupdate_plot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mmodel\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    253\u001b[0m         \u001b[0mpb\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mshow_fig\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    254\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/bmcs_utils/app_window.py\u001b[0m in \u001b[0;36mupdate_plot\u001b[0;34m(self, model)\u001b[0m\n\u001b[1;32m     44\u001b[0m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot_fig\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mclf\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     45\u001b[0m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0maxes\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mmodel\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0msubplots\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot_fig\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 46\u001b[0;31m         \u001b[0mmodel\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mupdate_plot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0maxes\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m     47\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     48\u001b[0m \u001b[0;32mclass\u001b[0m \u001b[0mK3DBackend\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mPlotBackend\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/bmcs_cross_section/pullout/pullout_sim.py\u001b[0m in \u001b[0;36mupdate_plot\u001b[0;34m(self, axes)\u001b[0m\n\u001b[1;32m    934\u001b[0m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mhistory\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot_geo\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0max_geo\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    935\u001b[0m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mhistory\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot_Pw\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0max_Pw\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 936\u001b[0;31m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mhistory\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot_G_t\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0max_energy\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    937\u001b[0m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mhistory\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot_dG_t\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0max_dG_t\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/bmcs_cross_section/pullout/pullout_sim.py\u001b[0m in \u001b[0;36mplot_G_t\u001b[0;34m(self, ax, label_U, label_W, color_U, color_W)\u001b[0m\n\u001b[1;32m    259\u001b[0m         \u001b[0max\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mt\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mW_t\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mcolor\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mcolor_W\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mlabel\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mlabel_W\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    260\u001b[0m         \u001b[0max\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mt\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mU_bar_t\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mcolor\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mcolor_U\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mlabel\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mlabel_U\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 261\u001b[0;31m         \u001b[0max\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mt\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mU_bar_t\u001b[0m\u001b[0;34m+\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mG_omega_t\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mcolor\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'black'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mlinestyle\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'dashed'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mlabel\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'G_t'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    262\u001b[0m         ax.fill_between(t, W_t, U_bar_t, facecolor='gray', alpha=0.3,\n\u001b[1;32m    263\u001b[0m                         label='G(t)')\n",
+      "\u001b[0;32m/opt/conda/lib/python3.8/site-packages/bmcs_cross_section/pullout/pullout_sim.py\u001b[0m in \u001b[0;36m_get_G_omega_t\u001b[0;34m(self)\u001b[0m\n\u001b[1;32m    178\u001b[0m         \u001b[0mA\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mxmodel\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mA\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    179\u001b[0m         \u001b[0meps_tEms\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget_eps_tEms\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mslice\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;32mNone\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 180\u001b[0;31m         _, D0_rs = mats.get_corr_pred(np.array([[0,0,0]]), 0,\n\u001b[0m\u001b[1;32m    181\u001b[0m                            kappa_n=np.array([0]), omega_n=np.array([0]))\n\u001b[1;32m    182\u001b[0m         \u001b[0mw_ip\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mfets\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mip_weights\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+      "\u001b[0;31mTypeError\u001b[0m: get_corr_pred() got an unexpected keyword argument 'kappa_n'"
+     ]
     }
    ],
    "source": [
@@ -374,7 +395,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.8.8"
+   "version": "3.8.10"
   },
   "toc": {
    "base_numbering": 1,
diff --git a/tour4_plastic_bond/4_2_BS_EP_SH_I_A.ipynb b/tour4_plastic_bond/4_2_BS_EP_SH_I_A.ipynb
index aec6732..409829b 100644
--- a/tour4_plastic_bond/4_2_BS_EP_SH_I_A.ipynb
+++ b/tour4_plastic_bond/4_2_BS_EP_SH_I_A.ipynb
@@ -67,13 +67,13 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 1,
+   "execution_count": 9,
    "metadata": {},
    "outputs": [
     {
      "data": {
       "application/vnd.jupyter.widget-view+json": {
-       "model_id": "aa2d73eaf48e4676b6688232447794fc",
+       "model_id": "0b02442c72fb4a6f98046ba52f7f69e0",
        "version_major": 2,
        "version_minor": 0
       },
@@ -89,7 +89,7 @@
     "%matplotlib widget\n",
     "from plastic_app.bs_model_explorer import BSModelExplorer\n",
     "bs = BSModelExplorer(name='plasticity explorer')\n",
-    "bs.bs_model.trait_set(E=100, K=3, gamma=9, tau_bar=20)\n",
+    "bs.bs_model.trait_set(E=100, K=0, gamma=0, tau_bar=20)\n",
     "bs.n_steps=200\n",
     "bs.interact()"
    ]
@@ -104,13 +104,13 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 2,
+   "execution_count": 10,
    "metadata": {},
    "outputs": [
     {
      "data": {
       "application/vnd.jupyter.widget-view+json": {
-       "model_id": "00feafdb5e594b2391ecc0b781276dc9",
+       "model_id": "2236eef753d545c1a37952c3883db077",
        "version_major": 2,
        "version_minor": 0
       },
@@ -240,7 +240,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 3,
+   "execution_count": 12,
    "metadata": {
     "slideshow": {
      "slide_type": "skip"
@@ -316,7 +316,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 4,
+   "execution_count": 13,
    "metadata": {
     "slideshow": {
      "slide_type": "fragment"
@@ -325,7 +325,7 @@
    "outputs": [
     {
      "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAGYAAAAXCAYAAAD5oToGAAAACXBIWXMAAA7EAAAOxAGVKw4bAAADfElEQVRoBe2Y603cQBCAj4gC8ugAOghJB6QDUCoIdAC6X/APkQ5ICREdQDoAOsiVgOiAfJ9v7djGj7Vv7xQRj7S3u+OZ2fG8dnxbz8/PswnGW+D8/Pwj3PfjJTRzbjejJ+wAC1xCu4eDHgbw9JK+6aWYCFotELJlltopHjg5ptXsUQ/mUJkxyWEqZSNNSpbswLrDfDtGBHxv4dOxgrKEb+CfXEyO0Qrj4BQ2x1i4xAnHOTPrK9Y2EbviplKmFQZCiPZPzKOyJRx3BP9+6WhLohlolzc5pmSYIUuNeDGEoYHWbLlrwGeordfyHROiz4j70fayKfDI9264Z85KTgqZykCezj7I5b6mUqbBHOuG5J0YzrB8HTD2cuWnyz+3BDMG0jhPzF13RxHVJdbRS86yIzNb/EjNOjKFJXEMAhV8xGiL2GNo1lpifJlVAP1O4DcbFowicssyA43dUyMMtQP0OuWU+YsCw955sbJjEKKievqQ8cjQSZU2Epqkf1cgPymgn6XEwPnAOHHforMB1ni3gB9kB+h1ijw6JuvEWNsQZLZbyTFB4A3zNQIzYD1jVBzBPs8oaS7Yf19SL3/Z59H6E4z9vVG7MeC8TF9mO61cFwOtAJ5ZERqzhWcattcOhbDlwm8WK0zlD1Bk6ZzZ1tnZmQ9/MdrKkHR1OERAxfgSBAXnzJWXCs90ji/9judmWAHsxd8yv5BZEIUFNBpnv45n/z7gHhuePcD3QqcGulmQrxN2WRcBwvo3uMo90MQvDlod1WiHNp46fhshGqmxptaJI/ZfoblpoStHY1HqOF8jXDMXRmjhz9DQZRFVpwHvxW27XMnGOl3EPs9udczOCrLVsRJQHbK67NDB9vdR6nZZ4zR+NIWXso6bHRmA0yl3zFFOCWxrnYIudmV+medVZM7ewIqFVjvECljpjikfwkt4mRmxXeUoi0ZodI7Rt+ihLx+xybV6Wi51jkFj8ERlC3Qxduh9l5QZY5R0Rn54SaPRCNQpXd8LvcqviyDoZYCpp0NHxUKvHWIEpXTMZw68ijg0j8B/0ikl/S1dljIDqDPgSjwuY+1QY6tuU5ayqK6H4y0RMQ6sarrhHc7wsjdrikYlRgV4Yu3QKS6ZYzpPCQ9R2gi0Bq8jW7wDou6BoE7vhL6putXes+oEG/t3mZe0AzOazBhb2itwQ0oELP8P/AFbrDhX0imFlwAAAABJRU5ErkJggg==\n",
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAGoAAAAhCAYAAAAxvj44AAADaklEQVR4nO3aSYgdRRzH8c+EuIAaRYMoIuhBCCLjqDm4oERzmIOMy0URXHIQgyISMAiCwiguOC4MCgEPYuICHp0oIi4oohAxhhFcg4FBMQdHFKOiMWPi4d9j2pfX25tuO2P6C0XR1f2v+r3+dS1d/YbGx8d1HPwsaVtARzmWti1gkXEcfmqj4a5HVeNKXIOh/zp1RlVjFK+30XBnVHkOT9IvbTTeGVWey/BuW413RpVnDJvbarwzqjyn4pua6robH2EXZvEKzsoL6IwqxzmYrrG+VdiAC8WQOoe3cHxWQPceVY4rxFNfF6M9xzfgZ1yU1U7Xo8pxHrY2WP8xwosfsy7ojCrmFOzEvgbbmBRD65asC7qhbz9LsLdP+Zh6h71eHsUluBh/ZV30f+lRa8QTv2rA+DF8of/9WI23c2KXC4P3FaTdOLIn9nHcmLTxdZ7AOnvUcnwv9qby+BPH4o8a214IK4XBZ4hV2Pupc0eJpzxP69G4P3V8Gm7Cx3g1VT7bU8+TuBaX4vMikXUaNajgttmapPPF6i5t1CjeKIifwXjq+Gbxu1/CYxkxG3A9rhILiJOS8l+TdAB1GjWjuuCDiZeF5rtSZZeLl9MqjCT5dM41tyZ575B6n3/fw39ocjExkuTTDbZRJ1OYwAp8KearE8RwXoWzk/yTnGuKpocD6J08ZxRPiun0woCCl2KHGAJ6mcB3YsumH/00Ppuce6fPuY05GtNsx1di+IML8GHJ2HmGMCyW87MVY3Pp7VE7VJs/dmaUFwmew8N4SnT1+af2dqwVy9VvM+qeFF9a04yIj3qbhJFppjPq6cdUUs+EMOy5CrFwOpbhg4pxhfQatbqmessI3oR7cQfuwdXinWJM/rAx2adsjbjBGy3sU8QU1uNEnInPKsavSPJPF6ChL029R5URvAeP4DaxunoRt4jNybbYgh+wTsF7TQbLknxXXYLmacqosoKfwe94DQ/g+Yb0lGWv2IVYb7DdiO1Jvg5P4Lp6ZDVnVFnBu8USdRseakhLVTbjN7w3QOw2MYzvEb99ZV2imjKqiuBh1VdXTfKmWEzMDRj/IE4W9/bOukQNtfxP2SPEn0XW2r/E7uhD25uywzhM9MCOHNo26lwxT1VdBh9ytG3U02Lrf9D54JChbaM6StIZtUjojFok/A1nysHf0522WAAAAABJRU5ErkJggg==\n",
       "text/latex": [
        "$\\displaystyle - \\tau_{Y} + \\sqrt{\\tau^{2}}$"
       ],
@@ -335,7 +335,7 @@
        "-τ_Y + ╲╱  \\tau  "
       ]
      },
-     "execution_count": 4,
+     "execution_count": 13,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -373,7 +373,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 5,
+   "execution_count": 14,
    "metadata": {
     "slideshow": {
      "slide_type": "fragment"
@@ -382,7 +382,7 @@
    "outputs": [
     {
      "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAGUAAAAWCAYAAADZylKgAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAER0lEQVRoBeWZ61EbMRDHjYcCPCnB6QBCBTEdAKmA0AF8hG8MdBCoIAkdQCrg0UHoIAwdOP+fIt3s+fSwfWcuwM4okndXq9U+dWQwnU4HfY/j4+Nx3zr0dX7s7sNBz3BycnIoFTZ6VqPP48feBpUOa0SIBc+wLdwnj7/T/Gx4RlqPNT5o3Ikf3qVAe3e0cUvz0VIC3sgmb/NnzRfuSqm0VVpNNe4z9EPRr1P0El57Rzn5pf1vjY4tsAn3ipYveWzig/AmE4x49TFDL5HOxPCtxPSO6NgCm8SdInwoSdcwZeB3hlYi7VXpWuJ8B3RvC2wyimaKbOAyRQy1TNFvmrIDrekzNbonFSftpZe0ybLiGa+UAZvsrSeU5zVUM7gMiaNswx8I95DYX0KTiTX5sQ2STzr/8bSPmo+Eq+kQ27cK3Avpgk22G07R4S5LROSp9tNfkBcXeAzTBfCyy/YTnU3pxAnO8V4vnHTQhQKLyHhBXWgHk4ZThAz9ZNcYhCfwtX43So5wZNWlxpPWYa9+ZgEnP6U4JIfzJpptJvJsDkGS2to5vitdJIc7Y6cdrdcSimKTccwpLlOsQbR+1LiygvR7pMHb+kEDZzUcZvln1nzjFMuQ5JItOOJG66LDxcPFf2kwzwtV8OU2LKrLrCztx1bfhSfgUoANRzGnNPqJl1CVG3/5r8KfexqO/OzXrSfJJwiQzcMiBMm5cNmPTNFx9GZrBYyAZXUxIuyy1EtdsNacIgWcASSFCK0ByhkEDnFfn8LjRCKTcsP+LY19rTFQCkjTbDRrPw6gpyCT9aHWpwW5YuseSrqIjj0IhntzOg5AX1uCuUuuJ2KTp6ERwjKUiOTLSIeQfvxpJBidg3DYhXBE962G+wjSnAL4o2ksGWca1feP1nOVrtRBbfDz6CKeEFyUWaoJ+hKwpxpVD/R8PJ6SthU/mfI46xS+HwbaaL0LyoHwOIBoqEqZ1jgSQwYnMfO6ygHyyagYcIaVjz7ohdPDGbF9q8DNowuG/KFBxbiSjqGioKsNPGRF7Sp8ALLtYV1CRlrwKkCAEyJc5WHhAHgwNDON3XobvK31CA6KaRkFGt7sGYFxXwtKISWB8wCaZC7t/3F1/29RF+nl7qqZ4LRlHyc9G5WgW7sZUrXEcQc4hY27FXrxBa8wGwF7EsFlkgC/xkBjg7Vl9L9rOEt/yfWCuswG5xfpSgYFcAaXTOxlneXowoWEuBmGHS1mXkouojUT3fw5v/Z8Tsju5UMwoUsrtDcoNnDG9r/JFFtBQj/BRjGA1/XiLpxCltFTeL4SBaRpEcRHM0RRFyHFDf83Q+gXfBhiBwy8qbXNCHoitEYV8DbAFu5Fu972rhLEIcvWexxKb5nLkW11XeF++iivLl6fURAtZyMeNhW9i0yJKjEPUooSSQeaiaDXDGTK7TIX8HfHBtXjqPHfwcsIfs97ZEz6AEFF6enkr9h/AXhx+5axk4GHAAAAAElFTkSuQmCC\n",
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAGEAAAAaCAYAAACn4zKhAAAEOUlEQVR4nO3ZW4hVZRQH8N+xm5ORlUpRIYqmRqVndLIL1YtWFhRFZARlZFBQYA8JUQkOkpVERQ8SRPWg3aQMzd4sohuUlZklhd0vZJkYCZlJaQ9rb2bPnn3mXObMcZTzh813zlr7W2t9e33r8u1d6u7u1saBxbAW61uB7RjRYr1DATOwH7fkGUVOWJ/c3N+1qAEjunADHsJfDcw/2PEx1uB+HJNlHF5w83T8i6X9CHylASMewC480cDcQwUP4gMsEM8DfZ0wASfgE3Q3UfkkzMZT+LuJcg82bMCXuA3L8B9901FXMn7YZOXzUcKqCvwLRXR9gz2ibmyQ2S0txGDb8iLGik2J1jlhtvD6+wW8e/G2KFxv4FGsw1GY02Q7qqEVtryXjBenhHw6Sp1wPk6tIORh9RXWESjji4J5J2IJ3sUs7M3xR9ehZ6BolS3pBr8oJWSdUEJn8vvmCgJ26l0rlmKmjFcLcAoOw7YC3pSEt1XfRcOOfuQ2G4NhyzJMxWUZ2p8izY1NCdl0NAkjxU4oVbhG5ZSU8WkVQ9I5fxTwtiRGzceruA7HV5GXxfeqt9PZ69l+ZA3UliKUsamAvlMmsrKRkKaiag81r+SFKvek3dDwAt4OXIDFuBxXiNqxHvdhYxXZafGsFb/0wxuoLUWYhmcK6B0yXeJAnDAGJ2OfKGLn4SvRfmUL8PZkzEdRis9xLY4UefLW5P/ZIpX9048Ns2q0tVbUasto/I67Ev4Z+BV3CKfBSaLObMrpGIbj8F2WkKJeJ6T1Y6EoaNPFTlult3O3JQZPriJvL17HXJESR4lFHAhUsyVd++3i7UEXNuM5PRHfid1iY2YxWaT2TSlhWGYsi139WY2GlhNjr8Zb4hBytyg44zL37Rdt32hMzC1kQoHciTgTP+LnGm0ZKOq1pSzeKswRWWCrSFljcFrmns3imWZxbjK+mRLSHXu6eJ+xSzzISnhST17txMv4IcNPW9D8+WM1rsGl+DqhLcBN4iC0RaSt8bgy4c8vWMBgoV5bylirZy307ajKiovyJaLWrE0JqRNmJOOxojAVYR8eySl5OnfPdNFhfJujr8ZvmIflCW1ton+myKvDhYOfF61dPowHE/XaUta3IenSO/2UxWEvi5G4Cq/hp5SYOmFFctWKo0XYHZGhlXAnVopQzWIvHhdH/07xbmpNcg0FrFG7LcNFXs+fsRYIp+0Rz2eivpEwL5mf3cwNf0+YKnL9jTgnUbhSnLIrRdJjIrcuaVDnUMFZyXi9aGmniLWPF3WBaE3pXV87cI/ICu9kBTbqhLLo0RfhJVGAOkTR2Vlhzh7htI8c3B91yiLlLBYv4zaK9czU045PE8V6d2beOFFTF+YFltqfN+vGctEFzW2WwFZ/3jwUUBaR3zS0nVAfSqImNNUJRZ8326iM/aKNbyrakTAE0HbCEEDbCUMA/wPPCgq3T3pJlwAAAABJRU5ErkJggg==\n",
       "text/latex": [
        "$\\displaystyle E_{b} \\left(s - s_{pl}\\right)$"
       ],
@@ -390,7 +390,7 @@
        "E_b⋅(s - sₚₗ)"
       ]
      },
-     "execution_count": 5,
+     "execution_count": 14,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -475,7 +475,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 6,
+   "execution_count": 15,
    "metadata": {
     "slideshow": {
      "slide_type": "fragment"
@@ -484,7 +484,7 @@
    "outputs": [
     {
      "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAAADMAAAAvCAYAAABOtfLKAAAACXBIWXMAAA7EAAAOxAGVKw4bAAADlklEQVRoBe2Y0XHbMAyG5V4HyF03cDdw2g3cDdJ6giQbJOcn+y2XbtCskGzgdoO0GzQjpNnA/T+F5ImWrFCkyCfhjgZJASCBH4Rkzvb7fVWSttvtQuut1K7GXvf92AYD7N1K5lpOXQfIDhJ5N0g6UdigUon/STTVqV4ambV2ATKDSQE4kRL60PyVVeeafzH9qpgzWpQNzMV/2sUH8lvpXlod9X+o/1vto53z0kwCJ2r3NCswIueMpJyTC+1r2dgPCBMcCkpNh84A2bnamYSi0uHVrP8rW6TIJ/FYVDAIKo90jpHnDEJaEIce1C4Yj0QE5ibFlvZ1Z/ZmzeDck+ZcMWk5YyTJR1LuzGrGctkAlaU4ARqFZIvUYm+nTYOdzkiYdAAhd+CaSgP70RWsax3tba55kD5Vnz066qtmd5K6kgIIeUpWW/NE50W87yxw/lzFsbqW6xkbI6VBsIsuJcNeKnEc4YX7pTFm/olxnzOkGp8cLPRdzSMZ4BlRx5AHtxU0MtjpJD3nGYH6qvashmNexZNMfSbEcQR5nLEVjMxx8n3OYBxCwXPGGCNaH9RAb6HmDqLmLBHVTlSMjZ24O0vqV2pddrDHOwX04I4k747CrOtDUwIo8a7Zqdn8bC1i5P5J5kF9outI4zp1xL1AOIGDjuSI9vrQzoFY77BVAGQMOH+psbk6V8Wd9+o7kiwpggznAr0mAb/Vb84f66/0gOBFk+eMNkR0MFh/85jNkgbfelYAOcjlrvQoDKCFs6GEzmOocJecc8Y4AiLkeTOlOHRH3zmSpQBQzfjcID2htdpN3Qv4kR6o8mnSXDdA0xepnZGRpaY5WCDilVkzZsOdqWbMWXRwqI6w+FBU6vJq7EUxiwwpAiKushxY4zlvcSLYIs0TAKIKIjTrnLpB9FlSZEASdVazGItyCESogK3KFmMvRsciE6Pr6RhUQccVAk+gwGA0ZArs9c0lRkPmzZUKCEzOFAhy1BITMlFhK6A022w2Ze9nMzo1leaMwU0yPRWApPBlVJ6QyRjcJNMTMknhy6g8IZMxuEmm+240kwzrnyf3AEF3yEkLNZSzfM7IES4nuJ3hDu5ZDce8v9OSSbpWkr02cT07ZtOH60LtrGlT411znKs/epqZiLuoa8wt6ZA7tHbEA2dKVLOV9pJ0hxzoS1XCmeQ75FBnshQAu7hSbK7+X/GZncvJcyMDKsl3yKEByO3MKHfIoc78ByMwvF+EqgZXAAAAAElFTkSuQmCC\n",
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAEsAAAAgCAYAAABafYVzAAADvElEQVR4nO3Za4hVVRQH8N9MDwzToqQHPTChqOhhZGAPJIs0CB0nqAh6QEQSQQ70AKFgKnrQGz/4pfpSEX2qpjdZERFlZDFBFkqg9MEoe5qSlVkf1pk8Xs/d59w75x5pmD8c9j37rL3WOv+119p739M3PDxsEtXQv7cd+D9h34bsHIyfG7LVir66FDU1swZwuXC86as2NEXWQrzZkK2eoQmy9s+u3xqw1VM0QdYFeK8BOz1HE2QtwssN2Ok5miDrGHxTk67l+ARbsBmv4JSadJei12SdgdEa9Z2PlThHpPcOvI1DarTRFr3eZy0W0a8LC1vur8avOLdmO4Uomlk34B/8gMfayFTFmVgzjvFlmCb8+6mHNv5DEREb8DD2wxAu6lL3UdgkiO8VHhdpvrqC7OH4Gyu6NVZE1irchgey+7ld6CBWwV6mxkOYh8sECWUYEL6+2K3BVIqNRev0hMwifNVGz4V4JzF2BnaKmZe6/sCUlrGP4JrMxtcJG3kM4ke8X1F+D6QK/IasPa3N8zlidTperE4f5J5NFdHentB/IO7O3c/EtfgUr+b6N7foWYErMB9fJvTncZBYPZ+zaxbOwPfKz49/ZuO3p8i6M2tniRfb2vJ8TXbNFatenqyFeKvEiY0Yzt1fL8h6XtTMIqzEVVgiivoRWf/WAv/yuEQcuV7I9XUcrHZkLcB12e8+nIqP2si+JF709hbnliecL8LsrB1NyNyYta3pfZfdiW/FILaJejyGjToMVlGtmYYn8AuezfrapSKM4AScmNN5qJjinWCsNn6ekGn3N8xwYswUXIw3pMvC7KwdbSdQRNaDOBbL8HrWlyry67FOpCKcjY8T8kXoEwHZJKZ9nVggUq5sFSwNVitZ87FU5OzTuYGpmUXMroHs92KdH5yPw/SUo+PAoCjSryVkKgUrT9ZUPCmOD0uzvnX4XdSs1KoxIgr9YTgZa5Pu74mxFP6iw3Fl2Edsb94V79UOlYKVJ+t+sfLdLBgmltm1maKZCT2rxfFoSPV9Tx7Ts3ZLF2NTmCfqZ1kKVgrWGFnn4Sax436mRWY0a1OpuDMbe6vudu3rs3YIj+LKLnQU4VLh20iJXKVg9eMAPGX39MtjbGqmijxRp7bpbof8Ge7AX4KwOV3oKMISfIjvSuQqBasf94ilfxm+LZAZzdqyIr9KrKQ7SuTa4V4cmfl0S5c68jgLR6t2FqwUrL4J/EX6PrExnmXX0W1cmMhfpAdFCamFKJr7Ir03cFLdCifyzKodk2R1gEmyOsC/pUK+aj/Gp7MAAAAASUVORK5CYII=\n",
       "text/latex": [
        "$\\displaystyle \\frac{\\lambda \\sqrt{\\tau^{2}}}{\\tau}$"
       ],
@@ -496,7 +496,7 @@
        "       \\tau       "
       ]
      },
-     "execution_count": 6,
+     "execution_count": 15,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -606,7 +606,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 7,
+   "execution_count": 16,
    "metadata": {
     "slideshow": {
      "slide_type": "fragment"
@@ -615,7 +615,7 @@
    "outputs": [
     {
      "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAGcAAAAaCAYAAACq/ULmAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAEtElEQVRoBe2Z23EVMQxAN5kUwKOD0AGPCoAOAh0AHcDwl1/oAKgASAeECnh0AB0At4Nwjsfe2Yd3vdlsNrkEzTiyZVuSJVnW3uycnJxU/2EdCxweHl5D0kfaA/o7Jam7pQVrzaPsU9qDteRdhBzOt0Hus65s6G+6NMeXwjkod4Aud8DHOSX/MdrvzHlecfZPXfrORac1lPKqfwbf6Sq3TWP030dfb4D4Fc1z3Yz4BfPemgos/Q+4ldYYP3Ut+LXrhMtwc8zB2WsdNNySPxj1J6q+oOmcY41Mc/yD9pk2Cqx9y4JnYJ0XYC91EmbyOf2HtLuR9hUcvB7HblaBG7SvrHftLGCvfO5O4cEaD+mhezl7lvAzbCrosmFeRwWgr5NMWxYBpbRtkL6jPXJzzjleKxlaxn0HZ40PPTlRPnPByDJiRgFZKSBSwIyuP8/JmboY3AZiCbSFjjS9bXZzq5lIVdOYp2VUR0iOzwSaebb3EHb3qSi06+ALf5dm6mJwFW0VeX9n7WNtkHUO9HRbSoYz1cwCFLkdN5o2ixAVL65bY0FBF6NeZwSgb4YxHY8FelwdkPYIaW3IOeHmdBlGQYEDfaN5qsCwp/NHGT8jn87U1g8fc66DaK9b4BDsYJ1mJVfRHyqCDPiQvntvjhsBo7pleJhpTB1SAzSv4Fy4x8YWvyFGyPFAv+L8LXBdmg7tOS/6BF0sCLLvKHTPa0EzVtRo03D7es6BgU4Q9ulb5gp6XLqGWQrkmfsga/FHB1OrzgiBEPXTWWMHbPFYajBRF891Fkg22c+ltfTePEKZ0JCkIUxBvUcN2m3aN1rpfeoqbCluJA0CPPeZtAQNjokLrfC+DW46p4mSLnH+peLph9Q1U5Vkkxu9mwPDcHMQUBuEvo45agpjHMo91znPXM9xzfWZ/qSb4z7463hvsQ9rCh6nesC8fP3oE08Fg7A+79imIV2ge/7wkI/tnzCXbk6Vc07vvYkM6wcsGsAyOP3UoEPvTxDcXGKEeHsGwQPTlGHFk4ImfXln97FevouX3HN0ySpYJiab/G6lNRQIBmB/L0WpXIOvjgmPHnSdaZSafp7TPtJOE7UNtv0uvHxv/B3KG2ORoozF+PclDlNW0iWdrfcRmlJGq1JrqouCvgP3wJtI16E67i00o/wLbUrOdU9ShG4b4OWXcv0dRb+Y0toclhutrEv+5nAcf7qvUCabf6HrCB/jOsXR16EaMjlLPOVnFp04tk5ZTTnqpX4GQZLFcBVYUxeD33J8s8efawze0STaKmiphHYouEZDii0AmjdLuhVUAvN9MwUmehcbAKGeh9+mO8n4Cc1UGX5Kj/MqvXoJvbIu2i/YV+domLNUGRq4edP8XUjDjgJ7TFPK1rlNZ4d9kWeT7yi/85xcWRftETLG7gKHsqLyRlVgo9x/I7TK7hEZH5hL79zIsqsxFe1ogaVdqr0Fju2tqx9vBJzG2EaI3yTNtLiASlvLwqxzhA1Dmr8M/6b2fXt/itu2tZYvKY4NrE4fgsObvURaK8kszfs+vSwt+tfncYhPwpvkGM974TdHJVDoAOS305VMb5zbKlnHtJ6Ey3BzKpSygPgB9nviKoJB2auY/wL14PT1wFRiCwAAAABJRU5ErkJggg==\n",
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAGEAAAAbCAYAAABsv+EEAAAEJklEQVR4nO3ZfayXYxjA8c/JW8kKZTZaO1aKDTu9IIax8jLDmOEPhIz+scywhjZnUZiIP9CMreV9ZIXNHwkLUxkxWub9bTVJTRvSUv647qfzePq9n9Pv/Oh8t3vP73dd98v1PNd9X9d9P09bZ2enPnqXfr1tQAlGpNIKdMeW+Xitlop7NzjA7mRpurb3phGJptjSiithj6PZTliA9RhYoU671lgFlLflbczDw9iUyhzsVaGvcdiBa4uKUk5YkipXKjNquYMC43EF7sXvDbRvNS4Xz+8kTMV1uKlC/Q+xCHfjgLyiVE4Yi22YVaHDl2u3dSezsRmPNdC2FVmHaWJSfo5RwglzKrS5BytSu9mZsOiEETgYq9DZY+aGgZPwBP7swX57k+XCARnv4y4MEpOtFCuFw6biPvzNruFofLp+0FOWJqagDS+U0Z8qVtfX2CLyxkq52dJEdrctz2O4mJRonhMmCa8vL6G7HctE4lqKB/Eq9sM5PWxHNeqx5UQxsTImYK3yqyDjvXQ9MxMUw1HmhJMxrEwn96svsQ5EB9aUaHcoZuJdTMTWgn5oHeN0l3ptOQwP4VEci1tF0q1GNsFPywR5J7RhTPp9TZkONqo/Vxwutm7rSuiOSrov7HrTsKHOsbpDvbY8k+qvELnhScytYZzfRJgbngny4WgUBouZ0FamDCl0OEtsaSuRtdlUQrc6GTUFr+AyHFT9Pnbynerb6Xx5ukJf9dqyDTfgwFTvZinRJq7GeWXabpRbWfmVkIWiTyoMXKSjhvrZbqh/Cd0GnII7cS7OFzeyBHfgoyp9Z8mzVtZW0HXXlnoYILdL7AknPFelzvp0La6ijM9wCfYVcfL69P94Ecr+qtD3xFoNrZHu2FIr/cTq+TYvyKjXCYeI5LRd7CT+SG0nFOqtwy8YXaW/rXgDl4qQOEQky96gmi2ni1DUCKNFaP84E/TLXTvEA/20xs6yJH6L2FWMFcv9Bf9eYTvEtm8oRhbal3pNPBLH4Af8VKMt3aWZtmST9K1MkD2so8X7jM2YXqGDx3XF1Q4xYy7C90k2XayGdnyVa7cQF+PsnHwarhIHodUibB2BC5J+ipgUzaCZtpwlcs3iTJA5YVy6DhKJqRTb8UDu/xi8pMsBdJ0DiofAhfgZk/FIki1O458g4m5/4eBnxZH+yxpuqKdoli2DcaH42PNjJsycsCCVeugQe+M8Y8U275uCfKt47TtbOG+VeKO4qM4xdxeLNMeWycLB+cnc8PeE/XEk9snJ2nAjnhJ76CJzRWyd2eCY/3UG4DYRFd7JKxr9vHmcSLhX4k38Kk7Sw3TF0SJbUv0zxKuM/8M3hXpoFzl1flHRqBM6xEFpBl4UO5/XRebfWKHdslT2RNYo88qnUSfMS4XGPvD0kaPvQ38L0OeEFqDPCS3AP0EBBUE8GnK3AAAAAElFTkSuQmCC\n",
       "text/latex": [
        "$\\displaystyle E_{b} \\left(\\dot{s} - \\dot{s}^\\mathrm{pl}\\right)$"
       ],
@@ -623,7 +623,7 @@
        "E_b⋅(\\dot{s} - \\dot{s}__\\mathrm{pl})"
       ]
      },
-     "execution_count": 7,
+     "execution_count": 16,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -649,7 +649,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 8,
+   "execution_count": 17,
    "metadata": {
     "slideshow": {
      "slide_type": "fragment"
@@ -658,7 +658,7 @@
    "outputs": [
     {
      "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAJUAAAA8CAYAAABvsfhpAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAJGUlEQVR4Ae2d73EUORDF15QDML4MfBmAiQBfBuaIwJABFJ/sby7IAMgAyIAjggMygAzORwa+95MllTSjGWl2vDvas7pq0IzUkp66n1qaP172rq+vVzXKxcXFW+H6rPRTjfjuEib54EDj/aLjsc5/5cZ+L6ewRLkl1GEj1BLW7/dpiXSmkm86h2CjUh2pBPqZEB8rfTKKvBVu1QLyx3d1+FrH+1zHVZFKwB9Y4I9zwFv59i0g/7xTrwdKT8d6r4pUAsoseCfQ2XV7bFCtbKMWYAV5Lx8NLoP7G+1+QuOW/USqRaKU+qfvbxMgF6mq3b0ixR1R0nh+6fgguASA5BalGlIJIOv1G0ArXULo/6H6Z+/QZNwC2OqHbHWk42dXtYrlT8CIEkc6LrsAt3Ft+18pbYQqMLjsBJH+0vE8pV5LpHoFSIFdKkrRP7Nvkggv+wrqIkwK5GzBcdwg2M6/H9UNNnvZ7a6KSCVQ3E0AcusiAkAGwjgzb6q8Vr2X9mB/caXj1vdlU0FtSR97Je8Ee5FKBnoh5T90HFtwX5WGEYTZiSMOdXyVPrpri+qz9CHFTlWdH+grTYZf01r5P8y03mwrrP5MGD7qcNiZueQ90LHIUqp+8Y95lqTz5Ea6cGyjamr7J4eU8H/01iNFqjdSYsPM+5vvSpOkUb4jn9RmiSGvBZhtSHqO1I702TpDCrYtHrSuS07qMemqEY2Fu7MzAfpXqYmkGwTHZOr5Ibn8CcyJBeJmYAoXD8J6O/+UYiYP0hbPaowm/ftKH2baLSkmslyWKKZ0hKH7TA2SMYOLx5Nqd26etRHRg7cTmxRWDLfS+H6SpFKpi06fvWb6hEbnCgSeRE5rtFn9qg0i3onSKHSv26jawbinOm6D7OvCCOvxQj655wmVZp4bv9mx+6aGSGUilZSjSKVrljwjOidiROW2aGqCcyeRamoHA/pr3fGl2pItjpRP1OM5F3ZZXIQD34Bl3aW9ZAzOb4zfS29PZUuYdRFhBBKiRQZT3m2F+eKIpz5x3j8W5+9Kufvq4iJisLeIxmDruORU5dTvie2DpeOgV3iT8Vw6LP8rpRgUDCa622vyncFvaizzLxhfCAsRK7IRcJRnbix0mh0n+glxYzwMy3qkUkcmSkmJ22x3m0+n5CedEDY45VztO4ZfldSTPssxDjRkVgomDONno/KIpkQhBpxciqwOy0NPVEY+DnCPCGg/ujuUjusf/OiDye0twBLp63opARv2YIJwA+Zlyjh9pc6J2mDikhuRcr+jx6XbTz1RhdB4fDDnmOmrKQ9jcgt7pXNX15dnTiIwY7pqGweyBzKYrC7Oc8RfqQwszM7fdDBDh27tiTS9CWLrRx8GKm+lI+xTTXvhmRRjiJ5NSd+T3Gsuc3JluwWPJ5XwYacp48yhx95eUqQ6oTQ0pM65o4k2tLo2IRU9ylWlRzjfy/DJoS0iMhSJ+iJaQSSeU0UkBguNKL1U4iJW9KxGZczaoShFfU8g6WL8QWwqv6/yyaJ6EPGLDtJS8ZO8pILtAzsx8Xi04CeYzieNs6C/aBwpUmHI1F7EO8ICDkMqRHxc0HlX5cpmRKC6SlyrT4jNbIMsjvg8T+stNcojLBOxnintvvREP7ksKr8rT5WRuwPu1sleg09KpRiy7XUVGLPyIJTzidsiDEXQWx3nvRCQwBhnKa9nSJWFkQhCuY0qJIQUJ9JhyeEJc5Yk0p8sape9C5+SEKEgvtmEDjSEIRFPOtVlA/9JKU4tEfSreriZA62x4Q/8Z95B2rGyyvw5UnfuOKMbrYhU6tQtJ6lIZTAJJLPgUeAYiAjheBBIJPlbh3OoTkfFOfdwTEvtEr49cJ33lr5ufemAiXEQrRzJX+mapTErqsM4iXJ+OcxWWlhBWCEUyyp7xhA3qwzbFcgTifJuY5zOj6btLqlMpx1AHoTyIRCbUr8U6hwi4nTXMOmxjhIpXf7oN+xzpf7A2n2i3e3TkRtioc+7Soezq9u9Rj+Mzt3yqq41LucbIlQUFOw1Y0ktf2uPU+1CSMT50Vzsq+BAZ+91GMaSqzzW41DQOdZBysY8BE2+X2J0zl6hyBlqh70PTu7diSkvlDNdsLyy7IIBoW7KSDel+lflRDRmLBEKTNGmXddj8kiFEZHHlCsowwdEqOiGKsBFOVuT7h5zzji9L4J+VpAKp04xdlifc8JqGGpZuyFBqeBsx/hkHdt+2EdSbyCT5Y5J4t6qD6jF2epzjk3ixrZwJbxu65LsTeWQrfdp88xxOr9F+857SQTTMnGWYaxSIglLzNBsSbUMIKLdRsRigZBhNN1IX3ewUaIc/icwedn3Z+ufMKP9RlodjM6YRDfs0cxmugsuobtWltrd2O37WoD+P5W4MeitILNJJYfR6OjeJmPDDypn78JGc0qEyzTbirdgAXzW8/1tLH+zsNvoxMZ/aoSb1W+rPM8C8huEQsKbNpOxOKkMipuN9NjDOavWkoosQBBwr+giWLWQiiWQu0jH/ghku6jSAtyUcWfdkypIZZdAnsa7h5U9oC2jHgvIX7x/XSlN7oGrIJUFyC0/D+a4o2hStwV4mHw2BLEaUlmAEIun+00qtYCNUjybSkYpYFdFKgF1Xz6wXjepzALyz5EgEaVG3zZURSprQ74B4hMXBtCkEgvIH7w1cZ/UjL7brY5UAs8jf25XN/ZdViV+2jUYvD/lzcngsucGVB2pACbg7ouCFq2cpxZM5Q+iFIQy25MclL1af504B7yV12uBKiNVveZqyEossHd+fl7nD6mXoG86VVqgLX9VumW3QbXlb7f9VyX6Rqoq3bLboBqpdtt/VaJvpKrSLbsNqpFqt/1XJfpGqirdstugGql2239Vom+kqtItuw2qkWq3/Vcl+v0qUVUCSm/l+WY+/P2GLjL/25/dgrt83V7TDHhfhOIPXPm2iw/TrnRAsOhP56XT++tc6TTh05d2xDbQS/YHOk5Du+j6c3jdzmObhfZoy18irtgI5KOQrvkLH6JWkwILtI16gZGk8lQHy2CTAgs0UhUYSSr82tzXMtWm1TbqGQ5o6TuSCv+1a+8HwzJV72xxi1R51xOlRv8kKd/E3dJopMr7m1+L4/FCk0IL/AcC2rvuV4tV4wAAAABJRU5ErkJggg==\n",
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAKEAAAAhCAYAAACx4ctuAAAG/UlEQVR4nO2caWwVVRiGn8smihQFFCOgFRTBBVoKBhAJS6RGRKgGiBE01C1GRUyIC4LWBQygYjAhxOUPGsXduoBaF1RUFIQaxSgKVDRVoaBUXECo/njPtMMwM3dmemduae6bTE7vPWfOe+7Md771pKmysjJyyCGbaJEl3jbAmVnizsEfvYDDkiTMlhA+DNTEzHEU8F+WrkMZvwL3JkmYDSGcBqwFfo6ZZxwwEUhl4TqUsQtYDdyaFGHSQtgVuAJ4LAGuYuCNBHiaI14AxqD3FTuSFsI70A+M22S1MdcfMfM0ZywDypIgSlIIOwOXA0sT4BoJrEyApznjKWAykB83UZJCOAkFI1UBxvY0V1SMBV5pxP05wG/ARuDSuImSFMILgHUBx75jrqjoDmxtxP123AasAWqB7cCrwBkZmrupYz16b7EiKSFMAUOAygS4CjPMMxxYjNY/EtgHvA10zCBHU0UlUAS0jZOkVZyT23A8kIdyUEGQ3wiuC5G2yhSKHZ+noDTG2RnmaYr4BWiNXKMNcZEkpQlPNO2uBLiKUB4yLrRHz21njBxBcDXKMtQAC4nnXdaa9oQY5q6H28IrSF8RmBWSJ8+0tb6jGo+uQDXxpoAeQmZqdYwcQbAFuB9pqunAuTFwWEqjvc+YLsB+YFFUEjdz3B/5PXN87nsxJE9r0+4NeZ8XWgB1Lt+PJV4TuQAYBpyDHnw2UWGuncBcYBDwZoY59pi2jc+Yceh9vBSVxCmEPZHDvZ7MJir/Mq3fjgqKsUgD9OFgQRyFfDYvdAa2kb60thfoAPxj++4BlDcbAXwfYr1xw9LI/WKY27Jgf/qMKQF2AB9EJXEK4QDTrok6oQe2mzbPd5Q0zE3ogXZF5rsKRaMzzfqGA6egaHWV7d52SDvZBceJI4G7bZ/zUQL9c+A1x3rt8yxCec4RwNdpfkPS2GLavh79jdl4HUy7/eDh9f0jUWLbsgyh+ZISws2m7eAzZiZyAbYis1KD/I0BwHmmf625BqEo2C6ExcBbadZRxYEa/kokhMuQdnXDYqQBxyPTd5z5fre5so3Zpu2BNplzTVE3HjQojU0e3GOQqba7Z6H5vIRwCNDNg3gB/urZDbtR9r2XR38XtPBVyKQ6fcfOjs8vIwG62fbdGJRYDoMC01b6jLnWtM7k+V0kVFv1wWig1PydQmc0P3GMqSL8xrNwKgr0vE48lSBZqGgMn10IUyjRCzDVg3Sng2AOcBbBIrP3UNDjht5ASySobsGL8+xhOTDf3PcNcow7ITMQBpYf9YXPmKZ6NKs98CjwO9Iwk5FJdgqhEwWmrQzAUYh3Db4tslAr8HeB0vLZUzS9kLlchfc5uU4uBH4v0I5n0U51i8g3oHRAKar5TgKO9plrI/AtMskAg4FPA67DQgq9tGq8fZ5MoYpwh2KfDDDnfJS/uxFYbr4LEpwE2XgWCtF7c8NoZHrTRcVp+ewCYZnioEIFEsKnA45diYSnGHjd0VcDDAXuBM5HEfB+pOZvx73mXI7SA/ORMIY9nXMS8nk+CnlfFGzCX1s4UZ2mfwRwDdKAS4HTzPdewYmFMBtvEPAvDQLuRAmyWs53GZqvMUJ4DCrH1SF/aTDwHXo4boncOuSzlXos/CtgAnJ0h6GKwARgIIqU9zjGlwMzgGPRSwhbVupt440bozI4Vzt0KHgXetYgq/A3sjQpvJP1YTbeVGAeEkQnWiJF8S7+VbBAfHZzHFYILf9xBgoq+iOJfwbvmvRypMK7+8y7F6VkJiLXoBMKXJxYjTTodKLl7ZKq4mQa96FIeBoNGnM/2oR5+Nfdg268jsDpwBKP/mHovaQzxYH4WtjaAqStvkwzsYUCJDAlwPsoQLgF+Sn5PvddRUNaASTMbmcHT0ZHprYCP7n016HqyAyiVUk2mnY68CBwSYQ5ksZQ4Dr0e59w9FWa1s8kB914cw3PPo/+i9DzL08zTyA+S2P1QRqqFgmSFx6hYfcVAs8DP9j6rdSNXzH9RzPPcOQnTkMh/GdoN29DatwKOkpxL9GBgpiLiZatX4dq4NcjQVxIcP82GzgceJwDzbAdlgXrh7dw2DdeR5QPdv7mIpRv9bOI44GPSX8qKghfvRAWmTYPBQduqEOlKwsF6KHY0R89pM34Yy3yK0APrBVK9UxAoX81ysLPQ36mFypQYOK1Y9NhDv418qaEe1AG4zLc83aVpvXThEE2XiVKLHthIMohL0yz3qB8pCL+B4YjkNacjXwUkEP8Iao73xBl0hwOCcxFAWYPGkqGjULUQ619UQQ2BUVIO1ASuxsNZjSH5okSZKozIoAQXQgLUO5rFvAcKqutQLmlbB/2zCFe9Mn0hFGFcAkN4XvYs4U55HAAsvW/aHLIoR45Icwh68gJYQ5Zx/+0N7B3QUpbRgAAAABJRU5ErkJggg==\n",
       "text/latex": [
        "$\\displaystyle \\frac{E_{b} \\left(\\dot{s} \\sqrt{\\tau^{2}} - \\lambda \\tau\\right)}{\\tau}$"
       ],
@@ -670,7 +670,7 @@
        "                  \\tau                 "
       ]
      },
-     "execution_count": 8,
+     "execution_count": 17,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -693,7 +693,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 9,
+   "execution_count": 18,
    "metadata": {
     "slideshow": {
      "slide_type": "fragment"
@@ -702,7 +702,7 @@
    "outputs": [
     {
      "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAAACcAAAAvCAYAAABg8NNYAAAACXBIWXMAAA7EAAAOxAGVKw4bAAADYElEQVRYCe2Y0VEbMRCG7UwKYCYdOB0A6cDpAIYKgA5geLLfMqSD0ELoAEpI0oHpIAkdON933Gnuzmf7tAcMk/HOyCtptatfu1qdpfFyuRxFaD6fL9C7h59H9PvovOszqD0GQHv0TSiHbdmztvVcpMxms72IXo7O2MFvlUJhfa3FvI9MxJ67Ru93qfsRfknfY2WrlJ/Rdm920TljbroE9b5scBi9w4BgfmkIPoUJtsha2t+oC/SY8oei7JKSqNJNHWsqWeAwaoZOW8ad+Lv26d+H3cFvbUvULcVCnnr6/2aBq8wymd4TkOfc51q/IBIQZIJN4a7G9eVZCcFkDxj+SjGUhm9Bn2FbRycIXEiIssA5A2Dcb2OqeuyeckF73cY/Qv6DEqO+hyKH7jVl0R5P35KyciDTN1HWHp/TzvFcFcrkBTymZ27gXftKmdsgTDkJccosU4DUz69H2us+/J8Y774M0+7zFXVdzp6LzhHW24GLum7nuf/Sc2M/MdGVvbTe7hCOevhNZ2vOh7/TAXz4/bfrn8qLzgEDOgeDY+7iAgPIxiVmAKakOiispddG8HRvSJafoTLUc1dg2HSHWAuRBfnXXn1p8sRGp/SnP65hcBjR4ATuPSJC1+imP6rU/WP6k+IlvaAhYXWPDdlnZwCaljhkRsDFmmAFhcBhwJAcwqNec3K9tvFmFg2rq/ziDFFiYe23EsE+0J+SK9tzpde86KQnhyjASg9bhtLb2kHVJ88Gh044Q+sTV3WATagbiQPqKVOVN8KKUPRe9zbtpSPkKaM0UidkTlS/PtbF1tPzF2MF5gtC8d5Stkfw4r6bwNHh50evKGi4l3ZB5RhTvpOQK3P1W5+/GCswxwuuylD3XToBCnCl0A36geLbxz4lbUz6KnLVnV4rbeQ8f3mmmfXyRNgRYEGN/3MIHPyXckvd1SeiXYQK7ivTVmKc3rhq29mqWBvQSAgMGRI96L7S7XXS3e30r8vb9RM6ws9fGmuAK627oaUUe4CaKHpT8H1JnY2H7DZDK+AAYEKYrX5eDLNkovQ+dNHT636Kuvat9nrRCrhSq/KeAAsPwHO9Nuj5Sxyd4ACi51y1HrNUYKn2osHPX87SyNb6tKXHfJReydz6uJesd3rOCQHnt1PvpcSw/zXpH0J8+gbdfQPeAAAAAElFTkSuQmCC\n",
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAB8AAAAqCAYAAABLGYAnAAACZklEQVR4nOXXS4hOYRzH8c870UyKhctiCiMLooxLs5EshMKGRHakocSCjctykpKUkr3EQkmZsRBC0pDkuplSrilmJpMs5DaMxXOmeTvzzpznnN7MNPOrp97T8/8/3/Nc3v/zO6WWlhYRmo0OrMSzmIQYTYiM+4gl+FAtcB54L15VEww11R5wTMFLOISX+IFuXE76WtE3RNueNXDMnh/ETuwV9r0ei5O+ZuxBLd5gG9qTvp5qwNfhGm4nz+/xMAVYKqziPWFlohSz7FdwALewD9MrxDSiMw84Fn4G84TZ78ZrLKgAf5EHHAuXAE+hSTiAjan+ucKBzKWsPT+MLjwSCs0O/MLdVFwNGjBTqIZ/Y+BZM6/FETzBA+GUr05eqFynhUP3FhNjwGTP/GjSsnRHmHkujfoKNzbh/XveN5Lw0kjAR8WyZ6kO36sNj535GgOltWotFr4MTyNjoxXrZPrdyX+HN+FxznFnCZdPh3DVbi4KXyXU7jzqxX4sxFrh4plUBF6HnznhnwyYi258UcEBZcEb8C4nOK0m4Zod9LVTDp+CXan+DYJ9KqppOC+43EEHthy+BccEA9GvGficymkV59VrBfN5XDAig1QOP4ubglOByfhaIadZ8O5zBLu0NXmux8UkpoRzwkG9UAmchkMbNia/1+NGhZwewSZPNeDVO5P2O4lZIXxAbMLzpC1KD5Su7ddxMnnz+bg01Fsb3qu3i/gnpQO+CYVhufBfHU6FvPpwcMLSn8D9jNxCXj0LflVY8ix4uVcv5AsqJXUJn0V/MnILefVyDWUm2iJyC3n1co1f6zx+4f8AjFJ/AGn1DasAAAAASUVORK5CYII=\n",
       "text/latex": [
        "$\\displaystyle \\frac{\\dot{s} \\tau}{\\sqrt{\\tau^{2}}}$"
       ],
@@ -714,7 +714,7 @@
        " ╲╱  \\tau   "
       ]
      },
-     "execution_count": 9,
+     "execution_count": 18,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -778,7 +778,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 10,
+   "execution_count": 19,
    "metadata": {
     "slideshow": {
      "slide_type": "fragment"
@@ -787,7 +787,7 @@
    "outputs": [
     {
      "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAEgAAAAaCAYAAAAUqxq7AAAACXBIWXMAAA7EAAAOxAGVKw4bAAADZklEQVRYCe2Zv1IUQRCHOSWwjFDfAN4AfYALztBMJDJUQjMpo7sU3gBCjSh8A7CKxEghNIM3UHgD/L5ldtjb3dvd2/vDldhVU9M73T39m76e6RnoXF9fL82bBoPBCj4PaT34zrz9j+OvUxYgQL9mki3aHvzXcSZsqsu8q+ie33WA6tb6IL8gDE4ZMzgb8DMJTvD5J+/7Lr7DGjfwvQnvD2Z2R1qOHAzCHbpV+ufZ8XF4bM2MPeehOZ8On4V+G/kV/EJRwGRCnAPsE207BZjPIBd1kQrb9DjRXgfOdcz3Ls1vnX+jLTKJfT0LMB8gZY1T34jTzJYyukIWgw2/i9I6fa9MeZZj+KzCmXVdyO6yAGUNRvI4deuYJS9GKhUFAtBmbtQSZ8Q3SYBc7BMAjHNeGdSYVRHFDBnwtcEZEbUOkDME53GyHLOC3IAkBP8RxjPpOAzNravBWYljqIpVarYTvgGcZ5rbag3+pdPQGzgrnLx3La8Vi0JPs0DyAYq/eFapjGdRLvB3kK3R50u4h/T+CFvT3qCMFRjmsxKahY3tGuAsgxjH4hZjIoNjhfEJUEnoHqFwQG8Jtzppk2RExrBxsDM2I9mAz0xsXBQa4sz6PODDSqufhJIAMWDt9wZtFpT+6jfqyZbQuIfeWTqmHU37JcaVe9mSzwfN4VbEXGZd46IQcIzEWQYCG18OZucRfHId6fT7/fcMuEBvktlFM1Sk4DhJdaRmjik/18pURFUcmQRnsHWX7D08OTk5pf3i4zP9o263+73o7nYE+SV6jxl5S3tF++A343OvTreoilxbnATHhPHyuwX/Jb7m+fDMuAyCym2GTkLYmIZmn73p7zZYOGqKEz3/iuGusOImuyI5g1wRAy7OLPBlW0ro7NDcXgnBu72S0p2OLUI/Ac5N8J9hH4+M5dyCDJKZNIrMlKG3F5MZ9X36mWYP84vLDHcBdbf3SXAOvUXzAcJ/Jb1DamVwn6aB9L7T+F5SOXuFEB/p4/dnhVoqmhrOsQIESKtcbaVLUU67x78X0lqaJs54BtV6vacK/wNU88OXBWjosVZj/6+J03M1risfoB9IVqP0/jGuPZZ4lz8UIA43H54X9F6W7hWFNXtV8eIbKd6k4wgMyt5tLN0z+79Y1t9d8qzVO5NBOYQvvCD+At4Wa74RwH92AAAAAElFTkSuQmCC\n",
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAEoAAAAbCAYAAADbAhkjAAADaElEQVR4nO3Zb6jeYxzH8ddBtgcLG5JJHVlCLP8a8qdjUqdWeMIDHjjmiQdroTzAkzvJqDX/SiS18ieFdRZZFnUnZVNEbVmkDpnJJlm0sz8cD77Xj9vv/v25//x0ztZ516/r3Od7ne/vc32v6/pe3/s6I61Wyzz1HNewv3PTMwgb8W5zUv7DMLrACQ0JyfgwtaMN+x2WoXU1vaKOWeoCdQe24Tf8iUtq+o8qnrU2XsAz+DU963F8r0KHpEwXtHAYP2ES5xV1qtp6F+JV/IE3sTs5G5Q7RR66GsvxUvK3fgifTdDGQrEIbsFSrMh3qgrUDRjBI3i2AUF7sBYz2CVm7n5zI1Dt9PNnuAKL8Htnp6qttyS1XzUkaJsIUsYnYvZOash/E+wSi2Nx3lAVqCx/HK5xfh024VtM42d8isf7ltk8/WrLxtqVO4c99R7GR7hcHMEb8A4WYDzX90oxWxlX4Ufs7+N9G8WqnGhYWy1VOeqU1B4osZ+BR/ExbsShnP203OeleBrP42I8iMd6l9oX/WrLmE7tyXlDWaBGxLKdwXclfc4XS/TrAiGwL/f5tdR/e/L7Mp4q8V3GQ3hCHAxV9KstYyq1Y/iy05AP1DhuwvW4DC8qLwl2ivpqNU4XgdgqaqQijmBNeoqYKPl9J3vUB2kQbRmv4F5xEl+Db0QZM5XPUeN4QByRO/BchdN9uBZvieX9BvZiiwjybDKott1YJ5L6bSLPjdKdzO8T+/NWLMP7qqvnHcnhYrESN4lgbxVJczYZRNsq8Q3ic1yEE6Uaq+jU24/NeBtn4YIeRB3CB7hdJNBTRULNGFO+5f5v6rR1sjK1LbF9/ymNqsqDLIkvKbBdqvjaYpmYie/xQ4XvQTlTJOquU6mDYbRlheZU3lBVHmTRLArmWtwliredopA7Bzcn+2r8VeF7UNal994taqoihtGWjfVI3lAVqJkK2+b0tytEHlgoisfX8aQ4LWaLJrR1jb0qUAdT2/W9R1xHTPbwwqaZUF9GTBpcWzbW6byhKkdlkb8HZ9f0PdpZIFbgmDjM9uY7VA3+PXHErhIJsJeLu6ORllhB28VNxgYFOaxq6x0QhedKUSIsMtzF3VylLZL3L+Lq54uiTnX/XDgoqtktDQqba7T9e3FXyrGcdxplPlA9Mh+oHvkbyALM94rPeAgAAAAASUVORK5CYII=\n",
       "text/latex": [
        "$\\displaystyle \\left\\{ \\dot{s}^\\mathrm{pl} : \\dot{s}\\right\\}$"
       ],
@@ -795,7 +795,7 @@
        "{\\dot{s}__\\mathrm{pl}: \\dot{s}}"
       ]
      },
-     "execution_count": 10,
+     "execution_count": 19,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -840,7 +840,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 11,
+   "execution_count": 21,
    "metadata": {
     "slideshow": {
      "slide_type": "fragment"
@@ -849,7 +849,7 @@
    "outputs": [
     {
      "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAAADkAAAAVCAYAAAD8dkbIAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAC0ElEQVRYCdWY61EbMRCAg4cCIOnA6SCPCkIJOKkgoYMw+ed/GdJBoIJM6ACoAEIHoYN46MD5vstJo9Pp0MXGM2ZnhFa7q31oV3syO8vl8lkO8/l8D9ob6eCXOX/b1vg4xadXjFvwu9y/SU5A6DO0q5Z+k/O3cd0GZnAn4BcMkxRhJ80kzEM4Pxn74PdRagBB5jesS+ajAZFHIaP/pFX0h/klw2B6GVMG+gXTPfPMtZBn8i20OwTGBOhpWSZNWTNvBPDlF4qvmY8Z38CPGWZL2yW4hWjpRsiDjIwaghEPwoy/rsmuykf3J/buMZ8HHa1d198DrTavHKSKW4M1G+vwLTkzk8M1hAPsW01VWCvIqvb1BQ5QsSioCfdRfhV2qxKJACdnA2hKKCGn6BEypymhhCNTbVgjs/S8pD+n5UGafu9aDzDqHZBnCXm6BmwTiIBMqbQiX6R1fgpaa1ghgKI/rVL9zcEOHPY2vBgkxjX8ntFxXCl4dis7WtoApFeDcn8K7LG9j/pEpfsG8BcFuj76iTlkNP7ughicWXGese69cKAZTAyItUE/dMKwh4H9Y/YuhjXETJm1DqDbT+A+xDPmL8xfJxJFGDcMv0WlEoDVgQ+s/OhuDPAjHETJn0ALDSj3w+BMmvGcT+SC+OYLr5YzaRXwZeShbBqsKp3NIdy5UtXpm0/Td8TU8JsgEw2WpKU4CGzU6NSDGRR6PIZPzFKD8gFiYkK2U4u+2rz3kZcHmQoP4Z7UUJkM7Yl0jPuCWTJ8rj0IyPg5WjBrswFwS9UG+fEfpfi3c59jdy2Klome1OgnVa4CJz1lD2lsuZs1u6V2bTTOluLoSvrvIFE+w8hagA5/SYwCZC270C9G7cmFVinXXMfWr/Mge6+FrY+g76B3NjYd2XmQvhBsDPGiK/RUAL9DU+r0jM5/BgwGwSmTTzu/RT9YG/hWQ+uz91bfez7/BYNb9kRn92zfAAAAAElFTkSuQmCC\n",
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAEEAAAAZCAYAAABuKkPfAAADVUlEQVR4nO3YX6hlcxTA8c81QrkYQ8qgLsMMEpIYTTRumW7Ng/HggZIRD6K4eMEoZ9TkZZqpEZGXCQ9KMkpE4b7492CiRjNMtw4aI5o0g67bxXhYv9097XN+e+9z59zunfKt3TqdtX6//futvX5rrb2HWq2WAbMiyclBTzxfnDgPc36Y5Mg8zD0vnLDQC1gM1DnhTnyOw/gHVzeYc8TiiYIWZvAzdmFlL6Oq43A5XsOfeAMH0mQLyfl4BmM4CwfF5jbjtx72EzhFPLxbsRzXlY2qnHAzhrAJO+a66gGyAp/iHLyNfWJDDwunrMGh0piJdMGXuBbD+KPTqOo4LEty75yXPVheEA54CBvwOEaxHauwpWb8PvFQzywrqpywJMmZjP5s/IujNde0CMlj4SKsQxvPl3RPiyN7F06tmKPYx5Ky4lhK5LA4nwUjuFuE3Tsd//+KvzJz7Exj7km/c4wm+YFwfCe/4xPhpNVmS3RjqpywNMmpjL4tsm/BfWJDr2NrvwupYVWS32X0+4UTVso7oXgQZ5QVueMwhBtFOH/faJmz5fOrhvbwBC7DWzV2xcIPZ/TF/0sr5mgnubasKEfCGG7BTbgGL2leFq9K8uuG9kSJO9iHfY6hJI9W2LyK+0WUrhHR8zLa5UgYw6OilOzBc30s4kr8JHLAoCmedFcoJ04v2fXiAJ4VCfJ2PCk1dWUnjKcbbcDFeF+PbNqDC9NC+omCfvg2yZ4dHy5JMpczYD1exG5cgZOkHqJXTjgimpE3cZ44s3VcmuSeBrZz4eMk1+le82kivKdEi5+jqDAtfKOj9Ff1CUVCXFZhU1CE45EGtp2cKxyYC/OCSVEeR/BgSbdZ9AeviH4hR9EktcuKKifMNLApKMJwHNtwR4MxxBndi9sa2D6AX0QLvyuN/QiPpPtvqhlf7OPvnKIXVZm2zG48JRw3LhLroJlM8+7E9XhMvE/swA263xtydO2rqlmaTrKr186wRX3/XmZjupryo+gu50Kxj67utSoS9id5Ly6osV3MnCzeNteKnNVVwqs29q7I9uvxg+YfVRYTLfHkvxDJe5vud4/K4zAlzuCoKJPDFv6jSr9MiER4CJ/JtPR1b5HTeC9dxyMTZj+qZDlez/lA+d8J+A/IOKnbkQwY0gAAAABJRU5ErkJggg==\n",
       "text/latex": [
        "$\\displaystyle \\left\\{ \\dot{\\tau} : 0\\right\\}$"
       ],
@@ -857,7 +857,7 @@
        "{\\dot{\\tau}: 0}"
       ]
      },
-     "execution_count": 11,
+     "execution_count": 21,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -909,7 +909,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 12,
+   "execution_count": 22,
    "metadata": {
     "slideshow": {
      "slide_type": "fragment"
@@ -918,7 +918,7 @@
    "outputs": [
     {
      "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAI8AAAAXCAYAAAAyVhy9AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAEgUlEQVRoBe2a7VHcMBCGD4YCCOkAOgiTCgIdkEkFgQ5g+AX/GNJBKCGhA6ADQgehBEIH5HmM7LHv5I87ZO7gvDOOvlar1erd1erIytPT02ig922B09PTT+zwT+pdrqUWOMhbSAuco9U2ILpLqd1qSmGDrMWzQIg6o9TAcacDeBbvvFNrdIxAI09yGq6t5CZdHIFEm0202aS8nkUr5q0zT/BJypK+0/9oZQCPVni/dMTW/Galc4BykE+m/pO6ifeWfSupX1ss8A+5IlZ03oePIiPRa+Z/D1+mwHP38G9qC2Bfz+CGcntW2cz1Kb5LmUUuSs/vL1+WfCeNPAjfQbBKH1C/oCyItv2iVlDNvKFC4FBps4B5zlkbU8u4Uee2jicpeFjkqx9AuYwseEPfBp+oFUBLR+xb5zIHqThWakMgX0fdoSyunFnWiOipPG+N7MmfGjwbCJ4ADn1XLOp1JXC8ypaVPFS/vin5C4tz8/z2+IpbIzV4JsIki5pk6XHenUl/pELm0hE29AAfKZteUHuMJ8spkWWuk//QWNwaScHDIhVw0HbBfT5zoKbNwtIvlXSp8/yJPK1fjaaXzh4OmWVUMXoXEaAsKfDosFGa1g7wC5wjyl0FhrblfVLwlLVFuKBxsy7c6x1fXjdWZ32NqceYkz3wCerKExaeCvAZXyhCP68N7fiR79B2jc46QTTq0D+VHeAXOM7xDF1fMu/JbNcLeFjI0OqiF9R/UM6NwqavKItcjPqIrwIW2nmUVNcz2hW9aede/4txf/941dyN9TJ9KU0Ncl10hoIY02G1+wQx5uG32mFsoq9jI3Xlj6rIEkCjlZOTEwd9CdWFc/nGyRdVxfg5A/071E2QL6mPb67OW/LpI+Yk1acQHCrI14jH47o5TJ8A8mA+UH+0Lyfa9l9TRved81nC4wFqh3HaCB0P4wO075hXsVeEJ+sK8gXKFvUCxNTz32AqusfkwFtrhxh/rG8NIS4UvT9jE5r6gkICRyPHDCFiM9TWyUmpT80a3+hXxxiVvbq41tDJg9IZioOKTc774IvukX4jsk/1SlTL501R5lFSHbO1gmx1bAVOWKfJDp1UWe3E1YEJpb0fjWB6UJZclafRpydWwl95/BXrHmD0h69gePMKo0xG9AmcW8pOwAnTei2CLj5A9qmvh8WOKQV/V6q1Q1cBa10Zm/jCBvRmw/GXGl5DeZIIVyO/tRs9Bbie33T1ZF4NjwDSi4sfxVoXeF0G9dQhBZDAFuCdog58XezQuptUkSf/9djfcioboG2eY8TxECpjrdqlZ9DbGiMIOjquV+vJ6jzXnxjQIUpBL51APf0EU1dqtUMXQS+OPGziNwuZfGn0n7QpMjI5FOHrz83mXCfw9F18ZgEjYBtlAFtU4JSU95rS/oK80SlKc6x2tcPYtGozBXhiiXF1lQVpYeCuunoddAHZXHfGfkyQjT5Fct9FoSns0CjuxeBplP4GBzGskdKI2cd15bWd9OpG37nlkcn/P88bxEuhMgfhy8roZOTxOe01PM11wJTlof8rY6gcjiQb1gAAAABJRU5ErkJggg==\n",
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAJAAAAAhCAYAAAA20CJoAAAENUlEQVR4nO3bW6hUVRzH8c8Ru0BlpV1EMZJu2uVkl4fuVD74EKcLQREUGVRSHMIHCaKiY/fsJglSD+Klwh6CPCYRqd0oMLI4QUlJkhT5kFFkQjfTHv7bHOfMnpk9Z89sR/cXFotZe13++39+s9b6r1mnZ2BgQElJq4wq2oCS7mZ00QaU5MJR+LWIgcsZaP/gGtyAnk6nUkD7BzPwdhEDlwLqfg5O0u9FDF4KqPu5Eu8XNXgpoO6nDyuLGrwUUPczCd/n1Nd9+BTbsBVv4sx6DUoBdTfnYCjH/i7HQlwklsYdWIOxaQ3Kc6Du5moxS+TFjKrPt+A3XJw2TjkDdTfnYX0b+z9CaOSXtAp5C+h+7GoybW/D+AcSE7FF+LJdzBdL5Lq0CnkvYRswt87zo9EvhDOInTmPvz8ySm0/9cl3+armaVyGS/FvWqW8BfRGkmoxDquFQ1bgtpzH3leZicW4Qvbzmj48g6mGi2i62KOkcQx+Ej851ONvHIk/K8qexc3C5m/rNe7UJvpYrMVZeB03iR1+pxiJM4vifBEVnSKioo8qnh0mZoV6dh6Ohys+n4hb8RlWVZRvrernBdwoxLOhkZGdENB4IZ7TsVx8a1KnxDbRqjOLZH2SLhDRVqWAZuCdBu03Y6Di8+3inV8Ts1otFoqZ51qxcR6flG9P0jDaLaCJeBenYplYtorY92yW3Zn7CiuEvfdWlF0lDv2yMC3Jh+rUuSvJ11aVz7W3//6nnQI6QYjnJCzCnfadTfO0JB8q0IZmGcQ8TMHXYg85TizJWTg7yb+oU6fREj+M6jB6s+bD8F14JaXfyfhAiOdF3CGbePKyI416zhyNTWI6r2YefhQ/HzRr9+Lk2Xs1ni1pwtaN+EYsY3AhPmmiXSU96BVh/9aMbetSPQNtkm0PsKVG2cli5pmEBbinBbvysCONRs7cgSeE7QP2fNP7MUuEtj+k9D1f3A6sZJq48LVUCKySoSZtHkz6mCeEtKzJdruZjDH4OGO7hlQLaPoI+ztNiGeCCAXntNjPSO2oRzPOXIoHhfgfwHXiXKRP/SVgfo2ymeKPv0Tr1y4GhS+PE8HIVxnbT0nyL1scP5U8T4LPEMvWBDypdfG0m2ac+Q+ewt0i4nlV7OHWtNe0VNbhZ8zW4FwmhTFJvi0vg3aTl4B6xRp/PB6RPULoJM06cxH+wFt4FC+306gG7BSnznO0dvq8Mcln4zlxDpcLeURhE4V4xuI78bIDdepvUuwfo9KZY8X9l+U16v0lwtmpeLwjltVnJa7Hhy20/Vwsxf3ivZ9X+50zk4eALrHnvshkPNSg/gLFCiiLM3vtfYBXJKvFJrrVE/zHkpQrPeV/pqZyiLioPsueULykivI6RTq9OEjMWCUplAJK51yxD8oaMh9QlAJK5yUcqrO3BrqOUkAlI6IUUMmIKAVUMiL+A3fI+3AHcVB6AAAAAElFTkSuQmCC\n",
       "text/latex": [
        "$\\displaystyle - Z - \\tau_{Y} + \\sqrt{\\tau^{2}}$"
       ],
@@ -928,7 +928,7 @@
        "-Z - τ_Y + ╲╱  \\tau  "
       ]
      },
-     "execution_count": 12,
+     "execution_count": 22,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -961,7 +961,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 13,
+   "execution_count": 23,
    "metadata": {
     "slideshow": {
      "slide_type": "fragment"
@@ -970,7 +970,7 @@
    "outputs": [
     {
      "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAO0AAAAmCAYAAAAoYhz4AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAKxElEQVR4Ae2d73UdNROHHR8XYPJ2YDogSQWYDkxSQUwHcPiUfMuBDiAVvCQdhFRASAdABTjuIDyPvBL7R7urvXvvNb7RnCNLK41Go9/OjLTavcm9jx8/Ho3R8+fPv6XtG9LXlN+P8dX6ikBFYD0C+NglUr4j/UD55zGJ93JOS4czOrwhXZN02D/HBNT6ikBFYHsI4GvnSPuJpO99ybV5hwZO2zjs73D9SvnrDne9qAhUBHaOAH53yiD64BHlz/sDHvcruP6BZKenmbZaVRGoCOwYARzV1dXH0jPKbpk7lHPaL+B433TsMNeLikBFYD8I4H+/MpLO+6A/4km/guv7pFWHTgz4BzLcXhstKlUEDh4BbN3d6SvSOeV7W5yw/tih3ErbYVh60Sh/Rr+HS/tW/sNEAJu4JHnAcrDE/OKWdjBH2jxY2hrtwmlV/jMUHSzrW9O6CrozCGAHFyj7gNzt3qHT1cgEfYXj25itUG57vFowCuq4lT5xBLADt4zfk9/pAI7+7hxdLc3jQe3/KDu/7+bsnfY/Sa9I35J+pM8q2vpKu0qb2vnQEPAZb6tbw9sASKdjXD960Gk9q/mR5LVnN29JswS/H0t8Q66jr6KtrLQo8gEtVOaa5ATbH2M4UU+kjTaDd07Urybk+uXWI5JbMSmcvFHfec/MtUYkj3q+IxklB4du1O3tIG3PY4nTV6R43iAGYhHJe+j98vDjHbrJuxHRVzkPS2TsE4O5yczock17sm3KOq9bXw+fSrb/BrCXpI5dzunUb1/ttCqMUG+2UaTz6RXX1vuSWMPY2RaJccKWg9xvMg0OY8b2f9qlp/Co04CoV+dgcIPGLVfscyxVZzxx0tDEydd6WZyoj85tt03JlahjDzlBjLU3vHPjt+s21EU70l5KSDx08lNS1v5KhByXMM3wGDX81DF3g9w6GLU9iNhYyZnxQzPyXc2l1zdZ9y/t4SU1ubqO6tK07eUgbZ9jRTQY0yArTa0M3su0ogTu5X/Ee/bw5TYwGJvKhrqcIq8Iq0a+O7vHYzqU1G/Dae+jzMBRqPOG6Uh+P1k0qRKFJ3iiMXYMhbGNah4e+Cwy0DMnD75Rp87xr6nb51iNnnF17eCUmYOPCBsRc4oB1O33LN0CBqM6zeiiLemkgSi7I9GupgJgw50yMVm1Pc45bVIqDTNdeNFvZhLu3XWirygPnhn7/Fu6fqKcNoCUNR5PL3123Ufg2NJUdiomBLc2To7GtQYYiLJBa4khNj1T5hg+puwt+KWRd194zLwuGrw+Jw9BkFy/cXE4ojx1+GYwjGcKsk+Rr5AGW++Tdg8GixFyLgqnbvTpOCXXKu7WyGfcRTce/jUHQOqedEGWOrgtL/4qC15leFglhZ8jUido5+S57X9gXPoHWRGjsQA5OB/oj4GMTbFyjp37giydrONg1CUs+2MXXD/qyxvr02Dxd9PuQeXsK5QxWWvrC3TxICprB9SLn7Y2Z2/iGlbsps+U2u4MfU0kf7o/HaeF4SXJA4qsYlPSbaOfjmLEFvhFMuDXgHWQ0igE6w3RV6OTgjFyrQ46RvHq2ozv/I2SvoPzwECA/Wh71XYGGYmQpXxvgDKvSOrpoU2iZtx03S/QvhFW9Is4OacYnJRl/TZP9pXp3CYJHVwc0gl+o594zBn+pNxNGgt1cV5rKeKirWtfo4ROYuPbDu9VfKw5Ck5Lg8rEm/jlqJSJhka4Bvkz5cUvkOljFPMA6HpimLGmOKG/6a+zOVmN8JLcrUzJs+x9+Du/X6SfBvR0bNCl9chzlXvT1ofyEWny5vXHgX9TrCJOYRehXGRpPOo0CHDUqa+B7Ipy7MvlLInlQF67VzPuOXl77gavaIdt9p2W53Rp2r9XCcrBvlYoFO1bjErIty5vGdcAFw5RT7gwyupsP1Fe7Gz0O2pkCPZryp0oyfUXpPaNsUuW4IsTyrZPVDoHSWOOc9DhXPkFe9Zp6ZeMjLJBzH4vVuhE9y4hSxwSFlzrFBvNeUO9Ak6NHkE5yj57dvDhOmzH5LMdxoRN6DT/R/yu5tmC7WiM2o4HOpOBQb3ge0syL6UUoOY6ID+rS4PBtnZbRbhEXRlb+/AxT1v+QP7ghII3xAYfqsPNolxM9NHwnKyg5yZWss8vHm+EUR0cP23JKWuMbpfPyZcEDlceV+fOlnVk3LXVHp6J3b4o4JQZzKAdiHnrEBpIDH46+tLdl/Z0nzRKjOP9cQwfZWIwiV8aZfvBHww427iichNdVgwXcSl2XvRzi6w9hkeJ40Zhl2AB+Ytrb2wRwauBG/myL+pp92aEX+AXCdyAqRnDnjnjd7WVXDVnqZGlk0eDPaKsEe+KvBnvdiW8LbeZm1UDnGhrr6Q6bAh+1GsLzl9MPBDx+9mt4YEsjdCfsbnCGmDDoQv53mmPukT8rksmiV7a8EuSO4Zgl8exo0pTFjidcJbgd3ANwIgxFomN4L+QdklxSxUOodoDoaN1GqQrZwSrzZLKtOtA7jTSat00CliiOTmJcaaAHAOeB0LvZ1i31TyKUxyg0ekReTQog674xXOK3yjHQBi75XL7jOKNfJ8LPf0ORHl2axx5t53fgi7FKy26ib87kc4W/6QHgo77B8xGvLTa9Hjipc6tAoOvnehrhNbYi9/V0ceb7DfMrtpLPnnU2Y7oM2b8Bg4NzdXW+Q2IvurrNt6xzQ1GniArO235aQs6OhZpiY6IGZCy2yvcgGGsIupB+xI9JnFCpgbis2WaL2UdXQeLTmz+kDRHOvfU7saxvC+JGEP9DA5xrNS248K+dTFYe/ZSMk/tNT7mJRg6Toug6GSPEkemAJ83V0PX6DzAIgukE6vU6c1l+dG9k3B8+s1uF+FTvkHB3PGOqFOn38hTsKGss7rdkwxEwTnJ+87yhLqwEpFrPCHgkPu+NAUDylFHV2RTCfCOnSMx7hhujilX19JjEit1pL9zCau6sqgTpzbJoyOaGwTaOxbr24HOQNXHjqoBidkURk9pd8vtvXFcSWw7h5g31Tv/u29dxLCN8dQExX94j/0nVNvp2bNnH0iv2nW1PMDoAoxODx0X5sgU/517YxsX7bqxcsN7Ptb+qdaDy++ky5L5Nxi+6fMeT7l5bRtFoP3cN8p0AA3uvMJK2KyK/lzvdeG8PMsIu5dC/oNna7B0t7fknOe6D0xne9xvrNdDBBrg42d3Q4bDqvH5Nh0aMfclTujW33OP9vb6sNBZPpvHdPFbhoEjLhFVnXYJWje8l4CenpuXd787PZinz6YbPWfal+SpcOkXaXcHmM01NYAtCXzZker2OAvLeCVG+Ek47DgCi1o85Jk6RV4k7C4zYzceunloW3KQNznV6rST8NTGNQhgoG4D/RS05N3umqH+032Z/xkKpo8j1ipbnXYtgrX/JAIYrAdXvvs/n2Q87Ea3xe3336tmW59pV8FXO5cggMP2vzIr6XYwPMx/o3OBMQByK+0VzOGYf6xTra8IVARuD4Gc0/oc4pdNlSoCFYFbQoDV2YXT5CLaoZzT+omVP2U763DWi4pARWCfCMRT98EhXs5pX6CZx9L+awZ1m7zP21THqgiAAH7nL3vir3sGr4gG/xO8qDXO6ofl4YNxrj/pgwQxqVQR2DUC+JmfOMYfdHR+jtceO+u0kQEhF5SfkHzX5tcxlSoCFYEdIdA4rf+O+OSvgP4Bb6z3LzRhhjkAAAAASUVORK5CYII=\n",
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAANUAAAAcCAYAAAD2p/r/AAAIcElEQVR4nO2ceYwURRTGf4tKwCWAgKJgyOAqqyiHK54RlYB4G41BEm9QMSgexCsoieuBtyIqalTMRlQkigIeeJ8QFcVbUZR4hlVAzHqgKO76x1ft9vR2VXfP9MzOrvMlncrW+bqm3qvvvareitraWlJClUlXptVhGWW0RXSIKN8cOAk4MkZfL5qnjDLaO2qBPW2FLqUaALwArAGeTFemMspo05gJjAHuAiqChS6lmgNcDiwCGmMMlDFPqaKKZoranscsFkrp3fKRpY7km8Ya4AKgOzA+WLippdEwYAjti85575Jp52MWC6X0bq0lyyvABGCWP9O2Uw0xZf8UVqYyymjTaAKGBjNtStUNccXOBRSo2MhQfEvWGmPmi/uB1UBlRL0MpfNuGcJleQX5PTOAn81zI7CJo6/dkLKcGmPcSqAjCuj9B5tSefmbWcovNQPHeX5zjJMEz5v+RoSU9QZeMuWLgB4pjFfK8ObC9UzNod9hwAnAtcDvqUja+jgerb+9gTOA04HJjvrLgPnAVUCXiL47mjRrfdt8Kg9NlvxPURDDhi2ASWawBcQLdEShxsizLJA/HJiLFOty4IqUxitl1AAbgWmOOo/l0O/VwC/AnbkIVaKoB85Ba+czFNWejHYsG64B3jLtrnbUC11nNqXaFPgVTXAYHjdPGHoiS9oBafw4h1BxUYV2nxUBmc5HVrUBOBR4NoWxSh3eXLyHzkvSwgBgFHAv8EeK/bY23iR7c3gDuBLoin19L0UKeAZwHfbYwiqTZumRjZbVYFcaF7YEXgZ2BR5Fsfy/cugniGEmfdukXZElvhF4F8kbVKgluOnRqynI5aEXslpRlGwD0CnPsYJzkRbGIz96rqPOcDTvK4E/ke+1FLc1LxQKLcvDQD9kaGx4DtHkIf7MsJ3qAKAaGJ1QiK1RaHMgOuM6kejoYR1wMtrN6hz1vIX0DjAYKewOyAk9l3DFfRztmEGMQ5P1coRsSdAF0U4PGfRey8g+A1mDFkAY6kg2F2kr1Sj0e71pKb8E0c1vkQFbiyj3MOBgU14sJJFlT2QsvN1qL7TD2HYpD0tMeiB2BrQa+WhXAIcA66GlUt2GNG8fkjmqfVGgYACKHo0jXb/GW0gDkSVqQkr7gKNNGGe+ASlUHdlKkC++JpuKnYYU5GGLHPnAm4t9gG0tdW4g2e9XiULDyy3teqP5WgyMpKUR65VgrHyRVJY+wC3AHcAg4EIUhIiCZ7T2i6g3B/jIpLcCLwaVaipy9ichLhkH/ZBCVaFDsAnEV6gpyCeqd9SpQHQSZBV+R9bm45hjeH3MBCaa9GzsQZg0MNSk7ydok3QubL7qOpL7Wn1RmNk29o6mfAXhrGBtwvHyQVJZHjT130K/+SxgeoxxGhCr6BdRbytk6O9BetDCp2pAkZHj0AKMQn/km1QhKnY6yXaoeuQQNjjqDEDnZiuQ/1QJHJVgjE3QzjQRuB4ZjEIqFDRz7A8StEkyF4uRgoU9PUPaTSOcCnvw2vxsKf/EyDUeWAiMRRHeOPia+McvTbjZRy6ybES/eXdT73yy3ZJTgMMtbdfh3oUrgGcQPXzCyB8aqGgywl7o6Axge6RQGUQbJ1KYxerRnTfQy3+PojfHx2i7GaJgJyHrfXEB5AuiAvl9q5APlSa8uUiirKCd09XGi/bZgihrgX2RLzsSzekadCZYEzH2SuDzBM+q8G5SkSUpOuOOhI5CzGG+P9MWUm9AytKJcMe6Gm11fYCb0OXCQsEfpKhHirUYuA/4DnjN0q4TmvjDjHw3FVBGP/qj6OSSqIo5IB+lmuMoX23SsF3Ow8comtsR+RkTzN+7I/q4wdJuZBJBYyJXWZKgA9rdvnLUGWjSLHbhuulQQbjl2hntUH2QD1BIhYJspQItqLGI1s1HHDuISuApdHZ1JsVTKHzyJPH54iIXpdoS/VaNKDq73rTfy1enHln76hj9/YU+CToWGbeeKHjQGoiS5QBE/XJBNdKB9x11Qq/x2ZTKo3F/B/IHo1B0b0TBpsQWMRzboEXYzVLeAVnZjWQvpKeB8xBHfho5ix66ofOD/RFfjnM7oA698ynxxHaiq0mjQrZBxJ2LRhRtigsvsHEBiprVIIo1l2am0oR2/F6I1gfbh31WsT2wCwprf59AnnxQTFk8o+M6evHiB1luj43+bTAN/OHVvmaAHmhLbMQdZVoJzHaUg66DuM5mdkJnQB/Sktvejs6qzkE+4AhT5yEUbl4KbGeR8RqyKYJnXDZGyBsHK0x6Hpqrt3FTL79McebiF9y+4d1k+yVDkUU/GvjG5F2MjFQG+NLkzQOOAQ7y5YHm92Q0n58gqtif5q/Bx1O8a2HFlGU0CmgscNRZb9IsF8mmVF/RcoHtS/NF1f7AZRFC3Ua0UkVhN5O+YymfbGQ5AkWNxtB8rrCHeYJYTUtFG4SuZT2Vh6we3kVHE5OQYk0nnlJFwZuLrtjnvpGWVNe73fKNL88zln6mMg/4EQV1ZvryF6B1sgea305IaR9Cxy5fJHmJPFEsWbqhCPOTyG+34Q9TnnXJocLyj1+6IGGrcZ+btAd0B35Ci/Gi1hWlIFiOzmb8h9Bj0LlKL7KN5xR05lKD7hb+X3E2OsjdD3jdUW8aWj9n+TNtPtVvKJ4/g/zvqpU6hiPf8ebWFqQA2BxRZP8nPBXoatdsWrKR6cgvSfO2SVtDZ2Rc5uFWqEFoN7syWOD69OMeRAdmmQFy+ZSgLeAJ2q/hGEzzla6X0I5ci643hf2HrD9N3REogtpevqlKggzyS+scdaYjQzUa+CFYGPU91ULzlNE2MRQFjKYCjyC6twhFttZZ2ryG/ezv/4DlRF/zcn3kGKlUZbRt3GUeaL9Mo+SQxmfuZZRRhg9lpSqjjJTxL5m+IQVqKuYZAAAAAElFTkSuQmCC\n",
       "text/latex": [
        "$\\displaystyle \\left\\{ \\dot{Z} : K \\dot{z}, \\  \\dot{\\tau} : E_{b} \\left(\\dot{s} - \\dot{s}^\\mathrm{pl}\\right)\\right\\}$"
       ],
@@ -978,7 +978,7 @@
        "{\\dot{Z}: K⋅\\dot{z}, \\dot{\\tau}: E_b⋅(\\dot{s} - \\dot{s}__\\mathrm{pl})}"
       ]
      },
-     "execution_count": 13,
+     "execution_count": 23,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -1036,7 +1036,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 14,
+   "execution_count": 25,
    "metadata": {
     "slideshow": {
      "slide_type": "fragment"
@@ -1045,7 +1045,7 @@
    "outputs": [
     {
      "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAAALYAAAA/CAYAAABNTzIUAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAKSElEQVR4Ae2d65XUOBCFmzkbALAZDBnwiGCHDHhEAGQAh3/840AGsBHsQgZABAtksGwEy5IBez+3SsduS7Ktdk+7u6vO0ci2nr66LpdKas+Vnz9/robkxYsX58rzTuGbju8P5fd0R2AXCIh7kPWtwjMd/yi1cVZKJE0VPFb0t8J3hUcKLo7AvhC4q4YfKPwjXt4sdeJKSWOr8FMVfqXwVsdPShV5miNwGQiIh1fVzhcFrIgbOv+WajdLbBW4UIEPCl91fCtV2K+VERBuaJWHCiiIWUV1X5m1wgOqLOAKuTGNb6S6/kvqYrj2JsRufhRAGkjibYc9+GwgnydPQEB4omxfq8hTxQSOO5K0sZURTYOqp4KvnRJ+MgqBgOHK8RsFV02ml6HQ81ThnMbm9Yn8sY78bwUCAI7Gnix6GLAjbcBQMMgjXS96AtbZTuMvWChgX58rJnRs7RyxDcxO5tOAbPu7BOgA+MfK2l6pjjhZ1zFmITZl0p6sbOMYijXE1o1gYXS4enYMd1dzDyLLVYV3hJryA2Wwqbexqx+rX0zeTdD8aKWii8syn1CcfYOdMrEBhYnxPRGmymRIEUh1YUbcVlyrrakWbf2ZA5c6BE6W2MAl8kHu9wosQs0lPCQ2samqU/1i3aCtjSA6ri2fyI9E9KSJHTDCfsUsuTcSs2w21YG2vlDMwzKLqC7MD/rmawkTED15Yos4mAxoxzhZm4DfZtZqT8hmRZyrb+eKeAPc0nFbg5PsUkAg5xUpFDnKJDbW4OhHcycJpOtoTVxMJdsZez3ruVAaJMXsQbOn5Iny0JeVYkjN4g77I+ycuDP7J82lj4ATe40J5gjL3pCut4olMpGGNoZUSZMg5KGepCidNB4adkeyoQySdzwnytPY0IohNfkhtnlCeKN08uvcJYOAE3sNDERDIE+H2IFYaNFfFdDqNxVSkzi0bVJbhzo+KI62t45XCql61Ezjs0ar47uOovxzmEuxvmM+OHliiywQCF822pCFkQ5xjXyKXyrdNHdnT7rS0PQlbQ2BI4lpQ+do76Qo/VoywS+ORuCkJ48iEK/8TwoQtbFtFSe1ovJCRPJgR1OuLTwUVr59PXf8UAnsnHTZEQKLIbbIwsSN1/XwT3pmAEPtoDUhV7MHQ+cQF1PhgUJOsIuRaOuqHJPK96F8kzjiD2U+j8jnWSoRWBKxIVZSW1beW7ZYIDWaGrs4mgg6x5zI+rSVl8kjXhGWvK8qRp4rYKaMEpVD27M83m53VFnPNB6BxRA7dPn7+K7X5RShLlSSSRmauuO6C+eQt/SAmdaG3I3mVcxDOVYo4y67sWhV5ruUyaMGHi2FNiSGGGg7vAzEgz/MVJ45BTMCTR09FBuVk87mqN5WSPLp+kcFtC2aGoJ2JpI6H5I7ygAWLjtEYHZia9D54S+DH7WejtnnAGHQlKQ1GksxXgZMgqRvWNdnF7XZLHjkKlY6hB/62RWmh/1qf5L2Vf1TH4RcV3d2XX3sjeEuGlM7KLbfqXtuXHKmCA1OltBRtPLtTGFW7iIRdIzPGPca5sHBiPoL+dHacRJ5MJ0f6KjujbEvjeFADeOT1dYP5Z59hyU9yBH7eujepAlO6Og1xVM0MDcHkAcl3KNCfEgPqvOFzlaOYaHGclJoD0VRs8PS8O/xJ0fsZgGhZuBCR8t3001FQ1gHuyl+thcEKsZw234y58h6owqV/xXSeuZlj9i6KTMLOh6DQuVTkug8RG5Ex9jY2Ny7aCu04tHSEQjjz5s7zsvG9Fnl0PSUM87GYqnJo7mzsH0miRqi7L+hEPsmUh6PB8qHW4/XBx88aZ42xRC+aVvHbxQGb1J5LmWSo36dhAhP3tT2U7n7OucrBYzTheLiyqrStx0L6i/usMwMAhN5tkJQNu7z6RBbCdg53BzuMJ6E0aL8rOJB5MYuV8xTBFHbBGXymARI1+2JbefPtq/8PAiAnpuoZst6Qh+BgCceCswCXLGQhbHE7Vn05Mw0FrTLGzy5w1LXk6K2XyvgQqW/rAA3Zu2Z5dYFFg6oHHImyWd5N2Plt6e6Pdk09147O2ScRdQmD8LUieosbR9pJTgMfhOuEAUOmN06+OaeaSx4iyOjFNs66/qv2ufBg3tsyYCLq+YTZzpBs/KksBqH3TJJQmXNq0gFeZVFXzUVhXTa4OFpgOO6yzIR0HihgJ4rvNQxCmSnEtqDN7z14Qkep7aSHNW+yhiP7xux/1PJP5Uw+WmxFkOlvEpMFkdg9fFSNlgZALuMdS9Di0hVzateNB47GKO9WlXRyEKhPUj9WygCF6s+ghrq4uH4aMTmCf0UKmbSUO1+U1lsa8wQYkyFnT/xod8ebYlAGDts6miK6hhP1k7GUPXaZLWZqNJ9XYPkTFYn7UlXfqwByjYK9UwHVMakzhZVvugYoo8S5cVoxwxpRMeYIWaf2WWPd4yAMB89ZqmuqDzEgMSR1CEfE8rZRe1AapTppB2WqY6EvkdSk6fjFdE5rwNeBTRoRNdhUdDMTDqjhIY2v40R0+c4UBsMJH09+c8cGxaKq7BQucYTFrDEHOV1jmcEshc9IkpfWfs6HNW+8sMZ2kBTd9YwOFfAYqAfg/M95cV0gtS0Hbc4NBpbFxtRAq8c27/BTY2RZtasso8V8CViZ19XXG2vj2lU9dNXAPg8Jv8x52lhgcat0dwPVe6uAsRASTFxhHyb2lSX+tJqf+xY0A5154jbmLJKh7RDYkq18wD2PvweKsO0yH5Ue6glT98PAho7lBEaj4f+6KXF1d6boqOxQUKZ0YIAwySi5umnGpf9IHBHY3YSpA7wYkIhHXOGCz1ic1Fir5Tb61P/u3QEghKy7QxL7+5c/TMnRXReWMU5YttT7xrbkFp+zBznUnzPC4LieuiLrVrGruWIHTP4wWEgcIKkLg6ME7sIjyceKgKbfuxDvY+99Vua0vYn5Mw23Fqbix576++pNOzE3mKkRVh8qMxH8KFi50HyuEig45XyTN7MQzmX7RBwYlfiJ8LiapryocnKlrxYDQJO7BrUVCZo4qiNA9HNm1RZqxebCwGfPM6F5PpfS7P/wWUBCDix5xsElrNtYWu+Wr2mKgSc2FWwdQvJDDnXFbYgRNOkm8PPLhsBJ/Y8iKOtq3+cMU8XvJY2Ak7sNhr1x3dU1LZP1tfiJWdDwL0iM0ApE6SzF3iGKr2KLRFwjb0lgF58mQg4sZc5Lt6rLRFwYm8JoBdfJgI5YrPvAbH9rusz/+sILAsB23jWW/HNEZv/PIDYLxTWZ/7XEVgWAvYLr97CWJLYmuXbNkv8sy6OwOIQEEcv1Ck0dvLHy0lih7tofmakCtqfLVvcDXqHThYBtggjnW3C60v5H/OuRGgKsETMl57s18BWzmNHYG8IiI+QGk7yObPkNoaSxqbjfBmKpWI+e4bqd3EE9opAIDVWBF8aS2prOtj7YE6q16qAT0hhb/vPnFIA+bVLQUA8xKmBph7k4ZDGbjqsClky5lt+vtGnQcT/7AkB9uPwBV9zbmS78T8S+YzYlA0H5wAAAABJRU5ErkJggg==\n",
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAALsAAAAgCAYAAABKMQnqAAAHrElEQVR4nO3cf7BVVRXA8c8DBQwNDIKyX0/tBzaZWEpUppkmJaE+K52mwl9hUUQ0kU5Ozbx+GCah6UzUTDnD9IPIMn2COWUZ40jipA6ZWGKO5KROipVIP0DC/ljnDveed+4598e57/HgfmfOnHv32Wef9e5ae++119rn9fT39+vSZV9gVIP1DsBrOilIly5t8mo8L69CkbGPxeW4FtNKEmokMxHPDcPRpZhpuA5XYFxWhSJjX4UxOA0DpYo2MjkdZ6FniI8uxdyId2ObMPpB7Jdz8yyciCm6o0uFWfjocAvRJZd+PIo5WF19IW9kfy8exI6OiTWyGJMczwy3IF1y2YWH8L70hTxjfy3275REI5B3YO1wC9GlIfYT9ltDnrE/X0RhGuXw5GiFFVjT4r1DxRzhF3appR29d4oDMSFdmGfso8S03Si/To69lZfhkRLa+Rx+h614UviVryuh3eFiT9T7WBm2XRSN6UYCgqOxoaS23o7leItwjXbiV3hBSe13qWO3edEY+HsTD+itU74Wf8J2zE3KrsHF+F8T7Q8np0mt7NtgVur7h/E03lriM4aS3uEWIIN/4OB0Yd7IPgW/L+nhH0ye9WYRupuHT5fUdj0uFCHTLbhS49niLN6Iu8oQKoODhGzNDCz7Cq3q8F5MThfWu/l40TOuaEHALB7HQjHCX4ulOm/sD+PrIqK0CO9ssZ2X4DGdyzV8Q7hI6xuoO1XMhld3SJY9jVZ1uExsHTi5ujDL2F8hFHAm/tCikGnWqzWWO3CIiPh0ilvwWVyWfJ9ZUL9exx+UnCiRpWJgeb/GXLrThZzXd0iePY1mdVjhfvSJjnJopTCt4HnihzzVyPQfs6iMmEfl1JmDP8o2+JPkRxsmi0RG0f6W7Wr3bCwTa5iT8OeiPyKhD0/htgbr7y00osM0N4n10U/wCQYvUL8jesUA5uOeBhp9m3BJjhJT/lZsFhGGS5I6bxIr5MroPlO4BlubEL5VHk7Or69z/RgRIXmViJDcXnVtvBhx/5vT/oH4UtX3XpyDu9XmDp6saudqnC22Y9xfIH+FCSJ6s9LuWWAynlAcNduR3J/3dzTDOvFb1eM2nFDSsyjWYRbThT1fhFvJjsaswxL8Jmn8LzkNXoJLRfz5F2IhMVUY0LvsNvZDhGu0HEeKqekrTQhOJJ7OwXnJ50b5QnI+TBjmttT1u5Jjpoi6VBv7LPyyoP3NYj9GhY8kcq4S02ia5fgQzhCL0hcl5dsyZKtmtsh7/KyqrJWOVgbXCxcjzXl4ubCdLFbojA7TvDSR4eMSQ6d+6PEGMTUvxifr1JkqfujbxVSc3kNTvRr+IUbjTjG6XyNW153mFJyffO4RHe2OOnVvEIZ6UVXZbJEEaobpyXlDnevzk3PaNfqi2k6Tpg//UmtkmzXX0coiq+2lwtBXqO2A7dKMDissFgPDj6sL80I5j2JGzvVpwoA3yd4stqXq804sEPvBD8Zn1C7IzsV7cp5FGN0RGl+cHSSmsX/iB0lZ3jQ4IF4AqOzbH4VJwk1ohopfWS9sW28bb39Om+PETHmz/BF6enLe0JCk5dAjZqvF+KYwzF116nZahxVmiAhgjRxFccspOdc2imTI+WLPyNkyAvkl8rgIXT7dYP3LxUjzKfw8Kctb4GzCA8KVIXICdzYpY49QxmPCdSiLU8T0XWQkRR2tbEaLkXy++L0XyA/RdlqHFV6YVdhOomULjsNPhRuzSij4ZryhjXbL4ESRvFqD79mt/KJRYUCE9wijb3bj16EinFq2sfWJ2fOmnDqd6mj12F/ofK6YlS4uuf1WdUgL2wWeU7youU/EiMeIePGFyfdjRWRmu4h0DCXj8V0xelRetHgA/xH+XnVUKM2AmI6niC2iG5t8dsUFuq/J+/IYLUKjt8ofESsdbV2Jz67HODHIzRa/17KS229Hh4Tb3NRGsG1iQdQIO0So8SyxYJ0kFrDDwRKxal8oRjlifbBRGENvzr3rxYy1SOOx72oqSbIyQ6rHi9+zyIXpREfLYryYYU4V0Y6yDZ32dEjY7aCITZ6xP6h+Vu9o2XuYXym2qz6CvxYI1CwvFgodtE+5iuNEAmE1vp+6tiE5502Du5J7F2stqbYpOS8SWy0+0EIbac5M5Cp6B7jRjrZCjIrntiDLBBGKPSG5/1tN3j8UOoR/C/utIc/Y14jFQRYLk8bWizDiEuG/3Sv2EuetyFtlichy9tW5fkAiS/XUV03F5yta4NwoRoZWspT34PN4Vhj8MS20keYM/BZ/K6jXaEer6HxnC7KsFMmku8XI259xjM25f6h02CvjZaA8n/06kRn9GL6dujaQ3DtD+OjjxHSzEl+T0auGgC+L0OFcsepPsyE5F40Kt4goQCvGQCTZLm3x3jTHigRJIzmJSkdbIAz+Svwoo96R4j3avMVuFqOES0XoPSss/YT8EGoRZehwXiLHoL+9p+A/gk0SC4VnRAy12VBcl/b4qohNH2Z3yrwdJoq9NcvUJs/2BqbjArGR8QIZEamilzeeElPO4fbMTfp7O31i6i7D0Il9TM8qb+v2nsREXCUnsFBk7BUeSo4uQ8sRJbe3Wp3/lrUXsLaoQjtJpS5dRhRdY++yz9A19i77DP8H54LIXoglhscAAAAASUVORK5CYII=\n",
       "text/latex": [
        "$\\displaystyle \\left\\{ \\dot{s}^\\mathrm{pl} : \\frac{\\lambda \\sqrt{\\tau^{2}}}{\\tau}, \\  \\dot{z} : \\lambda\\right\\}$"
       ],
@@ -1057,7 +1057,7 @@
        "⎩                             \\tau                         ⎭"
       ]
      },
-     "execution_count": 14,
+     "execution_count": 25,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -1091,7 +1091,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 15,
+   "execution_count": 26,
    "metadata": {
     "slideshow": {
      "slide_type": "fragment"
@@ -1100,7 +1100,7 @@
    "outputs": [
     {
      "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAAARUAAAA/CAYAAADOrlwGAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAOpklEQVR4Ae2d65XcRBOGxz4OYG0iwGQA6whYMsDeCGwygON//udjMsBEwCUDTAQsZABEwH7OYL/36e3SkUbdUksazWikqnM03epr9dvVpeqLNA/u7u52ffTmzZunSvOLrn/kf96X3uMdAUdgXQho3KMo3uv6Tv6PXa172BVJnAp4JedvXbe6XupycgQcge0h8JWa/ELXv9IJn3c1/0GXpaLM3yrzO13v5f+mqyCPcwSOgUAU6GvVhWwelFT2g4MWuLLChM+FmvSnLmYun+n+n1QTH6UCCVOGKzkolL/kd4UCKBkSPlhzTA0/ZJKcfbDahkAxBb6Svxp88v+g65jygUxign939qCeWQOE+UddLH+gWH7T9VmqCV3Tnx9iBp/ypJCLYQL5a3m/kLtahUJT1T7m0Snl8U5xCNjspHqC2S33r9kr8wqSCETsv1fkU/mT1mLSUomdh4mDleIdmIQ3DDSe3q+F0ReZJGsLZl2tQWo7Ftovur7VhbDNSa9VOJbKYBJvoa9iRmQbeqlwlKXTMATeKjkKhf5o9XlSqSghc1bop3vHfzMIMB0wiy6TZPnBGlgMMtqBy6BlAH4S3d7VfuVnze3v6M4ySCOPPB3HWoRYVJWlJT/txYxPmvAKd8ogIOyYBrGeQn9wNdZWckrFNHkjcaaOxQSrcWhPVqkvI1M3cutCzmChbU903Sg9aUcRYCrjZUkZSsPu2Qe5lVCPqvQAmVK8KAxrgzUKBhl8hn6XC56/6yqxxBikP+qa68gB/E1ZR3ml9mBRmVJCeRL2ua6TWeOqG4yf6WIaDcGfrV2EAH6UjgcYaZBn5Bpl3+JbYceStaBUxAdT0oaeyCkVpTs/EqCYYt/LZU+dqVtSaSjclM+URiLg7Nt3kuoyRWaKrjP9nJEFvNgTKLCh9GDJE57FWRuMORbBgrQXuuqKPJe+OJwylRgFPkUpk5fBuChSm8L0QS4yi3JPyqzibNaQnbIp7zFlLdvHDxeF8AGYEbBXsZiuQcAAaGjXEVWz49O7QCl+AP+x3JKn/Qg2yrOM5AX+sco6KZbNk/NFZ8JxkVgVb8dlvc8l/pii1QcCSoZB3HraT6lnTF7xYOc+fk3lVzyytpP7XFe9DY3kMe7ksrYqSyUibJq+b8BjJo4idZ4JQdGTr0sQRjEwIdMIXi5UXakCBg+mP70WXGkTxC/1YylNsVIa1aks+o+pxMkVfWTMHoQNmY1tf600bNsX9YHSZZVOrGt2Z3WWihALHSRwG5aK7pnyBIrAN+ItrtClDp5yJ+/AQn5LkzF1YRAHipixxlKKFYr60NM8BhWWykFIbcHqojyOASyl/65pXB1n+VF87CyydlKkUChjCbRGS4XOaAwCdQpKoCFACpti9rKw1igv15mqBwH+L8az09C7m5Ira2p4IS8vlO5WdTH4ODUZLD+5KJswuOXPHXgD06CYlKYXH6UJC49yG/2lMur0teKzOzSKgyemB/CXom+UJlhOcmkT+FubuN/p/tSDFpmt5FH80B6UXrF1prSUwWIuxDSJNUXah5Vn7S/GKpQy8mdVSkXgoTwgtrkMYISN8KxgkmEgUSYDr5PEA+ZstUof+aNji4Wls4IBkYW8sFCbnLooHCUB3128GyYIczVIUmyqPCxHrBAGdHIaEtOwq5QkxRMHX0y5qBtsGztEShP4kAtPpKc/GIAQbWmkD6FH/BEvJrNBseoeXGhHsaJTHuTxR1207xNdLJjTbsYB2OzkFmNF+im0KqUiIGw9JWhqgBGYCNNvcludpDCEi864ld/y6raXnihFq7x6rlgvT4n64EKATdnVk8/qH8ALwjmFGOAQ+GRJ/IA7yosBwKG53LYuVkbyYRDLoF+rxU35d7rqeNd5YLuc9uFWpPRdSrJKN6PH5O4/8YIyQOnR5ldysdKq9nXwAN5fKq3hv5MfxfSSPPKD9xCsyDaa1qZUgtYXiJVgyc/aR6NjdB+2PUlHvNDrVBAJdBHO20R4K0jlY62gSFibMAFqpSNA8ZTLuRDcUqoUaF+GHC8KR/G+Jr/8QbD7ysrEl2IS+kd1vVU5ZrGEJ6qVqzimADxdk6R4ygjlkED3DJxqUBFWJ8U/rt+X+pVv1j4RH0Fm5WIl2ulUFALtp08asqv7FilfJb+RX/K9lT/gIXcQVq0KBgasTakgWKn5eSWcEXQ6zDqQTv1yIG50Fk+HLKkelBl1MGhM2XHuI2tuK45yk1OBbEUFEX28EK9iGoO6oNhUEsPkNhW5H0Z7dWGx8FTeP5kJTkOwuFZ6FPhBCR5V4BA+htYfZDbiEPLKj+wgx1dyc1Zcqx6l5eGAdZOVsZhpFqyMoYfmmcNV4zi6XQ3oOeqwMlWPafyWYCmu0uRKj0KxhSs69EIXnYcZzolL7g9CKgtTljd6sVAQEuo4WPlDmDwSL9a2jwN446kMVQNBvLKA+6vcIeWQ54aCzoXUvqzMqg2GC1ZHL8WykGN7WO7kt/7Yzz8rVrMpldggNOflfotmurepRcpSCVWKJ/h5JteElU5F4XAwis74Q5d1prxZIk+uw3YqiylEdQ5G/t6pT7amiRFH5mWQpULTxB9Y0mdYK4ZpMN+JLyHlo1+xdKrpUEm+BaTJyqzaAiZgg+VhuCRZVjxKgin9/iI764UNUprZsXrUqPGAN2Ie05bTfTaAD1h6siiA3am+pGApHAXC2kbdzKdTUQDGI26JEkT5dD1BqKthoakO+Ns/1amg2emYvCCw9LvhWdo4FDl8olgYSLyXNaQMsCXfuVGnzKoxyBDYIGuVJSd/RcIJa5vFZtYHcbHUWQCn7Lqs6zbQ7Fg9sprmcAcKxmAWVP6FMqGNg/alAIWhOOpEmktduABft2QIr3cWc+cS4URx8WTgSgn/S8VjijLVol6IwXaKnYZj8gJ+dXxDw/t+hAuWHJgyeMA/NRi6inmmyIYS70p8yji1E3lAZnGR253CkNk/5NanLigT5AcKO2Ry2Q3bl89rhQWLRy4Kg7JJQ9rUA3Z2rGZVKmrYrCTQGNBDBbDOE0qhDvwLRTIIO0l5GATUjVJqDaJYZr3czvLmjDwyL+AxdnCzE8TgYpFyf+B0QqT0U2Sgs+xDR4rXIplVOh529QdekpWYLsTJz05R526R0syO1axKRQ1AwAARqra47m8X8YsAB2tDLk8FzO7OTqlx/bP8PCFaSqWWZjNecFRjMcXBZTCBuy4Uce9AGly4ZzgqAjmlgoBMIgkI5hsLdygWzl6wDlE6YJX0KITWrhZVxXMwIwtrtnb5ILgHDCtv6I5NA2rhP+fWbaMuv5kPgZxSsVX8KSY8T36mEpR1o2txT/T4ZBy1zkFeXUyDSk89CoJVE8p1iFJeNRgbaJxNUcO6UL29OaWCGcvComWs5+n1K9+VEmHt2NrD4hRKbyPKEqA0scKWZoGVcX+gVOpvHiDFr+cfqFov5rQIMPOAeJBUC8wEPOSnTlEhEDRFEVAR6xW2nlKvYjX+2D7WipjqbZLUdp5UvCrQEKxNgrGhRqu/eZAyvjEgGtRSKoq1AcJTeBCpIs4ZsIIftsLw6yreDVDao53AHdSwjsTiGXDhuwVuR7Y1RTHtmX1HYU2Aragt7NjtJPu8ilJR4x8KFRnMWMWyx71/Oq/K1OdRXr63yRH14qeX0jJd+p8u1ip8wa4PZI93BBaAgMYqRgTnY6p/LKwsFUUSgVWBMpiiUDCHoUHTJ9WJKXXy72sGzv3HEXAEihDQuMVKZUOHTyuEsR8sFd3Y6T2+1D1p0VH5UU5Me6q/xizizhM5Ao7A2SJQ0yHPbfeHac/PUxVKROSZ3FG7RkMRFb9Ms5wcAUfggAiMNAiY5WBQVErlU938rsL4KhYr+VOUAtvRU863KHsZjWx8WeGeyhFwBIoQ0DgMsxMlDt8LCmsqCuRMii2O/ik/i6Zj6VIZbQ97bBmezxFwBM4AgX2FAss2/TH2+QIaOzAc6DIlY3G9ripgoQaFNNhSiYrMd396UfYEjsAyEIjjnd0fdmyr11WCpWIsKoIdGLaB+YQdJs1QYuqzU95BOz8xD3Uz7brh3skRcAQWj4CdQWucU9q3VGgFCTnMwo7Q0J0gTtIOzaMs9yRllPxyusW76wg4AstAQGOVWQkHPrFSGmuwDUsFdmMCrAY+z8dUppOUhhfqMIEg1lN+Cj7/cQQcgTUjEGYlamBrVtJSKhEFm4KgJPoIi4ZvkgQFJHe0pdJXkcc7Ao7AYhBgVgJV32K+v20v1Fo4lgrUa6koTVAq0W3MrSjAyRFwBFaJwJPYqtv91qXWVPbTdN7LMhl9pL+zYI90BByBs0QgN/05y8Y4046AI3B6BCZbKqdvgnOwRARkwTIt5vWP3BR60pvwS2yz83SPgCsVl4SDIyCFwrEE1uVYY2POjYKpDkfJv1OawQckyee0fARcqSy/j86KQykLthp5Db7aBZR/p8uVyFn15HhmXamMx85zJhCIyqNSILpHydhuYiKHB60NAV+oXVuPLq8912KJv+J02ggCrlQ20tEnbCbvkNlhyhOy4VUfCwFXKsdCeoP1aOrzVM3mdY9qOrRBGDbXZFcqm+vyozYYK6XxstlRa/fKToKAK5WTwL6ZSvm0qL0ev5lGb72hvvuzdQmYsf2a9vi7YDPiu9Si3VJZas84X47AmSLgSuVMO87ZdgSWioArlaX2jPPlCJwpAjmlchvbY99MONPmOduOgCMwEwL2omjrtHROqfD/P5B93en+zn8dAUfAEbhHwL4K2TrYmFQqWrW3Dy9xzsDJEXAEHIEKAemHK91gqXyQv9hSoQD+qmOnTHxZ38kRcAQcAUOAT1lAjc9Z3AftdklLhUgpEzJwvPqd/Lxp6uQIOAIbR0C6AIWCPuAvTpOvX2SVSsSOfyzkmDV/hYrJ4+QIOAIbRSAqFGYu7+VPWilA8+Du7q4XIhXA//qwvuKfAOxFyxM4AutDQDqAzRsslF4d0GepBHRUIMet+W9lfzksIOI/jsDmEOAdrsfSBbaJkwXg/9DzLOA8voOcAAAAAElFTkSuQmCC\n",
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAARYAAAAhCAYAAAAcRYTuAAAL2klEQVR4nO2dfZgVVR3HP8tbvCiorJCitojAYmXrCkgWBL7RI668KFCPLyGpSCJhkj2QPS0qWoFRPkFmZTypJBUKaFJRaoXyjliShfG4kLqGvGcoBGx/fGe6c2fPmTszd+69Czuf57nP3b0zc865Z875nd/bnFtWW1tLSkpKSpK0yHF8DPAocGIR2mKjooR1p6Q0dT4ItCpynacAvwQG204IEiyTgKuAccDuJFsVgSnA8SWqOyU/TgAGAA0leDUnDgLfBMqKWGc9cAPwZaR8NMImWE4GpgE3A4cK0rTcjAH2A38pUf0p+TEcOB0N+GK/mhO7gFVokheTPcA1wDcwLP42wXIlsBY1OogezitpugATgR8WoOyU4jAU+HWpG9FM+AXQD/hYkevdDWwBavwHbIKlP/CvEAX/3nklzUxgHs1PrT1WaOO8/l3qhjQjZgNzS1DvbiQvsrAJlm5A64I2x0434NPA4hLVn5I/FwLPl7oRzYzVQAfgkiLX2xY4zf+hTbC0dy7IRQXJR21uApYA/0243JTiUQMsLXUjmiE/AKYWuc42SF5kYRMsZZROY/kssKxEdYelUL6lOOTTlvnA08k15f+cDmxLqKxpyN+3D3gHeAr4SEJll5JCjKFlwMUoBF0sWmGQI0Hh5lJ413sDPYEXS1B3FArlW4pDU2oLwLnAxgTLG4z8bRcgE+sQ8DvgpATrKAWFuG9bkfAdnnC5kQlKrPlP0VqRYSBa6UqVN5OSP1cgrSIphvr+vxbYC3wi4XqOFdYAg5BZVAz2I3MoC5vGchzwRohCOwIPAHUoUacB+Eq89gFa7bbkcX0S/BTYjhxhNipoOhnBFZjb8jzwIPBdJKh3o8hByxzlnYfu4+djtuc8YF3Ma8NwPBq3uVIhmjoVFGYMbUT3oFi8gWGumARLa2SShPFzPAbcCmwCvgXMID+nXQ/grRDnLSd39uWdMervSybppxQaW9Jcje7xx4EJwI3AbTmuWY8icvegBSYK3dD9K2SawHfQ5FlVwDrC4o7DIYZjXYFnnePLKJ7pthPojl1p6AocRgpBEjwDVOJznZhMoWuBFc4riErgcuA3wLAEGggamJtDnFeNbO2ZAec8EaP+e5GT8Psxrm2K1AOT0eD+G9ALCZbZOa67D4UvJ6M+8dMCOGL4vIbCmiezkJo/EE2OUlON+na97/OBwEI0iWcAd2Hur0KwB5km5Ujz9jMc3b8nE6pvGfAaMMJbplewtABGISHxmRAFXui8L0qmfYBWyPdznNMDSf+XgNoE6+6FPOo/At5LsNxSsops7WElcDcyYfcFXLcGCaIJ6DkU7ySuQYKpD40ny0VoYSoE9yNtcgjwjwLVEQV3HG4muy9vRxrvXuAytPAWE9c/aTPlRyKt5o8J1XcYyY15aDwsBRq86lJ7YCwSFDsDCroSDVY3y+8hMuZHZZ6NbEPulaiv8742z7r8jEfq3MKAcwYiTWgLEoDb0SQ0reqFptBteRw4Awlbl74oStMTRWm8dED3LmhhKEeDL5cZe4DsPKoHyAiVv8b5MgXAPw47ovsxG9iAtBm/UHmB4O/9B+L3kYtrwn/AcKwTUgieIjPP8q0P4G1kEo3C8eF5NZZ3gdFoxe4KzLEUUo/Uuy+g7RRcc6QBqUR+5gOfA653/g7iPfTlgyiUYLkYdbbNdp+Ovus2NGB2oH7qizKFpyfcniDCtuV8JCxdrWUA8oEEaSsuLzjvl5CZIOuc1wAU/fGay0OB3+Yo8zhkFrhUoLGxnux8mnfICKh5SKiMQA5bN0fjXedVKtxxuA44B20j0BM5zL+Ighl+nkR+GT/XIyH+HPH6yIsrUPYbjg1Di7fXTZBvfSB50M65DjD7WO5EA/bnwJuG4y8i+/sO4FWSNUf2ocftg3Bv6AUYUokdZhHN+doBqELfx3RdV9T5K5C67x805RHqypcobTkVOTvnAR9FT8DeE7IeV3APMhxbjB6bv8Pz2TCUzBZEHdnj5QY0GB/H7veZ6Lz7cz5mkOzYi4o7Ds9GWmIDMgMfDbjG9B1nIaEyn4wvptZzPEwfeWnnvJsWj5FofHuFW12e9VWicZWVlGcSLG8jgTIW+LalsLPRF9gQouJpyOasD3HuVoIFSxkKSYOkvIldRB9w3ZAKZ2tjpXN8M+aVaEfE+vIhSlsec85djQb+j7Fron72ohXqDMOxJSgKWIl8MS2AzpidhUFUOe8bA85pitsgeMfhjWiyDgBeiVjGXCQ456LoqimaVuW8bwxZbld07/b4Pm+LtNllBJurUeu7GlkqWYLMFpJ6EzkzbVQ772EESz0afHtDnLsZ7U5loxcylVZg34ujs++amZjVTy/uNbbEvE2o/eORc2os4XfVqyPaJkVBK17UthxCG3ad4JxzO9k+rHEosmdjF2ZtbDPwd2QOgcLZq3O024T7mP/LMa6NQh3J3gN3HG5Gc6ADMtXC0hJpKBORgJ6EPUQftY9OxeySuBSZPbmiQVHr64VhQQ7KvPVPUC+uYHkpZOVhWQl8CWlDpsiMq35GGYhVIc5367I5p3YAnwS+jjz9NWiCLge+SrCAdZ2rYcmVx5NPW6Jiuw8grWU4mhhXoMTCKJQh38RbyH4vJEnfA3ccrkQa+RoUbXsdaYlBtAYWoN0Za5FJZyNOH30Ys59wJNJwf5VwfSdhCLgECZYgFdSN328MWXlY/uTUW03GeeglrmD5WY5zXBU+SJi+gpzbbZDf4Sbn/37IlDpgue6isA2NQNy2RKEF0nRetxxfgp6k7YJM400Ry++OIimm+5w0Sd8Dr+O2Hml9K4CHgX9iD+W2RU7eYajv7s9RT5w+qqKxcGuJFqBnCbYc4tRnlBNBDyHawr5lSF16jeQ38tmJHjAbYDkeVbCcjFTDI8j5t9+51l9+PZLQvUOUedBp4xg0mDoju7YUBLVlMFKx49Ib3euNluOrkPY0hXh5JW5qQhS/RFPBK1hAY2osmsCLMadddEDawmUooppLqED0PjoTCXp/1vwgNDZymUFx7okbqs7CJlgOYlcde6HnNcKaQaegBucKI7s8jHmD3hZIGh8h/D64roNtKvK4VyM1byHZ2loDWmXKgbMMZZgebz8LPb6/jXDPVSVBMdviCt/nLMePoHyIqcTLtu3ovIcJfTcl3HF4iOwF7hkkZE90/u7iOdYJheI/hfxaYTO7o/bR5Uh4+CObo9D9WpJwfSBZ0UhDtplCO00nO0Rx3ILSw8PmsYAS9GaiSeQVXn2Q82kfwQ86PkTGRq5CX3wkijjhXPsyitd7V9pFKPlvqO/zyU771yB1fztSGV3H5XiKl65dzLZcirTWoMG4FPVZnCxO99GNKchOX0tuk7Up4I7DP9PY//Q9lMsyGfXNEOecBSg9Yg3SKmoN5d5H4zkXtY+uQ1EqPyNQmkiu7Wbj3JMDNI5AWQXLJuBDlmNRBUtUDiOH2AT0KwEu7hObHZHj0sQRslXMc5FNu9XzmSvN/draItTx15G9d+gS1E/9kR+jLRJcC1C6u8kDXyiK1ZZOaDA+jXwGNpYj522cX3LYgHKmJqGBPIejQ7C449D2BPdtSNjXoOjSaDK5QP0x7A+LFohaw+dR+mgA0lb9lkQ/lO8VJs0gzj1pjcF0KrP8EmI12iG/mI9f+1kEfI38UrhfRbkb3kSf0ei7ldN4QkxDyU7VJB/xOpq4FaXRD0IO9ZSmzxMoolrn+/xeNK7PxO6Ij0trp8wh+BY1m49lg3PBNQk3JAoTUAgv1/4hNtojtdS7xWYZSrd+BPMqOwf5Ke4yHGsutEMDcRGpUDlaGIt+AqTOcGwkMv2TFiogrWY5Bk05KCo0DvkbfoJhs9wisAOph1fFvP4cMmnW5yMH5yNILbSZUu87568jeKOnY5kK5Kcq9qbMKfHojkxXm8nSh0w2bVJ0RC6HKuAW0wk2U8jLaSgUm0RuRBxaEc+GvxlJ1OnoeZlyFIa7BT22kJJyLBB3fuRDe5TKsdV2Qpgfky5WKNVG3E570HlBvE2fUlKOBkrxE8j7CRAqEGwKpaSkpMQiFSwpKSmJ8z+AAhKjlqb21QAAAABJRU5ErkJggg==\n",
       "text/latex": [
        "$\\displaystyle \\left\\{ \\dot{f} : \\frac{E_{b} \\left(\\dot{s} - \\dot{s}^\\mathrm{pl}\\right) \\sqrt{\\tau^{2}} - K \\dot{z} \\tau}{\\tau}\\right\\}$"
       ],
@@ -1112,7 +1112,7 @@
        "⎩                                       \\tau                              ⎭"
       ]
      },
-     "execution_count": 15,
+     "execution_count": 26,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -1137,7 +1137,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 16,
+   "execution_count": 27,
    "metadata": {
     "slideshow": {
      "slide_type": "slide"
@@ -1146,7 +1146,7 @@
    "outputs": [
     {
      "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPwAAAA/CAYAAADE16V9AAAACXBIWXMAAA7EAAAOxAGVKw4bAAANp0lEQVR4Ae2d25UUtxaGm1kEAJwIjDPAEIHHGXCJAJyBWbzxxrIzsB0Bx87AOAKwM8AnAo/JYM7/qbXLUpVU9+ruabbW0khV2pK2/tK+SKrquXV9fb0bCq9fv74vml8U/1L+yRC9lzsCjsDhEJBMIsQ/Kb5U/lNfzxd9hZSpgRdKPipeKT5X9OAIOAKnhcA3Yuep4v8krw/6WLvVZ+FV+TtV/l7xJ+W/7WvIyxwBR+B4CEg+76j3PxTxxr/U9V8lbqoCrwqXqvCb4p/Kf1Wq7PeWIyBs0cjPFFGuqwa1fWvVBr2xk0YgziWEnqX3lyVmb5duxns/xtTd+B6QVijCg2Lt9XKFtryJzxgBzSGM8w+C4DulRPJZKK7hRYjVwTWggT+zGn6xGgIR551jvBqk3tBu9yaC8KoERs3C42IS3u4T/7sRAjwULPzkICXBms0eKsqZ8Fz3e3dp92T+91wR4Pkrsn6/r5SYreVrAm8TKCM+V5COMS4eRnwo72b2/73aaDZSlWcJxvqtuHab2YdXu5kIBIEX63jqmQzXBP5GDlOTno0vjigexgF8UJpaPKwignZP8YPooZ0UVIcjyndKG2Gb1MC/xKzZl6zbX4iHXxRNYeApcO+Boi/D/sX5c8ylcz4b/7kJPJsUP2jC8yIC+w9FgdZ9UwwZGEMXqmcKwxTKUJVieWznodIlSoO6KDQPjsBoBM5K4Bm1hOgyjt4sXwkM3kqa7PqqbdZHd0lLjU64hzV+M4G+QyoeGEMaUAAcx7h1T1HxfIbARXZ1Hhdm1XmHoC/gmk8OS4Vd9fESLpX+OrnzSgW1xVrtsaK/L1HByG/vEThHgQ8WXkKQWXhdNy+2KI+FzsoPOCFm78yXeNRY2JPAY/gqjqtE5vccgYDA2bn0GhXWLhNmCQJKIHPDdW+y66s6CNbfigSWBNnHCirHyuL2Z/1DnITHKq8uJ2IffL+AJ1AK34omuPNKEXZ4CF5NvN4pzXZmS434vc8TgbMSeE30YN31KDl/5Os+QnChlVaFLFAN/FF7LBEQrqAoYl8ogLDxpms8CKw3wlZ0rSMNx2fFoHLKUEx8kXilSPvZTr5orH+EHXp4QskR4CWjD3f9jyMQETgrgdeYbP3+pCUYv+m6Y/WioPyselfKW93O5FAZwsW6O/UKEKygVGI7WN3/KPJKY+1oDOtcVDyxDfhs1vbK7xTTPtV8EzhzR5mRNkH0QQE1NzzjCCQInJvABwufCony7Fw3QsTYdX1HEdebozsUQUcZJBg1WdFi5RFyzuEbBUE7ECl9o8Qsffa7ASrDTe+z7rTRCLfosdpY+2JQ+d1igd90BHoQuOgpW1ykSflRsTrJF3fQbaCzfo8kDQ/iB6uI8FlASby1i1KqOigEzvihpS3GhbudBd1DQLH0rNPxCtKAR9A+SkvL2/lnujF00tCu49eOQC8Cmwm8JjyCxaRf9JJKL/dJofoL1l23OkKistSCI+y26YWCgE/cdVxx3lzjuhN0n7Uyn5ti2dmUg75Ea4qgWUuLjs28X5V+Ujo2UMdfrBmLltONQmBLgWdy85JKcQNrFHfTiMzFru6QixcU0KNE8FASKAN+4AML/l7RBFbZ4KbzznpzZq985s4HouSPymkPHl4obwrhla5x90cF1YNPNh4bF39URSdyBAYQ2Ezg6VcTdopFG2B1sBiLSJ9FIdF9hJsNrsa9Vx4lgUAbn6Rtj8TceBXtg+jpCyVh9azIUlMaCH2w1D20VidNqZN6JWmZ5x2B2Qjcnl3zBCpKiO6IDXbZg0WEJd2z4zguCdAgxKRs0qUeAPcb11t5vJG2oD3XvUvVS8/G2fCr7obThyKKB8tOe9kGnq6HwiMRpIppiN7LHYFRCGwq8Jr0TFqzgm90bflRzA0RxfamClPaLLv1qUfwVIUIeBNieUrTlA1kcOFRPpwStJVIb1XRLxlTb9uHKtQYOK3Ag0KpEtiPSJ8/ChhFfY8y0UM7KoiWuij6nfIng5V4Ycwoazw0AsYF45DxqGvmBTTgAS7N+x3K71S+2fhqAk+Hi4KYxq3lYSL0vyuyPs6Ox3R97IAwIvQ8FCw4E28VHmlHEUWRehDHHu/B+tfYN/tyUW3zvFDM/yhlSXYSGIsPxrxTyteazK2aEnsLnULxB0tUb7Px1dbwCCphjmXb19wffTGwoMGVpq600Rw7RfMyYdDMCH7tAc3iU+3xfvsk6z6roxOtpLGz/0Hoe/acmEzGSG1jHVHO6RGrLo8bxNeDyEHRcKg88KuUl8MYQzHEsrnjMzzxoLJQs/AwjZaxilmloQvV40HjJbCWZVB9D3youc3KxRsKrboW36zjz6dhU6Cdo9IWBM0pSOv+0CXeI+88EIsCNtTABuWm5LIxiz/k4ZXij8qPlau548ObJoB/8DrClf50LLyYMYaXCCkd4dJUNZgx4OlZIxDmkuZBNpd0jUcVQpwjWbmVDaWxXebYKSntZ/AdeQtDUB4D+kopa/Wxwr4TLbhMHp/qofyoZ7Ks7D50BF637UjpuRGNTdXRC0U2JMxt4UUWtNSoINpDv5k3ii8nmo0AEz0TZj1jJmFmCHRvydKRJcGl2rgzm8t1KzLmZjziC1ngG4q5+wxzx8em8U79NsqV69v8sRCZg2EYzB6K0fSlqgNznE+zacGufOZODNTlgbHmsF3dPnIvO3EE9OzNuqz+5WJr6BgUJjWCNXq+tdpY5TIZc1ByuoYvDOhoq15gZNb4kD1FTgzYo2IDOfDQWHjd4JiAxnE7ENxZQXURWkKm2fe36n9VDwVzyDfz6sx4yRoI2PqdzakQ1Siud/GYUjR8YfiHYrb2HcHIVaQ5Bbfexvy3xoGgoYCQJZSeHdWNGFJGMnt86pNNabwNvsIMchksvC7QQmhIHszSzQ88hJ3aadwarscE1ZnsVYxp12m6CAhrPCqOS0nHBubH2OcaLHxKrzzCns0vXXM6wgbxpC8XYZi6SlhC4i4jYLXPklW8D7HOVuM2r4bxmLdhsvVKHGRjN55q6ZzxtdtSG5wUwQPK9EkQeF0g7P/VjUkMtRuP17gRS1yYSrPd2+KXpYOHCgLCp/q/5VSGct3yO4fO+j2yiRcZgnhAYFNXHIH5el/a/1d1sVgIu9EzqbHyvZZ+43GHMauPxkNWHiWHt3updFAhiS4E0c4an9VvpWCOh9EI/Be6+F2d8K45WnyJwDLosVZApPOD+KxO6Pmtes2lCOi5mKXruOetuYWwB+HQfeYNCgDBoD6Go/ZiCrQIO3MVxbVTirF6qtgr8NBuESLPNN0Zs+6hjBgTVh43uzeordXGp7YQdLBiTf/ygp6VwQUxbY/pB/i54aEqvp9b2eudBQK2lq3u42iOYcEeMffiiBEIDE31y0XoRI8w4JKzsZwaFqwYy4O5a2VVXxSqYxZP4MDYeF+gV7ZUvtr4IhaNsDO6IPBkYsA9giEAnRzUAQ+R+umDGNUOQCheK2Y/2TSqshOdGgJB6PQsi/NA9xHuyV8uJvWw/JkyidcI1VEsvPrtHbPKbSmDlS+GNcentpBFhJ29keZI8Hbaswpss4Efd5jz9hLaaae62cNI+6jlY988sA81Gr9/ugjo+aHof1ZkohGZB0y4NECDBxiMgsrTecL9ZmIqX/pykXIse22viXLe/WBXnLm0aVAfNmbSdMzvVWabdjvlbeMOfsJvHiplHG0e1xyfKZhsCXHr+jrf9xITMP5Rkc2G4g8uqqwYRE8n95RmnRSJ/aYjkCCgOYN31+zJKP+PirHkNeFOans2RUCYmQxj3W2pHkguUkLyIkDrfFJES6K5eoNo8ARMk6Ol3/ZW8EJHoIwABibMN6Vs5q325WK5u7O+GzxtjTD1oMKAOwIfYTC3GgEeCrgrrL/tYblGHkLMy0sI4BVu9uViqcMzvmcbiHjqWcjW8EkJFp4waOFFEwQ+pu7Kg5qHyQjIYLDBd6wNt8n8nniFe5G/qzafNYFv01Wv9aCalwyqRF7gCDgCJ4FAzaU/CeacCUfAEVgXgcUWfl12blZr8m7suKW29OHoxT2gm/VYz5pbF/iZj1eCzBEkex3sW7BWQvjTc+SdaIovnojOgyNwFARc4GfALkHm2GPKP36c0YtXcQTWR8AFfgam0XI31jsqADvZmNGiV3EEDoOAb9qtg/MzNVP6Smqd1r0VR2AlBFzg1wHysZqxl5XWadFbcQQ2QMAFfiGocufvqwleQ25c/IVNenVHYDMEXOCXQ4t1b3/1tLxVb8ER2AABF/jloD5SE/Yp4vLWvAVHYEMEfJd+Ibhy5f37gYUYevXDIeAW/nBYe0+OwNERcIE/+iNwBhyBwyHgAn84rL0nR+DoCNQEnnfDCfZd7f7K/zoCjsBNQMA+5uq8/VkTePvlWPvljJswSOfREXAE9gjYL1V1XgYrCrx2nu2TTs6YPTgCjsANQUCyeylWsfDvlB9t4Rle+JldVcr+3SwFHhwBR+BkEeAzbUL2qfb+VvcfUdj9nQSdCrwuGv5JX1PgGUfAEThJBCSzCDufbvNvpYqvehdd+mQ0/CcaXhvl30/hKnhwBByBE0QgCjveOP+qq2jdYbvzjyhKY1ED/O4863n/yaYSQH7PETgiApJPNtmx7IPyOWThwzDUIK+Plv71Tyj3P46AI3BUBPiW467k1Dbbq8z8H0IZxowXBNlgAAAAAElFTkSuQmCC\n",
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPQAAAAhCAYAAADu3IZiAAAJMklEQVR4nO2df7AWVRnHPxfxiiJXDQQtdW5IzoVRvIAkmjiCI5R4xasif2Q6Mv4Yy4iKdCwr0LRSi8bSaUibBqJyKvOKinaFaRz8BQhXk0zMBDTvJD8mqIwQuf3xPdu+73J29+y+u/e99333M/PO3nff3XPOs2ef5zznec7ubViwYAEFBQW1wYCY3z8EXNkbDSkoKIhlBHA+EXobpdBnAE8Ab2fcqFrgSGAS0FOFT0H9sgtoAZ5Cg+0BDAw58TBgGTAN+EsuTevfzATeAxqq3ZCCumIP8H3gDWAxcGnwgLARejawjnhlPtF86o3pyHspKKgGHcBEYGTwh7AReirwmkPBK822OVWz+ieN5vPPajekoK7ZCkwB/lq6M2yEHkXhToYxFfhDtRtRUPcMRHpaRphCNwGHOBTaTH2NzgBtwCPVbkRB3TMIGBLcGabQB5kTapFK5/3HI3cnC24G1gK7gW3AcuDkjMouqG0GYpkyR6WtatXlXok/90/KOKAru6ZwDnAfcCZy5fcRkZIoKCjBqp9hQTGAvTk1pD9zIRpFs2J64PtnUK7xExnXU1B7vG/bGTZCNwLdDoU2AfcAm5EB6AFuStG43qSZ9PP+CSidlxdDUJ/szLGOWuFadL9tBxYRv+qxr5JWjm6kp2XYTm4ETgCecyh0GfB5YCNwJ7CQygNGncSvlrqlwjrS8BHgHfJdrfUD5NI/73DsCOADZFD7Mnn155vA3cDBwDzgvAzaWimerFMsv40AVpnfV+BPq9LK8QyWKLfN5Z6JRqHVMQW2ABcATwIzHBvhwng0l7w94piHMqwvyABgv2V/G/m6wXcBZwOTkaLGMRO19Xc5tikL8urPTvPZCdyBluI+maKcLBmPFPbFwP7JwINIqRcCt+LfY2nl+ClwAxpo/ubtDCr0KcD1wCyHAqea7W8djnXlRGS5NgALMizXlTZkLUdzoFKfi+a4efA94HJk2V2X2rYDO4Cnc2pTFvRGf3rezKk5le+KJ+smlLXw+DLwHRQbOZ9wZU0qRzfwaeBe4Grkspe53EOBJcCNwFsRBV2CrNC95vtifNepxbExYZxmtmsrLMfGZDQSvIHWxL4LrEFW0av7HOBjKOpcymA0au6JKH8YMgJx7uV/KU8J3oOvzH9ylOUIZFCX44/maevPkzz70+NNsx2bYx0uBGVtQvfb3cB6NHpHjbxp5FiFliB7ulg2Qu9AN/Kj6CYPS+10I7fhs8BR+K5UD/C65fifoUcwrzJ/R5HXDfBV1M6t6KJuR+7PacAnze/rzGcSimaXTjmmA7+PqeNw5Ep5NCO5X0TX1GMbvmG4DynzRcjlOsbs/5f5hDEDxTpKXdU09edNbyj01812JLoGUdctTzxZ1yGl/A0aHH4MfIH4rFEaOb4BHIqevQAOdLn/A3wFaf3xyJoHeRZ4AY3kr5KtK+VdlDOB40KOuQv4d4IyR6AbfTVym4MXdljg+8PIhbmxZN8MtAgkis2UX4urkUL9CllpG9ebbdB4LiT6uraja9BZYf15k0d/ljINmGP+bkBTRpdgbh54so5BA2IPmqL93OHcNHJMAD6H9PT/2IJiG9DjkxcQPj8egyzDeofG3ozmEHFpsAa0cAM0mtvYSXID0oJWvm3CbiW3B753oIh9C/BnNC0Zilz0JLSabVfEMWkW7wxCXsUKokdal/rzJK/+9BgC/AT4B/JALkcjYzUUulTWa5CBmgS84nBuWjmuQt5A2T1tS1v1oKjZmIjCxputi0J3I8XYFXPcSWhuuBpdINtnaOCc2ykfpWxsNHXPQSm12WiqEMYm9KTZheb7GcgjSYoX3HgpxblRTEMuWVx0O039m0n2soWo0SdNf4Jbn4KM7gnInX3c7HMNKG0mOznBl3UT0onBaBrlQlo5xmB5+UjYSrEe4OiIwjyF3uBQsSuey5LkBmx1OH47cBbwTRRlbEOBpE7ga9iNUgdKC92JFHtJgjaBbtaxKG+9LeG5cbQjq/xYDvV7AUNX3on4LU1/glufTgGuQyPaEvzBxzWglKWc4Mv6HPJI1wC3oUDXsojzKpFjGDIgZUQt/YxaseLl27ocKnYlrUL/0uG4V1AqrhHleq813yeiPF4wVtABzAeGo4u8MUGbAD6KopzPJDwvjoOQQVpFtMeTtv5zU7bLRiUKHdWng4H7kfzXmX2vofjPKciYxS3+yVJOKA+IdaPp6mqUK34Le2qxUjms+hmltGGLGxqQS/A62T7kn/QGOBr4MErVrESvBHoJzV3C2IsefrgMXfChKGgW5Hk0ss8j3SuYvPSdyxwqCWejNse523nVn4Q0Cu3Sp99GkeC5+CPnB8joNlGdx3lLFRrU5tnIAD+MPZ1bqRxeirKMMIXeS7hLchKayLu628cigY6IOGYAssz7gT86lusFIeajKPZ4dGEexPc8xmF/VHIUekxxK/aXIO5HOd75pFsd1mS2uyOPSs7FqG0dVarflTT9CfF9ehaK7C4HlgbO7TLb3s5He7Luo9x4PY4GhKPM38NLfstCjvexBHnDXO4d2FNWkCwgBrJEcXno0SjQs5vohzsW41uzViRQO7DF7LsJXdRmNLLONXWvQZbvXeSOegGvOdiXeYICaJeQbiWWN7eZh1YPrcVtahDHRSht+Pcq1e9Kmv6E6D4dDTxAuYtaiqdMpxJv8LLEk/Vl5C6X8iOUi56L7idvjXcWcuxBelpGmEJvxPIkhyGpQrswwWybUPDKxn60RNJjHErebynZ5+UzPc+jA8n4cTRnHoRuoF8A38W+EMajEwXF9jlJUM569MDBDUipFlG5Qk1EudxFVao/CWn6E6L7dCHyDq/AngLtMtveHqE9WcOewvsiGkTaULR8C9nI0YglttMQ8p8zzgO+BHwqptBq8iqydKWLJmahnN4w0iliX+YOFEEdib9MsNaotz5Ny5FoMBpFIDgaNofuxI+o9kUOQ67MwSX7GlAubym12fHtyB2rVWWuxz5Ny63AD7FkOqKi3LPM54GcGlUJY/GX1p2OLNVS5JKGuXj9ndH4q79qkXrs06QMR3GoPcC3bAdE5aF3IT+/L77fqhUtDrgF+DVyx1ag9Ebxto/+SStFn8axG71QJDRdHDaHLigo6If01/cwFRQUWCgUuqCghvgfa8yEuTrOptoAAAAASUVORK5CYII=\n",
       "text/latex": [
        "$\\displaystyle \\left\\{ \\dot{f} : \\frac{E_{b} \\dot{s} \\sqrt{\\tau^{2}}}{\\tau} - E_{b} \\lambda - K \\lambda\\right\\}$"
       ],
@@ -1158,7 +1158,7 @@
        "⎩                  \\tau                                   ⎭"
       ]
      },
-     "execution_count": 16,
+     "execution_count": 27,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -1181,7 +1181,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 17,
+   "execution_count": 28,
    "metadata": {
     "slideshow": {
      "slide_type": "fragment"
@@ -1190,7 +1190,7 @@
    "outputs": [
     {
      "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAGQAAAA0CAYAAAB8bJ2jAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAIeUlEQVR4Ae2b63EUORDHx5QDML4MfBkYOwJMBjwiwGRgyp/gGwUZABFQkAFcBHBkYDKAcwbc/yfUsjSjeWhWy06x21WyNFKr1W9pR+O9nz9/NtsMz58/P5b8/9bWgejuzaG5P2fSHzbnpeS5IwV+XYJct5bAxKZ48NHRLMUY6GHbI+RSOiBCikAGPNAE5gJHv6rmsfqvfXt2tbUGkfJQ5JHqTzO091Lzntg8tV+rzT70t/XNrasZRExdiIl7KieemS+qY4/Bq1DCocoX4YNbBJpzpQmfVAdlFBFIkZ/qkTIHzsXDexUzJlFG37HKSntRTYO8ElOvxBDHtq+qswpXvxmuSBGaZwY1gxfNj5E9rRPVcw3LPByuOlQzCJxJwDPPoXlOjuE36iwObdG+VrlNnSNa2IdHvyicE9DFAzLEgIG+qX+l6IBg7VOWRcXHmNtMm9RTDDWMIRpE2pnqD8UMZCaIDr9j7qvcyQwXd9U2iIsQMZlEiJ5JUw7UxsOTcT/0u6pZJ6scc5KFPdF+x9SI3OrHXrwlUbaYxkgJs+orDm3NQfDvKgAp76n62nTxVFJbwgMTIriv8WzK9GucC/cgwo+bT4Tj0pVqjAEPLiv450b1t3hCabvaHiJGbP/gKPneM4Jg9GcVMJVZ0SMFIrwzpF8LA4VNWX1EId6PQrLpw+NwRO2AxujHwA9UfqhAPzmFCcfWxxjgwxNOCMBLgu96C/9UM4jWtf3jQYvxj3rueI0X5K3m/VDb5nbY1xjCk/PjqEJwM3qjMZSC5/6lcsFzC1/dDvDwjnP4+fAZ9hW1G5V4TaNBzW8OnC15Byb84CAgzYGaBiESEiHEICePIKQfP1AfaYWjMYbqGAu8NgiXKMEI/A5JDAgt8FW/UGWRgqcH0BipqC86mB+UL1wMTLRkQeO3swMVOm9VoGEkOvuHHwhKkCB4FYoxwIjv7CFXaw4G4zcOuNC6Uh/ppAPqR4lECvsEkRUDUdU+rsbjcfuRHsZOijF+tXYVg0h4Fx05ITQWRwDGsE0RA2Ig0hFphl++PHdA/eRqXmcTGWzY4GdxNWbGCvlcuGz2H1RjsCkA/lp++I0tXsUgWsRSSO/pRsrAY08jpWBEjPVGfUTAZxVTppouBfHOKPxmUbuTrhxi9Ec40ISPc7XNaJd6Jp2NgubAJweTkMJGJ1VEqGUQPKrpE0L9KJ8NMKQvtTEiCjevpT5RicHSVOgTPmthRJsXxqKGGRajOG8fwY+muh95cVTHY2tv789dQQLifZySnEdBR33h5MOzAByUTM0mHkcQ/SGtqM1Rta2Ix+ojpcW/DTgQDJ5mWEcFDycyoJls8HoeglMNxo4zhFt9bBWD4KElgraZ57QVp4WHQsAAAfx4jBPGJjRIUTgIJ722oXunC3cVmXrpTh2olbKmrhfjoSgip1FNBPBKPjkix8ilbU8LY8ZRWErmt+PvbeojBymMUxapx23aemZj33rYmEG2XvM9Cthkyuphabu7dwZZmP13BtkZZGEaWBg7e8+ePdvub0mXZpBNHXsXpofFsLPbQxZjil+M7AyyDQbRr+6jhcm5KHaG9FM9QrQYV6i8FtlBvwa4b0FPHaj66kSLcPfAJVT2hZ5ngnsQXr0D3Mrx1tiAl41E16HKrO9/jdC6ai8Dr+iRFeBKgSuB5C2xnnnTDA7yIWf4akbtxtNhXnKtXM0gIowy/1Gd/QQHJgyEY9//ZnE1jvfcU43xFgleBt5Yd75igWH1Y4xHKr3/piAcLu3uqg5OOfs+RITawC3d6MWOFucWEIgvq3713PzFa7KC3qCUt/zapIvEK0spab6l5Ox1gca5TuDzpiRqMuugL/QWLtxq7iEPJwpqXj/2VUe4S88IMreLKKasCuZUiQySn0s3FMyNZdZY8cJeX+gt8FTFICJIeE69lXPCaE4SIXoOm5zahHAyHguygDapqIllUJuouVTNXjFVF5AB9yENIKQsEcGy8d21Q4j+hO9aoz5r4vVTFQjjCa7WxkgYIYD65l7dBhprbCBD4E+8ojf+cTSknoK10QX6c2nUGUSEyGUohJz3QwXjJCelEQWdCL9k/6j+/a/W/y0gPbgI12LOqfRMZKOvkqiIeSU1G81mXwSxdsl3rTExa5MDMeQY2P5R9fvfsUUrj5sM36U7PmMiRXEAOVfNF5Oje0eLH/TGUd8BBiH0KA70jIGS9OGHhqrDiXOcJ/g1HT21q33/awyKJtEavM76VcNno/FcauEzpbFTEdONLr8h7DsAS/eXGi81CJGFQzvYt0ZUs2Elp4dobNUmxk72D08wpDsJCXPkZBMWBdz1eJOqHoU36ufwQbo02pPotZCcDKIRjs5q41TIdaa678v7FpnwmDhzziAwPcVTAkU1CLtg5XjA2mLUPKtjbAQyPNUYwwmLcGpDF0GZf6rS+0NLY2uFIRm0MFECj0RJif6QD/05SI69WpBchgeFFObxxioUytwhsNybixA3z69f9P3v0IJrGOuVQbwjF3pgHxl0zhZfREhwyMQgEIsHWxOHHjEg3jsE0G7EbNbY6se75nz/O7Rm7bFBGbSYpV6iZCrw+ijoZL81C6Ua0dbQ4OM7jba/60X5eMpbFRd5UFBfGw+cExXqOd//atr6IJIB/lwW8DJ8Vh32IrVtY4cZ/l2CdMtvt+D9DGQAR7w5ZHCFW6Pobv5K5bgGLaPBfb+1qfX8n8r9uK+kzVyVi5I568QVL0cqV/Ea7ZSVMeDkLjzkxtKTpw0i1v7+91qrUZYC/PhGbwGqvX6HosKTE9SUMA0MDDUs7IXzx33/K9lIf69V20HBqaK9hwzpZ8oYxz32iGSRKRNzOGKWza521OWW2kQfe3VHtqoRglTe8hz9woa3CWmXvKZ0w/sv/uexs+H/D8zCQpWNuobfAAAAAElFTkSuQmCC\n",
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAEAAAAAsCAYAAADVX77/AAAFJ0lEQVR4nO3aaYxdYxgH8N8omdopEgRFxZZ0GDM+WCKxK4qILYiSaoIQiZ36MBFLkAhFbEkJQWqtRlCisVVqa9HEhKhdq/Z9KVofnnM6556ee++5vXfuTFL/5M25533e93mf83/X8z+3o6+vz6qM1Yc6gAKMxB+D5Lsjn7HaIDXUDA5Arwi21WkFDEcCdsPcdjU23AjowLIktQXDjYBevNlgnS3xAt7DOzi6kcrtIGAr/IruEmX3xawG/f+Dc7EzDsRNWKts5XbsAguxKz4vUXYk/mrQ/6Ikwdf4ARvjszKVy4yAWQbmZTY9WTLAf/Ch+g82Gp+U9FkNvVhDObJRjoBuTMZmuXTiSgSYYj2cnss7FE814XMj3IuJGlhE6xEwBhvgRXyVSz/XqduBi/A+/hTD85HEdgyuRGem/Cb4tmzgOXTicVyDVxupWI+AHvyLeSsR1IU4DWdhBxyB5xLbVDyL/ZP7dfFTFT/TFU/BZThFEH2PmKr3NRpkGQJGiN77NZOmJfbDcVuVuoeIIf08PsUc3JGxP4Ejk9/jMLOKn4liym2NpTjWwDR8EHvheByFt5M0ts5zLUe9XaAHj+KSXP6PybVL9dHxOG7ELsnvaSqH+DO4XvTgDnioip/vkmu36LCXRIekeEUT23m9it2YLVbxbEofpAvbi8NLv0rmb05sT2ESFmCnjP03cXjZQ+wU9dAl1p6v6xVsBLUI2Aaj1J7/XWKf78W1YtHLYgFuMPBy05WzP5HUm10i1i5x0mspak2BnuT6FTbN2b5J6q4phjm8K+YiXIzFeF307gQsEUfWLGbgKuUI2FbsKC1FGQL6c/nLsKHYIvvFwkRMl7eT351i3RgthvocseIvzvlaLKbHvyViXS3xt4UYdUtrFy+HjiYUoVPFQ47F+nga48WIGQzsh7uxOdbR+JG5EM28C4wVO8QbiZ8LDN7DE/v86FY7bYaA85Pr5FYEMlQYbnpA2/E/AUMdwFAjXQPapsENN6QEFErGqwIanQIbisPLmEGIZTDxCM4rMjRKwGXiLW5BJq9ZyawVmIX7c3mT8Ls4rF2By8WBrQKNnAPWSpyOz+WnktnUXP7vDfguwj1CI+wrUbbbAOGduFXoA0caEGE+xsmJbTmyBEw3IFDkMUF8r1sq3r9T5CWzoUAaw1viG8Gj4rl6VYqsM4SOWUFAdgrUU172ThrJ7hjNSGatQo+IaZT4pNYvVKJPcuVew+7iDXY5sgR8J3pxlAHlJRVA/xbELCxovJZk1g70iA57WGgL6WjNY6GQzDfPZhatAdWUl5FC3c03XksyI3TDw3BmlQdIcVmSUnSKnr0gkzcOLxfE8Dx2FB9WqyElpWIEVCOgSHn5VmyDWXQLefvDGg3X0g2zuF2lLngtvsSUTN6XBfW6cbVQoWeLRbOvoNyo5PpNNrOIgGrKyzyhAaQoI5kRBHwhdMO1cRzmF5T7Pkkpfknua5GbxjBXKFIn4TF8gAdyZceKaVAhyhSdA7LKS9Y+U4iaGyf3ecksm0Zk6tXTDZtBGkPaCTPENJqKPXNl9xFnmAoUEXCTGFYfi0UjxXyxkp6Qa7zfwAfKReJh10lsnVbUDVMCW4EefKRyzblO7FrTxQiRxHAU7so7aFQSO1jI3Tspp+PtJtaIQ5P7idgOlzbSaAtwtvgydVDe0KgiNBO3iOnxaYnyXWJNWUMcQ8+w4kmyHViCc4oMKyOJTalfZDnarRtWw53VDIP9B4lhrxuu8orQ/wQMdQBDjf8AH/wvoAr7RYcAAAAASUVORK5CYII=\n",
       "text/latex": [
        "$\\displaystyle \\frac{E_{b} \\dot{s} \\sqrt{\\tau^{2}}}{\\tau \\left(E_{b} + K\\right)}$"
       ],
@@ -1202,7 +1202,7 @@
        "    \\tau⋅(E_b + K)    "
       ]
      },
-     "execution_count": 17,
+     "execution_count": 28,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -1233,7 +1233,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 18,
+   "execution_count": 29,
    "metadata": {
     "slideshow": {
      "slide_type": "fragment"
@@ -1242,7 +1242,7 @@
    "outputs": [
     {
      "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAAASQAAAA/CAYAAABJn7UAAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAQfklEQVR4Ae2d7ZUctRKGx3s2APDNwGQA6wi8ZIDtCLxkAMf//M8HMgAiMJABJgIbMoAbgfc6A9/30UpaaVqt/pjume7pqnN6pZZKaulVdalUUs8++PTp066LXr169Ug8v+n6V/GnXfyWbwgYAoZAQEA6AyXzs67vFf8Y0kvhRSkxTVMFN7r/R9etrhdpnsUNAUPAEOiBwNfieabrv9InX9b4H9QsJBX+ToV/0PWz4t/WKrI8Q8AQMATaEJD++Ex5f+litfWF7v8t8bYqJBW4VoE/dP2t+FelwpZmCBgCp0VA7yYWBy/6pKR6H0xaoSpL2orr54tS/ZelRJ/2kw9tmVYBybIMgRMjwArmK73gf5+4HZ2Pp426fhTjdwq5iGdU9CGJEa2LaUUFi+9o1iO7MQQ2goB/T3cre0df++F5WRqmNgvpuWd+UypkaYaAIbAIBHipsZAGkRQY/pygEDA8oBdKr+6A3bEd9pdn6MJ/9EghV+ZLalNIoZEZ82FNsdKGgCEwFQK8zKqLF/rtiDp/ULm4SaU47hn8UEW/zoj6u4o4hSQmVmKZjmlTSF0VWn6CgAaU3Ui2Nq988nuF6WzDjIQAPdT1Xvzw9ibxc+zircIoRL0LG2MvBOYeQxox8Th+ryq5xtCN2vKbrqDMsLJI+1LXMVw06buRtd8UUgbHuBsNIs65HxVyAAy/W1HhKD0ort4PUpmgzIKy613WGPsjIJxnG0NaMeU4+rquFI6doCjHpLk4ulhci1baIAnHtW96mHVKPeG0amailpjSNNXLbPK5Qjt6kQIzQ3yuMaSpE48jFs3rsRCoLZwrTK0UFBRb8cewjqrNNoVUhWdQZrCKOLtVI5Zfg2hPeAaVNeZBCMw2hrRiinFUHVjM1wp/H9SzFmbVgx/nG12LmPBMIbUM1IhkZyFpgDMLSfcs0xwpzqyU5fssC5aBwBrGcNTOWgleyeMjpWNtcY4ptZhK7EdJMx/SdDAz02TKRoOMgGcDrbRBZrH4EZgPvpnsgnR+oOh5LRiOwCxjSDP6jKN4sFTYFs/kaK8b3yi/uBvmn8G3p1hRJfpWPLgNdgpRRsiSswr9PemDXAqlhxySZgrpEPR8WQ2im1l1yzYsv4oAIRSkF4UHhi5SXSz/EBqnxPxzUFBjnZldj9xsvseW/k86hlSoujvHUTxY0lg/KITi8snzsEXfIOWRzuTHr3Hc6kJOsl048QQ5QhnBj2yhhCFkKuN3qUf+YwppGsCD7+Hp3qD/ofvGjOOF4Bc9+lbxUDZridIRGnwFqUWFwASFl/HbzcEIhHGYbAxpUZ9xFA9KAcvlP7r4pKJt+x0LpzHB+fLIWvQrKb7TlcqOqo7EmSMmzOwbOPGffKIzhRTH6KAIllAmABpcdi2igPj8z5SGSc7RABRVQ1nBl5L4mF1RQpxDCi9NymLxaRCYbQxpXm0clRcsl9diDZZS9rtj4mEp1mYdUT4qH/Gi4LCWiqT8z4sZC0i8WEAbzqEJDd+D71QUIAkBMxJCFYgX4E242Q/Fj7LibAx81POP0jDDjeZBYPIxpJlDxlG8KBEsJfxEWMgpYR07/0+a2BJ/rnQmstXRYhSSBgDrAbOz+ycsFwSz2utmVjWpIQDKSy0glFFwKCL8KKhr8WCic2qW+4yUxhqfn4HAMsLRCW+DLytkN4MREKazjSGNUf1DxjFMOtGfo/I4u39XiMLqQ/Av8uBjV+OXpJAA++Rr2C7ACvlhGdW6MyJBYrZ7nAgULwDKigNqWEHvdAVB3CmNb43ieSXFbbkmgGakyceQto4ZR5VBLpClG8XD5PNS9yznOkllkDUc83EJ11loQQyLUUgek9sFYdO3KcxGuzYBUDrKB+dhXL4pzguA0gkzHuGVrkCUSfmpn+fsn7AN/BYehsAcY0iLxo5jmJxQSrSN7x+DrHT1FP7UMu/iX1T+5TFaIzDR2rxghICN5mdHgRBzti/YYj89qb20m10yNxvRIqXt737Bc6WLECd2akGRHk1yxdnmTYXohe6vVSY9U4IzfI0WpLqyPBKWc48hnR41jsiKLiwcLCPkInNw675Gj5WZTWY15qXlTa6QBGTjy3SlsePEC4ilANju5VPIjsKfuornLpS+SFK7P6phQ4Rkvx/4y1KT+pkYEF5HPi/ND1mrCtWPhiwspQNzjyH9PHAcWaIxyfHupJMVVbeSeA+Ry9Z6j5Vx0fIgZo/BJDAoh9WABVAiZvkIruL4TzhzgWm7JULIHMYKsYIwybMjAmsHw/evJgtr7+KsY+jlgUkptaTXjlln+9sU0kNfctAsLRCxHIZ+mU4ZBHdLxCyGDwkLEWsJn9JZ0UhZWBMGs4+hMOQbsziBrwmcjraGPjXe+7YlG9vSmTXT8YCY7QUx3veIYCmEBvZgXz+LMELRn70/aIQsrGZwtzKGMw0Iu8oQEzGrpEgNC0lAh+VT6oSNBQ6MYA24pQr1KI6FgE9pjmcd2FQrbggYAnMgoPcd9wQro6Br4mNKFlLYcoxO1sjdEdGDKPvBs7V9mf5MfLfiwVzjH8a55YpCFJV7tuI/6eq0IMSzWKepx2CTgcYFCzvsOrpvw5TGeLNzWD1tbGO6GZHBae/cFhrzaCVl/yhSGThY2TKMP1PQFx6V5aQyW/jO76QQ7YcwOsWiECH9U+Ek39GoHhTY/3Sxpb6qXTq1+WzJjws7p290cbSDcUcmOKxX3QGyMRVKGyKNN5MW56bif7KNSzZlkoEyQqlUZzHxZCT+MPs5ZeQzwzZ/yhuXa2nimLieick31IE+5lFWZhgCD8X+ROPDb4wjS8Fh32lx25gOA3rt3BpvJih0Bp+MoUN2zkLSDUslrCP+N9Pg7WdfmVs+qQ60XjxrpPjO5/MMlJ4TVNKNzhsBjTsTEIf7XivOBGJkCDQQkGwE/fM0KCSWPr8qo9Nv06jNJ/hKcVIHWpziURtX9eFuAHLOUJhM/j/caa/qZcbjq/XoH5ijH6rfxnQOYEfWqfEYLE9eVnD5vA0KiZmMdT+E32f0NrzK4jtiuUbIkspmRgGxJfIygM8oLv0VZ4fVZGFLgtCjr5IJVk2sqpwB43xICIqu4Bj+S3EUVC8Sr32Z3gupdTANGftSj1QeAUP5RGXk+fj2z8gQiAh4WYnKiIz9bf8nSmP5hrUUFJSiVcISyj7m8w+a9ct0PQOlabts1aEZlhkwVThq51Ll2FFj2U95QsxwdtpQUtUdNuXvVMbGFCA2QBprlvQoI2Qlfh7jLKTQf2VgUofvyxCiPuR2T1T2Rhc/IIYf6aHC0f6oPg9V/bSVpeUqf4iqTx+PzZNgmh1gHdCO56rja10IGJMUDm0mLI6RpDuwSmpS8vyTjqnawctiVEBgQmyCEZNNVNk5JJ7vH8iOGR8PNn5QvNBGSzozBDTuTEbslKL0N0XqMxMqsj94t3kLQAkXJhg+iB+9WaGyKHx0TMMSbygkQFUBlkKYz2ftlFY/ET7OyVzpgpiZ05cQDACPszV8kR/O1Oh2GeT78FitCRYtn+HgE8xmHt1jHsND/+hnPMSqeEbixS8Yzegs84xv1Gfw4Zc9i31XusmLAPI4IGP7fsJe0qFy4Jz5jkLBfR9SSEdg0YRXus72OzMBg5bHu8/WMdq6qHCUHgRRbMsi34edQvrAzF7sg/I4OQ1x1ixVunep/q/yUMLh858s75xvfL9fKmz1nSrP5EVCAA662Pz6VVerLFXkJcgoVlJGbQopPAThPGsSoCheqKZ4mQkmXb7652Zb464VI/6oLpzJUHGZoXwOvd4qzKwmV6L550Z8o83xZnWrSeFwXvBrtDZa2Ji83KEDVmA2xlfMigO6vQvu/2ZO7fvkTcWCtmZHqEYNbV5j7pGHsp9K4YeXJOuDXh6c0wgN/qCistpvp/i2qIyAgY+++yxBTF4ElscKzKaSYSeKppDulqYAnFlIumeZ5kjxj4pk+T5rKcFzGpL2QXGsJpYg+IpGH3RdSgfnbIfwwafRFyOn/FUmkwfdb1FewOzZlGNzOWVlK62LF3dfuBA6lFAkCVzntnVkPn6EPsT2qa0s0fi1wTHm9PFbP8ET1VcsQfrdNmPXfsECqyeTgUqTTF7uwQEzsOtjWd6XqsQ2rZAkxG62Ez74cvD6Qwg06ZP6jKh4Dkr64F4o3TNT83L2nfHnaNZR61Sf8WcwgeAjwy9B/7OdMvFEha28fWLzZoj/yOTlDkHcGOEd2sd01P2mFZIQC/4A9yNiICjBZZufn0NovNBKY3bkEwgcxKGsbk9KoR0f1Ca3Xa8QZXqjkI9be/mOTtqDAx6u/jEmjFfsp+I7XTUFtP9EJiEUWRcFrE1e7pACM96XyWjrCslp91R4FW8cilMazmHOXXA0AEXVUFZtIyJ+Zl73nD2eh9wrv7Ss4jl9dsSoItRN+4JDOixfXio/vqgwnxupzyieqHx0j4LCWhpCjEWfMiYvOaq8ByjzyehysprWWRHCW/IdRPNdAg7g+CbCy45QPtHVi1S+pHB2SseRiukf6u1VX4HJ9UH1xHW84ihV+nWtkFO18YUtlD+3JBz82W7jhB00ecnB7KvI81KVu80qJL2kbrYTNg3hVV5qAaGM3MvOy604CooXnfKPdVUPGip/Nqr1QQ/FSqKNWEl9rS2xrp5Q9EP7y9KjOtPXsFbeVuUFzMBuMrqYrKb1VRT8ASULyfVGgsb6mE8JPvru8YIjfPySAZbNO128+Kei1j6offSLtuJHqr5sp2r81M9VPxkvrM6hFiE4UbZGrViHQv75W5IXLKRUGQcoRodbVkjMpLs24VU6yod//R2Xb4ojlDiOg4IiZIfmVFTtgxoV2o6VtAUCjzEvCAoMa7dGVaw3Ki98ZjNU+dcw3l1Wc88sU0KDpcAumZtJ6Z7SwnY/txA8KBlCnMupBUV6up3MgIx5AVRsHCV9oH1uVvd9eKcw+qMUD45tHsTPwrDc5CzOUdvLw49IKJWghIc89o2Y9+VgJ6xMXurywqRd9JEOAT/l3ZpC+qjOD/UvpHix25bOCJxSfZEyzB3X83v1QXwozlR5zt20k9evPo8aW8ZU105XtgGg+15YVzp+tvIibMJkmE7YFSj6ZW15ydYPoZyL3StmzZ1CnN38JMnYbXWEnctoGQhgUU4626u+c5YXJrvJ/aebspAmkHtm4Pgb4lJGwdE5uGqVnXRmGdwAK5AhoPFgo8L9fzCFUy1rz1JehI9zeSicWoFvy4eUSeCIGw0Ay7XJB2FEU6zIPAigQPAljZ5o0madsbzgp5vlPWhbst16YNnWMzIENoGAFAhLaBz/8cv9TXR8QCc9Nodujji3hx7bcFm0LdnY7oaYKeIJYJdifwyBM0ZALxzLtbhbecZdHdU14TMFNuxWQ+/vgvu/RQtJDw1KyJ29uGe3mCFgCBgC4xGQbrlWaSyk4j+RKCok/zinCVWBma/j8beShoAhkCMQduaKR1JaFZIUEQVw4rKrxKE6I0PAEDAERiMgPYIyQpfwTwLS83yxzlaF5Dn4qp01Nf9hAFPLyBAwBAyBwQh4ZcRqi+MVReuISov/l23/aaqArVD8SbWfAd0vZveGgCFgCOykP9gkC58uBf90EZkuC8kVUoWczzj6d1vFFluiIWAIrA0Bzi3xT2eryohO/R8zqE/02su74AAAAABJRU5ErkJggg==\n",
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAANoAAAAsCAYAAAAD1qbfAAAKeElEQVR4nO2dfZQXVRnHPwusu2Am4kutGS9aYC+7te5GGGCvdDRf0jLqZNIax1JPURaE4MtZPCWiQlJ0tJPQlqFoGG9CYUlaUGiKJZqlbvRiIEilokWIS39877TD7Lzc+c3L7zc0n3Pumd2ZO3Pv3pnn3uc+93nu1nV2dlJSUpIt/WLmH5RJLUpKKqcR2JdBimIA0N+2kraCNg74BXCr7YNLSnLifUA7UJdyimIY8CdgAXBYVGYbQTsDuBOYBXzQIn9JSZ6cAGyqQrndQDMa2e4HXh2WOUrQjgRuBr4I/DSN2pWUpEgd9qpeFjwLXACsB74XljFK0L4A7AS+k0atSkpSph14IOY9rwXuAX4H/Bb4UAr1uAKYAIwNyhAlaGejkaxaPUaRGAq8ALRWuyIFJm4bvhtYF7OMvcDngTci4ZhPciPfX4HtwEeCMoQJ2mBgJBoeS6LZCrwV9ZQllRG3DRuB/8QsYxsayQB2AP8Ejoj5DD92AGOCLoYJ2khzHJBCJYrAOvzNvHda3r8XeJL4L/5AIs82dKx+SWgH6tGIlJSDgNcHXQwTosHm2JBCJYpAK3ApsMhz/l9VqEtRyaoNXwlMBG5ynfsA8IMEzzwcGTAmk87UqBE4NOhi2Ig20PWAuBxnUiV0Yd8DpsVxqGO5F3jak56PuLcO+DLwB2A3UiGWZlXRCkjyLuKWM5hs2vBs4Cvs3+kfiQx1ldAALANmA7+s8Bl+z+xPwMAUJmjOol0lI9rdJhWFNuBl4KEK7p0GnAdcBIxC644/Sa9qicnrXWTZhouAu4D3mt8PAZ4LeNZygr09JqHvugupuTdXUNcgHO3Qd7HbZsHaZpW86LSh3mgHsno56TZXntOAG3zuPRlYgz7mPwMbgW9lWdkaJes2XEGvw8QpwNqAekwGmoDhQA+yBDaZdCsywX8UOBP4jUnN0X9eJKFykpWhY3jA+XuA36PJ7iRzbiEwHfWG1aINuAO4xHP+WdfPLfj31suA64G3mJ9vo3KVJguG51RO1m34Y+Ba9EGPAm4PqMffzbEVDSQ/R8LvsJ74Pr6JsSkwSr+Oyzmm3BOBzwDnAxenXEZcWoENyOLlTu6X3YIssQ8Aj9HbC37DnF+D/pZu4A251Lq2yLoNX0Rm/xORdTKKFjQ/3BGVMSV2hV0MEzTn2t/SqwugdYwpaGS7HfVS1RS0EcAQoucWLWidpx2YgybvDt3APHqdW1vSr2ZNk1cbrjD3bbCoUwu962V54MiJr0yFCZrjJPlIqtWR/u02p/4KOBqZcKtBmzk+jf5md3LCIBqQFfZ68/vDaJFzOtCBvAxGIsfrPUhF/n8irzZcidRGG0E7Flkx82KzOR7ldzFsjnYSGvbXpF2jGsP5SB7znN+Hwh+eA95krveYa61oEt2A5iTDkGqzEVnGtmda49ojrzbcjlRLm/l8P/PMY9Ao2hOePTHfRQ7G4/FZSA8a0YYiq8ws4CXLgsYDP0QqgLMWcj9wlSff29nfQjMGNUTac0FbZuAfj9SPXhNyC+oh61EvfAHykbsSOB711Ecgq1qe6koQGwgParw35fLybMMVlnWaj4R5iykzazaiEXcKPnLlJ2iHAIvRWsMCy0JmIutOGzLRzgNWod7qZE/eo5H6MAotRE4DvmZZjkMX+mA6Yt5XKc3IovZrpNJcjtSkWmUZ6iS96S/m+s8C7usiu3bNuw3XoRGtnvzc4s5FWuI13gte1XECEq4r2H/9I4xXoV5pPRry93iuex02FyO9/T70UhcSX9Dy5kvmeGlVa2HPdT7nrkWaShd6X3lTtDashOeRxjYHWVZPR8a/PoL2KJqTnYu2Lthq8fDjkeA8Tl8hg77rIXuBz5rkR4dFmTOAqzF/REkodcA3gQvN8XME+/aV7ZqcMWjrjwX0run1EbStyNR+OdLjR6MwgjAeRXr4p5D/2WLkLhN1XxK2UX4MNvRH7kuTkDozPSJ/2a7JGAf8CHgP0tj+R5Ax5Crk9n+jxcN3mgKWItVxCfCMKfCEyupbkgL16F1MAjqJFrKSZAxCMvB9PEIGweb9l9Fi8sXIRyyql3sE+ZQdhJYFPm1+fxvwGnono++KVfVkFC0qPE2f0kb00k8FpgJzU3x20do1TcLe0URkr/DdKS5swfphpHqcEaMie9DWBxORceRwU3g1SHv7saxTWhwMrEbxWheRrpBB9dupVt/RWebo6+ARJmiO79awkDyt+Mc6vQ54MzInPxVRwUpoQkaYwEA7F0mjfpOyDs1b3ZyPgiG9DriQzMR+KJofv9Pc7+cpH0acdo3iMLTAnEcsXJosRbu+xWWoOb7gd9HGqbgp5NoU4Am0WLcQBdItQaNhAzKQZLEiPxt5GZwVlZHeqN8mT/p4gvK70LzHhlbgQfNzA4oSno1CPq72ye+8ExvHWS+3AO8w5R1r6uhNYfGFcdo1ipnI477bc77WO74rgcuI39mE7utoEyYTtqq+wjxjNJqTNSLL5S1oLeEJuzpmhjfqt1rlP4i2ObsDtVc7wftdNCNtYnXMsvqh+THofYz2ybMD+w4iCYPQx3u6zzWn40tzu4Mu1J6dFnlb6RXqBrTkcSbq+Jxg0y3AJ8w1W0K9T5LGoy03KW86sFOtkkT9pkEb6q2HIJVkDXI9+ndA/sHIVWku8ZdHepBXTxI6sGvX5QTvWv1J9Pf1oHm6m6J0fCuRxhNH0EKxUR2rGZCZFJuo36zL70GbyMyh9yMMYjzyLZ2XfdUSERXFPA59zF4LZS11fJuQmjyWvtrFfchiPhB7QqdIYSOaM0fYHaOwWsMm6vc0ZAa/MOQ5M01yaEAvbKrr3CnIm8Zb/t3IwGCzpriKyjZDypuoKObh+HsVuTs+N6vR9gJZ4+74LsHfVQ1U93rkl+udYwbhON/7zq3DBO1FcyyyoLWi3ZOeDMkTFF7v5kb2D52fgwL9vu465xcg24oW/6chj/pO8pkj5UVQFHMj/t9NUTo+R+uIM6LtRuvFsQXNWaT2818sAiOwj/p9CjmBHozWADd78vzDJIdd5vcwAXbK34SssOegMKLHkbHoQCAoinkn/v/KqCgd3xBzfCaiHm52E+LYETZH60ZClvtGJilhE/UL4eH1aZTvfDQrUS+8CJngDwSCopgfQhHTbuJ0fH77irhxOjkn7fI5550L+3V8l+G/zNOMvok4Abz7CNnKPEyIXkJxQ8NjFFZLuKN+t7nSVuAV5lpQeH1a5f+R/dWia5CxYDl68UXHHcXs/pbWos113G1ZpI7vJLQGaMsA1A53BWWIGq1uQo7Cg2MUWivYRP0GhddH0UH0XGsG/l4R56F9JbZYlFPrzMc/inkzstx9zHWuKB3fQLSu9u0Yz/4w+q6WBGWIErSlaAi/gQPzn10EhdeX2BEWxTwLeQ45o1VROr7JyNNpo0W5mHuvA75KiKoZJWj70BrJUWiRcYJl4UWhaFsUFIm1KPjxmBj31ELHtwcFx0ZRD7wffTerkPtaIHWdnZ02hdchL/6xpKczl5R4mYusd6ciDWoq8eZKeTIC1W8JfZcR+mAraCUlJQkoqum+pKRQlIJWUpIDpaCVlOTAfwGw1+kazeuL8wAAAABJRU5ErkJggg==\n",
       "text/latex": [
        "$\\displaystyle \\left\\{ \\dot{s}^\\mathrm{pl} : \\frac{E_{b} \\dot{s}}{E_{b} + K}, \\  \\dot{z} : \\frac{E_{b} \\dot{s} \\sqrt{\\tau^{2}}}{\\tau \\left(E_{b} + K\\right)}\\right\\}$"
       ],
@@ -1254,7 +1254,7 @@
        "⎩                        E_b + K                 \\tau⋅(E_b + K)    ⎭"
       ]
      },
-     "execution_count": 18,
+     "execution_count": 29,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -1278,12 +1278,12 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 19,
+   "execution_count": 30,
    "metadata": {},
    "outputs": [
     {
      "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAIAAAAAyCAYAAACUPNO1AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAJa0lEQVR4Ae2c7XHUPBDHHSYF8FBC6ICXCggdBKgA6ACGT/CNgQ4CFTDQAVABgQ6ACoB0kOf/02k10lnWyT77cufcziiSpdVK2l3trmRfDi4uLpoSvHr16kjtp+CofL+Eu2/bDg5ITrc0k69Kr1V+W5rVtVKjOj9T+0+lX0oPSrj7tu3hgOT2w8vrhco/ldjEWTjosgDq9EY9UICnKr/L9t6RSs2fdWC97vgpnyk/92Wy60ow6YbSmfCrLZ2nfVf9TpSAL0rnqk82jJ4/qh4cxmX856pDUC1QPZvui/KnrcYeFerPurAErO22ntnICWQVQIjHwvqs9G7dSSSjXfKD1oK/+6H8dm4qqneKorxaAYyOp/1L+U2ri3PVI/xHSo9VjpUvRmvUhtD+KXXOM+mw4sHT+y20rGJ3uYDn0FXntTRwxdw22qy1oNQAO7QLsHStXdKFbPWijc8FPi2y9K/an1Cj/IFSp/A9Du3/CS+rpOD0AT8e6zpW2eYZSByGUlqAWSVGpdi78WS7GstWAsxvXzDlSmiL4ezmF0qnKlcrlhda3zmU8D+oEevGPBO307IAftIQK2oqCDsGTkhaX6LYeoYxDjzjk3ZrW5Fj2puYtsrsNoIwfH218FeMM7TZxm+5p8MCxb+Ftl1sQiCJcCUYlCJRdNUlO6RyodAO/UQDk0/Q1cuFCp/A+48fE2GhPMn8fNtoWUkBRhvksgmJiW73ax5HKhONA5hn6lu7gsZaiGg75dIzFgVB2q6rIqV+uI9wMvB0odNLiaoGi5BaLiBqm1PR/D9BmEtaHIwlam8JSnW3lL4rJT69gyFG+4/w3yhx8ULQhbLZ0bCj66JaeEcqEaQFK6JnAvHvC4zp/l4VBXAWIGawygg+idpVh1VoPB7tsUBoyoFZF87+7vQkJHYuQABYDer/WemJEspzX2ny+5erogAt/++l4q64KYvZCN8d13wbgv3gy6XM0Y6FpTLKg0vAktBeBI+P5WBM5sTtnSlRse+6jbNXADHSdmjLnHvGGw8RvttxqkdoLkZQ+ZnSRyVnHQyZXHWdtNVsAqyyAqKF/z9QP1wKysO4rTFVPyrMXgHELfPRyQkg5qIYfaTnu8rPfT2CZRdzE8rO/KZkAlUxQCdt9WM8aJyo3ClItRE3hLsH+ikZ3TDQVIWroAAuEBNTs/5c9QibYCu4A5URAIIxhSC/o7QMRdpCNpolK8D4hufoa1zoonw2vquf4s/hFEQvm6YYx457r8TOJjWqs+MfjwA4CJWce/fYQlBvAZ2KDdey7GboGG3ymPY3tYVXrypjMSymwJzjVnix5uiobPBYhWPVgwtNgIBy0uPfYpimGVUBNOlR3mLZ5Ibmmse5+iZv43rSui4ascV4qP4IqqmlLTwUKFYiurfAjxOP1cKZsuLaWMS1ENsR7J5dB+4H3G5Uzs7kTVpyZNz1Bdr8x1QAdt1ob7FsgpeUYz2IAbjVwxpsLCibeL03lunnXIDza0JEoL1AjOrdp9cAG0LWOjDJG/HBm1gSclFiKGfV4jFzFuDYI7TOzXHHfXnnOIBSm2zD5HMWgCMLPjCOikOH5YLwLNptaZfHXflJmWhsRfC4vLaZPb/WerjQeqIUrpgTC6AGO49W+bwIH5/JUQmlIQ8pHkz1LVA7ioPbmUPw2FrftlSIzwSxCJ6PU8zNNw0fhZJevnx5qvRP6ZbVlXLwlE5iHD1/jp9ry+p3vRZ3j7eQ11A+iNeJnA9gvrSCmzB8RPGDRbV3grSKiw6+gFnn/N1JP27QGOUfM8TIV6QsnvAeoQqEeyJEXPdpHAP8rerdjfRITRsJHPsstnu6V77FHQmviZkcEeyrmN8qs5OHAFp1NqTjvs/mOCD5ciXONfk9ld8mvwtQBUEg157ZHxGoPgvqR1DBO+xqM5QltK+clAOSD2afyy3ki8tvlk8BXH7gCpZfnIBbAnb/8kuOEn5o00S4abtQmvzzpzDoFSyIv9wBIHyO5U74sCFRACoEWAG+ZGldGrjW/J+7vl++tVCrcc7VjPLs3UeBTyM0uRdT4ne4A4BmHATaGHYBxF2Ala0tm4voWpG/+lsMkqW/qlL90Wzma3cJKBOKZWB3DQQ+2Z9IGeJl5X4NbCSsKQDvic8S3uoZ6wwO62Od4UtilUvAhg473xBzCmCmHKbtBIgpvId/q9x++5e9yFK7KcrWrcuvofFr4CY2uwZN/IOffPWRXbRMlibbsP6cAoTGXSpokWg4ULJamL+1rI0bIfrjx+Ur3sS0RijVRdGwE9inXCe182r6r/LEKuRwO+qI7xKYjQJoVbZjVt1FhO/vEk4Mf2B32Q4bTmXR05Q4WYMEDn3e0fT6jWHNZHJBYE2/bcRxzBOzEgugZ8y+A5XPVUjafdO2ZFymNfEaVMYqTPYbwzlZABiVCFfMQykQegDVtQKh0Hj5BdYQ5qe5YvJ7/8awzzJmoQBe0KwbX2x3GJhNFGBUn88gU0C0BqfEesZycXHTCtzGHH8uLsD8/xS//RuT3yVatobBvzEsEe9qm4UF0OLY6Y12TWw+OUol0bSeuXXkbM1n4Oys6t0lfC7I3DiMFcENymrPfULGOLURu9FmfvZ5ORYAN0AAmKxFz6PAXBSg5f89dxCaAzEVlwAzjbkw/J5rrPjTIeBG9VzK4HqMbgW1LIpbg+iE46TK9mXWscrczgYFz1IYULnzCiCm2M5Jjk7wQm3xDkf4jrmqh9kuRvD9uYGrvlgR7qjg5wDN1hpUhxVgjViBWmsi1DqYQwxgvjM5AcTLF4OP9Dzkt38xmSnLnWvQ3FkXilz8jeHQyc1BAdzduRiVNY9+dw397d9QvvbtV1yDiJkrwwqMCiUX4IKbUUcbiZiEivl+r8TOJjWqs+MfjwA4d3xe/ds/Om4CojUwz3gNQ39jWJp2pyxbCqCJEYVCjIltJTBHTWwdf8hpILYYD0Xv8SYXW7sG4fEa173KXWN+Jkv4lkCXC8DvEHlax6TTDB6Irt3alBMcrvPbP5jaYuyW8YggErA3iYsn/W1ZAN+CxuE38Tnrap8nuVUZ1iP8YwYpgQVhvSepvp3BZ29iE3Twio4ccYOx1XOjJd8ExuMLmeMH15Hcrk1yCRGPty+PzwEv/K+iTIyR/c6zywU06my+h58TWRQ6/iz3FCfhgGSG2ceK4+qywmfgTgtAIyBCaI9TAJUHm0pHbP9nIxyQnLjoYue/Vrl4Q/k/XtnNam14LSUAAAAASUVORK5CYII=\n",
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAGAAAAAmCAYAAAA7mZ5JAAAFxElEQVR4nO3bV4we1RXA8d9ijEMv8gMGATYlAYQNZi0sQPDCC8VAQCjIQQEsiBR4QEhUQ0ALInQhihBGIMoDCAgtCQGJElEEpGCHllhGVMlgShQI1dgJ5OHMaGdnp33fN9/uAvuXrubbmTn33rnnzjn3njM7MDQ0ZJJKtsU/sR/+3nbl6zS8b1Ns0Hbj3xHewx5CCZ0yve6GOgVsjNuTskUXHfg+8F+8jq+7kD0Ff8ZeZTdUKWAd/A5v4qdY2UUH+s2f8G1Beaih7B25c7/Elzi7pf6djxNwJ/YsumHdCuGF2Aq/aakz/WAuzsUtufNfNpRNFTUN14uJdjgea6l/8A+cgyUK3oSqN2AR7hGvYCfskJR+swM2w1N4P1c+bSi7FNvgGTFD5xke/AGciRVYjQ9xb5d9/UNS9+z8hTIFDGAfvN1FY08kpd8M4n+6W5kMClO1BZZhOfY18nnPEJPwZPwEh+n+zfgKn2D//IUyBWyN9VWbqPFmEFPEzPw8U+5Ori/ADRWy3+C3uAzHiUHKciAeFpPpHeFMb+yhv1MVWIayAd4kOf6oi4ZmdiHTDYO4z2iH+UlynKP87RgUA7uzEueIB3A1dk9+341/dd1b1hPL+RGUvQHrJcdpPTTYb+biWbFEzJZ0kObgx3hBmJjZOdlHcWhShgrqvy6Rf1isjt7ALj30d6qC8azbB0zpocF+MkvY7yr7P0dsouYJM3NmTnYZXsYx+DV+XlDHG7gqqWMgqbNbBopOTmQbX8VgcnwfW+aufSSea31hQoiBPjonmyrv92KZeItwws/hLHyAv4pV4HFYgydbe4KEOgV8XnFtunCAhZrNsEbYvtUd9KuOdBCX585/i82Fs1suHC1hcl7MyL5p2FfA5cK8PIj5wlScje3whXDABwildEvh3qRMAen5dysq3AgXZv6eKWbKUiN3oh8pH/zbEplFye+mLE5KGXOwvbC7m+JXwtZXyS7K/L7QyGdrg3cVjHeZAmaJV+/5igrfNtJ5nSgG8y5c2U0PW2S2WCH9TTzj6cJcjSdPK9iIlSngSNyvs07vkRxf7EBmMS7Fqg5kmnBacjy35Xp7YYlYtc2U2fAVKeBA4fXnd9jA7snxpQ5kVml/8CcqL+A8oYhDsZbRCrhI2M698HEHladLtPeEzZ+kmMvFXuUJsbdYkd8H/FGsbvbusOJZYvfcyez/IbI5jhLBuZWM3og9j1/gJhGcasrOyfHVHjtYRFG8f6KWKqaI2NNLuEIsbwt9wAe4VdiqUV67hDR2VBcGzjNDLBNX4T8l99TtM74rLBR7kAXZk2WhiKewm+GVTR2vJcdTxdZ9YUO5S8SG6YgG9/aS/WqDXjNox+IvcnuiMgWkAa1dG3ZumYinrBVKmNdQrhPS7NeMXCmK4TTlNsWBuLL2lya/p+FmMYEOF0vpOnZVEE0t2wek9mxGw84RqctO05fHJ6WOfPZrrMln0O4TYzdP86RVPmaF+mjo1IaV95tesl9ttV+XQaujMLJcp4A6zz5W1GW/xqL9ugxaV9SZoLVtNNICddmvBTgEJ9XUc05SUqaJZz09c+4gkaTPt1+XQavjGwUTukwBqXa7+RipH8wVu/TXS65XpR+zLBFfeqRcJqKU12bOFUWA5+Jikah/VjjuoQbtZVkjWftnKVNAGk6YCApomv1aKeItG+JneKXgvn8nJeWz5O8yxWbbz2bQ7hdL7zsbPUHwtYIwTZkP+FisNiZCSjKf/cqWtH9l6cc22y/KoO3TQT3rio+0Rp0s4xHs1EED/aIu+7VaefqxrfbrMmhv1dSxo/A3o74rGqj4PH23RGCXXOMTjT2Ffzg4+fsE8cBVGbOx5iYRUjkxf6FqGfoqLsDjevsaoN9k04/TRfrxmnHt0TADYmU233CSaAR1SfklInq3WHj9FS12ri0mYvox5RLhfPcVDn8UVSZokjGg6X/ITNInJhUwzkwqYJz5P4/udg2SMK7rAAAAAElFTkSuQmCC\n",
       "text/latex": [
        "$\\displaystyle \\left\\{ \\dot{\\tau} : \\frac{E_{b} K \\dot{s}}{E_{b} + K}\\right\\}$"
       ],
@@ -1293,7 +1293,7 @@
        "⎩               E_b + K   ⎭"
       ]
      },
-     "execution_count": 19,
+     "execution_count": 30,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -1627,7 +1627,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 20,
+   "execution_count": 33,
    "metadata": {},
    "outputs": [
     {
@@ -1637,7 +1637,7 @@
        "       1.76641039, 2.05694336, 2.31490712, 2.53621714, 2.71736924])"
       ]
      },
-     "execution_count": 20,
+     "execution_count": 33,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -1653,13 +1653,13 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 21,
+   "execution_count": 34,
    "metadata": {},
    "outputs": [
     {
      "data": {
       "application/vnd.jupyter.widget-view+json": {
-       "model_id": "5413ad5ab7174be9be8f2d585b3b1675",
+       "model_id": "4ce9295f06e742c890468b129c0e9ffb",
        "version_major": 2,
        "version_minor": 0
       },
@@ -1673,10 +1673,10 @@
     {
      "data": {
       "text/plain": [
-       "[<matplotlib.lines.Line2D at 0x7ff9f7a9aac0>]"
+       "[<matplotlib.lines.Line2D at 0x7f4c3d69fee0>]"
       ]
      },
-     "execution_count": 21,
+     "execution_count": 34,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -1710,7 +1710,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 22,
+   "execution_count": 35,
    "metadata": {
     "scrolled": true,
     "slideshow": {
@@ -1746,13 +1746,13 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 23,
+   "execution_count": 36,
    "metadata": {},
    "outputs": [
     {
      "data": {
       "application/vnd.jupyter.widget-view+json": {
-       "model_id": "cbe6af1675a5402d9220714879882fe4",
+       "model_id": "2ebcbcb77ee245a39b472cc87eb08e13",
        "version_major": 2,
        "version_minor": 0
       },
@@ -1885,7 +1885,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.8.8"
+   "version": "3.8.10"
   },
   "toc": {
    "base_numbering": 1,
-- 
GitLab