l02c01_celsius_to_fahrenheit.ipynb

{
  "nbformat": 4,
  "nbformat_minor": 0,
  "metadata": {
    "colab": {
      "name": "Copy of l02c01_celsius_to_fahrenheit.ipynb",
      "provenance": [],
      "collapsed_sections": [],
      "toc_visible": true,
      "include_colab_link": true
    },
    "kernelspec": {
      "display_name": "Python 3",
      "name": "python3"
    }
  },
  "cells": [
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "view-in-github",
        "colab_type": "text"
      },
      "source": [
        "<a href=\"https://colab.research.google.com/github/guowenbin90/AI-TensorFlow/blob/main/Copy_of_l02c01_celsius_to_fahrenheit.ipynb\" target=\"_parent\"><img src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open In Colab\"/></a>"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "HnKx50tv5aZD"
      },
      "source": [
        "##### Copyright 2018 The TensorFlow Authors."
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "cellView": "form",
        "id": "IwtS_OXU5cWG"
      },
      "source": [
        "#@title Licensed under the Apache License, Version 2.0 (the \"License\");\n",
        "# you may not use this file except in compliance with the License.\n",
        "# You may obtain a copy of the License at\n",
        "#\n",
        "# https://www.apache.org/licenses/LICENSE-2.0\n",
        "#\n",
        "# Unless required by applicable law or agreed to in writing, software\n",
        "# distributed under the License is distributed on an \"AS IS\" BASIS,\n",
        "# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.\n",
        "# See the License for the specific language governing permissions and\n",
        "# limitations under the License."
      ],
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "YHI3vyhv5p85"
      },
      "source": [
        "# The Basics: Training Your First Model"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "_wJ2E7jV5tN5"
      },
      "source": [
        "<table class=\"tfo-notebook-buttons\" align=\"left\">\n",
        "  <td>\n",
        "    <a target=\"_blank\" href=\"https://colab.research.google.com/github/tensorflow/examples/blob/master/courses/udacity_intro_to_tensorflow_for_deep_learning/l02c01_celsius_to_fahrenheit.ipynb\"><img src=\"https://www.tensorflow.org/images/colab_logo_32px.png\" />Run in Google Colab</a>\n",
        "  </td>\n",
        "  <td>\n",
        "    <a target=\"_blank\" href=\"https://github.com/tensorflow/examples/blob/master/courses/udacity_intro_to_tensorflow_for_deep_learning/l02c01_celsius_to_fahrenheit.ipynb\"><img src=\"https://www.tensorflow.org/images/GitHub-Mark-32px.png\" />View source on GitHub</a>\n",
        "  </td>\n",
        "</table>"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "F8YVA_634OFk"
      },
      "source": [
        "Welcome to this Colab where you will train your first Machine Learning model!\n",
        "\n",
        "We'll try to keep things simple here, and only introduce basic concepts. Later Colabs will cover more advanced problems.\n",
        "\n",
        "The problem we will solve is to convert from Celsius to Fahrenheit, where the approximate formula is:\n",
        "\n",
        "$$ f = c \\times 1.8 + 32 $$\n",
        "\n",
        "\n",
        "Of course, it would be simple enough to create a conventional Python function that directly performs this calculation, but that wouldn't be machine learning.\n",
        "\n",
        "\n",
        "Instead, we will give TensorFlow some sample Celsius values (0, 8, 15, 22, 38) and their corresponding Fahrenheit values (32, 46, 59, 72, 100).\n",
        "Then, we will train a model that figures out the above formula through the training process."
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "fA93WUy1zzWf"
      },
      "source": [
        "## Import dependencies\n",
        "\n",
        "First, import TensorFlow. Here, we're calling it `tf` for ease of use. We also tell it to only display errors.\n",
        "\n",
        "Next, import [NumPy](http://www.numpy.org/) as `np`. Numpy helps us to represent our data as highly performant lists."
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "-ZMgCvSRFqxE"
      },
      "source": [
        "import tensorflow as tf"
      ],
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "y_WQEM5MGmg3"
      },
      "source": [
        "import numpy as np\n",
        "import logging\n",
        "logger = tf.get_logger()\n",
        "logger.setLevel(logging.ERROR)"
      ],
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "AC3EQFi20buB"
      },
      "source": [
        "## Set up training data\n",
        "\n",
        "As we saw before, supervised Machine Learning is all about figuring out an algorithm given a set of inputs and outputs. Since the task in this Codelab is to create a model that can give the temperature in Fahrenheit when given the degrees in Celsius, we create two lists `celsius_q` and `fahrenheit_a` that we can use to train our model."
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "gg4pn6aI1vms",
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "outputId": "64a1024e-7941-4731-aa57-6a03544cdf19"
      },
      "source": [
        "celsius_q    = np.array([-40, -10,  0,  8, 15, 22,  38],  dtype=float)\n",
        "fahrenheit_a = np.array([-40,  14, 32, 46, 59, 72, 100],  dtype=float)\n",
        "\n",
        "for i,c in enumerate(celsius_q):\n",
        "  print(\"{} degrees Celsius = {} degrees Fahrenheit\".format(c, fahrenheit_a[i]))"
      ],
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "-40.0 degrees Celsius = -40.0 degrees Fahrenheit\n",
            "-10.0 degrees Celsius = 14.0 degrees Fahrenheit\n",
            "0.0 degrees Celsius = 32.0 degrees Fahrenheit\n",
            "8.0 degrees Celsius = 46.0 degrees Fahrenheit\n",
            "15.0 degrees Celsius = 59.0 degrees Fahrenheit\n",
            "22.0 degrees Celsius = 72.0 degrees Fahrenheit\n",
            "38.0 degrees Celsius = 100.0 degrees Fahrenheit\n"
          ]
        }
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "wwJGmDrQ0EoB"
      },
      "source": [
        "### Some Machine Learning terminology\n",
        "\n",
        " - **Feature** — The input(s) to our model. In this case, a single value — the degrees in Celsius.\n",
        "\n",
        " - **Labels** — The output our model predicts. In this case, a single value — the degrees in Fahrenheit.\n",
        "\n",
        " - **Example** — A pair of inputs/outputs used during training. In our case a pair of values from `celsius_q` and `fahrenheit_a` at a specific index, such as `(22,72)`.\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "VM7_9Klvq7MO"
      },
      "source": [
        "## Create the model\n",
        "\n",
        "Next, create the model. We will use the simplest possible model we can, a Dense network. Since the problem is straightforward, this network will require only a single layer, with a single neuron.\n",
        "\n",
        "### Build a layer\n",
        "\n",
        "We'll call the layer `l0` and create it by instantiating `tf.keras.layers.Dense` with the following configuration:\n",
        "\n",
        "*   `input_shape=[1]` — This specifies that the input to this layer is a single value. That is, the shape is a one-dimensional array with one member. Since this is the first (and only) layer, that input shape is the input shape of the entire model. The single value is a floating point number, representing degrees Celsius.\n",
        "\n",
        "*   `units=1` — This specifies the number of neurons in the layer. The number of neurons defines how many internal variables the layer has to try to learn how to solve the problem (more later). Since this is the final layer, it is also the size of the model's output — a single float value representing degrees Fahrenheit. (In a multi-layered network, the size and shape of the layer would need to match the `input_shape` of the next layer.)\n"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "pRllo2HLfXiu"
      },
      "source": [
        "l0 = tf.keras.layers.Dense(units=1, input_shape=[1])"
      ],
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "_F00_J9duLBD"
      },
      "source": [
        "### Assemble layers into the model\n",
        "\n",
        "Once layers are defined, they need to be assembled into a model. The Sequential model definition takes a list of layers as an argument, specifying the calculation order from the input to the output.\n",
        "\n",
        "This model has just a single layer, l0."
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "cSp-GpLSuMRq"
      },
      "source": [
        "model = tf.keras.Sequential([l0])"
      ],
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "t7pfHfWxust0"
      },
      "source": [
        "**Note**\n",
        "\n",
        "You will often see the layers defined inside the model definition, rather than beforehand:\n",
        "\n",
        "```python\n",
        "model = tf.keras.Sequential([\n",
        "  tf.keras.layers.Dense(units=1, input_shape=[1])\n",
        "])\n",
        "```"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "kiZG7uhm8qCF"
      },
      "source": [
        "## Compile the model, with loss and optimizer functions\n",
        "\n",
        "Before training, the model has to be compiled. When compiled for training, the model is given:\n",
        "\n",
        "- **Loss function** — A way of measuring how far off predictions are from the desired outcome. (The measured difference is called the \"loss\".)\n",
        "\n",
        "- **Optimizer function** — A way of adjusting internal values in order to reduce the loss.\n"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "m8YQN1H41L-Y"
      },
      "source": [
        "model.compile(loss='mean_squared_error',\n",
        "              optimizer=tf.keras.optimizers.Adam(0.1))"
      ],
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "17M3Pqv4P52R"
      },
      "source": [
        "These are used during training (`model.fit()`, below) to first calculate the loss at each point, and then improve it. In fact, the act of calculating the current loss of a model and then improving it is precisely what training is.\n",
        "\n",
        "During training, the optimizer function is used to calculate adjustments to the model's internal variables. The goal is to adjust the internal variables until the model (which is really a math function) mirrors the actual equation for converting Celsius to Fahrenheit.\n",
        "\n",
        "TensorFlow uses numerical analysis to perform this tuning, and all this complexity is hidden from you so we will not go into the details here. What is useful to know about these parameters are:\n",
        "\n",
        "The loss function ([mean squared error](https://en.wikipedia.org/wiki/Mean_squared_error)) and the optimizer ([Adam](https://machinelearningmastery.com/adam-optimization-algorithm-for-deep-learning/)) used here are standard for simple models like this one, but many others are available. It is not important to know how these specific functions work at this point.\n",
        "\n",
        "One part of the Optimizer you may need to think about when building your own models is the learning rate (`0.1` in the code above). This is the step size taken when adjusting values in the model. If the value is too small, it will take too many iterations to train the model. Too large, and accuracy goes down. Finding a good value often involves some trial and error, but the range is usually within 0.001 (default), and 0.1"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "c-Jk4dG91dvD"
      },
      "source": [
        "## Train the model\n",
        "\n",
        "Train the model by calling the `fit` method.\n",
        "\n",
        "During training, the model takes in Celsius values, performs a calculation using the current internal variables (called \"weights\") and outputs values which are meant to be the Fahrenheit equivalent. Since the weights are initially set randomly, the output will not be close to the correct value. The difference between the actual output and the desired output is calculated using the loss function, and the optimizer function directs how the weights should be adjusted.\n",
        "\n",
        "This cycle of calculate, compare, adjust is controlled by the `fit` method. The first argument is the inputs, the second argument is the desired outputs. The `epochs` argument specifies how many times this cycle should be run, and the `verbose` argument controls how much output the method produces."
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "lpRrl7WK10Pq",
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "outputId": "115b808e-7647-4e76-ca89-06d7ee64c485"
      },
      "source": [
        "history = model.fit(celsius_q, fahrenheit_a, epochs=500, verbose=False)\n",
        "print(\"Finished training the model\")"
      ],
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "Finished training the model\n"
          ]
        }
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "GFcIU2-SdCrI"
      },
      "source": [
        "In later videos, we will go into more detail on what actually happens here and how a Dense layer actually works internally."
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "0-QsNCLD4MJZ"
      },
      "source": [
        "## Display training statistics\n",
        "\n",
        "The `fit` method returns a history object. We can use this object to plot how the loss of our model goes down after each training epoch. A high loss means that the Fahrenheit degrees the model predicts is far from the corresponding value in `fahrenheit_a`.\n",
        "\n",
        "We'll use [Matplotlib](https://matplotlib.org/) to visualize this (you could use another tool). As you can see, our model improves very quickly at first, and then has a steady, slow improvement until it is very near \"perfect\" towards the end.\n"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "IeK6BzfbdO6_",
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 296
        },
        "outputId": "8850239f-a0fa-4849-ee93-b0fc8b15125b"
      },
      "source": [
        "import matplotlib.pyplot as plt\n",
        "plt.xlabel('Epoch Number')\n",
        "plt.ylabel(\"Loss Magnitude\")\n",
        "plt.plot(history.history['loss'])"
      ],
      "execution_count": null,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "[<matplotlib.lines.Line2D at 0x7f5c962fd150>]"
            ]
          },
          "metadata": {},
          "execution_count": 8
        },
        {
          "output_type": "display_data",
          "data": {
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAEGCAYAAACUzrmNAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3deXwX1b3/8dcne0JWIGwJO1gEZI2KaKvWVrH1grXWrSpavGirLW29bW177229bX9tbSt1X1oXtC61tlZq3QBxqaISFBc2WQQBQcIWEghZP78/ZgJf02wm+eab5f18PObxnTkz35nPYMwnc86Zc8zdERERaUxcrAMQEZGOT8lCRESapGQhIiJNUrIQEZEmKVmIiEiTEmIdQDT07t3bhwwZEuswREQ6lWXLlu1099z69nXJZDFkyBAKCwtjHYaISKdiZpsa2qdqKBERaZKShYiINEnJQkREmqRkISIiTVKyEBGRJilZiIhIk5QsRESkSUoWEUoOVjJ3wXss37w31qGIiHQoShYRqqqdGxat5Y1Ne2IdiohIh6JkEaFHcvBC+/7yqhhHIiLSsShZREhKiCMpIY7SCiULEZFIShZ1pCcnUHpQyUJEJJKSRR3pyQmqhhIRqUPJoo4eyQmUllfHOgwRkQ4lasnCzO42sx1m9m5E2W/MbLWZvW1mj5lZdsS+H5rZOjNbY2anRZRPC8vWmdk10Yq3VnpyPKXlldG+jIhIpxLNJ4t7gWl1yhYAY919HPAe8EMAMxsNnAeMCb9zq5nFm1k8cAtwOjAaOD88NmqCaig9WYiIRIpasnD3F4HddcqedffaBoFXgfxwfQbwsLuXu/v7wDrgmHBZ5+4b3L0CeDg8NmqCaii1WYiIRIplm8XXgKfC9Txgc8S+LWFZQ+X/xsxmm1mhmRUWFRW1OKiMFCULEZG6YpIszOzHQBXwQFud093vdPcCdy/Iza13Ctlm6ZGk3lAiInW1+xzcZnYJcAZwirt7WLwVGBhxWH5YRiPlUdEjOYEDFdVU1zjxcRbNS4mIdBrt+mRhZtOA7wPT3f1AxK75wHlmlmxmQ4GRwOvAUmCkmQ01sySCRvD50YwxIyUc8kNvcYuIHBK1Jwszewg4CehtZluAnxD0fkoGFpgZwKvufoW7rzCzR4CVBNVTV7p7dXieq4BngHjgbndfEa2Y4ePjQ2WmJEbzUiIinUbUkoW7n19P8V2NHP8L4Bf1lD8JPNmGoTWqNlmUHqyCrPa6qohIx6Y3uOvIqE0WauQWETlEyaKOw9VQejFPRKSWkkUd6YeeLDTkh4hILSWLOg4nCz1ZiIjUUrKoo0dyPKDZ8kREIilZ1JGeogZuEZG6lCzqSE6IJzHelCxERCIoWdSjh2bLExH5GCWLeqQnJ1CiebhFRA5RsqhHVmoi+8rUdVZEpJaSRT2yUhMpVrIQETlEyaIe2WmJ7FWyEBE5RMmiHlmpiew9oGQhIlJLyaIeWalJ7Cur5PDcTCIi3ZuSRT2yUhOpqK6hrFJDfoiIgJJFvbLTgkmP1MgtIhJQsqhHVmqQLNRuISISULKoR7aShYjIxyhZ1KNXejIAu/aXxzgSEZGOQcmiHr3SkwDYVVoR40hERDoGJYt65KQlEWewq1RPFiIioGRRr/g4o2ePJHbu15OFiAgoWTSoV49kdpboyUJEBKKYLMzsbjPbYWbvRpT1NLMFZrY2/MwJy83MbjSzdWb2tplNivjOzPD4tWY2M1rx1tUrPYlderIQEQGi+2RxLzCtTtk1wCJ3HwksCrcBTgdGhsts4DYIkgvwE+BY4BjgJ7UJJtp6pydTpCcLEREgisnC3V8EdtcpngHMC9fnAWdGlN/ngVeBbDPrD5wGLHD33e6+B1jAvyegqOiflcL2fQepqdH4UCIi7d1m0dfdt4Xr24G+4XoesDniuC1hWUPl/8bMZptZoZkVFhUVtTrQAdmpVFTVqCpKRIQYNnB7MKRrm/3Z7u53unuBuxfk5ua2+nz9s1IA2FZc1upziYh0du2dLD4Kq5cIP3eE5VuBgRHH5YdlDZVH3YDsVAA+3KtkISLS3sliPlDbo2km8HhE+cVhr6gpQHFYXfUMcKqZ5YQN26eGZVGXFyaLrXsPtsflREQ6tIRondjMHgJOAnqb2RaCXk2/Ah4xs1nAJuCc8PAngS8A64ADwKUA7r7bzH4GLA2P+z93r9toHhXZaYlkpSayvqi0PS4nItKhRS1ZuPv5Dew6pZ5jHbiygfPcDdzdhqE1i5lxRN901n5U0t6XFhHpcPQGdyOO6JvBmu0lml5VRLo9JYtGjOqXwb6DVWzerUZuEenelCwacdzw3gC8uLb1722IiHRmUWuz6AqG5/YgPyeVv7+5lU+P7M2Dr39ARVUN3zhpBLkZybEOT0Sk3ShZNMLMuOLE4fz339/lxN88T5xBjcOS9bt4/KrjSU6Ij3WIIiLtQsmiCRccM4ik+Dh27a/grEl5vLOlmMvuK+T+JZu47NPDYh2eiEi7ULJoQlyccc7Rh18i7zs6hWOG9mTeko187fihxMVZ7IITEWknauBugQunDGbz7jJeWrcz1qGIiLQLJYsWOG1MX3r2SOKh1z6IdSgiIu1CyaIFkhPi+dLEPBat/oi9BzSEuYh0fUoWLfSliXlUVjv/fGdb0weLiHRyShYtNGZAJiP7pPPYG+0yYrqISEwpWbSQmXHmxDwKN+3hg10HYh2OiEhUKVm0wowJAwB4fLmeLkSka1OyaIX8nDSOGdqTx5Zv1ci0ItKlNZkswtnrLjSz/w23B5nZMdEPrXP40sQ8NhTtp3DTnliHIiISNc15srgVOA6oncyoBLglahF1MtPHD6BnjyRuem5drEMREYma5iSLY939SuAggLvvAZKiGlUn0iM5gdmfGcaL7xWxbFO7zPgqItLumpMsKs0sHnAAM8sFaqIaVSdz0ZTB9M1M5nuPvs2BiqpYhyMi0uaakyxuBB4D+pjZL4B/Af8vqlF1Mj2SE5h7zgTe37mfH/3tHTV2i0iX0+Sos+7+gJktA04BDDjT3VdFPbJOZuqI3nz3c0fwuwXvMSA7le9PGxXrkERE2kyDycLMekZs7gAeitzn7qqgr+Oqz47gw+KD3Pr8evpmpjBz6pBYhyQi0iYae7JYRtBOYcAgYE+4ng18AAxt6UXN7DvAZeH53wEuBfoDDwO9wmtf5O4VZpYM3AdMBnYB57r7xpZeO5rMjJ/NGENRSTk//ccK+mamMG1sv1iHJSLSag22Wbj7UHcfBiwE/sPde7t7L+AM4NmWXtDM8oBvAQXuPhaIB84Dfg3MdfcRBIlpVviVWcCesHxueFyHlRAfx03nT2R8fjbfe/QtthWXxTokEZFWa04D9xR3f7J2w92fAqa28roJQKqZJQBpwDbgs8Cj4f55wJnh+oxwm3D/KWbWoaenS02K54bzJlBV7fzgr2rwFpHOrznJ4kMz+28zGxIuPwY+bOkF3X0r8FuCqqxtQDFBtdNed6/td7oFyAvX84DN4XerwuN7tfT67WVwrx5cc/ooXnyviPlvtfifS0SkQ2hOsjgfyCXoPvsY0IfDb3N/YmaWQ/C0MBQYAPQAprX0fBHnnW1mhWZWWFRU1NrTtYkLpwxmfH4WP3tiFcVllbEOR0SkxZpMFu6+293nuPvEcJnTyp5QnwPed/cid68E/gYcD2SH1VIA+UDtUK5bgYEA4f4sgobuunHe6e4F7l6Qm5vbivDaTnyc8YsvHcXu/eX85pnVsQ5HRKTFmjOQ4GIze67u0oprfgBMMbO0sO3hFGAlsBg4OzxmJvB4uD4/3Cbc/5x3okaAsXlZzJw6hAde+4A3P9BggyLSOVlTv3fNbHLEZgrwZaDK3b/f4ouaXQucC1QBbxJ0o80j6DrbMyy70N3LzSwFuB+YCOwGznP3DY2dv6CgwAsLC1saXpsrLa/ic797gZweSfzjquNJiNfI8CLS8ZjZMncvqHdfS/5IN7PX3b3DDlPe0ZIFwFPvbOPrD7zBNaeP4ooTh8c6HBGRf9NYsmhONVTPiKW3mZ1G0G4gn8C0sf2YNqYfv31mDcs094WIdDLNqQ9ZBhSGn0uAqzn8wpw0k5nx67PH0T87hasefIPd+ytiHZKISLM1J1kc6e7Dwje6R7r7qcDSaAfWFWWlJnLrBZPZVVrBnIffpLqm07TTi0g315xk8Uo9ZUvaOpDu4qj8LP5vxhheWruT6xesiXU4IiLN0tios/0IeiilmtlEgkEEATIJhuiQFjrvmEEs37yXWxavZ1x+NqeN0WCDItKxNTbq7GnAJQQvyF0fUV4C/CiKMXULP50+hpXb9nH1I28x4qp0huemxzokEZEGNTbq7Dx3Pxm4xN1Pjlimu/vf2jHGLiklMZ7bLpxMUkIcl9+/jP3lmo5VRDquBpOFmV0Yrg4xs+/WXdopvi4tLzuVm86fyIaiUn70mEanFZGOq7EG7h7hZzqQUc8ibeD4Eb35zueO4PHlH/Lg6x/EOhwRkXo12Gbh7neEn9e2Xzjd05Unj2Dppj1cO38l4/OzGZundx5FpGNpzhvcuWb2IzO708zurl3aI7juIi7O+P25E+iVnsQ3HnhDw5mLSIfTnPcsHicY3mMh8M+IRdpQzx5J3HzBJD7cW8b3/vKW2i9EpENprOtsrTR3/0HUIxEmD87hmtNH8fN/ruK2F9bzjZNGxDokERGgeU8WT5jZF6IeiQAw64ShTB8/gN88s4ZnV2yPdTgiIkDzksUcgoRRZmb7zKzEzPZFO7Duysy47uxxjMvL4tt/Xs6qbfqnFpHYa860qhnuHufuqe6eGW5ntkdw3VVKYjx3XlxAZkoil80rZGdpeaxDEpFurjm9oSbVswyPmC9boqBvZgp/uLiAXfvLueL+ZZRXVcc6JBHpxppTDXUr8Crwh3B5FfgLsMbMTo1ibN3eUflZ/PYr4ynctIcfP/auekiJSMw0J1l8CEx098nuPhmYAGwAPg9cF83gBM4YN4A5p4zk0WVb+MNLjU49LiISNc2pSjrC3VfUbrj7SjMb5e4bzKyx70kbmXPKSNbtKOWXT61meG46pxzZN9YhiUg305wnixVmdpuZnRgutwIrzSwZ0KvG7SAuzvjtV8YzZkAm33roTdZsL4l1SCLSzTQnWVwCrAO+HS4bwrJK4ORoBSYfl5oUzx8uLiAtOYHL7lvKLvWQEpF21Jyus2Xu/jt3/1K4/NbdD7h7jbuXtkeQEuiflcofLi5gx75yZs0rpKxCPaREpH00p+vsSDN71MxWmtmG2qU1FzWz7PCcq81slZkdZ2Y9zWyBma0NP3PCY83MbjSzdWb2tplNas21O7sJA7O58fyJvLVlL9986A2qqmtiHZKIdAPNqYa6B7gNqCKodroP+FMrr3sD8LS7jwLGA6uAa4BF7j4SWBRuA5wOjAyX2WEs3dppY/px7fQxLFy1g5/+Y4W61IpI1DUnWaS6+yLA3H2Tu/8U+GJLL2hmWcBngLsA3L3C3fcCM4B54WHzgDPD9RnAfR54Fcg2s/4tvX5XcfFxQ7jixOH86dUPuPX59bEOR0S6uOZ0nS03szhgrZldBWwlmD2vpYYCRcA9ZjYeWEYw/lRfd98WHrMdqO0fmgdsjvj+lrBsW0QZZjab4MmDQYMGtSK8zuP7p32K7cVl/OaZNfTPSuGsSfmxDklEuqjmDiSYBnwLmAxcBMxsxTUTgEnAbe4+EdjP4SonADyoV/lEdSvufqe7F7h7QW5ubivC6zzi4ozrzh7P1OG9+P6jb/OvtTtjHZKIdFHN6Q211N1L3X2Lu1/q7meF1UEttQXY4u6vhduPEiSPj2qrl8LPHeH+rcDAiO/nh2UCJCXEcftFkxnRJ50r/rSMFR8WxzokEemCGkwWZja/saWlF3T37cBmM/tUWHQKsBKYz+EnlpkEM/QRll8c9oqaAhRHVFcJkJmSyD2XHk1mSgIz736dDUXq0Swibcsa6kljZkUEbQUPAa8BHxvbw91faPFFzSYAfwSSCF7yu5QgcT0CDAI2Aee4+24LxhS5GZgGHAAudffCxs5fUFDghYWNHtIlrS8q5Zzbl5CSGM9frjiOAdmpsQ5JRDoRM1vm7gX17mskWcQTDBZ4PjCOYN7thyLHieqoumuyAHh3azHn3/kqfTKTeeTy4+iVnhzrkESkk2gsWTRYDeXu1e7+tLvPBKYQDPnxfNgjSjqosXlZ3HXJ0WzZU8Yl9yyl5KCG7xKR1mu0gdvMks3sLIKX8K4EbgQea4/ApOWOGdqT2y+czKpt+5g1r5CDlRoWRERap7EG7vuAJQQ9la5196Pd/Wfurp5IncDJo/pw/bkTWLpxN1//0zIqqjQsiIi0XGNPFhcSDLExB3jFzPaFS4mZ7Wuf8KQ1po8fwC/OPIrFa4q46sE3qNQ4UiLSQo21WcS5e0a4ZEYsGe6e2Z5BSstdcOwgrp0+hmdXfsSch9/UwIMi0iLNGe5DOrmZU4dQWV3Dz/+5ivi4t5h7zngS4pvz8r6ISEDJopu47NPDqK5xfvnUahLCmffi4zQtrog0j5JFN3L5icOpqnF+88wa4uOM6748jjglDBFphiaThZn1AMrcvcbMjgBGAU+5uzrwd0JXnjyCyuoafr9wLe5w3dnj9IQhIk1qzpPFi8Cnw5nrngWWAucCX41mYBI9c04ZiWHMXfge5VXVzD13AolqwxCRRjQnWZi7HzCzWcCt7n6dmS2PdmASPWbGnM+NJCUxjl8+tZqKqhpuumAiyQnxsQ5NRDqo5vw5aWZ2HMGTxD/DMv1W6QIuP3H4oW61s+9bpje9RaRBzUkW3wZ+CDzm7ivMbBiwOLphSXuZOXUIvzrrKF5cW8Sl9yxlf3lVrEMSkQ6oOZMfveDu09391+H0qjvd/VvtEJu0k/OOGcT154zntfd3MfPu19mnwQdFpI4mk4WZPWhmmWGvqHeBlWb2veiHJu3pSxPzufmCSSzfvJfz7niVHSUHYx2SiHQgzamGGu3u+4AzgaeAoQTzcEsX84Wj+vPHmQW8v3M/Z9+2hE279sc6JBHpIJqTLBLNLJEgWcwP36+of8Yk6fRO+lQfHvzPYyk5WMmXb3uFd7dqTm8RaV6yuAPYCPQAXjSzwYBGne3CJg7K4S9XTCUpPo7z7nyVV9bvjHVIIhJjzWngvtHd89z9Cx7YBJzcDrFJDI3ok85fvzGVAdkpXHL3Up58Z1usQxKRGGpOA3eWmV1vZoXh8juCpwzp4vpnpfKXy6cyLj+LKx98g3mvbIx1SCISI82phrobKAHOCZd9wD3RDEo6jqy0RO6fdSynjOrLT+av4KfzV1BdoyYrke6mOcliuLv/xN03hMu1wLBoByYdR2pSPHdcNJlZJwzl3lc2Mvu+Qr28J9LNNCdZlJnZCbUbZnY8UBa9kKQjio8z/ueM0fzszLEsXrODr9y+hG3F+jEQ6S6akyyuAG4xs41mthG4Gbi8tRc2s3gze9PMngi3h5rZa2a2zsz+bGZJYXlyuL0u3D+ktdeWlrtoymDuvuRoPth9gDNveVlda0W6ieb0hnrL3ccD44Bx7j4R+GwbXHsOsCpi+9fAXHcfAewBZoXls4A9Yfnc8DiJoZM+1YdHv34c8WZ85fYlPLtie6xDEpEoa/YkBu6+L3yTG+C7rbmomeUDXwT+GG4bQQJ6NDxkHsFLgAAzwm3C/aeEx0sMjeqXyd+vPJ6RfdOZff8ybli4lho1fIt0WS2d8aa1v6x/D3wfqAm3ewF73b221XQLkBeu5wGbAcL9xeHxHw/IbHZt996ioqJWhifN0SczhUcuP46zJuYxd+F7XPGnZZRoEEKRLqmlyaLFf0Ka2RnADndf1tJz1Mfd73T3AncvyM3NbctTSyNSEuP53Tnj+d8zRrNo9Q6+dOsrbCgqjXVYItLGGkwWZlZiZvvqWUqAAa245vHA9LCx/GGC6qcbgGwzq525Lx/YGq5vBQaGMSUAWcCuVlxf2piZ8bUThnL/rGPYVVrOjFteZvHqHbEOS0TaUIPJwt0z3D2zniXD3ZszHWtD5/2hu+e7+xDgPOA5d/8qwYRKZ4eHzQQeD9fnh9uE+59zd1WOd0BTh/dm/lUnMDAnja/NW8pNi9SOIdJVtLQaKhp+AHzXzNYRtEncFZbfBfQKy78LXBOj+KQZBvZM469fn8r08QP43YL3uPTepezeXxHrsESklawr/pFeUFDghYWFsQ6jW3N3Hnz9A679x0p6piVx8wUTKRjSM9ZhiUgjzGyZuxfUt68jPVlIF2JmfPXYwfzt61NJTozj3Dtf5fYX1qtaSqSTUrKQqBqbl8U/vnkCp43py6+eWs1l9xWyR9VSIp2OkoVEXWZKIrdcMIlrp4/hpbVFfOHGlzShkkgno2Qh7cLMmDl1CH/9+lRSE+P56h9f45dPraKiqqbpL4tIzClZSLsal5/NE986gfOOHsQdL2zgS7e+zLodJbEOS0SaoGQh7S4tKYFfnnUUd140mW3FB/nijf/i/iUb6Yo980S6CiULiZlTx/Tj6W9/minDevE/j6/g0nuXsr34YKzDEpF6KFlITPXJSOHeS4/mp/8xmlc37OLzc1/gkcLNesoQ6WCULCTmzIxLjh/K03M+w5H9M/n+o28z856lbN2rmfhEOgolC+kwhvTuwcP/OYVrp4+hcONuTpv7Ig++9oGeMkQ6ACUL6VDi4oIuts98+zMclZfFjx57hwvveo2NO/fHOjSRbk3JQjqkgT3TeOCyY/n5mWN5e3Mxp/7+RW5YuJbyqupYhybSLSlZSIcVF2dcOGUwC68+kVNH92Xuwvc4/fcv8co6vf0t0t6ULKTD65uZws0XTGLe146hqsa54I+v8Z0/L6eopDzWoYl0G0oW0mmceEQuz37nM3zzsyN44u0POeV3z3PPy+9TWa0hQ0SiTclCOpWUxHiuPvVTPDXnM4zLz+baf6zk9Bte4oX3imIdmkiXpmQhndKIPuncP+sY/nBxAVXVNcy8+3W+du9SNhSVxjo0kS5JyUI6LTPj86P78sx3PsMPTx/F6+/v5tS5L/LzJ1ZSfKAy1uGJdClKFtLpJSfEc/mJw1n8Xyfx5Un53PXy+3z6uue4/YX1HKxUV1uRtqBkIV1GbkYyvz57HE9+69NMHpzDr55azUm/eZ4/L/2AKjWCi7SKkoV0OUf2z+SeS4/h4dlT6JeVwg/++g7TbniJZ1Zs19AhIi2kZCFd1pRhvXjsG1O5/cJJ1Lhz+f3LmHHLyyxa9ZGShsgn1O7JwswGmtliM1tpZivMbE5Y3tPMFpjZ2vAzJyw3M7vRzNaZ2dtmNqm9Y5bOy8yYNrY/z377M1z35XHsOVDBrHmFTL/5ZRauVNIQaa5YPFlUAVe7+2hgCnClmY0GrgEWuftIYFG4DXA6MDJcZgO3tX/I0tklxMdxztEDee7qk7ju7HEUl1Vy2X2F/MfN/+JZVU+JNKndk4W7b3P3N8L1EmAVkAfMAOaFh80DzgzXZwD3eeBVINvM+rdz2NJFJMbHcU7BQBZdfSK/OXscJQermH3/Mr5447+Y/9aHaggXaUBM2yzMbAgwEXgN6Ovu28Jd24G+4XoesDnia1vCsrrnmm1mhWZWWFSkt3mlcYnxcXylYCCLvnsiv/3KeA5WVfOth97kpN8GQ4gcqKiKdYgiHUrMkoWZpQN/Bb7t7vsi93lQJ/CJ6gXc/U53L3D3gtzc3DaMVLqyhPg4zp6cz8LvnMidF02mX2YK1/5jJVN/9RzXP7uGnaUarFAEICEWFzWzRIJE8YC7/y0s/sjM+rv7trCaaUdYvhUYGPH1/LBMpM3ExRmnjunHqWP6sWzTbu54YQM3LV7HHS9u4KxJecycOoRR/TJjHaZIzMSiN5QBdwGr3P36iF3zgZnh+kzg8Yjyi8NeUVOA4ojqKpE2N3lwT+68uICF3z2Rsybl8bc3tjLt9y9x7h1LePKdbRrlVrola+9eIGZ2AvAS8A5Q+3/djwjaLR4BBgGbgHPcfXeYXG4GpgEHgEvdvbCxaxQUFHhhYaOHiDTb3gMVPFK4mfuWbGLLnjL6Zabw1WMHcd4xg8jNSI51eCJtxsyWuXtBvfu6YpdBJQuJhuoa5/k1O7j3lY28tHYnifFB1dW5BQM5YURv4uIs1iGKtEpjySImbRYinVF8nHHKkX055ci+rC8q5U+vbuKxN7fyz7e3kZedytmT8/lKQT75OWmxDlWkzenJQqQVDlZWs2DlRzxSuJl/hXODnzCiN+cePZDPHdmXlMT4GEco0nyqhhJpB5t3H+DRZVt4dNkWtu4tIyM5gdPG9mPGhAFMHd6beFVTSQenZCHSjqprnCXrd/H48q08/e52SsqryM1I5oxx/TlzQh7j8rMI+m2IdCxKFiIxcrCymsWrd/D35VtZvLqIiuoahvbuwReP6s+0sf0YMyBTiUM6DCULkQ6g+EAlT6/YxuPLP+TVDbuoccjPSWXamH5MG9uPSYNy1KNKYkrJQqSD2b2/goUrP+LpFdv519qdVFTXkJuRzKmj+/L50X2ZMqyXGsel3SlZiHRgJQcrWbymiKff3cbi1UWUVVaTmhjP8SN6cdKn+nDyqD7kZafGOkzpBvSehUgHlpGSyPTxA5g+fgAHK6tZsn4Xi9fs4LnVO1i4Khgi7Yi+6Zw8qg8nHdGHSYOzSU7QU4e0Lz1ZiHRQ7s76olIWry5i8ZodvP7+bqpqnJTEOI4e0pPjR/Tm+OG9GT0gU91ypU2oGkqkCyg5WMkr63exZP0uXl63k7U7SgHISk1kyrAgeUwd3ovhuenqYSUtomookS4gIyWR08b047Qx/QDYse8gSzYEiePldbt4ZsVHAOSkJTJ5cE8KhuRw9JAcxuZlqdpKWk3JQqST6pOZwowJecyYkIe7s3l3GUs27KRw4x4KN+1h4aogeSQlxDE+PytIIINzGDcwiz4ZKTGOXjobVUOJdFE7S8sp3LiHZZt2s3TjHlZ8WExldfD/e/+sFI7Ky2L8wGzG5WcxLi+brLTEGEcssaZqKJFuqHd6MtPGBi/8QfA2+Ttbi3lr817e2VrM21uKeXblR4eOH9wrjXH52YwZkMmofhkc2T+TPhnJav8QQMlCpPrx4XgAAAruSURBVNtISYzn6CE9OXpIz0NlxQcqeffDYt7aspe3NxfzxqY9/OOtDw/tz0lLZFS/TEb1z+DI8POIvhl6YbAbUrIQ6cay0hKDLrgjeh8qKz5Qyert+1i9vYTV2/exalsJD7++mbLKagDiDAb1TGN4bjrD+6QzrHcPhvdJZ3huOj17JMXqViTKlCxE5GOy0hI5dlgvjh3W61BZTY3zwe4DrN6+j5XbSli/o5T1RaW8tG4nFVWH5yTPSUtkWG46w3N7MDw3ncG9ejCwZyqDeqaRkaI2kc5MDdwi0mLVNc6He8tYV1QaJpD9bCgKPneWln/s2Jy0RAb1TGNguAyKWPpnpZAQHxeju5BaauAWkaiIj7NDv/xP/lSfj+0rLqtk8+4DfBAutevvbi3m6Xe3U1Vz+A/VOIPcjGT6ZaXSPzOF/tkp9M9KCbazUuiXmULfzBSSEpRQYkXJQkSiIis1kay8LMbmZf3bvuoaZ1txGZt3l/HB7v1s3XuQ7cVlbCs+yLqiUl5aW8T+iuqPfccs6OHVNzOZ3umHl9yMZHqnJ5GbnkzvjGRy05PJSk3UcO9tTMlCRNpdfJyRn5NGfk4axw3vVe8xJQcr2V58kG3FB9kWJpJtew+yo+QgO0srWLO9hJ2l5YfeHYmUEGf0Sk+id3oyPXskkZ2WRHZqIjlpicF6WiI54Wd2WhI5aYlkpijBNEbJQkQ6pIyURDJSEhnZN6PBY9yd4rJKdpaWs6OknJ2lFewsKWdnaTlF4eeeA0F12N6ySorLKmmomdYseBrKSUsiIyWB9OQEMlISyEhJJD05gcyUBNLD7cP7Ez9WnpoY32UHdew0ycLMpgE3APHAH939VzEOSURizMzCJ4UkRvRpOKnUqq5x9pVVsudABXvLKtl7oII9+ysPrx+oYO+BSkoOVlFaXsXOnfspPVgVbFdUNZhoIiUnxJGWFE9qYjypSfGkJSWEn/FhecLh9UOfCaQmxpOcEBcs4XpS7XZCuC/x8HpSfFy7Pgl1imRhZvHALcDngS3AUjOb7+4rYxuZiHQm8XFGTo8kclrwPkhNjbO/IkgiJQdrl8qI7UoOVFRTVllNWUV1sF5RzYGKKg5UVLN7fwVb9ny8rDyi23FLJMXHHUoiSfFBkhmbl8VN509s1Xnr0ymSBXAMsM7dNwCY2cPADEDJQkTaRVycHaoa6//vbfYtUl3jlFUGyaOsopqKqhrKDy1BMimvrKGiuobyyupD+yr+bX815ZXBvoE9ozOrYmdJFnnA5ojtLcCxkQeY2WxgNsCgQYPaLzIRkRaKjzPSk4P2j46uy3Radvc73b3A3Qtyc3NjHY6ISJfSWZLFVmBgxHZ+WCYiIu2gsySLpcBIMxtqZknAecD8GMckItJtdPyKMsDdq8zsKuAZgq6zd7v7ihiHJSLSbXSKZAHg7k8CT8Y6DhGR7qizVEOJiEgMKVmIiEiTlCxERKRJXXLyIzMrAja14hS9gZ1tFE5noXvuHnTP3UNL73mwu9f7olqXTBatZWaFDc0W1VXpnrsH3XP3EI17VjWUiIg0SclCRESapGRRvztjHUAM6J67B91z99Dm96w2CxERaZKeLEREpElKFiIi0iQliwhmNs3M1pjZOjO7JtbxtBUzu9vMdpjZuxFlPc1sgZmtDT9zwnIzsxvDf4O3zWxS7CJvOTMbaGaLzWylma0wszlheZe9bzNLMbPXzeyt8J6vDcuHmtlr4b39ORy5GTNLDrfXhfuHxDL+1jCzeDN708yeCLe79D2b2UYze8fMlptZYVgW1Z9tJYtQxDzfpwOjgfPNbHRso2oz9wLT6pRdAyxy95HAonAbgvsfGS6zgdvaKca2VgVc7e6jgSnAleF/z6583+XAZ919PDABmGZmU4BfA3PdfQSwB5gVHj8L2BOWzw2P66zmAKsitrvDPZ/s7hMi3qeI7s+2u2sJGvmPA56J2P4h8MNYx9WG9zcEeDdiew3QP1zvD6wJ1+8Azq/vuM68AI8Dn+8u9w2kAW8QTD+8E0gIyw/9nBMM+X9cuJ4QHmexjr0F95of/nL8LPAEYN3gnjcCveuURfVnW08Wh9U3z3dejGJpD33dfVu4vh3oG653uX+HsKphIvAaXfy+w+qY5cAOYAGwHtjr7lXhIZH3deiew/3FQK/2jbhN/B74PlATbvei69+zA8+a2TIzmx2WRfVnu9PMZyHR4+5uZl2yD7WZpQN/Bb7t7vvM7NC+rnjf7l4NTDCzbOAxYFSMQ4oqMzsD2OHuy8zspFjH045OcPetZtYHWGBmqyN3RuNnW08Wh3W3eb4/MrP+AOHnjrC8y/w7mFkiQaJ4wN3/FhZ3+fsGcPe9wGKCKphsM6v9wzDyvg7dc7g/C9jVzqG21vHAdDPbCDxMUBV1A137nnH3reHnDoI/Co4hyj/bShaHdbd5vucDM8P1mQR1+rXlF4c9KKYAxRGPtp2GBY8QdwGr3P36iF1d9r7NLDd8osDMUgnaaFYRJI2zw8Pq3nPtv8XZwHMeVmp3Fu7+Q3fPd/chBP/PPufuX6UL37OZ9TCzjNp14FTgXaL9sx3rhpqOtABfAN4jqOf9cazjacP7egjYBlQS1FfOIqinXQSsBRYCPcNjjaBX2HrgHaAg1vG38J5PIKjXfRtYHi5f6Mr3DYwD3gzv+V3gf8PyYcDrwDrgL0ByWJ4Sbq8L9w+L9T208v5PAp7o6vcc3ttb4bKi9ndVtH+2NdyHiIg0SdVQIiLSJCULERFpkpKFiIg0SclCRESapGQhIiJNUrKQLs/MqsPROWuXNhtR2MyGWMRovo0c91MzOxC+cVtbVtqeMYi0hob7kO6gzN0nxDoIgkHrrgZ+EOtAIplZgh8eR0mkXnqykG4rnBPgunBegNfNbERYPsTMngvH/l9kZoPC8r5m9lg4X8RbZjY1PFW8mf3Bgjkkng3fnq7P3cC5ZtazThwfezIws/8ys5+G68+b2VwzKzSzVWZ2tJn9LZyz4OcRp0kwswfCYx41s7Tw+5PN7IVwwLlnIoaDeN7Mfm/BXAhzWv+vKV2dkoV0B6l1qqHOjdhX7O5HATcTjF4KcBMwz93HAQ8AN4blNwIveDBfxCSCt2chmCfgFncfA+wFvtxAHKUECeOT/nKu8GDOgtsJhnC4EhgLXGJmtSOmfgq41d2PBPYB3wjHxroJONvdJ4fX/kXEeZPcvcDdf/cJ45FuSNVQ0h00Vg31UMTn3HD9OOCscP1+4Lpw/bPAxXBodNdiC2Yje9/dl4fHLCOYO6QhNwLLzey3nyD+2jHK3gFWeDiuj5ltIBggbi+w2d1fDo/7E/At4GmCpLIgHG03nmDYl1p//gQxSDenZCHdnTew/kmUR6xXAw1VQ+Hue83sQYKng1pVfPwpP6WB89fUuVYNh/8frhu7E4wJtMLdj2sgnP0NxSlSl6qhpLs7N+JzSbj+CsEIpgBfBV4K1xcBX4dDkwxltfCa1wOXc/gX/UdAHzPrZWbJwBktOOcgM6tNChcA/yKYES23ttzMEs1sTAtjlm5OyUK6g7ptFr+K2JdjZm8TtCN8Jyz7JnBpWH4Rh9sY5gAnm9k7BNVNLZqj3d13EsxBkBxuVwL/RzAK6gJgdcPfbtAagnnGVwE5wG3uXkEwDPevzewtgpF3pzZyDpEGadRZ6bbCCXMKwl/eItIIPVmIiEiT9GQhIiJN0pOFiIg0SclCRESapGQhIiJNUrIQEZEmKVmIiEiT/j/Vfw0EtEsDIQAAAABJRU5ErkJggg==\n",
            "text/plain": [
              "<Figure size 432x288 with 1 Axes>"
            ]
          },
          "metadata": {
            "needs_background": "light"
          }
        }
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "LtQGDMob5LOD"
      },
      "source": [
        "## Use the model to predict values\n",
        "\n",
        "Now you have a model that has been trained to learn the relationship between `celsius_q` and `fahrenheit_a`. You can use the predict method to have it calculate the Fahrenheit degrees for a previously unknown Celsius degrees.\n",
        "\n",
        "So, for example, if the Celsius value is 100, what do you think the Fahrenheit result will be? Take a guess before you run this code."
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "oxNzL4lS2Gui",
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "outputId": "8e69a6df-9fde-43ea-c243-c8724c3cd818"
      },
      "source": [
        "print(model.predict([100.0]))"
      ],
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "[[211.33829]]\n"
          ]
        }
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "jApk6tZ1fBg1"
      },
      "source": [
        "The correct answer is $100 \\times 1.8 + 32 = 212$, so our model is doing really well.\n",
        "\n",
        "### To review\n",
        "\n",
        "\n",
        "*   We created a model with a Dense layer\n",
        "*   We trained it with 3500 examples (7 pairs, over 500 epochs).\n",
        "\n",
        "Our model tuned the variables (weights) in the Dense layer until it was able to return the correct Fahrenheit value for any Celsius value. (Remember, 100 Celsius was not part of our training data.)\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "zRrOky5gm20Z"
      },
      "source": [
        "## Looking at the layer weights\n",
        "\n",
        "Finally, let's print the internal variables of the Dense layer. "
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "kmIkVdkbnZJI",
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "outputId": "9100c911-2cff-4beb-e66e-ba861d67e48c"
      },
      "source": [
        "print(\"These are the layer variables: {}\".format(l0.get_weights()))"
      ],
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "These are the layer variables: [array([[1.8203623]], dtype=float32), array([29.302042], dtype=float32)]\n"
          ]
        }
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "RSplSnMvnWC-"
      },
      "source": [
        "The first variable is close to ~1.8 and the second to ~32. These values (1.8 and 32) are the actual variables in the real conversion formula.\n",
        "\n",
        "This is really close to the values in the conversion formula. We'll explain this in an upcoming video where we show how a Dense layer works, but for a single neuron with a single input and a single output, the internal math looks the same as [the equation for a line](https://en.wikipedia.org/wiki/Linear_equation#Slope%E2%80%93intercept_form), $y = mx + b$, which has the same form as the conversion equation, $f = 1.8c + 32$.\n",
        "\n",
        "Since the form is the same, the variables should converge on the standard values of 1.8 and 32, which is exactly what happened.\n",
        "\n",
        "With additional neurons, additional inputs, and additional outputs, the formula becomes much more complex, but the idea is the same.\n",
        "\n",
        "### A little experiment\n",
        "\n",
        "Just for fun, what if we created more Dense layers with different units, which therefore also has more variables?"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "Y2zTA-rDS5Xk",
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "outputId": "2cd961f2-05a8-44e2-9e43-f79c3e319f48"
      },
      "source": [
        "l0 = tf.keras.layers.Dense(units=4, input_shape=[1])\n",
        "l1 = tf.keras.layers.Dense(units=4)\n",
        "l2 = tf.keras.layers.Dense(units=1)\n",
        "model = tf.keras.Sequential([l0, l1, l2])\n",
        "model.compile(loss='mean_squared_error', optimizer=tf.keras.optimizers.Adam(0.1))\n",
        "model.fit(celsius_q, fahrenheit_a, epochs=500, verbose=False)\n",
        "print(\"Finished training the model\")\n",
        "print(model.predict([100.0]))\n",
        "print(\"Model predicts that 100 degrees Celsius is: {} degrees Fahrenheit\".format(model.predict([100.0])))\n",
        "print(\"These are the l0 variables: {}\".format(l0.get_weights()))\n",
        "print(\"These are the l1 variables: {}\".format(l1.get_weights()))\n",
        "print(\"These are the l2 variables: {}\".format(l2.get_weights()))"
      ],
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "Finished training the model\n",
            "[[211.74744]]\n",
            "Model predicts that 100 degrees Celsius is: [[211.74744]] degrees Fahrenheit\n",
            "These are the l0 variables: [array([[-0.78645104, -0.22119656,  0.21469362,  0.0202496 ]],\n",
            "      dtype=float32), array([-3.6797    , -3.7694812 , -0.27469736,  3.3668094 ], dtype=float32)]\n",
            "These are the l1 variables: [array([[-0.3104039 ,  0.1838275 ,  0.65251106,  1.0428878 ],\n",
            "       [-1.1052334 , -0.4231577 , -0.02089227,  0.4438429 ],\n",
            "       [ 0.29398388, -0.18603238,  0.43652624,  0.4234833 ],\n",
            "       [ 1.1874797 , -0.81879693,  0.15731238,  0.04826059]],\n",
            "      dtype=float32), array([ 3.6194685, -1.7177329, -1.6768664, -3.5703118], dtype=float32)]\n",
            "These are the l2 variables: [array([[ 1.2847482 ],\n",
            "       [-0.2723774 ],\n",
            "       [-0.33395264],\n",
            "       [-1.078701  ]], dtype=float32), array([3.5507586], dtype=float32)]\n"
          ]
        }
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "xrpFFlgYhCty"
      },
      "source": [
        "As you can see, this model is also able to predict the corresponding Fahrenheit value really well. But when you look at the variables (weights) in the `l0` and `l1` layers, they are nothing even close to ~1.8 and ~32. The added complexity hides the \"simple\" form of the conversion equation.\n",
        "\n",
        "Stay tuned for the upcoming video on how Dense layers work for the explanation."
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "VM0UzFUVsfsz",
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "outputId": "62d2d5ce-b08d-41a1-aebf-4bbd935c0b7f"
      },
      "source": [
        "l0 =tf.keras.layers.Dense(units=1,input_shape=[1])\n",
        "model = tf.keras.Sequential([l0])\n",
        "model.compile(loss='mean_squared_error',optimizer=tf.keras.optimizers.Adam(.1))\n",
        "history = model.fit(celsius_q,fahrenheit_a,epochs=500,verbose=False)\n",
        "model.predict([100.0])"
      ],
      "execution_count": null,
      "outputs": [
        {
          "output_type": "execute_result",
          "data": {
            "text/plain": [
              "array([[211.29646]], dtype=float32)"
            ]
          },
          "metadata": {},
          "execution_count": 12
        }
      ]
    }
  ]
}