3_Hermite_and_Newtons_Divided_Difference_Interpolation.ipynb

{
  "cells": [
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "view-in-github",
        "colab_type": "text"
      },
      "source": [
        "<a href=\"https://colab.research.google.com/github/Mouly22/Numerical-Methods/blob/main/3_Hermite_and_Newtons_Divided_Difference_Interpolation.ipynb\" target=\"_parent\"><img src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open In Colab\"/></a>"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "NQ8skUBePlvZ"
      },
      "source": [
        "---"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "mUa8n1tfPlva",
        "nbgrader": {
          "cell_type": "markdown",
          "checksum": "9b7261b005274f20aab5b5157815cc4f",
          "grade": false,
          "grade_id": "title",
          "locked": true,
          "schema_version": 3,
          "solution": false,
          "task": false
        }
      },
      "source": [
        "# Part 1: Hermite Interpolation\n",
        "---\n",
        "Hermite Interpolation is an example of a variant of the interpolation problem, where the interpolant matches one or more **derivatives of $f$**  at each of the nodes, in addition to the function values."
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "t6HuHqs8Plva",
        "nbgrader": {
          "cell_type": "markdown",
          "checksum": "2ac7b29da925009a6e6e93ce16cd1871",
          "grade": false,
          "grade_id": "intro1",
          "locked": true,
          "schema_version": 3,
          "solution": false,
          "task": false
        }
      },
      "source": [
        "## Importing the necessary libraries"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "zaOqCWsePlvb",
        "nbgrader": {
          "cell_type": "code",
          "checksum": "b3ce57ac45d5516f69d56b4a7c02a1d6",
          "grade": false,
          "grade_id": "import",
          "locked": true,
          "schema_version": 3,
          "solution": false,
          "task": false
        }
      },
      "outputs": [],
      "source": [
        "import numpy as np\n",
        "import matplotlib.pyplot as plt\n",
        "from itertools import combinations\n",
        "from numpy.polynomial import Polynomial"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "wR9qR23NPlvb",
        "nbgrader": {
          "cell_type": "markdown",
          "checksum": "9368014a5a45cd31c846853230dae411",
          "grade": false,
          "grade_id": "intro2",
          "locked": true,
          "schema_version": 3,
          "solution": false,
          "task": false
        }
      },
      "source": [
        "## Creating the components for Hermite interpolation\n",
        "\n",
        "For the case of Hermite Interpolation, we look for a polynomial that matches both $f'(x_i)$ and $f(x_i)$ at the nodes $x_i = x_0,\\dots,x_n$. Say you have $n+1$ data points, $(x_0, y_0), (x_1, y_1), x_2, y_2), \\dots, (x_n, y_n)$ and you happen to know the first-order derivative at all of these points, namely, $(x_0, y_0 ^\\prime ), (x_1, y_1 ^\\prime ), x_2, y_2 ^\\prime ), \\dots ,(x_n, y_n ^\\prime )$. According to hermite interpolation, since there are $2n + 2$ conditions; $n+1$ for $f(x_i)$ plus $n+1$ for $f'(x_i)$; you can fit a polynomial of order $2n+1$. \n",
        "\n",
        "General form of a $2n+1$ degree Hermite polynomial:\n",
        "\n",
        "$$p_{2n+1} = \\sum_{k=0}^{n} \\left(f(x_k)h_k(x) + f'(x_k)\\hat{h}_k(x)\\right), \\tag{1}$$\n",
        "\n",
        "where $h_k$ and $\\hat{h}_k$ are defined using Lagrange basis functions by the following equations:\n",
        "\n",
        "$$h_k(x) = (1-2(x-x_k)l^\\prime_k(x_k))l^2_k(x_k), \\tag{2}$$\n",
        "\n",
        "and\n",
        "\n",
        "$$\\hat{h}_k(x) = (x-x_k)l^2_k(x_k), \\tag{3}$$\n",
        "\n",
        "where the Lagrange basis function being:\n",
        "\n",
        "$$l_k(x) = \\prod_{j=0, j\\neq k}^{n} \\frac{x-x_j}{x_k-x_j}. \\tag{4}$$\n",
        "\n",
        "**Note** that, we can rewrite Equation $(2)$ in this way,\n",
        "\n",
        "\\begin{align}\n",
        "h_k(x) &= \\left(1-2(x-x_k)l^\\prime_k(x_k) \\right)l^2_k(x_k) \\\\\n",
        "&= \\left(1 - 2xl^\\prime_k(x_k) + 2x_kl^\\prime_k(x_k) \\right)l^2_k(x_k) \\\\\n",
        "&= \\left(1 + 2x_kl^\\prime_k(x_k) - 2l'_k(x_k)x \\right) l^2_k(x_k) \\tag{5}\n",
        "\\end{align}\n",
        "Replacing $l^\\prime_k(x_k)$ with $m$, we get:\n",
        "$$h_k(x) = (1 - 2xm + 2x_km)l^2_k(x_k). \\tag{6}$$\n",
        "\n",
        "# Tasks:\n",
        "\n",
        "* The functions: `l(k, x)`, `h(k, x)` and `h_hat(k, x)` calculate the corresponding $l_k$, $h_k$, and $\\hat{h}_k$, respectively.\n",
        "\n",
        "* Function `l(k, x)` has already been defined for you. Your task is to complete the `h(k, x)`, `h_hat(k, x)`, and `hermit(x, y, y_prime)` functions.\n",
        "\n",
        "* Later we will draw some plots to check if the code is working.\n",
        "\n",
        "---"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "PPIUqb63Plvc",
        "nbgrader": {
          "cell_type": "markdown",
          "checksum": "b73facf56ebe52af4d34ce6aed213204",
          "grade": false,
          "grade_id": "part01_intro",
          "locked": true,
          "schema_version": 3,
          "solution": false,
          "task": false
        }
      },
      "source": [
        "### Part 1: Calculate $l_k$\n",
        "This function uses the following equation to calculate $l_k(x)$ and returns a polynomial:\n",
        "\n",
        "$$l_k(x) = \\prod_{j=0, j\\neq k}^{n} \\frac{x-x_j}{x_k-x_j}.$$"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "hLF2Df6LPlvd",
        "nbgrader": {
          "cell_type": "code",
          "checksum": "193776fb462b953e4cc901203895dce4",
          "grade": false,
          "grade_id": "l_ans",
          "locked": true,
          "schema_version": 3,
          "solution": false,
          "task": false
        }
      },
      "outputs": [],
      "source": [
        "# Already written for you.\n",
        "\n",
        "def l(k, x):\n",
        "    n = len(x)\n",
        "    assert (k < len(x))\n",
        "    \n",
        "    x_k = x[k]\n",
        "    x_copy = np.delete(x, k)\n",
        "    \n",
        "    denominator = np.prod(x_copy - x_k)\n",
        "    \n",
        "    coeff = []\n",
        "    \n",
        "    for i in range(n):\n",
        "        coeff.append(sum([np.prod(x) for x in combinations(x_copy, i)]) * (-1)**(i) / denominator)\n",
        "    \n",
        "    coeff.reverse()\n",
        "    \n",
        "    return Polynomial(coeff)"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "WDJyus6iPlve",
        "nbgrader": {
          "cell_type": "markdown",
          "checksum": "2b2776844ce83abcacb6b0ba6028ccd0",
          "grade": false,
          "grade_id": "part02_intro",
          "locked": true,
          "schema_version": 3,
          "solution": false,
          "task": false
        }
      },
      "source": [
        "### Part 2: Calculate $h_k$\n",
        "This function calculates $h_k(x)$ using the following equation:\n",
        "$$h_k(x) = \\left(1 + 2x_kl^\\prime_k(x_k) - 2l'_k(x_k)x \\right) l^2_k(x_k).$$\n",
        "\n",
        "This equation is basically a multiplication of two polynomials.\n",
        "\n",
        "First polynomial: $1 + 2x_kl^\\prime_k(x_k) - 2l'_k(x_k)x$.\n",
        "\n",
        "Second polynomial: $l^2_k(x_k)$.\n",
        "\n",
        "The `coeff` variable should contain a python list of coefficient values for the **first** polynomial of the equation. These coefficient values are used to create a polynomial `p`."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "deletable": false,
        "id": "QFOHvxn4Plve",
        "nbgrader": {
          "cell_type": "code",
          "checksum": "b40e0245554b77b7aeb6e6fada5c28f9",
          "grade": false,
          "grade_id": "h_ans",
          "locked": false,
          "schema_version": 3,
          "solution": true,
          "task": false
        }
      },
      "outputs": [],
      "source": [
        "def h(k, x):\n",
        "    # initialize with None. Replace with appropriate values/function calls\n",
        "    l_k = l(k, x)\n",
        "    l_k_sqr = l(k, x) * l(k, x)\n",
        "    l_k_prime = l_k.deriv(1)\n",
        "    coeff = [1 + 2*x[k]*l_k_prime(x[k]), -2*l_k_prime(x[k])]\n",
        "    p = Polynomial(coeff)\n",
        "  \n",
        "   \n",
        "    # --------------------------------------------\n",
        "    # YOUR CODE HERE\n",
        "    \n",
        "    # --------------------------------------------\n",
        "    \n",
        "    return p * l_k_sqr"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "GzGTUWJ1Plvf",
        "nbgrader": {
          "cell_type": "code",
          "checksum": "084dac03acc384f3952f1decde3eda17",
          "grade": true,
          "grade_id": "h_test",
          "locked": true,
          "points": 1.5,
          "schema_version": 3,
          "solution": false,
          "task": false
        }
      },
      "outputs": [],
      "source": [
        "# Test case for the h(k, x) function\n",
        "x = [3, 5, 7, 9]\n",
        "k = 2\n",
        "h_test = h(k, [3, 5, 7, 9])\n",
        "h_result = Polynomial([-2.5, 0.5]) * (l(k, x) ** 2)\n",
        "\n",
        "assert Polynomial.has_samecoef(h_result, h_test)\n",
        "assert h_result == h_test"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "M724n62yPlvf",
        "nbgrader": {
          "cell_type": "markdown",
          "checksum": "e95ffdba5dfcaf981dc8b884538cb7fa",
          "grade": false,
          "grade_id": "part03_intro",
          "locked": true,
          "schema_version": 3,
          "solution": false,
          "task": false
        }
      },
      "source": [
        "### Part 3: Calculate $\\hat{h}_k$\n",
        "This function calculates $\\hat{h}_k(x)$ using the following equation:\n",
        "\n",
        "$$\\hat{h}_k(x) = (x-x_k)l^2_k(x_k).$$\n",
        "\n",
        "This equation is also a multiplication of two polynomials.\n",
        "\n",
        "First polynomial: $x-x_k$.\n",
        "\n",
        "Second polynomial:  $l^2_k(x_k)$.\n",
        "\n",
        "The `coeff` variable should contain a python list of coefficient values for the **first** polynomial of the equation. These coefficient values are used to create a polynomial `p`."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "deletable": false,
        "id": "R0VeC9AIPlvf",
        "nbgrader": {
          "cell_type": "code",
          "checksum": "526567f27fdc6b0c693944715388aeb3",
          "grade": false,
          "grade_id": "h_hat_ans",
          "locked": false,
          "schema_version": 3,
          "solution": true,
          "task": false
        }
      },
      "outputs": [],
      "source": [
        "def h_hat(k, x):\n",
        "    # Initialize with none\n",
        "    l_k = l(k, x)\n",
        "    l_k_sqr = l(k, x)* l(k, x)\n",
        "    coeff = [-x[k], 1]\n",
        "    p = Polynomial(coeff)\n",
        "    \n",
        "    # --------------------------------------------\n",
        "    # YOUR CODE HERE\n",
        "    # --------------------------------------------\n",
        "    \n",
        "    return p * l_k_sqr"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "URbFerVBPlvg",
        "nbgrader": {
          "cell_type": "code",
          "checksum": "dab48dc3e3425ad92aac1fdc643320cd",
          "grade": true,
          "grade_id": "h_hat_test",
          "locked": true,
          "points": 1.5,
          "schema_version": 3,
          "solution": false,
          "task": false
        }
      },
      "outputs": [],
      "source": [
        "x = [3, 5, 7, 9]\n",
        "k = 2\n",
        "h_test = h_hat(k, [3, 5, 7, 9])\n",
        "h_result = Polynomial([-7, 1]) * (l(k, x) ** 2)\n",
        "\n",
        "assert Polynomial.has_samecoef(h_result, h_test)\n",
        "assert h_result == h_test"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "yJzRxYHuPlvg",
        "nbgrader": {
          "cell_type": "markdown",
          "checksum": "cd0d1d10b688aec869dfd36a17eda759",
          "grade": false,
          "grade_id": "part04_intro",
          "locked": true,
          "schema_version": 3,
          "solution": false,
          "task": false
        }
      },
      "source": [
        "### Part 4: The Hermite Polynomial\n",
        "This function uses the following equation:\n",
        "\n",
        "$$p_{2n+1} = \\sum_{k=0}^{n} \\left(f(x_k)h_k(x) + f'(x_k)\\hat{h}_k(x)\\right).$$\n",
        "\n",
        "The polynomial denoted by the equation is calculated by the variable `f`."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "deletable": false,
        "id": "CVWuSAIqPlvg",
        "nbgrader": {
          "cell_type": "code",
          "checksum": "238ae82cc84e0f1db6ee42946346ae4a",
          "grade": false,
          "grade_id": "part04_ans",
          "locked": false,
          "schema_version": 3,
          "solution": true,
          "task": false
        }
      },
      "outputs": [],
      "source": [
        "def hermit(x, y, y_prime):                      #x = nodes\n",
        "    assert len(x) == len(y)\n",
        "    assert len(y) == len(y_prime)\n",
        "    \n",
        "    f = Polynomial([0])                        #main formula of summation\n",
        "    # --------------------------------------------\n",
        "    # YOUR CODE HERE\n",
        "    for i in range(len(x)):                  ##as the node no is len(x)\n",
        "      f += y[i]* h(i, x) + y_prime[i]* h_hat(i, x)            #summation\n",
        "\n",
        "\n",
        "   \n",
        "    # --------------------------------------------\n",
        "    return f"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "BZyouHLCPlvg",
        "nbgrader": {
          "cell_type": "markdown",
          "checksum": "be144e5dff21332b4b93a63e59a035c2",
          "grade": false,
          "grade_id": "plot_intro",
          "locked": true,
          "schema_version": 3,
          "solution": false,
          "task": false
        }
      },
      "source": [
        "## Testing our methods by plotting graphs.\n",
        "\n",
        "**Note:** \n",
        "\n",
        "* For each of the 5 plots, there will be 2 curves plotted: one being the original function, and the other being the interpolated curve. \n",
        "\n",
        "* The original functions are displayed in orange color, while the hermite interpolated curves are in blue.\n",
        "\n",
        "* `x`, `y`, and `y_prime` contain $x_i$, $f(x_i)$, and $f'(x_i)$ of the given nodes of the original function $f$.\n",
        "\n",
        "Upon calling the `hermit()` function, it returns a polynomial `f`. For example, for plot 1, it is called `f3`.\n",
        "\n",
        "In general, a polynomial may look like the following: $f = 1 + 2x + 3x^2$. Next, we pass in a number of $x$ values to the polynomial by calling the `.linspace()` function on the polynomial object using `f.linspace()`. This function outputs a tuple, which is stored in a variable called `data`. First element of `data` contains a 1D numpy array of $x_i$ values generated by `linspace()`, and the second element of `data` contains a 1D numpy array of the corresponding $y_i$ values outputted by our example polynomial:\n",
        "$f = 1 + 2x + 3x^2$. \n",
        "\n",
        "Using `test_x`, we generate a range of $x_i$ values to plot the original function, and `test_y` contains the corresponding $y_i$ values of the original function. For the first plot, our original function is the *sine curve*.\n",
        "\n",
        "For all the plots:\n",
        "\n",
        "`plt.plot(test_x, test_y)` plots the original function.\n",
        "\n",
        "`plt.plot(data[0], data[1])` plots the interpolated polynomial."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "I0AAKqRiPlvh",
        "nbgrader": {
          "cell_type": "code",
          "checksum": "9e619e399993ba9b52b9a5acfc6d952a",
          "grade": false,
          "grade_id": "samples",
          "locked": true,
          "schema_version": 3,
          "solution": false,
          "task": false
        }
      },
      "outputs": [],
      "source": [
        "pi      = np.pi\n",
        "x       = np.array([0.0, pi/8.0,  pi, 4.0*pi/2.0])\n",
        "y       = np.array([0.0,    1.0, 0.0,       -1.0])\n",
        "y_prime = np.array([1.0,    0.0, 1.0,        0.0])"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "CE1RRR6pPlvh",
        "nbgrader": {
          "cell_type": "markdown",
          "checksum": "62a455aacd036963d695b540c90f7cd0",
          "grade": false,
          "grade_id": "plot01",
          "locked": true,
          "schema_version": 3,
          "solution": false,
          "task": false
        }
      },
      "source": [
        "**Plot 1:** trying to interpolate a sine curve (`np.sin()`) using first 2 nodes in `x` and `y`, and their corresponding derivative in `y_prime`."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "QQPpgFHEPlvh",
        "nbgrader": {
          "cell_type": "code",
          "checksum": "6d5bb5d8680003a78893653e98f6a087",
          "grade": true,
          "grade_id": "plot01_test",
          "locked": true,
          "points": 1,
          "schema_version": 3,
          "solution": false,
          "task": false
        },
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 286
        },
        "outputId": "62a1f835-2587-4671-dc77-5851180b2949"
      },
      "outputs": [
        {
          "output_type": "display_data",
          "data": {
            "text/plain": [
              "<Figure size 432x288 with 1 Axes>"
            ],
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAD4CAYAAAAD6PrjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAcwUlEQVR4nO3dfXRU933n8fdXEhKP5kE8mmewwMF2bWMFY5o4fkhs7N2GOJtkcbNrkvgsSdbebLvdNs76NEmTk3OSbpue+MRxDm1onW4d20lKTBNsjB0nPm5iDJgngQEJgZEEEnoASQiQkPTdP+YiD8PoCc3VnYfP65zR3PndO5rv/OZKn3t/986MuTsiIiIAeVEXICIi6UOhICIiPRQKIiLSQ6EgIiI9FAoiItKjIOoChmry5Mk+b968qMsQEckoO3bsaHD3KYntGR8K8+bNY/v27VGXISKSUczs3WTtGj4SEZEeCgUREemRklAws/VmdtLMyuLaJpnZFjMrD64nBu1mZk+YWYWZ7TGzpXH3WRMsX25ma1JRm4iIDFyq9hT+CViZ0PYY8Kq7lwCvBrcB7gNKgsta4CmIhQjwNeBWYBnwtYtBIiIiwyMloeDurwNNCc2rgKeD6aeBj8W1/9hj3gQmmNkM4F5gi7s3ufspYAuXB42IiIQozGMK09z9RDBdC0wLpmcCVXHLVQdtvbVfxszWmtl2M9teX1+f2qpFRHLYsBxo9thHsabs41jdfZ27l7p76ZQpl51mKyIiVyjMUKgLhoUIrk8G7TXA7LjlZgVtvbWLyDD57pZD/P5wY9RlSITCDIWNwMUziNYAL8S1PxSchbQcaA6GmTYD95jZxOAA8z1Bm4gMg9bzF3ji1XJ2Vp2KuhSJUEre0WxmPwHuACabWTWxs4i+DTxvZg8D7wKfChbfBNwPVABngc8CuHuTmX0T2BYs9w13Tzx4LSIhOVR3BoDF08ZFXIlEKSWh4O4P9jLr7iTLOvBIL79nPbA+FTWJyOAcqmsFYJFCIafpHc0iAsDB2lbGFOYzc8KoqEuRCCkURASIhULJtHHk5VnUpUiEFAoiAsSGj3Q8QRQKIkLDmXYa2zpYNF2hkOsUCiLCodrYQWbtKYhCQUQ4ePHMo+ljI65EoqZQEBEO1bUyaUwhU8YWRV2KREyhICIcrG1l0bSxmOnMo1ynUBDJcd3dzoFanXkkMQoFkRz3btNZznZ0cd3V46MuRdKAQkEkx+0/3gLAkquvirgSSQcKBZEct/9EMwV5Rsk0nXkkCgWRnLf/eAvXTB1LUUF+1KVIGlAoiOS4/SdaNHQkPRQKIjms4Uw7dS3tLJmhUJAYhYJIDnvnhA4yy6UUCiI5rOfMI+0pSEChIJLD9h1vYeaEUUwYXRh1KZImFAoiOaysppnrNHQkcRQKIjmq+dwFKhvauHH2hKhLkTSiUBDJUWU1zQDcMFMfbyHvUSiI5Kjd1acB+INZCgV5j0JBJEftqWpmbvFoHWSWSygURHLU3ppmDR3JZRQKIjmo4Uw7NafPceMsHWSWSykURHLQHh1PkF4oFERy0O6qZszgeg0fSQKFgkgO2ll1mpKpYxlTVBB1KZJmFAoiOaa729n57ilumTsp6lIkDSkURHJM+ckztLZ3Ujp3YtSlSBpSKIjkmB3vngLgFoWCJKFQEMkxO949RfGYQuYWj466FElDCgWRHPP2sVMsnTsRM4u6FElDCgWRHNJ4pp0jDW0aOpJeKRREcoiOJ0h/FAoiOWTb0SYK8/P0mUfSK4WCSA55s7KJm+ZMYOSI/KhLkTQVeiiY2VEz22tmu8xse9A2ycy2mFl5cD0xaDcze8LMKsxsj5ktDbs+kVzRfO4C+443c9uC4qhLkTQ2XHsKd7r7Te5eGtx+DHjV3UuAV4PbAPcBJcFlLfDUMNUnkvW2HWmi22G5QkH6ENXw0Srg6WD6aeBjce0/9pg3gQlmNiOKAkWyzZuVjRQW5HHzHH1ctvRuOELBgZfNbIeZrQ3aprn7iWC6FpgWTM8EquLuWx20XcLM1prZdjPbXl9fH1bdIlnlzSONLNXxBOnHcITCB9x9KbGhoUfM7Pb4me7uxIJjwNx9nbuXunvplClTUliqSHZqPnuBfcdbNHQk/Qo9FNy9Jrg+CWwAlgF1F4eFguuTweI1wOy4u88K2kRkCLYeacR1PEEGINRQMLMxZjbu4jRwD1AGbATWBIutAV4IpjcCDwVnIS0HmuOGmUTkCr1eXs/ownyWztGb1qRvYX/DxjRgQ/AZKwXAM+7+kpltA543s4eBd4FPBctvAu4HKoCzwGdDrk8k67k7vz1Uz4qFxRQW6K1J0rdQQ8HdK4Ebk7Q3AncnaXfgkTBrEsk1RxvPUtV0jrUfXBB1KZIBtNkgkuV+ezB2yO72RTopQ/qnUBDJcq+XNzCveDRzi8dEXYpkAIWCSBZr7+zi94cb+ZD2EmSAFAoiWezNyibOXejiQ4sVCjIwCgWRLPbyvlpGF+azYuHkqEuRDKFQEMlS3d3Olv113LF4ij7aQgZMoSCSpXZXn+Zkazv3LJkedSmSQRQKIlnq5f11FOQZdy6eGnUpkkEUCiJZavO+WpYvKGb86BFRlyIZRKEgkoUO1LZQWd/GvddN639hkTgKBZEs9MKu4+TnGfffoO+oksFRKIhkme5uZ+Ou43ywZDLFY4uiLkcyjEJBJMvsOHaKmtPnWHXT1VGXIhlIoSCSZX6xs4aRI/J0KqpcEYWCSBY5f6GLX+09wUeWTGdMUdhflyLZSKEgkkVeKqvl9NkL/OfS2f0vLJKEQkEkizyz9Rhzi0ezYqG+i1mujEJBJEtUnGzlraNNPLhsDnl5FnU5kqEUCiJZ4pmtVYzINz5xy6yoS5EMplAQyQKt5y/w0+1VrLx+BpP13gQZAoWCSBb4yVvHaG3vZO0HF0RdimQ4hYJIhuvo7Gb9G0dZsbCYG2aNj7ocyXAKBZEMt3H3cWpbzrP2du0lyNApFEQy2IWubr7/63LeN+MqPrRI38MsQ6dQEMlgz22r4mjjWf783kWY6TRUGTqFgkiGOtfRxfdeLef98ybq29UkZRQKIhlq3euV1Le28+WV12ovQVJGoSCSgY40tPHkbyr4DzfMoHTepKjLkSyiUBDJMO7OX/6ijKL8PL76R0uiLkeyjEJBJMP8dEc1b1Q08OcrFzPtqpFRlyNZRqEgkkEqTrbytRf2sXzBJD5969yoy5EspFAQyRDnOrp49JmdjC7M53urbyZfn4QqIdBXM4lkgM6ubv7HT3ZysK6Vf/zM+zVsJKHRnoJImnN3vrZxH6+8U8fX/+g67tB7EiRE2lMQSWOdXd08vqGM57ZX8YUPLWTNinlRlyRZTqEgkqaaz13gz57fzSvv1PGlu67hTz+yKOqSJAcoFETS0FtHmvjT53ZR13Keb6y6jodumxd1SZIj0i4UzGwl8D0gH/gHd/92xCWJDJsjDW38zcsH+dWeE8yZNJqffXEFN82eEHVZkkPSKhTMLB94EvgIUA1sM7ON7r4/2spEwtPU1sFvD51kw87jvH6onpEj8vjS3SV8/vYFjClKqz9RyQHptsYtAyrcvRLAzJ4FVgGpD4VDm6G9FcwAi11bXjCdF1wS24i7bXHt1k8bsem4q1558MP94o33pt3Buy+d790J7cRNJ87rvvx+kNDey3VPfXHTiS75UDZL3geJfZPYZ32+Bv31e9yyPfXEze+t/z1houc59tEfcf3o7rSd76Dl3AXa2jvouNBFV1cX7hfnx5bx7i7aO7tp77hA87kOms60c/x0Gw2t7eDOwlEF/Jcbilm+YCJXjWyE/W9dWsuAJT7nPtbLpH2cl/A6JP599PMaJLuOK633/o/v8/76P9m6Td/r+mXzGPi639d63/PcEv/W+1tXr/R/T9zvm7sC8vL7r20Q0i0UZgJVcbergVsTFzKztcBagDlz5lzZI738l9Bw8MruKxLHgLHB5YqMCK67gPLgIjIQj9dlfSgMiLuvA9YBlJaWDmYT6j2ffh46O0i+VZG4ZZG4lZ5sizthCyPpVkbc7T4/6jh+KyK+LXErO8mWQ7Jlet3Ki59HL1t5JEwHty+TuDfR195Osj2fXrbs+t3aS/b7gnou2etKqDFZn/f0Qfzt957/qXMXeLGsjtcPNdJ8/gJ5eXnMnzyWOcVjmTFhNONHFzJ2VCFFBQUU5OeRl58PxPrY8oz8vDxGjihkVFEBV40cEcyP32rvo/8vqa0Pl/RJ/PNP0ve9vQa9bX13dyXct6+t7MS+T+j/y/4G4vcmEvfs+tqqTrI+x2+JX7bHk5fk98T178W/hb5eh+QdH9f/yfq9+/LX55L/LwnTl93u5e8gv7CPmq5MuoVCDTA77vasoC31Js4L5ddK9jl/oYsnXi1n/b8foaOzmI8sWcJnls7i9pIpjCpM7VaaSNTSLRS2ASVmNp9YGKwG/jjakiSX7Tx2ij97fjeVDW08cPNMvnR3CfMnj4m6LJHQpFUouHunmT0KbCZ2Sup6d98XcVmSozbsrObLP9vLlHFF/L+Hb+UDJZOjLkkkdGkVCgDuvgnYFHUdktt++NvDfPvFAyxfMImnPn0LE8ekfuxWJB2lXSiIRO1Hbxzh2y8e4KM3Xs3ffPJGCgv0uZGSOxQKInE27Kzmm7/cz33XT+e7n7qRgnwFguQWrfEigbKaZh77+V6WL5jE91bfrECQnKS1XgQ4fbaDz//zDorHFPL9P16qISPJWRo+EgH+6t/2U9dynp9/cQWTxxZFXY5IZLQ5JDlvy/46Nuys4ZE7r+FGfSKp5DiFguS0lvMX+D8b9vK+GVfxyJ3XRF2OSOQ0fCQ57cnXKmg4086P1pTqOIII2lOQHHas8Sz/+MZRPn7zLP5gloaNREChIDnsOy8dID/P+IuVi6MuRSRtKBQkJ+0/3sKv9p7gv92+gGlXjYy6HJG0oVCQnPT918oZV1TAwx+YH3UpImlFoSA5p7yulRfLalmzYh7jR43o/w4iOUShIDnnB785zKgR+XxOewkil1EoSE6paznPv+0+zur3z2GSPg5b5DIKBckpz2w9Rpc7D902N+pSRNKSQkFyRkdnN8+8dYw7Fk1hnr5SUyQphYLkjBfLTlDf2s6aFfOiLkUkbSkUJGc8s/UYc4tHc3vJlKhLEUlbCgXJCVVNZ9l6pIlP3jKLvDyLuhyRtKVQkJzw87erMYMHls6KuhSRtKZQkKzn7vzr2zXctqCYmRNGRV2OSFpTKEjW23b0FMeazvKftJcg0i+FgmS9X+yqYXRhPiuvnx51KSJpT6EgWa2zq5vNZbXcde1UxhTpO6VE+qNQkKz21tEmGts6uP+GGVGXIpIRFAqS1V7cW8vIEXncsVjvTRAZCIWCZK2ubuelfbXcuXgqows1dCQyEAoFyVo73j1FfWs792noSGTAFAqStbbsr2VEvnHXtVOjLkUkYygUJGv9+sBJbp1fzFiddSQyYAoFyUrHGs9yuL6NO7WXIDIoCgXJSr8+UAegoSORQVIoSFZ67WA98yePYb6+TEdkUBQKknXOdnTy+8pG7lysvQSRwVIoSNZ5s7KRjs5uDR2JXAGFgmSdN8obKSrIo3TexKhLEck4oYWCmX3dzGrMbFdwuT9u3lfMrMLMDprZvXHtK4O2CjN7LKzaJLv97nADpfMmMnJEftSliGScsPcU/s7dbwoumwDMbAmwGrgOWAn8wMzyzSwfeBK4D1gCPBgsKzJgDWfaOVDbyoqFk6MuRSQjRfGunlXAs+7eDhwxswpgWTCvwt0rAczs2WDZ/RHUKBnqd4cbAfjDaxQKIlci7D2FR81sj5mtN7OLA7wzgaq4ZaqDtt7aRQbsdxUNjBtZwA0zx0ddikhGGlIomNkrZlaW5LIKeApYCNwEnAD+NgX1XnzctWa23cy219fXp+rXShb498MNLF9QTH6eRV2KSEYa0vCRu394IMuZ2d8Dvwxu1gCz42bPCtrooz3xcdcB6wBKS0t9ECVLFqtqOktV0zke/sP5UZcikrHCPPso/vOKHwDKgumNwGozKzKz+UAJ8BawDSgxs/lmVkjsYPTGsOqT7LPtaBMAty4ojrgSkcwV5oHmvzazmwAHjgKfB3D3fWb2PLEDyJ3AI+7eBWBmjwKbgXxgvbvvC7E+yTLbjjYxbmQBi6aNi7oUkYwVWii4+3/tY963gG8lad8EbAqrJslu246eonTuRB1PEBkCvaNZskJTWwcVJ89QOm9S1KWIZDSFgmSF7cHxhGXzFQoiQ6FQkKyw7WgThfl5en+CyBApFCQrbDt6ihtnj9fnHYkMkUJBMt75C12U1TRzy1wNHYkMlUJBMl5ZTTOd3c7SOROiLkUk4ykUJOPtqjoNwE0KBZEhUyhIxttVdZqZE0YxddzIqEsRyXgKBcl4u6pOc9Ns7SWIpIJCQTJaw5l2qk+dUyiIpIhCQTLarmOx4wk3KhREUkKhIBltV9Vp8vNMb1oTSRGFgmS0XVWnWTxtHKMK9aY1kVRQKEjGcnf2VJ/W0JFICikUJGNVnzpHy/lODR2JpJBCQTJWWU0zANfPvCriSkSyh0JBMlbZ8Wby80zftCaSQgoFyVhlNS2UTB2rT0YVSSGFgmQkd6esppnrdTxBJKUUCpKR6lraaWzr4PqrdTxBJJUUCpKR3jvIrD0FkVRSKEhGKjvejBm8b4b2FERSSaEgGamspoX5k8cwpqgg6lJEsopCQTLSgdoWlmgvQSTlFAqScVrPX6D61DkNHYmEQKEgGedQXSsAi/WmNZGUUyhIxnnnRCwUrp2hUBBJNYWCZJyDta2MKypg5oRRUZciknUUCpJxDtS2sHj6OMws6lJEso5CQTKKu3OgtlVDRyIhUShIRjnefJ7W851cO11nHomEQaEgGeXAiRYArp2uPQWRMCgUJKMcqI2debRIoSASCoWCZJRDda3MnDCKq0aOiLoUkaykUJCMUl53hpJpY6MuQyRrKRQkY3R1O4frz1AyVaEgEhaFgmSM6lNnae/spmSqjieIhEWhIBmjvO4MAAu1pyASmiGFgpl90sz2mVm3mZUmzPuKmVWY2UEzuzeufWXQVmFmj8W1zzezrUH7c2ZWOJTaJPuUn4yFwjUKBZHQDHVPoQz4OPB6fKOZLQFWA9cBK4EfmFm+meUDTwL3AUuAB4NlAb4D/J27XwOcAh4eYm2SZSpOnmHaVUWMH6Uzj0TCMqRQcPd33P1gklmrgGfdvd3djwAVwLLgUuHule7eATwLrLLYh9jcBfwsuP/TwMeGUptkn4qTrTqeIBKysI4pzASq4m5XB229tRcDp929M6E9KTNba2bbzWx7fX19SguX9OTulJ88o6EjkZD1+wW3ZvYKMD3JrMfd/YXUl9Q/d18HrAMoLS31KGqQ4XW8+TxnO7r0HgWRkPUbCu7+4Sv4vTXA7Ljbs4I2emlvBCaYWUGwtxC/vAjlwbetafhIJFxhDR9tBFabWZGZzQdKgLeAbUBJcKZRIbGD0Rvd3YHXgE8E918DRLIXIumpQmceiQyLoZ6S+oCZVQO3Ab8ys80A7r4PeB7YD7wEPOLuXcFewKPAZuAd4PlgWYAvA//LzCqIHWP40VBqk+xyuL6NCaNHMGmMzlQWCVO/w0d9cfcNwIZe5n0L+FaS9k3ApiTtlcTOThK5zJGGMyyYPCbqMkSynt7RLBmhsr6NBVM0dCQSNoWCpL0z7Z2cbG1nvvYUREKnUJC0d6S+DYCFUxQKImFTKEjaq2yInXk0f7KGj0TCplCQtFdZ34YZzC0eHXUpIllPoSBp70hDG7MmjmLkiPyoSxHJegoFSXuVDWc0dCQyTBQKktbcnSP1bXqPgsgwUShIWjvZ2k5bR5fOPBIZJgoFSWuVwemoGj4SGR4KBUlrI0fkcc+SafrIbJFhMqTPPhIJ281zJrLuodL+FxSRlNCegoiI9FAoiIhID4WCiIj0UCiIiEgPhYKIiPRQKIiISA+FgoiI9FAoiIhID3P3qGsYEjOrB969wrtPBhpSWE6qqK7BUV2Do7oGJ1vrmuvuUxIbMz4UhsLMtrt72r1dVnUNjuoaHNU1OLlWl4aPRESkh0JBRER65HoorIu6gF6orsFRXYOjugYnp+rK6WMKIiJyqVzfUxARkTgKBRER6ZH1oWBmnzSzfWbWbWalCfO+YmYVZnbQzO7t5f7zzWxrsNxzZlYYQo3Pmdmu4HLUzHb1stxRM9sbLLc91XUkebyvm1lNXG3397LcyqAPK8zssWGo6/+a2QEz22NmG8xsQi/LDUt/9ff8zawoeI0rgnVpXli1xD3mbDN7zcz2B+v//0yyzB1m1hz3+n417LqCx+3zdbGYJ4L+2mNmS4ehpsVx/bDLzFrM7E8SlhmW/jKz9WZ20szK4tommdkWMysPrif2ct81wTLlZrbmigpw96y+AO8DFgO/AUrj2pcAu4EiYD5wGMhPcv/ngdXB9A+BL4Zc798CX+1l3lFg8jD23deB/93PMvlB3y0ACoM+XRJyXfcABcH0d4DvRNVfA3n+wH8HfhhMrwaeG4bXbgawNJgeBxxKUtcdwC+Ha30a6OsC3A+8CBiwHNg6zPXlA7XE3tw17P0F3A4sBcri2v4aeCyYfizZOg9MAiqD64nB9MTBPn7W7ym4+zvufjDJrFXAs+7e7u5HgApgWfwCZmbAXcDPgqangY+FVWvweJ8CfhLWY4RgGVDh7pXu3gE8S6xvQ+PuL7t7Z3DzTWBWmI/Xj4E8/1XE1h2IrUt3B691aNz9hLu/HUy3Au8AM8N8zBRaBfzYY94EJpjZjGF8/LuBw+5+pZ+UMCTu/jrQlNAcvw719n/oXmCLuze5+ylgC7BysI+f9aHQh5lAVdztai7/oykGTsf9A0q2TCp9EKhz9/Je5jvwspntMLO1IdYR79FgF359L7usA+nHMH2O2FZlMsPRXwN5/j3LBOtSM7F1a1gEw1U3A1uTzL7NzHab2Ytmdt0wldTf6xL1OrWa3jfMougvgGnufiKYrgWmJVkmJf1WMPja0o+ZvQJMTzLrcXd/YbjrSWaANT5I33sJH3D3GjObCmwxswPBVkUodQFPAd8k9kf8TWJDW58byuOloq6L/WVmjwOdwL/08mtS3l+ZxszGAj8H/sTdWxJmv01siORMcLzoF0DJMJSVtq9LcMzwo8BXksyOqr8u4e5uZqG9lyArQsHdP3wFd6sBZsfdnhW0xWsktutaEGzhJVsmJTWaWQHwceCWPn5HTXB90sw2EBu6GNIf00D7zsz+HvhlklkD6ceU12VmnwH+I3C3BwOqSX5HyvsriYE8/4vLVAev83hi61aozGwEsUD4F3f/18T58SHh7pvM7AdmNtndQ/3wtwG8LqGsUwN0H/C2u9clzoiqvwJ1ZjbD3U8EQ2knkyxTQ+y4x0WziB1LHZRcHj7aCKwOzgyZTyzx34pfIPhn8xrwiaBpDRDWnseHgQPuXp1sppmNMbNxF6eJHWwtS7ZsqiSM4z7Qy+NtA0osdpZWIbFd740h17US+Avgo+5+tpdlhqu/BvL8NxJbdyC2Lv26tyBLleCYxY+Ad9z9u70sM/3isQ0zW0bs/0GoYTXA12Uj8FBwFtJyoDlu6CRsve6tR9FfceLXod7+D20G7jGzicFQ7z1B2+CEfSQ96guxf2bVQDtQB2yOm/c4sTNHDgL3xbVvAq4OphcQC4sK4KdAUUh1/hPwhYS2q4FNcXXsDi77iA2jhN13/wzsBfYEK+WMxLqC2/cTO7vl8DDVVUFs7HRXcPlhYl3D2V/Jnj/wDWKhBTAyWHcqgnVpwTD00QeIDfvtieun+4EvXFzPgEeDvtlN7ID9imGoK+nrklCXAU8G/bmXuLMGQ65tDLF/8uPj2oa9v4iF0gngQvC/62Fix6BeBcqBV4BJwbKlwD/E3fdzwXpWAXz2Sh5fH3MhIiI9cnn4SEREEigURESkh0JBRER6KBRERKSHQkFERHooFEREpIdCQUREevx/ztCfF7cAX9wAAAAASUVORK5CYII=\n"
          },
          "metadata": {
            "needs_background": "light"
          }
        }
      ],
      "source": [
        "n      = 1\n",
        "f3     = hermit(x[:(n+1)], y[:(n+1)], y_prime[:(n+1)])\n",
        "data   = f3.linspace(n=50, domain=[-3, 3])\n",
        "test_x = np.linspace(-10, 10, 50, endpoint=True)\n",
        "test_y = np.sin(test_x)\n",
        "\n",
        "plt.plot(data[0], data[1])\n",
        "plt.plot(test_x, test_y)\n",
        "plt.show()"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "W7v9mEcLPlvh",
        "nbgrader": {
          "cell_type": "markdown",
          "checksum": "d8f7bf96ae4f35419f880d25e72cff09",
          "grade": false,
          "grade_id": "plot02",
          "locked": true,
          "schema_version": 3,
          "solution": false,
          "task": false
        }
      },
      "source": [
        "**Plot 2:** trying to interpolate a sine curve (`np.sin()`) using first 3 nodes in `x` and `y` and their corresponding derivative in `y_prime`."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "qnNNOyC5Plvh",
        "nbgrader": {
          "cell_type": "code",
          "checksum": "25f490548e99af7d5fccb9fd5c653e80",
          "grade": true,
          "grade_id": "plot02_ans",
          "locked": true,
          "points": 1,
          "schema_version": 3,
          "solution": false,
          "task": false
        },
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 286
        },
        "outputId": "4caa6a28-5f70-48a1-fa53-c0d9304ca567"
      },
      "outputs": [
        {
          "output_type": "display_data",
          "data": {
            "text/plain": [
              "<Figure size 432x288 with 1 Axes>"
            ],
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAD4CAYAAAAJmJb0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3deXAk53nf8e+DGdw3BrMg9iKwu1gey4gXtLyXsiWbtKWY8iXTdhSmQhfjKtlRyq44kpMqO4kVW1Zsy7GtOIyliH9YkllOGDE+ZFO0LHJ5Y0VyucslF9gb2AP3fc7gyR89WIFL7IGdHjTQ+H2qUD3TPTP99ADzwztvv91t7o6IiMRTUdQFiIhI4SjkRURiTCEvIhJjCnkRkRhTyIuIxFgy6gIWa2xs9JaWlqjLEBFZU/bt29fv7umllq2qkG9paaGjoyPqMkRE1hQzO3GxZequERGJMYW8iEiMKeRFRGJMIS8iEmOh7Hg1s+PAGJAFMu7ebmYNwF8ALcBx4BPuPhTG+kRE5MqE2ZL/AXe/xd3bc/c/Azzr7m3As7n7IiKyggrZXfMQ8ETu9hPAxwu4LhERWUJYIe/A35vZPjN7LDevyd3P5G6fBZqWeqKZPWZmHWbW0dfXF1I5Istw+nV49j/D9GjUlYiELqyQv9fdbwN+BPiUme1ZvNCDk9YveeJ6d3/c3dvdvT2dXvKALZHCOnsAnv+vMD0SdSUioQsl5N29JzftBZ4CdgPnzKwZIDftDWNdIqFLlATT7Gy0dYgUQN4hb2aVZla9cBv4YeAA8DTwSO5hjwDfzHddIgWRyA0ym89EW4dIAYQxhLIJeMrMFl7va+7+LTN7DXjSzB4FTgCfCGFdIuErKg6m2blo6xApgLxD3t2PAjcvMX8A+HC+ry9ScImFkFd3jcSPjngVWQh5dddIDCnkRdRdIzGmkBc535JXyEv8KORF1JKXGFPIiyQU8hJfCnkRdddIjCnkRdRdIzGmkBdZOOJVIS8xpJAXWTh3jbprJIYU8iJFOuJV4kshL3J+dI2OeJX4UciL6Nw1EmMKeRH1yUuMKeRFNIRSYkwhL1JUBJZQd43EkkJeBIIuG7XkJYYU8iKgkJfYUsiLQDDCRt01EkMKeRHIteQV8hI/oYW8mSXM7HUz+6vc/VYze8XMuszsL8ysJKx1iYQukVR3jcRSmC35TwOHFt3/PPAH7r4DGAIeDXFdIuFSS15iKpSQN7PNwEeBP8vdN+AHgb/MPeQJ4ONhrEukIBIlOhhKYimslvwXgV8D5nP3U8Cwuy+cDKQb2LTUE83sMTPrMLOOvr6+kMoRWaZEsbprJJbyDnkz+xjQ6+77rub57v64u7e7e3s6nc63HJGro+4aialkCK9xD/BjZvajQBlQA/whUGdmyVxrfjPQE8K6RAqjSC15iae8W/Lu/ll33+zuLcDDwD+4+88D3wF+KvewR4Bv5rsukYJJFENmJuoqREJXyHHy/w74FTPrIuij/3IB1yWSn2SpdrxKLIXRXXOeu/8j8I+520eB3WG+vkjBJEogoz55iR8d8SoCOq2BxJZCXgQgUaqQl1hSyIuAWvISWwp5EdA4eYkthbwIBKNrtONVYkghLwLqrpHYUsiLwPd3vLpHXYlIqBTyIhD0yeMwn7nsQ0XWEoW8CATdNaAuG4kdhbwIBDteQeevkdhRyIuAWvISWwp5EQh2vIJCXmJHIS8Ci7prFPISLwp5EciNrgGy6pOXeFHIi4B2vEpsKeRFYFFLXt01Ei8KeRFQS15iSyEvAhpdI7GVd8ibWZmZvWpmb5rZQTP7j7n5rWb2ipl1mdlfmFlJ/uWKFEgy9+eplrzETBgt+RngB939ZuAW4EEzuxP4PPAH7r4DGAIeDWFdIoWx0JLPTEdbh0jI8g55D4zn7hbnfhz4QeAvc/OfAD6e77pECiapHa8ST6H0yZtZwszeAHqBZ4AjwLC7L5zSrxvYdJHnPmZmHWbW0dfXF0Y5IsuX0I5XiadQQt7ds+5+C7AZ2A1cv4znPu7u7e7enk6nwyhHZPmSZcFULXmJmVBH17j7MPAd4C6gzsySuUWbgZ4w1yUSKg2hlJgKY3RN2szqcrfLgR8CDhGE/U/lHvYI8M181yVSMAp5iank5R9yWc3AE2aWIPin8aS7/5WZvQ18w8x+C3gd+HII6xIpjKIkWJFG10js5B3y7r4fuHWJ+UcJ+udFVj+zoF9eJyiTmNERryILkqXqrpHYUciLLEiUqrtGYkchL7JALXmJIYW8yIJkmVryEjsKeZEFyVJd/k9iRyEvskAteYkhhbzIgqR2vEr8KORFFqglLzGkkBdZUFwGcwp5iReFvMgCteQlhhTyIgsU8hJDCnmRBcXlMDcVdRUioVLIiyzQEa8SQwp5kQXJcshMgXvUlYiERiEvsqA4dwlAteYlRhTyIgsWrvOqna8SIwp5kQULIa+drxIjCnmRBcUVwTSjkJf4CONC3lvM7Dtm9raZHTSzT+fmN5jZM2bWmZvW51+uSAEVlwdTteQlRsJoyWeAX3X3G4E7gU+Z2Y3AZ4Bn3b0NeDZ3X2T1Oh/y6pOX+Mg75N39jLt/L3d7DDgEbAIeAp7IPewJ4OP5rkukoM6H/GS0dYiEKNQ+eTNrAW4FXgGa3P1MbtFZoOkiz3nMzDrMrKOvry/MckSWZ6FPXt01EiOhhbyZVQH/G/g37j66eJm7O7DkESbu/ri7t7t7ezqdDqsckeU7P4RSIS/xEUrIm1kxQcD/ubv/n9zsc2bWnFveDPSGsS6RgtGOV4mhMEbXGPBl4JC7//6iRU8Dj+RuPwJ8M991iRTU+e4a9clLfCRDeI17gE8Cb5nZG7l5vw78DvCkmT0KnAA+EcK6RApHLXmJobxD3t33AnaRxR/O9/VFVkxJZTCdVUte4kNHvIosSBRDUTHMTURdiUhoFPIiixVXqCUvsaKQF1mspEIteYkVhbzIYsUV2vEqsaKQF1msRN01Ei8KeZHFiithdjzqKkRCo5AXWaykUgdDSawo5EUWK6mEWe14lfhQyIssVlKlkJdYUciLLFZSCTNjUVchEhqFvMhipWrJS7wo5EUWK6mE+TnIzEZdiUgoFPIii5VUBVMNo5SYUMiLLFZaHUzVLy8xoZAXWWyhJa+Ql5hQyIssttCSV3eNxIRCXmSx0ppgqpa8xIRCXmQx9clLzCjkRRY7H/Kj0dYhEpJQQt7MvmJmvWZ2YNG8BjN7xsw6c9P6MNYlUlBlue6aaYW8xENYLfmvAg9eMO8zwLPu3gY8m7svsrqVVIEVqSUvsRFKyLv7c8DgBbMfAp7I3X4C+HgY6xIpKLNg5+v0SNSViISikH3yTe5+Jnf7LNC01IPM7DEz6zCzjr6+vgKWI3KFymoV8hIbK7Lj1d0d8Isse9zd2929PZ1Or0Q5IpemkJcYKWTInzOzZoDctLeA6xIJT3kdTA1HXYVIKAoZ8k8Dj+RuPwJ8s4DrEglPeT1MDUVdhUgowhpC+XXgJeA6M+s2s0eB3wF+yMw6gY/k7ousfgp5iZFkGC/i7j97kUUfDuP1RVbUQsi7B6NtRNYwHfEqcqHy+uDCITpJmcSAQl7kQhWpYDo5EG0dIiFQyItcaCHkJxTysvYp5EUuVNEYTCf7o61DJAQKeZELVS605BXysvYp5EUuVLkhmE7o+D1Z+xTyIhcqrQrORjmukJe1TyEvspSqDTB2NuoqRPKmkBdZSnWzQl5iQSEvspSajTDaHXUVInlTyIsspWYTjJ6B+fmoKxHJi0JeZCl1W4JTG4yry0bWNoW8yFLqW4Lp0PEoqxDJm0JeZCn1rcF08Gi0dYjkSSEvspS6ayFRAn3vRl2JSF4U8iJLSSQh1Qa9h6KuRCQvCnmRi2n+AJx5M7h4iMgaFcqVoeTyJmYynBmZomd4mnOj08zMZZnJzDObnaeiOEF1WTHVZUk211dwbaqCytK1/6sZnZ6jZ2iKc6PT9I7OMD6TYd6dzLxTkiiiobKEuopimmrKaG2spKw4EXXJ77Xpdnjz68HO14bW788/9jyceQN2PgiNbe95yvDkLMcHJhkYn2FgYpaRyTkguMBUkRmpqhIaq0rZUF3K1lQFpclVts3L5O4MTc5xrH+cvrFZRqfmGJ2eIzPvJIuMkmQRNWXFXFNbRnNtGdfUlsVim/vGZzg1OMnp4WkmZjJMzGaZnstSmiyioiRJZWmC5tpyNteX01RTRqIouiuMFTxJzOxB4A+BBPBn7h77a71OzWb53skh9p0Y4uDpEQ6eHqV7aGpZr5GuLuXG5hpu3lLHLVtquW1rPXUVJQWqOH8jU3N87+QQHccHefv0KO+cHePMyPQVP98MNtWVc11TNbdurePWrfXcvKWOqij/2bXeH0yPPAsNvwDzWfibfwsdXwbA/+G3ePXeL/MPU9t5q3uEzt5x+sZmrvjlE0VGS6qCnU3V3Lyljg+21HPTptpVHYLdQ5N87+Qwb5wc5s3uYbp6xxmZmrvi5yeLjO3pKm5oruamTbXc0ZrihuZqkonV2ang7hzrn2DfiSHe6hlhf/cI754dY2oue8WvUZwwdjZVc9PGWm7aXMvulgbaNlRRtELBb17Ar6JmlgAOAz8EdAOvAT/r7m8v9fj29nbv6OgoWD2FksnO8/qpYZ4/3McLRwbY3z3MXNYxg9ZUJTdsrOGGa6rZ0lDBxrpyrqkpo7wkQWmyiOJEERMzGcamM4xMzdE9NMXxgQmO9U9woGeEw+fGmM9danTXxhru3t7IPTsauaO1IdKW7/RclteOD7K3s5/nO/s5dHYU9yC42jZUcf011Vx3TQ0tqQo21JSyobqM6rIkRUVGwozZzDxDk7MMTc5xeniKo30THOkb5+DpEY70TQDBa928uZZ7dwTbfNu19RSvZBi4w5/cASWV8OgzzP3fX6b4ra/xwoaf50uTH+K3Rv8DCeb5kezvs6O5gZ1N1bQ1VdHaWEW6upRUZQm1FcUUmTHvTjbrDEzM0jc2Q+/YNF294xw+N8Y7Z8c4MTAJQEmyiNu21nFfW5p7dzRy06baSFuBgxOz7O3q58Wufl440s+pwaCxUlZcxE0ba7m+uZqWVCWtjZU01ZRRW15MTXkxxQljLuvMZuYZmZrj3Og0p4eDv+23T49y6MwYZ0eDRkBlSYL2lgbua2vk/p1pdmyowiK8tu650Wm+e7iPF7r6efnoAOdGg3/cVaVJdm2s4caNNbSkKtmaqmBjbTnVZUkqS5KUFhcxk5lnei7L2PQcp4en6R6a4sTABG+fGeWtnhGGc9/s6iuKuaM1xT1tjdzflmZrqiKvms1sn7u3L7mswCF/F/Cb7v5A7v5nAdz9t5d6/NWG/Osnh/jtv32H+3em2dOWZtfGmoL/lzw1OMlznX08d7iPF7sGGJvJUGRw85Y67tyWYndrA+3X1lNdVpzXeiZnM+zvHuGVo4O8eKSf108OM5udpzRZxB3bUuxpa2TPzjRtBf5guDtH+iZ47nAf3z3cx8tHB5jJzFOcMG7bWs9d21Psbmnglq11VJTk1/oemZzjje5hXj02wAtdwT/NeQ8+ZHdvT7En93vO94NxOfPzzpnvPM6m53+N/qI0jfN9fDHzE/ypfYLdrSl+pu4dPrr/l5n72B9T3P7JvNbVPz7DvhNDvHZskL1d/bxzdgwIwuC+tjT370xzX1sjG2rKwti0i5rLzvPmqeHg99zZz/7uYdyhuizJndtS3LM9RXtLA9ddU533P9yzI9O8enyQV48N8OKRAY7m/rlfU1PGnp3B3/W9OxoL/g12ocHyfGc/zx3uO//eN1aVctf2FHdtS7G7tZ5tjfm1vt2dU4NTvHp8kJePDvDSkQF6hoN/mi2pCh7bs52fu2PrVb12lCH/U8CD7v4LufufBO5w919a9JjHgMcAtm7devuJEyeWvZ69nf389t8e4uDpUQBSlSXcsa2BO7eluKM1xY4NVXm1htyds6PTvHZ8iJeODPDy0QGO9Qd/kJvqyoM/yLY0d+9opLY8v1C/nKnZLK8cG+C7h4N/MEcWfTDu3pHijtYGdremaElV5BX67s7xgUn2nRjixSP9vNg1cL7ltS1dyZ5c8OxubSj4/oORqTleOtLPdw8HH8KFD8bWhorz32raW+rZVFee1zZn553D58bYd2KIl44O8GJXP0OTs/zrxFN8uPxdTl77k6Tu/iS3X1sfdKm4w5fuhNJq+IVvh7W5APSNzfBCV7C9z3X20T8+C8DOpqrcNqe4bWtd3qE/l53n7dOjvHZ8kJePBuEzvqjBcv/ONHt2prl5c13Bv1F0D02yt7Of5zr72NvZz+h0UMeujbXcvT3FndtT3LqlLu/Qn8lkcw2nAV45NshrxweZnpunJFFEe0s9e3YGf9vXX1Nd8IbTsf6J3O+4nwd2NfEzH4xhyC+Wb3dN39gMe7v6eP5w8DXrdK5PuKIkwa6NNezaWEtrYyWb68vZXF9BbXkx5cUJSouLyMw7EzMZxmcynBuZ5uTgJCcHJ3nn7Bj7u0foHw++slWXJrljWwN3bW/k/p2NbE9H+9Vy8QfjlaODDEwEYVBfUcwNzTVcf00NbU1V53d6pSpLKUkUUZw03IOdoyNTcwyMz3KsP+gm6uwdZ3/38Pmvlg2VJdy1PcU92xu5r62RLQ2FbUFfysI3ir2dfezt6uflo4OMz2QAaKop5YbmGq5rqmbHhio21pXTWFVKY1UJJckiinK/p7HpDP25HaMnBiboPDdOZ+8YB3pG3/Na9+xo5N7cz0XDdO8X4du/AZ/eD/XXFmSb5+edt8+Msrernxe6+nn12CAzmeCcOpvqyrlpUw07NlSxPV3FloYK6iuKqasoobIkSdad7LwzNZuldyzYAd49NMm758Y4dGbsPf3LC/8097Q1cvf2RmorCttguZRMdp43u4d57nA/Lx0d4PWTQ8xlg6xqSVXwTzbXsT1dSUuqki0NFTRUllBVmqS6LMm8OzNz80xnsvSOznB6eIqe4SnePTvGwdOjdPaOnX+t65qquWt7ij07G7lzWyrvb6FRiX13zVLcne6hKV45NsiBnhHe6hnh7dOjy9phkiwyWhsr+cDmOj6wuZZbt9axa2O0faSXshCArxwb4K3uEQ6dCXaALgTClSgrLqK1sYqbN9dyy5ZgB+hK7iRaruy8887ZUfadGOJ7J4Z499w4R/rGmV3GNleXJtnRVMWNzTW0t9Rz+9YGtjRc4beCgSPwR7fBg5+HO38xjy25cjOZLAd6Rnn95BCvnxzm0NlRTg5Mkpm/8s9yXUUx1zVVc0NzDbdfW88HWxq4prawXUH5mJrN8vrJId7sHuHNU8O81TPC6ZGpZY1ubawq4caNtdzYXMMtW+rY3dpAQ+XqHcywHFGGfJJgx+uHgR6CHa8/5+4Hl3p8oXe8zs87/RMzdA9N0TM0xfhMhqnZLFNzWUoSRVSUJqgoSdBUXcaWhgqaa8tW7V7/K5XJznNmJBi2eXZ0msGJWeayzlw2CMGasmJqy4uprywOdp5Vl63aQL9Smew8p3JDN/vHZ+gfm2Eu6ziOO1SVJUlVlpKqKmFLfQVNNaX5fRv7b7dBwzb4Z38Z3kYs01x2npODk5wenmJoco7hyVkmZrIkiiBRVERpsoimmjI2VJfSXFdGuirPbV4FpueydA8F37hHpuYYmw4GMCSLjNJkEaXFCdJVpWysK2djXbBTeK1v88VEFvK5lf8o8EWCIZRfcffPXeyxa3V0jaxzf/2r8MbX4TMnIBFdF4esX5cK+YI3U939b9x9p7tvv1TAi6xZLffC3ERwdKzIKrO2+yJEVoOtdwXTky9HW4fIEhTyIvmqvgZqt0L3a1FXIvI+CnmRMGy+HXr2RV2FyPso5EXCsPFWGDkFE/1RVyLyHgp5kTA03xJMtfNVVhmFvEgYmj8QTM/uj7YOkQso5EXCUF4PtVvg3JLH+YlERiEvEpamXXD2QNRViLyHQl4kLBtuhIFOyMxGXYnIeQp5kbA07YL5DPQfjroSkfMU8iJhSV8fTPveibYOkUUU8iJhaWwDS0DvoagrETlPIS8SlmQpNLRC/7tRVyJynkJeJEzp66FPIS+rh0JeJEyNbTB4FLJzUVciAijkRcLVeF0wwmbwWNSViAAKeZFwNe4MpgOd0dYhkqOQFwlT445gqrHyskrkFfJm9tNmdtDM5s2s/YJlnzWzLjN718weyK9MkTWirBYqN0B/V9SViACQzPP5B4CfAP7H4plmdiPwMLAL2Ah828x2uns2z/WJrH6NbequkVUjr5a8ux9y96XGiz0EfMPdZ9z9GNAF7M5nXSJrRmoH9CvkZXUoVJ/8JuDUovvduXnvY2aPmVmHmXX09fUVqByRFdTYBlODMDkYdSUilw95M/u2mR1Y4uehMApw98fdvd3d29PpdBgvKRKtVFswHVC/vETvsn3y7v6Rq3jdHmDLovubc/NE4i+1MMKmE7aol1KiVajumqeBh82s1MxagTbg1QKtS2R1qb8WipJqycuqkO8Qyh83s27gLuCvzezvANz9IPAk8DbwLeBTGlkj60aiGOpbFPKyKuQ1hNLdnwKeusiyzwGfy+f1Rdas1A6FvKwKOuJVpBBSO4ITlc3PR12JrHMKeZFCSO2AzDSMnLr8Y0UKSCEvUgiNuWGUg0eirUPWPYW8SCGcH0apfnmJlkJepBCqmqCkWuewkcgp5EUKwSw47bDOYSMRU8iLFErjTg2jlMgp5EUKpbEtGF0zOxF1JbKOKeRFCmXhUoBX0mXz5Qfg+d8rbD2yLinkRQrlfMhf5lKAmRk49QpkZgtfk6w7CnmRQmnYDpa4fMgPnwQcGlpXpCxZXxTyIoWSLIGGbdB76NKPGzoeTOtbCl2RrEMKeZFCSl8HfUtdIXORgdxRsQ3bC1+PrDsKeZFC2nBDcKKyzMzFHzN4NDhwqrJx5eqSdUMhL1JIG24Az166X36gC1LbgwOoREKmkBcppA27gum5gxd/zEDn9891IxIyhbxIIaV2QKIUzr619PLZSRg+9f3hliIhU8iLFFIiCRuuv3jID3QCDmmFvBSGQl6k0JpvhrP7wf39y3rfCabpG1a2Jlk38r2Q9xfM7B0z229mT5lZ3aJlnzWzLjN718weyL9UkTWq+RaYGsod9HSB3oOQKAl2vIoUQL4t+WeAm9z9A8Bh4LMAZnYj8DCwC3gQ+JKZJfJcl8jatPHWYNqz7/3Lzr4F6eshUbyyNcm6kVfIu/vfu3smd/dlYHPu9kPAN9x9xt2PAV3A7nzWJbJmNd0EyTLo7njvfHc48yY0fyCaumRdCLNP/l8Cf5u7vQlYfAXj7ty89zGzx8ysw8w6+vr6QixHZJVIlsCm2+HkS++dP3QcJgdg422RlCXrw2VD3sy+bWYHlvh5aNFj/j2QAf58uQW4++Pu3u7u7el0erlPF1kbrr07aLVPj3x/3qlXg+nmD0ZTk6wLycs9wN0/cqnlZvYvgI8BH3Y/P3ygB9iy6GGbc/NE1qfW++G5L8Cx5+GGjwXzTuyF0lpo2hVtbRJr+Y6ueRD4NeDH3H1y0aKngYfNrNTMWoE24NV81iWypm25A0pr4N1cj6Y7dD0LrfdBkcYkSOHk2yf/x0A18IyZvWFmfwrg7geBJ4G3gW8Bn3L3bJ7rElm7kiVw/Ufh0P8LjnI99SqM9gTzRArost01l+LuFz3hhrt/DvhcPq8vEiu3/XN48+uw738FV4IqqYIb/mnUVUnM5RXyIrIMW++CHR+Bv/v14P6Hfh1Kq6OtSWJPIS+yUszgJ/8MvvsFqEzB3Z+OuiJZBxTyIiupvB4e/C9RVyHriE5QJiISYwp5EZEYU8iLiMSYQl5EJMYU8iIiMaaQFxGJMYW8iEiMKeRFRGLMfKmLC0fEzPqAEwV46UagvwCvu5LW+jas9fph7W+D6o9eobbhWndf8oIcqyrkC8XMOty9Peo68rHWt2Gt1w9rfxtUf/Si2AZ114iIxJhCXkQkxtZLyD8edQEhWOvbsNbrh7W/Dao/eiu+DeuiT15EZL1aLy15EZF1SSEvIhJjsQ55M/tpMztoZvNm1r5ofouZTeUuPn7+AuSrzcXqzy37rJl1mdm7ZvZAVDUuh5n9ppn1LHrffzTqmq6EmT2Ye5+7zOwzUddzNczsuJm9lXvfO6Ku53LM7Ctm1mtmBxbNazCzZ8ysMzetj7LGy7nINqz4ZyDWIQ8cAH4CeG6JZUfc/Zbczy+ucF1Xasn6zexG4GFgF/Ag8CUzS6x8eVflDxa9738TdTGXk3tf/wT4EeBG4Gdz7/9a9AO5930tjDX/KsHf9mKfAZ519zbg2dz91eyrvH8bYIU/A7EOeXc/5O7vRl3H1bpE/Q8B33D3GXc/BnQBu1e2unVjN9Dl7kfdfRb4BsH7LwXk7s8BgxfMfgh4Inf7CeDjK1rUMl1kG1ZcrEP+MlrN7HUz+66Z3Rd1Mcu0CTi16H53bt5a8Etmtj/3VXZVf93OWcvv9WIO/L2Z7TOzx6Iu5io1ufuZ3O2zQFOUxeRhRT8Daz7kzezbZnZgiZ9LtbbOAFvd/VbgV4CvmVnNylT8XldZ/6p1me3578B24BaC38HvRVrs+nKvu99G0O30KTPbE3VB+fBg7PdaHP+94p+BZKFXUGju/pGreM4MMJO7vc/MjgA7gRXfIXU19QM9wJZF9zfn5kXuSrfHzP4n8FcFLicMq/a9Xg5378lNe83sKYJuqKX2Va1m58ys2d3PmFkz0Bt1Qcvl7ucWbq/UZ2DNt+SvhpmlF3ZUmtk2oA04Gm1Vy/I08LCZlZpZK0H9r0Zc02XlPpgLfpxgx/Jq9xrQZmatZlZCsMP76YhrWhYzqzSz6oXbwA+zNt77Cz0NPJK7/QjwzQhruSpRfAbWfEv+Uszsx4E/AtLAX5vZG+7+ALAH+E9mNgfMA7/o7pHvILnQxep394Nm9iTwNpABPuXu2ShrvUK/a2a3EHzNPg78q2jLuTx3z5jZLwF/BySAr7j7wYjLWq4m4Ckzg+Az/zV3/1a0JV2amX0d+BDQaGbdwG8AvwM8aWaPEpyS/BPRVXh5F9mGD630Z0CnNRARibF12V0jIrJeKORFRGJMIS8iEmMKeRGRGLFxOvQAAAAXSURBVFPIi4jEmEJeRCTGFPIiIjH2/wHPseadIxmtYAAAAABJRU5ErkJggg==\n"
          },
          "metadata": {
            "needs_background": "light"
          }
        }
      ],
      "source": [
        "n      = 2\n",
        "f5     = hermit(x[:(n+1)], y[:(n+1)], y_prime[:(n+1)])\n",
        "data   = f5.linspace(n=100, domain=[-0.85, 2.5])\n",
        "test_x = np.linspace(-5*pi, 5*pi, 100, endpoint=True)\n",
        "test_y = np.sin(test_x)\n",
        "\n",
        "plt.plot(test_x, test_y) # 25-\n",
        "plt.plot(data[0], data[1]) # 10-33\n",
        "plt.show()\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "InK1vx-bPlvi",
        "nbgrader": {
          "cell_type": "markdown",
          "checksum": "aae683a514e99dfecaa8672e1d6a4cb3",
          "grade": false,
          "grade_id": "plot03",
          "locked": true,
          "schema_version": 3,
          "solution": false,
          "task": false
        }
      },
      "source": [
        "**Plot 3:** trying to interpolate a sine curve (`np.sin()`) using first 4 nodes in `x` and `y` and their corresponding derivative in `y_prime`."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "eondZ70xPlvi",
        "nbgrader": {
          "cell_type": "code",
          "checksum": "29fc16454184db75748c0381eff20504",
          "grade": true,
          "grade_id": "plot03_ans",
          "locked": true,
          "points": 1,
          "schema_version": 3,
          "solution": false,
          "task": false
        },
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 286
        },
        "outputId": "34d1c17f-202f-438f-e11b-503b2b602e52"
      },
      "outputs": [
        {
          "output_type": "display_data",
          "data": {
            "text/plain": [
              "<Figure size 432x288 with 1 Axes>"
            ],
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAD4CAYAAAAJmJb0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3deZxcZZ3v8c+v9yTdnYV0FpaQgIBGkK1BEBCuBFlGRUS9qDOXGXHQEUYZZbxwvRLUQWEQXBEBZUQH9c4ICAqyhGFABMEOa0IICUs2kk4n3Z3eu6urfvePpyppYlXSVXVqa77v16teXVWn+pynq09961fPec5T5u6IiMjEVFXqBoiISOEo5EVEJjCFvIjIBKaQFxGZwBTyIiITWE2pGzDWzJkzff78+aVuhohIRVm6dOkWd29Jt6ysQn7+/Pm0tbWVuhkiIhXFzNZkWqbuGhGRCUwhLyIygSnkRUQmMIW8iMgEppAXEZnAFPIiIhOYQl5EZAJTyIsUSHvPELc/tZ5EQtN5S+mU1clQIhPJBbc+RduaLl7vHuTC9xxQ6ubIm5QqeZECeG1LP21rugC485nXS9waeTNTyIsUwFNrQ8CffcTerNrcR3vPUIlbJG9WCnmRAnh6bTeN9TV8/J3zAHhmXXeJWyRvVgp5kQJYsbGHhXs2c+DsRgBe6egvcYvkzUohL1IAr20dYMEeU2hqqGVWUz0vd/SVuknyJpV3yJvZPmb2kJm9YGbLzezzyftnmNkDZrYq+XN6/s0VKX/9w6Ns6Rtm3h6TAdi/pZFXFPJSIlFU8qPAF919IXAMcIGZLQQuAR509wOAB5O3RSa8NVsHANg3GfLzZ05mbedAKZskb2J5h7y7b3T3p5LXe4EVwF7AmcAtyYfdAnww322JVIJUoO87YwoAc6dOYkvfCMOj8VI2S96kIu2TN7P5wOHAE8Bsd9+YXLQJmJ3hd843szYza+vo6IiyOSIlsWnbIABzpzUAMGdq+Lm5Z7hkbZI3r8hC3swagduAi9y9Z+wyd3cg7bnd7n6ju7e6e2tLS9qvKBSpKO29w9RWGzMm1wGw59RJALzePVjKZsmbVCQhb2a1hIC/1d1vT97dbmZzk8vnApuj2JZIuWvfNsSspgaqqgzYUclv0glRUgJRjK4x4CfACne/dsyiu4Bzk9fPBe7Md1silaC9d4jZzfXbb89Nhvzr3Qp5Kb4oKvnjgL8B3mNmzyQvZwBXAqeY2SpgUfK2yITX3jPM7OaG7ben1Ncwua6aLX3qk5fiy3sWSnd/FLAMi0/Od/0ilaa9Z4jj3zLzDffNbKxXyEtJ6IxXkQgNj8bpHRplZmPdG+7fo7FOIS8loZAXiVBn/wgAezTWv+H+mY31bO0bKUWT5E1OIS8SoVSQz5jyxkpe3TVSKgp5kQhtTVbyO3fXzGyso7N/hLi+ClCKTCEvEqGtyWp9xpS/7K5JOHQNqMtGikshLxKhVHfNHjtV8tOT3Tdd/Qp5KS6FvEiEtvaPUFttNNW/cXTy9Mm1AHQNxErRLHkTU8iLRKirf4Tpk+sIJ4LvMD05j426a6TYFPIiEeoaGNke6GNNS1by3Qp5KTKFvEiEugdj2wN9rB2VvLprpLgU8iIR6h4YSRvyk+uqqauuUneNFJ1CXiRC3QOxtN01Zsa0ybV096uSl+JSyItExN3pHogxNU0lD6HLRpW8FJtCXiQig7E4I/FE2koeYOrkWroHVclLcSnkRSKSOqg6bVL6Sn7qpFp6FPJSZAp5kYikhkemO/AKIeS3KeSlyBTyIhFJBfjUSRm6a1TJSwko5EUikgrw5knpv3Bt6qRa+kfixOKJYjZL3uQU8iIR6RkcBUKYp5O6X9W8FJNCXiQiO7pr0od8qsJXv7wUk0JeJCI9QzGqDKbUZe6uAYW8FJdCXiQi2wZjNE+qparK0i5XyEspRBLyZnazmW02s2Vj7rvczDaY2TPJyxlRbEukXG0bjGXsqgGFvJRGVJX8T4HT0tz/bXc/LHm5J6JtiZSlnsEYzQ2ZQz61rHdotFhNEokm5N39EaAzinWJVKrdVfJNCnkpgUL3yV9oZs8lu3OmF3hbIiXVMzSacYw8QENtFTVVRu+QumukeAoZ8tcD+wOHARuBa9I9yMzON7M2M2vr6OgoYHNECmt33TVmRvOkWnoU8lJEBQt5d29397i7J4CbgKMzPO5Gd29199aWlpZCNUek4HqHRmlqyFzJAzQ11Ki7RoqqYCFvZnPH3DwLWJbpsSKVLhZPMBiL77KSh3DwVWe8SjHtuuwYJzP7JXASMNPM1gOLgZPM7DDAgdeAT0exLZFylKrOVclLuYkk5N39Y2nu/kkU6xapBDsmJ9t9Jf/Klr5iNEkE0BmvIpFIHUxt2k13jSp5KTaFvEgEUsHdvJvummbNKS9FppAXiUAquHdXyTfW19A/Eiee8GI0S0QhLxKF7ZX8Lk6Ggh0HZvuG1WUjxaGQF4nAePvkd8xfoy4bKQ6FvEgEepKVfGP9riv5xmQlr4OvUiwKeZEI9A7FaKyvoTrDXPIp6q6RYlPIi0Sgd2h0t1U8jJ2JUt01UhwKeZEI9A7Fdnu2K+zozlF3jRSLQl4kAn3Du5+cDHaMo1fIS7Eo5EUiEGag3PXIGtCBVyk+hbxIBMYzzTDApNpqqvXFIVJECnmRCIQ++d1X8mam+WukqBTyIhHoGRrd7bw1KY31NfRrCKUUiUJeJE/Do3FGRhPjGkIJIeR7VMlLkSjkRfLUN84vDElpaqihb1h98lIcCnmRPO34Vqjd98mnHqczXqVYFPIieRrvV/+lNNbXbK/+RQpNIS+Sp97h8c1AmdLYUKNKXopGIS+Sp2wr+SYdeJUiUsiL5CmX7pqR0QTDo/FCNksEUMiL5K1vnF8YkpJ6M+gfVshL4SnkRfKUdSWffDPQwVcphkhC3sxuNrPNZrZszH0zzOwBM1uV/Dk9im2JlJve4VEaaquorR7fyyl10lSP5q+RIoiqkv8pcNpO910CPOjuBwAPJm+LTDjhC0PG11UD+nYoKa5IQt7dHwE6d7r7TOCW5PVbgA9GsS2RctM7FBv3vDWwo5JXd40UQyH75Ge7+8bk9U3A7HQPMrPzzazNzNo6OjoK2ByRwhjvNMMpquSlmIpy4NXdHfAMy25091Z3b21paSlGc0QiNd5phlO2f3GIQl6KoJAh325mcwGSPzcXcFsiJZN1JV+v0TVSPIUM+buAc5PXzwXuLOC2REom25BvqK2ipso0E6UURVRDKH8JPA4cZGbrzew84ErgFDNbBSxK3haZcHqHYlmNrjEzGvXtUFIk4y8/dsHdP5Zh0clRrF+kXMUTTv9IPKtKHjQTpRSPzngVyUNqhEwuIa8Dr1IMCnmRPPQmz1ptzmJ0TerxquSlGBTyInnIdt6alMaGmu3z0IsUkkJeJA/ZfvVfivrkpVgU8iJ56N0+zXB2lXyTvh1KikQhL5KHVCXfPCnLSr5B3w4lxaGQF8lDT46VfHNDrb4dSopCIS+Sh57B3EJeM1FKsSjkRfLQOzRKfU0V9TXVWf1e6k1BZ71KoSnkRfLQMxTLuj8exlTyOvgqBaaQF8lDT5aTk6WkhlzqKwCl0BTyInnoGYxlfbYrjPniEHXXSIEp5EXykO00wynqk5diUciL5CHXPvlUd02vumukwBTyInnoHRrN6ku8U1IHXlXJS6Ep5EXysG0wt0q+rqaKhtoqTTcsBaeQF8nRUCzOyGiCqTmEPIQuG3XXSKEp5EVytG0wt7nkU5obaugZVCUvhaWQF8lRakqDfCp5jZOXQlPIi+RoW54h3zypVjNRSsEp5EVytL27JudKvkZ98lJwCnmRHOVdyTfUaAilFFz2A3yzZGavAb1AHBh199ZCb1OkGPLtk29uqN2+DpFCKXjIJ/0Pd99SpG2JFMW25MiYXE6GgtBdM5z84pBspyoWGS9114jkqHtwhMb6Gmqqc3sZpfryNYxSCqkYIe/A/Wa21MzO33mhmZ1vZm1m1tbR0VGE5ohEo3sgxrTJuXXVwI5uHg2jlEIqRsgf7+5HAKcDF5jZu8cudPcb3b3V3VtbWlqK0ByRaHQPjDB9cl3Ov5+q5LepX14KqOAh7+4bkj83A3cARxd6myLF0BVRJa+Ql0IqaMib2RQza0pdB94LLCvkNkWKZdtgjGn5VPKpb4dSyEsBFXp0zWzgDjNLbesX7n5vgbcpUhRdAyNMVyUvZa6gIe/urwCHFnIbIqUQT3jelfz2kB9QyEvhaAilSA56h2K4w7QcT4SCMKf8pNpqVfJSUAp5kRx09o8AMH1K7iEPoZrXEEopJIW8SA5SIb/HlPq81jN1Ui3d6q6RAlLIi+RgS18y5Btz75OH8ElAIS+FpJAXyUFUlfz0yXV0DoxE0SSRtBTyIjnY2jcMwIwp+VbydXQr5KWAFPIiOdjaP0JTQw11Nfm9hKZPrqVrIIa7R9QykTdSyIvkYGv/CDMb8+uqgdBdE0+4vgZQCkYhL5KDrX3DeXfVANsnOFOXjRSKQl4kB+09Q8xujqCST46z79IIGykQhbxIDtp7hpnV1JD3elKVfGf/cN7rEklHIS+Spb7hUfqGR5kzNf+QT/Xrb+lVd40UhkJeJEvtPUMAkXTXtDSFdXT0qZKXwlDIi2RpR8jnX8k31FbT1FBDR69CXgpDIS+SpShDHqClsV6VvBSMQl4kS2u3DmIGe02bFMn6ZjbVq5KXglHIi2RpTWc/c5obaKitjmR9LU31bFHIS4Eo5EWytHbrAPvMmBzZ+mY3NbBx25CmNpCCKPR3vJYfd9j4DKxaAg3NsODd0PJWCN9DW1k6XoJV90PTnPA3zDwAavIf8VF0fZvh9aehYSpMmgGTpodLdXnunms7BzjxwJa/XJBIwEgfxEd2XEaTPxtnhUsae0+fxGAsTmf/CHc9+zq/eGItv/3H4yP7pJCX2BBsWw/b1kLP6zDnEJhbgd/o6Q5dr8L6pbB1NSw4Aea9C6omfp1bnq+iqI0Ow6t/gJV3w8p7off1Ny5vnB3CfsGJsN+JMG1eado5HiMD8MKd8NTPYO1jb1xm1TBjP5j1Vpi1EN7xP2GP/UvTzvHofBUe+z48/e8QT9NdUT8VmufCcZ+Hd5xTFi/IbYMxNvcOs6Blyo47E3F49pfw0DehZ336X6yqhSP+F5zwRZi61xsWpT4VrOsa5KX2Prb2j5Qu4DtegkevhS0vQfc66N/8l4+Zdyy88zPw1veV7Rsx8Ri8+gisb4MNbbBhKQxs3bH84SuhaU84+ENwyIdh7mGVWeiNQ5n+hyKybT08cBm8dF+osGqnwFveAwd9BQ44Ndz36sPwSvLy/H+G39v7aDjzB9ByUGnbP9bGZ2HpLaGNwz0wY39Y9NWwgw5tg80roOPF8HPzCnjxbnj0O3Dil+Bdn4Oa/OdZiUz7cnj027DsdqiqhkM/Ft6QRgdhoAsGu2CwEwY6Yd0T8Jt/gCdvgtOuhHnvLGnTl7++DYCD95waqsNV98OSy2HzC7DXkXDMZ6CmAaproboueakN+9dTP4Onfw5H/h2c8IXwCYxQyQOs7xpgbWc/8yLsChq3kX545FvhTbd2UvhbDjoNpu4TLtP2gcY5sOo+eOIG+M9zoXlvOPpTcMS5MHlG8ducyevPwJ0XQvvzgIXX8YGnw95Hwl6tMH1fWPUALLst/C2P/wD2eAsc/GE46lPQmOZTWgWzcuoHbG1t9ba2tmhW9vyv4XdfgMQovOOjcNAZoVqvzTDszT2E5OoH4Q/XhJ3+5MvgmM+WtoIc6oG7vxDCvboeFp4JR54L+x6368qjdxPcewksvyNU9e//HuxzVPHanc7rz8B/fxNeuje84bb+HRx7ATTvmfl3Egl4/j9CkPZuhEM+Aosuh6l7F6nRb3TjIy/zjXte5NnzpjP10StgzaPh09PJl8HCD+76f9K1Bv7wLXj61hD8refB8RfRWzOdQy6/n/992lv5xZNrOGLedL57zuHF+YPcQ0Fw7yWwbR0c+nE45Wu7DrpEPPwPn/hRqJZrJsFR58F7vpL59VUMo8Pw8L+GAmJKC5x6BRxwSugGzGSgE1b8Nry+Xns0/N4Hr4cDFhWv3REws6Xu3pp2WaFD3sxOA74LVAM/dvcrMz02kpAf7IK7L4Zlv4Z93gln3QAzFmS3jt52+O3n4aXfhzA987rs1xGFDUvh1+dB9xp49z/DMf8Q+qqzsfL34fno2RCqlJMvC8ciiskdHr8OliyG+ubwdxz1qeyqv+G+8OJ97PtgVXD8RXDcRUUPlc/f/CDv33Ati+KPwuSZcNIlcOTfhtAer85XQtX87C+hdjK87zsce/cMDp83jfuWt/PZk/bni+8twqfIzlfh918Kn0ZmLYS/ugb2fVd262hfDo/9AJ79Bcx5B3zkp6XpIlzfBndeEAq1wz4RAj7b10r7C3DbeeFT2TGfDcVEhRzjKlnIm1k18BJwCrAe+DPwMXd/Id3j8w75Vx+BOz4Dfe1w4iVw/D/l3mfoDs/8IlQ4iTic+i/hY3Yx+u0SifAR8sGvho/IZ/8Y9j029/UN98J/XREqr6a58P7vwIGnRtfeXRnoDC++lfeEPtwzf5D9i2+srjWhC+6F38DsQ+CjtxQtVPpWP862n/81c6q6qT7hC3Dc56C+KfcVblkNd10Iax/nD83v41ObP8wwdVz70UP50BEF/KTiHrq/HvgKVNXASZfCOz+d3RvVzl68J3SrJeJh/zrkw9G1d1dig/Bf/wJ/+mFy3/5uqN7zWd8Di+HJG8L+dfaPwzGuMrerkMfdC3YBjgXuG3P7UuDSTI8/8sgjPSexIfd7/4/74qnu3zvCff3S3NaTTtda95++331xs/vPznLv2RjdutPp3ez+8w+F7f3y4+79W6Nb97o29x++K6z795e6x4ajW3c6a590v/bt7l/dw/3x690TiejWvfJe9yv3db9iL/fnb4tuvekkEr7y9m96bPF0X3vZ/r7q6YejW/foiPv9l7kvbvblXznET7zkJm/vGYxu/Tsb6HL/1SfCPvDvH3HftiG6dXetdf/xKWHdd/6j+3B/dOtOZ9My9+8fFbZ31+fcB7ujW/fKe92v2s/967Pcn7wp2n030/a61+X860CbZ8rhTAuiuAAfJnTRpG7/DfCDnR5zPtAGtM2bNy+3v3DpLeEf/dt/ch/uy20duxKPuz9xo/vXZ7tftcD9xd9Hvw1391UPuF99gPvXWgq3Y8WG3O/+5/B83XCi+9ZXot9GIuH+x++7f3WG+7cPdl/fFv023EOo3HRy+Ft+98Xwt0VtoMs33nC2++JmX3LZe/zWh56Jfhvu3vX0Xd69eE8fuHxO4d601i91//Yh4f/yx+8VZv8aHXF/4PLwP7nuGPf2FdFvI5EIr4+vtYTXy6ol0W/D3b1nUyjsFje73/rRcDtqw33ud30+mV8X5byasg75sZecK/l43H3N47n9bjbaV7j/8LjwD7n7YveRgWjWO9Dl/pvPhvV+/yj3jc9Hs95deeEu92/u4/6Nvd2X3R7dervWhgox9UlkoCu6dacTG05+imt2/9G7o33TSoZibPF0v+6Ki3xwOBbdutMY7Vzj8bFvWlFVwomE+59+FD5RXbPQfe0T0ax3V1Y9sKMSfuSaEP5R6N8a9qvFzeETb+/maNabSTzu/th14Q3lyvnRvlbWtbl/9/DQA3Hfl/MqUnYV8oXukz8WuNzdT03evjTZRfTNdI+PdHRNoYwOw5Kvwp+uCwerzv4JzF6Y+/peug9+exH0bQrjwU+8pHgHE7vWwK8/GcYRt34STv1GGD6Xi9GRcBzhkavD7ZMXh37eYo09fvHu0CfshINuh34s9+Mxg11hhNUTN+BTWvho52c49JhF/N/35fF/Hq94LIwkevwHYYjiKV+Fg8/O/Xkc6AyDCFbcBQeeFkaOFGu4Y+8muOfiMHpl9sFhhNfeR+a+vjWPw22fCq+VRZfDMRcUb+Rbx8pwvO/1p8L/44xv5f48xkfD/vXwVeE4wlk/Cidn5aGUB15rCAdeTwY2EA68ftzdl6d7fEWEfMqqJSFUhraFIWetn8xuLPpAZzio+9z/C28WZ14Hex1RuPZmEo/Bg1+Dx74Xxgofe0E48agui7HarzwcXsxbXgoHV0+7MoyrLrau1+C2v4f1T4bzCE66JLwgq8Z5YlFsKBxw+8M1Yejqoefw3Nu/xAduXsH1nziC0w+ZW9Dmv8Frfwz7x6bnYJ9j4PQrYc8shlX2bgojmtr+LZx/sOhyOPbC0pzws+J3cM8/hyGw7/w0vOf/ZnfAetsGePLGsI9OmwcfvjmM4y+2+GgY4fXwlTB5D/jA97MfwND5Ctx+Pqz/MxzyUTjjapg0Le+mlXoI5RnAdwhDKG929ysyPbaiQh6grwPu/GwYgtYwDRZ+IIzj3ve49MHiDltfhtf+AA99I5zwc8IX4YSLS3+y0qolYTTPpufC33LkuXDU3+86rHs2wv1fDieVTJ8Pp18NB763aE1Oyz2M5HnoG9C+DGYeGMJ+4VmZq75EHJ77jzBKo2c9vGVRONFszsH8eul6Lv7PZ3no4pNYMHNK+t8vlEQcnrk1vAn3bwlDA0++DJpmZ/6dzlfgj98NI8MSo/D2s+D4L8Ccg4vX7nSGesLf8ecfh/MiTr8qPM+ZPjnGY+FT7lO3wOol4IlwwtwZ3yr+EOCdbXwuVPWbl4fzVlLn4GQ636N/azivYOU94W+pqYe/ujbSEUglDflsVFzIQwiV1UvCyRQv3h3Oom2cEyrIhWeGSn9DW/L06qUw1B1+b8474IM/DHOBlAt3WPt4GGq54rfhvre+Dw49JwzD7F4X5jDpXhdOnOleC1gYqnr8Rbl39RRCIhG6KP77m2Hs9KyF4X/iHqZQiI+EIImPwNo/hTeEuYeFT2X7nbh9Ndc9tJqr71vJi18/rXRTDQz1hG6wP10fAmKfo5Pz+yTn+Zk8I5zws3pJOPmtqia8IRz3uXCiVjlZ92ToPtr8AmChOJj1tjD30qy3hbNrV90X3qT62kN3xmGfgMP/ujTnqmQyOhy6W9r+LRRrAHscEPadBSeGtr78UAj2dU+EN6mmPeGg08PZzhGfzKeQL5aRgfCOvey2UN3Hk9/baVXQ8rYdp1Xv3Rpul8FcLBl1r4M/3xSmUki9MQFMmRWq+6n7hNPDjzi3vOfHScRD8D18VehOAsBCWKamHGiaG154aar9r/xmGXc9+zrPLi7xJxQInwIfuRq2rNox9cNgN+FABFDXBEd9MpzIk5wyoSyNjoQgb1++YzqOravDJw8IczAdeGqY6+ctp5Tv/DgQion2ZTumR1nzGMT6dyyfc0io9A86vaDz4yjkS2GwG17+r3Ca9J6HQ31jqVuUm5GBMENk4+wwsVY5VevZcA8nulTXha60cb7Y/v5nbazdOsB9//TuAjcwR4l4+LQ42BVmuczn5KxSGh2BzpdDd9OeR4SJ6SrR6Eg4ONv5Csw/oWjHpnYV8mX8FlnhJk0LM9xVurrJMP+4Urcif2bZHUxOau8ZYs7UEs7HsjtV1aG7ppwmCMtFTV3orpn1tlK3JD81dTDvmHApE2XcXyBSepu2DTEnou9yFSkFhbzILnQPxJg+pYymaRbJkkJeJIOhWJyReIKmBvVqSuVSyItk0DsURns0K+SlginkRTLoHYoB0KiQlwqmkBfJIFXJN9XnMc+6SIkp5EUy2B7yquSlginkRTJIddc0NaiSl8qlkBfJQJW8TAQKeZEMepKVfLMqealgCnmRDFKVvEbXSCVTyItk0D88yuS6aqqrSvBFGyIRUciLZDAYizOpVHPIi0REIS+SwVAsUbovChGJiEJeJIOh0Tj1tXqJSGXTHiySwdBInIYaVfJS2RTyIhkMjcZpUCUvFU57sEgGQ7EEk+pUyUtlU8iLZDAUU3eNVD6FvEgGQ7G4RtdIxStYyJvZ5Wa2wcyeSV7OKNS2RAphKJbQ6BqpeIU+X/vb7v6tAm9DpCBUyctEoDJFJAP1yctEUOiQv9DMnjOzm81seroHmNn5ZtZmZm0dHR0Fbo7I+A2NJphUpzpIKltee7CZLTGzZWkuZwLXA/sDhwEbgWvSrcPdb3T3VndvbWlpyac5IpGJxRPEE65KXipeXn3y7r5oPI8zs5uA3+WzLZFiGorFAdQnLxWvkKNr5o65eRawrFDbEonaUCwBoDNepeIVcnTNv5rZYYADrwGfLuC2RCKVquTrVclLhStYyLv73xRq3SKFNjyq7hqZGPRZVCSN4dHQXVNXrZeIVDbtwSJpxOIOQF2NvvpPKptCXiSNWDxVyau7RiqbQl4kjViyu6a2WpW8VDaFvEgaI8lKvrZGLxGpbNqDRdLY3ievA69S4bQHi6SR6pOvVchLhdMeLJLGjpBXn7xUNoW8SBojo6rkZWLQHiySxo5x8nqJSGXTHiyShvrkZaLQHiyShvrkZaJQyIukMaJKXiYI7cEiacRGQ5+8Ql4qnfZgkTRi8QTVVUZ1lbprpLIp5EXSiMUT6o+XCUEhL5LGSDyhrhqZELQXi6QRiyc0b41MCNqLRdKIjboqeZkQtBeLpBGLJ6jVt0LJBKCQF0lDffIyUWgvFklDffIyUeS1F5vZR8xsuZklzKx1p2WXmtlqM1tpZqfm10yR4orF1ScvE0NNnr+/DPgQcMPYO81sIXAO8HZgT2CJmR3o7vE8tydSFCOjGicvE0NepYq7r3D3lWkWnQn8yt2H3f1VYDVwdD7bEikm9cnLRFGovXgvYN2Y2+uT9/0FMzvfzNrMrK2jo6NAzRHJTiye0FzyMiHstrvGzJYAc9Is+rK735lvA9z9RuBGgNbWVs93fSJROHa/PZhSn29vpkjp7XYvdvdFOax3A7DPmNt7J+8TqQhfOu2tpW6CSCQK9Xn0LuAcM6s3swXAAcCTBdqWiIhkkO8QyrPMbD1wLHC3md0H4O7Lgf8AXgDuBS7QyBoRkeLLq9PR3e8A7siw7ArginzWLyIi+dHwARGRCUwhLyIygSnkRUQmMIW8iMgEppAXEZnAzL18TjI1sw5gTZpFM4EtRW5ONsq9fVD+bSz39kH5t1Hty1+5tzFT+/Z195Z0v1BWIQmyOtIAAARBSURBVJ+JmbW5e+vuH1ka5d4+KP82lnv7oPzbqPblr9zbmEv71F0jIjKBKeRFRCawSgn5G0vdgN0o9/ZB+bex3NsH5d9GtS9/5d7GrNtXEX3yIiKSm0qp5EVEJAcKeRGRCaysQ97MrjazF83sOTO7w8ymjVl2qZmtNrOVZnZqidr3ETNbbmYJM2sdc/98Mxs0s2eSlx+VU/uSy0r+/O3MzC43sw1jnrczSt0mADM7Lfk8rTazS0rdnnTM7DUzez75vLWVQXtuNrPNZrZszH0zzOwBM1uV/Dm9zNpXVvufme1jZg+Z2QvJ1/Hnk/dn9zy6e9legPcCNcnrVwFXJa8vBJ4F6oEFwMtAdQna9zbgIOC/gdYx988HlpXB85epfWXx/KVp7+XAxaVux05tqk4+P/sBdcnnbWGp25Wmna8BM0vdjjHteTdwxNjXAfCvwCXJ65ekXs9l1L6y2v+AucARyetNwEvJ125Wz2NZV/Lufr+7jyZv/onwNYIAZwK/cvdhd38VWA0cXYL2rXD3lcXe7njton1l8fxViKOB1e7+iruPAL8iPH+yC+7+CNC5091nArckr98CfLCojRojQ/vKirtvdPenktd7gRXAXmT5PJZ1yO/kk8Dvk9f3AtaNWbY+eV85WWBmT5vZw2Z2Qqkbs5Nyfv4uTHbP3VzKj/NjlPNzNZYD95vZUjM7v9SNyWC2u29MXt8EzC5lYzIot/0PCF3AwOHAE2T5PJb86+jNbAkwJ82iL7v7ncnHfBkYBW4tZtuS295t+9LYCMxz961mdiTwGzN7u7v3lEn7SmZX7QWuB75OCKyvA9cQ3txl94539w1mNgt4wMxeTFarZcnd3czKbfx2We5/ZtYI3AZc5O49ZrZ92Xiex5KHvLsv2tVyM/tb4H3AyZ7shAI2APuMedjeyfuK3r4MvzMMDCevLzWzl4EDgcgPiOXSPor4/O1svO01s5uA3xW4OeNRsucqG+6+Iflzs5ndQehmKreQbzezue6+0czmAptL3aCx3L09db1c9j8zqyUE/K3ufnvy7qyex7LurjGz04AvAR9w94Exi+4CzjGzejNbABwAPFmKNqZjZi1mVp28vh+hfa+UtlVvUJbPX3KHTTkLWJbpsUX0Z+AAM1tgZnXAOYTnr2yY2RQza0pdJwxYKIfnbmd3Aecmr58LlNUnzXLb/yyU7D8BVrj7tWMWZfc8lvoI8m6OLq8m9Ic+k7z8aMyyLxNGPawETi9R+84i9NEOA+3Afcn7zwaWJ9v8FPD+cmpfuTx/adr7c+B54Lnkjjy31G1KtusMwsiGlwndYCVv007t248w6ufZ5H5X8jYCvyR0W8aS++B5wB7Ag8AqYAkwo8zaV1b7H3A8oevouTEZeEa2z6OmNRARmcDKurtGRETyo5AXEZnAFPIiIhOYQl5EZAJTyIuITGAKeRGRCUwhLyIygf1/5m4fGKuzFIEAAAAASUVORK5CYII=\n"
          },
          "metadata": {
            "needs_background": "light"
          }
        }
      ],
      "source": [
        "n      = 3\n",
        "f7     = hermit(x[:(n+1)], y[:(n+1)], y_prime[:(n+1)])\n",
        "data   = f7.linspace(n=50, domain=[-0.45, 3.3])\n",
        "test_x = np.linspace(-6*pi, 6*pi, 50, endpoint=True)\n",
        "test_y = np.sin(test_x)\n",
        "\n",
        "plt.plot(data[0], data[1])\n",
        "plt.plot(test_x, test_y)\n",
        "plt.show()\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "qR1mUG58Plvi",
        "nbgrader": {
          "cell_type": "markdown",
          "checksum": "d8edc2554aa6f2946081a47711e0f556",
          "grade": false,
          "grade_id": "plot04",
          "locked": true,
          "schema_version": 3,
          "solution": false,
          "task": false
        }
      },
      "source": [
        "**Plot 4:** trying to interpolate an exponential curve (`np.exp()`) using all nodes in `x` and `y` and their corresponding derivatives in `y_prime`."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "tMgsKsC9Plvi",
        "nbgrader": {
          "cell_type": "code",
          "checksum": "4978eef584205c7563997760c20241c6",
          "grade": true,
          "grade_id": "plot04_test",
          "locked": true,
          "points": 1,
          "schema_version": 3,
          "solution": false,
          "task": false
        },
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 286
        },
        "outputId": "6cb66e38-1d50-4bee-baf1-728162c7876c"
      },
      "outputs": [
        {
          "output_type": "display_data",
          "data": {
            "text/plain": [
              "<Figure size 432x288 with 1 Axes>"
            ],
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAD4CAYAAAAJmJb0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3de3Qc9X338fd3d7WSLNuybMt3GV8wGBuMbRQDARJSEzCUYqCES9pA0qSQFJ4mPUlb0vQ0aVKe0KbkQp+EhADBpAHiQgKGAA02EAjYGAPGd/Dd8lWybPmm6+5+nz92ZRYj25K1u7NafV7nzJnZ38zOfEer89Hot3Mxd0dERApTKOgCREQkexTyIiIFTCEvIlLAFPIiIgVMIS8iUsAiQReQbvDgwT5mzJigyxAR6VHefPPN3e5e2dG8vAr5MWPGsGTJkqDLEBHpUcxs89HmqbtGRKSAKeRFRAqYQl5EpIAp5EVECphCXkSkgCnkRUQKmEJeRKSAKeRFRAL2w/nv8eq63VlZt0JeRCRAB5rb+NGCtSzZtDcr61fIi4gEaPnWfbjD1NEDsrJ+hbyISIDermkAYOoohbyISMFZWtPAuMFllPcpysr6FfIiIgFxd5bWNDC1KjtH8dCFkDezKjN70cxWmdlKM/tyqv1bZrbNzJamhsvS3vN1M1tnZu+a2SXZ2AERkZ5q+75m6g60ZK0/Hrp2q+EY8FV3f8vM+gFvmtnzqXk/cPf/TF/YzCYB1wOTgRHAfDM7xd3jmShcRKSnW7ol1R+fD0fy7r7D3d9KTR8AVgMjj/GW2cCj7t7i7huBdcCM7hQrIlJI3t6yl2gkxMRh/bO2jRPqkzezMcA04PVU021mtszMHjCzilTbSKAm7W1b6eCPgpndbGZLzGxJXV3diZQjItIjLa1p4PQR/YlGsvf1aJfXbGZ9gceBr7j7fuAeYDwwFdgB3NWV9bn7ve5e7e7VlZUdPr1KRKTgtMUTLN+2j6lVFcdfuBu6FPJmVkQy4H/l7r8BcPdd7h539wTwc97vktkGVKW9fVSqTUSk13t35wFaYomsfukKXTu7xoD7gdXu/v209uFpi10FrEhNzwOuN7NiMxsLTAAWd79kEZGer/0iqGlZ/NIVunZ2zXnAZ4DlZrY01fZPwA1mNhVwYBNwC4C7rzSzucAqkmfm3Koza0REkpZuaWBQWZRRFaVZ3U6nQ97d/whYB7OeOcZ77gDuOIG6REQK2tKavUytGkCykyR7dMWriEiO7WtqY33doayeH99OIS8ikmPLtqYugsryl66gkBcRybn2K12nZOnOk+kU8iIiOba0poHxlWWUl2bnzpPpFPIiIjn0/p0ns3sRVDuFvIhIDm3d20T9odac9MeDQl5EJKdydRFUO4W8iEgOLd3SQElRiFOH9cvJ9hTyIiI5tLRmL2eMLKconJv4VciLiORIayzBiu37c3IRVDuFvIhIjqzZuZ/WWCJnZ9aAQl5EJGeW1uTuStd2CnkRkRxZuqWByn7FjCgvydk2FfIiIjnydk1DTu48mU4hLyKSAw2NrWzcnZs7T6ZTyIuI5MCSTXsBmD46d1+6gkJeRCQnFm6oJxoJMS2HX7qCQl5EJCcWrq9n+ugBlBSFc7pdhbyISJY1NLayeud+zh03OOfbVsiLiGTZ6xv34A7njh+U820r5EVEsmzh+npKikKcWVWe820r5EVEsmzRhnqqTxpIcSS3/fHQhZA3syoze9HMVpnZSjP7cqp9oJk9b2ZrU+OKVLuZ2d1mts7MlpnZ9GzthIhIvqo/2MKanQcC6aqBrh3Jx4Cvuvsk4BzgVjObBNwOLHD3CcCC1GuAS4EJqeFm4J6MVS0i0kO8vnEPAOeMy/OQd/cd7v5WavoAsBoYCcwG5qQWmwNcmZqeDTzkSYuAAWY2PGOVi4j0AIs21NMnGmbKqNz3x8MJ9smb2RhgGvA6MNTdd6Rm7QSGpqZHAjVpb9uaajtyXTeb2RIzW1JXV3ci5YiI5K2F6+upHjMwZw8JOVKXt2pmfYHHga+4+/70ee7ugHdlfe5+r7tXu3t1ZWVlV8sREclbdQdaWFt7kHMD6qqBLoa8mRWRDPhfuftvUs272rthUuPaVPs2oCrt7aNSbSIivcKiDfVAMOfHt+vK2TUG3A+sdvfvp82aB9yUmr4JeDKt/cbUWTbnAPvSunVERArewg319C2OcPqI/oHVEOnCsucBnwGWm9nSVNs/AXcCc83s88Bm4NrUvGeAy4B1QCPwuYxULCLSQyxaX8+MsQOJBNQfD10IeXf/I3C0O93P7GB5B249wbpERHq0Xfub2bD7EDfMGB1oHbriVUQkC9r744M6P76dQl5EJAsWrq+nf0mESQH2x4NCXkQkKxZuqGfG2EGEQ7l7nmtHFPIiIhm2vaGJzfWNgZ462U4hLyKSYQvXp86PD7g/HhTyIiIZt2hDPRV9ipg4rF/QpSjkRUQybeGGes4eO4hQwP3xoJAXEcmomj2NbN3bxDnjBgZdCqCQFxHJqNfW7wbg3PG5f2h3RxTyIiIZtGB1LSPKSzhlaN+gSwEU8iIiGdPcFueVtbv5k9OGkLynY/AU8iIiGbJwQz1NbXFmnjb0+AvniEJeRCRDFqzeRZ9oOC/Oj2+nkBcRyQB354XVtZx/8mBKisJBl3OYQl5EJANW7djP9n3NXJRHXTWgkBcRyYgFq2sxg09MHBJ0KR+gkBcRyYAFq3dx5qgBVPYrDrqUD1DIi4h0U+3+Zt7Zuo+LTsuvo3hQyIuIdNsLa2oB8urUyXYKeRGRbpq/upaRA0rz4q6TR1LIi4h0Q3NbnD+uq2NmHl3lmq7TIW9mD5hZrZmtSGv7lpltM7OlqeGytHlfN7N1ZvaumV2S6cJFRPLBa+t309yWyMuuGujakfyDwKwO2n/g7lNTwzMAZjYJuB6YnHrPT8wsf64OEBHJkPmraymLhvPm1sJH6nTIu/vLwJ5OLj4beNTdW9x9I7AOmHEC9YmI5K32q1wvmFBJcSQ/j2Mz0Sd/m5ktS3XnVKTaRgI1actsTbV9iJndbGZLzGxJXV1dBsoREcmNldv3s3N/MzPz8NTJdt0N+XuA8cBUYAdwV1dX4O73unu1u1dXVlZ2sxwRkdyZv3pXXl7lmq5bIe/uu9w97u4J4Oe83yWzDahKW3RUqk1EpGC8sKaWqVUDGNw3v65yTdetkDez4WkvrwLaz7yZB1xvZsVmNhaYACzuzrZERPLJrv3NLNu6L+9uSHakSGcXNLNHgAuBwWa2FfgmcKGZTQUc2ATcAuDuK81sLrAKiAG3uns8s6WLiATn/atc87erBroQ8u5+QwfN9x9j+TuAO06kKBGRfPfsip1UDSzl1KH5d5VrOl3xKiLSRbUHmvnj2jpmnzkyL69yTaeQFxHponlLt5NwuHLaiKBLOS6FvIhIFz2xdBtnjCzn5CH53VUDCnkRkS5Zu+sAK7bt56ppHV7fmXcU8iIiXfDE0m2EQ8afnZn/XTWgkBcR6bREwnni7e2cf/LgvHvM39Eo5EVEOmnJ5r1sa2jqMV01oJAXETk2d1jzDLQc4Ldvb6NPNMzFk/P7Ktd0CnkRkWOpWwOP3kDb0kf53bLtXDJ5GH2inb6ONHAKeRGRY1n1JGC8Ej6H/c0xruxBXTWgkBcRObZV82D0ucxd3crgvsWcN35Q0BV1iUJeRORodq+D2pU0nfynvLCmlivOHEEk3LNis2dVKyKSS6ufBOB/EzNojSd61Fk17RTyIiJHs2oejKzm4TVxxleWcfrI/kFX1GUKeRGRjuzdBDuW0jDmUhZv2sPV00fl/R0nO6KQFxHpyOqnAHiqrRqAK3rIbQyOpJAXEenIqnn4sCk8sAo+MqaCqoF9gq7ohCjkRUSOtG8bbF3MxiEz2bj7EH95zklBV3TCFPIiIkda8zQAP607gyH9irn09OEBF3TiFPIiIkdaNY+Wgacyd2MJnznnJKKRnhuVPbdyEZFsOFgLm1/l1aKPEo2E+PTZo4OuqFs6HfJm9oCZ1ZrZirS2gWb2vJmtTY0rUu1mZneb2TozW2Zm07NRvIhIxq15GnB+uH0Ss88cwaC+PeO+8UfTlSP5B4FZR7TdDixw9wnAgtRrgEuBCanhZuCe7pUpIpIjq56koXQ0y9pG8LnzxgZdTbd1OuTd/WVgzxHNs4E5qek5wJVp7Q950iJggJn13G8uRKR3aNyDb3yFp1qrOWfcICaN6HlXuB6pu33yQ919R2p6J9B+J/2RQE3acltTbSIi+WvN05jH+XXj9II4iocMfvHq7g54V99nZjeb2RIzW1JXV5epckREusYd3rifreEq9g2YxEWn9ZynPx1Ld0N+V3s3TGpcm2rfBlSlLTcq1fYh7n6vu1e7e3VlZWU3yxEROUFbl8COpfy0+SJu+uhYwqGed5+ajnQ35OcBN6WmbwKeTGu/MXWWzTnAvrRuHRGR/LP4XppCZTwXvpBPVVcdf/keotMPKjSzR4ALgcFmthX4JnAnMNfMPg9sBq5NLf4McBmwDmgEPpfBmkVEMuvALnzlb5nbNpPLqk+mvLQo6IoyptMh7+43HGXWzA6WdeDWEy1KRCSn3pqDJdp4MPZJ7vvomKCryaie88hxEZFsiLfhbzzAQpvKSadMYXxl36Aryijd1kBEerfVT2EHd3Bfy0V88ePjg64m43QkLyK9WmzRz9jJEBLjL+KccYOCLifjdCQvIr3XzuVEti7iwbZP8rVZk4KuJit0JC8ivVbTqz8Fj7J/4rWcPrI86HKyQkfyItI7Ne4hvOJ/mJc4jy9d+pGgq8kahbyI9EoNrz1I1FvYOfFGxg4uC7qcrFHIi0jv09YEi37CGz6R6y6/LOhqskohLyK9Tt38HzIgVse7p/0tw8pLgi4nqxTyItK7HKqnbPF/8SJncfkVnwq6mqxTyItIr7Lrd9+hONHI9rP+kQF9okGXk3UKeRHpNbx+A4NW/ZInQzO58uIP3XarICnkRaTX2PnEN2j1MLEL/pGy4t5xmZBCXkR6hYPrX2d4zTM8UXoVV328OuhycqZ3/CkTkd7NnV2P/wMt3p+p1/0LReHec3zbe/ZURHqtFS/NZXzjUt4cewuTxo4MupycUsiLSEE71NhE2cvfocZG8LHrvxZ0OTmnkBeRgvbGQ7cz1mto/sS3KCkp7AufOqKQF5GCtXrx81ywYw5vD7yMCR+7LuhyAqGQF5GC1HywgfJnb2VXaAinfPbHQZcTGIW8iBSk9+bcytBELXUX3U1Z/4FBlxMYhbyIFJxNLz/MlLqneXHIjZx53qygywlURs6TN7NNwAEgDsTcvdrMBgK/BsYAm4Br3X1vJrYnInI0e3ZsouKFv2eVncxHbroz6HICl8kj+U+4+1R3b7+U7HZggbtPABakXouIZE1rW4yaX3yOIm8jcs19lPftE3RJgctmd81sYE5qeg5wZRa3JSLCH+7/R85sfYv3pn2dUyZPC7qcvJCpkHfg92b2ppndnGob6u47UtM7gaEdvdHMbjazJWa2pK6uLkPliEhvs/CxH/DJnfexYvClTJ39laDLyRuZunfN+e6+zcyGAM+b2Zr0me7uZuYdvdHd7wXuBaiuru5wGRGRY3n3D3OZsfxfWV5azaRbHgKzoEvKGxk5knf3balxLfBbYAawy8yGA6TGtZnYlohIutpVLzP6xVtZGx7PSX/zGOGiwn8QSFd0O+TNrMzM+rVPAxcDK4B5wE2pxW4CnuzutkRE0jVtX0nJ/3yaWgZSctPj9O9fEXRJeScT3TVDgd9a8t+jCPCwuz9nZm8Ac83s88Bm4NoMbEtEBICWPTU03n8liUSI7Vf8inNPGhN0SXmp2yHv7huAMztorwd6x/O1RCSnWmrXs+feK+gbO8AfznuQy8/qPQ8B6So9NEREepTmjYto+eW1lMTjvH7uz7j84t59RevxKORFpMdoWj6P0ONfYJ/3Z/Wf/IJLPn5B0CXlPYW8iPQIzX/8McXzv8GyxDh2Xj6HWTPOCLqkHkEhLyL5Ld5G8zPfoOTNn/F8ohq/+ufMmjou6Kp6DIW8iOSvundpmfsFSuqW8VB8FiOv/wEzJ48IuqoeRbcaFpH8k0jAonuI33MBjbUb+Zp9lXE3/lgBfwJ0JC8i+WXfVhK//RtCm/7AS/FpPDz0a3znL2cyYkBp0JX1SAp5EckPsVZ48xckFvwbLa2t/GvbF+j/0c/z01kTKQqr0+FEKeRFJFjusOoJmP+vsHcjizmdb3MzX/n0JVw8eVjQ1fV4CnkRCc7mhfD7f4ZtS9haNJZ/bv0H9gz/GD/7i7OoGqgHfmSCQl5EcssdNrwIC38C657nQFEl/zd2C08lLuRLnzyFL1wwluJIOOgqC4ZCXkRyo60Jlv0aFt0DdWtojg7igdCnufvAJ7ls2jjmz5rIsPKSoKssOAp5Ecmu2jWw7FF4cw407eHAgNO4v99X+UndFE4bNZiHPzuZ6aN1i+BsUciLSOY11MCKx2H5Y7BrOW4haio/zve5iCd2jmFY/1LuuOYU/nz6KEIhPcUpmxTyItJ97lC3BtbNhzXPwJbXAGgZOp2Xqv6O7245jU1b+jK1agB3zxrLpacP02mROaKQF5ET09QAm16Btc/DugWwfysAbYMmsvLkW5lzsJonNxcTMuPSM4bz/fPGqFsmAAp5ETk+d9izAWoWQ82i5Lh2NeBQ3J/GUeezZPhn+eXuCTy/rQi2wbjKMm77xHBuOHs0w8t1tWpQFPIi8kHxGNSvhZ3L04Zl0FgPgBf3p2noWaw/dSavtp7Cb+pG8N7KFgBOH9mfr108jFmnD+PkIf2C3AtJUciL9Ebu0NyQPDrfvQ52v5cM9t3roH4dxJOh7eEorQNPZfeQj7O2aCIvNI7jdzv6U/9eDIABfYo4a3Q51549iEsmD9MFTHlIIS9SiNqa4cAOOLATDmxPjhtqoGELNGxOjlv2H17cLUxT2Sj2lp7EjqFTWRkfzWuHRvDK3goaa97/gnTs4DI+cVoF1SdVUD2mgnGD++rsmDynkBfJd4kEtB6Apr3JoXHPB6cP1ZE4VEfiQC1+qI7QoTrCLQ0fWk1LqA97osOoDQ1le/hCNhUPZnXzQFa2DmWLD6WtKRkHZjB6YB9OruzLZyb35eTKvkwY2o/xlWX0KynK9d5LN2U95M1sFvAjIAzc5+53ZnubIkHweIyW5sbk0NRIW8sh2pobaWtuJNbaRLzlIN5yCG89RKK1EVobsbZDhNoOEmo7RKTtEJHYISKxRqLxAxTHD1EaP0hxopEQftTt7vW+7PF+7Kaceq+g3k9il1ewiwp2+kB2eQU7vYKmcF8GhIupKC5iQJ8og8qiDCsv4bryEoaVlzKivIRh5SUM6VdCNKLTGwtFVkPezMLAj4FPAluBN8xsnruvyuZ2u8UdNr8GFjpi4P1pLK3d0trap0lOt7cddcwH2z70Pt5f5vDkB/81diDhRjzhJDw5xN1JJJKlhDHCISNkHB5/YA3uaWs6Spv7h6c7NQY8kdaWmvbEh18fHqfmJeLvv/YEeDy5TCKemk68P51IzT/8Om2ciEMiRiIRJxZrJdbWRizWRqytjXgsOR2PtZGItRKPteHxNjzWSiLeBvFWLN4GiTYs3kYo0UYo0UrYY4S9lYi3EfEYUVqJehtFFqcE6MqF+Qk3DlHCIUpo9BIarTQ1DKIxdBLN4TKaov1oCZfREulPU1E5sWg5seIKEsUDSJRWUFxcTEkkTJ9omL4lEcqLI4wqjlBWHKGsOEz/kiIqyqKURcOYqWult8n2kfwMYJ27bwAws0eB2UB+h/yDlwVdRacZyX+RdDunYwsB0dQQdyNGmDhhYoSIESZGhDbCxD1EGxHiRIhbmLgVEbciPBQhYX1JhKMkQlE8XISHoxCO4uESiBRDUQkWKSYUKYFoCaGiPoSjfQhFSwkXlxIpLiNc0pdIcV+KSvtSVFpGtLiM4miEgeEQQ8OmEJaMy3bIjwRq0l5vBc5OX8DMbgZuBhg9enSWy+kEM7jxyY6PMDs8Am1vSz+KTT967Xh8oCXGrn1N7NzXzK59Teza30ztgSZiccfSjqrDBn2LI/QtDlMcCRGNhChODUXhENGwETbDQqSO0pNBEbLUga87jpPw5FGjuxNPOHGHeCJBLAHxhBNL+BHT77fF4qnl3InFnZhDWzw5L7n3yWByDPdk9Y6RwFLzkm0JQofb2vfSMeKH20NHtIVIeHI9CUKEwmFCoRAWihAKh4mEI0QiYcLhCOGiIiLhCEWRCEVFEaJFUYqKIhQXRYlGi4gWRYlGo5QURymNFlEaDVNalDz6LYkmx6VFYUqjYQZGI5REQkR0RaYUgMC/eHX3e4F7Aaqrq4/e8ZgrZjDuwoyvdu+hVl5dv5tX3tvNK2vr2L6v+fC8QWVRJg7vx6kT+zOxsoxh/UsY2r+EoeXFDCorJpzHZy8kDv8xSCTHcScWT6T+qCS7j/zwdHJscPiINTmd/OMUTv3BCoc+OETDIYp0lCtyQrId8tuAqrTXo1JtBc/deWtLAy+9W8vLa3ezbGsD7tCvJML5Jw/mr86vYOKw/pw6rB+V/YqDLveEhUJGNGRE9Ux4kbyU7ZB/A5hgZmNJhvv1wKezvM1A7T3UyuNvbeWRxVtYX3eIcMiYWjWAL8+cwMdOqWTKyHJ1A4hIzmQ15N09Zma3Af9L8rvBB9x9ZTa3GQR3Z/HGPTy8eAvPLt9JazzB9NED+N41U7h48jDKS3VusYgEI+t98u7+DPBMtrcTBHfnqWU7+NH891hfd4h+JRFumFHFDWePZuKw/kGXJyIS/BevPdWSTXv4zu9W805NAxOH9eN710zh8ikjKI3qZEYRyR8K+S7atPsQ//7cGp5dsZOh/Yv53jVTuHr6qLw+A0ZEei+FfCfta2rjR/PX8stFmygKh/i7i07hrz82lj5R/QhFJH8poTphaU0Dtz38FtsbmrjuI1X83UWnMKS/niovIvlPIX8M7s4vXt3Ed59dzZB+JTz2pY/q8WUi0qMo5I9iX2Mbf//YO/x+1S4uOm0o//mpKQzoEw26LBGRLlHId+CdmgZuffgtdu5r5p//9DQ+f/5YXVIvIj2SQv4Ijyzewr88uYIh/Ur4ny+eyzR1z4hID6aQT/PAHzfy7adXceGplfzwuqnqnhGRHk8hn/LzlzdwxzOrufT0Ydx9wzSKdH8ZESkACnngJy+t4z+ee5c/nTKcH143VQEvIgWj14f8fy1Yy13Pv8fsqSO461Nn6g6RIlJQem3Iuzs/nL+WHy1Yy9XTR/K9a87UrQlEpOD02pD/wfPvcfcL67i2ehTfvXqKAl5EClKvDPnfvr2Vu19Yx3XVVXz36jMIKeBFpED1ug7oldv38fXfLOfssQP5t6tOV8CLSEHrVSHf0NjKF//7TQaURvnxX0zXWTQiUvB6TXdNPOH8n0feZte+FuZ+8VwG9+25D88WEemsXhPyd/3+XV5Zu5s7rz6DqVUDgi5HRCQnekV/xXMrdvCTl9Zzw4zRXD9jdNDliIjkTMGH/NpdB/jq3HeYNnoA37piUtDliIjkVLdC3sy+ZWbbzGxpargsbd7XzWydmb1rZpd0v9Sua2yNccsv36Q0GuGevziL4ogesi0ivUsm+uR/4O7/md5gZpOA64HJwAhgvpmd4u7xDGyv84U9/x4bdh/ikb8+h2HlelyfiPQ+2equmQ086u4t7r4RWAfMyNK2OrRi2z4eeHUTN8wYzbnjB+Vy0yIieSMTIX+bmS0zswfMrP0JGyOBmrRltqbaciKecL7+m+VU9Ily+6yJudqsiEjeOW7Im9l8M1vRwTAbuAcYD0wFdgB3dbUAM7vZzJaY2ZK6urou70BHHnxtE8u37eObfzaJ8j5FGVmniEhPdNw+eXe/qDMrMrOfA0+nXm4DqtJmj0q1dbT+e4F7Aaqrq70z2zqWbQ1N3PX7d/nEqZVcPmV4d1cnItKjdffsmvQUvQpYkZqeB1xvZsVmNhaYACzuzrY6w9355pMrcIdvzz5dD98WkV6vu2fX/IeZTQUc2ATcAuDuK81sLrAKiAG35uLMmudW7GT+6lq+cdlpVA3sk+3NiYjkvW6FvLt/5hjz7gDu6M76u2J/cxvfnLeSySP687nzxuRqsyIiea1g7l3zvefeZffBFu67qVqP8BMRSSmINHxz817++/XN3PTRMUwZpZuPiYi0K4iQj4SMCyZU8tWLTw26FBGRvFIQ3TVnVg3gob/K6QW1IiI9QkEcyYuISMcU8iIiBUwhLyJSwBTyIiIFTCEvIlLAFPIiIgVMIS8iUsAU8iIiBczcu30L94wxszpgc9B1pAwGdgddRDeo/uD19H1Q/cHr7D6c5O6VHc3Iq5DPJ2a2xN2rg67jRKn+4PX0fVD9wcvEPqi7RkSkgCnkRUQKmEL+6O4NuoBuUv3B6+n7oPqD1+19UJ+8iEgB05G8iEgBU8iLiBQwhXyKmQ00s+fNbG1qXHGU5eJmtjQ1zMt1nR3UM8vM3jWzdWZ2ewfzi83s16n5r5vZmNxXeXSdqP+zZlaX9jP/QhB1Ho2ZPWBmtWa24ijzzczuTu3fMjObnusaj6UT9V9oZvvSfv7/kusaj8XMqszsRTNbZWYrzezLHSyTt59BJ+vv3mfg7hqS30v8B3B7avp24N+PstzBoGtNqyUMrAfGAVHgHWDSEcv8DfDT1PT1wK+DrruL9X8W+H9B13qMffgYMB1YcZT5lwHPAgacA7wedM1drP9C4Omg6zxG/cOB6anpfsB7HfwO5e1n0Mn6u/UZ6Ej+fbOBOanpOcCVAdbSWTOAde6+wd1bgUdJ7ke69P16DJhpZpbDGo+lM/XnNXd/GdhzjEVmAw950iJggJkNz011x9eJ+vOau+9w97dS0weA1cDIIxbL28+gk/V3i0L+fUPdfUdqeicw9CjLlZjZEjNbZGZB/yEYCdSkvd7Kh39BDi/j7jFgHzAoJ9UdX2fqB/jz1L/Zj5lZVW5Ky5jO7mM+O9fM3jGzZ81sctDFHE2qK3Ia8PoRs3rEZ3CM+qEbn0FBPMi7s8xsPjCsg1nfSH/h7m5mRzu39CR332Zm44AXzGy5u6/PdK1y2FPAIwLkhd8AAAHxSURBVO7eYma3kPyv5E8Crqk3eYvk7/xBM7sMeAKYEHBNH2JmfYHHga+4+/6g6+mq49Tfrc+gVx3Ju/tF7n56B8OTwK72f+FS49qjrGNbarwBeInkX96gbAPSj2xHpdo6XMbMIkA5UJ+T6o7vuPW7e727t6Re3geclaPaMqUzn1Hecvf97n4wNf0MUGRmgwMu6wPMrIhkQP7K3X/TwSJ5/Rkcr/7ufga9KuSPYx5wU2r6JuDJIxcwswozK05NDwbOA1blrMIPewOYYGZjzSxK8ovVI8/4Sd+va4AXPPVtTh44bv1H9J1eQbLPsieZB9yYOsPjHGBfWrdg3jOzYe3f4ZjZDJKZkS8HCaRqux9Y7e7fP8piefsZdKb+7n4Gvaq75jjuBOaa2edJ3u74WgAzqwa+6O5fAE4DfmZmCZI/6DvdPbCQd/eYmd0G/C/JM1UecPeVZvZtYIm7zyP5C/RLM1tH8gu264Oq90idrP9vzewKIEay/s8GVnAHzOwRkmc/DDazrcA3gSIAd/8p8AzJszvWAY3A54KptGOdqP8a4EtmFgOagOvz6CABkgdanwGWm9nSVNs/AaOhR3wGnam/W5+BbmsgIlLA1F0jIlLAFPIiIgVMIS8iUsAU8iIiBUwhLyJSwBTyIiIFTCEvIlLA/j8Kxk7StdbmZQAAAABJRU5ErkJggg==\n"
          },
          "metadata": {
            "needs_background": "light"
          }
        }
      ],
      "source": [
        "#defining new set of given node information: x, y and y'\n",
        "x       = np.array([0.0, 1.0,          2.0       ])\n",
        "y       = np.array([1.0, 2.71828183,  54.59815003])\n",
        "y_prime = np.array([0.0, 5.43656366, 218.39260013])\n",
        "\n",
        "\n",
        "f7      = hermit( x, y, y_prime)\n",
        "data    = f7.linspace(n=50, domain=[-0.7, 2.45])\n",
        "test_x  = np.linspace(-0.5, 2.2, 50, endpoint=True)\n",
        "test_y  = np.exp(test_x**2)\n",
        "\n",
        "plt.plot(data[0], data[1])\n",
        "plt.plot(test_x, test_y)\n",
        "plt.show()\n",
        "\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "DYVL-_XcPlvi",
        "nbgrader": {
          "cell_type": "markdown",
          "checksum": "5419db887410ebe73b42fa5f21ffc50e",
          "grade": false,
          "grade_id": "plot05",
          "locked": true,
          "schema_version": 3,
          "solution": false,
          "task": false
        }
      },
      "source": [
        "**Plot 5:** trying to interpolate $y = (x-3)^2 + 1$ using all nodes in `x` and `y` and their corresponding derivatives in `y_prime`.\n",
        "\n",
        "For this plot you might be able to see only one curve due to the two curves overlapping. This means that our polynomial is accurately interpolating the original function.\n"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "tuHX_czoPlvj",
        "nbgrader": {
          "cell_type": "code",
          "checksum": "92cc6293f29ebf42544cedf49e4c1fc5",
          "grade": true,
          "grade_id": "plot05_test",
          "locked": true,
          "points": 3,
          "schema_version": 3,
          "solution": false,
          "task": false
        },
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 286
        },
        "outputId": "49276e20-a484-4904-cb9d-49ee9f9af943"
      },
      "outputs": [
        {
          "output_type": "display_data",
          "data": {
            "text/plain": [
              "<Figure size 432x288 with 1 Axes>"
            ],
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3deXhU5f3+8fcnGwlhCxB2whr2JUCAqLiCCm6grai17i1ad639qrX91VZr3fdWK4q44S6i1Yq4VFzYEhJD2AOyhSUJS4CQhSTP7w/GNmKQkMzkTGbu13XlysxzZubczBVuhifnnMecc4iISGiJ8DqAiIj4n8pdRCQEqdxFREKQyl1EJASp3EVEQlCU1wEA2rZt67p37+51DBGRRiUjI6PQOZdY07agKPfu3buTnp7udQwRkUbFzNYfapumZUREQpDKXUQkBB223M2sq5l9bmbLzGypmd3gG29tZnPMbLXve4Jv3MzscTPLNbNsMxse6D+EiIj8UG0+uVcAv3XODQDSgGvMbABwG/Cpcy4Z+NR3H2ACkOz7mgI85ffUIiLykw5b7s65Lc65xb7be4DlQGdgIvCC72EvAJN8tycCL7oD5gOtzKyj35OLiMghHdGcu5l1B4YBC4D2zrktvk1bgfa+252BjdWetsk3dvBrTTGzdDNLLygoOMLYIiLyU2pd7mbWDHgbuNE5t7v6Nnfg0pJHdHlJ59wzzrlU51xqYmKNh2mKiEgd1arczSyaA8X+inPuHd/wtu+nW3zf833jeUDXak/v4hvzu/w9pfzl/WUU7dsfiJcXEWm0anO0jAHPAcudcw9X2/QecInv9iXArGrjF/uOmkkDiqpN3/jVjuJypn39HVO/XBuIlxcRabRq88n9GOAi4CQzy/J9nQbcC5xsZquBcb77AB8Ca4FcYCpwtf9jH9CvQwtOH9KR57/+jh3F5YHajYhIo3PYyw84574C7BCbx9bweAdcU89ctXbj2GQ+XLKFf85dw+0T+jfUbkVEglqjP0M1uX1zJg7txIvfrKdgT5nXcUREgkKjL3eA68cmU1ZRyT+/WON1FBGRoBAS5d4zsRnnDO/CS/PXs213qddxREQ8FxLlDnD9SclUVDme+o8+vYuIhEy5J7VpyrkjujBjwQY27yrxOo6IiKdCptwBrj2pNw7H3z/P9TqKiIinQqrcuyQ05byRXXkjfSMbd+zzOo6IiGdCqtwBrjmxN2bGk5/p07uIhK+QK/eOLeP4xagk3lq8ifXbi72OIyLiiZArd4CrT+hFVITx6CervY4iIuKJkCz3di1iufTo7ryblcfKrXu8jiMi0uBCstwBrjq+F81ionjw45VeRxERaXAhW+4J8TFMOa4nc5ZtI3PDTq/jiIg0qJAtd4DLxvSgTXwMD8zWp3cRCS8hXe7NmkRxzYm9+WbNdr5aXeh1HBGRBhPS5Q5wYVoSnVrG8sDsFRy41LyISOirzTJ708ws38xyqo29Xm1VpnVmluUb725mJdW2PR3I8LXRJCqSG8f14dtNRcxeus3rOCIiDaI2n9ynA+OrDzjnznPOpTjnUjiwcPY71Tav+X6bc+4q/0Wtu3OGd6ZXYjwPfbySyip9eheR0HfYcnfOzQV21LTNt3j2ZOBVP+fyq6jICH57Sl9W5+/l3cw8r+OIiARcfefcjwW2Oeeqnwraw8wyzewLMzv2UE80sylmlm5m6QUFBfWMcXjjB3ZgUOcWPPLJKsorqgK+PxERL9W33C/gh5/atwBJzrlhwM3ADDNrUdMTnXPPOOdSnXOpiYmJ9YxxeBERxu9O7cemnSW8unBDwPcnIuKlOpe7mUUB5wCvfz/mnCtzzm333c4A1gB96hvSX45LbsvoHq154rPV7C2r8DqOiEjA1OeT+zhghXNu0/cDZpZoZpG+2z2BZGBt/SL6j5lx24R+FO4tZ+rcoIklIuJ3tTkU8lVgHtDXzDaZ2RW+Tefz41+kHgdk+w6NfAu4yjlX4y9jvTIsKYHTBndg6ZczKdi26fBPEBFphCwYTuxJTU116enpDba/jWtX0OmFNBYk/oyjr32uwfYrIuJPZpbhnEutaVvIn6Fak649+7G47VmMLJjJ+tVLvI4jIuJ3YVnuAD0n381+oih89w6vo4iI+F3Ylnub9klkJ13MiOIvWJ7+mddxRET8KmzLHWDw5D9QSCuY/QdclU5sEpHQEdblHt+8FWsGXkf//UvJ+mSG13FERPwmrMsdYMSk61kf0YU28++horzM6zgiIn4R9uUeFR1DYdrvSarKI3PWE17HERHxi7Avd4Dh4y5gWfQgei59nOLdWm9VRBo/lTtgERHYqXfThiKy37zb6zgiIvWmcvfpn3oiGc1OYOiGl9iW953XcURE6kXlXk3Hn91LFJWse/N2r6OIiNSLyr2aTj36k9npPEbu/IjVWV96HUdEpM5U7gcZcN5dFFlzyj+8XSc2iUijpXI/SPNWbVg94DoGli8hc87LXscREakTlXsNhp99I+sikmg3/6+Ule7zOo6IyBFTudcgKjqG3cfdSRe3lay37vc6jojIEavNSkzTzCzfzHKqjd1pZnlmluX7Oq3attvNLNfMVprZqYEKHmhDTvgZ38aOZEDu0+ws2Ox1HBGRI1KbT+7TgfE1jD/inEvxfX0IYGYDOLD83kDfc/7x/ZqqjVHLSfcT58rIff33XkcRETkihy1359xcoLbroE4EXnPOlTnnvgNygVH1yOep7v2Gk544ieEF77JhRcMtAygiUl/1mXO/1syyfdM2Cb6xzsDGao/Z5BtrtPpM/ivFFseud2/VoZEi0mjUtdyfAnoBKcAW4KEjfQEzm2Jm6WaWXlBQUMcYgde6XSeWJ/+GIaXpZH32htdxRERqpU7l7pzb5pyrdM5VAVP539RLHtC12kO7+MZqeo1nnHOpzrnUxMTEusRoMMPP/T82RHQh8es/UVqiQyNFJPjVqdzNrGO1u2cD3x9J8x5wvpk1MbMeQDKwsH4RvRcdE8vuE+6mi9tK5ht/9TqOiMhh1eZQyFeBeUBfM9tkZlcA95vZEjPLBk4EbgJwzi0F3gCWAR8B1zjnKgOWvgENOu5sspoezdC1U8nXVSNFJMiZc87rDKSmprr09OA/GmXzd8tpM/1YlrQ8ntSb3/Y6joiEOTPLcM6l1rRNZ6gegU49+pPZ9Zek7v6EFQs+9jqOiMghqdyP0JDz/0w+rYn6+DYqKyq8jiMiUiOV+xFq2qwlG1J/T+/KNWS8+7jXcUREaqRyr4MRp13BsuhBJOc8QtHO4D1GX0TCl8q9DiwigpizHqSF28PyGbd5HUdE5EdU7nXUe/BRpCeezcj8t8nN/sbrOCIiP6Byr4d+F95PkTWn4r2bqaoMicP5RSREqNzroWVCImtSbqVfxXIWvfuk13FERP5L5V5PqWddzfLogfRZ8gC7Crd6HUdEBFC515tFRBA76RGau2JWzvid13FERACVu1/0GDia9A7nMXL7+6zM+NzrOCIiKnd/GXTh3yi0BCI//K3OXBURz6nc/aRZiwQ2jrqD3pVrSH/7Qa/jiEiYU7n70fDxl7OkyTD6L3+Mwq0bD/8EEZEAUbn7kUVE0PJnjxLrylg34yav44hIGFO5+1lSnxQWd72U1N1zWDJ3ptdxRCRMqdwDIOXCu9honUj4/DZK9+31Oo6IhKHaLLM3zczyzSyn2tgDZrbCzLLNbKaZtfKNdzezEjPL8n09HcjwwSo2Lp6icfcfWHP15d97HUdEwlBtPrlPB8YfNDYHGOScGwKsAm6vtm2Ncy7F93WVf2I2PoOOOZNFrcaTmvcy3y1d4HUcEQkzhy1359xcYMdBYx87574/mHs+0CUA2Rq95AsfZa81pWzmDbqwmIg0KH/MuV8O/Lva/R5mlmlmX5jZsYd6kplNMbN0M0svKAjNBS9aJXZkzbDbD1xY7K2HvI4jImGkXuVuZncAFcArvqEtQJJzbhhwMzDDzFrU9Fzn3DPOuVTnXGpiYmJ9YgS1EWf+hpwmwxiw7GHy89Z5HUdEwkSdy93MLgXOAC50zjkA51yZc26773YGsAbo44ecjZZFRJAw+QliqGDjjOu8jiMiYaJO5W5m44H/A85yzu2rNp5oZpG+2z2BZGCtP4I2Zp17DSazx68ZUTyXrDkzvI4jImGgNodCvgrMA/qa2SYzuwJ4EmgOzDnokMfjgGwzywLeAq5yzu2o8YXDzIhf/Im1Ed3p9PUd7N5V6HUcEQlx5ptR8VRqaqpLT0/3OkbArVz8Bb1nTSSjzRmMuv5lr+OISCNnZhnOudSatukM1QbUd/jxLOr4C0bteJ+cL2d5HUdEQpjKvYGlXHw/G60TrT/7HcV7iryOIyIhSuXewGKbNmPvKQ/TyW1jyYu3eB1HREKUyt0D/Y+awMK25zAq/02WL5zjdRwRCUEqd48MvPhh8q0tTf99I6UlxV7HEZEQo3L3SHyLBApPvJ9ubhOLX9KVI0XEv1TuHhp0/DmktxrPqLwXWZX5pddxRCSEqNw9lnzJk+ywVsS8fzVlpfsO/wQRkVpQuXusZUIim4+7j+5VG8h88Vav44hIiFC5B4GUkyazMOEMRua9xAodPSMifqByDxL9L32CfGtL/L+vZ1/xbq/jiEgjp3IPEs1btmb7yY/Q1W3m2+k3ex1HRBo5lXsQGXTMmSxody5HFbzJki/f9zqOiDRiKvcgM/TSR9honWj76U3sLtLVkkWkblTuQSa2aXNKTn+Cdq6QZc9r5SYRqRuVexDqkzqOjC4XkbbrX2TOeeXwTxAROUityt3MpplZvpnlVBtrbWZzzGy173uCb9zM7HEzyzWzbDMbHqjwoSzl4gdYE9mT7l/fxvZtG7yOIyKNTG0/uU8Hxh80dhvwqXMuGfjUdx9gAgfWTk0GpgBP1T9m+IlpEkvEz58jzpWwefrluKoqryOJSCNSq3J3zs0FDv7t3kTgBd/tF4BJ1cZfdAfMB1qZWUd/hA03PfoPJ7P/bxlcsoiFb97vdRwRaUTqM+fe3jm3xXd7K9Ded7szsLHa4zb5xn7AzKaYWbqZpRcUFNQjRmhLm3wrWbGjGLrsQdYtz/A6jog0En75hao7sMr2Ea207Zx7xjmX6pxLTUxM9EeMkGQREXS5dBr7LI6Kt35FaYkuLiYih1efct/2/XSL73u+bzwP6FrtcV18Y1JHbTt0ZeOx99O7ci2Lp2tpPhE5vPqU+3vAJb7blwCzqo1f7DtqJg0oqjZ9I3U0dOwFLGozkbStM8ia+57XcUQkyNX2UMhXgXlAXzPbZGZXAPcCJ5vZamCc7z7Ah8BaIBeYClzt99RhavDlT5IX2YlOn93AjvzNXscRkSBmB6bLvZWamurS09O9jtEorF0yjy5vncHypiMY8ruPsAidhyYSrswswzmXWtM2NUMj03PwUWT0u4WhJQtIf/2vXscRkSClcm+ERk++lYy4Yxi64hFys+Z6HUdEgpDKvRGKiIyg5xXPs90SiJv1K/bo6pEichCVeyOV0LY9u057mvZVBax89gpdnkBEfkDl3oj1H3Uy6T1/Q+qez5j/zmNexxGRIKJyb+RG/fIucmKHk7LkHnJzFnodR0SChMq9kYuIjKTTZS+wz5oS/fZl7N290+tIIhIEVO4hoHX7JLad8g+6VOWxcuplmn8XEZV7qBhw9Oks6HENI/Z8Tvqb93kdR0Q8pnIPIWkX/YXFcWkMXfYAqzM+8zqOiHhI5R5CIiIj6fmrlyiIaEuL93/FjnxdjFMkXKncQ0yrNu3YN/F5Wrnd5D33Syr27/c6koh4QOUegpJTjiF7yB0MLlvMgum3eh1HRDygcg9RI8+5kYyECRy1aRoZc17zOo6INDCVe6gyY9CUqXwX3ZM+X93I+pVZXicSkQakcg9hTeKa0+zi16iwaHj9Qop36wJjIuGizuVuZn3NLKva124zu9HM7jSzvGrjp/kzsByZ9kl9yBv3FJ0rN5P7zwtxVZVeRxKRBlDncnfOrXTOpTjnUoARwD5gpm/zI99vc8596I+gUneDxpzBgj63MLT4GxbqF6wiYcFf0zJjgTXOufV+ej3xs6MvuJ2FLcczesNUFs9+yes4IhJg/ir384FXq92/1syyzWyamSXU9AQzm2Jm6WaWXlBQ4KcYcigWEcGQq6axKqoPfb+5hdylWrNWJJTVu9zNLAY4C3jTN/QU0AtIAbYAD9X0POfcM865VOdcamJiYn1jSC3ExsXT+rI3KLVYmrz1S7YXbPE6kogEiD8+uU8AFjvntgE457Y55yqdc1XAVGCUH/YhftK2cw92nvkc7asK2DJ1MuVlpV5HEpEA8Ee5X0C1KRkz61ht29lAjh/2IX7Ue8Q4ckbew6DybL596nJdIlgkBNWr3M0sHjgZeKfa8P1mtsTMsoETgZvqsw8JjOFnXMm8zlcwctcHLJrxZ6/jiIif1avcnXPFzrk2zrmiamMXOecGO+eGOOfOcs5pYjdIjb78ARbFn0Dq6sf4ds7LXscRET/SGaphLCIykgG/eZnV0ckkf3Uzq7O+8jqSiPiJyj3MxTdrTpsr3ma3NafluxexZdNaryOJiB+o3IW2HZMonTyDePaxd9rPKNqla9CINHYqdwGg+4DRrD/xSXpUrmP9Uz+nvKzM60giUg8qd/mvAcefS1bKnQwpyyD7HxfpEEmRRkzlLj+QevYNzEuaQmrRbBY9/1uv44hIHanc5UfSLr2P+a3OYNTGaWS89aDXcUSkDlTu8iMWEcHwq58nK3Y0KUvuZvHHOgZepLFRuUuNYmJi6HPtm6yJTmbA1zeyZN7HXkcSkSOgcpdDatqsJe2unMX2yLYkfXQpq7Pnex1JRGpJ5S4/qVViJ6IvnUWZxZLwznlsWJ3tdSQRqQWVuxxWu6S+lF3wNpFUEf3KOeTrLFYRv3j6izXMWbYtIK+tcpda6dp3GIWTXqWZ20vptDMpKtzsdSSRRu3ZL9dy779XMHvp1oC8vspdai05ZQzrTn2edpXbKHjqDHbv2u51JJFG6eX567n7g+WcNrgD954zOCD7ULnLERl89ASWH/cPulWsY+Pfz2Tv3t1eRxJpVN5M38gf3s1hXP92PHreMKIiA1PDKnc5YsPGTmZp2oP0L1/G2sfPYl/xHq8jiTQK7327mVvfzubY5LY8+YvhxEQFroL9sUD2Ot/KS1lmlu4ba21mc8xste97Qv2jSjBJmXA536bew6CyLHIfP4vSfXu9jiQS1D7K2cpNr2eR2r01z1yUSmx0ZED3569/Nk50zqU451J9928DPnXOJQOf+u5LiBl25tUsHnYXg0ozWfX4REpLir2OJBKUPl+Zz3WvLmZIl5ZMu3QkcTGBLXYI3LTMROAF3+0XgEkB2o94LHXSdWQM/TNDStNZ/fhEykv3eR1JJKhkfvkvbnvpc/p2aM70y0bRrElUg+zXH+XugI/NLMPMpvjG2ldbO3Ur0P7gJ5nZFDNLN7P0goICP8QQr4w85wbmD/wTg0sWseKxSZSp4EUAyPzkNQZ8cin3Np3BS5ePpmVcdIPt2x/lPsY5NxyYAFxjZsdV3+iccxz4B4CDxp9xzqU651ITExP9EEO8lHbuzcwb8EeGlCxg+aOTNEUjYW/xRy8y8Mur2RjdjRFXPkNCfEyD7r/e5e6cy/N9zwdmAqOAbWbWEcD3Pb+++5Hgd9TkW1g08I+klC5g1aNnULynyOtIIp5I/+BZhsy7ge9ikml37WxatPnR5EXA1avczSzezJp/fxs4BcgB3gMu8T3sEmBWffYjjcfIc29hccrdDCzNZP3jE3Sik4SdRe8+ybCFt7CqyQA6X/cRLVq19SRHfT+5twe+MrNvgYXAB865j4B7gZPNbDUwzndfwsTwSdeRnfYwyeUr2PbkKRRtD8zp1SLBZtHbDzEi8w8si02h+/Uf0qyFd0eB1+vXts65tcDQGsa3A2Pr89rSuA2bcDnfxjSl39xr2fz3U6i48l+0aZ/kdSyRgJk/46+krbqfrLhR9LvhXWLj4j3NozNUJWCGjj2fVeOm0a5yKyVPn8KW9au8jiTid66qim+m3kjaqvvJjB9D/xtneV7soHKXABt87FmsP30GLVwR9vwE1i5L9zqSiN9U7C9n4RMXc3Te8yxsfSZDbnqXJrFNvY4FqNylAQwYNY4dP59JFBW0feMsls77yOtIIvVWWlLMkkcmMXrn+8zvchkjr32RyKiGO479cFTu0iC6D0qj4rLZFFlLen/0SzJnv3D4J4kEqaKdhax95FSG7fua+X1vJe1Xj2IRwVWnwZVGQlqHbv1ofvVnrIvuxdBvbmDBG/d5HUnkiBVuXs/2J8bSu2wZGakPknbB772OVCOVuzSoVokdSbrpE7Lj0xi97B4WTL0BV1XldSyRWlmzZAEVz5xEh8otrDzpWUac8WuvIx2Syl0aXFx8cwbd9B4LEs5kdN500h+drMsVSND79rPX6fDWWURQxeZz3mHw8ed4HeknqdzFE1HRMYy67kUW9LiakbvnsO6hsRRs3eB1LJEfcc4x/9W/MeiLK9kS1Rl+/Rm9h47xOtZhqdzFMxYRwehL/kZm2mN027+GyqdPZM2SBV7HEvmviv3lzP/7r0hbeS9L4o+i002f0a5zD69j1YrKXTw3bPyl5J39DhFU0fGtM8maM8PrSCLsKdpBzkOnc1ThWyzocAFDbn6fps1aeR2r1lTuEhR6pxxLxJTPyYtOYshXVzP/xT/qF63imfUrs9j52BgGlaSzaOAfGX3V00RENcwiG/6icpeg0bZTd7re/DmZzY8nbe3jZDxyri4bLA1u8eyXaTNjPM2q9rDi1JcZee4tXkeqE5W7BJXYps0ZfvNM5ne/muG7PyX/kTGsX/Wt17EkDFRWVPDN1BsZPu8aNkd3pfyK/zDo6NO9jlVnKncJOhYRQdqlf2PZ2Om0qtpJm1dOJeMjndEqgbNrez45D57K0XnPsyjhDLrd8h86dO3ldax6UblL0Bp03CTKL/+czdFdGTH/euY9fTX795d7HUtCTG7WXIqfHEP/kkwWDfojI69/mSax3l/Vsb5U7hLU2icl0/2WuSxsezZHbX2FVfefROHmdV7HkhDgqqpYOOMvJM2cRJSrYP2ZbzDy57eAmdfR/ELlLkEvJjaOUddOJ2P4vfQoX0XkM2N0uKTUy478zWQ/MJ5Rqx5iadPRxFw7j+TUcV7H8qs6l7uZdTWzz81smZktNbMbfON3mlmemWX5vk7zX1wJZyPO+g0Fv/iY7ZHtSPn6N8x/4lJKivd6HUsamSVff0DFP46h/74MFvS7lZTffUBC24ZfwDrQzDlXtyeadQQ6OucW+xbJzgAmAZOBvc65B2v7WqmpqS49XYs4SO2Ule4jc/otpG19hXURXak4+1l6D07zOpYEuf3lZSx84TbSNj3P5oiO7D/7OXoOOdrrWPViZhnOudSattX5k7tzbotzbrHv9h5gOdC5rq8nUltNYpuSdtU/yDnpBZpV7aHrW2fwzSt3U1VZ6XU0CVIbli9kw31HcUzeNL5tfSptfzu/0Rf74fhlzt3MugPDgO8vDHKtmWWb2TQzq3H5bzObYmbpZpZeUFDgjxgSZgYdN4moa+axIn4ER69+gBX3nUDemhyvY0kQqdhfzvwX/kCH1yaQUFlARtoTDL/hdeKatfQ6WsDVeVrmvy9g1gz4Avirc+4dM2sPFAIOuIsDUzeX/9RraFpG6sNVVbHo3Sfo9+29RFNBVp/rGXXe7UQ2stPFxb++W5FF+VtT6FuxksXxx9H14qdIbN/F61h+FZBpGd8LRwNvA684594BcM5tc85VOueqgKnAqPrsQ+RwLCKCUefcQMmUr1nVdBhHrX6QVfeNYd2KLK+jiQcqysv4+sX/R8dXx9GhIo/FIx9k2G9nhVyxH059jpYx4DlguXPu4WrjHas97GxA/0+WBtG+c0+G/O4jMkbcS6f9G+j46ji+nn4H5WWlXkeTBrJm8edsvG80x6x9jFXNRlJ51TcMP/3XQbe+aUOoz9EyY4AvgSXA95fv+z1wAZDCgWmZdcCVzrktP/VampYRf9u+bQMbX7qGlL1z2RDRhb1j/8aAY87yOpYESNGOAla+cguphbMosAQ2pf2ZEadeFDInJB3KT03L1HvO3R9U7hIomZ+8Rruv/x+d3TYymp9I5/MepkOXnl7HEj+pqqxi4fvPkJx1D63cbha2n8yAX9xLy1atvY7WIFTuEtZK9hWT+eqdDN/wPJVEktXzSkac93tiY2O9jib1sCZrLmX/uo0BFUtZHdWHyLMeC/nDGw+mchcBNn+3nMI3b2TIvvlstE7sSLuNISdfFJbzsY1Z/sZc8t6+jWG75rCdlqwfciPDJl6PRYbf0VEqd5Fqlv7nTZrO/Qs9qjawKrofleP+TP/R472OJYdRtHMHOW/eyYi8GRiwsOMvGHL+nWEzBVMTlbvIQfbv38+iWX+nV85jtGcHmXFpNDv9bpIHjfQ6mhykdN9eFr/zMH1zn6UNRWS0GEfHc+6hU/e+XkfznMpd5BBKiveQ/fa9DFj7HE1dKRktxtLxjDvo2ne419HCXlnJXnLee5xuy/9JW3axvEkKsaf9hR5Dj/c6WtBQuYscxu7tW1n51l0M3PwmsZTzbYvjaHHy7fQacpTX0cJOcfFesmc9Su9Vz5LITpZEDyHixNsZeLQuMHswlbtILRVu28yKd+9j6ObXaW4lLI5LI+r43zF49FgsxI+Z9lrR9q3kvPcYfdbPIJFdLI8ZzP7jbmXwMWfovT8ElbvIESraUciK9x6g77qXacVeVkT1o2T4rxly8sVERsd4HS+k5OVms2X2IwzM/xdxVk5OXCoxx99Mn9GnhfxJSPWlchepo9K9u8j54O90WPESXdwWCkkgN2kyvSZcS2LHJK/jNVoVFRVkzZ1FxKJnGF4ynzIXRWarU2h36k30HKDLUdWWyl2kniorK8n87E2i06cytCydMhdFVvPjaZJ6IYPHTNQVKGupMC+XNbP/SdKGmXSkgB20YFXX8+h12vX6x7IOVO4ifrRp9bds/eRJ+mz7gBYUk08b1nc+nXbHXkq3fiO8jhd09hbvZdnnrxGXM4OBJYuJMEdOk1r6giMAAAarSURBVGFUpFzEoJMuIKpJU68jNloqd5EAKC/dx5LPXseyX2VIySKirIpVkb3J73YG3Y45j669Bngd0TNlJXtYNncmFTkzGbD7G+KtlK20JbfzRLqNnULXnv28jhgSVO4iAVawZQPfff4CbdfOpGfFGgC+i+xOYeexJAw/m56DjyEiMrQvc7C9MJ/cee8Rvep9+u2ZR1PK2EkLctucQMsRP6f3qNOJ0PSVX6ncRRpQ/voVrP3qDZqvm02/8qVEmmMrbVjXcjQRPY+j58jxtO3Uw+uY9Va5v5zczM/ZmT2bhK1f0Xv/KiLNsZ2WrEw4gaYp5zDw6NOI1tFFAaNyF/FIwbY81s97h+jc2fTYu5gWFAOwKaITm1uNoCppDO0HHE1Sr4FERkZ6nPanFe0o4Lvsuexds4D4gkz6lGQTb6VUOiM3pi9FHY+lzZBT6ZFyoj6hNxCVu0gQqKqoYO3SBWzL/oSmeV/TuzSb5pQAsMfFsbFJMsWtB2CdhtKs2zA69ehPixatGjxnZUUFm9evpHBdDvs2Lye2MIf2e5bSpWrzgT+HMzZFdmFr61TodSJ9R59Gy9aJDZ5TPCp3MxsPPAZEAs865+491GNV7hKOqioqyFuVTv7KBVTkZdFy1zK67V9LnJX/9zGFtKIwuhPF8Unsb9EN17ILTVok0rRVO5q3bk/Lth2Jb55Qq8sWu6oqyspKKd1bxK7CPPYU5FGyczOVRVtwe7cRu28zrfeto2PlFprY/v8+r4AENsQNoLx9Cs17p9Fj8BjiW4bvlRiDSYOXu5lFAquAk4FNwCLgAufcspoer3IXOaC8vJyta7PZuS6b0m25ROz8jvjiDbTZv5n27KjxOftdJGXEsN+iqCSSCqKotEgqiSLK7acJZcS6MmIpI9Jq/vte4mIojGhLYWw3ylr2JLJdX1p2HUDHXoNp3rpDIP/IUg8/Ve6BmhgbBeQ659b6ArwGTARqLHcROSAmJoakfqkk9fvx39fykmKKCjaxe8c2indupbSogIo9hUSUbieiohQq90PVfqyqAqvaD64SImKojIqDqFhcdBxExUFMPNEtOxCX0IkWiV1o3aELTZsl0NWMrh78mSUwAlXunYGN1e5vAkZXf4CZTQGmACQl6cw0kcOJiYsnMakviUm6jrkcnmcH3jrnnnHOpTrnUhMT9csYERF/ClS558EP/ofXxTcmIiINIFDlvghINrMeZhYDnA+8F6B9iYjIQQIy5+6cqzCza4HZHDgUcppzbmkg9iUiIj8WsNPInHMfAh8G6vVFROTQQvtKRiIiYUrlLiISglTuIiIhKCguHGZmBcB6D3bdFij0YL/BTu/Loem9OTS9N4cWqPemm3OuxhOFgqLcvWJm6Ye6LkM40/tyaHpvDk3vzaF58d5oWkZEJASp3EVEQlC4l/szXgcIUnpfDk3vzaHpvTm0Bn9vwnrOXUQkVIX7J3cRkZCkchcRCUFhV+5mdq6ZLTWzKjNLPWjb7WaWa2YrzexUrzIGAzO708zyzCzL93Wa15m8ZmbjfT8buWZ2m9d5gomZrTOzJb6flbBeM9PMpplZvpnlVBtrbWZzzGy173tCoHOEXbkDOcA5wNzqg2Y2gAOXJh4IjAf+4VsLNpw94pxL8X2F9UXgfD8LfwcmAAOAC3w/M/I/J/p+VsL9WPfpHOiQ6m4DPnXOJQOf+u4HVNiVu3NuuXNuZQ2bJgKvOefKnHPfAbkcWAtWBKqtC+ycKwe+XxdY5Aecc3PhR6uZTwRe8N1+AZgU6BxhV+4/oaZ1Xzt7lCVYXGtm2b7/Zgb8v5FBTj8fP80BH5tZhm99ZPmh9s65Lb7bW4H2gd5hwK7n7iUz+wToUMOmO5xzsxo6T7D6qfcJeAq4iwN/ae8CHgIub7h00siMcc7lmVk7YI6ZrfB9gpWDOOecmQX8GPSQLHfn3Lg6PC3s1n2t7ftkZlOBfwU4TrALu5+PI+Gcy/N9zzezmRyYxlK5/882M+vonNtiZh2B/EDvUNMy//MecL6ZNTGzHkAysNDjTJ7x/QB+72wO/CI6nGld4EMws3gza/79beAU9PNysPeAS3y3LwECPoMQkp/cf4qZnQ08ASQCH5hZlnPuVOfcUjN7A1gGVADXOOcqvczqsfvNLIUD0zLrgCu9jeMtrQv8k9oDM80MDnTKDOfcR95G8o6ZvQqcALQ1s03An4B7gTfM7AoOXN58csBz6PIDIiKhR9MyIiIhSOUuIhKCVO4iIiFI5S4iEoJU7iIiIUjlLiISglTuIiIh6P8DOFhL2HgzXRwAAAAASUVORK5CYII=\n"
          },
          "metadata": {
            "needs_background": "light"
          }
        }
      ],
      "source": [
        "#defining new set of given node information: x, y and y'\n",
        "x       = np.array([1.0, 3.0, 5.0])\n",
        "y       = np.array([5.0, 1.0, 5.0])\n",
        "y_prime = np.array([-4.0, 0.0, 4.0])\n",
        "\n",
        "f7      = hermit( x, y, y_prime)\n",
        "data    = f7.linspace(n=50, domain=[-11, 10.5])\n",
        "test_x  = np.linspace(-10, 10, 50, endpoint=True)\n",
        "test_y  = (test_x-3)**2 + 1\n",
        "\n",
        "plt.plot(data[0], data[1])\n",
        "plt.plot(test_x, test_y)\n",
        "plt.show()"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "UgP_TWwOPlvj",
        "nbgrader": {
          "cell_type": "markdown",
          "checksum": "c07a0720672d1b686830ce5cef78046c",
          "grade": false,
          "grade_id": "cell-580b64f506c61583",
          "locked": true,
          "schema_version": 3,
          "solution": false,
          "task": false
        }
      },
      "source": [
        "## Part 2: Polynomial Interpolation Using Newton's Divided Difference Form\n",
        "---\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "AgAbn67kPlvj",
        "nbgrader": {
          "cell_type": "markdown",
          "checksum": "030540dd5ee10acfca978c35d1e76cb8",
          "grade": false,
          "grade_id": "cell-95fa6d3532e84eff",
          "locked": true,
          "schema_version": 3,
          "solution": false,
          "task": false
        }
      },
      "source": [
        "### Newton's Divided Difference Form\n",
        "\n",
        "Newton form of a $n$ degree polynomial:\n",
        "\n",
        "$$p_n(x) = \\sum_{k=0}^{n} a_kn_k(x),$$\n",
        "where the basis is:\n",
        "$$n_k(x) = \\prod_{j=0}^{k-1}(x-x_j),$$\n",
        "$$ n_0(x)=1,$$\n",
        "\n",
        "and the coefficients are: $$a_k = f[x_0, x_1, ..., x_k],$$\n",
        "\n",
        "where the notation $f[x_0, x_1,\\dots,x_k]$ denotes the divided difference.\n",
        "\n",
        "By expanding the Newton form, we get:\n",
        "\n",
        "$$p(x) = f [x_0] + (x-x_0) f[x_0,x_1] + (x-x_0) (x-x_1) f[x_0,x_1,x_2] + \\dots + (x-x_0) (x-x_1) \\dots (x-x_{k-1}) f[x_0, x_1, \\dots, x_k]$$\n",
        "\n",
        "### Tasks: \n",
        "1. Complete the `calc_div_diff(x,y)` function which takes input `x` and `y`, and calculates all the divided differences. You may use the lambda function `difference()` inside the `calc_div_diff(x,y)` function to calculate the divided differences.\n",
        "\n",
        "2. Complete the `__call__()` function which takes an input `x`, and calculates `y` using all the difference coefficients. `x` can be a single value or a numpy. In this case, it is a numpy array.\n",
        "\n",
        "`res` variable must contain all results (corresponding y for x)."
      ]
    },
    {
      "cell_type": "code",
      "source": [
        "class Newtons_Divided_Differences:\n",
        "  \n",
        "    def __init__(self, differences):\n",
        "        self.differences = differences\n",
        "\n",
        "    def __call__(self, x):\n",
        "        '''\n",
        "        this function is for calculating y from given x using all the difference coefficients\n",
        "        x can be a single value or a numpy\n",
        "        the formula being used:\n",
        "        f(x) = f [x0] + (x-x0) f[x0,x1] + (x-x0) (x-x1) f[x0,x1,x2] + . . . + (x-x0) (x-x1) . . . (x-xk-1) f[x0, x1, . . ., xk]\n",
        "\n",
        "        work on this after implementing 'calc_div_diff'. Then you should have\n",
        "        f[x0], f[x0,x1]. . . . . ., f[x0, x1, . . ., xk] stored in self.differences\n",
        "\n",
        "        'res' variable must return all the results (corresponding y for x)\n",
        "        '''\n",
        "\n",
        "        res = np.zeros(len(x))                            #Initialization to avoid runtime error.\n",
        "        \n",
        "        for p in  range(len(x)):\n",
        "          for q in  range(len(data_x)):                   #0,1,2,3\n",
        "\n",
        "            coefff=1\n",
        "            for r in range(q):                             #(x-x0) (x-x1) being created\n",
        "              coefff  = coefff*(x[p]-data_x[r]) \n",
        "          \n",
        "            res[p] +=  coefff*self.differences[q] #for each of the value of the elements of x passes in the list is, x\n",
        "            \n",
        "        return res"
      ],
      "metadata": {
        "id": "Ids8f_Y3Mzt9"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "deletable": false,
        "id": "iLUOVj3APlvj",
        "nbgrader": {
          "cell_type": "code",
          "checksum": "50cdde63e49e009a199201f86744564a",
          "grade": false,
          "grade_id": "cell-9e446c528b62c87f",
          "locked": false,
          "schema_version": 3,
          "solution": true,
          "task": false
        }
      },
      "outputs": [],
      "source": []
    },
    {
      "cell_type": "code",
      "source": [],
      "metadata": {
        "id": "EPdRXfPdskLp"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "deletable": false,
        "id": "HZhqt-6JPlvj",
        "nbgrader": {
          "cell_type": "code",
          "checksum": "abd6bbf6334a986a608c38206f702989",
          "grade": false,
          "grade_id": "cell-fbe4c16afb6b3b32",
          "locked": false,
          "schema_version": 3,
          "solution": true,
          "task": false
        }
      },
      "outputs": [],
      "source": [
        "# basic rule for calculating the difference, implanted in the lambda function. \n",
        "# You may use it if you wish\n",
        "difference = lambda y2, y1, x2, x1: (y2-y1)/(x2-x1) #a er formula\n",
        "\n",
        "def calc_div_diff(x,y):                 #ekhane 'a' er list generate korbo\n",
        "    assert(len(x)==len(y))\n",
        "    #write this function to calculate all the divided differences in the list 'b'\n",
        "    b = [y[0]]  #initializing\n",
        "    #----------------------------------------------\n",
        "    # YOUR CODE HERE\n",
        "    y_var = y\n",
        "    for m in range(1, len(x)):\n",
        "      store = []\n",
        "      for n in range(len(x) - m):\n",
        "        store.append(difference(y_var[n+1],y_var[n], x[n+m], x[n] ))\n",
        "      y_var = store.copy()\n",
        "      b.append(y_var[0])\n",
        "    return b\n",
        "\n",
        "\n",
        "\n",
        "\n"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "FancUC_xPlvk",
        "nbgrader": {
          "cell_type": "code",
          "checksum": "a23dd58fcab86885a166b1e138b9ea08",
          "grade": true,
          "grade_id": "cell-72b9329c6c32517e",
          "locked": true,
          "points": 5,
          "schema_version": 3,
          "solution": false,
          "task": false
        },
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "outputId": "7c2552ad-6a85-45ef-cd65-f4bcb5c38c5b"
      },
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "[-60.0, -20.0, 53.0, -32.833333333333336, 14.041666666666668, -3.465277777777778, 0.7232142857142857]\n"
          ]
        }
      ],
      "source": [
        "data_x = [-3.,-2.,-1.,0.,1.,3.,4.]\n",
        "data_y = [-60.,-80.,6.,1.,45.,30.,16.]\n",
        "\n",
        "test = calc_div_diff(data_x, data_y)\n",
        "print(test)\n",
        "\n",
        "assert len(test) == len(data_x)\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "affgKuuNPlvk",
        "nbgrader": {
          "cell_type": "markdown",
          "checksum": "e2b6e9e71eaf31b857393044c2768ac0",
          "grade": false,
          "grade_id": "cell-f6a9f264b02b2642",
          "locked": true,
          "schema_version": 3,
          "solution": false,
          "task": false
        }
      },
      "source": [
        "### Plotting the polynomial\n",
        "* `data_x` and `data_y` are the coordinates of the given nodes.\n",
        "\n",
        "* `differences` is a list which contains the divided differences as each of its elements: $f[x_0], f[x_0,x_1], f[x_0,x_1,x_2], \\dots$\n",
        "\n",
        "* `obj` is an object of type `Newtons_Divided_Differences`. Creating the object runs the constructor of the class where the `difference` are stored in `self.differences`.\n",
        "\n",
        "* `X` contains $x_i$ values through which we want to plot our polynomial.\n",
        "\n",
        "* Calling the object using `obj(X)` executes the `__call__()` function of the class, which returns a numpy array containing the corresponding $y_i$ values, and storing them in variable `F`.\n",
        "\n",
        "* Using `plt.plot(X,F)`, we plot the $(x_i, y_i)$ pairs of the polynomial."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "deletable": false,
        "editable": false,
        "id": "akcjxODFPlvk",
        "nbgrader": {
          "cell_type": "code",
          "checksum": "3f9ba22ae7061e1d74101e40a2d5a7a4",
          "grade": true,
          "grade_id": "cell-a16b9abefbcdc519",
          "locked": true,
          "points": 5,
          "schema_version": 3,
          "solution": false,
          "task": false
        },
        "colab": {
          "base_uri": "https://localhost:8080/",
          "height": 288
        },
        "outputId": "56b57c2c-b887-48ef-f09a-6c38f85f2e2f"
      },
      "outputs": [
        {
          "output_type": "display_data",
          "data": {
            "text/plain": [
              "<Figure size 432x288 with 1 Axes>"
            ],
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYMAAAD7CAYAAACIYvgKAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3deXhc9X3v8fd3RstoG1mbJVmyLW+AFcJiHGyWpARSsH2TQEouDzQNbkLjLrht7r19WpK0IU0uaZab5An3NvTS1g30kpA0KYWkbI4bnhYoBpuAjW3A8gJeJFm2jBZrnZnf/WOO7EEeWduMzmjm83qeeebMb86c+c7RaD5zzu93zphzDhERyW0BvwsQERH/KQxERERhICIiCgMREUFhICIiKAxERIQJhIGZzTezX5rZbjPbZWZ/7LVXmtlmM9vrXVd47WZm95pZi5ntMLMVCcta782/18zWJ7RfZmY7vcfca2aWjhcrIiLJTWTLIAL8D+dcM7AauNPMmoG7gC3OuWXAFu82wFpgmXfZANwH8fAA7gZWAZcDd48EiDfPZxIet2b6L01ERCYqb7wZnHOtQKs33WNme4AG4EbgGm+2B4BngD/z2h908aPZXjCzOWZW78272TnXCWBmm4E1ZvYMEHbOveC1PwjcBDxxrrqqq6tdU1PTJF6qiIhs3779uHOuZnT7uGGQyMyagEuBrUCtFxQAbUCtN90AHEp42GGv7Vzth5O0n1NTUxPbtm2bTPkiIjnPzN5K1j7hDmQzKwV+CnzWOdedeJ+3FZD281qY2QYz22Zm2zo6OtL9dCIiOWNCYWBm+cSD4CHn3D97ze3e7h+862Ne+xFgfsLDG722c7U3Jmk/i3PufufcSufcypqas7ZyRERkiiYymsiAvwf2OOe+nXDXY8DIiKD1wKMJ7bd7o4pWA13e7qSngOvNrMLrOL4eeMq7r9vMVnvPdXvCskREZAZMpM/gKuCTwE4ze8Vr+zzwNeDHZnYH8BZwi3ff48A6oAXoAz4F4JzrNLOvAC958315pDMZ+APg+0AR8Y7jc3Yei4hIatlsPYX1ypUrnTqQRUQmx8y2O+dWjm7XEcgiIqIwEBGRHAyDB//zID979ajfZYiIZJScC4OHXzzEv/wq6chVEZGclXNhUFceoq17wO8yREQySm6GQZfCQEQkUe6FQTjEiVNDDEaifpciIpIxci8MykMAHOse9LkSEZHMkXNhUO+FgfoNRETOyLkwqAvHw6BV/QYiIqflXhiMbBl09ftciYhI5si5MCgL5VNSEKStS30GIiIjci4MYORYA20ZiIiMyN0wUJ+BiMhpuRkG4SKFgYhIgpwMg/ryEO09g0Rjs/O3HEREUi0nw6C2PEQ05jjRq05kERHI0TCo17EGIiLvkpNhUKejkEVE3iW3w0BbBiIiQI6GQWVxAflB05aBiIgnJ8MgEDBqwzrWQERkRE6GAcSHl7bq/EQiIkAOh0FtOES7ftNARATI4TAY2TJwTgeeiYjkbBjUhkMMDMfo7o/4XYqIiO9yNgzqy4sAaNXZS0VEcjcM6soLAR2FLCICOR0G8S2DdoWBiEjuhsHcskLMtGUgIgI5HAb5wQDVpYW06yhkEZHcDQMYGV6qMBARyekwiB94pjAQEcnpMNCWgYhIXE6HQW04RFf/MP1DUb9LERHxVU6HQb1+5EZEBMjxMBj5kRudvVREcl1uh4H3W8jqRBaRXJfbYXB6y0BhICK5bdwwMLNNZnbMzF5LaPuSmR0xs1e8y7qE+z5nZi1m9oaZ3ZDQvsZrazGzuxLaF5nZVq/9R2ZWkMoXeC7FBXmEQ3n6xTMRyXkT2TL4PrAmSft3nHOXeJfHAcysGbgVeI/3mO+ZWdDMgsBfA2uBZuA2b16Ar3vLWgqcBO6YzguarPryIoWBiOS8ccPAOffvQOcEl3cj8LBzbtA5dwBoAS73Li3Ouf3OuSHgYeBGMzPgWuAn3uMfAG6a5GuYltrykEYTiUjOm06fwUYz2+HtRqrw2hqAQwnzHPbaxmqvAt5xzkVGtc+Y+nBIWwYikvOmGgb3AUuAS4BW4Fspq+gczGyDmW0zs20dHR0pWWZdeYiO3kGGo7GULE9EZDaaUhg459qdc1HnXAz4W+K7gQCOAPMTZm302sZqPwHMMbO8Ue1jPe/9zrmVzrmVNTU1Uyn9LHXlIZyDjp7BlCxPRGQ2mlIYmFl9ws2PASMjjR4DbjWzQjNbBCwDXgReApZ5I4cKiHcyP+biv0b/S+Dj3uPXA49Opaap0vBSERHIG28GM/shcA1QbWaHgbuBa8zsEsABB4HfBXDO7TKzHwO7gQhwp3Mu6i1nI/AUEAQ2Oed2eU/xZ8DDZvY/gV8Bf5+yVzcBIweeqd9ARHLZuGHgnLstSfOYH9jOuXuAe5K0Pw48nqR9P2d2M804nZ9IRCTHj0AGKC/KpzAvQJvOTyQiOSznw8DMqC8P0datDmQRyV05HwYQ70TWloGI5DKFAfFOZI0mEpFcpjAA6sqLONY9SCzm/C5FRMQXCgOgLlzIUDRGZ9+Q36WIiPhCYUB8ywB0rIGI5C6FAWeOQlYYiEiuUhigA89ERBQGQHVpIcGAactARHKWwgAIBoy5ZYUaXioiOUth4KkrD9Gu3UQikqMUBp66cIijOgpZRHKUwsDTWFHEkZP9OvBMRHKSwsCzsKqEwUiM9h7tKhKR3KMw8DRVlQBw8Hifz5WIiMw8hYFnYVUxAG+dOOVzJSIiM09h4Jk3p4j8oHHwhLYMRCT3KAw8wYAxv7JYWwYikpMUBgmaqkq0ZSAiOUlhkGBhVXzLwDkNLxWR3KIwSNBUVULfUJSOXv0esojkFoVBgjMjirSrSERyi8IgwZljDdSJLCK5RWGQoKGiiGDAtGUgIjlHYZAgPxhgfkURBzW8VERyjMJglIVVJdoyEJGcozAYpamqmIMaXioiOUZhMMrCqhJ6BiKc7Bv2uxQRkRmjMBilqTo+vFT9BiKSSxQGoyz0hpfqHEUikksUBqM0VhQRMP2ugYjkFoXBKIV5QebNKdKWgYjkFIVBEjp7qYjkGoVBEiNnLxURyRUKgySaqko42TdMl4aXikiOUBgkcfrspZ3aOhCR3KAwSKKpOj689IDOXioiOUJhkMSCSv2ugYjklnHDwMw2mdkxM3stoa3SzDab2V7vusJrNzO718xazGyHma1IeMx6b/69ZrY+of0yM9vpPeZeM7NUv8jJCuUHqS8P6ShkEckZE9ky+D6wZlTbXcAW59wyYIt3G2AtsMy7bADug3h4AHcDq4DLgbtHAsSb5zMJjxv9XL6IjyjSloGI5IZxw8A59+9A56jmG4EHvOkHgJsS2h90cS8Ac8ysHrgB2Oyc63TOnQQ2A2u8+8LOuRdc/DShDyYsy1dNVSUaXioiOWOqfQa1zrlWb7oNqPWmG4BDCfMd9trO1X44SbvvFlaVcLx3iJ4BDS8Vkew37Q5k7xv9jJz838w2mNk2M9vW0dGR1udqqlInsojkjqmGQbu3iwfv+pjXfgSYnzBfo9d2rvbGJO1JOefud86tdM6trKmpmWLpE3Pm7KUKAxHJflMNg8eAkRFB64FHE9pv90YVrQa6vN1JTwHXm1mF13F8PfCUd1+3ma32RhHdnrAsX40ceKYRRSKSC/LGm8HMfghcA1Sb2WHio4K+BvzYzO4A3gJu8WZ/HFgHtAB9wKcAnHOdZvYV4CVvvi8750Y6pf+A+IilIuAJ7+K7ksI8asoK1YksIjlh3DBwzt02xl3XJZnXAXeOsZxNwKYk7duAC8erww/x30PWbiIRyX46AvkcFmp4qYjkCIXBOTRVFdPePUjfUMTvUkRE0kphcA4jI4re7tSuIhHJbgqDc2jywkC/hywi2U5hcA4LTh94pn4DEcluCoNzKC/Kp7KkQCOKRCTrKQzGod9DFpFcoDAYR/zspdoyEJHspjAYx8KqYo529TMwHPW7FBGRtFEYjGNJTSnOwb6OXr9LERFJG4XBOJrnhQHYfbTb50pERNJHYTCOpqoSiguC7G5VGIhI9lIYjCMYMC6oK2OXtgxEJIspDCageV6YPUe7iZ+UVUQk+ygMJqC5vpyewQiHT/b7XYqISFooDCZgpBNZu4pEJFspDCbggroyAoY6kUUkaykMJiCUH2RJTSm7j3b5XYqISFooDCaoeV5YxxqISNZSGExQc32Yo10DnDw15HcpIiIppzCYoJFO5D3qNxCRLKQwmKDmeu+0FAoDEclCCoMJqiotpC4c0vBSEclKCoNJUCeyiGQrhcEkNNeHaeno1W8biEjWURhMQvO8MNGYY2+7fttARLKLwmASznQi6+AzEckuCoNJWFBZTGlhnjqRRSTrKAwmIRAwlteXqRNZRLKOwmCSmuvD7GntJhbTbxuISPZQGExS87wwp4aivN3Z53cpIiIpozCYpOb6ckBHIotIdlEYTNKy2lLyAsYunc5aRLKIwmCSQvlBls4tVSeyiGQVhcEUNNeHtZtIRLKKwmAKmueFae8e5HjvoN+liIikhMJgCkaORNZvG4hItlAYTMHID93oSGQRyRYKgymYU1xAw5widSKLSNaYVhiY2UEz22lmr5jZNq+t0sw2m9le77rCazczu9fMWsxsh5mtSFjOem/+vWa2fnovaWYsVyeyiGSRVGwZfNA5d4lzbqV3+y5gi3NuGbDFuw2wFljmXTYA90E8PIC7gVXA5cDdIwGSyZrnhdnf0Uv/kH7bQERmv3TsJroReMCbfgC4KaH9QRf3AjDHzOqBG4DNzrlO59xJYDOwJg11pVRzfZiYg9fbtHUgIrPfdMPAAU+b2XYz2+C11TrnWr3pNqDWm24ADiU89rDXNlZ7Rrtk/hwAtr910udKRESmb7phcLVzbgXxXUB3mtkHEu90zjnigZESZrbBzLaZ2baOjo5ULXZK6spDLK4u4fl9J3ytQ0QkFaYVBs65I971MeAR4vv8273dP3jXx7zZjwDzEx7e6LWN1Z7s+e53zq10zq2sqamZTukpceXSKrbuP8FwNOZ3KSIi0zLlMDCzEjMrG5kGrgdeAx4DRkYErQce9aYfA273RhWtBrq83UlPAdebWYXXcXy915bxrlxSzamhKDsO66R1IjK75U3jsbXAI2Y2spwfOOeeNLOXgB+b2R3AW8At3vyPA+uAFqAP+BSAc67TzL4CvOTN92XnXOc06poxqxdXAfCf+45z2cKMHwAlIjKmKYeBc24/cHGS9hPAdUnaHXDnGMvaBGyaai1+qSwpoLk+zHMtJ9h47TK/yxERmTIdgTxNVy6pYvvbJxkY1vEGIjJ7KQym6aql1QxFYhpiKiKzmsJgmt63qJJgwHh+33G/SxERmTKFwTSVFuZxcWO5jjcQkVlNYZACVy2tZsfhLnoGhv0uRURkShQGKXDFkiqiMceLB2bFiFgRkbMoDFJgxYIKCvMCPNeiXUUiMjspDFIglB9kZVOFOpFFZNZSGKTIlUuqeb2thxO9g36XIiIyaQqDFLlyiXdqiv3aVSQis4/CIEXe21BOWWGehpiKyKykMEiRvGCAVYsr+U+FgYjMQgqDFLpiSTUHjp/i6Dv9fpciIjIpCoMUumppvN9Au4pEZLZRGKTQeXPLqCop4PkWDTEVkdlFYZBCgYCxekkVz+87QfznG0REZgeFQYpdtaSatu4BDhw/5XcpIpJlnHNp+6KpMEixkeMNntOuIhFJscd3tvGx7z3Pse6BlC9bYZBiC6uKWVxdws93tPpdiohkkVjM8d0tb9IzMExVaWHKl68wSDEz4+bLGtl6oJNDnX1+lyMiWeLJXW282d7LH123jGDAUr58hUEafOzSBszgpy8f9rsUEckCsZjj3i17WVxTwocvmpeW51AYpMG8OUVcuaSKf375iEYVici0Pb27jdfbevjDa5emZasAFAZpc/OKRt7u7OOlgyf9LkVEZrF4X0ELi6tL+EiatgpAYZA2ay6so6QgyE+3a1eRiEzd5j3t7GntZuO1S8kLpu8jW2GQJsUFeax9bz3/urOV/qGo3+WIyCzknOO7v9hLU1UxH704fVsFoDBIq5tXNNI7GOHp3W1+lyIis9Dm3e3sbu1m47XL0rpVAAqDtFq1qJKGOUX8RLuKRGSSnHN8d8teFlYVc9Ml6d0qAIVBWgUC8WMOnm05TmuXTmstIhO3Zc8xdh3t5s4PprevYITCIM1uXtGAc/DIr474XYqIzBIjWwXzK4v42KUNM/KcCoM0W1hVwvuaKvjp9sM65kBEJuSXbxxj55Eu/vCDy8ifga0CUBjMiJtXNLKv4xSvHu7yuxQRyXDdA8N8+We741sFK2ZmqwAUBjNi3UX1FOYFdMyBiJyTc44//acdHDrZz7dvuWTGtgpAYTAjwqF8bnhPHY+9epTBiI45EJHk/u4/DvDkrjY+t/YC3tdUOaPPrTCYITdf1khX/zD/tueY36WISAZ68UAnX3vyddZeWMcdVy+a8edXGMyQq5dWUxsu5P9tfcvvUkQkwxzrGWDjD15mQWUx3/j4RZil52R056IwmCHBgPGZ9y/muZYTPPOGtg5EJC4SjfFHP/wV3QPD3PdbKygL5ftSh8JgBt1+RRNNVcXc8697iERjfpcjIhngW5vf5IX9ndxz03u5oC7sWx0KgxlUkBfgrrXL2Xusl4dfOuR3OSLis82727nvmX3cdvkCbr6s0ddaFAYz7Ib31LJqUSXf2fwm3QPDfpcjMmNiMcdwNKaDLz0/3naIjT94mQsbwtz9kWa/yyHP7wJGmNka4LtAEPg759zXfC4pLcyMv/hwMx/5P8/yvV/u4661F/hdksi4BoajHO8d5HjvECd6B09Pnzw1RM9AhJ7BYbr7I/QMDNM9EOHUYIThaIzhaDwAIjFHNHYmBPKDRkEwQEHemUs4lE9lSQEVxQVnrksLqC0rZEFVMfMriikpzJiPrCnrH4ryxUdf45+2H+aqpVV899ZLCeUH/S4rM8LAzILAXwO/DhwGXjKzx5xzu/2tLD0ubCjnNy5tZNOzB/jEqgXMryz2u6TZ6aGH4AtfgLffhgUL4J574BOf8G85mWYCrysWcxw/NUhb1wCtXQO0dw/Q1jVAW3d8ur17kPauAXoGI0mfoig/SLgoj7JQPmWhPOYUFzC/spjSwjzygwHvYuQFjfxggKAZkZhjKBpjKHLmMhiJ0j0QofPUEG939tHZO5T0OatKCmisLGZBZTGLqoo5r66MC+rKaKoqeffJ3DLpb5pQy3BDI9/9tfX8ZP5q/ujapfzxh85L289YTpZlwiabmV0BfMk5d4N3+3MAzrm/GusxK1eudNu2bZuhClOvrWuAb/72F/n8s/9I1Yl2/9+wmWa8f+aHHoING6Cv70xbcTHcf//k1mGqlpNpHnoIt2EDlvC6hgtD/GzjX7JlxYdoT/jwj8Te/RmQFzDmlhVSWx6itixEbbiQueEQ1aUFVJUUUl1WSHVpAdWlhWn9RjsUiXGyb4jWrgEOdfbxdmcfh0/Grw919nP4ZB8jpRcEAyyZW8oFdWX8l52/5IPf+gLBgYQzBfv1N03y/urPL+TAV79D85/8/szW4jGz7c65lWe1Z0gYfBxY45z7He/2J4FVzrmNYz1mtocBDz3E8B2/Q/7gwJm2bPgQSoUk/0CxoiJe/Ytvsu2qtXT2DbHh9mup6Gg966Gt5XNZ+9//kUjC7gnnHHmBAMGAkReIf0sNBuLfWB/5xieoe6f9rOWcrK7nf2/6BSWFQUoK8ygpCFJckBefLgxSPHK7II/iwiBF+UFC+cGUf8tzzjEcdfQPR+kdjNDVN0xX/zDdA951/zAn+4Y40TvE8d4hjvcOcuLUIA9/9TYaus4ewnykfC6f/PMfUVceil/CIerLQ9SVF1EXjrdVlRQQyJBvq+cyMByl5Vgvb7T18EZ7D6+39fBGWzc//fpv0tjdcdb8gw2NRPcfoLhgBneINDXBW0mOLVq4EA4enLk6EmRFGJjZBmADwIIFCy57K9lKni0y8E3ip+FojH0dvexp7ebaGy6nPMkH/eFwDVf//j+QFzDe/KsPE+Ds964z40uP7CAvGDj9wW8YURffZz0cjRGNOSIxx3AkxjduuRRLspwYxkVffIJTQxEm8y+SHzQK84KE8gMU5gUpzIuH0FkXi9cUizmvNk5PD0ai9A/FGBiO0j8cfde+9mSCAaOypICqkgJqygqpKingO795GZascDOIZfewZhcIJH3tMYyld/2MpXNLuahxDhc1lvPehnKW14dTvoXTPxTlyV2t3LRiQdL3l59/h0wPg5zbTUQgQNJPmVz4Z3WO/cdPsXV/Jy+/fZLdR7tpOdbLkHfsxf6vf2TMD/ru3kHCRXnYokWpCdNxQtm5+LfyU4NRTg1G6BuK0jcU4dRQlD7v9qmhCP1DUQYj8Q/wgeEYA5EoA8NRhiIxYs4Ribr4tdeRGnOOgBkBi4dD/Dr+wR7KCxIqiG9tFOUHKSqIh0o4lE+4KI9wUT7hUD7lRfmEi/IpK8w7+5t8Ln/ZGOO199c38jcPbGHnkS52HH6H471DQHy32OKaEpbNLWPp3FKW1ZaybG4ZTdXFFOZNLCScc3T3R9jT1s0jLx/hX3e20jsY4YX/+2nq3klykGkGbhlkRAcy8BKwzMwWAUeAW4Hf9LekNFuwIPk/64IFM19Lmjnn2Husl637T/DCgU627u/keO8gEO8QbJ4X5v3LmmieF2Z5fRj74fx4X8EotmAB5cXe0Zn33JN8X/8990yuuHGWY2YUF+RRXJBHTVnh5Jbtp1Stn9lojNde9M2v8d9+/Twg/p5s7Rpgx+Eudh55hzfaetl1tIsnXms93Q8R9PpORkK3POESDNjpzvb27niH+8Bw/MtMcUGQde+t5+OXNTL3wm/D786Ov0NGhIFzLmJmG4GniA8t3eSc2+VzWemV5A07kF9I5It/SamPZaXKwHCU51qO84s9x9iyp51jPfEP//ryEFcvrWLV4ipWLapkUXXJ2edh+epXx/8gG+lXme6IkVQtJ9Nk6+uaiAm8djNj3pwi5s0pYs2FdafbB4ajHDh+ir3Hetnb3kNb1wBd/fH+mUOdfezypoejjrnhQurCId7TUM6HltdSVx6isaKI9y+rOTMEdvEnwM5dS6bIiN1EUzHrdxPBu0bMDNY38Pn33cab132UH3xmlW/nJ5mOYz0DbPE+/J9tOc7AcIySgiAfOK+Ga86v4YrF1cyvLJrYSbgyaWigSBbJ6D6DqciKMBjl315vZ8OD27lsYQUPfPryjDgQZTydp4Z44rVWfvbqUbYe6MQ5aJhTxIeWz+W65bWsWlw54f2uIpJ+md5nIMC1F9TyrVsu5rM/eoWNP3iZ+37rshn9paOJ6uof5qldbfx8RyvPtRwnGnMsrinhD69dxtoL67igrsyXU/CKyNQpDDLMjZc00D0Q4S/+5TX+9Cc7+NZ/vTgjxnwPDEfZsucYj75yhGfe6GAoGmN+ZRG/+4HFfPiieSyvVwCIzGYKgwz0ydUL6eob4n89/SZloTzu/sh7fDlkPRKN8fy+Ezz6ylGe2tVG72CEuWWFfPKKhXz04nlc1FiuABDJEgqDDHXnB5fS1T/M3/7HAV7Yf4LPrV3ONefXpP3DNxKNsfVAJ0+81sqTr7VzvHeQslAe695bx02XNLBqcVXGnEtFRFJHYZChzIzPr1vOpQsq+PqTr/Op77/ElUuq+Py65VzYUJ7S5xqKxHhu33Ge3NnG07vbONk3TFF+kA9eUMNHL57HNefPnRWd2SIydRpNNAsMRWI8tPUt7t2yl5N9w3zs0gb+5IbzaZhTNKXlRWOOPa3dbD3QyYsHTvD8vhP0DEQoLczjuuVzWXthHb923lyKChQAItlGQ0uzQPfAMPc9s49Nzx7AObiwIcz5dWWcX1vmnco3TGVJARA/wrJnMEJn7xAnTg3ReWqIN9t7ePFAJ9vfOkmvd3rgBZXFrF5cyZoL67hqabWGgYpkOYVBFjnyTj//8OwBdh7p4o32Ht7pO/OLadWlhQQDcPLU8Olz/SQ6r7aUyxdVcvmiKi5vqqSuPDSTpYuIz3ScQRZpmFPEn384/jN5zjk6ega90/f28GZ7D2ZQWVJIZUk+lSXxs1hWlhSwoLKYCm/LQUQkkcJgljMz5oZDzA2H+MB5NX6XIyKzVOYd3ioiIjNOYSAiIgoDERFRGIiICAoDERFBYSAiIigMREQEhYGIiDCLT0dhZh3AW0nuqgaOz3A50zXbap5t9YJqnimqOf2mW+9C59xZR6jO2jAYi5ltS3bejUw222qebfWCap4pqjn90lWvdhOJiIjCQEREsjMM7ve7gCmYbTXPtnpBNc8U1Zx+aak36/oMRERk8rJxy0BERCZp1oeBmf3IzF7xLgfN7JUx5jtoZju9+Xz9iTQz+5KZHUmoe90Y860xszfMrMXM7prpOhPq+KaZvW5mO8zsETObM8Z8vq/j8daZmRV675kWM9tqZk0zX+W76plvZr80s91mtsvM/jjJPNeYWVfC++WLftQ6qqZz/q0t7l5vPe8wsxV+1OnVcn7CunvFzLrN7LOj5vF9HZvZJjM7ZmavJbRVmtlmM9vrXVeM8dj13jx7zWz9lApwzmXNBfgW8MUx7jsIVPtdo1fLl4A/GWeeILAPWAwUAK8CzT7Vez2Q501/Hfh6Jq7jiawz4A+Av/GmbwV+5PN7oR5Y4U2XAW8mqfka4Od+1jnZvzWwDngCMGA1sNXvmhPeI23Ex9pn1DoGPgCsAF5LaPsGcJc3fVey/z2gEtjvXVd40xWTff5Zv2UwwswMuAX4od+1pMjlQItzbr9zbgh4GLjRj0Kcc0875yLezReARj/qmICJrLMbgQe86Z8A13nvHV8451qdcy970z3AHqDBr3pS6EbgQRf3AjDHzOr9Lgq4DtjnnEt2wKqvnHP/DnSOak58vz4A3JTkoTcAm51znc65k8BmYM1knz9rwgB4P9DunNs7xv0OeNrMtpvZhhmsaywbvc3nTWNs+jUAhxJuHyYzPiQ+TfwbXzJ+r+OJrLPT83gB1wVUzUh14/B2WV0KbE1y9xVm9qqZPWFm75nRwpIb72+dqe/fWxn7C2OmrSaWsA0AAAKrSURBVGOAWudcqzfdBtQmmScl63pW/Aaymf0CqEty1xecc49607dx7q2Cq51zR8xsLrDZzF73kjgtzlUzcB/wFeL/UF8hvnvr0+mqZSImso7N7AtABHhojMXM6DrOJmZWCvwU+KxzrnvU3S8T363R6/Uv/QuwbKZrHGXW/a3NrAD4KPC5JHdn4jp+F+ecM7O0Df+cFWHgnPvQue43szzgN4DLzrGMI971MTN7hPguhbS9ecereYSZ/S3w8yR3HQHmJ9xu9NrSYgLr+LeBDwPXOW9HZZJlzOg6TmIi62xknsPe+6YcODEz5SVnZvnEg+Ah59w/j74/MRycc4+b2ffMrNo559v5dCbwt57R9+8ErQVeds61j74jE9exp93M6p1zrd5utmNJ5jlCvM9jRCPwzGSfKFt2E30IeN05dzjZnWZWYmZlI9PEO0RfSzbvTBi17/RjY9TyErDMzBZ532huBR6bifpGM7M1wJ8CH3XO9Y0xTyas44mss8eAkdEWHwf+baxwmwlef8XfA3ucc98eY566kX4NM7uc+P+tbwE2wb/1Y8Dt3qii1UBXwu4Ov4y59yDT1nGCxPfreuDRJPM8BVxvZhXeLufrvbbJ8bP3PFUX4PvA741qmwc87k0vJj6y5FVgF/FdH37W+4/ATmCH98euH12zd3sd8dEl+/ysGWghvk/yFe8yMhon49ZxsnUGfJl4kAGEgH/yXtOLwGKf3wtXE99duCNh/a4Dfm/kPQ1s9Nbpq8Q78K/0ueakf+tRNRvw197fYSew0ueaS4h/uJcntGXUOiYeVK3AMPH9/ncQ78/aAuwFfgFUevOuBP4u4bGf9t7TLcCnpvL8OgJZRESyZjeRiIhMg8JAREQUBiIiojAQEREUBiIigsJARERQGIiICAoDEREB/j/suPVRcs3s7AAAAABJRU5ErkJggg==\n"
          },
          "metadata": {
            "needs_background": "light"
          }
        }
      ],
      "source": [
        "import numpy as np\n",
        "import matplotlib.pyplot as plt\n",
        "\n",
        "data_x = np.array([-5.,-1,-0.5,0.5,3.5,7.,9])\n",
        "data_y = np.array([-30., -50., 36., 31., 75., 60., 46.])\n",
        "differences = calc_div_diff(list(data_x), list(data_y))\n",
        "p = Newtons_Divided_Differences(list(differences))\n",
        "test_x = np.linspace(-7, 10, 50, endpoint=True)\n",
        "test_y = p(test_x)\n",
        "\n",
        "#generating 50 points from -3 to 4 in order to create a smooth line\n",
        "plt.plot(test_x, test_y)\n",
        "plt.plot(data_x, data_y, 'ro')\n",
        "plt.show()\n"
      ]
    },
    {
      "cell_type": "code",
      "source": [
        "# ans hobe [-0.86]"
      ],
      "metadata": {
        "id": "3akMZQcpq2e3"
      },
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "source": [
        "**Problem related Newton's Divided Difference interpolation**\n",
        "\n",
        "Suppose, you have three nodes (-0.5, 1.87), (0, 2.20), (0.5, 2.44). Using Newton's Divided Difference method, print out the value of the interpolating polynomial at x = 6.\n",
        "\n",
        "You have to solve the problem using Newtons_Divided_Differences class."
      ],
      "metadata": {
        "id": "Cy_hrOvvz7B2"
      }
    },
    {
      "cell_type": "code",
      "source": [
        "\n",
        "data_x = np.array([-0.5, 0, 0.5])\n",
        "data_y = np.array([1.87, 2.20, 2.44])\n",
        "differences = calc_div_diff(list(data_x), list(data_y))\n",
        "p = Newtons_Divided_Differences(list(differences))\n",
        "temp = p([6])\n",
        "print(temp)"
      ],
      "metadata": {
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "id": "XdjKvzeh1sY3",
        "outputId": "bd986212-2732-4543-d2ec-7fb576178925"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "[-0.86]\n"
          ]
        }
      ]
    },
    {
      "cell_type": "markdown",
      "source": [
        "**Problem related Hermite interpolation**\n",
        "\n",
        "Suppose, consider the following data set:\n",
        "\n",
        "![Capture.JPG](data:image/jpeg;base64,/9j/4AAQSkZJRgABAQEAeAB4AAD/4RDaRXhpZgAATU0AKgAAAAgABAE7AAIAAAAFAAAISodpAAQAAAABAAAIUJydAAEAAAAKAAAQyOocAAcAAAgMAAAAPgAAAAAc6gAAAAgAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAERlbGwAAAAFkAMAAgAAABQAABCekAQAAgAAABQAABCykpEAAgAAAAMxOAAAkpIAAgAAAAMxOAAA6hwABwAACAwAAAiSAAAAABzqAAAACAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAMjAyMjoxMDoxNCAxOTozNjoyNQAyMDIyOjEwOjE0IDE5OjM2OjI1AAAARABlAGwAbAAAAP/hCxdodHRwOi8vbnMuYWRvYmUuY29tL3hhcC8xLjAvADw/eHBhY2tldCBiZWdpbj0n77u/JyBpZD0nVzVNME1wQ2VoaUh6cmVTek5UY3prYzlkJz8+DQo8eDp4bXBtZXRhIHhtbG5zOng9ImFkb2JlOm5zOm1ldGEvIj48cmRmOlJERiB4bWxuczpyZGY9Imh0dHA6Ly93d3cudzMub3JnLzE5OTkvMDIvMjItcmRmLXN5bnRheC1ucyMiPjxyZGY6RGVzY3JpcHRpb24gcmRmOmFib3V0PSJ1dWlkOmZhZjViZGQ1LWJhM2QtMTFkYS1hZDMxLWQzM2Q3NTE4MmYxYiIgeG1sbnM6ZGM9Imh0dHA6Ly9wdXJsLm9yZy9kYy9lbGVtZW50cy8xLjEvIi8+PHJkZjpEZXNjcmlwdGlvbiByZGY6YWJvdXQ9InV1aWQ6ZmFmNWJkZDUtYmEzZC0xMWRhLWFkMzEtZDMzZDc1MTgyZjFiIiB4bWxuczp4bXA9Imh0dHA6Ly9ucy5hZG9iZS5jb20veGFwLzEuMC8iPjx4bXA6Q3JlYXRlRGF0ZT4yMDIyLTEwLTE0VDE5OjM2OjI1LjE4MDwveG1wOkNyZWF0ZURhdGU+PC9yZGY6RGVzY3JpcHRpb24+PHJkZjpEZXNjcmlwdGlvbiByZGY6YWJvdXQ9InV1aWQ6ZmFmNWJkZDUtYmEzZC0xMWRhLWFkMzEtZDMzZDc1MTgyZjFiIiB4bWxuczpkYz0iaHR0cDovL3B1cmwub3JnL2RjL2VsZW1lbnRzLzEuMS8iPjxkYzpjcmVhdG9yPjxyZGY6U2VxIHhtbG5zOnJkZj0iaHR0cDovL3d3dy53My5vcmcvMTk5OS8wMi8yMi1yZGYtc3ludGF4LW5zIyI+PHJkZjpsaT5EZWxsPC9yZGY6bGk+PC9yZGY6U2VxPg0KCQkJPC9kYzpjcmVhdG9yPjwvcmRmOkRlc2NyaXB0aW9uPjwvcmRmOlJERj48L3g6eG1wbWV0YT4NCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgPD94cGFja2V0IGVuZD0ndyc/Pv/bAEMABwUFBgUEBwYFBggHBwgKEQsKCQkKFQ8QDBEYFRoZGBUYFxseJyEbHSUdFxgiLiIlKCkrLCsaIC8zLyoyJyorKv/bAEMBBwgICgkKFAsLFCocGBwqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKioqKv/AABEIAFoBZgMBIgACEQEDEQH/xAAfAAABBQEBAQEBAQAAAAAAAAAAAQIDBAUGBwgJCgv/xAC1EAACAQMDAgQDBQUEBAAAAX0BAgMABBEFEiExQQYTUWEHInEUMoGRoQgjQrHBFVLR8CQzYnKCCQoWFxgZGiUmJygpKjQ1Njc4OTpDREVGR0hJSlNUVVZXWFlaY2RlZmdoaWpzdHV2d3h5eoOEhYaHiImKkpOUlZaXmJmaoqOkpaanqKmqsrO0tba3uLm6wsPExcbHyMnK0tPU1dbX2Nna4eLj5OXm5+jp6vHy8/T19vf4+fr/xAAfAQADAQEBAQEBAQEBAAAAAAAAAQIDBAUGBwgJCgv/xAC1EQACAQIEBAMEBwUEBAABAncAAQIDEQQFITEGEkFRB2FxEyIygQgUQpGhscEJIzNS8BVictEKFiQ04SXxFxgZGiYnKCkqNTY3ODk6Q0RFRkdISUpTVFVWV1hZWmNkZWZnaGlqc3R1dnd4eXqCg4SFhoeIiYqSk5SVlpeYmZqio6Slpqeoqaqys7S1tre4ubrCw8TFxsfIycrS09TV1tfY2dri4+Tl5ufo6ery8/T19vf4+fr/2gAMAwEAAhEDEQA/APpGiq19f2+m23n3bssZdYxtRnJZmCqAFBJJJAqv/btr/wA8b/8A8F8//wARQBo0Vnf27a/88b//AMF8/wD8RR/btr/zxv8A/wAF8/8A8RQBo0Vnf27a/wDPG/8A/BfP/wDEUf27a/8APG//APBfP/8AEUAaNFZ39u2v/PG//wDBfP8A/EUf27a/88b/AP8ABfP/APEUAaNFZ39u2v8Azxv/APwXz/8AxFH9u2v/ADxv/wDwXz//ABFAGjRWd/btr/zxv/8AwXz/APxFH9u2v/PG/wD/AAXz/wDxFAGjRWVN4jsLeLzJ0vo03BdzafPjJIAH3O5IFSf27a/88b//AMF8/wD8RQBo0Vnf27a/88b/AP8ABfP/APEUf27a/wDPG/8A/BfP/wDEUAaNFZ39u2v/ADxv/wDwXz//ABFH9u2v/PG//wDBfP8A/EUAaNFZ39u2v/PG/wD/AAXz/wDxFH9u2v8Azxv/APwXz/8AxFAGjRWd/btr/wA8b/8A8F8//wARR/btr/zxv/8AwXz/APxFAGjRWd/btr/zxv8A/wAF8/8A8RR/btr/AM8b/wD8F8//AMRQBo0Vnf27a/8APG//APBfP/8AEUf27a/88b//AMF8/wD8RQBo0Vnf27a/88b/AP8ABfP/APEVHN4jsYEDzJfIpZUBOnz/AHmIUD7nckD8aANWis7+3bX/AJ43/wD4L5//AIij+3bX/njf/wDgvn/+IoA0aKzv7dtf+eN//wCC+f8A+Io/t21/543/AP4L5/8A4igDRorO/t21/wCeN/8A+C+f/wCIo/t21/543/8A4L5//iKANGis067aDrDfj/uHz/8AxFL/AG7a/wDPG/8A/BfP/wDEUAaNFZ39u2v/ADxv/wDwXz//ABFH9u2v/PG//wDBfP8A/EUAaNFZ39u2v/PG/wD/AAXz/wDxFH9u2v8Azxv/APwXz/8AxFAGjRWd/btr/wA8b/8A8F8//wARR/btr/zxv/8AwXz/APxFAGjRWd/btr/zxv8A/wAF8/8A8RU1hqlrqTTratJvt3EcqSRPGyMVDDIYA9CDQBbooooAxfFP/IPsv+wnZ/8Ao9K2qxfFP/IPsv8AsJ2f/o9K2qACiiigAooooAydY8SWOhqz33m+VHs86SNNwhDttUt35bjAye+Mc1a0zUk1S2eVILiAxytE8dxHsYMpwfYj3BIrD8S6RruszC3i/s86dFdW10sbs6yS+U4dkY4IAJAIIH8OCDuyOitBcC1X7ayGc5LiPlVyc7QcDIHTOBnGcCgCaiiigAooooAyfEv/ACBD/wBfNv8A+j0rWrJ8S/8AIEP/AF82/wD6PStagAooooAKKKKAMvxLrS+HfDOoau8LzrZwPMyJjOFBPcjjiuYsrrUdE+Hp1K8vNQfUtQFtCqXrKwguJdkWUHOFLtvxn8hxVzx7d29/pjeFomZ77VGhiMKox/cPKBI2cY4RZDjOeDWh4006S+8LuLZGeSzngvUjQZL+TKku0D1IQge5oAsXV/Z+GNJCuk8sVvC0shQbmEa8vKxJGfU9yTwCarz3yeItHvf7OFzbXFrh7eaWIp8+zejrnqpBGR6EggdKq+J7DV/E+jT2ejXFgNM1Cy2+c5fzDk5IBGRsZfl9RuJ5xg6V3e3GmeH7y81BYw6qfKt4PmwcBUjBwCzM2Ow5bFAFjQdWj13w7p+rQKVjvraOdVP8IZQcfrV+srwtpLaD4R0nSZGDvZWcUDsOhZVAJ/OtWgAooooAKy/EH/IMi/6/rT/0pjrUrL8Qf8gyL/r+tP8A0pjoA1KKKKACiiigArC8a67N4b8GapqttA809vbO0SqAcPtO0nJHGcVu1yPjS7i1JrTw1aky3dze2j3EQU/LbiQyMxOMYIhYfiPUUAVYfDd/qGk6PaX9xqU0T3sd9em9lUtH5SgrGMHO1pVVsc8bumQK2oLmTUvGt7C+Ra6TFGEAYgPNICWJHfam3H++3tW9XN2qT2HjbWY/3ajU4Irm1Zzw0iL5bqfoBGfox9KAFufG+m20xiNveynyRcxmCDzPNh3bTIuDkqCR6E5GAamMzad41htlb/RdUt5JNh/gmiK8j/eV+f8Acz3NZem6H4kXxNJrGpPpqzXNjBZu1s74txGzs5RWXneWGASMYGd2OdK5U33jyxEa7o9MtZZJn7LJKVVF+u1XJ9Pl9aAN+iiigAooooAKwtE/5GbxJ/19Q/8ApPHW7WFon/IzeJP+vqH/ANJ46AN2iiigDC8Xq76PbpFJ5UjahaBJMZ2Hz0wcHr9Kd/Zeuf8AQxN/4BR0vin/AJB9l/2E7P8A9HpW1QBif2Xrn/QxN/4BR0f2Xrn/AEMTf+AUdbdFAGJ/Zeuf9DE3/gFHVe+i1PTrfzrzxOyIWCqBYIzOx6KqjJY+wBNdHXIeKWmHj/wSBu+zm7ut/pv+zPt/HG+gC5YxanqNv59n4oMibirf6CgZGHVWU8qR3BAIqKaW6t9Ut9Nm8YRLe3LFYrf7LGXYhS547fKCcmoPCchTxZ40LMFtv7UhCEnA3/ZYdw+ucUeKv+R/8Ef9fd1/6SvQBqf2Xrn/AEMTf+AUdH9l65/0MTf+AUdbdFAGJ/Zeuf8AQxN/4BR0f2Xrn/QxN/4BR1t0UAcl4g03WU0cmXX2kXz4Bt+xxjnzkwfwPNaX9l65/wBDE3/gFHUviX/kCH/r5t//AEela1AGJ/Zeuf8AQxN/4BR0f2Xrn/QxN/4BR1t0UAYn9l65/wBDE3/gFHR/Zeuf9DE3/gFHW3VPUr+TT7dZYtPu78s23y7UIWHHU7mUY/GgDB1D7Zps0X2zxM3nurGNI9NWWTaMbm2qCdo4yegyM9quQ2Or3MEc1v4nEsUih0dLOMqwPIII6isvTbmS6+LNzLcQS22/QIGihnA3p+/l3jgkZ+5nBPak+Fkkg+Hdis7YZprryFY8mMTybMe23H4YoAkE8ttqCaanikRytJ5YVNOXyxIctsLgbQ55O0nJ9KSW6dtVh0+fxWv2hpNsPmacoR5Bn5UkI2lxz8oO4c8VyPhTULrSvh34W1HUIo7yW71wxXFvJErFJZZ5VMinGRIrHJPpke9dFrWh2l19g0TT5Gg02w1RNQvbqaVpCsgk8xYYySTuZ2HThV47igDRa6uV1D7EfFn74OIif7PXYsh6IXxtDnI+UnPI45rQ/svXP+hib/wCjrgb1ro/BXXCm7+0f7anxt+9539o/Jj3+7j8K9ZoAxP7L1z/AKGJv/AKOj+y9c/6GJv/AACjrbooAxP7L1z/AKGJv/AKOs7XNM1ldPiMmvs4+2WoA+xxjBM8YB/A4P4V1lZfiD/kGRf9f1p/6Ux0AQf2Xrn/AEMTf+AUdH9l65/0MTf+AUdbdFAGJ/Zeuf8AQxN/4BR0f2Xrn/QxN/4BR1t0UAYn9l65/wBDE3/gFHWZqVzNpE27UfFXlMqZkl/s1WWFCfvSMARGvB5YgcH0rrq4Dx0dPv8AR/FSaXLafa7Sx2a2ghKzy2/ls6xrL0QlS+CVcDPQUAdGum606B08SblYZDCzjIIrD0jWbbxLcpHovja11KaLL7YLON2h7Zb+4eSOcE810Gkarp//AAiFjqEAa3sv7OjuUibl44fLDDIHJwOK52DWz4efwrYrN9s07VLOXe5g8p0EcAkVwvUAgEENk5I545ANCC4u7m9FrF4qfzWz5ZbTlVJcDJ2ORtfAHO0nFUdN1e2v9Xm03SvGtrPqAmdJ7aCzjMiOnDF1HK/dxubAPHPSsLw9ba5pFz4aOpf8THwzJIH0s7wLnTy8TCNZsDEqhGKgjoTzuwDW1H4h/su28NXelz/brHXtUaAh7cxMEk811dQQGyu0ZznIyeOKAOg/svXP+hib/wAAo6P7L1z/AKGJv/AKOtuigDE/svXP+hib/wAAo6P7L1z/AKGJv/AKOtuigDE/svXP+hib/wAAo6r+GIpoNZ8Qx3dx9pmW7i3TFAm7/R48cDgYro6wtE/5GbxJ/wBfUP8A6Tx0AbtFFFAGH4tkSLSrWWVlSOPUbRndjgKBOmST2FW/+Ei0X/oMWH/gUn+NVfFQDadZhhkHUrQEHv8Av0rW+zQf88I/++BQBT/4SLRf+gxYf+BSf40f8JFov/QYsP8AwKT/ABq59mg/54R/98Cj7NB/zwj/AO+BQBT/AOEi0X/oMWH/AIFJ/jVe81Pw5qESR3ep6e4RxIhF2qsjDowIbIPuK1Ps0H/PCP8A74FYuualFp+qaPpltBb/AGrVbh40aSPcqIkbSO2ARnhQBz1bPagCNj4OeyFpLcaTLB9oF0VkuEfMwbcJCSclsgHJ5p15N4V1C+t7y71GxkuLUkwSfbgDESMErhuCRwcdRT9B1FNTvNXsbq2t1udLvPs7mNMK6tGsiMAc4+VwCMnkGsbxL4jv9G0XxDq0WlQRW+iMCsd3bkfbUCqzGOQNgfeKj5TyvvwAdN/wkWif9Biw/wDApP8AGj/hItF/6DFh/wCBSf41Ygjt57aOYW6KJEDBWQZGRnBqT7NB/wA8I/8AvgUAU/8AhItF/wCgxYf+BSf40f8ACRaL/wBBiw/8Ck/xq59mg/54R/8AfAo+zQf88I/++BQBz/iLXtIl0YrHqti7faIDhblCcCZCe/pWp/wkWi/9Biw/8Ck/xqt4lt4BopxDGP8ASLf+Ef8APZK1fs0H/PCP/vgUAU/+Ei0X/oMWH/gUn+NH/CRaL/0GLD/wKT/Grn2aD/nhH/3wKPs0H/PCP/vgUAU/+Ei0X/oMWH/gUn+NH/CRaL/0GLD/AMCk/wAaufZoP+eEf/fAo+zQf88I/wDvgUAY15eeGL64juLjUrHz4kaNJo70Ruqtjcu5WBwcDI6cD0pFn8JLcWM63OkiXTo2itHE0eYEIClV54BAA/Cm3F+0vjA6Fp8FpG0VgLyaaaHfjc5SNQoI/uOSc9h61L4V1W18TeGrXVEtI4ml3JJHtB2SIxRwPbcpx7UAQRN4QgvFuYrzTVkWZp0Au12JK2dzhN20MdzZYDJ3H1NZa+GvhqmuDWVttC/tIT/aPtRmQv5mc785655plx4nvbKPQb280qKG31bUxp72U1qY7mEsXCPncQwGzJ4GVOe2KmtfFtlfahJa2j6fJeW+pfYp9M8vFwi7yu/k9NoMmduNoIz3oA0i/hI3jXJvdO8xphcMBeKEaUAASFN20sMDnGeB6Vo/8JFov/QYsP8AwKT/ABrnNX13U9KuLXzNItVF1rKWEFsUDSTwN1mVlY4wu5yGUYCnPrXX/ZoP+eEf/fAoAp/8JFov/QYsP/ApP8aP+Ei0X/oMWH/gUn+NXPs0H/PCP/vgUfZoP+eEf/fAoAp/8JFov/QYsP8AwKT/ABrN17XtHk06IR6rYsRe2rELcoeBcRknr2AJre+zQf8APCP/AL4FZniC3gGmxYhjH+nWn8I/5+I6AJ/+Ei0X/oMWH/gUn+NH/CRaL/0GLD/wKT/Grn2aD/nhH/3wKPs0H/PCP/vgUAU/+Ei0X/oMWH/gUn+NH/CRaL/0GLD/AMCk/wAaufZoP+eEf/fAo+zQf88I/wDvgUAU/wDhItF/6DFh/wCBSf41nag/hLVFmW+vtPlS4XZOn20Ksy9MOoYBxjjBzxxW79mg/wCeEf8A3wK5fxXq2oaFpus39tploLews1mtnmQP9rlJb90ArBlOQoHByXGM9KANAXfhRdRS/W70kXaW/wBlWYTR7liznYDn7uRnFR203hS0ukuYb/T/ADo4zFGz3gfykOMqgZjtB2jgYHA9K1rNEuLCCaaySCWSJXeJlGY2IyVP06VieH9SudR8Q67pupWdjGdMkgVDbqTkSRhzknrjOOgoAktZPCdm0Jt7/T1FvzBGb0FIOCPkQttTgkfKBwSKWOXwnFeQ3S32nmWAsYd14rLCWGCUUthCQSOAOCRUL64n/CcQaRFZ25smsrmaS4K/MZIniUqO2B5hyeeRjjBrNs/Fxn0bw9rclla/YNcvVt0jEWJIVk3eUxOSCchQRgfe9uQDqP8AhItF/wCgxYf+BSf40f8ACRaL/wBBiw/8Ck/xq59mg/54R/8AfAo+zQf88I/++BQBT/4SLRf+gxYf+BSf40f8JFov/QYsP/ApP8aufZoP+eEf/fAo+zQf88I/++BQBT/4SLRf+gxYf+BSf41n+HLiG717xFPaypNC93FtkjYMrYt4wcEcVufZoP8AnhH/AN8CsbQlVPEniRUUKou4cADA/wCPeOgDeooooAxfFP8AyD7L/sJ2f/o9K2qy/EFjdX+molgIWuIbmG4RZ3KI3lyK5BYBiMgdcGoPtPir/oE6P/4NJf8A5HoA26KxPtPir/oE6P8A+DSX/wCR6PtPir/oE6P/AODSX/5HoA265rxNpM8/iDw7rdtG8w0m4lM0UYyxjkiZCyjuQdpx1IzjJwDa+0+Kv+gTo/8A4NJf/kej7T4q/wCgTo//AINJf/kegDN0Wz1LTr3XdZXTZJpNY1OJ47cyrGyQLHHF5jbuhAVn29cYHXioPE82q3urfZn8I6nqmm2rLLEIbi0SK5lHKlw86ttU/wAJXkjPOBWz9p8Vf9AnR/8AwaS//I9H2nxV/wBAnR//AAaS/wDyPQBsxFzCnmgB9o3AdAe9OrE+0+Kv+gTo/wD4NJf/AJHo+0+Kv+gTo/8A4NJf/kegDborE+0+Kv8AoE6P/wCDSX/5Ho+0+Kv+gTo//g0l/wDkegCXxL/yBD/182//AKPStauZ1OPxVqNibf8AszR4/wB5G+7+05T9x1bH/Hv324q39p8Vf9AnR/8AwaS//I9AG3RWJ9p8Vf8AQJ0f/wAGkv8A8j0fafFX/QJ0f/waS/8AyPQBt1R1bT59StVittVvNMZX3GWzERZhg/KfMRxjn0zx1ql9p8Vf9AnR/wDwaS//ACPR9p8Vf9AnR/8AwaS//I9AGWml3uheNP7Wf7ZqtvPpSWck+1GmMscrupZVCjDCQjIAA284zmn+DNL1LwxoOk6RNYmbzTcT3lwsy7bZ3cyBMdW5crkf3c960ftPir/oE6P/AODSX/5Ho+0+Kv8AoE6P/wCDSX/5HoAwUm1u/wDFcF5feENQR4J/KtZ7i4tGt7WEth5QFmLl2Qf3eMhRxuLU7jwk17d2l22n3EHim1vRKdXQ4V4/NG4FwfmjMY2iM8j0HWuq+0+Kv+gTo/8A4NJf/kej7T4q/wCgTo//AINJf/kegDmfEtj4k1mzS2jso4Ndi1RXs9St49sdvarKDuZyxJJjBDIOpb7uOa9BrE+0+Kv+gTo//g0l/wDkej7T4q/6BOj/APg0l/8AkegDborE+0+Kv+gTo/8A4NJf/kej7T4q/wCgTo//AINJf/kegDbrL8Qf8gyL/r+tP/SmOoPtPir/AKBOj/8Ag0l/+R6q6gniq+tUi/szR02zwzZ/tOU/6uRXx/x799uPxoA6WisT7T4q/wCgTo//AINJf/kej7T4q/6BOj/+DSX/AOR6ANuisT7T4q/6BOj/APg0l/8Akej7T4q/6BOj/wDg0l/+R6ANuuK8XwazqNhr+nTaXHeJPb7dFktk/eJKUILPIWwhD4IYYwB1Jra+0+Kv+gTo/wD4NJf/AJHo+0+Kv+gTo/8A4NJf/kegB+inVrS203TtUt2uHj06P7VqIlXa9wAFZdv3jnls4xWRodlqqeL/ABRc3GnXNjBqZgNtcs8LbdkIQkqrk5yOOK1PtPir/oE6P/4NJf8A5Ho+0+Kv+gTo/wD4NJf/AJHoA506Dr0HjvSmVZJ9LttPntZLoRxAAyPEcYMu852HLYzn1zVTTvDGpr4Z8JeGZbZ1fQr+Ga6uSMRtHBuKMh/iLHZx1GTnGOet+0+Kv+gTo/8A4NJf/kej7T4q/wCgTo//AINJf/kegDborE+0+Kv+gTo//g0l/wDkej7T4q/6BOj/APg0l/8AkegDborE+0+Kv+gTo/8A4NJf/kej7T4q/wCgTo//AINJf/kegDbrC0T/AJGbxJ/19Q/+k8dO+0+Kv+gTo/8A4NJf/kenaFYX9vd6peaqltFNfTpIIraZpVRViRPvMqnJ2k9KANmiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKAP/2Q==)\n",
        "\n",
        "Using Hermit basis, print out the interpolating polynomial and find the value at x = [0.15,0.30,0.50].\n",
        "\n",
        "You have to solve this problem using hermit function."
      ],
      "metadata": {
        "id": "1FlrGq6b0JSk"
      }
    },
    {
      "cell_type": "code",
      "source": [
        "x = np.array([0.1, 0.2])\n",
        "y = np.array([-0.620, -0.283])\n",
        "y_prime = np.array([3.585, 3.140])\n",
        "h8 = hermit(x, y, y_prime)\n",
        "print(h8(0.15))\n",
        "print(h8(0.3))\n",
        "print(h8(0.5))"
      ],
      "metadata": {
        "id": "PVqdwknilBpq",
        "colab": {
          "base_uri": "https://localhost:8080/"
        },
        "outputId": "52a25421-0e58-47eb-9314-2e8ec4c8eaa9"
      },
      "execution_count": null,
      "outputs": [
        {
          "output_type": "stream",
          "name": "stdout",
          "text": [
            "0.005562500000001691\n",
            "0.4370000000000047\n",
            "35.342000000000006\n"
          ]
        }
      ]
    }
  ],
  "metadata": {
    "colab": {
      "provenance": [],
      "include_colab_link": true
    },
    "kernelspec": {
      "display_name": "Python 3",
      "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.9.1"
    }
  },
  "nbformat": 4,
  "nbformat_minor": 0
}