Skip to content

Instantly share code, notes, and snippets.

@vuddameri
Created March 24, 2024 08:42
Show Gist options
  • Select an option

  • Save vuddameri/76afc03ba6cf3af3be6ecb1ffbe8ae87 to your computer and use it in GitHub Desktop.

Select an option

Save vuddameri/76afc03ba6cf3af3be6ecb1ffbe8ae87 to your computer and use it in GitHub Desktop.
Rectangular-Load.ipynb
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "markdown",
"metadata": {
"id": "view-in-github",
"colab_type": "text"
},
"source": [
"<a href=\"https://colab.research.google.com/gist/vuddameri/76afc03ba6cf3af3be6ecb1ffbe8ae87/rectangular-load.ipynb\" target=\"_parent\"><img src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open In Colab\"/></a>"
]
},
{
"cell_type": "markdown",
"id": "1c98c236",
"metadata": {
"id": "1c98c236"
},
"source": [
"<h1><font color='blue'><center> CVEN 2320 Differential Equations for Civil Engineers </center></h1>"
]
},
{
"cell_type": "markdown",
"id": "cdee0e99",
"metadata": {
"id": "cdee0e99"
},
"source": [
"<center><font color='teal'><h2> Venki Uddameri, Ph.D. P.E., </h2>"
]
},
{
"cell_type": "markdown",
"id": "daa40516",
"metadata": {
"id": "daa40516"
},
"source": [
"<hr>"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "d16437c5",
"metadata": {
"id": "d16437c5"
},
"outputs": [],
"source": [
"import numpy as np\n",
"from matplotlib import pyplot as plt"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "f44fc4b5",
"metadata": {
"id": "f44fc4b5"
},
"outputs": [],
"source": [
"r = 1\n",
"a = 0.2\n",
"k = 0.35\n",
"eta = 0.05/24\n",
"td = 3\n",
"So = 20"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "1fad74a4",
"metadata": {
"id": "1fad74a4"
},
"outputs": [],
"source": [
"\n",
"u = lambda t,td: np.where(t<td,0,1)"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "d1be9e73",
"metadata": {
"id": "d1be9e73"
},
"outputs": [],
"source": [
"t = np.arange(0,2400,0.25)\n",
"num1 = r/eta*u(t,0)\n",
"num2 = r/eta*np.exp(-eta*t)\n",
"num3 = r*u(t,td)*(1/eta - 1/eta*np.exp(-eta*(t-td)))\n",
"num4 = a*np.exp(-k*t)\n",
"num5 = a*np.exp(-eta*t)\n",
"num6 = So * np.exp(-eta*t)\n",
"S = num1 - num2 - num3-num4 - num5 + num6"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "d8cdfb09",
"metadata": {
"id": "d8cdfb09",
"outputId": "a4b41e28-724a-4cc6-93df-130af6513ee8"
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAEWCAYAAACT7WsrAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAsR0lEQVR4nO3dd5hU5fn/8fdN79J0pSlgAVERha+9AEqwV4wtsURFo4k1X+s3Br1siT32+kOjghobsTcQg7GAoiBKUERBQEBpa0HK/fvjedYdl53Z2WVnz87M53Vd55o5/X7mnJl7znPOeY65OyIiIpVpkHQAIiJSfylJiIhIWkoSIiKSlpKEiIikpSQhIiJpKUmIiEhaShJSJTPbxcxmmFmpmR2cdDx1wcyON7N/19KyupuZm1mj2ljeuoqxbJrg+kvNrGeW02a975nZQDObk9I/y8z2WsdwK66jXm3LuqAkUUHcsX6IO+V8MxtpZq2SjithlwG3uHsrd3+qpgsxs07xC1aSMuziNMNeyGJ5tf4jkIR1LYeZjTWzhWa2zMw+MLODajO+StY3wsxWxu/IEjN708x2ynb+uB/NzHLyWtn3Korf659iGb41s5fNrHdtLb+QKElU7gB3bwX0A7YFLkw2nMRtDHy0rgtx93nAp8DuKYN3Bz6pZNj4dV1fJhYUyv5/JtDJ3dsAw4EHzaxTbSw4wz/mR+J3pCMwFnisNtZXiVrZ99L4WyxDV2ABMDJH68lrhfIlyQl3nw+8SEgWAJjZjvGf05L4r21gyrjjzWymmS03s8/N7JiU4RPM7GYzW2pmn5jZninzdTazMfEfzadmdnLKuBFm9qiZPRCX+5GZDUgZf76ZfRXHTS9brpk1MLMLzOwzM/smLqN9urKa2clx3d/GWDrH4Z8BPYF/xX9dTSuZt2w9y81smpkdkuFjHU9MCGbWkJCEb6owbCdgvJltYmavxfgXmdlDZtY2TvcPYKOUuM7LYvuMM7MrzGwC8D3QM902S5nnWjNbHMftkzJ8PTO718zmxc//8hg7ZtYwzrfIzGYC+2X43NOV48C4rZfEuLdItwx3/9DdV5X1Ao2Bbhm2AcBeFqpxFpvZrWZmcb1l++oNZvYtMCLTQuJ6HwK6mNn6cRnbm9l/YuzzzOwWM2uSUuafq7ss/KO/1cyejdvgbTPbJI5ba98zsxPM7OM47UwzO6WKclbJ3b8HHga2iuv9xZFd/A4+WNm8mfYfM/tdjHWxmb1oZhuva6yJcHd1KR0wC9grvu8KTAFuiv1dgG+AfQkJdkjsXx9oCSwDesVpOwFbxvfHA6uAswlf4COApUD7OP514DagGSEhLQT2jONGAD/GdTYErgLeiuN6AbOBzrG/O7BJfH8W8FYsQ1PgTmBUmjIPBhYB28VpbwbGV/aZpJn/cKBz/EyOAL4j/LOtbNrjgA/i+wGEpLFZhWE/AE2ATeNn3DR+xuOBG9PFlWn7xPHjgC+BLYFGwHpVbLOVwMnxc/89MBewOP6p+Jm2BDYA3gFOieNOJRwddQPaE/5pO9Coqn0u9m8eP8MhhP3lPMIRWJMM2+CZuJ848ALQIMO0HqdvS0hQC4G9K+yrf4yfUfNK5h8BPBjfNwGujvtPozisP7BjnL878DFwVoX1bxrfjwS+BbaP0z8EjM7w2ewHbAIYsAch2W8Xxw0E5mSz38b1Xh7ftyIkiTfSrDO1vN3LtiWZv/MHx222RZz2/4A3k/59q0mXeAD1rYs7SCmwPO4MrwJt47jzgX9UmP5Fwg9fS2AJcFjFL1b84v38AxOHvQP8lvBDshponTLuKmBkyg76Ssq4PsAP8f2mhMPkvYDGFdb5MTHRxP5OhB+9tX6ogHsJh95l/a3itN1TPpO0SaKS5U0GDkozrnssbztC0rwiDv8qZdjYNPMeDLxfYVulfpnTbp/4fhxwWcq4qrbZpyn9LeL+sCFQAqxInQc4qixu4DXg1JRxv6J6SeLPwKMp/Q3i5zOwis+9MbAPcHYV0zmwa0r/o8AFKeX+sor5RwA/xc9uNSERp42N8IflyQrrT00S96SM2xf4JN1nU8mynwLOjO8HUr0k8WMsw3xgDOV/sCpujxGkTxLp9p/ngRMrbMPvgY2z/R7Vl07VTZU72N1bE3a63oR6Vwj1o4fHw+glZrYE2JXwr/k7wr/oU4F58fA59UTYVx73lugLwr/vzsC37r68wrguKf3zU95/DzQzs0bu/inhCzgCWGBmo8uqiWKsT6bE+THhC13C2jrHdQLg7qWEL36XSqZdi5kda2aTU9a1FeWf2S+4+yxgDuFz2x14I476T8qw8XG5G8QyfWVmy4AH0y03pcyVbp+UaWanxFLVNpufMu338W2ruJ7GcZ6y9dxJOKKA8HnOTlnOF1RPxe2xJi4v4/Zw95Xu/jww1MwOBIhVVqWx262yshH2qdSLM1JjT+dRd29L2J+mEo4eiOvc3MyesXDhxzLgSjJvt0yx/IKZ7WNmb1moFl1CSCqZlp3Jte7e1t03dPcD3f2z6sxcxf6zMXBTyv7xLeHoJ6vvVH2iJJGBu79O+MdxbRw0m/BPtW1K19Ldr47Tv+juQwg/Sp8Ad6csrktZvW+0EeHoYi7Q3sxaVxj3VZYxPuzuuxJ2Sgf+mhLrPhVibebulS13bpwfADNrCXTIJoZYz3o38AegQ/zhmEr4QqTzBiEZ7AS8WWHYrpSftL4qlqmvh5Oyv6mw3NSkC1Vsn8rmqWKbpTObcCTRMWU9bdx9yzh+Hr88J7BRFcurWI6K28Pi8rLaJwj/cjcBcPctPVwZ1Mrd36hivnTxpJ/QfRFwCjDCyk+W3074LDeL2+0iMu8PWbFwPuxxwvexJO5rz9XGsiv4jnDkWGbDdBNm2H9mE6ofU/fF5u7+Zrpl1VdKElW7ERhiZv0I/2QPMLOh8eRkMwvXZnc1s5J4srEl4QeklPDPvcwGwBlm1tjMDifUVT7n7rMJP5RXxeX1BU4k1M1mZGa9zGxw/PL8SKjLL1vnHcAVZSfLzGx9S39p5MPACWbWLy7rSuDt+K+/Ki0JPyoL43pOIJ4AzGA8cCww192XxWH/jsPWIxxVALQmfI5LzKwL8L8VlvM14cRmmbTbp7IgsthmlfJwldZLwHVm1sbCRQKbmNkecZJHCdu6q5m1Ay6oYpEVy/EosJ+Z7WlmjYFzY3xr/cCYWe/477p53Ld+Q0i2r1dVjtri7p8QqvXOi4NaE+rqS+M/69/X0qqaEM5PLQRWWbiQ4Fe1tOxUk4Ej4+c5ABhW2URV7D93ABea2ZZx2vXi9z7vKElUwd0XAg8Af44/6AcR/hktJPxb+F/C59iA8GWeSzi03AM4LWVRbxNO0C4CrgCGufs3cdxRhLrOucCTwF/c/eUswmtK+UnD+YREdFEcdxOhnvUlM1tOOIm9Q5oyvkqoB3+c8C94E+DILNaPu08DriP8sH8NbA1MqGK212OsqTerTQaaA5NSqnYuJZxMXwo8CzxRYTlXAf8XD+n/VMX2qUxV2yyTYwk/WtOAxcA/Ka/Wupvwo/kB8F4lcVdUsRzTCUdNNxO27QGEy7J/qmReI1Y3Esp8JnCEu7+XZTlqyzXAcDPbAPgTcDThvN7dwCO1sYJYJXsGIYkujusYUxvLruDPhO/AYsI++HCa6dLuP+7+JOGofnSscptKOF+Ud8qu1JAcMrPjgZNitZCISN7QkYSIiKSlJCEiImmpuklERNLSkYSIiKSVF83dduzY0bt3716jeb/77jtatmxZuwHlkWIuv8penGWH4i5/atknTZq0yN3XX5fl5UWS6N69OxMnTqzRvOPGjWPgwIG1G1AeKebyq+wDkw4jMcVc/tSym1l17/Zfi6qbREQkLSUJERFJS0lCRETSUpIQEZG0lCRERCQtJQkREUlLSUJERNIq6CSxejU880wnSkuTjkREJD8VdJJ47z247rpeHHNM0pGIiOSngk4Sq1aF1zFjQO0YiohUX0EniebNy9/PmJFcHCIi+aqgk0S/fvDss28wfz5svnnS0YiI5J+CThIALVqspqQk6ShERPJTQSeJCRNg+PD+jBoFQ4fCf/+bdEQiIvmloJPEkiUwY0ZrmjaFl16CJ55IOiIRkfxS0Emi7OqmHj1g++3h8ceTjUdEJN8URZJo2BCGDYOJE2HWrERDEhHJKwWdJFavDq+NGsFhh4X3OpoQEcleQSeJjh2hX7/FtGoFPXvCscdCp05JRyUikj/y4hnXNTV4MNxwwwdstNFAAO6/P9l4RETyTUEfSVRm+XL49NOkoxARyQ8FnSRGj4bf/GZ7vv66fNigQXDSScnFJCKSTwo6SSxeDF991eIXjfvtvz+MH88vEoeIiFSuoJNE2SWwjVLOvBx2WGgR9qmnEglJRCSvFHSSSL0EtsxWW4XG/v75z2RiEhHJJwWdJCo7kjALRxNjx8K33yYTl4hIvijoJNGzJ+y88yIaN/7l8N//HiZPhnbtEglLRCRvFPR9EoceCu3bT6Vp04G/GN6tW+hERCSznB1JmFk3MxtrZh+b2UdmdmYc3t7MXjazGfE1kf/z06fDb38L8+YlsXYRkfyQy+qmVcC57r4FsCNwupn1AS4AXnX3zYBXY39O/PWvMGzYTpWOW7MGHnwQHnssV2sXEcl/OUsS7j7P3d+L75cDHwNdgIOAsgYy7gcOzlUMS5fC0qWNKx23xRawzTYwalSu1i4ikv/q5JyEmXUHtgXeBkrcfR6ERGJmG6SZZzgwHKCkpIRx48ZVe70zZ/akQYMuaefdYYdu3HXXJowa9RadOv1Y7eXng9LS0hp9doVAZR+XdBiJKeby13rZ3T2nHdAKmAQcGvuXVBi/uKpl9O/f32vi7LPdmzdfmXb855+7g/uVV9Zo8Xlh7NixSYeQGJW9eBVz+VPLDkz0dfwNz+klsGbWGHgceMjdyx4e+rWZdYrjOwELcrX+1auhYUNPO7579/AwohYtchWBiEh+y1l1k5kZcC/wsbtfnzJqDHAccHV8fTpXMfTvDzNnLgQ6p51GJ65FRNLL5TmJXYDfAlPMbHIcdhEhOTxqZicCXwKH5yqAY4+FjTb6L5mSBIQjjtmzw5GFiIiUy1mScPd/A5Zm9J65Wm9NHH00TJoEM2aEZjtERCQo6GY5sjV0KHz2Gbz9dtKRiIjUL0oShJPXzZrBAw8kHYmISP2iJAG0aQOHHBKeZLdiRdLRiIjUH0oS0XHHhSfZPfdc0pGIiNQfBd0KbHXstRe89BIMHpx0JCIi9YeSRNSwIQwZknQUIiL1i6qbUqxaBeedByNHJh2JiEj9oCSRolEjeO01+Pvfk45ERKR+UJKo4Ljj4P33YcqUpCMREUmekkQFRx0Vjih0z4SIiJLEWjp2hP32C0+tW7Uq6WhERJKlq5sqMXw4tGwJS5aEpCEiUqyUJCqx776hExEpdqpuymDKFFiQs0ciiYjUf0oSacyZA337wp13Jh2JiEhylCTS6NoV9twT7r0X1qxJOhoRkWQoSWRw8snwxRfw8stJRyIikgwliQwOPhg6dIB77kk6EhGRZChJZNC0abgD+5VX4Pvvk45GRKTuKUlU4cILQ5VTixZJRyIiUvd0n0QVym6mcw+vZsnFIiJS13QkkYXPPoN+/cJDiUREiomSRBa6dYOvv4Zbbkk6EhGRuqUkkYUmTcLlsM8+C59/nnQ0IiJ1R0kiS6ecAg0awO23Jx2JiEjdUZLIUteu4b6Je++FH35IOhoRkbqhq5uq4YIL4MADoWHDpCMREakbShLVMGBA6EREioWqm6qptBSuvRbeey/pSEREck9Joprc4bLL4Prrk45ERCT3lCSqqXVr+N3v4JFHwjMnREQKmZJEDZx5ZnjGhG6uE5FCpyRRAz16wKGHhqfWlZYmHY2ISO7o6qYaOvdcWLIEFi6EVq2SjkZEJDeUJGpoxx31xDoRKXw5q24ys/vMbIGZTU0ZNsLMvjKzybHbN1frrytz5sCHHyYdhYhIbuTynMRIYO9Kht/g7v1i91wO159z7rDnnnDaaUlHIiKSGzlLEu4+Hvg2V8uvD8xCgpgwAd5+O+loRERqn3nZI9dysXCz7sAz7r5V7B8BHA8sAyYC57r74jTzDgeGA5SUlPQfPXp0jWIoLS2lVQ7PLH//fUN+/eud2G67xVx22Uc5W09N5br89ZnKXpxlh+Iuf2rZBw0aNMnd160xIXfPWQd0B6am9JcADQlHMFcA92WznP79+3tNjR07tsbzZuvii93N3KdNy/mqqq0uyl9fqezFq5jLn1p2YKKv4+94nd4n4e5fu/tqd18D3A1sX5frz5Wzzgp3Yr/xRtKRiIjUrjq9BNbMOrn7vNh7CDA10/T5omNH+PJLWG+9pCMREaldWSUJM2sHdAZ+AGbFI4Gq5hkFDAQ6mtkc4C/AQDPrBzgwCzilRlHXQ2UJYt486NQp2VhERGpL2iRhZusBpwNHAU2AhUAzoMTM3gJuc/ex6eZ396MqGXzvuoVbv/3973D++eE52BtumHQ0IiLrLtM5iX8Cs4Hd3L2Xu+/q7gPcvRtwNXCQmZ1YJ1HmiX32gZ9+ghtuSDoSEZHakfZIwt2HZBg3CZiUk4jy2Gabwa9/DbffHh512q5d0hGJiKybrK5uMrMuZrazme1e1uU6sHx14YWwfHmoehIRyXdVnrg2s78CRwDTgNVxsAPjcxhX3urbFw45BO67Dy6+GBqpCUURyWPZ/IQdDPRy9xU5jqVg3HQTtGypBCEi+S+bn7GZQGNASSJL3bqFV/dwIrtp02TjERGpqWySxPfAZDN7lZRE4e5n5CyqArBiBey6K+y1F1x1VdLRiIjUTDZJYkzspBqaNg1XO918M5x9NmywQdIRiYhUX5VJwt3vr4tACtEll8Ajj8A114RORCTfpL0E1sweja9TzOzDil3dhZi/eveGY46BW2+F+fOTjkZEpPoy3SdxZnzdHzigkk6ycMkl4eT1rbcmHYmISPVlqm6aD+DuX6SbwMwstlkuaWy6KbzyCuy8c9KRiIhUX6YjibFm9kcz2yh1oJk1MbPBZnY/cFxuwysMAwdCkybhiEJEJJ9kShJ7E+6wHmVmc81smpnNBGYQWoa9wd1H1kGMBWHSJOjZE959N+lIRESyl6mBvx+B24DbzKwx0BH4wd2X1FFsBWXzzWHlSjjvPHjtNTBLOiIRkapl1cCfu69093lKEDXXunU4iT1uHDz/fNLRiIhkp06fcV3shg8PJ7LPPx9Wr656ehGRpClJ1KHGjeGKK2DqVHjuuaSjERGpWrbPuN4Y2MzdXzGz5kAjd1+e29AK0+GHQ8eOMGhQ0pGIiFStyiMJMzuZ8CjTO+OgrsBTOYypoJnB4MHhtbQ06WhERDLLprrpdGAXYBmAu88A1FzdOhozBrp2hRkzko5ERCS9bJLECnf/+TYwM2tEeDKdrIPtt4c1a+Ccc5KOREQkvWySxOtmdhHQ3MyGAI8B/8ptWIVvww3hz3+GZ56BF15IOhoRkcplkyQuABYCU4BTgOeA/8tlUMXijDPCJbFnnx1utBMRqW+qTBLuvsbd73b3w919WHyv6qZa0LQpXH89TJ8Or7+edDQiImur8hJYM5vC2ucglgITgcvd/ZtcBFYs9t8fPv4YevVKOhIRkbVlc5/E84SG/h6O/UfG12XASPRsiXViVp4gpk2DPn2SjUdEJFU2SWIXd98lpX+KmU1w913M7De5CqzYvPAC7LMPPPUUHHRQ0tGIiATZnLhuZWY7lPWY2fZAq9i7KidRFaE994SttoI//lE32YlI/ZFNkjgJuMfMPjezWcA9wMlm1hK4KpfBFZPGjeGOO2D2bLj00qSjEREJsrm66V133xroB/Rz977u/o67f+fuj+Y8wiKyyy5w0klwww3wwQdJRyMikn0Df/sBWwLNLD4tx90vy2FcRevqq+HVV8Nlsdtsk3Q0IlLssrkE9g6gBTCIUNU0DHgnx3EVrQ4dQoJo3DjpSEREsjsnsbO7HwssdvdLgZ2AbrkNq7g1bgzu8NBD8NFHSUcjIsUsmyTxY3z93sw6AyuBHrkLSQCWLoWzzoITToBVuoZMRBKSTZL4l5m1Ba4B3gNmAaOqmsnM7jOzBWY2NWVYezN72cxmxNd2NYy74LVtCzffDO++G5ruEBFJQsYkYWYNgFfdfYm7Pw5sDPR290uyWPZIYO8Kwy6Iy9sMeDX2SxpHHAEHHwyXXAKffJJ0NCJSjDImCXdfA1yX0r/C3Zdms2B3Hw98W2HwQcD98f39wMFZR1qEzOD226FlS/jNb2D16qQjEpFiY1U16GpmlwIfAk9Ut/VXM+sOPOPuW8X+Je7eNmX8YnevtMrJzIYDwwFKSkr6jx49ujqr/llpaSmtWrWqesJ67M03O7B6tbHbbouqPW8hlL+mVPbiLDsUd/lTyz5o0KBJ7j5gnRbo7hk7YDmwhnDCelnsX1bVfHHe7sDUlP4lFcYvzmY5/fv395oaO3Zsjeetj1asqN70hVb+6lDZi1cxlz+17MBEz+I3NlOXzR3Xrd29gbs3dvc2sb9NDXPS12bWCSC+LqjhcorS//t/oX2npVlV+ImIrLtsrm7CzA40s2tjt/86rG8McFx8fxzw9Dosq+j06QMzZ8LppycdiYgUiyqThJldDZwJTIvdmXFYVfONAv4D9DKzOWZ2InA1MMTMZgBDYr9kaYcdwpVODz0UOhGRXMum7aZ9CQ37rQEws/uB96ni8lV3PyrNqD2rFaH8wkUXwcsvwymnQP/+0Lt30hGJSCHLqroJaJvyfr0cxCFZatQIRo2CVq3gjTeSjkZECl02RxJXAe+b2VjAgN2Bi3IalWTUtWtoBHA9pWsRybEqk4S7jzKzccD/EJLE+e4+P9eBSWZlCeKVV2DRIjjyyMzTi4jURDYnrl9193nuPsbdn3b3+Wb2al0EJ5m5h+dP/O53MHVq1dOLiFRX2iRhZs3MrD3Q0czaxcb52se7qDvXWYSSlhn84x/Qpg0ceigsWZJ0RCJSaDIdSZwCTAJ6x9ey7mng1tyHJtno1Akeeww+/xyOPlrtO4lI7UqbJNz9JnfvAfzJ3Xu6e4/YbePut9RhjFKF3XYLzYo//zw88kjS0YhIIUl74trM/geY7e43x/5jgcOAL4AR7l6xhVdJ0KmnQvfuMHRo0pGISCHJVN10J/ATgJntTrg7+gFgKXBX7kOT6tp773Ce4tNP4cMPk45GRApBpktgG6YcLRwB3OXhwUOPm9nknEcmNeIOhx0GixfDW28lHY2I5LtMRxINzawsiewJvJYyLpub8CQBZjByJHz7LRxwAPzwQ8OkQxKRPJYpSYwCXjezp4EfgDcAzGxTQpWT1FPbbguPPgqTJ8Nll/Vh1aqkIxKRfJXp6qYrgHMJz6reNT7AomyeP+Y+NFkX++4Lt94Kb73VgZtuSjoaEclXGauN3H2tWm13/2/uwpHadOqpMHv2NE47rU/SoYhInsq2FVjJU0OGLKB583A39gsvJB2NiOQbJYkicfHFsP/+8MwzSUciIvlESaJIXH01bLcdHH44vP560tGISL5QkigSrVvDc89Bjx7h0tj33ks6IhHJB0oSRaRjR3jpJWjfHoYPDzfeiYhkopviikzXruFBRU2bhhvvREQy0ZFEEdp0U+jWDdasgfPOg08+SToiEamvlCSK2Ny58MADMHCgEoWIVE5Jooh17Qpjx4b3ShQiUhkliSK3xRbliWKPPUJ7TyIiZZQkhC22CPdOrLdeaD1WRKSMrm4SAHr1go8+gsaNQ//s2eHktogUNx1JyM/KEsTTT4croEaPTjYeEUmekoSsZdAg2HFHOPpouPFG3XQnUsyUJGQtbdqEFmMPOQTOPhvOOAM9uEikSClJSKWaN4fHHoNzz4VbboEXX0w6IhFJgpKEpNWgAVx7Lbz5Juy3Xxi2enWyMYlI3VKSkCrttFN4nTQJ+vSBd95JNh4RqTtKEpK1xo3hp59gt93gvvuSjkZE6oKShGStb1+YODEkiRNPhNNPD0lDRAqXkoRUS4cO4cqnP/0JbrsN7rgj6YhEJJcSuePazGYBy4HVwCp3H5BEHFIzjRrBNdfAkCEweHAYVloKrVolG5eI1L4kjyQGuXs/JYj89atfhYSxaFE4of3nP+t+CpFCo+omWWfNm4ejissvD0cWc+YkHZGI1BbzBNpcMLPPgcWAA3e6+12VTDMcGA5QUlLSf3QNGxIqLS2lVRHXg9Rl+V9+uYTrr9+cJk3WcM4509ljj0V1st50innbF3PZobjLn1r2QYMGTVrn2hp3r/MO6BxfNwA+AHbPNH3//v29psaOHVvjeQtBXZd/+nT3/v3d99nHfc2aOl31Wop52xdz2d2Lu/ypZQcm+jr+Xidy4trd58bXBWb2JLA9MD6JWKR2bb45/Oc/4US2GcyaBf/9bzh/ISL5p87PSZhZSzNrXfYe+BUwta7jkNxp3BjatQvvL78chg6F006D5cuTjUtEqi+JE9clwL/N7APgHeBZd38hgTikDtx8M5xzTrifok8fGDMm6YhEpDrqPEm4+0x33yZ2W7r7FXUdg9Sd5s3huutCI4Ft28JBB8E99yQdlYhkS5fASp3YcUd47z24/nr49a/DsLlz1aqsSH2nJCF1pnHj8BCjNm1g5cpwMnvAABivSxZE6i0lCUlEo0bhDu1vvoE99oAjjoAvv0w6KhGpSElCEmEWEsMnn8Bf/hJOaPfqBe++m3RkIpJKSUIS1aIFjBgB06fDH/4A224bhk+ZEqqkRCRZShJSL2y0UWhZtlEj+O472HNP2GILGD0a1qxJOjqR4qUkIfVOixYwciS0bAlHHQX9+8Ozz0ICzYyJFD0lCal3zGDffeH99+Ghh2DZMth/fz1bWyQJShJSbzVoAEcfHU5uP/007LBDGH7ddaEaSvdYiOSekoTUe40bw4EHhverV8OoUaEaqk8fuPNO+P77ZOMTKWRKEpJXGjYM1U7//Ce0bg2nnhpOej/9dNKRiRQmJQnJOw0awGGHhXsqxo+H3XYL91gAfPxxOJchIrVDSULylllIEE8+Cb17h2GXXw7bbRfainr++Q1VFSWyjpQkpKDccgvccAMsXQp/+1tvOneGSy9NOiqR/KUkIQWlXTs46yyYNg1uvPF99tuvfNyqVXDvvbB4cWLhieQdJQkpSGawzTZLeeih0DYUwBtvwEknwYYbwrBh4WT3Tz8lG6dIfackIUVj4ECYOBF+//uQMA4+GDp3Ds/gFpHKKUlI0TALTXzceCN89VVo6mPYMNhkkzD+8svDs7hfey1UTYmIkoQUqUaNQtMfd9wR7r0AWLgQ7r8/NC7YqROcfDKMHZtsnCJJU5IQiW66KSSKxx8PT8175JGQNCA0LjhyJMyenWiIInWuUdIBiNQnLVrAoYeG7scfw6W0ENqPOuGE8H7rrUODg3vvHe7HaNIkuXhFck1HEiJpNGsGJSXhfe/e8NFH8Le/Qfv24XWPPcJ5DQjnOCZNUqODUniUJESyYBYaFPzf/4Vx42DRInjqKRg8OIz/xz9gwABYf/3QZMitt8IHH+iBSZL/VN0kUgNt28JBB5X3n3ACdOsGr74auieeCNVQS5eGI5JXXgmJZocdoFWrxMIWqTYlCZFaUFICxxwTOneYNSucx2jWLIwfMQImTAhXUvXrF85lDBwYLsEVqc+UJERqmRn06BG6Ms88A2+9FRLFhAnhqqkvvihPEkcdFRLNgAGh23zz0NqtSNKUJETqQNu24WqovfcO/WvWwJIl4f3KlTBnDowZEy7DhXCV1cUXw0UXhZPhr78erqpaf/0kopdipiQhkoAGDcJVUhCevPfGG+Eu708+CU2HTJ4MW24Zxs+cGW7wg9Du1NZbh+7YY2GbbUL1FoQjGJHapiQhUk80agRbbRW6VF26wEsvwZQpofvwQ7jtNthll5Akxo4NV1T16lXerVzZkb59yxORSE0pSYjUcy1awJAhoSuzenX55bUdO8KRR8L06eEqqgceANiKffaB7bcPd5DfeGP5eZIePaBnzzCu7MS6SDpKEiJ5qGHD8jan+vaF228vH7d8OTz88ES23noAEKq2GjQI93c8+GB59dS8eaH66vbbw9P9unb9ZTd0aPk6pHgpSYgUmNatoVevUpo3D/2HHBI6CM/P+PLLcJ6j7G5ys3A/x7RpIXGsWRPOk/z4Yxh/xhnwwgvlyaNLF9hoo9DkOsDcueGekPbtdUVWIVKSECkiTZrAppuGrsypp4YOwsnz+fPh66/Lf/D79oUFC8IVWOPHhyZIunYtTxInnBDOmTRqFK6+KikJTbLfc08Y/9hj8P330KFDSCQdOsAGG4SnCEr9pyQhIj9r1Kj8iKHMSSeFrox7qNIqc845odn1r78u71KvtLr00tDuVarBg8Od6RDOtSxbFpJHWSLZccdw7wjA889D06aw3nq/7NSwYt1QkhCRajGDNm3K+4cODV06EybAN9/8sks9iujZM9yhvmBBuAT4m2/g22/Lk8SwYeFIJNVJJ8Hdd4eEtcUWoYotNYH06NGRgQND9dptt4WmUFq1gpYtw+umm4ZmVFavDs88b9UqJCJdRry2RJKEme0N3AQ0BO5x96uTiENEcq/sh7tnz8rH33nn2sPKTq4D/Pvf4ZxJWbdsWWiVF8KNiP36lQ+fPz+8Hzq0BRBuWDz77LWXf+WVcOGF4fxMWVwNGpQnkiuvhOOPh88+gxNPDFeBNW9e/nrSSbDzzmH+++//5bhmzUILwV27hgQ0fXpIQGVdkyahuq1p0/JWg+vzBQJ1niTMrCFwKzAEmAO8a2Zj3H1aXcciIvVT6j/6bbdNP12TJjB69NrDx437EuhJx47hyOS770JXWhpeN944TNe2Lfz97+XDy6YpG796dUhYixeHE/Q//gg//BCeJwLhAoBLLll7/U8/HZLEhAlwwAFrj3/tNRg0KJyvOeqokKDKEkjTpvDiiyH51QdJHElsD3zq7jMBzGw0cBCgJCEitarszvZ0NxW2awd//GP6+TffPDSJks7AgeFo5ocfyhPIjz+Gx99CaPX3+efDsBUrQvXXihXlR0JbbQWXXfbLcT/9VL9ugjRPPa6rixWaDQP2dveTYv9vgR3c/Q8VphsODAcoKSnpP7qyvwtZKC0tpVURt81czOVX2Yuz7FDc5U8t+6BBgya5+4B1WV4SRxKVnRpaK1O5+13AXQADBgzwgQMH1mhl48aNo6bzFoJiLr/KPjDpMBJTzOWv7bIncevLHKBbSn9XYG4CcYiISBWSSBLvApuZWQ8zawIcCYxJIA4REalCnVc3ufsqM/sD8CLhEtj73P2jKmYTEZEEJHKfhLs/BzyXxLpFRCR7ao5LRETSUpIQEZG0lCRERCStOr+ZribMbCHwRQ1n7wgsqsVw8k0xl19lL17FXP7Usm/s7uuvy8LyIkmsCzObuK53HOazYi6/yl6cZYfiLn9tl13VTSIikpaShIiIpFUMSeKupANIWDGXX2UvXsVc/lote8GfkxARkZorhiMJERGpISUJERFJq6CThJntbWbTzexTM7sg6XhywcxmmdkUM5tsZhPjsPZm9rKZzYiv7VKmvzB+HtPNLMPj6+sfM7vPzBaY2dSUYdUuq5n1j5/Zp2b2dzOr7Bkn9U6a8o8ws6/i9p9sZvumjCuY8ptZNzMba2Yfm9lHZnZmHF7w2z9D2etm27t7QXaEFmY/A3oCTYAPgD5Jx5WDcs4COlYY9jfggvj+AuCv8X2f+Dk0BXrEz6dh0mWoRll3B7YDpq5LWYF3gJ0ID8B6Htgn6bKtQ/lHAH+qZNqCKj/QCdguvm8N/DeWseC3f4ay18m2L+QjiZ+fpe3uPwFlz9IuBgcB98f39wMHpwwf7e4r3P1z4FPC55QX3H088G2FwdUqq5l1Atq4+388fGseSJmnXktT/nQKqvzuPs/d34vvlwMfA10ogu2foezp1GrZCzlJdAFmp/TPIfMHm68ceMnMJsXnggOUuPs8CDsYsEEcXoifSXXL2iW+rzg8n/3BzD6M1VFl1S0FW34z6w5sC7xNkW3/CmWHOtj2hZwksnqWdgHYxd23A/YBTjez3TNMWyyfCaQva6F9BrcDmwD9gHnAdXF4QZbfzFoBjwNnufuyTJNWMiyvy19J2etk2xdykiiKZ2m7+9z4ugB4klB99HU8tCS+LoiTF+JnUt2yzonvKw7PS+7+tbuvdvc1wN2UVx8WXPnNrDHhR/Ihd38iDi6K7V9Z2etq2xdykij4Z2mbWUsza132HvgVMJVQzuPiZMcBT8f3Y4AjzaypmfUANiOcyMpn1SprrJJYbmY7xis7jk2ZJ++U/UBGhxC2PxRY+WOs9wIfu/v1KaMKfvunK3udbfukz9zn+KqAfQlXAnwGXJx0PDkoX0/CVQwfAB+VlRHoALwKzIiv7VPmuTh+HtOp51d1VFLeUYTD6pWEf0Un1qSswID4hfoMuIXY8kB979KU/x/AFODD+OPQqRDLD+xKqBr5EJgcu32LYftnKHudbHs1yyEiImkVcnWTiIisIyUJERFJS0lCRETSUpIQEZG0lCRERCQtJQkpCmbW1sxOS+nvbGb/zNG6DjazS+L7kWY2LIt5mpjZeDNrlIuYRGpKSUKKRVvg5yTh7nPdvcof7xo6D7gt24nNrKGHRihfBY7IUUwiNaIkIcXiamCT2O7+NWbW3eJzGczseDN7ysz+ZWafm9kfzOwcM3vfzN4ys/Zxuk3M7IXYmOIbZta74krMbHNghbsvShm8u5m9aWYzy44qzGxgfEbAw4QbogCeAo7J4WcgUm06tJVicQGwlbv3g59b00y1FaF1zWaEppXPd/dtzewGQvMFNxIeMH+qu88wsx0IRwuDKyxnF+C9CsM6Ee6a7U24M7asmmv7GNPnsX8q8D81L6JI7VOSEAnGemirf7mZLQX+FYdPAfrGFjh3Bh5LeZhX00qW0wlYWGHYUx4aYZtmZiUpw99JSRC4+2oz+8nMWsdYRBKnJCESrEh5vyalfw3he9IAWFJ2JJLBD8B6GZad2lzzd5XM3xT4sapgReqKzklIsVhOePRjjXhov/9zMzscQsucZrZNJZN+DGxak3WYWQdgobuvrGmcIrVNSUKKgrt/A0wws6lmdk0NF3MMcKKZlbW6W9njcMcD22b1gPm1DQKeq2FsIjmhVmBFapmZ3QT8y91fqeZ8TwAXuvv03EQmUn06khCpfVcCLaozQ3ww1lNKEFLf6EhCRETS0pGEiIikpSQhIiJpKUmIiEhaShIiIpKWkoSIiKT1/wFjm5JtHGtdWwAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.plot(t,S,'b--')\n",
"plt.xlabel('time (hr)')\n",
"plt.ylabel('Storage (in)')\n",
"plt.grid()\n",
"plt.title('Response of a Watershed to 3-hr Rainfall Pulse')\n",
"plt.show()"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "cc263014",
"metadata": {
"id": "cc263014",
"outputId": "89f8f853-38bb-42e0-df58-6ad940dd98e4"
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAEGCAYAAACQO2mwAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAyH0lEQVR4nO3dd3xUVdrA8d+TAgQINYrRgIAIRBADicCuqAkYRWRFyirKrqyiwYINy6Ls6sorioW2ioVVXoFFWDuoiIUimHdRQUAQRIqhS5GWUENy3j/OpJDMJJMyc6c838/nfmbm3Dtzn5MJebjnniLGGJRSSil3IpwOQCmlVODSJKGUUsojTRJKKaU80iShlFLKI00SSimlPIpyOgBvxMXFmebNm1fqvUeOHKFOnTrVG1AQCef6a93Ds+4Q3vUvXvfly5fvM8acUZXPC4ok0bx5c5YtW1ap9y5atIjU1NTqDSiIhHP9te6pTofhmHCuf/G6i8iWqn6eNjcppZTySJOEUkopjzRJKKWU8kiThFJKKY80SSillPJIk4RSSimPNEkopZTyKKSTRF4efPxxPDk5TkeilFLBKaSTxPffw9ixbRg0yOlIlFIqOIV0kjh1yj7OmQO6tpJSSlVcSCeJmJii5xs2OBeHUkoFq5BOEklJ8MknS/j1V2jd2ulolFIq+IR0kgCoXTuPJk2cjkIppYJTSCeJzEzIyEhm5ky46ir4+WenI1JKqeAS0kni4EHYsCGWmjXh88/h/fedjkgppYJLSCeJgt5NLVpA587w3nvOxqOUUsEmLJJEZCQMGADLlkFWlqMhKaVUUAnpJJGXZx+joqB/f/tcryaUUsp7IZ0k4uIgKekAdetCy5Zw880QH+90VEopFTyCYo3ryureHcaPX0WzZqkATJ3qbDxKKRVsQvpKwp3sbNi40ekolFIqOIT0lcSsWfDQQ51ZvpzCAXVpaVC3Lixa5GhoShWKi4ujefPmToehHJSVlcW+ffucDsOtkE4SBw7Ajh21T5vcr3dvGDUKdu9GR2KrgNC8eXOWLVvmdBh+N2P1DEbOH8nWQ1tpVr8Zo3uMBihVNujC0J/GOSUlxekQPArpJFHQBTaqWC3794cnn4QPP4ShQx0JS6mwN2P1DDI+yuBo7lEAthzawi0f3oKIcDLvZGFZxkcZZG7NZO6GueUmk6qUhUMiqqyQThLFu8AWaN/eTvb37ruaJJRyysj5IwsTRIHc/NxSxx3NPcory14pfO0pmVSlrCqJqNf5vUq9L9QSTkgnCXdXEiL2auK552D/fmjUyJnYlApnWw9trfR73SWTqpRVJRGVfF9FEk6wJJOQThItW8Lvf7+P6Oi408rvvBNuugkaNnQoMKXCXLP6zdhyaIvTYXjkbYIpyW3C+eAWJKL01QsQFIkipLvA9usHo0evoWbN08ubNrXNTiLOxKVUuBvYfmCpsuiIaGpE1jitTAj+f6S5JrcwQRQ4mnuUkfNHOhRRxfgsSYhIUxFZKCLrRORHEbnPVd5IRL4QkQ2uR0f+P79+Pfz5z7BrlxNnVyp8bTu0jde/f52E2ASa1muKIJxb/1z+97r/ZUqfKZxb/9zCsjtS7qB2dO3T3u8umVSlzKlEVJUmN3/yZXPTKeBBY8z3IhILLBeRL4C/APONMWNEZAQwAvirLwJ49ll4/vnf4a77cX4+/PvfcPHFcO+9vji7UqpA8e6u0ZHRCELmrZm0iWtT6tiSTTCXNLukWnsylSzrdX4vpq6aetqN9OiI6NPuP3gqK0kQDMbj/uKa1W/m1XFO81mSMMbsAna5nmeLyDrgHKAPkOo6bCqwCB8liUOH4NChaLf7EhPhootg5kxNEkr5UsnurifzTlIjsgbLdi1zmyRKGnThILdt99VZVtlEVLJ3k7cJp3Z07cLPC3RijHdZr0onEWkOLAbaA1uNMQ2K7TtgjCnV5CQiGUAGQJMmTZJnzZpV4fO++mpLPvjgHD77bInb/TNnNmXy5PN4662lxMcfr/DnB4OcnBzq1q3rdBiOCJa6P/TQQyE9mK75hOZub1KfW/9csu7P8n9APlaZQYJt2rThtddeq5bzF/+9T0tLW26MqdJIPZ8nCRGpC3wFjDbGvC8iB71JEsWlpKSYyvwjGj4cXn31FEePur9gysqyCxI9/TQ8+miFPz4oLFq0iNTUVKfDcESw1D0lJSWkk0TEkxFum2AEIf+JfAciCjzV+TtQ/PdeRKqcJHzau0lEooH3gBnGmILFQ3eLSLxrfzywx1fnz8uDyEjPSbB5c7sYUe3aHg9RSlVRQr0Et+XB0iYf7nx2T0JEBHgDWGeMGVds1xxgMDDG9TjbVzEkJ8PmzXuBsz0e8847vjq7UgqgRYMWbDu87bSyYGqTD3e+7N10CfBnYLWIrHSVPYZNDm+LyBBgK/BHXwVw883QrNnPlJUkwF5xbNtmryyUUlVX0C5fcC8i+axk9h3bF5QjjsOdL3s3fQ0eOyD38NV5K+Omm2D5ctiwQQfYKVVVJXszAaz7bR2T/zBZE0MQCukR19666irYtAm++cbpSJQKfu4m7wumEcbqdJoksDeva9WCadOcjkQFi23btpGWlkZiYiLt2rVj4sSJADz88MO0bduWDh060LdvXw4ePOhsoA7wNJI4WEYYq9NpkgDq1YO+fe1KdidOOB2NCgZRUVGMHTuWdevWsXTpUiZNmsTatWtJT09nzZo1/PDDD7Ru3ZpnnnnG6VD9rkGtBm7LtTdTcNIk4TJ4sF3Jbu5cpyNRwSA+Pp5OnToBEBsbS2JiIjt27ODKK68kyjU3fdeuXdm+fbuTYfrdN9u/4dDxQ0RK5Gnl2pspeGmScLniCvj8c7j2WqcjUcEmKyuLFStW0KVLl9PKp0yZwtVXX+32PZMnTyYlJYWUlBT27t3rjzB97rejv3H9u9fTrEEzXr7m5dMm6tOb1sErpNeTqIjISEhPdzoKFWxycnLo378/EyZMoF69eoXlo0ePJioqikGD3P9hzMjIICPDrikQyOsbe2PG6hk8Nv+xwnsOo9JGkZGcQUZyhsORqeqgVxLFnDoFjzwCb77pdCQqGOTm5tK/f38GDRpEv379CsunTp3Kxx9/zIwZM5AQ71Nd0N21+E3pMV+PYcbqGQ5GpaqTJolioqJgwQL45z+djkQFOmMMQ4YMITExkeHDhxeWz5s3j2effZY5c+ZQOwzme9HurqFPk0QJgwfDihWwerXTkahAlpmZyfTp01mwYAFJSUkkJSUxd+5chg0bRnZ2Nunp6SQlJXHHHXc4HapPaXfX0Kf3JEq48UY7e+y0afD8805HowJVt27dcDeDcq9evRyIxhl5+XnUiKzBibzS/ca1u2vo0CuJEuLi4Jpr7Kp1p045HY1SgWvUV6M4kXei1HKg2t01tGiScCMjA7p3hzAcLKtUmWasnkHzCc2JeDKCUYtHcVmzy0qtS63dXUOLNje50auX3ZRSRdxN3Pfdzu/IICMkV5hTll5JlGH1atjjsyWRlAou7noyHTt1THsyhThNEh5s3w4dOkA1LTurVNDTnkzhSZOEBwkJ0KMHvPEG5OsyvEoRVzvObbn2ZAptmiTKcPvtsGULfPGF05Eo5awtB7dw5OQRpMQ6YtqTKfTpjesyXHcdNG4Mr79uFyZSoWfZsmUsWbKEnTt3EhMTQ/v27bniiito1KiR06EFjJN5J7n+3euJioxibPexTPxmoi5DGkb0SqIMNWvaEdhffglHj5Z/vAoeb775Jp06deKZZ57h2LFjtGnThjPPPJOvv/6a9PR0Bg8ezNat4d3WXtDdteZTNfl2x7f85aK/8MDvHiDr/izyn8gn6/4sTRBhQK8kyvHoo/DEExAG0/CElSNHjpCZmUlMTIzb/StXrmTDhg00axae7e3uuru+vuJ1Oid01sQQZjRJlCPOda+uYAaGEJ/UM2zcfffdZe5PSkryTyABqqyJ+zRJhBdNEl7YtAn69YPnntN7E6Fm7969/Otf/yIrK4tTxeZhmTJlioNROU+7u6oCmiS80LQp7N4NL72kSSLU9OnTh0svvZQrrriCyMjI8t8QJurUqEPOyZxS5drdNfxokvBCjRq2O+zo0fDLL9CihdMRqepy9OhRnn32WafDCChv//g2OSdziIqI4lR+0dWVdncNT9q7yUtDh0JEBLzyitORqOrUu3dv5s6d63QYjis+cd/AdwdyXsPzeOPaN3TiPqVXEt5KSLDjJt54A558Ejx0ilFBZuLEiTz99NPUrFmT6OhojDGICIcPH3Y6NL9x15NpZ/ZOIiMideI+pVcSFTFiBIwfD9p0HTqys7PJz8/n2LFjHD58mOzs7LBKEKAT96myaZKogJQUuPlme49CBbeffvoJgO+//97tVp5t27aRlpZGYmIi7dq1Y+LEiQC88847tGvXjoiICJYtW+bTOlQX7cmkyqLNTRWUkwOvvmoXJerUyeloVGWNGzeOyZMn8+CDD5baJyIsWLCgzPdHRUUxduxYOnXqRHZ2NsnJyaSnp9O+fXvef/99hg4d6qvQq118bDw7s3eWKteeTAo0SVSYMTBqFKxcaZc4VcFp8uTJACxcuLBS74+Pjyc+Ph6A2NhYEhMT2bFjB+np6dUWoz8cP3WcSCndfqo9mVQBbW6qoNhYuPVW+M9/7JoTKjh9/fXXZe4/fPgwa9as8eqzsrKyWLFiBV26dPH6/JMnTyYlJYWUlBT27t3r9fuq2wPzHmDb4W08+LsHtSeTckuvJCrhvvvgxRft4LoxY5yORlXGe++9xyOPPELPnj1JTk7mjDPO4Pjx42zcuJGFCxeyZcsWxo4dW+7n5OTk0L9/fyZMmEC9evW8Pn9GRgYZGRkApKSkVLoelTFj9QxGzh/JlkNbAOh9fm9euPIFXrjyBb/GoYKDJolKaNHCTtPx2mvwt79B3bpOR6Qqavz48Rw4cIB3332Xd955h127dhETE0NiYiJDhw6lW7du5X5Gbm4u/fv3Z9CgQfTr188PUVedu+6uC7IWMGP1DL1yUG5pkqikBx+Egwdh715NEsGqYcOG3H777dx+++0Vfq8xhiFDhpCYmMjw4cN9EJ1v6MR9qqI0SVRS1666Yl04y8zMZPr06Vx44YWFM8Y+/fTTnDhxgnvuuYe9e/dyzTXXkJSUxGeffeZssMVod1dVUT5LEiIyBegN7DHGtHeV/QO4HSi4U/eYMSao50TYvh3274cOHZyORPlTt27dMAXzx5fQt29fP0fjvca1G7Pv6L5S5drdVXniy95NbwI93ZSPN8YkubagThDGQI8ecNddTkeiVPk27d9E9olsIuT0f/ba3VWVxWdJwhizGNjvq88PBCI2QWRmwjffOB2NqojvvvuOX3/9tfD1tGnT6NOnD/feey/794fOr23xifsSJyUSKZGMu3KcdndVXnPinsQwEbkZWAY8aIw54O4gEckAMgCaNGnCokWLKnWynJycSr/XG+efH0mdOr/j4YcPMGrUjz47T2X5uv6BrKy6Z2Rk8MILL1CvXj1WrVrF//zP/3DPPfewadMm+vbty5NPPunfYH2gZE+m3PxcIiSCuDpxOnFfgMnOzq62f6fV/m/eGOOzDWgOrCn2ugkQib2CGQ1M8eZzkpOTTWUtXLiw0u/11siRxogYs3atz09VYf6of6Aqq+4dOnQofH7XXXeZJ554ovD1RRdd5Lug3KjK73dZzh1/ruEflNrOHX+uT86nKq86fweK/94Dy0wV/477dcS1MWa3MSbPGJMP/Avo7M/z+8r999uR2EuWOB2J8lZeXl7hcqXz58+ne/fuhfuKL2MazLQnk6oOfm1uEpF4Y8wu18u+gHfzHgS4uDjYuhXq13c6EuWtG2+8kcsvv5y4uDhiYmK49NJLAdi4cSP1Q+SLbFq/qduEoD2ZVEV4lSREpCFwNnAMyHJdCZT3nplAKhAnItuBJ4BUEUkCDJAFBM9UmeUo+Luyaxe45n1TAWzkyJH06NGDXbt2ceWVVyIiAOTn5/Piiy86HF31aH9G+1JJQnsyqYrymCREpD5wN3AjUAM7tqEW0ERElgIvG2M8TqFpjLnRTfEbVQs3sP3zn/DXv9p1sM86y+loVHm6du1aqqx169YORFL9Pv75Y+ZunEt6i3R+3v8zWw9tpVn9ZozuMVp7MqkKKetK4l1gGnCpMeZg8R0ikgz8WURaGmNC+g9/RVx9NTzwgF297tlnnY5GhZuCifu2HtqKiNCsXjPm3DSHWlG1nA5NBTGPN66NMenGmOklE4Rr33JjzP2aIE53/vlw/fXwyitwwG3HXqV8o6C765ZDWzAY8k0+e47u4b117zkdmgpyXvVuEpFzROT3InJZwebrwILVo49CdrZtelLKX9xN3Hf81HFdp1pVWbk3rkXkWeAGYC2Q5yo2wGIfxhW0OnSAvn1hyhQYORKidApF5Qfa3VX5ijd/wq4D2hhjTvg4lpAxcSLUqaMJQvnPObHnsD279FKJ2t1VVZU3zU2bgWhfBxJKmjaFRo3sBIAnNLUqH8s3+dSvVXpsh3Z3VdXBm//rHgVWish8oPBPnjHmXp9FFQJOnIBu3eCKK+CZZ5yORoWa4j2Z6tesz8ETBxnScQhfbv5Su7uqauVNkpjj2lQF1Kxpezu9+KLtFnvmmU5HpEJFyYn7Dp44SKREktY8jdevfd3h6FSoKbe5yRgz1d3mj+CC3eOPw7Fj8PzzTkeiQom7nkx5Jo+RC7Qnk6p+HpOEiLztelwtIj+U3PwXYvBq2xYGDYJJk6DY0gUqBGzbto20tDQSExNp164dEydOBGD//v2kp6dz/vnnk56ezgEfDJjRnkzKn8q6krjP9dgb+IObTXnh8cfh5EmbKFToiIqKYuzYsaxbt46lS5cyadIk1q5dy5gxY+jRowcbNmygR48ejBkzptrP7anHkvZkUr5QVpL4FcAYs8XdBiAFs6Ipj1q1gi+/hL//3elIVHWKj4+nU6dOAMTGxpKYmMiOHTuYPXs2gwcPBmDw4MF8+OGH1X7u2zvdXqpMezIpXykrSSwUkXtE5LT/nohIDRHpLiJTgcG+DS80pKZCjRr2ikKFnqysLFasWEGXLl3YvXs38a5pgOPj49mzZ4/b90yePJmUlBRSUlLYu3ev1+fKPpHN1FVTaVCzAQn1EnQJUuVzZfVu6gncCswUkRbAQewssJHA58B4Y8xKXwcYKpYvhz594IMP4OKLnY5GVZecnBz69+/PhAkTqFevntfvy8jIICMjA4CUlJQyjy3e3bV2dG2O5B5h0eBFXN788irFrpQ3PCYJY8xx4GXgZRGJBuKAY+4m/FPla90acnPhkUdgwQLQhrrgl5ubS//+/Rk0aBD9+vUD7Hrsu3btIj4+nl27dnFmFfs+l+zueiT3CNER0W5HVyvlC15N8GeMyTXG7NIEUXmxsfYm9qJF8OmnTkejqsoYw5AhQ0hMTGT48OGF5ddeey1Tp9oe4lOnTqVPnz5VOo+77q65+bk6cZ/yG7+ucR3uMjLsjey//hXy8so/XgWuzMxMpk+fzoIFC0hKSiIpKYm5c+cyYsQIvvjiC84//3y++OILRowYUaXzaHdX5TSdgs6PoqNh9Gi44QaYOxf+oB2Jg1a3bt0wxrjdN3/+/Go7T7P6zdhyaIvbcqX8wdv1JM4VkStcz2NEJNa3YYWuP/4R5s+H3r2djkQFg/6J/UuVaXdX5U/lJgkRuR27lOlrrqIE4EMfxhTSRKB7d/uYk+N0NCqQ7Ti8g6mrptKsXjOa1Wum3V2VI7xpbrob6Ax8A2CM2SAiOl1dFc2ZAzffDN99ZycCVApO7+5aI7IGxhgyb82kTVwbp0NTYcqb5qYTxpjCYWAiEoVdmU5VQefOkJ8PxTrGqDBXcp3qE3knQGDZrmVOh6bCmDdJ4isReQyIEZF04B3gI9+GFfrOOstO1fHxxzBvntPRqEDgrrvrybyT2t1VOcqbJDEC2AusBoYCc4G/+TKocHHvvbZL7AMP2IF2Krxpd1cViLxZTyLfGPMvY8wfjTEDXM+1uaka1KwJ48bB+vXw1VdOR6OcprO7qkDkTe8md+tJLBGR8SLS2B9BhrLevWHdOrvMqQpvo3uMpnZ07dPKtLurcpo3vZs+BfKAt1yvB7oeDwNvomtLVIkItHF1XFm7Fi64wNl4lHMKurUW9G7SdapVIPAmSVxijLmk2OvVIpJpjLlERP7kq8DCzbx5cPXV8OGHdrZYFT6ysrJOmwk2jjjMXkPcGXGM//d4xjO+wp+ZnZ1NbKzvxrzu3buXM844w2efXxV79+6lVq1aPq1/ZXn6uWVlZfk/GC95kyTqikgXY8w3ACLSGajr2nfKZ5GFmR49oH17uOce+7xu3fLfo0LDvn37SpWlpKSwbFnlu74uWrSI1NTUKkRVtqrG50spKSm88MILPq1/ZQXyz80Tb3o33Qa8LiK/iEgW8Dpwu4jUAZ7xZXDhJDoaXn0Vtm2DJ590OhqllLK86d30nTHmQiAJSDLGdDDGfGuMOWKMedvnEYaRSy6B226D8eNh1Sqno1FKKS9ngRWRa4B2QK2CZa2NMaN8GFfYGjPGTgC4fj1cdJHT0SinFKxaF6gCOT6NrXp50wX2VeAG4B5AgD8C5/o4rrDVuLFNENdf73QkykmB/sckkOPT2KqXN/ckfm+MuRk4YIx5Evgd0NS3YYW36GgwBmbMgB9/dDoapVQ48yZJHHc9HhWRs4FcoIXvQlIAhw7B/ffDLbfAKe1DppRyiDdJ4iMRaQA8D3wPZAEzy3uTiEwRkT0isqZYWSMR+UJENrgeG1Yy7pDXoAG8+KKdSnzcOKejUf40b9482rRpQ6tWrRgzZozT4RTatm0baWlpJCYm0q5dOyZOnOh0SKXk5eXRsWNHHn30UadDKeXgwYMMGDCAtm3bkpiYyH//+1+nQ/JKmUlCRCKA+caYg8aY97D3ItoaYx734rPfBHqWKBvh+rzzgfmu18qDG26A666Dxx+Hn35yOhrlD3l5edx99918+umnrF27lpkzZ7J27VqnwwIgKiqKsWPHsm7dOpYuXcqkSZMCJrYCEydOJDEx0ekw3Lrvvvvo2bMnP/30E6tWrQrYOEsqM0kYY/KBscVenzDGHPLmg40xi4H9JYr7AFNdz6cC13kdaRgSgVdegTp14E9/grw8pyNSvvbtt9/SqlUrWrZsSY0aNRg4cCCzZ892OiwA4uPj6dSpEwCxsbEkJiayY8cOh6Mqsn37dj755BNuu+02p0Mp5fDhwyxevJghQ4YAUKNGDRo0aOBsUF7ypgvs5yLSH3i/GmZ/bWKM2QVgjNlV1gp3IpIBZAA0adKERYsWVeqEOTk5lX5voHjwwcbk5QlLlpQemVueUKh/ZQVj3b/66iuioqIK487OzmbdunUVroev6/7rr7+ydOlSMjIyAuZn/MQTT3DTTTfxww8/cOrUqYCJC2Djxo3UqlWLq6++mk2bNtG6dWuGDRtGTExMtZ+r2r97Y0yZG5AN5GNvWB92vT5c3vtc720OrCn2+mCJ/Qe8+Zzk5GRTWQsXLqz0ewPRiRMVOz7U6l8RwVj3t99+2wwZMqTw9bRp08ywYcMq/Dm+rHt2drbp1KmTee+993x2jor66KOPzJ133mmMsXXv2rWrwxGd7rvvvjORkZFm6dKlxhhj7r33XvO3v/3NJ+cq/t0Dy4wXf2PL2rwZcR1rjIkwxkQbY+q5XterZE7aLSLxAK7HPZX8nLD0v/9r53c65FWDnwpGCQkJbNu2rfD19u3bOfvssx2M6HS5ubn079+fQYMG0a9fP6fDKZSZmcmcOXNo3rw5AwcOZMWKFfzpT4Ez/2hCQgIJCQl06dIFgAEDBvD99987HJV3vOndhIhcKyIvuLbeVTjfHGCw6/lgIDAaW4PEBRfA5s1w991OR6J85eKLL2bDhg388ssvnDx5klmzZnHttdc6HRZgWx2GDBlCYmIiwwNscfZnnnmG7du3k5WVxaxZs+jYsSP//ve/nQ6r0FlnnUXTpk1Zv349APPnz+eCIFkXwJsR12OA+4C1ru0+V1l575sJ/BdoIyLbRWQIMAZIF5ENQLrrtfJSly62p9OMGXZToScqKoqXXnqJq666isTERK6//nratWvndFiA/d/69OnTWbBgAUlJSSQlJTF37lynwwoaL774IoMGDaJDhw6sXLmSxx57zOmQvOLNjete2In98gFEZCqwgnK6rxpjbvSwq0eFIlSneewx+OILGDoUkpOhbVunI1LVrVevXvTq1cvpMErp1q1bwb3EgJaamsozzwTeBNVJSUlBN004eNncBDQo9ry+D+JQXoqKgpkz7XoTS5Y4HY1SKtR5cyXxDLBCRBZiJ/i7DAiO66QQlZBgJwGsr+laKeVj5SYJY8xMEVkEXIxNEn81xvzq68BU2QoSxJdfwr59MHBg2ccrpVRleHPjer4xZpcxZo4xZrYx5lcRme+P4FTZjLHrT9x6K6xZU/7xSilVUR6ThIjUEpFGQJyINHRNztdIRJoDgdNxO4yJwPTpUK8e9OsHBw86HZFSKtSUdSUxFFgOtHU9FmyzgUm+D015Iz4e3nkHfvkFbrpJ53dSgePgwYO8/PLLHvcfO3aMyy+/nLy8PBYtWkTv3u6HYA0cOJANGzb4KkxVDo9Jwhgz0RjTAnjIGNPSGNPCtV1kjHnJjzGqclx6qZ1W/NNP4T//cToapazyksSUKVPo168fkZGRZX7OnXfeyXPPPVfd4SkvldXcdLGInGWMedH1+mYRmS0i/3Q1Q6kAcscdNknc6Gl0ilJ+NmLECDZt2kRSUhIPP/xwqf0zZsygT58+ha9zcnIK11sYNGhQ4ZiMSy+9lC+//JJTuvqWI8pqbnoNOAkgIpdhR0dPAw4Bk30fmqqonj3tfYqNG+GHH5yORoW7MWPGcN5557Fy5Uqef/750/adPHmSzZs307x588KyFStWMGHCBNauXcvmzZvJzMwEICIiglatWrFq1Sp/hq9cykoSkcaYgvUgbgAmG2PeM8b8HWjl+9BUZRgD/ftD796wc6fT0Sjl3r59+0qtp9C5c2cSEhKIiIggKSmJrKyswn1nnnkmO/UX2hFlJgkRKRhH0QNYUGyfN4PwlANE4M03Yf9++MMf4Nixstt7lXJCTEwMx48fP62sZs2ahc8jIyNPa146fvy4T9ZeUOUrK0nMBL4SkdnAMWAJgIi0wjY5qQDVsSO8/TasXAmjRl2ANuUqJ8TGxpKdne12X8OGDcnLyyuVKDz5+eefA2aiw3BTVu+m0cCD2LWqu5mimb0igHt8H5qqil69YNIkWLq0MQG4Xr0KA40bN+aSSy6hffv2bm9cX3nllXz99dflfs7u3buJiYkhPj7eF2GqcpTZbGSMWeqm7GffhaOq0x13wLZta7nrruCYt16FnrfeesvjvmHDhjFu3DiuuOIKUlNTSU1NLdz30ktFvezfeusthg4d6sswVRm8nQVWBan09D3ExNjR2PPmOR2NUkU6duxIWloaeeWMAG3QoAGDBw8u8xjlO5okwsTIkbbH08cfOx2JUkVuvfXWcgfT3XLLLURFaV8Zp2iSCBNjxkCnTvDHP8JXXzkdjVIqWGiSCBOxsTB3LrRoYbvGBska7Eoph2mSCCNxcfD559CoEWRk2IF3SilVFm3oCzMJCXahopo17cA7pZQqi15JhKFWraBpU8jPh0cegZ9+cjoipVSg0iQRxnbuhGnTIDVVE4VSyj1NEmEsIQEWLrTPNVEopdzRJBHmEhOLEsXll9v5npRSqoAmCUVioh07Ub++nT1WKaUKaO8mBUCbNvDjjxAdbV9v22ZvbiulwpteSahCBQli9mzbA2rWLGfjUUo5T5OEKiUtDbp2hZtuggkTdNCdUuFMk4QqpV49O2Ns377wwANw773owkVKhSlNEsqtmBh45x148EF46SX47DOnI1JKOUGThPIoIgJeeAH+7//gmmtsWTlT/yulQowmCVWu3/3OPi5fDhdcAN9+62w8Sin/0SShvBYdDSdPwqWXwpQpTkejlPIHTRLKax06wLJlNkkMGQJ3322ThlIqdGmSUBXSuLHt+fTQQ/Dyy/Dqq05HpJTyJUdGXItIFpAN5AGnjDEpTsShKicqCp5/HtLToXt3W5aTA3XrOhuXUqr6OXklkWaMSdIEEbyuvNImjH377A3tv/9dx1MoFWq0uUlVWUyMvap46il7ZbF9u9MRKaWqixgH5lwQkV+AA4ABXjPGTHZzTAaQAdCkSZPkWZWcSCgnJ4e6YdwO4s/6f/FFE8aNa02NGvkMH76eyy/f55fzehLO33041x3Cu/7F656Wlra8yq01xhi/b8DZrsczgVXAZWUdn5ycbCpr4cKFlX5vKPB3/devNyY52ZirrzYmP9+vpy4lnL/7cK67MeFd/+J1B5aZKv69duTGtTFmp+txj4h8AHQGFjsRi6perVvDf/9rb2SLQFYW/PyzvX+hlAo+fr8nISJ1RCS24DlwJbDG33Eo34mOhoYN7fOnnoKrroK77oLsbGfjUkpVnBM3rpsAX4vIKuBb4BNjzDwH4lB+8OKLMHy4HU9xwQUwZ47TESmlKsLvScIYs9kYc5Fra2eMGe3vGJT/xMTA2LF2ksAGDaBPH3j9daejUkp5S7vAKr/o2hW+/x7GjYPrr7dlO3fqrLJKBTpNEspvoqPtIkb16kFurr2ZnZICi7XLglIBS5OEckRUlB2h/dtvcPnlcMMNsHWr01EppUrSJKEcIWITw08/wRNP2BvabdrAd985HZlSqjhNEspRtWvDP/4B69fDsGHQsaMtX73aNkkppZylSUIFhGbN7MyyUVFw5Aj06AGJiTBrFuTnOx2dUuFLk4QKOLVrw5tvQp06cOONkJwMn3wCDkwzplTY0yShAo4I9OoFK1bAjBlw+DD07q1rayvlBE0SKmBFRMBNN9mb27NnQ5cutnzsWNsMpWMslPI9TRIq4EVHw7XX2ud5eTBzpm2GuuACeO01OHrU2fiUCmWaJFRQiYy0zU7vvguxsXDHHfam9+zZTkemVGjSJKGCTkQE9O9vx1QsXgyXXmrHWACsW2fvZSilqocmCRW0RGyC+OADaNvWlj31FHTqZOeK+vTTs7QpSqkq0iShQspLL8H48XDoEDz3XFvOPhuefNLpqJQKXpokVEhp2BDuvx/WroUJE1ZwzTVF+06dgjfegAMHHAtPqaCjSUKFJBG46KJDzJhh54YCWLIEbrsNzjoLBgywN7tPnnQ2TqUCnSYJFTZSU2HZMrjzTpswrrsOzj7brsGtlHJPk4QKGyJ2io8JE2DHDjvVx4ABcN55dv9TT9m1uBcssE1TSilNEipMRUXZqT9efdWOvQDYuxemTrWTC8bHw+23w8KFzsaplNM0SSjlMnGiTRTvvWdXzfvPf2zSADu54JtvwrZtjoaolN9FOR2AUoGkdm3o189ux4/brrRg54+65Rb7/MIL7YSDPXva8Rg1ajgXr1K+plcSSnlQqxY0aWKft20LP/4Izz0HjRrZx8svt/c1wN7jWL5cJx1UoUeThFJeELETCj78MCxaBPv2wYcfQvfudv/06ZCSAmecYacMmTQJVq3SBZNU8NPmJqUqoUED6NOn6PUtt0DTpjB/vt3ef982Qx06ZK9IvvzSJpouXaBuXcfCVqrCNEkoVQ2aNIFBg+xmDGRl2fsYtWrZ/f/4B2Rm2p5USUn2XkZqqu2Cq1Qg0yShVDUTgRYt7Fbg449h6VKbKDIzba+pLVuKksSNN9pEk5Jit9at7Wy3SjlNk4RSftCgge0N1bOnfZ2fDwcP2ue5ubB9O8yZY7vhgu1lNXIkPPaYvRn+1Ve2V9UZZzgRvQpnmiSUckBEhO0lBXblvSVL7Cjvn36yU4esXAnt2tn9mzfbAX5g55268EK73XwzXHSRbd4CewWjVHXTJKFUgIiKgvbt7VbcOefA55/D6tV2++EHePlluOQSmyQWLrQ9qtq0Kdpyc+Po0KEoESlVWZoklApwtWtDerrdCuTlFXWvjYuDgQNh/Xrbi2raNID2XH01dO5sR5BPmFB0n6RFC2jZ0u4ruLGulCeaJJQKQpGRRXNOdegAr7xStC87G956axkXXpgC2KatiAg7vuPf/y5qntq1yzZfvfKKXd0vIeH07aqris6hwpcmCaVCTGwstGmTQ0yMfd23r93Arp+xdau9z1EwmlzEjudYu9Ymjvx8e5/k+HG7/957Yd68ouRxzjnQrJmdch1g5047JqRRI+2RFYo0SSgVRmrUgFat7FbgjjvsBvbm+a+/wu7dRX/wO3SAPXtsD6zFi+0UJAkJRUnillvsPZOoKNv7qkkTOyX766/b/e+8A0ePQuPGNpE0bgxnnmlXEVSBT5OEUqpQVFTRFUOB226zWwFjbJNWgeHD7bTru3cXbcV7Wj35pJ33qrju3e3IdLD3Wg4ftsmjIJF07WrHjgB8+inUrAn165++6cSK/qFJQilVISJQr17R66uuspsnmZnw22+nb8WvIlq2tCPU9+yxXYB/+w327y9KEgMG2CuR4m67Df71L5uwEhNtE1vxBNKiRRypqbZ57eWX7VQodetCnTr2sVUrO41KXp5d87xuXZuItBtxaY4kCRHpCUwEIoHXjTFjnIhDKeV7BX+4W7Z0v/+110qXFdxcB/j6a3vPpGA7fNjOygt2IGJSUlH5r7/a51ddVRuwAxYfeKD05z/9NDz6qL0/UxBXRERRInn6afjLX2DTJhgyxPYCi4kperztNvj97+37p049fV+tWnaG4IQEm4DWr7cJqGCrUcM2t9WsWTRrcCB3EPB7khCRSGASkA5sB74TkTnGmLX+jkUpFZiK/4++Y0fPx9WoAbNmlS5ftGgr0JK4OHtlcuSI3XJy7OO559rjGjSAf/6zqLzgmIL9eXk2YR04YG/QHz8Ox47Z9UTAdgB4/PHS55892yaJzEz4wx9K71+wANLS7P2aG2+0CaoggdSsCZ99ZpNfIHDiSqIzsNEYsxlARGYBfQBNEkqpalUwst3ToMKGDeGeezy/v3VrOyWKJ6mp9mrm2LGiBHL8uF3+Fuysv59+astOnLDNXydOFF0JtW8Po0advu/kycAaBCmm+HWdP04oMgDoaYy5zfX6z0AXY8ywEsdlABkATZo0SZ7l7r8LXsjJyaFuGM/NHM7117qHZ90hvOtfvO5paWnLjTEpVfk8J64k3N0aKpWpjDGTgckAKSkpJjU1tVInW7RoEZV9bygI5/pr3VOdDsMx4Vz/6q67E0NftgNNi71OAHY6EIdSSqlyOJEkvgPOF5EWIlIDGAjMcSAOpZRS5fB7c5Mx5pSIDAM+w3aBnWKM+bGctymllHKAI+MkjDFzgblOnFsppZT3dDoupZRSHmmSUEop5ZEmCaWUUh75fTBdZYjIXmBLJd8eB+yrxnCCTTjXX+sevsK5/sXrfq4x5oyqfFhQJImqEJFlVR1xGMzCuf5a9/CsO4R3/au77trcpJRSyiNNEkoppTwKhyQx2ekAHBbO9de6h69wrn+11j3k70kopZSqvHC4klBKKVVJmiSUUkp5FNJJQkR6ish6EdkoIiOcjscXRCRLRFaLyEoRWeYqayQiX4jIBtdjw2LHP+r6eawXkTKWrw88IjJFRPaIyJpiZRWuq4gku35mG0XknyLibo2TgOOh/v8QkR2u73+liPQqti9k6i8iTUVkoYisE5EfReQ+V3nIf/9l1N0/370xJiQ37Ayzm4CWQA1gFXCB03H5oJ5ZQFyJsueAEa7nI4BnXc8vcP0cagItXD+fSKfrUIG6XgZ0AtZUpa7At8DvsAtgfQpc7XTdqlD/fwAPuTk2pOoPxAOdXM9jgZ9ddQz577+Muvvluw/lK4nCtbSNMSeBgrW0w0EfYKrr+VTgumLls4wxJ4wxvwAbsT+noGCMWQzsL1FcobqKSDxQzxjzX2P/1Uwr9p6A5qH+noRU/Y0xu4wx37ueZwPrgHMIg++/jLp7Uq11D+UkcQ6wrdjr7ZT9gw1WBvhcRJa71gUHaGKM2QX2Fww401Ueij+Titb1HNfzkuXBbJiI/OBqjipobgnZ+otIc6Aj8A1h9v2XqDv44bsP5STh1VraIeASY0wn4GrgbhG5rIxjw+VnAp7rGmo/g1eA84AkYBcw1lUekvUXkbrAe8D9xpjDZR3qpiyo6++m7n757kM5SYTFWtrGmJ2uxz3AB9jmo92uS0tcj3tch4fiz6Sidd3uel6yPCgZY3YbY/KMMfnAvyhqPgy5+otINPaP5AxjzPuu4rD4/t3V3V/ffSgniZBfS1tE6ohIbMFz4EpgDbaeg12HDQZmu57PAQaKSE0RaQGcj72RFcwqVFdXk0S2iHR19ey4udh7gk7BH0iXvtjvH0Ks/q5Y3wDWGWPGFdsV8t+/p7r77bt3+s69j3sF9ML2BNgEjHQ6Hh/UryW2F8Mq4MeCOgKNgfnABtdjo2LvGen6eawnwHt1uKnvTOxldS72f0VDKlNXIMX1D2oT8BKumQcCffNQ/+nAauAH1x+H+FCsP9AN2zTyA7DStfUKh++/jLr75bvXaTmUUkp5FMrNTUoppapIk4RSSimPNEkopZTySJOEUkopjzRJKKWU8kiThAorItJARO4q9vpsEXnXR+e6TkQedz1/U0QGuDnmDBGZ54vzK1UdNEmocNMAKEwSxpidxphSf7yrySPAy2UdYIzZC+wSkUt8FINSVaJJQoWbMcB5rvn3nxeR5uJan0FE/iIiH4rIRyLyi4gME5HhIrJCRJaKSCPXceeJyDzXpIpLRKRtyZOISGvghDFmX7Hiy0Tk/0Rkc4mrig+BQb6rslKVp0lChZsRwCZjTJIx5mE3+9sDN2HnwRkNHDXGdAT+i53GAOxC8/cYY5KBh3B/tXAJ8H2Jsnjs6Nne2GRVYBlwaeWqo5RvRTkdgFIBZqGxc/Zni8gh4CNX+Wqgg2smzt8D7xRb1Kumm8+JB/aWKPvQ2MnY1opIk2Lle4Czq6sCSlUnTRJKne5Esef5xV7nY/+9RAAHjTFJ5XzOMaB+GZ9dfNrmWq7jlQo42tykwk02dgnISjF2Hv9fROSPYGfoFJGL3By6Dmjl5ce2pmgGT6UCiiYJFVaMMb8BmSKyRkSer+THDAKGiEjB7LvulsVdDHT0aqF5SAM+qWQsSvmUzgKrlI+IyETgI2PMl+UctxjoY4w54J/IlPKeXkko5TtPA7XLOkBEzgDGaYJQgUqvJJRSSnmkVxJKKaU80iShlFLKI00SSimlPNIkoZRSyiNNEkoppTz6fxNUcvjRfnSOAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"\n",
"fig, ax1 = plt.subplots()\n",
"\n",
"# These are in unitless percentages of the figure size. (0,0 is bottom left)\n",
"left, bottom, width, height = [0.55, 0.5, 0.3, 0.3]\n",
"ax2 = fig.add_axes([left, bottom, width, height])\n",
"\n",
"ax1.plot(t,S,'b--')\n",
"ax1.set_xlabel('time (h)')\n",
"ax1.set_ylabel('Storage (in)')\n",
"ax1.grid()\n",
"ax2.plot(t[0:25],S[0:25],'go-')\n",
"ax2.set_xlabel('t (h)')\n",
"ax2.set_ylabel('S (in)')\n",
"#plt.show()\n",
"plt.savefig('Storage.pdf')"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "76ce144f",
"metadata": {
"id": "76ce144f"
},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.10.12"
},
"colab": {
"provenance": [],
"include_colab_link": true
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment