lesson6.ipynb

{
  "nbformat": 4,
  "nbformat_minor": 0,
  "metadata": {
    "colab": {
      "name": "lesson6.ipynb",
      "provenance": [],
      "collapsed_sections": [],
      "include_colab_link": true
    },
    "kernelspec": {
      "name": "python3",
      "display_name": "Python 3"
    }
  },
  "cells": [
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "view-in-github",
        "colab_type": "text"
      },
      "source": [
        "<a href=\"https://colab.research.google.com/github/onolab-tmu/audio-processing-100-knocks/blob/master/lesson6.ipynb\" target=\"_parent\"><img src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open In Colab\"/></a>"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "oakad_2x5FF2",
        "colab_type": "text"
      },
      "source": [
        "# Lesson 6: スコープ、クラスとオブジェクト指向\n",
        "\n",
        "## スコープ\n",
        "\n",
        "プログラムの中の変数の有効範囲は、スコープと読んでます。\n",
        "\n",
        "パイソンでは、スコープは二つあります：\n",
        "* **ローカルスコープ** 関数の中\n",
        "* **グローバルスコープ** プログラム全体\n",
        "\n",
        "### 最初はローカルスコープの例"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "SGkDLcDE5V9I",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "def my_func():\n",
        "  # `a` is a local variable of function `my_func`\n",
        "  a = 10\n",
        "  print(\"a value is\", a)"
      ],
      "execution_count": 0,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "7Qmku-5G18sK",
        "colab_type": "text"
      },
      "source": [
        "`my_func`を呼ぶと、想像通りに行きます。"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "q2SYVx4swuB7",
        "colab_type": "code",
        "outputId": "b700f883-e8df-48d5-b9e9-0b22292f5028",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 36
        }
      },
      "source": [
        "my_func()"
      ],
      "execution_count": 0,
      "outputs": [
        {
          "output_type": "stream",
          "text": [
            "a value is 10\n"
          ],
          "name": "stdout"
        }
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "ZHv4Dfw62GKi",
        "colab_type": "text"
      },
      "source": [
        "けれども、`my_func`の外から、`a`にアクセスしようとするとエラー発生です。"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "3GHwJ4_4wvhg",
        "colab_type": "code",
        "outputId": "216ea4d1-2033-47a9-982e-de8cfec6d3c5",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 174
        }
      },
      "source": [
        "print(a)"
      ],
      "execution_count": 0,
      "outputs": [
        {
          "output_type": "error",
          "ename": "NameError",
          "evalue": "ignored",
          "traceback": [
            "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
            "\u001b[0;31mNameError\u001b[0m                                 Traceback (most recent call last)",
            "\u001b[0;32m<ipython-input-3-bca0e2660b9f>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0mprint\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0ma\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
            "\u001b[0;31mNameError\u001b[0m: name 'a' is not defined"
          ]
        }
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "YLT82p0b9B0_",
        "colab_type": "text"
      },
      "source": [
        "### グローバルスコープのれい"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "m_Qf6eQ02meQ",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "b = 7\n",
        "\n",
        "def my_func2():\n",
        "  print(\"b value is\", b)"
      ],
      "execution_count": 0,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "rbk5wbrD2siO",
        "colab_type": "text"
      },
      "source": [
        "`my_func2`を走らせると、また思う通りに行きます。"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "SORpQnG03M4l",
        "colab_type": "code",
        "outputId": "5f07f69a-b2fd-43a2-a2e9-c9e1aefd99a0",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 36
        }
      },
      "source": [
        "my_func2()"
      ],
      "execution_count": 0,
      "outputs": [
        {
          "output_type": "stream",
          "text": [
            "b value is 7\n"
          ],
          "name": "stdout"
        }
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "xHOAQTiJ3Qq-",
        "colab_type": "text"
      },
      "source": [
        "そして、今回は"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "yUIuzr7S3Vn6",
        "colab_type": "code",
        "outputId": "00aa510f-5b8c-4054-9d8a-a2b3ce946e33",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 36
        }
      },
      "source": [
        "print(b)"
      ],
      "execution_count": 0,
      "outputs": [
        {
          "output_type": "stream",
          "text": [
            "7\n"
          ],
          "name": "stdout"
        }
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "dmPR7aaQ3Wdi",
        "colab_type": "text"
      },
      "source": [
        "エラー発生しません。\n",
        "\n",
        "という事で、\n",
        "* ローカルスコープの変数はグローバルスコーデアクセス**出来ません**.\n",
        "* グローバルスコープの変数はローカルスコープでアクセス**出来ます**.\n",
        "\n",
        "なう、グローバルで`b`を編集するとどうなるのでしょうか？\n",
        "\n",
        "\n"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "AZbVzWlm35SQ",
        "colab_type": "code",
        "outputId": "6773329c-a9b8-43cb-b52d-9c4906dd56ed",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 36
        }
      },
      "source": [
        "b = 36\n",
        "my_func2()"
      ],
      "execution_count": 0,
      "outputs": [
        {
          "output_type": "stream",
          "text": [
            "b value is 36\n"
          ],
          "name": "stdout"
        }
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "sPi30mXX4nML",
        "colab_type": "text"
      },
      "source": [
        "`my_func2`の中は本当にグローバルスコープの`b`を使ってるので、`b`を編集すると、関数の結果も変わります！"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "vkC8-gsp35vq",
        "colab_type": "text"
      },
      "source": [
        "### 変数名が被った時"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "U-tIFId237aW",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "c = 21\n",
        "\n",
        "def my_func3():\n",
        "  c = \"a string\"\n",
        "  print(\"c value is\", c)"
      ],
      "execution_count": 0,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "XjsFocRh5CgR",
        "colab_type": "code",
        "outputId": "748f8903-43a7-4965-b48e-1d6bcec32dce",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 36
        }
      },
      "source": [
        "my_func3()"
      ],
      "execution_count": 0,
      "outputs": [
        {
          "output_type": "stream",
          "text": [
            "c value is a string\n"
          ],
          "name": "stdout"
        }
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "eFxjzm1S5WgX",
        "colab_type": "code",
        "outputId": "79a8b62e-054c-4316-ce30-79447a71953b",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 36
        }
      },
      "source": [
        "print(c)"
      ],
      "execution_count": 0,
      "outputs": [
        {
          "output_type": "stream",
          "text": [
            "21\n"
          ],
          "name": "stdout"
        }
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "FdA5LkYz5E5k",
        "colab_type": "text"
      },
      "source": [
        "グローバルも、ローカルスコープも、どちらも同じ変数名が存在しますと、\n",
        "1. ローカルスコープで作成した変数は**別の変数**になっています。例のコードで、関数の中`c = \"a string\"`を書いた時、グローバルスコープの`c`の中身は変わっていませんでした。\n",
        "2. ローカルスコープの変数があると、グローバルスコープの同じ名前の変数は**アクセス出来なくなります**。(shadowing)\n",
        "\n",
        "もうし、関数の中にグローバルスコープの変数を変えたい場合は、`global`キーワードを使えます。"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "iM_zdARE590T",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "d = 73\n",
        "\n",
        "def my_func4():\n",
        "  global d\n",
        "  print(\"my_func4: d value was\", d)\n",
        "  d = \"I am a string now\"\n",
        "  print(\"my_func4: d value is\", d)"
      ],
      "execution_count": 0,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "jIoRdUxq6gQo",
        "colab_type": "code",
        "outputId": "7f1e9520-1876-4d2b-eb0b-957d67219b2e",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 54
        }
      },
      "source": [
        "my_func4()"
      ],
      "execution_count": 0,
      "outputs": [
        {
          "output_type": "stream",
          "text": [
            "my_func4: d value was 73\n",
            "my_func4: d value is I am a string now\n"
          ],
          "name": "stdout"
        }
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "WTqGAZ6L6hg8",
        "colab_type": "code",
        "outputId": "94fe5cc0-876d-49cf-d0b5-80d76059c81f",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 36
        }
      },
      "source": [
        "print(d)"
      ],
      "execution_count": 0,
      "outputs": [
        {
          "output_type": "stream",
          "text": [
            "I am a string now\n"
          ],
          "name": "stdout"
        }
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "Hbp5Zw-Q6icL",
        "colab_type": "text"
      },
      "source": [
        "## クラスとオブジェクト指向\n",
        "\n",
        "オブジェクト指向は複雑データタイプと、このタイプに対するオペレーションのプログラムの書き方です。\n",
        "\n",
        "例えば、簡単な物理シミュレーションを考えましょう。ボールがあって、投げた時、ニュートン力学に関して、どこまで行くか、プログラムで書きたい。\n",
        "\n",
        "ボールの必要なプロパーティは\n",
        "* 半径\n",
        "* 位置\n",
        "\n",
        "投げた時は速度のベクトルと飛ぶ時間で、新しい位置を計算するのが目的です。\n",
        "\n",
        "とりあえず、オブジェクトを使わず書いて見ましょう。\n"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "UGBg_hPi7168",
        "colab_type": "code",
        "outputId": "ce8be4d8-45a0-4b0d-f185-878679d3a43c",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 36
        }
      },
      "source": [
        "import math\n",
        "\n",
        "my_ball_position = 10.\n",
        "my_ball_radius = 0.5\n",
        "a_earth = 9.81\n",
        "\n",
        "def compute_fall_time(position, radius, acceleration):\n",
        "  \"\"\"\n",
        "  compute the ball new location after fly time, when thrown at speed\n",
        "  \"\"\"\n",
        "  \n",
        "\n",
        "  fall_time = math.sqrt(2 * (position - radius) / acceleration)\n",
        "\n",
        "  return fall_time\n",
        "\n",
        "t_fall = compute_fall_time(my_ball_position, my_ball_radius, a_earth)\n",
        "print(f\"The ball falls for {t_fall} seconds\")"
      ],
      "execution_count": 0,
      "outputs": [
        {
          "output_type": "stream",
          "text": [
            "The ball falls for 1.3916893275819882 seconds\n"
          ],
          "name": "stdout"
        }
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "tBfE5NDvDVTJ",
        "colab_type": "text"
      },
      "source": [
        "こう書くと、色々不便なところはあります。\n",
        "* 位置と半径は別の変数が入ります。\n",
        "* ボールがいっぱいあると、ちょっと変数の扱いが大変になります。\n",
        "* `compute_fall_time`関数は`position`と`radius`と対応してます。\n",
        "* `a_earth`はどう扱えばいいのか？\n",
        "\n",
        "オブジェクト指向使えば、もっとまとまった書き方があります。\n",
        "その時、新しいタイプを作成します。\n"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "-rklfqLMC5-U",
        "colab_type": "code",
        "outputId": "b32f4885-eeea-45bd-f81d-5b73f377be23",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 110
        }
      },
      "source": [
        "class Ball:\n",
        "\n",
        "  def __init__(self, pos, rad, accel):\n",
        "    self.position = pos\n",
        "    self.radius = rad\n",
        "    self.acceleration = accel\n",
        "\n",
        "  def __repr__(self):\n",
        "    return f\"This is a ball with position {self.position}, radius {self.radius}, and acceleration {self.acceleration}\"\n",
        "\n",
        "  def fall_time(self):\n",
        "    a_earth = 9.81\n",
        "    fall_time = math.sqrt(2 * (self.position - self.radius) / self.acceleration)\n",
        "    return fall_time\n",
        "\n",
        "class EarthBall(Ball):\n",
        "\n",
        "  def __init__(self, pos, rad):\n",
        "    super().__init__(pos, rad, 9.81)\n",
        "\n",
        "  def __repr__(self):\n",
        "    return f\"This is an Earth ball with position {self.position}, radius {self.radius}\"\n",
        "\n",
        "class MoonBall(Ball):\n",
        "\n",
        "  def __init__(self, pos, rad):\n",
        "    super().__init__(pos, rad, 1.635)\n",
        "\n",
        "  def __repr__(self):\n",
        "    return f\"This is a Moon ball with position {self.position}, radius {self.radius}\"\n",
        "\n",
        "my_ball = Ball(10, 0.5, 9.81)\n",
        "print(\"This is the ball radius:\", my_ball.radius)\n",
        "print(\"This is the ball position:\", my_ball.position)\n",
        "\n",
        "my_earth_ball = EarthBall(27, 1.2)\n",
        "my_moon_ball = MoonBall(27, 1.2)\n",
        "\n",
        "t_fall = my_ball.fall_time()\n",
        "print(f\"The ball falls for {t_fall} seconds\")\n",
        "print(f\"The ball falls on Earth for {my_earth_ball.fall_time()} seconds\")\n",
        "print(f\"The ball falls on Moon for {my_moon_ball.fall_time()} seconds\")"
      ],
      "execution_count": 0,
      "outputs": [
        {
          "output_type": "stream",
          "text": [
            "This is the ball radius: 0.5\n",
            "This is the ball position: 10\n",
            "The ball falls for 1.3916893275819882 seconds\n",
            "The ball falls on Earth for 2.293455654230203 seconds\n",
            "The ball falls on Moon for 5.6177961005649655 seconds\n"
          ],
          "name": "stdout"
        }
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "_5a42crWQphn",
        "colab_type": "code",
        "outputId": "c178b650-a99b-4b92-cd5e-35a2de0ea74a",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 54
        }
      },
      "source": [
        "print(my_moon_ball)\n",
        "print(my_ball)"
      ],
      "execution_count": 0,
      "outputs": [
        {
          "output_type": "stream",
          "text": [
            "This is a Moon ball with position 27, radius 1.2\n",
            "This is a ball with position 10, radius 0.5, and acceleration 9.81\n"
          ],
          "name": "stdout"
        }
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "JJkRis8Sa2aC",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "class Ball2:\n",
        "\n",
        "  def __init__(self, radius, position=0, speed=0, acceleration=9.81):\n",
        "    self.position = position\n",
        "    self.speed = speed\n",
        "    self.radius = radius\n",
        "    self.acceleration = acceleration\n",
        "\n",
        "  def fall_time(self, x):\n",
        "    discriminant = self.speed ** 2 - 4 * (0.5 * self.acceleration) * (self.positioin - x)\n",
        "\n",
        "    if discriminant < 0:\n",
        "      return math.isinf\n",
        "    elif discriminant == 0:\n",
        "      return -self.speed / self.acceleration\n",
        "    else:\n",
        "      sol1 = (math.sqrt(discriminant) - self.speed) / self.acceleration\n",
        "      sol2 = (math.sqrt(discriminant) + self.speed) / self.acceleration\n",
        "\n",
        "      return min(sol1, sol2)\n",
        "\n",
        "  def fly(self, t):\n",
        "\n",
        "    # update position and speed\n",
        "    self.position += self.speed * t + 0.5 * self.acceleration * t ** 2\n",
        "    self.speed += self.acceleration * t\n",
        "    "
      ],
      "execution_count": 0,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "m7JfnOW0iQMh",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "import numpy as np\n",
        "\n",
        "time = []\n",
        "position = []\n",
        "\n",
        "my_ball = Ball2(0.5, position=0, speed=-30, acceleration=9.81)\n",
        "\n",
        "step = 0.1\n",
        "\n",
        "for t in np.arange(0, 10, step):\n",
        "\n",
        "  position.append(-my_ball.position)\n",
        "  time.append(t)\n",
        "\n",
        "  my_ball.fly(step)"
      ],
      "execution_count": 0,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "yCYplwzri6Xy",
        "colab_type": "code",
        "outputId": "c2c667fd-006f-4b69-cdc0-d104080bad44",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 284
        }
      },
      "source": [
        "import matplotlib.pyplot as plt\n",
        "\n",
        "plt.plot(time, position)"
      ],
      "execution_count": 0,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "[<matplotlib.lines.Line2D at 0x7faecba933c8>]"
            ]
          },
          "metadata": {
            "tags": []
          },
          "execution_count": 51
        },
        {
          "output_type": "display_data",
          "data": {
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX8AAAD4CAYAAAAEhuazAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0\ndHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nO3deXxU9b3G8c83GyEJBEJCgIQlkICE\nsGkEFMQVRVRAbXHXuqF1bbVatVpbWpd6W723bhWtWy0idQOLguKGIoJBtrCHnRAgBEggIfvv/sHc\n3tSCAsnMmcw879drXpmc2Z7j8uTMOb/zO+acQ0REwkuE1wFERCTwVP4iImFI5S8iEoZU/iIiYUjl\nLyIShqK8DnC4kpOTXbdu3byOISLSbCxYsGCncy7lYI81m/Lv1q0beXl5XscQEWk2zGzjoR7Tbh8R\nkTCk8hcRCUN+3+1jZhuAvUAdUOucyzWzJOANoBuwARjnnNvt7ywiInJAoLb8T3XODXDO5fp+vwf4\n2DmXBXzs+11ERALEq90+Y4BXfPdfAcZ6lENEJCwFovwd8KGZLTCz8b5lqc65It/9bUDqwV5oZuPN\nLM/M8oqLiwMQVUQkPARiqOcw51yhmbUHPjKzlQ0fdM45Mzvo1KLOuYnARIDc3FxNPyoi0kT8Xv7O\nuULfzx1m9g4wCNhuZh2dc0Vm1hHY4e8c0jT2V9dRUl5Fyb5qSvfXsK+qlr2VNeyvrqOqtp6q2npq\n6x34pgo3M2KiIoiJjKBFdAQJLaKIbxFF69hokuJjSIqPoW1cNFGRGngmEkh+LX8ziwcinHN7fffP\nBCYA04CrgEd9P6f6M4ccvsqaOtYVl7OhpJyNJRVsLCmncM9+ikor2VZayb6q2sN+L7N//Q34weel\nJLSgQ2IsHVrH0iUpjq7t4ujSLp7M9gl0SozFzBqxViLyXf7e8k8F3vH9jxsFTHLOzTCzb4ApZnYt\nsBEY5+ccchDbyyrJLywlv7CM5UWlrN6+j40l5dQ3KOx28TGkt21Jj5R4hmUm0751C9rFx5AU34K2\ncdEkxEbRKjaauOhIWkQf2MJvuBXvnKO6rp7q2noqa+opr6plX1UtZftr2FVRze7yaor3VbO9tJKi\nskrW7yzn89XFVNXW/+s94mIiyWyfQO8OremT1prsjq3p0ymRljGRgfzHJRJSrLlcySs3N9dpeoej\nV1fvWFFUxrz1u/h2424WbtrN1tJK4MCWd0a7eHp1aEVWaiuy2ifQPSWeLklxtIqNDnjW+npH8b4q\nNuwsp6B4H2u272PNjr0s31rG7ooaACIjjN4dWzGwc1tyu7VlSPd2pLaODXhWkWBmZgsaDLH/98dU\n/qHJOcfa4nJmry5mTsFO5m/Yxd7KA7ts0tq05NiubRnQuQ390hPp3bE1CS2Cf5on5xxFpQe+rSze\nsoeFm/awePMeyqvrAOjWLo4TeiQzPCuZE3skkxgX+D9cIsFE5R8mqmrrmLu2hFkrtvPpymIK9+wH\n/r8Uh3RPYnBGOzokhs4Wcm1dPSuK9jJvfQlfr9vFvHUl7K2qJcJgYJe2nN67PWf0TiWrfYKOG0jY\nUfmHsP3VdXy6agfTlxbx2codlFfXERcTybDMZE7ulcLwrBQ6J8V5HTNgaurqWbx5D7NXF/PJqh3k\nF5YB0LVdHGfndGRU3w70TUvUHwIJCyr/EFNTV88Xa4p5Z+FWZi3fzv6aOpITYhiR3YEz+6RyQvd2\nxEbrYCjAttJKPl65nRn52/hqbQl19Y7OSS0Z0z+NsQM7kdm+ldcRRfxG5R8iVhSVMSVvM9MWbaWk\nvJq2cdGM6tuRc/p1ZHBGOyIjtDX7fXaXV/PR8u28t2Qrcwp2Uu+gb1oi43LTGd0/TccIJOSo/Jux\n8qpapi7ayuRvNrFkSykxkRGMyE7l/IFpDO+ZQkyUTo46GjvKKnlvSRFvLtjCiqIyYqIiODunA5cN\n7srx3dpqt5CEBJV/M1SwYy+vfb2JtxZsYW9VLcd0aMVFx3dm7IA02sbHeB0vpOQXljIlbzPvLCxk\nb2UtPVMTuHxIVy48Np34ZjAKSuRQVP7NhHOOL9bs5K9frufz1cXEREZwTr+OXD6kK8d2aaOtUT+r\nqK7ln4uLeG3eRpZsKaVVbBSXDOrClSd0Jb1t+Bw0l9Ch8g9yNXX1/HPJVv7y2TpWbd9LSqsWXDmk\nK5cO7kK7hBZexwtL327azYtfrueD/G0AnNuvIzcM70F2p9YeJxM5fCr/IFVZU8eUvM089/k6Cvfs\np2dqAuOH9+C8/h1pEaXROsFg6579vDRnPZPmbaK8uo7hPVO47bRMcrsleR1N5Aep/INMZU0dr8/f\nxF8+X8v2siqO69qWm07pwam92hOhETtBqbSihtfmbeTFL9dTUl7NiT3acdvpWQzp3s7raCKHpPIP\nEtW19byRt5knP17Djr1VDM5I4vYzsjihezvtz28mKqprmTRvE8/NXkfx3iqGZSbzi7N6MaBzG6+j\nifwHlb/H6uodUxcV8sSs1WzetZ/ju7XljhG9OKGHthqbq8qaOl77eiPPfraWkvJqzuidyt0je9Ez\nVSeNSfBQ+XvoizXFPPz+SlYUldGnU2vuOqsXJ/dM0ZZ+iCivquWlOet5bvY6yqtq+fFxnfn5iJ4h\nNX+SNF8qfw+s2b6X301fwezVxaS3bcldZ/XivH6dtE8/RO0ur+apTwt4de4GIiOMG4b34MaTe+ia\nA+IplX8AlVbU8MSs1fzt643Ex0Ry62lZXHliV43eCRObSir4w8yVTF9SRMfEWO45+xhG9++kb3ri\nCZV/ANTXO97I28xjM1ZSur+GSwZ14Y4RPTVOP0zNX7+LCf9cRn5hGYO6JTFhbB+O6aBzBCSwVP5+\nll9Yyv3v5rNo8x4GdUviN6P76GQgoa7e8Y+8zfxhxkrKKmu56oRu/HxElidXR5Pw9H3lr4lLGqG8\nqpY/fbial79aT1J8DI+P68/5A9P0FV+AA5eavHhQF0bmdOC/Zq7ipa/WM33pVn47OoeROR28jidh\nTlv+R+nTlTu4/918Cvfs57LBXbh75DEkttQWnRzaos17uPftpawoKmNEdioTxvShY2JLr2NJCNNu\nnya0u7ya37y3jKmLtpLVPoFHLuirU/3lsNXU1fPil+t5YtZqoiMi+NU5vbno+M76tih+ofJvIh8s\nLeKBqfnsqajhltMyuemUTM2nL0dlY0k5v3xrCV+v28XQzHY8ekG/sLrcpgSGyr+R9lRU88DUZby3\neCs5aa35rx/1p3dHHdCVxqmvd0yav4lH3l8BwK/Py2Zcrr4FSNPRAd9G+GzVDu5+cwm7yqu5Y0RP\nfnpKD6IjtbUvjRcRYVw+pCun9Erhrn8s4ZdvLWXmsu08emFf2rfSGcLiX2qxQ6isqePXU/P5yUvf\nkNgymndvHsptp2ep+KXJpbeN4+/XDebX52Yzp2AnI//7Cz5avt3rWBLi1GQHsaKojNFPfcmrczdy\nzdAM3rt1GDlpiV7HkhAWEWFcMyyD6bcNo0PrWK5/NY9fvbOU/dV1XkeTEKXyb8A5xytfbWDM03PY\nVV7Dq9cM4tfnZRMbrakZJDAy27finZtP5Ibh3Zk0fxPnPvkFK7eVeR1LQpDK36e0ooYbX1vAg9OW\nMSwzmRk/O4nhPVO8jiVhqEVUJPeO6s1r1w6mrLKWMU/NYdK8TTSXwRnSPKj8gYWbdjPqz1/w8Yod\n3H9Ob/56VS7JmpNHPDY0M5n3bzuJQRlJ3PfOUm59fSH7qmq9jiUhIqzL//9284x7bi5m8OZPT+S6\nk7prqJ0EjZRWLXjl6kHcPbIX7y8tYvRTX7Jq216vY0kICNvyL6+q5bbJi3hw2jKGZ6Uw/daTdCk+\nCUoREcZNp2Qy6fohlO2vZczTX/L2t1u8jiXNXFiW//qd5Zz/zBymL9nKXWf14vkrc0mM07w8EtyG\ndG/H+7cNo196G+6YspjfTFtGTV2917GkmQq78v905Q5GP/UlxXurePWawdx8aqauriXNRvvWsUy6\nbjDXDsvg5a82cNnz8yjeW+V1LGmGwqb8nXM89ckarnnlGzq3jWPaLcMYlpXsdSyRIxYVGcED52bz\nPxcPYEnhHs578kuWbNnjdSxpZjwrfzMbaWarzKzAzO7x52ftr67j1tcX8scPVzOmfyfe+umJmkRL\nmr0xA9J4+6dDiYwwfvyXuUxdVOh1JGlGPCl/M4sEngbOBrKBS8ws2x+fta20knHPzWX60iLuOfsY\nnrhogC6qLSEju1Nrpt4ylP7pbbh98iIem7GS+nqdDyA/zKst/0FAgXNunXOuGpgMjGnqD9lTUc3o\np75kXfE+XrgylxtP7qFhnBJykhNa8Np1g7lkUGee+WwtN/39W00LIT/Iq/JPAzY3+H2Lb9m/MbPx\nZpZnZnnFxcVH/CFt4mK4ZlgGb980lNN7px59WpEgFxMVwcPn9+X+c3ozc/k2Lpo4lx1llV7HkiAW\n1Ad8nXMTnXO5zrnclJSjm2rhxpN70KtDqyZOJhJ8zIzrTurOxCtyKdixj7FPz9G8QHJIXpV/IdC5\nwe/pvmUi0kgjslOZcsMJ1DnHj5+dy5yCnV5HkiDkVfl/A2SZWYaZxQAXA9M8yiIScnLSEnnnpqF0\natOSq16crzOC5T94Uv7OuVrgFmAmsAKY4pxb5kUWkVDVqU1Lptx4AoMykrhjymKe+axAM4PKv3i2\nz985975zrqdzrodz7iGvcoiEssSW0bx89SDGDOjEYzNW8dv3lmsoqAC6hq9IyIuJiuCJcQNITmjB\nX79cz859VfxpXH9aROl8l3Cm8hcJAxERxv3n9KZ9qxY88sFKSvfX8NwVxxEXowoIV0E91FNEmo6Z\nccPJPXjswn7MKdjJ5S/Mo7SixutY4hGVv0iYGXd8Z5657FjyC8sOnAy2VyeDhSOVv0gYGpnTkRd/\ncjybdlVw0XNfs3XPfq8jSYCp/EXC1LCsZF69ZhA791Yx7rm5bCqp8DqSBJDKXySM5XZL4u/XD2Zf\nVS3jnpvL2uJ9XkeSAFH5i4S5fultmDx+CLX19Vz03NcU7NAF4sOByl9EOKZDayaPH4IZXDzxa1Zv\n1x+AUKfyFxEAMtu3YvL4IUSYcfHErzUjaIhT+YvIv/RISeCNG04gJjKCy56fp28AIUzlLyL/JiM5\nnknXDyYywrj0eR0DCFUqfxH5D91TEph0/RDAuOT5eRoFFIJU/iJyUJntE5g8fjDOOS57fp7OAwgx\nKn8ROaTM9q3427WDqayt49IXdCZwKFH5i8j36t2xNa9eM4jSihouf2Ge5gIKESp/EflB/dLb8PI1\nx7OtrJIrXpjPnopqryNJI6n8ReSwHNc1iReuzGX9znJ+8tI3lFfVeh1JGkHlLyKH7cTMZJ68dCBL\nC0sZ/7c8KmvqvI4kR0nlLyJH5Kw+HXwXhCnh9skLqa2r9zqSHAWVv4gcsQuPS+fB87KZuWw797+b\nj3O6KHxzowt4ishRuXpoBiX7qnnq0wKSE1rwi7N6eR1JjoDKX0SO2p1n9qSkvIqnPi2gXUIMVw/N\n8DqSHCaVv4gcNTPjd2NyKNlXzYR/Lqd9q1jO6dfR61hyGLTPX0QaJSoygj9fMpDjurTl528sYt66\nEq8jyWFQ+YtIo8VGR/LCVbl0TmrJ9a/maSroZkDlLyJNok1cDK9cM4jY6Eh+8uJ8tpdpGohgpvIX\nkSaT3jaOl64+ntL9NVzzss4CDmYqfxFpUn06JfLUZceycttebpn0rU4CC1IqfxFpcqf2as/vxuTw\n6apiHpy2TCeBBSEN9RQRv7h0cBc27irnuc/XkZEcz3Undfc6kjSg8hcRv/nlWcewcWcFD72/gm7t\n4jkjO9XrSOKj3T4i4jcREcYTFw2gb1oit01eyLKtpV5HEh+Vv4j4VcuYSF64MpfEltFc+3IeOzQE\nNCj4rfzN7DdmVmhmi3y3UQ0eu9fMCsxslZmd5a8MIhIc2reO5YWrcindX8P1f1ug6wAEAX9v+T/h\nnBvgu70PYGbZwMVAH2Ak8IyZRfo5h4h4rE+nRJ64aACLN+/hnreWaASQx7zY7TMGmOycq3LOrQcK\ngEEe5BCRABuZ04FfnNmTdxdt5dnP13odJ6z5u/xvMbMlZvaimbX1LUsDNjd4zhbfsv9gZuPNLM/M\n8oqLi/0cVUQC4eZTMxndvxP/NXMVs5Zv9zpO2GpU+ZvZLDPLP8htDPAs0AMYABQBfzrS93fOTXTO\n5TrnclNSUhoTVUSChJnx2I/6kdMpkZ+9sYiCHZoEzguNKn/n3BnOuZyD3KY657Y75+qcc/XA8/z/\nrp1CoHODt0n3LRORMBEbHclzVxxHbHQE17+6gNL9NV5HCjv+HO3T8IoO5wP5vvvTgIvNrIWZZQBZ\nwHx/5RCR4NSpTUuevfw4tuyu4LbXF1JXrwPAgeTPff6PmdlSM1sCnAr8HMA5twyYAiwHZgA3O+c0\n7kskDB3fLYnfjO7D56uL+dOHq7yOE1b8Nr2Dc+6K73nsIeAhf322iDQflw3uytItpTzz2Vr6pScy\nMkeXgQwEneErIp777Zg+9O/chjunLNYB4ABR+YuI51pERfKXy4+lZUwk4/+2gL2VOgDsbyp/EQkK\nHRNb8vSlx7KxpIK7/qEzgP1N5S8iQWNw93bce/YxzFi2jYmz13kdJ6Sp/EUkqFw7LINRfTvwhxkr\nmbu2xOs4IUvlLyJB5cAZwP3plhzPra9/y3ZNAe0XKn8RCToJLaJ47vLjKK+q49ZJC3UReD9Q+YtI\nUMpKbcUjF/Rl/oZd/PHD1V7HCTkqfxEJWmMHpnHp4C785fO1fLxCM4A2JZW/iAS1X5+bTZ9Orblj\nymI276rwOk7IUPmLSFCLjY7kmcuOpb7ecevrC6mu1f7/pqDyF5Gg17VdPH/4UT8Wbd7DHzUBXJNQ\n+YtIszCqb0euGNKVibPX8clK7f9vLJW/iDQbvzqnN9kdW3PnlMUUle73Ok6zpvIXkWYjNjqSpy87\nluraem6fvEgXgGkElb+INCsZyfH8bmwO89fv4slP1ngdp9lS+YtIs3PBselcMDCNP3+8hnnrNP/P\n0VD5i0izNGFsDl2S4vjZG4vYXV7tdZxmR+UvIs1SQosonrzkWHbuq+KetzX//5FS+YtIs9U3PZG7\nzurFzGXbmfzNZq/jNCsqfxFp1q4b1p1hmclMeG85a4v3eR2n2VD5i0izFhFh/Glcf2KjI7h9sqZ/\nOFwqfxFp9lJbx/KHC/uRX1jGnz7S9A+HQ+UvIiHhzD4duHRwFybOXsfXGv75g1T+IhIy7j+nN93a\nxXPnlMWU7q/xOk5QU/mLSMiIi4niiYsGsK2skgen5nsdJ6ip/EUkpAzo3IbbTsvi3UVbmbZ4q9dx\ngpbKX0RCzs2n9mBglzY88G4+20orvY4TlFT+IhJyoiIjeHzcAKpr67n7LZ39ezAqfxEJSRnJ8dw3\n6hhmry7mtXmbvI4TdFT+IhKyLh/SleE9U3h4+grW7yz3Ok5QUfmLSMgyMx67sB/RkcadU3Txl4ZU\n/iIS0jokxjJhTA7fbtrD81+s8zpO0FD5i0jIGzOgE2f1SeXxD1ezZvter+MEhUaVv5n92MyWmVm9\nmeV+57F7zazAzFaZ2VkNlo/0LSsws3sa8/kiIofDzHjo/L4kxEZx5z8WU1unyd8au+WfD1wAzG64\n0MyygYuBPsBI4BkzizSzSOBp4GwgG7jE91wREb9KTmjB78fmsGRLKc9+ttbrOJ5rVPk751Y45w42\nhd4YYLJzrso5tx4oAAb5bgXOuXXOuWpgsu+5IiJ+N6pvR87r34k/f7KGldvKvI7jKX/t808DGl5W\nZ4tv2aGWH5SZjTezPDPLKy4u9ktQEQkvvx3dh8SW0fziH4upCePdPz9Y/mY2y8zyD3Lz+xa7c26i\ncy7XOZebkpLi748TkTCQFB/D78bkkF9YxsTZ4Tv6J+qHnuCcO+Mo3rcQ6Nzg93TfMr5nuYhIQJzd\ntyPn9OvI/8xaw4jsVHqmtvI6UsD5a7fPNOBiM2thZhlAFjAf+AbIMrMMM4vhwEHhaX7KICJySBNG\n9yEhNoq7wnT0T2OHep5vZluAE4DpZjYTwDm3DJgCLAdmADc75+qcc7XALcBMYAUwxfdcEZGAapfQ\ngt+O7sPiLaW8OGe913ECzprLbHe5ubkuLy/P6xgiEkKcc4z/2wJmry5mxs+Gk5Ec73WkJmVmC5xz\nuQd7TGf4ikjYMjN+PzaHmKgI7nlrCfVhNPePyl9Ewlpq61juP6c389bvYtL88Jn6WeUvImFvXG5n\nhma249EPVlJUut/rOAGh8heRsGdmPHJ+P2rr63ng3fywuPKXyl9EBOjSLo47R/Ri1oodvL90m9dx\n/E7lLyLic/XQbvRNS+TBafnsqaj2Oo5fqfxFRHyiIiN49MK+7K6o4eH3V3gdx69U/iIiDfTplMj4\n4d2ZkreFuWtLvI7jNyp/EZHvuP30LLokxfGrd5ZSWVPndRy/UPmLiHxHbHQkD52fw7qd5TwTohd+\nUfmLiBzESVkpjB3QiWc/K6BgR+hd91flLyJyCPefm018iyjuezs/5KZ+UPmLiBxCckIL7ju7N/M3\n7OLNBVu8jtOkVP4iIt/jR8elM6hbEg9/sIJd5aEz9l/lLyLyPSIijN+fn8O+yloeCaGx/yp/EZEf\n0DO1FdcP784/Fmxh3rrQGPuv8hcROQy3nZZFetuW/OrdfKprm/9lH1X+IiKHoWVMJBPG9KFgxz7+\n+mXzv+yjyl9E5DCddkwqZ2an8ueP11C4p3nP+6/yFxE5Ar8+LxuACe8t8zhJ46j8RUSOQHrbOG49\nPZOZy7bz6codXsc5aip/EZEjdN2w7vRIiefBacua7cRvKn8RkSMUExXB78bmsGlXBc8204nfVP4i\nIkfhxB7JnNe/E89+vpZNJRVexzliKn8RkaP0q1G9iY4wJvyz+R38VfmLiBylDomx3HZ6FrNW7OCT\nldu9jnNEVP4iIo1w9dAMeqTE85tpy5vVwV+Vv4hII8RERTBhzIGDv8/PXud1nMOm8hcRaaShmcmM\n6tuBpz8raDZn/qr8RUSawH2jegPwcDOZ9lnlLyLSBNLbxvHTkzOZvqSIr9bu9DrOD1L5i4g0kRtO\n7k5am5b8dtpyauuCe9pnlb+ISBOJjY7kgXN7s2r7Xl77eqPXcb6Xyl9EpAmd1acDQzPb8cSsNewO\n4mv+Nqr8zezHZrbMzOrNLLfB8m5mtt/MFvluf2nw2HFmttTMCszsz2ZmjckgIhJMzIxfn9uHvZU1\nPDFrtddxDqmxW/75wAXA7IM8ttY5N8B3u7HB8meB64Es321kIzOIiASVXh1acfmQrrz29UZWbivz\nOs5BNar8nXMrnHOrDvf5ZtYRaO2c+9o554BXgbGNySAiEox+fkZPWsVGM+G95Ryou+Diz33+GWa2\n0Mw+N7OTfMvSgC0NnrPFt+ygzGy8meWZWV5xcbEfo4qINK228THcMaInX60tYeay4Jv35wfL38xm\nmVn+QW5jvudlRUAX59xA4A5gkpm1PtJwzrmJzrlc51xuSkrKkb5cRMRTlw3uQs/UBB75YAVVtcE1\n788Plr9z7gznXM5BblO/5zVVzrkS3/0FwFqgJ1AIpDd4arpvmYhIyImKjOD+c7LZWFLBK19t8DrO\nv/HLbh8zSzGzSN/97hw4sLvOOVcElJnZEN8onyuBQ/4RERFp7ob3TOG0Y9rz5McF7NxX5XWcf2ns\nUM/zzWwLcAIw3cxm+h4aDiwxs0XAm8CNzrldvsduAl4ACjjwjeCDxmQQEQl2943qzf6aOh7/KHiG\nfkY15sXOuXeAdw6y/C3grUO8Jg/Iaczniog0J5ntE7jihK688tUGrhjSld4dj/gQaJPTGb4iIgFw\n++lZtG4ZzUPTVwTF0E+Vv4hIALSJi+H207P4smAnn63yfui6yl9EJEAuG9yVjOR4fj99OTUez/qp\n8hcRCZCYqAjuPfsY1haXM3n+Jk+zqPxFRAJoRHYqgzOSeGLWGsoqazzLofIXEQkgM+OBc7PZXVHN\n058WeJZD5S8iEmA5aYmcPzCNl+ZsYMvuCk8yqPxFRDzwizN7YcAfZx72xMhNSuUvIuKBTm1acu2w\nDN5dtJWlW0oD/vkqfxERj9x4Sg+S4mN4+P3An/il8hcR8Ujr2Gh+dkYWc9eV8OmqHQH9bJW/iIiH\nLhnUhe7J8Tzy/kpqA3jil8pfRMRD0ZER3D2yF2t27OOtb7f88AuaiMpfRMRjZ/XpwLFd2vD4R6vZ\nXx2YK36p/EVEPGZm3DuqN9vLqnhxzvqAfKbKX0QkCBzfLYkR2ak8+9laSgJwxS+Vv4hIkPjlyF5U\nVNfy5Cf+n/ZB5S8iEiQy27fiouM78/d5G9lU4t9pH1T+IiJB5PbTexIZYTz+kX+nfVD5i4gEkQ6J\nsVw9NIOpi7eyfGuZ3z5H5S8iEmRuPLkHrWOjeWzmSr99hspfRCTIJLaM5qZTevDZqmLmri3xy2eo\n/EVEgtBVJ3ajY2Isj85Y6ZdJ31T+IiJBKDY6kjvP7EW/tESqapt+zp+oJn9HERFpEj86Lp0fHZfu\nl/fWlr+ISBhS+YuIhCGVv4hIGFL5i4iEIZW/iEgYUvmLiIQhlb+ISBhS+YuIhCHzx2nD/mBmxcDG\no3x5MrCzCeM0B+G4zhCe6x2O6wzhud5Hus5dnXMpB3ug2ZR/Y5hZnnMu1+scgRSO6wzhud7huM4Q\nnuvdlOus3T4iImFI5S8iEobCpfwneh3AA+G4zhCe6x2O6wzhud5Nts5hsc9fRET+Xbhs+YuISAMq\nfxGRMBTS5W9mI81slZkVmNk9XucJBDPrbGafmtlyM1tmZrd7nSlQzCzSzBaa2T+9zhIoZtbGzN40\ns5VmtsLMTvA6k7+Z2c99/23nm9nrZhbrdSZ/MLMXzWyHmeU3WJZkZh+Z2Rrfz7ZH+/4hW/5mFgk8\nDZwNZAOXmFm2t6kCoha40zmXDQwBbg6T9Qa4HVjhdYgA+x9ghnPuGKA/Ib7+ZpYG3AbkOudygEjg\nYm9T+c3LwMjvLLsH+Ng5lwV87Pv9qIRs+QODgALn3DrnXDUwGRjjcSa/c84VOee+9d3fy4EySPM2\nlf+ZWTpwDvCC11kCxcwSgdWB2XEAAAINSURBVOHAXwGcc9XOuT3epgqIKKClmUUBccBWj/P4hXNu\nNrDrO4vHAK/47r8CjD3a9w/l8k8DNjf4fQthUIINmVk3YCAwz9skAfHfwN1A01/pOnhlAMXAS77d\nXS+YWbzXofzJOVcI/BHYBBQBpc65D71NFVCpzrki3/1tQOrRvlEol39YM7ME4C3gZ865Mq/z+JOZ\nnQvscM4t8DpLgEUBxwLPOucGAuU0YjdAc+Dbxz2GA3/4OgHxZna5t6m84Q6M0z/qsfqhXP6FQOcG\nv6f7loU8M4vmQPH/3Tn3ttd5AmAoMNrMNnBg995pZvaat5ECYguwxTn3f9/s3uTAH4NQdgaw3jlX\n7JyrAd4GTvQ4UyBtN7OOAL6fO472jUK5/L8Bsswsw8xiOHBQaJrHmfzOzIwD+4BXOOce9zpPIDjn\n7nXOpTvnunHg3/MnzrmQ3xp0zm0DNptZL9+i04HlHkYKhE3AEDOL8/23fjohfpD7O6YBV/nuXwVM\nPdo3imqSOEHIOVdrZrcAMzkwIuBF59wyj2MFwlDgCmCpmS3yLbvPOfe+h5nEf24F/u7bwFkHXO1x\nHr9yzs0zszeBbzkwsm0hITrNg5m9DpwCJJvZFuBB4FFgipldy4Ep7scd9ftregcRkfATyrt9RETk\nEFT+IiJhSOUvIhKGVP4iImFI5S8iEoZU/iIiYUjlLyIShv4XmTXtHMXhnGQAAAAASUVORK5CYII=\n",
            "text/plain": [
              "<Figure size 432x288 with 1 Axes>"
            ]
          },
          "metadata": {
            "tags": []
          }
        }
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "Af6gABJfjEaV",
        "colab_type": "code",
        "colab": {}
      },
      "source": [
        "\n",
        "\n"
      ],
      "execution_count": 0,
      "outputs": []
    }
  ]
}