EQE511-Week00.ipynb

{
  "metadata": {
    "kernelspec": {
      "display_name": "Python 3.6",
      "name": "python36",
      "language": "python"
    },
    "language_info": {
      "mimetype": "text/x-python",
      "nbconvert_exporter": "python",
      "name": "python",
      "pygments_lexer": "ipython3",
      "version": "3.6.6",
      "file_extension": ".py",
      "codemirror_mode": {
        "version": 3,
        "name": "ipython"
      }
    }
  },
  "nbformat_minor": 2,
  "nbformat": 4,
  "cells": [
    {
      "cell_type": "markdown",
      "source": "<img src=\"./a-figures/gtuinsaat_1.jpg\" width=\"100\" height=\"100\" align=\"middle\" >\n\n### \"GTU Department of Civil Engineering - Graduate Program\"\n## \"EQE511 FINITE ELEMENT METHODS AND COMPUTER APPLICATIONS\"\n#### Lecturer: _Dr. Ahmet Anıl Dindar (adindar@gtu.edu.tr)_\n#### Period: _2019-Fall_\n\n---\n\n#### Table of the Content of the Course\n\n\n|Week|Date|Topic|Tasks|\n|:---: |:---:|:---|:---:|\n| 1 | 11.2 | Introduction to finite elements method , Review of matrix theory and mathematical concepts|-|\n| 2 | 18.2 | Global and local coordinate systems, Stiffness and local - global transformation matrix for truss element|-|\n| 3 | 25.2 | Constructing and solving 2D truss system using Python | HW1| \n| 4 |  4.3 | Stiffness and local - global transformation matrix for 2D beam and frame elements | - |\n| 5 | 11.3 | Constructing and solving 2D frame system using Python | HW2|\n| 6 | 18.3 | Stiffness and local - global transformation matrix for 3D frame elements| -|\n| 7 | 25.3 | Introduction to two dimensional elements (Plane stress, Plan strain), Shape functions numerical integration| HW3|\n| 8 |  1.4 | Mid-term exam | Term Project Starts|\n| 9 |  8.4 | Constructing and solving plane stress and plane strain systems using Python | - |\n|10 | 15.4 | Formulation for plate, membrane and shell elements |  HW4 |\n|11 | 22.4 | Formulation for three dimensional solid elements |- |\n|12 | 29.4 | Dynamic time history and dynamic modal analyses of structure| HW5| \n|13 |  6.5 | Case Study 1 (Hackathon)| -|\n|14 | 13.5 | Case Study 2 (Hackathon)| -|\n\n---\n \n ",
      "metadata": {}
    },
    {
      "cell_type": "code",
      "source": "print(\"Galatasaray is cool\")",
      "metadata": {
        "trusted": true
      },
      "execution_count": 1,
      "outputs": [
        {
          "text": "Galatasaray is cool\n",
          "name": "stdout",
          "output_type": "stream"
        }
      ]
    },
    {
      "cell_type": "markdown",
      "source": "This is a text cell",
      "metadata": {}
    },
    {
      "cell_type": "code",
      "source": "import matplotlib.pyplot as plt",
      "metadata": {
        "trusted": true
      },
      "execution_count": 2,
      "outputs": []
    },
    {
      "cell_type": "code",
      "source": "plt.plot(range(100))",
      "metadata": {
        "trusted": true
      },
      "execution_count": 3,
      "outputs": [
        {
          "execution_count": 3,
          "output_type": "execute_result",
          "data": {
            "text/plain": "[<matplotlib.lines.Line2D at 0x7fd24b2c54e0>]"
          },
          "metadata": {}
        },
        {
          "output_type": "display_data",
          "data": {
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAD8CAYAAAB5Pm/hAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvqOYd8AAAIABJREFUeJzt3XlYlXX+//HnRwUFd1QUQcQFRQTccGsby5o0KzNzpmXKVqfZaraUTCebbNKmqalf2zhZYzNNG+CSZpvZXqZWsom7IoqCGyg7nM/vD05z+S1K5ZzDgfu8HtfFBefmPtzv2xtf3t7c54Wx1iIiIs7Vwt8DiIiIbynoRUQcTkEvIuJwCnoREYdT0IuIOJyCXkTE4RT0IiIOp6AXEXE4Bb2IiMO18vcAAF27drUxMTH+HkNEpFnZuHHjIWttt1Ot1ySCPiYmhg0bNvh7DBGRZsUYs+d01tOlGxERh1PQi4g4nIJeRMThThn0xpjnjDGFxpisk5aFGWPeMcZsc7/v7F5ujDGPG2O2G2MyjDHDfTm8iIic2umc0f8LmPCtZSnAGmttLLDG/RhgIhDrfpsBPO2dMUVEpKFOGfTW2g+BI99aPBlY4v54CXDFSctfsHU+BzoZYyK8NayIiJy5hl6j726tLQBwvw93L48E9p60Xr572XcYY2YYYzYYYzYUFRU1cAwRETkVb/8w1tSzrN7fVWitXWStTbbWJnfrdsr7/UVEHKW8qpYHV28m/2iZz7fV0BdMHTTGRFhrC9yXZgrdy/OBXietFwXs92RAERGn+XTHIVLSMsk7UkZU51CuH9Pbp9tr6Bn9CmC6++PpwPKTlt/gvvtmDFD8zSUeEZFAV1JRzd3pGVz7z3W0MPDyjDE+D3k4jTN6Y8xLwDigqzEmH7gXWAC8aoy5BcgDprlXfwO4BNgOlAE3+WBmEZFm552cg8xZlknR8Up+fl5ffnvhAEKCWzbKtk8Z9Nbaa77nU+PrWdcCv/J0KBERpzh0opJ5K7JZmVFAXI/2/POGZJKiOjXqDE2i1ExExGmstSz7eh/3vZ5DaWUNv79oALf/qB/BrRq/kEBBLyLiZfuPlXPP0kzWbiliWHQnHpqaRGz39n6bR0EvIuIlLpflxS/yWLg6l1qXZe6l8dx4VgwtW9R353njUdCLiHjBrkOlzErL4ItdRzinf1cevDKRXmGh/h4LUNCLiHikptbFsx/v4tF3thLcqgUPTU1iWnIUxvj3LP5kCnoRkQbK2V/CrLQMMvcV8+P47tx/RQLdO7Tx91jfoaAXETlDlTW1PPHedp5+fwedQoN48trhXJLYo0mdxZ9MQS8icgY27jnKrLQMthee4MrhkcydFE/ntsH+HusHKehFRE5DaWUNf31rC0s+203PjiH866aRjBsYfsrnNQUKehGRU/hoWxF3p2eSf7Sc6WN7c9eEONq1bj7x2XwmFRFpZMVl1cxflcNrG/Pp27Utr90+lpExYf4e64wp6EVE6vFW9gHmLMviSGkVvxzXjzvGx9ImqHFKyLxNQS8icpKi43UlZKsyC4iP6MDzN44kIbKjv8fyiIJeRIS6ErL0L/fx55U5lFfXctfFA5lxXl+CWjZ+CZm3KehFJODlHy1j9tIsPtxaxIjenVk4NYn+4e38PZbXKOhFJGC5XJb/rNvDwtW5WOC+ywdz/ZjetPBzCZm3KehFJCDtKDpBSloG63cf5dzYuhKyqM5No4TM2xT0IhJQqmtdLPpwJ4+t2UZIUEsenjaEqcMjm2x9gTco6EUkYGTtK2ZWWgbZ+0uYmNCD+yYPJrx90ysh8zYFvYg4XkV1LY+v2cY/PtxJ59BgnvnZcCYkRPh7rEajoBcRR1u/+wiz0jLYWVTKVSOimDspno6hQf4eq1Ep6EXEkU5U1vDXN3N54fM9RHYK4d+3jOLc2G7+HssvFPQi4jgfbC1idnom+4vLmT42hrsuHkjbZlRC5m2Bu+ci4jjHyqq4f+Vm0r7Mp1+3tqTePpYRvZtfCZm3KehFxBFWZxYwd3k2R8uq+PX5/fn1Bf2bbQmZtynoRaRZKzxewZ+WZfNm9gESIjuw5OaRDO7ZvEvIvE1BLyLNkrWW1I353L8yh4oaF7MmxHHbuX1o5YASMm9T0ItIs7P3SBmzl2by0bZDjIoJY8HURPp2c04Jmbcp6EWk2ah1WV74bDd/fWsLBrh/8mCuG+28EjJvU9CLSLOwvfA4s9Iy2bjnKOMGduOBKYlEdgrx91jNgoJeRJq06loX//hgB4+v2U5o65Y8+tMhXDHU2SVk3uZR0BtjfgfcClggE7gJiABeBsKAL4HrrbVVHs4pIgEoM7+YmWkZbC4oYVJSBPddPpiu7Vr7e6xmp8E/njbGRAJ3AMnW2gSgJXA1sBB41FobCxwFbvHGoCISOCqqa1mwOpcrnvqEwycq+cf1I3jy2uEK+Qby9NJNKyDEGFMNhAIFwAXAte7PLwHmAU97uB0RCRDrdh4mJT2TXYdKuXpkL+6+ZBAdQwKrhMzbGhz01tp9xpiHgTygHHgb2Agcs9bWuFfLByLre74xZgYwAyA6OrqhY4iIQxyvqGbhm7n85/M8eoWF8OKtozm7f1d/j+UIDQ56Y0xnYDLQBzgGvAZMrGdVW9/zrbWLgEUAycnJ9a4jIoFh7ZZC7knPpKCkglvO6cMffjyA0GDdK+ItnvxJXgjsstYWARhj0oGzgE7GmFbus/ooYL/nY4qIEx0treLPK3NY+tU+YsPbkfaLsxge3dnfYzmOJ0GfB4wxxoRSd+lmPLABWAtcRd2dN9OB5Z4OKSLOYq1lVWYB9y7Ppri8mjvGx/Kr8/vRupVKyHzBk2v064wxqdTdQlkDfEXdpZhVwMvGmPnuZYu9MaiIOMPBkgrmLsvi7ZyDJEV15D+3jmZQRAd/j+VoHl0Es9beC9z7rcU7gVGefF0RcR5rLa9u2Mv8VZupqnEx+5I4bj5bJWSNQT/tEBGfyztcRkp6Bp/uOMzoPmEsnJpETNe2/h4rYCjoRcRnal2Wf326m4ff2kLLFoYHpiRwzcholZA1MgW9iPjE1oPHmZmawdd7j3FBXDgPTEkgoqNKyPxBQS8iXlVV4+Lp93fwxNpttG8TxGNXD+XyIT1VQuZHCnoR8ZpNe48xKy2D3APHuWxIT+ZdFk8X9dP4nYJeRDxWXlXLo+9u5dmPdhLevg3P3pDMhfHd/T2WuCnoRcQjn+04TEp6BnsOl3HNqLoSsg5tVELWlCjoRaRBSiqqefCNXF76Io/eXUL5722jOaufSsiaIgW9iJyxNZsPcs/SLAqPV3DbuX34/UUDCQlWfUFTpaAXkdN2+EQl972ew4pN+xnYvT3PXD+Cob06+XssOQUFvYickrWWFZv2c9/rORyvqOa3F8byy3H9CW6l+oLmQEEvIj+ooLicucuyeHdzIUN6deKhqUkM7NHe32PJGVDQi0i9XC7Ly+v38uAbm6l2uZgzaRA3nd2HlqovaHYU9CLyHbsPlZKSnsHnO48wtm8XFkxNpHcXlZA1Vwp6EfmfmloXz32yi7+9vZXgli148MpErh7ZS/UFzZyCXkQAyD1QwqzUDDblF3PhoO7MvyKBHh3b+Hss8QIFvUiAq6yp5cm1O3hq7XY6hgTx/64ZxqVJETqLdxAFvUgA+yrvKDNTM9hWeIIpwyL506XxdG4b7O+xxMsU9CIBqKyqhr+9vZXnPtlFjw5teP7GkZwfF+7vscRHFPQiAeaT7YdISc9g75Fyrh/Tm5kTBtJeJWSOpqAXCRDF5dX8ZdVmXtmwlz5d2/LKjDGM7tvF32NJI1DQiwSAt7MPMGdZFodLq7j9R/347YWxtAlSCVmgUNCLONihE5XcuyKbVRkFxPVoz+LpI0mM6ujvsaSRKehFHMhay/Kv93Pf69mUVtbyxx8P4Oc/6kdQS5WQBSIFvYjD7D9Wzj1LM1m7pYhh0XUlZLHdVUIWyBT0Ig7hclle/CKPhatzqXVZ7r0snhvGxqiETBT0Ik6ws+gEKWmZfLH7COf078qDVybSKyzU32NJE6GgF2nGampdPPvxLh59ZyutW7XgoauSmDYiSvUF8n8o6EWaqZz9JcxM20TWvhIuHtyd+ycnEN5BJWTyXQp6kWamsqaWJ97bztPv76BTaDBPXzeciYkR/h5LmjCPgt4Y0wl4FkgALHAzsAV4BYgBdgM/sdYe9WhKEQFg454jzEzNYEdRKVcOrysh6xSqEjL5YZ7eVPsY8Ka1Ng4YAmwGUoA11tpYYI37sYh4oLSyhnkrsrnqmc+oqHax5OZRPPKToQp5OS0NPqM3xnQAzgNuBLDWVgFVxpjJwDj3akuA94FZngwpEsg+3FrE3emZ7C/+poQsjnatddVVTp8n3y19gSLgeWPMEGAjcCfQ3VpbAGCtLTDGqPtUpAGKy6qZvyqH1zbm07dbW179+VhGxoT5eyxphjwJ+lbAcOA31tp1xpjHOIPLNMaYGcAMgOjoaA/GEHGeN7MKmLs8myOlVfxyXD/uGK8SMmk4T4I+H8i31q5zP06lLugPGmMi3GfzEUBhfU+21i4CFgEkJydbD+YQcYzC4xXcuzyb1VkHiI/owPM3jiQhUiVk4pkGB7219oAxZq8xZqC1dgswHshxv00HFrjfL/fKpCIOZq0l7ct93L8yh/LqWu66eCAzzuurEjLxCk9/ovMb4EVjTDCwE7iJujt5XjXG3ALkAdM83IaIo+UfLeOepVl8sLWIEb07s3BqEv3D2/l7LHEQj4LeWvs1kFzPp8Z78nVFAoHLZfnPuj0sXJ2LBea5S8haqIRMvEz3aIn4wY6iE6SkZbB+91HOG9CNv0xJIKqzSsjENxT0Io2outbFog938tiabYQEteThaUOYOjxSJWTiUwp6kUaSta+YWWkZZO8v4ZLEHsy7fDDh7VVCJr6noBfxsYrqWh5fs41/fLiTsLbBPPOz4UxIUAmZNB4FvYgPrd99hFlpGewsKmXaiCjmTIqnY2iQv8eSAKOgF/GBE5U1PPRmLi98tofITiH8+5ZRnBvbzd9jSYBS0It42Qdbi5jtLiG76ewY/vjjgbRVCZn4kb77RLzkaGkV96/KIf3LffTr1pbU28cyordKyMT/FPQiHrLWsjrrAH9ansWxsmp+c0F/fn1Bf1q3UgmZNA0KehEPFJZUMHd5Fm9lHyQhsgMv3Dya+J4d/D2WyP+hoBdpAGstr23MZ/7KHCprXKRMjOPWc/rQSiVk0gQp6EXO0N4jZcxemslH2w4xKiaMBVMT6dtNJWTSdCnoRU5Trcvywme7+etbWzDA/VckcN2oaJWQSZOnoBc5DdsLjzMzNYMv844xbmA3HpiSSGSnEH+PJXJaFPQiP6C61sU/PtjB42u2E9q6JY/+dAhXDFUJmTQvCnqR75GZX8xdqZvIPXCcSUkR3Hf5YLq2a+3vsUTOmIJe5Fsqqmt59N2tPPvRLrq0DeYf14/g4sE9/D2WSIMp6EVO8vnOw9ydnsmuQ6X8NLkXsycNomOISsikeVPQiwDHK6pZsDqXF9fl0SsshBdvHc3Z/bv6eywRr1DQS8Bbm1vI7KWZHCip4Oaz+/DHiwcQGqy/GuIc+m6WgHWktIo/v57Nsq/3ExvejvRfnMWw6M7+HkvE6xT0EnCstbyeUcC8FdmUlFdzx/hYfnV+P5WQiWMp6CWgHCiuYM6yLN7dfJCkqI48dNto4nqohEycTUEvAcFay8vr9/KXVZupqnVx98Q4blEJmQQIBb04Xt7hMlLSM/h0x2FG9wlj4dQkYrq29fdYIo1GQS+OVeuyPP/JLh5+ewtBLVrwwJQErhmpEjIJPAp6caStB+tKyL7ee4zxceHMn5JAREeVkElgUtCLo1TVuHjq/e08uXY77dsE8djVQ7l8SE+VkElAU9CLY2zae4yZqRlsOXicyUN78qdL4+miEjIRBb00f+VVtTzyzhYWf7yL8PZtWDw9mfGDuvt7LJEmQ0EvzdqnOw5xd3omew6Xce3oaFImxtGhjUrIRE7mcdAbY1oCG4B91tpLjTF9gJeBMOBL4HprbZWn2xE5WUlFNQ++kctLX+TRu0soL902hrH9uvh7LJEmyRuvFrkT2HzS44XAo9baWOAocIsXtiHyP+/mHOSiRz7glfV5zDivL2/eeZ5CXuQHeBT0xpgoYBLwrPuxAS4AUt2rLAGu8GQbIt84fKKSO176iltf2EDn0GCW/vJsZl8yiJBgddSI/BBPL938HZgJtHc/7gIcs9bWuB/nA5EebkMCnLWWFZv2M29FNicqa/jdhQP4xbh+BLdSfYHI6Whw0BtjLgUKrbUbjTHjvllcz6r2e54/A5gBEB0d3dAxxOEKisuZszSLNbmFDO3ViYeuSmJA9/anfqKI/I8nZ/RnA5cbYy4B2gAdqDvD72SMaeU+q48C9tf3ZGvtImARQHJycr3/GEjgcrnqSsgefGMzNS7L3EvjufGsGFqqvkDkjDU46K21dwN3A7jP6P9orb3OGPMacBV1d95MB5Z7YU4JILsPlZKSnsHnO49wVr8uLLgyieguof4eS6TZ8sV99LOAl40x84GvgMU+2IY4UE2ti+c+2cXf3t5KcKsWLJyayE+Se6m+QMRDXgl6a+37wPvuj3cCo7zxdSVw5B4oYVZqBpvyi7kovjvzr0ige4c2/h5LxBH0yljxq8qaWp5cu4On1m6nY0gQT1w7jEmJETqLF/EiBb34zZd5R5mVmsG2whNMGRbJny6Np3PbYH+PJeI4CnppdGVVNfzt7a0898kuIjq04fkbR3J+XLi/xxJxLAW9NKpPth8iJT2DvUfK+dmYaGZNiKO9SshEfEpBL42iuLyav6zazCsb9tKna1temTGG0X3VTyPSGBT04nNvZx9gzrIsDpdWcfuP+vHbC2NpE6R+GpHGoqAXnyk6Xsm817NZlVHAoIgOLJ4+ksSojv4eSyTgKOjF66y1LPt6H/e9nkNZZS1/uGgAt4/rR1BLlZCJ+IOCXrxq37Fy7lmayftbihgeXVdC1j9cJWQi/qSgF69wuSwvfpHHgjc247Lwp0vjma4SMpEmQUEvHttZdIKU9Ey+2HWEc/p35cErE+kVphIykaZCQS8NVlPr4tmPd/HoO1tp3aoFD01NYlpylOoLRJoYBb00SM7+EmalZZC5r5iLB3fn/skJhKuETKRJUtDLGamoruWJ97bzzAc76BQazNPXDWdiYoS/xxKRH6Cgl9O2cc8RZqZmsKOolKnDo5h76SA6haqETKSpU9DLKZVW1vDXt7aw5LPd9OwYwpKbR/GjAd38PZaInCYFvfygj7YVcXd6JvuOlXPDmN7cNSGOdq31bSPSnOhvrNSruKya+1flkLoxn77d2vLqz8cyMibM32OJSAMo6OU73swqYO7ybI6UVvHLcf24Y7xKyESaMwW9/E/h8QruXZ7N6qwDxEd04PkbR5IQqRIykeZOQS9Ya0n7ch/3r8yhvLqWuy4eyIzz+qqETMQhFPQBLv9oGbOXZvHh1iKSe3dmwdQk+oe38/dYIuJFCvoA5XJZ/v35Hha+mQvAnycP5meje9NCJWQijqOgD0A7ik4wKzWDDXuOcm5sXQlZVGeVkIk4lYI+gFTXulj04U4eW7ONkKCW/G3aEK4cHqkSMhGHU9AHiKx9xcxMzSCnoIRLEntw3+UJdGvf2t9jiUgjUNA7XEV1LY+t2caiD3cS1jaYZ342nAkJKiETCSQKegdbv/sIs1Iz2HmolGkjopgzKZ6OoUH+HktEGpmC3oFOVNbw0Ju5vPDZHqI6h/DvW0ZxbqxKyEQClYLeYdZuKeSe9EwKSiq48awY7rp4IG1VQiYS0JQADnG0tIr7V+aQ/tU++oe3I/X2sxjRu7O/xxKRJqDBQW+M6QW8APQAXMAia+1jxpgw4BUgBtgN/MRae9TzUaU+1lreyDzAvSuyOFZWzW8u6M+vL+hP61YqIROROp6c0dcAf7DWfmmMaQ9sNMa8A9wIrLHWLjDGpAApwCzPR5VvO1hSwdxlWbydc5DEyI68cPNo4nt28PdYItLENDjorbUFQIH74+PGmM1AJDAZGOdebQnwPgp6r7LW8uqGvcxftZmqGhcpE+O49Zw+tFIJmYjUwyvX6I0xMcAwYB3Q3f2PANbaAmNMuDe2IXXyDpcxe2kmH28/xKg+YSycmkSfrm39PZaINGEeB70xph2QBvzWWltyui+nN8bMAGYAREdHezqG49W6LP/6dDcPv7WFli0M869I4NpR0SohE5FT8ijojTFB1IX8i9badPfig8aYCPfZfARQWN9zrbWLgEUAycnJ1pM5nG7bwePMTMvgq7xjnD+wGw9MSaRnpxB/jyUizYQnd90YYDGw2Vr7yEmfWgFMBxa43y/3aMIAVlXj4pkPdvDEe9tp27olf//pUCYP7akSMhE5I56c0Z8NXA9kGmO+di+bTV3Av2qMuQXIA6Z5NmJgyswv5q7UTeQeOM5lQ3py72XxdG2nEjIROXOe3HXzMfB9p5bjG/p1A11FdS2PvruVf364k27tW/PPG5K5KL67v8cSkWZMr4xtQj7feZiUtAx2Hy7jmlG9SJk4iI4hKiETEc8o6JuA4xXVLFidy4vr8ogOC+W/t47mrP5d/T2WiDiEgt7P3ss9yD1LszhYUsGt5/ThDz8eSEiw6gtExHsU9H5ypLSKP7+ezbKv9zOgezueuu4shkWrhExEvE9B38istazMKGDeimxKKqq5c3wsvzq/P8GtVF8gIr6hoG9EB4ormLMsi3c3H2RIVEcWXjWauB4qIRMR31LQNwJrLS+v38tfVm2m2uXinksGcfM5fWip+gIRaQQKeh/bc7iUlLRMPtt5mDF9w1hwZRIxKiETkUakoPeRWpfl+U928fDbWwhq0YIHr0zk6pG9VF8gIo1OQe8DWw7UlZBt2nuM8XHhzJ+SQERHlZCJiH8o6L2oqsbFU+9v58m122nfJojHrxnGZUkROosXEb9S0HvJ13uPMSs1gy0HjzN5aE/uvWwwYW2D/T2WiIiC3lPlVbU88s4WFn+8i/D2bVg8PZnxg1RCJiJNh4LeA5/uOERKWiZ5R8q4dnQ0KRPj6NBGJWQi0rQo6BugpKKaB9/I5aUv8ojpEspLt41hbL8u/h5LRKReCvoz9G7OQe5ZlknR8UpmnNeX3104QCVkItKkKehP0+ETlcx7PYfXN+0nrkd7/nlDMklRnfw9lojIKSnoT8Fay4pN+5m3IpsTlTX8/qIB3P6jfiohE5FmQ0H/A/YfK2fOsizeyy1kaK9OPHRVEgO6t/f3WCIiZ0RBXw+Xy/LfL/JYsDqXWpdlzqRB3HS2SshEpHlS0H/LrkOlpKRlsG7XEc7u34UHpyQR3SXU32OJiDSYgt6tptbF4o938cg7Wwlu1YKHpiYxLTlK9QUi0uwp6IHNBSXMSssgI7+Yi+K7M/+KBLp3aOPvsUREvCKgg76yppYn39vOU+/voGNIEE9cO4xJiSohExFnCdig/zLvKLNSM9hWeIIrh0Uy99J4OquETEQcKOCCvqyqhr++tYV/fbqbiA5teP7GkZwfF+7vsUREfCaggv7jbYdISc8g/2g514/pzcwJA2mvEjIRcbiACPri8moeWJXDqxvy6dO1La/MGMPoviohE5HA4Pigfyv7AHOXZXG4tIpfjOvHneNjaROkEjIRCRyODfqi45XMW5HNqswCBkV0YPH0kSRGdfT3WCIijc5xQW+tZelX+/jzyhzKKmu56+KBzDivL0EtVUImIoHJJ0FvjJkAPAa0BJ611i7wxXa+bd+xcmanZ/LB1iJG9O7MwqmJ9A9XCZmIBDavB70xpiXwJHARkA+sN8assNbmeHtb33C5LP9Zt4eFq3OxwLzL4rlhbAwtVEImIuKTM/pRwHZr7U4AY8zLwGTAJ0G/o+gEKWkZrN99lHNju/KXKYn0ClMJmYjIN3wR9JHA3pMe5wOjfbAdXl2/lznLswgJasnD04YwdXik6gtERL7FF0FfX9La76xkzAxgBkB0dHSDNtSnW1vGx4Vz3+TBhLdXCZmISH18EfT5QK+THkcB+7+9krV2EbAIIDk5+Tv/EJyOkTFhjIwJa8hTRUQChi/uOVwPxBpj+hhjgoGrgRU+2I6IiJwGr5/RW2trjDG/Bt6i7vbK56y12d7ejoiInB6f3EdvrX0DeMMXX1tERM6MXi4qIuJwCnoREYdT0IuIOJyCXkTE4RT0IiIOZ6xt0GuVvDuEMUXAngY+vStwyIvjNBeBuN+BuM8QmPsdiPsMZ77fva213U61UpMIek8YYzZYa5P9PUdjC8T9DsR9hsDc70DcZ/DdfuvSjYiIwynoRUQczglBv8jfA/hJIO53IO4zBOZ+B+I+g4/2u9lfoxcRkR/mhDN6ERH5Ac066I0xE4wxW4wx240xKf6exxeMMb2MMWuNMZuNMdnGmDvdy8OMMe8YY7a533f296zeZoxpaYz5yhiz0v24jzFmnXufX3HXYDuKMaaTMSbVGJPrPuZjA+RY/879/Z1ljHnJGNPGacfbGPOcMabQGJN10rJ6j62p87g72zKMMcM92XazDfqTfgn5RCAeuMYYE+/fqXyiBviDtXYQMAb4lXs/U4A11tpYYI37sdPcCWw+6fFC4FH3Ph8FbvHLVL71GPCmtTYOGELd/jv6WBtjIoE7gGRrbQJ19eZX47zj/S9gwreWfd+xnQjEut9mAE97suFmG/Sc9EvIrbVVwDe/hNxRrLUF1tov3R8fp+4vfiR1+7rEvdoS4Ar/TOgbxpgoYBLwrPuxAS4AUt2rOHGfOwDnAYsBrLVV1tpjOPxYu7UCQowxrYBQoACHHW9r7YfAkW8t/r5jOxl4wdb5HOhkjIlo6Labc9DX90vII/00S6MwxsQAw4B1QHdrbQHU/WMAhPtvMp/4OzATcLkfdwGOWWtr3I+deLz7AkXA8+5LVs8aY9ri8GNtrd0HPAzkURfwxcBGnH+84fuPrVfzrTkH/Wn9EnKnMMa0A9KA31prS/w9jy8ZYy4FCq21G09eXM+qTjverYDhwNPW2mFAKQ67TFMf93XpyUAfoCfQlrpLF9/mtOP9Q7z6/d6cg/60fgm5Exhjgqgprn/bAAABaElEQVQL+RettenuxQe/+a+c+32hv+bzgbOBy40xu6m7JHcBdWf4ndz/tQdnHu98IN9au879OJW64HfysQa4ENhlrS2y1lYD6cBZOP94w/cfW6/mW3MO+oD4JeTua9OLgc3W2kdO+tQKYLr74+nA8saezVestXdba6OstTHUHdf3rLXXAWuBq9yrOWqfAay1B4C9xpiB7kXjgRwcfKzd8oAxxphQ9/f7N/vt6OPt9n3HdgVwg/vumzFA8TeXeBrEWtts34BLgK3ADuAef8/jo308h7r/smUAX7vfLqHumvUaYJv7fZi/Z/XR/o8DVro/7gt8AWwHXgNa+3s+H+zvUGCD+3gvAzoHwrEG7gNygSzg30Brpx1v4CXqfgZRTd0Z+y3fd2ypu3TzpDvbMqm7I6nB29YrY0VEHK45X7oREZHToKAXEXE4Bb2IiMMp6EVEHE5BLyLicAp6ERGHU9CLiDicgl5ExOH+P9/tdpAkncOaAAAAAElFTkSuQmCC\n",
            "text/plain": "<Figure size 432x288 with 1 Axes>"
          },
          "metadata": {
            "needs_background": "light"
          }
        }
      ]
    },
    {
      "cell_type": "code",
      "source": "",
      "metadata": {},
      "execution_count": null,
      "outputs": []
    }
  ]
}