From 3f0a418a3ce997a9614be6bd9d9a2917f6b5e607 Mon Sep 17 00:00:00 2001
From: jsh9 <25124332+jsh9@users.noreply.github.com>
Date: Sun, 8 Jan 2023 18:16:41 -0800
Subject: [PATCH] Update example and docstring; v0.5.1

---
 PySeismoSoil/__init__.py      |   2 +-
 PySeismoSoil/class_curves.py  |   2 +-
 docs/source/conf.py           |   2 +-
 docs/source/index.rst         |   2 +-
 examples/Demo_04_Curves.ipynb | 203 ++++++++++++++++++++--------------
 setup.py                      |   2 +-
 6 files changed, 128 insertions(+), 85 deletions(-)

diff --git a/PySeismoSoil/__init__.py b/PySeismoSoil/__init__.py
index 8de3123..034e128 100644
--- a/PySeismoSoil/__init__.py
+++ b/PySeismoSoil/__init__.py
@@ -1,3 +1,3 @@
 # Author: Jian Shi
 
-__version__ = 'v0.5.0'
+__version__ = 'v0.5.1'
diff --git a/PySeismoSoil/class_curves.py b/PySeismoSoil/class_curves.py
index 384f689..486d7f2 100644
--- a/PySeismoSoil/class_curves.py
+++ b/PySeismoSoil/class_curves.py
@@ -211,7 +211,7 @@ class Damping_Curve(Curve):
     ----------
     data : numpy.ndarray
         A 2D numpy array with 2 columns. Its 0th column contains the strain
-        array, and the 1st column contains the G/Gmax values.
+        array, and the 1st column contains the damping values.
     strain_unit : {'1', '%'}
         The unit of the strain.
     damping_unit : {'1', '%'}
diff --git a/docs/source/conf.py b/docs/source/conf.py
index 842402f..3987f82 100644
--- a/docs/source/conf.py
+++ b/docs/source/conf.py
@@ -23,7 +23,7 @@
 author = 'Jian Shi'
 
 # The full version, including alpha/beta/rc tags
-release = 'v0.5.0'
+release = 'v0.5.1'
 
 
 # -- General configuration ---------------------------------------------------
diff --git a/docs/source/index.rst b/docs/source/index.rst
index a94e964..d4c9802 100644
--- a/docs/source/index.rst
+++ b/docs/source/index.rst
@@ -73,7 +73,7 @@ To report bugs and submit suggestions, please use the `"Issues" section <https:/
 Current version
 ---------------
 
-* v0.5.0
+* v0.5.1
 
 
 Indices and tables
diff --git a/examples/Demo_04_Curves.ipynb b/examples/Demo_04_Curves.ipynb
index f7c1176..f9eafb5 100644
--- a/examples/Demo_04_Curves.ipynb
+++ b/examples/Demo_04_Curves.ipynb
@@ -48,17 +48,6 @@
     "from PySeismoSoil.class_curves import GGmax_Curve, Damping_Curve, Stress_Curve"
    ]
   },
-  {
-   "cell_type": "code",
-   "execution_count": 3,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "strain = np.geomspace(1e-2, 10, num=20)\n",
-    "y = strain ** 2\n",
-    "data = np.column_stack((strain, y))"
-   ]
-  },
   {
    "cell_type": "markdown",
    "metadata": {},
@@ -68,26 +57,28 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 4,
+   "execution_count": 3,
    "metadata": {},
    "outputs": [],
    "source": [
+    "strain = np.geomspace(1e-2, 10, num=20)\n",
+    "y = strain ** 2\n",
     "y_ = np.flipud(y) / np.max(y)\n",
-    "data_ = np.column_stack((strain, y_))"
+    "fictitious_GGmax_data = np.column_stack((strain, y_))"
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": 5,
+   "execution_count": 4,
    "metadata": {},
    "outputs": [],
    "source": [
-    "ggmax = GGmax_Curve(data_)"
+    "ggmax = GGmax_Curve(fictitious_GGmax_data)"
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": 6,
+   "execution_count": 5,
    "metadata": {},
    "outputs": [
     {
@@ -109,7 +100,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 7,
+   "execution_count": 6,
    "metadata": {},
    "outputs": [
     {
@@ -121,7 +112,7 @@
        "        2.33572147,  3.35981829,  4.83293024,  6.95192796, 10.        ])"
       ]
      },
-     "execution_count": 7,
+     "execution_count": 6,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -132,7 +123,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 8,
+   "execution_count": 7,
    "metadata": {},
    "outputs": [
     {
@@ -145,7 +136,7 @@
        "       8.85866790e-06, 4.28133240e-06, 2.06913808e-06, 1.00000000e-06])"
       ]
      },
-     "execution_count": 8,
+     "execution_count": 7,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -161,13 +152,46 @@
     "### 2.2. A `Damping_Curve` class object."
    ]
   },
+  {
+   "cell_type": "code",
+   "execution_count": 8,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "strain = np.geomspace(1e-2, 10, num=20)\n",
+    "y = strain ** 2\n",
+    "fictitious_damping_data = np.column_stack((strain, y))"
+   ]
+  },
   {
    "cell_type": "code",
    "execution_count": 9,
    "metadata": {},
    "outputs": [],
    "source": [
-    "damping = Damping_Curve(data)"
+    "damping_curve = Damping_Curve(fictitious_damping_data)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 10,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAASIAAAEmCAYAAADGL52gAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA3JElEQVR4nO2dfXgcV32o39/uypJsS/JH4o8ojpPYJE6wIU4iQgxOwIEYAoXQ75u2EL7aQChNaSHNbcotH+2l9LmQtgR6S4EC5bYUAjEQJVaDSORYjuUkchQRmciOrMiKbEeW9WFpV7sz87t/zK68Wu9Yu/KOdjV73ufZZ8/OzM6e887op3POnDkjqorBYDAUk1CxM2AwGAwmEBkMhqJjApHBYCg6JhAZDIaiYwKRwWAoOiYQGQyGomMCkcFgKDomEBkMhqITKXYGZoOICHABMFbsvBgMBk9qgJc1h1HT8zIQ4QahI8XOhMFgmJELgf6ZNpqvgWgMoK+vj9ra2rNu6DgOR44c4cILLyQUMi3RFMaLN8ZNdvLxMjo6ypo1ayDHVst8DUQA1NbWzhiILMtidHSUxYsXE4nM6+IWFOPFG+MmO356kfl406uI1AIjIyMjMwYig8Ew94yOjlJXVwdQp6qjM20f+HqnbdscPHgQ27aLnZWSwnjxxrjJjp9eAh+IVJWTJ08yH2t+fmK8eGPcZMdPL3k3zUTkBuCTwDXAauA9qvpg2noBPgN8GFgC7AY+oqrdadssA/4J+DXAAR4A/kRVT+WYB9M0MxhKmLlomi0CngXu9Fj/KeDjwB3AdcA4sFNEqtK2+R7wauCtwDuBG4B/mUVeZsS2bQ4cOGCq2RkYL94YN9nx00veXd+q+jDwMIBb+TlNsjZ0F/B5Vd2RXPZe4BhwK/CfInIF8DagQVWfSm7zx0CjiPy5qr4828J4EY1GC73LQGC8eGPcnMkLx8Y4evIU631omhW6j+gSYBXwaGqBqo4Ae4Hrk4uuB4ZTQSjJo7hNtOuy7VREKkWkNvXCHbE5FZlt286atiwLgM2bN6OqOI4ztTxbOpFITEunmq2ptKqekU6WcSrtOM60dCoPXmnbtqelcylTenq2ZXIch6uuuopQKBSYMhXqOKkqmzdvRkQCU6ZzPU6qyrvub+X9O45x7FQipzLlQ6ED0ark+7GM5cfS1q0CjqevVFULGErbJpN7gJG01xGAzs5OALq6uujq6gKgo6OD7m63O6q9vZ0XX3yRzs5O2tra6OvrA6C1tZWBgQEAWlpaGBwcBKC5uZnh4WEAmpqaGBtzx2I1NjYSi8WwLIvGxkYsyyIWi9HY2AjA2NgYTU1NAAwPD9Pc3AzA4OAgLS0tAAwMDNDa2gq4AzHb2toA6Onpob29HYDu7m46OjpmLFNPTw/AOZfp2WefZXJyMlBlKsRx2r17N52dnfT29gamTOd6nIbG4yRsNzguXxjJqUx5kYqws3kBCtya9nlLctnqjO3+C/h+Mv0/gV9l2ddx3E7tbL9TCdSmveoBHRoaUlVVy7LUsqwz0olEQicnJ/W5557TWCymtm1PLc+Wjsfj09KO40xLO45zRlpVp6Vt256WTiQSZ01bljUtna0cmWVKT8+2TNFoVDs6OjSRSASmTIU6TrFYTJ977jmNx+OBKdO5HqfO/mFde/fP9LX/q3Hqe2cr08jIiCZjQa3mEEsKPWz0aPJ9JTCQtnwlsD9tmxXpXxKRCLAs7fvTUNVJYDJtewDC4fC098x0avTnxo0bp+0vfVRoerqiomJWaRGZSodCoanh77mkvfI+U5lyTXvlvaqqik2bNpHJfC5ToY5TZWXlGefMfC/TuR6nY6MxAOqXLZ62/mxlyodCN816cIPJTakFyT6d64A9yUV7gCUick3a97Yl87K3wPnBtm3a29vNFZAMjBdvjJszOTbq1gOqmSyNq2YishhYn7boEhG5ChhS1ZdE5D7gXhHpxg1MnwNeBh4EUNUuEXkE+LqI3AFUAF8B/lN9uGIGUF1d7cdu5z3GizfGzXSOjrg1ohWLK33Z/2zqUdcCv0j7/KXk+7eB24Ev4o41+hfcAY1PAG9T1Vjad34PN/j8nNMDGj8+i7zMSDgcZsOGDX7sel5jvHhj3JxJqml2+UWrzmiaFYK8m2aq+piqSpbX7cn1qqqfVtVVqlqlqm9R1Rcy9jGkqrepao2q1qnqBzTHUdX5YlkW+/btm7pMaXAxXrwxbs7kaDIQRYcGfPES+HvNRISlS5eeMfiy3DFevDFuziTVR3TR+XW+eAn8ZCvhcJj169fPvGGZYbx4Y9ycSappdvUVl5RG02y+YVkWra2tppqdgfHijXEznUnLZmg8DkDvgQ7TNJsNoVCI+vp6M+VnBsaLN8bNdI4nm2ULIiEuv9if6XMD3zQLhUKsXbu22NkoOYwXb4yb6aQ6qlfVVnHxxRf78huBD/mWZdHS0mKq2RkYL94YN9NJ9Q+trKn0zUtZ1IjWrVtnqtkZGC/eGDfTSQ1mXFlbxbp1K03TbDak2vuG6Rgv3hg300nViFbVVfnmJfAh37IsmpubTTU7A+PFG+NmOkeTndXn1yzwzUtZ1Ig2btxoqtkZGC/eGDfTSdWIVtdVs/ECf7yURSBasWLFzBuWGcaLN8bNdKYC0ZJqVqxY5stvBD7kJxIJdu7cOTWFpsHFePHGuDmNqk51Vi+rDvvmJfA1onA4TENDgy/D0uczxos3xs1pRqIJJi13DuoLlixkqU9eAh+IQqEQy5b5U52czxgv3hg3p0nd7Lp0YQXVlRVUV5qm2axIJBI89NBDppqdgfHijXFzmtSo6pW1Vb56CXwgikQibN26ddZz6QYV48Ub4+Y0x9IGM/rpJfCmRcQ8ljoLxos3xs1p0u8z89NL4GtEiUSCHTt2mGp2BsaLN8bNaabuM6ur8tVL4ANRJBLh5ptvNtXsDIwXb4yb0xxLqxH56SXwgQhm/6yloGO8eGPcuJzurHaf3uGXl8AHovRH9RpOY7x4Y9yc5uiIe/l+ZW2Vr15E3Uc6zyuSD20cGRkZmbHzTFWxLItIJGImQ0/DePHGuHFJ2A6X3fswqvDUvW9h+aIFOXsZHR2lrq4OoE5VR2f6rcDXiADzn80D48Ub4wZeGZtEFSrCwrKFCwD/vAQ+EFmWRVNTkzmxMjBevDFuXFL9QytqqgiFxFcvgW+aGQyG2fHwcwN85HvPcPVFS/jRR9+Q13dN0ywDVWV0dJT5GHD9xHjxxrhxOZo2MyP46yXwgciyLHbt2lX21exMjBdvjBuX1A2vK2vdQOSnl8APlqioqOAd73hHsbNRchgv3hg3LsfSbngFf70EvkbkOA5DQ0M4jlPsrJQUxos3xo1LakK0VclA5KeXwAci27bZt28ftm0XOyslhfHijXHjklkj8tNLwQORiIRF5HMi0iMiURE5JCJ/JWkjoMTlsyIykNzmURF5VaHzAm51cvv27VRUVPix+3mL8eKNceNyLKOz2k8vftSI7gY+AnwMuCL5+VPAH6dt8yng48AdwHXAOLBTRKoKnRnHcTh+/HjZV7MzMV68MW5gLJZgPO7WfFL3mfnpxY9AtAXYoaoPqephVf0h0AS8DtzaEHAX8HlV3aGqHcB7gQuAWwudGcdx6OzsLOuTKhvGizfGzenaUE1VhIUL3GtafnrxIxC1AjeJyGUAIvJa4I3Aw8n1lwCrgEdTX1DVEWAvcH22HYpIpYjUpl5ADTDVVrVtO2vasixEhG3btgFMCbQsK2s6kUhMS6fGS6TSqnpGOpn/qbTjONPSqUudXmnbtqelcylTenq2ZVJV3vzmNxMOhwNTpkIdJ4Bt27YRCoUCU6Z8j9PLJycAt6M6vU/oxhtvJBKJ5FSmfPAjEH0B+E/ggIgkgHbgPlX9XnL9quT7sYzvHUtbl8k9wEja6whAZ2cnAF1dXXR1dQHQ0dFBd3c3AO3t7bz44ov09/fT1tZGX18fAK2trQwMDADQ0tLC4OAgAM3NzQwPDwPQ1NTE2NgYAI2NjcRisWl3H8diMRobGwEYGxujqakJgOHhYZqbmwEYHBykpaUFgIGBAVpbWwHo6+ujra0NgJ6eHtrb2wHo7u6mo6NjxjL19PQAnHOZDh8+TDweD1SZCnWc+vv76e3tDVSZ8jlO7QfcvC+UxFSZnnnmGfbv34/jODmVKS9SEbZQL+B3gb7k+ybgD4ATwPuS67cACqzO+N5/Ad/32GclUJv2qgd0aGhIVVUty1LLss5IJxIJjcVi+vjjj2s0GlXbtqeWZ0vH4/FpacdxpqUdxzkjrarT0rZtT0snEomzpi3LmpbOVo7MMqWnZ1umiYkJfeyxxzQejwemTIU6TtFoVB9//HGdnJwMTJnyPU7/9PNf6dq7f6af+H771PJYLKaPPfaYJhKJGcs0MjKiyb/zWs0hbhT8XjMR6QO+oKr3py27F/h9Vd0gIpcCh4DNqro/bZvHgf2q+ic5/Ia518xg8JFP7+jkO3t6ufPN6/jk9g15f78U7jVbCGQ2Eu203+oBjgI3pVYmA8t1wJ5CZ8ZxHHp7e8u64zEbxos3xs2ZgxnBXy9+BKKfAn8pIu8QkYtF5D3AJ4AfA6hbBbsPuFdE3iUim4DvAC8DDxY6M47j0N/fX9YnVTaMF2+MmzMHM4K/XvxomtUAnwPeA6zADTD/AXxWVePJbQT4DPCHwBLgCeCjqvpCjr9hmmYGg4+8/m9/ztHRGD/52Bt4zYVL8v5+0Ztmqjqmqnep6lpVrVbVdap6byoIJbdRVf20qq5S1SpVfUuuQShfbNvm4MGDZT9cPxPjxZtyd2M7yiunpt95D/56Cfy9ZqrKyZMny35umUyMF2/K3c3gqUlsRwmHhPMWV04t99NL4KcBiUQiNDQ0FDsbJYfx4k25u0l1VJ+/uJJw6PQk+X56CXyNyLZtDhw4ULbVbC+MF2/K3U36013T8dNL4AMRQDQaLXYWShLjxZtydjMViGoqz1jnl5fAN83C4TCbN28udjZKDuPFm3J3kzlXdQo/vQS+RmTbNp2dnWVbzfbCePGm3N2kP901HT+9BD4QGQyG/Dg+duaoar8pi6bZxo0bi52NksN48abc3aSummXWiPz0EvgakW3btLe3l2012wvjxZtyd3O6j2h6Z7WfXgIfiACqq6uLnYWSxHjxplzdTMQtxmLuZGmZNSLwz0tZNM02bMh/GoOgY7x4U85uUs2yRQvC1FRNnyTfTy+BrxFZlsW+ffvK/qmdmRgv3pSzm8ynu6bjp5fAByIRYenSpaQ9zciA8XI2ytlNtuk/UvjppSyaZuvXry92NkoO48WbcnbjNZgR/PUS+BqRZVm0traWZTX7bBgv3pSzG69L9+Cvl8AHolAoRH19PaFQ4IuaF8aLN+XsJjWYMfVQxXT89BL4plkoFGLt2rXFzkbJYbx4U85uss1VncJPL4EP+ZZl0dLSUpbV7LNhvHhTzm6mrppl6SPy00vgA1EoFGLdunVlWc0+G8aLN+XqxnF06qqZV43ILy9l0TSrr68vdjZKDuPFm3J1MzQRx3IUETg/y1xEfnoJfMi3LIvm5uayrGafDePFm3J1k+ofWr6okorwmaHBTy+BD0ShUIiNGzeWXTV7JowXb8rVzTGPm11T+OmlLJpmK1asKHY2Sg7jxZtydXP0LP1D4K+XwIf8RCLBzp07SSQSxc5KSWG8eFOubs52nxn46yXwgSgcDtPQ0EA4HC52VkoK48WbcnVz7CyjqsFfL2XRNFu2bFmxs1FyGC/elKubXJpmfnkJfI0okUjw0EMPlV01eyaMF2/K1Y3X88xS+Okl8IEoEomwdetWIpHAV/7ywnjxplzdzFQj8tNL4E2LCLW1tcXORslhvHhTjm5iCZvhCbemk+2GV/DXS+BrRIlEgh07dpRdNXsmjBdvytHN8eQVs8pIiLrqiqzb+OnFl0AkIvUi8u8ickJEoiLynIhcm7ZeROSzIjKQXP+oiLzKj7xEIhFuvvnmsqtmz4Tx4k05ukmfEM1rBkY/vRQ8EInIUmA3kADeDlwJ/BlwMm2zTwEfB+4ArgPGgZ0i4ssT3crphMoH48WbcnNz9CxTxKbjlxc/akR3A32q+n5VbVPVHlVtUtVD4NaGgLuAz6vqDlXtAN4LXADcWujMWJZFY2Nj2d03NBPGizfl6OZ4DoHITy9+BKJ3AU+JyA9E5LiItIvIh9PWXwKsAh5NLVDVEWAvcH22HYpIpYjUpl5ADTD1oDfbtrOmLctCRLjlllsAcBxnanm2dCKRmJZW1WlpVT0jncz/VNpxnGnp1EHzStu2PS2dS5nS07Mtk6ry9re/nXA4HJgyFeo4Adxyyy2EQqHAlGmm4zQwHAVgVW2lZ5kAtm/fTiQSyalM+eBHILoU+AjQDWwHvgb8o4i8L7l+VfL9WMb3jqWty+QeYCTtdQSgs7MTgK6uLrq6ugDo6Oigu7sbgPb2dg4fPjz1GJS+vj4AWltbGRgYAKClpYXBwUEAmpubGR4eBqCpqYmxsTEAGhsbicVi0/4jxGIxGhsbARgbG6OpqQmA4eFhmpubARgcHKSlpQWAgYEBWltbAejr66OtrQ2Anp4e2tvbAeju7qajo2PGMvX09ADQ1tY26zI9/PDDjI+PB6pMhTpOe/bswbKsQJVppuN0sN/Nx8raqoKVKS9SEbZQLyAOtGYs+0dgTzK9BVBgdcY2/wV832OflUBt2qse0KGhIVVVtSxLLcs6I51IJDQajeqDDz6oExMTatv21PJs6Xg8Pi3tOM60tOM4Z6RVdVratu1p6UQicda0ZVnT0tnKkVmm9PRsyzQ+Pq4PPvigTk5OBqZMhTpOExMT+uCDD2osFgtMmWY6Tr/x1d269u6f6U+f7fcsU+pvKR6Pz1imkZERTf6d12oOcUM0WQUsFCLSC/y3qn4obdlHgHtVtV5ELgUOAZtVdX/aNo8D+1X1T3L4jVpgZGRkpOzGexgMfnDDF3/BS0MT/OCO62m4+Nxv4xgdHaWurg6gTlVHZ9rej6bZbuDyjGWXAb3JdA9wFLgptTIZWK4D9hQ6M6rK6OgohQ648x3jxZtyc6OqM46qTm3nlxc/AtGXgdeLyP8UkfUichvwh8D9AOqW4j7gXhF5l4hsAr4DvAw8WOjMWJbFrl27yuoKSC4YL96Um5vhiQRxy+1gXuExqhr89VLwphmAiLwT+N/Aq3BrQF9S1a+nrRfgM7gBagnwBPBRVX0hx/2bppnBUCC6BkZ5+z/sYtmiBTzzV28tyD5LoWmGqv5MVTepapWqXpEehJLrVVU/raqrktu8JdcglC+O4zA0NDSrS4pBxnjxptzc5DqY0U8vgb/XzLZt9u3bN20shMF4ORvl5ub0YEbvZhn46yXw49grKirYvn17sbNRchgv3pSbm6Mj7g2vZ+uoBn+9BL5G5DgOx48fL5tqdq4YL96Um5t8mmZ+eSmLQNTZ2Vk2J1WuGC/elJubY2l33p8NP70EvmkWiUTYtm1bsbNRchgv3pSbm2M59hH56aUsakT9/f1l898tV4wXb8rNzbE8mmZ+eSmLQHTo0KGyOalyxXjxppzcxC2HwVNxYObOaj+9+DKg0W/MgEaDoTAcOTnBG//uF1SEhRc+/3bP2RnzpSQGNJYSjuPQ29tbFv/d8sF48aac3KSe7rqixnuK2BR+eimLQFRO7f1cMV68KSc3uV4xA3+9lMVVsy1bthQ7GyWH8eJNObk5OjLzXfcp/PQS+BqRbdscPHiwbIbr54rx4k05ucn1ihn46yXwgUhVOXnyZNnMLZMrxos35eQm1zFE4K+XsmiaNTQ0FDsbJYfx4k05uTmaRx+Rn14CXyOybZsDBw6URTU7H4wXb8rJTeqqWa5NM7+8BD4QAUSj0WJnoSQxXrwpBzeqmldnNfjnJfBNs3A4zObNm4udjZLDePGmXNyMxiyiCbd2k0uNyE8vga8R2bZNZ2dnWVSz88F48aZc3BweHAdg2aIFVC8Iz7i9n14CH4gMBkN29vcNA/CaC+uKmxHKpGm2cePGYmej5DBevCkXN6lAdNWaJTlt76eXwNeIbNumvb098NXsfDFevCkXN6lAtPmipTlt76eXwAcigOrq6mJnoSQxXrwJupuT43F6kn1EV124JOfv+eWlLJpmGzZsKHY2Sg7jxZtycLP/yDAAl563iLqFFTl9x08vga8RWZbFvn37yuapnblivHhTDm7aXxoG4KqLluT8HT+9BD4QiQhLly4t2IRPQcF48aYc3Ez1D+XYUQ3+eimLptn69euLnY2Sw3jxJuhuHEd5Ns+OavDXS+BrRJZl0draGuhq9mwwXrwJupueE+OMRBNURkJcvqom5+/56SXwgSgUClFfX08oFPii5oXx4k3Q3exP9g9tqq+jIpx7Gf30EvimWSgUYu3atcXORslhvHgTdDftfScB2JxHRzX46yWYIT8Ny7JoaWkJbDV7thgv3gTdzekR1bn3D4G/XnwPRCLyFyKiInJf2rIqEblfRE6IyCkReUBEVvrx+6FQiHXr1gW2mj1bjBdvguwmGrc5MDAG5HfpHvz14qtpEWkA/gjoyFj1ZeDXgN8CbgQuAH7kRx6C3t6fLcaLN0F20/nyCJajrKip5IIcZmVMx08vvpkWkcXA94APAyfTltcBHwQ+oarNqvo08H5gi4i8vtD5sCyL5ubmwFazZ4vx4k2Q3aQ6qq9asyTv8UB+evEz5N8PPKSqj2YsvwaoAKaWq+oB4CXg+mw7EpFKEalNvYAaYOrmO9u2s6Yty0JV2bhxI47jTD2PybKsrOlEIjEtnZokPJVW1TPSyfxPpR3HmZZOHTSvtG3b09K5lCk9Pdsy2bbNq1/9akQkMGUq1HFyHGfqLvOglCmVbn8p1VG9NO8yOY7DlVdeSSgUyqlM+eBLIBKR3wWuBu7JsnoVEFfV4Yzlx5LrsnEPMJL2OgLQ2dkJQFdXF11dXQB0dHTQ3d0NQHt7O729vaxYsYKnnnqKvr4+AFpbWxkYGACgpaWFwcFBAJqbmxkedrPV1NTE2Jjblm5sbCQWi2FZFo2NjViWRSwWo7GxEYCxsTGampoAGB4eprm5GYDBwUFaWloAGBgYoLW1FYC+vj7a2toA6Onpob29HYDu7m46OjpmLFNPTw8AbW1tsy7TI488Qm1tLY7jBKZMhTpOTz75JCtWrKC/vz8wZUodp6cOnwDcGlG+ZXr22Wc5deoUoVAopzLlRSrCFuoFrMENKq9JW/YYcF8yfRswmeV7bcDfeeyzEqhNe9UDOjQ0pKqqlmWpZVlnpBOJhEajUX3kkUd0YmJCbdueWp4tHY/Hp6Udx5mWdhznjLSqTkvbtj0tnUgkzpq2LGtaOls5MsuUnp5tmcbHx/Xhhx/WycnJwJSpUMdpYmJCH3nkEY3FYoEpUyKR0KMjUV1798/0kr/4mZ6KJfIuUzQa1Ycffljj8fiMZRoZGVFAgVrNIW6IFvgZRSJyK/BjIH3SknAyUw6wHbdZtlTTakUi0psMVl/O4TdqgZGRkRFqa2vPuq3jOAwPD7NkyZJAdj7OFuPFm6C6eaTzKHf8+9NsWFXDI3fdkPf38/EyOjpKXV0dQJ2qjs60bz8GNP4c2JSx7FvAAeDvgD4gAdwEPAAgIpcDFwF7Cp2ZUCjEsmXLCr3beY/x4k1Q3ZyeCG3JrL7vp5eCh3tVHVPVzvQXMA6cSH4eAb4BfElE3iwi1+AGqj2q+mSh85NIJHjooYemOvEMLsaLN0F1sz81ojrPgYwp/PRSrFs8/hS3mfYAbv/PTuCjfvxQJBJh69atRCKBv5slL4wXb4LoxnaUjiMjQP4DGVP46WVOTKvqmzI+x4A7ky9fEZEZ+5HKEePFmyC6eeHYGBNxm8WVEdadv3hW+/DTS3B64jxIJBLs2LEjcNXsc8V48SaIblIzMr52TR3h0OwmNvPTS+ADUSQS4eabbw5UNbsQGC/eBNFNqn8o10cHZcNPL4EPRECgTqhCYrx4EzQ3s73jPhO/vAQ+EKWPSDWcxnjxJmhuxmIJuo+fAs6tRuSnl4IPaJwL8hnQqKpYlkUkEgn0ZOj5Yrx4EzQ3uw8O8nv/upcLl1bzxN3bZr2ffLzkO6Ax8DUiIDD/2QqN8eJNkNzk+0TXs+GXl8AHIsuyaGpqCtSJVQiMF2+C5iZ1x/25NMvAXy+Bb5oZDOWMqtLwN48yeCrOAx/ZwjVrz71WlAumaZaBqjI6Osp8DLh+Yrx4EyQ3R05GGTwVpyIsvPqCc/un7aeXwAciy7LYtWtXYKrZhcJ48SZIbtqT/UNXrq6lqiJ8Tvvy04tpmhkMAeazP32eb+7u4fYtF/PX73r1nP2uaZpl4DgOQ0NDs5q+MsgYL94EyU17AUZUp/DTS+ADkW3b7Nu3b2reXYOL8eJNUNxMWja/fNmtjBQiEPnpxTTNDIaAsr9vmFvv383ShRU881dvndPBmaZploHjOBw/fjwQ1exCYrx4ExQ3+9PGDxUiCPnppSwCUWdn57w/qQqN8eJNUNwU6kbXFH56MU0zgyGg3Pj3v6D3xATf+cDruOGy8+f0t03TLAPHcejv75/3/90KjfHiTRDcDI3H6T0xAcBrC9BRDf56KYtAdOjQoXl9UvmB8eJNENykJkJbd/4i6qorCrJPP70Ea/anLEQiEW64If9nOAUd48WbILg5/Yz7wt1b5qeXsqgR9fb2zuv/bn5gvHgTBDepWztm+8SObPjppSwC0Xxv7/uB8eLNfHfjOHp6DqIC9Q+5+/XPi7lqZjAEjIPHT/GWLz1OVUWIzr/eTiQ89/UNc9UsA9u2OXjw4Lwfrl9ojBdv5rubVG3oNfVLChqE/PQS+ECkqpw8eTIQc8sUEuPFm/nuZmpGxgL2D4G/XsriqllDQ0Oxs1FyGC/ezHc3p0dULynofv30EvgakW3bHDhwYN5Ws/3CePFmPruJxm0OHB0DYHOBa0R+egl8IAKIRqPFzkJJYrx4M1/d7O05ge0oK2srWV1XXfD9++Ul8E2zcDjM5s2bi52NksN48WY+u/l262EAtr96VcH37aeXwNeIbNums7NzXlaz/cR48Wa+ujl4fIxf/OoVROADb7ik4Pv300vgA5HBUC5844keAG6+ciUXn7eoyLnJj4IHIhG5R0T2iciYiBwXkQdF5PKMbapE5H4ROSEip0TkARFZWei8gFud3LhxI+HwuT3BIGgYL97MRzeDpyZ54Jl+AD609VJffsNPL37UiG4E7gdeD7wVqACaRCQ9RH8Z+DXgt5LbXwD8yIe8YNs27e3t866a7TfGizfz0c139/QStxxeu2YJ1/r0EEU/vRS8s1pV35b+WURuB44D1wAtIlIHfBC4TVWbk9u8H+gSkder6pOFzlN1deGvHgQB48Wb+eQmlrD59yd7Afjw1kt8nZvaLy9z0UdUl3wfSr5fg1tLejS1gaoeAF4Crs+2AxGpFJHa1AuoAaYis23bWdOpB8Ft2LABVZ26Wc+yrKzpRCIxLZ0aQZpKq+oZ6WT+p9KO40xLp/LglbZte1o6lzKlp2dbJsdxuPzyywmFQoEpU6GOk6qyYcMGRGRelOmHT73EifE49Uuquemy5b4dJ1XlVa96FeFwOKcy5YOvgUhEQsB9wG5V7UwuXgXEVXU4Y/NjyXXZuAcYSXsdAejsdHfZ1dVFV1cXAB0dHXR3dwPQ3t7OoUOH2LdvH3v37qWvrw+A1tZWBgYGAGhpaWFwcBCA5uZmhofdbDU1NTE25g4Ma2xsJBaLYVkWjY2NWJZFLBajsbERgLGxMZqamgAYHh6mubkZgMHBQVpaWgAYGBigtbUVgL6+Ptra2gDo6emhvb0dgO7ubjo6OmYsU0+P2ynZ1tZ2TmXau3fvVDmCUqZCHKfdu3ezb98+Dh8+XPJlGhkd5Z+angfgdzavoOXxx3w7Tk8//TQtLS1YlpVTmfLB17vvReRrwNuBN6rqkeSy24BvqWplxrZtwC9U9e4s+6kE0revAY4MDQ2xdOnSqYgdDoenpVPR+aWXXmLNmjVUVFQQCoWwLItQKHRGOpFIEA6Hp9KRSAQRmUqDG/HT0xUVFajqVNpxHGzbnko7jkMkEvFM27aNqk6ls5Ujs0wiMpXOVo5cyhSLxejr6+OSSy6Z+v35XqZCHadEIkFfXx9r165FREq6TM1dx/jAt5+ipjLCE3e/iYUVId+OUzwe5/Dhw6xbtw5VPWuZTp06ldfd974FIhH5CvBu4AZV7Ulbvg34ObA0vVYkIr3Afar65Rz2baYBMRiA277+JK2HTvDhrZfwl++4stjZmaLo04CIy1eA9wDb0oNQkqeBBHBT2ncuBy4C9hQ6P5Zl0draOtUWNrgYL97MFze/fHmE1kMnCIeE230YwJiJn178uMXjfuA23NrQmIik+n1GVDWqqiMi8g3gSyIyBIwC/wTs8eOKWSgUor6+nlDIjN1Mx3jxZr64+cYu93/8LZtWU7/E/6t8fnopeNNMRLx2+H5V/bfkNlXA/wH+B27fz07go6p6NMffME0zQ1lzdCTGG/+uGctRfvKxN/CaC5cUO0vTKHrTTFXF4/VvadvEVPVOVV2mqotU9ddzDUL5YlnWVE+/4TTGizfzwc239xzGcpTXXbxszoKQn15Ku+5ZAEKhEOvWrSv5avZcY7x4U+puxictvpccwPihrf73DaXw00vgpwFJtWsN0zFevCl1Nz94qo/RmMXFyxfylit8uUUzK356Kc2QX0Asy6K5ubmkq9nFwHjxppTd2I7yzd2HAfjgGy8hFPLvdo5M/PQS+EAUCoXYuHFjyVazi4Xx4k0pu/nv54/y0tAESxZW8BvXXDinv+2nl7Jomq1YsaLY2Sg5jBdvStnNvyYv2f/edRexcMHc/vn66aX0Qn6BSSQS7Ny5c+omQYOL8eJNqbppf+kkT/WeZEE4xPuuv3jOf99PL4EPROFwmIaGhnk1ydVcYLx4U6puUrWhd111AStqq+b89/30UhZNs2XLlhU7GyWH8eJNKbrpG5rg4U73DvcPvnHuLtmn46eXwNeIEokEDz30UMlVs4uN8eJNKbr51u7DOApbX3UeV6wuzt0EfnoJfCCKRCJs3bp1avoEg4vx4k2puRmNJfj+vpeA4tWGwF8vpWHaR0TE3I+WBePFm1JyY9kOf/5fzzIet7ls5WJuvOz8ouXFTy+BrxElEgl27NhRUtXsUsB48aZU3Kgq9/zoOZqeP8aCcIjPvnujr/NRz4SfXnydodEv8rn7XlWJxWJUVVUV9SCWGsaLN6XgRlX528Yuvr6rh5DA137/Gl+e3ppvnnL1UvS770uRUmnrlxrGizfFdvPVxw7x9eTl+i/8xmuKHoRS+OUl8IEofdJxw2mMF2+K7eZ7e3v5+52/AuAvb7mC3752TVHykYmfXsqiaZaadNw0QU5jvHhTTDc/63iZP/6PdlTho29ax6fetmFOf/9s5OPFNM2yYP7rZ8d48aYYbh5/4RX+9Pv7UYXbrruIT26/fOYvzTF+eQl8ILIsi6amJvNHl4Hx4k0x3Dzde5I7vvs0CVt552tW87kiXyHLhp9eAt80MxhKnV8dHeO3/+8eRqIJbrjsfP71vdeyIDK/6wimaZaBqjI6Osp8DLh+Yrx4M5duXjoxwR98Yy8j0QRXX7SEf/79q0s2CPnppTRLXEAsy2LXrl2mCZKB8eLNXLk5PhbjD765l+Njk1y+soZv3t4w53MM5YOfXkzTzGAoAsfHYrz3G20cODrGmmXVPHDHlqJM7eEXpmmWgeM4DA0N4ThOsbNSUhgv3vjpZmQiwRcfOcCNX3yMA0fHOL+mkn//4HXzIgj56SXwgci2bfbt24dt28XOSklhvHjjh5vxSYuvNHfzxi8289XHDhFN2Fy1Zgnf+9B1rF2+qGC/4yd+njOmaWYw+EgsYfO9vS/x1V8c5MR4HIANq2r485sv56YrVpTcJfpCkW/TrHR7xgqE4zgMDg5y3nnnleRTGYqF8eJNIdwkbIcfPn2Ef/x5NwMjMQAuXr6QT9x8Oe/ctHpOHwNUKPw8ZwJ/BjqOQ2dnp+kLycB48eZc3DiOsmN/P2/90uPc86PnGBiJsbquii/8+ib++xM38q7XXjAvgxD4e86YppnBcI7YjvL8y6M8+eIJHnjmCAeOjgGwfNEC7nzzem677iKqKkprIn6/MU2zDBzHYWBggNWrV5smSBrGizczuXEc5fkBN/A8+eIJ9vYMMRY7PbampirCHTeu4/YtF7OoMjh/Yn6eM8Gx5IHjOBw6dIiVK1eaP7g0jBdvMt04jnLg6BhPvniCPS+eoK1niJHo9FkKayojvO6SZVy/bjm/dc0a6hZWFCn3/uHnOVO0ppmI3Al8ElgFPAv8saq25fhd0zQzFJRYwmZgJMaRkxP0n4zSPxyl/2SUI8NRXjg2xvDE9MCzuDJCw8VLef2ly7l+3XJefUEd4Xna9+MH86JpJiK/A3wJuAPYC9wF7BSRy1X1eCF/y3Ec+vr6WLNmjfnPn0aQvagqsYTDWCzBaMxiLJZgLGYlX256cHySIyejU0HnlbHJs+5z4YIw1168jOsvXc7rL13Gpvo6IuFgeZsJP8+ZYjXNPgF8XVW/BSAidwDvAD4AfKFQP3Li1CQ9r4zx/PO9XOksKuqTO/Opd+ZXST1zY6/va9p627Z5vuswV8SrCYfD09a52+rUFzS5XNHke2pbddPJdY4DTnKZquKo+z1HFSe5Y0cV2wHbcdLeFctxt7EcxXFOv8dtZdKyiVsOccthMvket1Of7anlE3F7KtBYTv41/eqKMPVLq7mgroqK+CibL1vLmuWLWLt8Ea++oJaKMgs8mTiOQ39/P/X19fO/aSYiC4AJ4DdV9cG05d8Glqjqu7N8pxKoTFtUAxwZGhpi6dKlUyM9w+HwtPT323q5+0edvpXFUNqExG1C1VRFqKmqoKYqwuLKMHXVC1hSXcHqukouWr6Y1XVVXFBXyfm11W4AdRwikQiO40ylbdtGVafScOb5lp62LAsRmUqHQiFCoZBnOpFIEA6Hp9KpWRBTaWBqdsRUuqKiYmrWxIqKChzHwbbtqXS2csxVmU6dOlXyTbPzgDBwLGP5McBrXsx7gP+VubCzs5OtW7fS1dUFwMaNG+no6KC6upoNGzZwrP8lVtW4nYa2bRMJh4lURIhFY1RUVBCOhIlGYyxYsIBwOER0IkplZSWhcIiJiShVVVWEQsLE+ATV1dUgQnRiguqFC0GVaDTKwkULcRz36QYLF1bj2A6Tk5NUL6zGth3i8TjV1VXYlo1lJaiqqsayEliWTVVVFYlEAsexqaysIpGI4zhKZWUlibg7CndB5QImJycRERYsWEAsNkk4HKKiooJYNEYkEqGiIkI0GqWiooJIJMJENErlggVEIhHGx8epqqoiHA4zPj5OdXU14eSJsmDBAioqKpiYGGfRosWAMj4+Tk1NDY5tMxGNUluzGNu2iUaj1NbUYFkWk5MxahbXkEgkiMcnqalZTCKeIJGIU1tTQ3wyhm3b1NYsJhaNoupQs7iGiYlxBKirreHU2CgLImGW1NUyPHyS6soFLKmt5cTgK9QsXkRdbQ3Hjg6wfOkSli+ppb+vlwtXr2JpXQ0vdv+K9ZdewrK6Gp57dj9XbbqSJbWLeerJVrbdsIWli6rY09LMe975dmzbpqmpiXe/+92Mjo6ya9cu3vGOdzA0NMS+ffvYvmk7x48fp/PpdrZt28bLL7/MwYMHqa+vJxwOMzAwwJYtW+jp6eHkyZM0NDTQ3d1NNBpl8+bNnudee3s7S5cuZf369bS1tVFfX8/atWtpbW1l3bp11NfX09LSwsaNG1mxYgXNzc00NDSwbNkympqa2Lp1K7W1tTQ2NnLzzTcTiURobGzklltumZqg7N3vfjdjY2NTZRoeHnbLtH07g4ODdHZ2sm3bNgYGBjh06BA33HADfX199Pf3z6pMzzzzDI7j0NDQMGOZqqryu3euGDWiC4B+YIuq7klb/kXgRlW9Lst3ZlUjsiwL27bp6Ohg06ZNLFiwoKz+K52tTNFolOeee47NmzcDBKJMhTpO8Xic5557jte+9rWEQqFAlKkQx2lycpJnn32Wq6++GqCgNaJ50TTLsg9z1cxgKGFKfhoQVY0DTwM3pZaJSCj5eY/X92aLbdscOHDA3GWegfHijXGTHT+9FOuq2ZeAb4vIU0Ab7uX7RcC3/PixaDTqx27nPcaLN8ZNdvzyUswBjR/j9IDG/cDHVXVvjt81TTODoYQp+aZZClX9iqquVdVKVb0u1yCUL7Zt09nZaarZGRgv3hg32fHTS3mP0DIYDCWBmQbEYDAUnHlxr1mhGB2dsXxT1cmNGzcW9RaPUsN48ca4yU4+XnL520xnvtaI6oEjxc6HwWCYkQtVtX+mjeZrIBLgAmAsbXEb8Los6RrcoHVhxvbnSvpvFGr7s22TbV0uy7w+l4qXXL/jtU0+y8/mxpwzM3/O10sN8LLmEGTmZdMsWbBpUVZEnFRbNCOd2mQsl7ZqrqT/RqG2P9s22dblsszrc6l4yfU7Xtvks/xsbsw548s5k3NZg3TV7H6P9Fz8XqG2P9s22dblsmymz4VmNvs/Fzf5LD+bC3PO5P654MzLplk+pK6wkWPvfblgvHhj3GTHTy9BqhF5MQl8JvluOI3x4o1xkx3fvAS+RmQwGEqfcqgRGQyGEscEIoPBUHRMIDIYDEXHBCKDwVB0TCAyGAxFxwSiJCKyRkQeE5HnRaRDRH6r2HkqJUTkxyJyUkR+WOy8FBMReaeI/EpEukXkQ8XOTylxLueIuXyfRERWAytVdb+IrMKdV/syVR0vctZKAhF5E+69Q+9T1d8sbm6Kg4hEgOeBN+MO7Hsa92k0J4qasRLhXM4RUyNKoqoDqro/mT4KDALLipqpEkJVH6OwN4DOR14H/FJV+1X1FPAwcHOR81QynMs5Mm8CkYjcICI/FZGXRURF5NYs29wpIodFJCYie0Uk37vAU/u5Bgirat+55nsumEs385kCeEo9ky9FP1Dvc7bnhGKfQ/MmEOE+5eNZ4M5sK0Xkd3CfDvIZ4OrktjtFZEXaNvtFpDPL64K0bZYB3wH+0MeyFJo5cRMAztlTgCmuG1Wddy9AgVszlu0FvpL2OYT7H+sv8thvJdAC/EGxy1hqbpLfexPww2KXsViegC3Aj9PW3wfcVuyylIKbcz1H5lONyBNxnx57DfBoapmqOsnP1+e4DwH+DWhW1e/6kM2iUAg35UCOntqAjSJSLyKLgbcDO+c6r3PNXJxDgQhEwHlAGDiWsfwY7nPTcuENwO8AtyabKftFZFMB81gsCuEGEXkU+AFwi4gcEZGgBbEZPamqBfwZ8AvcZ/H9Hy2PK2Y5nUPnco7Myxka/UBVnyA4gbngqOpbip2HUkBVfwL8pNj5KEXO5RwJyh/eIGADKzOWrwSOzn12SgrjJjeMJ298dxOIQKSqcdzBZTellolIKPl5T7HyVQoYN7lhPHkzF27mTdMs2Tm4Pm3RJSJyFTCkqi/hXlr8tog8hdupeBfuJclvzXFW5xzjJjeMJ2+K7qbYlwrzuKT4JtzLipmvf0vb5mNAL+5UlnuB64qdb+OmdF7GU+m6MfeaGQyGohOIPiKDwTC/MYHIYDAUHROIDAZD0TGByGAwFB0TiAwGQ9ExgchgMBQdE4gMBkPRMYHIYDAUHROIDAZD0TGByDCniMjFyTmRryrAPlRE9p9jfm5P29d957Ivw+wxgajMEZHzReRrIvKSiEyKyFER2Skib0jbJutk6rOkD1gNdBZgX29h+h3hbxWRF0RkVES+m5xZMLWuLrlubcY+vp/MT1nfYV9sTCAyPABsBt4HXAa8C3gMWJ7PTtL/6M+GqtqqelTd2Q7PlROanCExOS3F/wP+GXf60muZ/gCELwD/rKq9GfmJqvv4qHgB8mOYJSYQlTEisgTYCtytqr9Q1V5VbVPV/63uTISIyOHk5j9O1owOJ5f/dXI63Q+JSA8QSy5/m4g8ISLDInJCRH4mIuvSfnNa00xE3pT8fJOIPCUiEyLSKiKX51mc85Kvr6rqL3FnUbwi+RtbgAbgH2ahyTAHmEBU3pxKvm4VkUqPbRqS7+/HbcI0pK1bD/wG8OvAVclli3DnrrkWt9nk4Aaxmc61v8GdD/pawAK+mU9BgFeAAeBmEVmIG2A7RKQC+BrwR6pq57lPwxxhAlEZk2we3Y7bLBsWkd0i8rci8pq0bV5JJoeTTapX0naxAHivqrarakdy+wdU9UeqelDdJ+d+ANgEXDlDdv5SVR9X1edxm1FbRKQqj7Io8NvAXwG/BNpxg9lf4E52H0uW71ci8rFc92uYG0wgKnNU9QHcJ5i+C3gEd4KsZ0Tk9hy+3psRmBCRV4nIf4jIiyIyChxOrrpohn11pKUHku95PbxPVZ9Q1QZVvURV7wQuAd4L3At8F/gX3JrSp9ODraH4mEBkQFVjqvrfqvo5Vd2C+3y3z+Tw1fEsy34KLAM+DFyXfIFbezobifQsJd/P9fz8v7jNvRBuh/wPVPU48Dhw4znu21BATCAyZON53L6eFAnc51qdFRFZDlwOfF5Vf66qXcBSf7I4Y14+iDvf8k84nfeKtPcZy2OYO+bN5PmGwpMMHD/A7UvpAMZwO4s/BexI2/QwcJOI7AYmVfWkxy5PAieAPxSRAdzm2Bf8yb03yeex34v70ExU9aSIdAF3iUgTbif638x1vgzemBpReXMKdxL0PwVacAcZfg74Ou5E6Sn+DHgr7mDEdq+dqfsY4t/FfTxxJ/Bl4JN+ZHwG/gH3Kawvpy27HTdvPwP+XlX3FSFfBg/M5PmGeYeIXAz0AJuTV+YKsc/HgP2qelch9mfID1MjMsxnWkWk9Vx2ICK/JyKncK+mGYqEqREZ5h0iEgEuTn6cVNW+c9hXDacfpTysqoPnmD3DLDCByGAwFB3TNDMYDEXHBCKDwVB0TCAyGAxFxwQig8FQdEwgMhgMRccEIoPBUHRMIDIYDEXHBCKDwVB0/j/NfIyV9cAxHwAAAABJRU5ErkJggg==\n",
+      "text/plain": [
+       "<Figure size 300x300 with 1 Axes>"
+      ]
+     },
+     "metadata": {
+      "needs_background": "light"
+     },
+     "output_type": "display_data"
+    }
+   ],
+   "source": [
+    "damping_curve.plot();"
    ]
   },
   {
@@ -179,7 +203,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 10,
+   "execution_count": 11,
    "metadata": {},
    "outputs": [
     {
@@ -196,12 +220,16 @@
     }
    ],
    "source": [
-    "HH_x_param = damping.get_HH_x_param(pop_size=30, n_gen=3, show_fig=True);"
+    "HH_x_param = damping_curve.get_HH_x_param(\n",
+    "    pop_size=30,\n",
+    "    n_gen=3,\n",
+    "    show_fig=True,\n",
+    ");"
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": 11,
+   "execution_count": 12,
    "metadata": {},
    "outputs": [
     {
@@ -210,7 +238,7 @@
        "PySeismoSoil.class_parameters.HH_Param"
       ]
      },
-     "execution_count": 11,
+     "execution_count": 12,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -228,7 +256,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 12,
+   "execution_count": 13,
    "metadata": {},
    "outputs": [
     {
@@ -247,7 +275,7 @@
        "}"
       ]
      },
-     "execution_count": 12,
+     "execution_count": 13,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -265,7 +293,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 13,
+   "execution_count": 14,
    "metadata": {},
    "outputs": [
     {
@@ -282,12 +310,16 @@
     }
    ],
    "source": [
-    "H4_x_param = damping.get_H4_x_param(pop_size=30, n_gen=3, show_fig=True);"
+    "H4_x_param = damping_curve.get_H4_x_param(\n",
+    "    pop_size=30,\n",
+    "    n_gen=3,\n",
+    "    show_fig=True,\n",
+    ");"
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": 14,
+   "execution_count": 15,
    "metadata": {},
    "outputs": [
     {
@@ -296,7 +328,7 @@
        "PySeismoSoil.class_parameters.MKZ_Param"
       ]
      },
-     "execution_count": 14,
+     "execution_count": 15,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -314,7 +346,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 15,
+   "execution_count": 16,
    "metadata": {},
    "outputs": [
     {
@@ -328,7 +360,7 @@
        "}"
       ]
      },
-     "execution_count": 15,
+     "execution_count": 16,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -346,16 +378,27 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 16,
+   "execution_count": 17,
    "metadata": {},
    "outputs": [],
    "source": [
-    "stress_curve = Stress_Curve(data)"
+    "strain = np.geomspace(1e-2, 10, num=20)\n",
+    "y = strain ** 2\n",
+    "fictitious_stress_data = np.column_stack((strain, y))"
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": 17,
+   "execution_count": 18,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "stress_curve = Stress_Curve(fictitious_stress_data)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 19,
    "metadata": {},
    "outputs": [
     {
@@ -368,7 +411,7 @@
        "        335.98182863,  483.29302386,  695.19279618, 1000.        ])"
       ]
      },
-     "execution_count": 17,
+     "execution_count": 19,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -379,20 +422,20 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 18,
+   "execution_count": 20,
    "metadata": {},
    "outputs": [
     {
      "data": {
       "text/plain": [
-       "array([0.00000000e+00, 1.06913808e-04, 3.28133240e-04, 7.85866790e-04,\n",
-       "       1.73298071e-03, 3.69269019e-03, 7.74759970e-03, 1.61377674e-02,\n",
-       "       3.34981829e-02, 6.94192796e-02, 1.43744989e-01, 2.97535144e-01,\n",
-       "       6.15748211e-01, 1.27417499e+00, 2.63655090e+00, 5.45549478e+00,\n",
-       "       1.12882789e+01, 2.33571147e+01, 4.83292024e+01, 9.99999000e+01])"
+       "array([1.00000000e-04, 2.06913808e-04, 4.28133240e-04, 8.85866790e-04,\n",
+       "       1.83298071e-03, 3.79269019e-03, 7.84759970e-03, 1.62377674e-02,\n",
+       "       3.35981829e-02, 6.95192796e-02, 1.43844989e-01, 2.97635144e-01,\n",
+       "       6.15848211e-01, 1.27427499e+00, 2.63665090e+00, 5.45559478e+00,\n",
+       "       1.12883789e+01, 2.33572147e+01, 4.83293024e+01, 1.00000000e+02])"
       ]
      },
-     "execution_count": 18,
+     "execution_count": 20,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -410,7 +453,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 19,
+   "execution_count": 21,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -435,7 +478,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 20,
+   "execution_count": 22,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -453,7 +496,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 21,
+   "execution_count": 23,
    "metadata": {},
    "outputs": [
     {
@@ -462,7 +505,7 @@
        "n_layers = 3, type: <class 'PySeismoSoil.class_curves.Damping_Curve'>"
       ]
      },
-     "execution_count": 21,
+     "execution_count": 23,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -473,7 +516,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 22,
+   "execution_count": 24,
    "metadata": {},
    "outputs": [
     {
@@ -509,7 +552,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 23,
+   "execution_count": 25,
    "metadata": {},
    "outputs": [
     {
@@ -518,7 +561,7 @@
        "3"
       ]
      },
-     "execution_count": 23,
+     "execution_count": 25,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -529,7 +572,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 24,
+   "execution_count": 26,
    "metadata": {},
    "outputs": [
     {
@@ -538,7 +581,7 @@
        "True"
       ]
      },
-     "execution_count": 24,
+     "execution_count": 26,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -549,7 +592,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 25,
+   "execution_count": 27,
    "metadata": {},
    "outputs": [
     {
@@ -558,7 +601,7 @@
        "False"
       ]
      },
-     "execution_count": 25,
+     "execution_count": 27,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -569,7 +612,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 26,
+   "execution_count": 28,
    "metadata": {},
    "outputs": [
     {
@@ -593,7 +636,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 27,
+   "execution_count": 29,
    "metadata": {},
    "outputs": [
     {
@@ -623,7 +666,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 28,
+   "execution_count": 30,
    "metadata": {},
    "outputs": [
     {
@@ -653,7 +696,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 29,
+   "execution_count": 31,
    "metadata": {},
    "outputs": [
     {
@@ -683,7 +726,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 30,
+   "execution_count": 32,
    "metadata": {},
    "outputs": [
     {
@@ -705,7 +748,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 31,
+   "execution_count": 33,
    "metadata": {},
    "outputs": [
     {
@@ -714,7 +757,7 @@
        "PySeismoSoil.class_parameters.HH_Param_Multi_Layer"
       ]
      },
-     "execution_count": 31,
+     "execution_count": 33,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -732,7 +775,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 32,
+   "execution_count": 34,
    "metadata": {},
    "outputs": [
     {
@@ -782,7 +825,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 33,
+   "execution_count": 35,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -791,7 +834,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 34,
+   "execution_count": 36,
    "metadata": {},
    "outputs": [
     {
@@ -809,7 +852,7 @@
        "       [   0.2637, -999.    ,    0.2637,    0.198 ,    0.2637, -999.    ,    0.2637,    0.1033]])"
       ]
      },
-     "execution_count": 34,
+     "execution_count": 36,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -835,7 +878,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 35,
+   "execution_count": 37,
    "metadata": {},
    "outputs": [
     {
@@ -872,7 +915,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 36,
+   "execution_count": 38,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -888,7 +931,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 37,
+   "execution_count": 39,
    "metadata": {},
    "outputs": [
     {
@@ -924,7 +967,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 38,
+   "execution_count": 40,
    "metadata": {},
    "outputs": [
     {
@@ -942,7 +985,7 @@
        "       [   3.    ,    0.0078,    3.    , -999.    ,    3.    ,    0.0138,    3.    , -999.    ]])"
       ]
      },
-     "execution_count": 38,
+     "execution_count": 40,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -968,7 +1011,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 39,
+   "execution_count": 41,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -984,7 +1027,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 40,
+   "execution_count": 42,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -993,7 +1036,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 41,
+   "execution_count": 43,
    "metadata": {},
    "outputs": [
     {
@@ -1002,7 +1045,7 @@
        "5"
       ]
      },
-     "execution_count": 41,
+     "execution_count": 43,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -1020,7 +1063,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 42,
+   "execution_count": 44,
    "metadata": {},
    "outputs": [
     {
@@ -1063,7 +1106,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 43,
+   "execution_count": 45,
    "metadata": {},
    "outputs": [
     {
@@ -1081,7 +1124,7 @@
        "       [ 3.    ,  0.0078,  3.    , 21.052 ,  3.    ,  0.0138,  3.    , 20.947 ]])"
       ]
      },
-     "execution_count": 43,
+     "execution_count": 45,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -1100,7 +1143,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 44,
+   "execution_count": 46,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -1109,7 +1152,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 45,
+   "execution_count": 47,
    "metadata": {},
    "outputs": [
     {
diff --git a/setup.py b/setup.py
index 1cfd719..630b90e 100644
--- a/setup.py
+++ b/setup.py
@@ -16,7 +16,7 @@ def load_requirements(req_filename):
 
 setup(
     name='PySeismoSoil',
-    version='v0.5.0',
+    version='v0.5.1',
     description='PySeismoSoil',
     long_description=long_description,  # will get passed to PyPI.org
     long_description_content_type="text/markdown",  # otherwise the default is RST