Example-2.1-case-1.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "b16ccb14",
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "from utils import *\n",
    "%matplotlib notebook"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "6a06026c",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "r = k/h**2 = 0.1\n",
      "Number of discretization on x: 11\n",
      "Number of discretization on t: 101\n"
     ]
    }
   ],
   "source": [
    "L = 1\n",
    "T = .1\n",
    "\n",
    "h = .1\n",
    "k = 0.001\n",
    "\n",
    "r = k / h**2\n",
    "print('r = k/h**2 = %g' % r)\n",
    "\n",
    "Nx = int(np.ceil(L / h)) + 1\n",
    "Nt = int(np.ceil(T / k)) + 1\n",
    "\n",
    "print('Number of discretization on x:', Nx)\n",
    "print('Number of discretization on t:', Nt)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "ae6a00e3",
   "metadata": {},
   "outputs": [],
   "source": [
    "def exact(x, t, N=20):\n",
    "    s = sum(1 / n**2 * np.sin(np.pi / 2 * n) * np.sin(np.pi * x * n) *\n",
    "            np.exp(-n**2 * np.pi**2 * t) for n in range(1, N + 1))\n",
    "    return 8 / np.pi**2 * s"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "3e255551",
   "metadata": {},
   "outputs": [],
   "source": [
    "@np.vectorize\n",
    "def initial(x):\n",
    "    if 0 <= x <= 1 / 2:\n",
    "        return 2 * x\n",
    "    elif 1 / 2 <= x <= 1:\n",
    "        return 2 * (1 - x)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "1093e809",
   "metadata": {},
   "outputs": [],
   "source": [
    "x = np.array([i * h for i in range(Nx)])\n",
    "t = np.array([i * k for i in range(Nt)])\n",
    "\n",
    "u0 = initial(x)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "2e15648b",
   "metadata": {},
   "outputs": [],
   "source": [
    "u = np.zeros((Nt, Nx))\n",
    "u[0, :] = u0\n",
    "\n",
    "# for j in range(Nt - 1):\n",
    "#     for i in range(1, Nx - 1):\n",
    "#         u[j + 1, i] = r * u[j, i - 1] + (1 - 2 * r) * u[j, i] + r * u[j, i + 1]\n",
    "\n",
    "# vectorized form\n",
    "for j in range(Nt - 1):\n",
    "    u[j + 1, 1:-1] = r * u[j, :-2] + (1 - 2 * r) * u[j, 1:-1] + r * u[j, 2:]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "be837689",
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "╒═══════════════╤═══════╤════════╤════════╤════════╤════════╤════════╤════════╤════════╤════════╤════════╤═══════╕\n",
      "│       x       │  0.0  │  0.1   │  0.2   │  0.3   │  0.4   │  0.5   │  0.6   │  0.7   │  0.8   │  0.9   │  1.0  │\n",
      "╞═══════════════╪═══════╪════════╪════════╪════════╪════════╪════════╪════════╪════════╪════════╪════════╪═══════╡\n",
      "│   j=0(t=0)    │   0   │  0.2   │  0.4   │  0.6   │  0.8   │   1    │  0.8   │  0.6   │  0.4   │  0.2   │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=1(t=0.001)  │   0   │  0.2   │  0.4   │  0.6   │  0.8   │  0.96  │  0.8   │  0.6   │  0.4   │  0.2   │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=2(t=0.002)  │   0   │  0.2   │  0.4   │  0.6   │ 0.796  │ 0.928  │ 0.796  │  0.6   │  0.4   │  0.2   │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=3(t=0.003)  │   0   │  0.2   │  0.4   │ 0.5996 │ 0.7896 │ 0.9016 │ 0.7896 │ 0.5996 │  0.4   │  0.2   │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=4(t=0.004)  │   0   │  0.2   │  0.4   │ 0.5986 │ 0.7818 │ 0.8792 │ 0.7818 │ 0.5986 │  0.4   │  0.2   │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=5(t=0.005)  │   0   │  0.2   │ 0.3998 │ 0.5971 │ 0.7732 │ 0.8597 │ 0.7732 │ 0.5971 │ 0.3998 │  0.2   │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=6(t=0.006)  │   0   │  0.2   │ 0.3996 │ 0.595  │ 0.7643 │ 0.8424 │ 0.7643 │ 0.595  │ 0.3996 │  0.2   │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=7(t=0.007)  │   0   │ 0.1999 │ 0.3992 │ 0.5924 │ 0.7551 │ 0.8268 │ 0.7551 │ 0.5924 │ 0.3992 │ 0.1999 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=8(t=0.008)  │   0   │ 0.1999 │ 0.3986 │ 0.5893 │ 0.746  │ 0.8125 │ 0.746  │ 0.5893 │ 0.3986 │ 0.1999 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=9(t=0.009)  │   0   │ 0.1998 │ 0.3978 │ 0.5859 │ 0.737  │ 0.7992 │ 0.737  │ 0.5859 │ 0.3978 │ 0.1998 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=10(t=0.01)  │   0   │ 0.1996 │ 0.3968 │ 0.5822 │ 0.7281 │ 0.7867 │ 0.7281 │ 0.5822 │ 0.3968 │ 0.1996 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=11(t=0.011) │   0   │ 0.1993 │ 0.3956 │ 0.5783 │ 0.7194 │ 0.775  │ 0.7194 │ 0.5783 │ 0.3956 │ 0.1993 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=12(t=0.012) │   0   │ 0.199  │ 0.3942 │ 0.5741 │ 0.7108 │ 0.7639 │ 0.7108 │ 0.5741 │ 0.3942 │ 0.199  │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=13(t=0.013) │   0   │ 0.1986 │ 0.3927 │ 0.5698 │ 0.7025 │ 0.7533 │ 0.7025 │ 0.5698 │ 0.3927 │ 0.1986 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=14(t=0.014) │   0   │ 0.1982 │ 0.391  │ 0.5653 │ 0.6943 │ 0.7431 │ 0.6943 │ 0.5653 │ 0.391  │ 0.1982 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=15(t=0.015) │   0   │ 0.1977 │ 0.3892 │ 0.5608 │ 0.6863 │ 0.7333 │ 0.6863 │ 0.5608 │ 0.3892 │ 0.1977 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=16(t=0.016) │   0   │ 0.197  │ 0.3872 │ 0.5562 │ 0.6784 │ 0.7239 │ 0.6784 │ 0.5562 │ 0.3872 │ 0.197  │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=17(t=0.017) │   0   │ 0.1963 │ 0.3851 │ 0.5515 │ 0.6708 │ 0.7148 │ 0.6708 │ 0.5515 │ 0.3851 │ 0.1963 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=18(t=0.018) │   0   │ 0.1956 │ 0.3828 │ 0.5468 │ 0.6632 │ 0.706  │ 0.6632 │ 0.5468 │ 0.3828 │ 0.1956 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=19(t=0.019) │   0   │ 0.1948 │ 0.3805 │ 0.542  │ 0.6559 │ 0.6975 │ 0.6559 │ 0.542  │ 0.3805 │ 0.1948 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=20(t=0.02)  │   0   │ 0.1939 │ 0.3781 │ 0.5373 │ 0.6487 │ 0.6891 │ 0.6487 │ 0.5373 │ 0.3781 │ 0.1939 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=21(t=0.021) │   0   │ 0.1929 │ 0.3756 │ 0.5325 │ 0.6416 │ 0.681  │ 0.6416 │ 0.5325 │ 0.3756 │ 0.1929 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=22(t=0.022) │   0   │ 0.1919 │ 0.373  │ 0.5277 │ 0.6346 │ 0.6731 │ 0.6346 │ 0.5277 │ 0.373  │ 0.1919 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=23(t=0.023) │   0   │ 0.1908 │ 0.3704 │ 0.5229 │ 0.6278 │ 0.6654 │ 0.6278 │ 0.5229 │ 0.3704 │ 0.1908 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=24(t=0.024) │   0   │ 0.1897 │ 0.3677 │ 0.5182 │ 0.6211 │ 0.6579 │ 0.6211 │ 0.5182 │ 0.3677 │ 0.1897 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=25(t=0.025) │   0   │ 0.1885 │ 0.3649 │ 0.5134 │ 0.6144 │ 0.6505 │ 0.6144 │ 0.5134 │ 0.3649 │ 0.1885 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=26(t=0.026) │   0   │ 0.1873 │ 0.3621 │ 0.5087 │ 0.608  │ 0.6433 │ 0.608  │ 0.5087 │ 0.3621 │ 0.1873 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=27(t=0.027) │   0   │ 0.186  │ 0.3593 │ 0.5039 │ 0.6016 │ 0.6362 │ 0.6016 │ 0.5039 │ 0.3593 │ 0.186  │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=28(t=0.028) │   0   │ 0.1848 │ 0.3564 │ 0.4992 │ 0.5953 │ 0.6293 │ 0.5953 │ 0.4992 │ 0.3564 │ 0.1848 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=29(t=0.029) │   0   │ 0.1835 │ 0.3535 │ 0.4946 │ 0.5891 │ 0.6225 │ 0.5891 │ 0.4946 │ 0.3535 │ 0.1835 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=30(t=0.03)  │   0   │ 0.1821 │ 0.3506 │ 0.4899 │ 0.583  │ 0.6158 │ 0.583  │ 0.4899 │ 0.3506 │ 0.1821 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=31(t=0.031) │   0   │ 0.1808 │ 0.3477 │ 0.4853 │ 0.5769 │ 0.6092 │ 0.5769 │ 0.4853 │ 0.3477 │ 0.1808 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=32(t=0.032) │   0   │ 0.1794 │ 0.3448 │ 0.4807 │ 0.571  │ 0.6028 │ 0.571  │ 0.4807 │ 0.3448 │ 0.1794 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=33(t=0.033) │   0   │ 0.178  │ 0.3418 │ 0.4761 │ 0.5651 │ 0.5964 │ 0.5651 │ 0.4761 │ 0.3418 │ 0.178  │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=34(t=0.034) │   0   │ 0.1766 │ 0.3389 │ 0.4716 │ 0.5594 │ 0.5902 │ 0.5594 │ 0.4716 │ 0.3389 │ 0.1766 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=35(t=0.035) │   0   │ 0.1751 │ 0.3359 │ 0.4671 │ 0.5537 │ 0.584  │ 0.5537 │ 0.4671 │ 0.3359 │ 0.1751 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=36(t=0.036) │   0   │ 0.1737 │ 0.333  │ 0.4626 │ 0.5481 │ 0.5779 │ 0.5481 │ 0.4626 │ 0.333  │ 0.1737 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=37(t=0.037) │   0   │ 0.1723 │  0.33  │ 0.4582 │ 0.5425 │ 0.572  │ 0.5425 │ 0.4582 │  0.33  │ 0.1723 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=38(t=0.038) │   0   │ 0.1708 │ 0.327  │ 0.4538 │ 0.537  │ 0.5661 │ 0.537  │ 0.4538 │ 0.327  │ 0.1708 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=39(t=0.039) │   0   │ 0.1693 │ 0.3241 │ 0.4495 │ 0.5316 │ 0.5603 │ 0.5316 │ 0.4495 │ 0.3241 │ 0.1693 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=40(t=0.04)  │   0   │ 0.1679 │ 0.3212 │ 0.4451 │ 0.5263 │ 0.5545 │ 0.5263 │ 0.4451 │ 0.3212 │ 0.1679 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=41(t=0.041) │   0   │ 0.1664 │ 0.3182 │ 0.4409 │ 0.521  │ 0.5489 │ 0.521  │ 0.4409 │ 0.3182 │ 0.1664 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=42(t=0.042) │   0   │ 0.165  │ 0.3153 │ 0.4366 │ 0.5157 │ 0.5433 │ 0.5157 │ 0.4366 │ 0.3153 │ 0.165  │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=43(t=0.043) │   0   │ 0.1635 │ 0.3124 │ 0.4324 │ 0.5106 │ 0.5378 │ 0.5106 │ 0.4324 │ 0.3124 │ 0.1635 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=44(t=0.044) │   0   │ 0.162  │ 0.3095 │ 0.4282 │ 0.5055 │ 0.5323 │ 0.5055 │ 0.4282 │ 0.3095 │ 0.162  │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=45(t=0.045) │   0   │ 0.1606 │ 0.3066 │ 0.4241 │ 0.5004 │ 0.527  │ 0.5004 │ 0.4241 │ 0.3066 │ 0.1606 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=46(t=0.046) │   0   │ 0.1591 │ 0.3038 │  0.42  │ 0.4955 │ 0.5217 │ 0.4955 │  0.42  │ 0.3038 │ 0.1591 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=47(t=0.047) │   0   │ 0.1577 │ 0.3009 │ 0.4159 │ 0.4905 │ 0.5164 │ 0.4905 │ 0.4159 │ 0.3009 │ 0.1577 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=48(t=0.048) │   0   │ 0.1562 │ 0.2981 │ 0.4119 │ 0.4857 │ 0.5112 │ 0.4857 │ 0.4119 │ 0.2981 │ 0.1562 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=49(t=0.049) │   0   │ 0.1548 │ 0.2953 │ 0.4079 │ 0.4808 │ 0.5061 │ 0.4808 │ 0.4079 │ 0.2953 │ 0.1548 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=50(t=0.05)  │   0   │ 0.1534 │ 0.2925 │ 0.4039 │ 0.4761 │ 0.5011 │ 0.4761 │ 0.4039 │ 0.2925 │ 0.1534 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=51(t=0.051) │   0   │ 0.1519 │ 0.2897 │  0.4   │ 0.4714 │ 0.4961 │ 0.4714 │  0.4   │ 0.2897 │ 0.1519 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=52(t=0.052) │   0   │ 0.1505 │ 0.287  │ 0.3961 │ 0.4667 │ 0.4911 │ 0.4667 │ 0.3961 │ 0.287  │ 0.1505 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=53(t=0.053) │   0   │ 0.1491 │ 0.2842 │ 0.3922 │ 0.4621 │ 0.4862 │ 0.4621 │ 0.3922 │ 0.2842 │ 0.1491 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=54(t=0.054) │   0   │ 0.1477 │ 0.2815 │ 0.3884 │ 0.4575 │ 0.4814 │ 0.4575 │ 0.3884 │ 0.2815 │ 0.1477 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=55(t=0.055) │   0   │ 0.1463 │ 0.2788 │ 0.3846 │ 0.453  │ 0.4766 │ 0.453  │ 0.3846 │ 0.2788 │ 0.1463 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=56(t=0.056) │   0   │ 0.1449 │ 0.2762 │ 0.3809 │ 0.4485 │ 0.4719 │ 0.4485 │ 0.3809 │ 0.2762 │ 0.1449 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=57(t=0.057) │   0   │ 0.1436 │ 0.2735 │ 0.3772 │ 0.4441 │ 0.4672 │ 0.4441 │ 0.3772 │ 0.2735 │ 0.1436 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=58(t=0.058) │   0   │ 0.1422 │ 0.2709 │ 0.3735 │ 0.4397 │ 0.4626 │ 0.4397 │ 0.3735 │ 0.2709 │ 0.1422 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=59(t=0.059) │   0   │ 0.1409 │ 0.2683 │ 0.3699 │ 0.4354 │ 0.458  │ 0.4354 │ 0.3699 │ 0.2683 │ 0.1409 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=60(t=0.06)  │   0   │ 0.1395 │ 0.2657 │ 0.3663 │ 0.4311 │ 0.4535 │ 0.4311 │ 0.3663 │ 0.2657 │ 0.1395 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=61(t=0.061) │   0   │ 0.1382 │ 0.2631 │ 0.3627 │ 0.4269 │ 0.449  │ 0.4269 │ 0.3627 │ 0.2631 │ 0.1382 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=62(t=0.062) │   0   │ 0.1369 │ 0.2606 │ 0.3592 │ 0.4227 │ 0.4446 │ 0.4227 │ 0.3592 │ 0.2606 │ 0.1369 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=63(t=0.063) │   0   │ 0.1355 │ 0.2581 │ 0.3556 │ 0.4185 │ 0.4402 │ 0.4185 │ 0.3556 │ 0.2581 │ 0.1355 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=64(t=0.064) │   0   │ 0.1342 │ 0.2556 │ 0.3522 │ 0.4144 │ 0.4359 │ 0.4144 │ 0.3522 │ 0.2556 │ 0.1342 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=65(t=0.065) │   0   │ 0.133  │ 0.2531 │ 0.3487 │ 0.4103 │ 0.4316 │ 0.4103 │ 0.3487 │ 0.2531 │ 0.133  │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=66(t=0.066) │   0   │ 0.1317 │ 0.2507 │ 0.3453 │ 0.4063 │ 0.4273 │ 0.4063 │ 0.3453 │ 0.2507 │ 0.1317 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=67(t=0.067) │   0   │ 0.1304 │ 0.2482 │ 0.342  │ 0.4023 │ 0.4231 │ 0.4023 │ 0.342  │ 0.2482 │ 0.1304 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=68(t=0.068) │   0   │ 0.1291 │ 0.2458 │ 0.3386 │ 0.3983 │ 0.4189 │ 0.3983 │ 0.3386 │ 0.2458 │ 0.1291 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=69(t=0.069) │   0   │ 0.1279 │ 0.2434 │ 0.3353 │ 0.3944 │ 0.4148 │ 0.3944 │ 0.3353 │ 0.2434 │ 0.1279 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=70(t=0.07)  │   0   │ 0.1267 │ 0.2411 │ 0.332  │ 0.3906 │ 0.4107 │ 0.3906 │ 0.332  │ 0.2411 │ 0.1267 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=71(t=0.071) │   0   │ 0.1254 │ 0.2387 │ 0.3288 │ 0.3867 │ 0.4067 │ 0.3867 │ 0.3288 │ 0.2387 │ 0.1254 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=72(t=0.072) │   0   │ 0.1242 │ 0.2364 │ 0.3256 │ 0.3829 │ 0.4027 │ 0.3829 │ 0.3256 │ 0.2364 │ 0.1242 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=73(t=0.073) │   0   │ 0.123  │ 0.2341 │ 0.3224 │ 0.3792 │ 0.3987 │ 0.3792 │ 0.3224 │ 0.2341 │ 0.123  │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=74(t=0.074) │   0   │ 0.1218 │ 0.2318 │ 0.3192 │ 0.3754 │ 0.3948 │ 0.3754 │ 0.3192 │ 0.2318 │ 0.1218 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=75(t=0.075) │   0   │ 0.1206 │ 0.2296 │ 0.3161 │ 0.3718 │ 0.391  │ 0.3718 │ 0.3161 │ 0.2296 │ 0.1206 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=76(t=0.076) │   0   │ 0.1195 │ 0.2273 │ 0.313  │ 0.3681 │ 0.3871 │ 0.3681 │ 0.313  │ 0.2273 │ 0.1195 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=77(t=0.077) │   0   │ 0.1183 │ 0.2251 │  0.31  │ 0.3645 │ 0.3833 │ 0.3645 │  0.31  │ 0.2251 │ 0.1183 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=78(t=0.078) │   0   │ 0.1172 │ 0.2229 │ 0.3069 │ 0.3609 │ 0.3796 │ 0.3609 │ 0.3069 │ 0.2229 │ 0.1172 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=79(t=0.079) │   0   │ 0.116  │ 0.2207 │ 0.3039 │ 0.3574 │ 0.3758 │ 0.3574 │ 0.3039 │ 0.2207 │ 0.116  │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=80(t=0.08)  │   0   │ 0.1149 │ 0.2186 │ 0.301  │ 0.3539 │ 0.3721 │ 0.3539 │ 0.301  │ 0.2186 │ 0.1149 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=81(t=0.081) │   0   │ 0.1138 │ 0.2165 │ 0.298  │ 0.3504 │ 0.3685 │ 0.3504 │ 0.298  │ 0.2165 │ 0.1138 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=82(t=0.082) │   0   │ 0.1127 │ 0.2143 │ 0.2951 │ 0.347  │ 0.3649 │ 0.347  │ 0.2951 │ 0.2143 │ 0.1127 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=83(t=0.083) │   0   │ 0.1116 │ 0.2123 │ 0.2922 │ 0.3436 │ 0.3613 │ 0.3436 │ 0.2922 │ 0.2123 │ 0.1116 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=84(t=0.084) │   0   │ 0.1105 │ 0.2102 │ 0.2894 │ 0.3402 │ 0.3578 │ 0.3402 │ 0.2894 │ 0.2102 │ 0.1105 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=85(t=0.085) │   0   │ 0.1094 │ 0.2081 │ 0.2865 │ 0.3369 │ 0.3543 │ 0.3369 │ 0.2865 │ 0.2081 │ 0.1094 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=86(t=0.086) │   0   │ 0.1083 │ 0.2061 │ 0.2837 │ 0.3336 │ 0.3508 │ 0.3336 │ 0.2837 │ 0.2061 │ 0.1083 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=87(t=0.087) │   0   │ 0.1073 │ 0.2041 │ 0.2809 │ 0.3303 │ 0.3473 │ 0.3303 │ 0.2809 │ 0.2041 │ 0.1073 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=88(t=0.088) │   0   │ 0.1062 │ 0.2021 │ 0.2782 │ 0.3271 │ 0.3439 │ 0.3271 │ 0.2782 │ 0.2021 │ 0.1062 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=89(t=0.089) │   0   │ 0.1052 │ 0.2001 │ 0.2755 │ 0.3239 │ 0.3406 │ 0.3239 │ 0.2755 │ 0.2001 │ 0.1052 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=90(t=0.09)  │   0   │ 0.1042 │ 0.1982 │ 0.2728 │ 0.3207 │ 0.3372 │ 0.3207 │ 0.2728 │ 0.1982 │ 0.1042 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=91(t=0.091) │   0   │ 0.1031 │ 0.1962 │ 0.2701 │ 0.3176 │ 0.3339 │ 0.3176 │ 0.2701 │ 0.1962 │ 0.1031 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=92(t=0.092) │   0   │ 0.1021 │ 0.1943 │ 0.2675 │ 0.3145 │ 0.3307 │ 0.3145 │ 0.2675 │ 0.1943 │ 0.1021 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=93(t=0.093) │   0   │ 0.1011 │ 0.1924 │ 0.2648 │ 0.3114 │ 0.3274 │ 0.3114 │ 0.2648 │ 0.1924 │ 0.1011 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=94(t=0.094) │   0   │ 0.1002 │ 0.1905 │ 0.2623 │ 0.3083 │ 0.3242 │ 0.3083 │ 0.2623 │ 0.1905 │ 0.1002 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=95(t=0.095) │   0   │ 0.0992 │ 0.1887 │ 0.2597 │ 0.3053 │ 0.321  │ 0.3053 │ 0.2597 │ 0.1887 │ 0.0992 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=96(t=0.096) │   0   │ 0.0982 │ 0.1868 │ 0.2571 │ 0.3023 │ 0.3179 │ 0.3023 │ 0.2571 │ 0.1868 │ 0.0982 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=97(t=0.097) │   0   │ 0.0972 │ 0.185  │ 0.2546 │ 0.2994 │ 0.3148 │ 0.2994 │ 0.2546 │ 0.185  │ 0.0972 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=98(t=0.098) │   0   │ 0.0963 │ 0.1832 │ 0.2521 │ 0.2964 │ 0.3117 │ 0.2964 │ 0.2521 │ 0.1832 │ 0.0963 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=99(t=0.099) │   0   │ 0.0954 │ 0.1814 │ 0.2497 │ 0.2935 │ 0.3086 │ 0.2935 │ 0.2497 │ 0.1814 │ 0.0954 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=100(t=0.1)  │   0   │ 0.0944 │ 0.1796 │ 0.2472 │ 0.2907 │ 0.3056 │ 0.2907 │ 0.2472 │ 0.1796 │ 0.0944 │   0   │\n",
      "├───────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│       x       │   0   │  0.1   │  0.2   │  0.3   │  0.4   │  0.5   │  0.6   │  0.7   │  0.8   │  0.9   │   1   │\n",
      "╘═══════════════╧═══════╧════════╧════════╧════════╧════════╧════════╧════════╧════════╧════════╧════════╧═══════╛\n"
     ]
    }
   ],
   "source": [
    "print_table(u, x, t)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "4be231e8",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "╒═════╤══════════════╤══════════════╤════════╤═══════════╕\n",
      "│  x  │  FD (x=0.3)  │  EX (x=0.3)  │  DIFF  │  Error %  │\n",
      "╞═════╪══════════════╪══════════════╪════════╪═══════════╡\n",
      "│ j=0 │    0.5971    │    0.5966    │ 0.0005 │  0.0812   │\n",
      "├─────┼──────────────┼──────────────┼────────┼───────────┤\n",
      "│ j=1 │    0.5822    │    0.5799    │ 0.0023 │  0.3986   │\n",
      "├─────┼──────────────┼──────────────┼────────┼───────────┤\n",
      "│ j=2 │    0.5373    │    0.5334    │ 0.0039 │  0.7346   │\n",
      "├─────┼──────────────┼──────────────┼────────┼───────────┤\n",
      "│ j=3 │    0.2472    │    0.2444    │ 0.0028 │   1.156   │\n",
      "├─────┼──────────────┼──────────────┼────────┼───────────┤\n",
      "│  x  │  FD (x=0.3)  │  EX (x=0.3)  │  DIFF  │  Error %  │\n",
      "╘═════╧══════════════╧══════════════╧════════╧═══════════╛\n"
     ]
    }
   ],
   "source": [
    "t_indices = [5, 10, 20, 100]\n",
    "x_index = 3\n",
    "u_tmp = u[t_indices, x_index]\n",
    "exact_tmp = exact(x[x_index], t[t_indices])\n",
    "diff_tmp = u_tmp - exact_tmp\n",
    "error_tmp = diff_tmp/exact_tmp * 100\n",
    "\n",
    "header = np.array(['FD (x=0.3)', 'EX (x=0.3)', 'DIFF', 'Error %'])\n",
    "table_2_3 = np.vstack([u_tmp, exact_tmp, diff_tmp, error_tmp]).T\n",
    "\n",
    "print_table(table_2_3, header)\n",
    "# table_2_3"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "66ec12c2",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "/* global mpl */\n",
       "window.mpl = {};\n",
       "\n",
       "mpl.get_websocket_type = function () {\n",
       "    if (typeof WebSocket !== 'undefined') {\n",
       "        return WebSocket;\n",
       "    } else if (typeof MozWebSocket !== 'undefined') {\n",
       "        return MozWebSocket;\n",
       "    } else {\n",
       "        alert(\n",
       "            'Your browser does not have WebSocket support. ' +\n",
       "                'Please try Chrome, Safari or Firefox ≥ 6. ' +\n",
       "                'Firefox 4 and 5 are also supported but you ' +\n",
       "                'have to enable WebSockets in about:config.'\n",
       "        );\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure = function (figure_id, websocket, ondownload, parent_element) {\n",
       "    this.id = figure_id;\n",
       "\n",
       "    this.ws = websocket;\n",
       "\n",
       "    this.supports_binary = this.ws.binaryType !== undefined;\n",
       "\n",
       "    if (!this.supports_binary) {\n",
       "        var warnings = document.getElementById('mpl-warnings');\n",
       "        if (warnings) {\n",
       "            warnings.style.display = 'block';\n",
       "            warnings.textContent =\n",
       "                'This browser does not support binary websocket messages. ' +\n",
       "                'Performance may be slow.';\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.imageObj = new Image();\n",
       "\n",
       "    this.context = undefined;\n",
       "    this.message = undefined;\n",
       "    this.canvas = undefined;\n",
       "    this.rubberband_canvas = undefined;\n",
       "    this.rubberband_context = undefined;\n",
       "    this.format_dropdown = undefined;\n",
       "\n",
       "    this.image_mode = 'full';\n",
       "\n",
       "    this.root = document.createElement('div');\n",
       "    this.root.setAttribute('style', 'display: inline-block');\n",
       "    this._root_extra_style(this.root);\n",
       "\n",
       "    parent_element.appendChild(this.root);\n",
       "\n",
       "    this._init_header(this);\n",
       "    this._init_canvas(this);\n",
       "    this._init_toolbar(this);\n",
       "\n",
       "    var fig = this;\n",
       "\n",
       "    this.waiting = false;\n",
       "\n",
       "    this.ws.onopen = function () {\n",
       "        fig.send_message('supports_binary', { value: fig.supports_binary });\n",
       "        fig.send_message('send_image_mode', {});\n",
       "        if (fig.ratio !== 1) {\n",
       "            fig.send_message('set_dpi_ratio', { dpi_ratio: fig.ratio });\n",
       "        }\n",
       "        fig.send_message('refresh', {});\n",
       "    };\n",
       "\n",
       "    this.imageObj.onload = function () {\n",
       "        if (fig.image_mode === 'full') {\n",
       "            // Full images could contain transparency (where diff images\n",
       "            // almost always do), so we need to clear the canvas so that\n",
       "            // there is no ghosting.\n",
       "            fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
       "        }\n",
       "        fig.context.drawImage(fig.imageObj, 0, 0);\n",
       "    };\n",
       "\n",
       "    this.imageObj.onunload = function () {\n",
       "        fig.ws.close();\n",
       "    };\n",
       "\n",
       "    this.ws.onmessage = this._make_on_message_function(this);\n",
       "\n",
       "    this.ondownload = ondownload;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._init_header = function () {\n",
       "    var titlebar = document.createElement('div');\n",
       "    titlebar.classList =\n",
       "        'ui-dialog-titlebar ui-widget-header ui-corner-all ui-helper-clearfix';\n",
       "    var titletext = document.createElement('div');\n",
       "    titletext.classList = 'ui-dialog-title';\n",
       "    titletext.setAttribute(\n",
       "        'style',\n",
       "        'width: 100%; text-align: center; padding: 3px;'\n",
       "    );\n",
       "    titlebar.appendChild(titletext);\n",
       "    this.root.appendChild(titlebar);\n",
       "    this.header = titletext;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function (_canvas_div) {};\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function (_canvas_div) {};\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = (this.canvas_div = document.createElement('div'));\n",
       "    canvas_div.setAttribute(\n",
       "        'style',\n",
       "        'border: 1px solid #ddd;' +\n",
       "            'box-sizing: content-box;' +\n",
       "            'clear: both;' +\n",
       "            'min-height: 1px;' +\n",
       "            'min-width: 1px;' +\n",
       "            'outline: 0;' +\n",
       "            'overflow: hidden;' +\n",
       "            'position: relative;' +\n",
       "            'resize: both;'\n",
       "    );\n",
       "\n",
       "    function on_keyboard_event_closure(name) {\n",
       "        return function (event) {\n",
       "            return fig.key_event(event, name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    canvas_div.addEventListener(\n",
       "        'keydown',\n",
       "        on_keyboard_event_closure('key_press')\n",
       "    );\n",
       "    canvas_div.addEventListener(\n",
       "        'keyup',\n",
       "        on_keyboard_event_closure('key_release')\n",
       "    );\n",
       "\n",
       "    this._canvas_extra_style(canvas_div);\n",
       "    this.root.appendChild(canvas_div);\n",
       "\n",
       "    var canvas = (this.canvas = document.createElement('canvas'));\n",
       "    canvas.classList.add('mpl-canvas');\n",
       "    canvas.setAttribute('style', 'box-sizing: content-box;');\n",
       "\n",
       "    this.context = canvas.getContext('2d');\n",
       "\n",
       "    var backingStore =\n",
       "        this.context.backingStorePixelRatio ||\n",
       "        this.context.webkitBackingStorePixelRatio ||\n",
       "        this.context.mozBackingStorePixelRatio ||\n",
       "        this.context.msBackingStorePixelRatio ||\n",
       "        this.context.oBackingStorePixelRatio ||\n",
       "        this.context.backingStorePixelRatio ||\n",
       "        1;\n",
       "\n",
       "    this.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
       "\n",
       "    var rubberband_canvas = (this.rubberband_canvas = document.createElement(\n",
       "        'canvas'\n",
       "    ));\n",
       "    rubberband_canvas.setAttribute(\n",
       "        'style',\n",
       "        'box-sizing: content-box; position: absolute; left: 0; top: 0; z-index: 1;'\n",
       "    );\n",
       "\n",
       "    // Apply a ponyfill if ResizeObserver is not implemented by browser.\n",
       "    if (this.ResizeObserver === undefined) {\n",
       "        if (window.ResizeObserver !== undefined) {\n",
       "            this.ResizeObserver = window.ResizeObserver;\n",
       "        } else {\n",
       "            var obs = _JSXTOOLS_RESIZE_OBSERVER({});\n",
       "            this.ResizeObserver = obs.ResizeObserver;\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.resizeObserverInstance = new this.ResizeObserver(function (entries) {\n",
       "        var nentries = entries.length;\n",
       "        for (var i = 0; i < nentries; i++) {\n",
       "            var entry = entries[i];\n",
       "            var width, height;\n",
       "            if (entry.contentBoxSize) {\n",
       "                if (entry.contentBoxSize instanceof Array) {\n",
       "                    // Chrome 84 implements new version of spec.\n",
       "                    width = entry.contentBoxSize[0].inlineSize;\n",
       "                    height = entry.contentBoxSize[0].blockSize;\n",
       "                } else {\n",
       "                    // Firefox implements old version of spec.\n",
       "                    width = entry.contentBoxSize.inlineSize;\n",
       "                    height = entry.contentBoxSize.blockSize;\n",
       "                }\n",
       "            } else {\n",
       "                // Chrome <84 implements even older version of spec.\n",
       "                width = entry.contentRect.width;\n",
       "                height = entry.contentRect.height;\n",
       "            }\n",
       "\n",
       "            // Keep the size of the canvas and rubber band canvas in sync with\n",
       "            // the canvas container.\n",
       "            if (entry.devicePixelContentBoxSize) {\n",
       "                // Chrome 84 implements new version of spec.\n",
       "                canvas.setAttribute(\n",
       "                    'width',\n",
       "                    entry.devicePixelContentBoxSize[0].inlineSize\n",
       "                );\n",
       "                canvas.setAttribute(\n",
       "                    'height',\n",
       "                    entry.devicePixelContentBoxSize[0].blockSize\n",
       "                );\n",
       "            } else {\n",
       "                canvas.setAttribute('width', width * fig.ratio);\n",
       "                canvas.setAttribute('height', height * fig.ratio);\n",
       "            }\n",
       "            canvas.setAttribute(\n",
       "                'style',\n",
       "                'width: ' + width + 'px; height: ' + height + 'px;'\n",
       "            );\n",
       "\n",
       "            rubberband_canvas.setAttribute('width', width);\n",
       "            rubberband_canvas.setAttribute('height', height);\n",
       "\n",
       "            // And update the size in Python. We ignore the initial 0/0 size\n",
       "            // that occurs as the element is placed into the DOM, which should\n",
       "            // otherwise not happen due to the minimum size styling.\n",
       "            if (fig.ws.readyState == 1 && width != 0 && height != 0) {\n",
       "                fig.request_resize(width, height);\n",
       "            }\n",
       "        }\n",
       "    });\n",
       "    this.resizeObserverInstance.observe(canvas_div);\n",
       "\n",
       "    function on_mouse_event_closure(name) {\n",
       "        return function (event) {\n",
       "            return fig.mouse_event(event, name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mousedown',\n",
       "        on_mouse_event_closure('button_press')\n",
       "    );\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseup',\n",
       "        on_mouse_event_closure('button_release')\n",
       "    );\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'dblclick',\n",
       "        on_mouse_event_closure('dblclick')\n",
       "    );\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mousemove',\n",
       "        on_mouse_event_closure('motion_notify')\n",
       "    );\n",
       "\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseenter',\n",
       "        on_mouse_event_closure('figure_enter')\n",
       "    );\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseleave',\n",
       "        on_mouse_event_closure('figure_leave')\n",
       "    );\n",
       "\n",
       "    canvas_div.addEventListener('wheel', function (event) {\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        on_mouse_event_closure('scroll')(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.appendChild(canvas);\n",
       "    canvas_div.appendChild(rubberband_canvas);\n",
       "\n",
       "    this.rubberband_context = rubberband_canvas.getContext('2d');\n",
       "    this.rubberband_context.strokeStyle = '#000000';\n",
       "\n",
       "    this._resize_canvas = function (width, height, forward) {\n",
       "        if (forward) {\n",
       "            canvas_div.style.width = width + 'px';\n",
       "            canvas_div.style.height = height + 'px';\n",
       "        }\n",
       "    };\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    this.rubberband_canvas.addEventListener('contextmenu', function (_e) {\n",
       "        event.preventDefault();\n",
       "        return false;\n",
       "    });\n",
       "\n",
       "    function set_focus() {\n",
       "        canvas.focus();\n",
       "        canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    window.setTimeout(set_focus, 100);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var toolbar = document.createElement('div');\n",
       "    toolbar.classList = 'mpl-toolbar';\n",
       "    this.root.appendChild(toolbar);\n",
       "\n",
       "    function on_click_closure(name) {\n",
       "        return function (_event) {\n",
       "            return fig.toolbar_button_onclick(name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    function on_mouseover_closure(tooltip) {\n",
       "        return function (event) {\n",
       "            if (!event.currentTarget.disabled) {\n",
       "                return fig.toolbar_button_onmouseover(tooltip);\n",
       "            }\n",
       "        };\n",
       "    }\n",
       "\n",
       "    fig.buttons = {};\n",
       "    var buttonGroup = document.createElement('div');\n",
       "    buttonGroup.classList = 'mpl-button-group';\n",
       "    for (var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            /* Instead of a spacer, we start a new button group. */\n",
       "            if (buttonGroup.hasChildNodes()) {\n",
       "                toolbar.appendChild(buttonGroup);\n",
       "            }\n",
       "            buttonGroup = document.createElement('div');\n",
       "            buttonGroup.classList = 'mpl-button-group';\n",
       "            continue;\n",
       "        }\n",
       "\n",
       "        var button = (fig.buttons[name] = document.createElement('button'));\n",
       "        button.classList = 'mpl-widget';\n",
       "        button.setAttribute('role', 'button');\n",
       "        button.setAttribute('aria-disabled', 'false');\n",
       "        button.addEventListener('click', on_click_closure(method_name));\n",
       "        button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n",
       "\n",
       "        var icon_img = document.createElement('img');\n",
       "        icon_img.src = '_images/' + image + '.png';\n",
       "        icon_img.srcset = '_images/' + image + '_large.png 2x';\n",
       "        icon_img.alt = tooltip;\n",
       "        button.appendChild(icon_img);\n",
       "\n",
       "        buttonGroup.appendChild(button);\n",
       "    }\n",
       "\n",
       "    if (buttonGroup.hasChildNodes()) {\n",
       "        toolbar.appendChild(buttonGroup);\n",
       "    }\n",
       "\n",
       "    var fmt_picker = document.createElement('select');\n",
       "    fmt_picker.classList = 'mpl-widget';\n",
       "    toolbar.appendChild(fmt_picker);\n",
       "    this.format_dropdown = fmt_picker;\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = document.createElement('option');\n",
       "        option.selected = fmt === mpl.default_extension;\n",
       "        option.innerHTML = fmt;\n",
       "        fmt_picker.appendChild(option);\n",
       "    }\n",
       "\n",
       "    var status_bar = document.createElement('span');\n",
       "    status_bar.classList = 'mpl-message';\n",
       "    toolbar.appendChild(status_bar);\n",
       "    this.message = status_bar;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.request_resize = function (x_pixels, y_pixels) {\n",
       "    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n",
       "    // which will in turn request a refresh of the image.\n",
       "    this.send_message('resize', { width: x_pixels, height: y_pixels });\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.send_message = function (type, properties) {\n",
       "    properties['type'] = type;\n",
       "    properties['figure_id'] = this.id;\n",
       "    this.ws.send(JSON.stringify(properties));\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.send_draw_message = function () {\n",
       "    if (!this.waiting) {\n",
       "        this.waiting = true;\n",
       "        this.ws.send(JSON.stringify({ type: 'draw', figure_id: this.id }));\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_save = function (fig, _msg) {\n",
       "    var format_dropdown = fig.format_dropdown;\n",
       "    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n",
       "    fig.ondownload(fig, format);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_resize = function (fig, msg) {\n",
       "    var size = msg['size'];\n",
       "    if (size[0] !== fig.canvas.width || size[1] !== fig.canvas.height) {\n",
       "        fig._resize_canvas(size[0], size[1], msg['forward']);\n",
       "        fig.send_message('refresh', {});\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function (fig, msg) {\n",
       "    var x0 = msg['x0'] / fig.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / fig.ratio;\n",
       "    var x1 = msg['x1'] / fig.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / fig.ratio;\n",
       "    x0 = Math.floor(x0) + 0.5;\n",
       "    y0 = Math.floor(y0) + 0.5;\n",
       "    x1 = Math.floor(x1) + 0.5;\n",
       "    y1 = Math.floor(y1) + 0.5;\n",
       "    var min_x = Math.min(x0, x1);\n",
       "    var min_y = Math.min(y0, y1);\n",
       "    var width = Math.abs(x1 - x0);\n",
       "    var height = Math.abs(y1 - y0);\n",
       "\n",
       "    fig.rubberband_context.clearRect(\n",
       "        0,\n",
       "        0,\n",
       "        fig.canvas.width / fig.ratio,\n",
       "        fig.canvas.height / fig.ratio\n",
       "    );\n",
       "\n",
       "    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_figure_label = function (fig, msg) {\n",
       "    // Updates the figure title.\n",
       "    fig.header.textContent = msg['label'];\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_cursor = function (fig, msg) {\n",
       "    var cursor = msg['cursor'];\n",
       "    switch (cursor) {\n",
       "        case 0:\n",
       "            cursor = 'pointer';\n",
       "            break;\n",
       "        case 1:\n",
       "            cursor = 'default';\n",
       "            break;\n",
       "        case 2:\n",
       "            cursor = 'crosshair';\n",
       "            break;\n",
       "        case 3:\n",
       "            cursor = 'move';\n",
       "            break;\n",
       "    }\n",
       "    fig.rubberband_canvas.style.cursor = cursor;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_message = function (fig, msg) {\n",
       "    fig.message.textContent = msg['message'];\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_draw = function (fig, _msg) {\n",
       "    // Request the server to send over a new figure.\n",
       "    fig.send_draw_message();\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_image_mode = function (fig, msg) {\n",
       "    fig.image_mode = msg['mode'];\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_history_buttons = function (fig, msg) {\n",
       "    for (var key in msg) {\n",
       "        if (!(key in fig.buttons)) {\n",
       "            continue;\n",
       "        }\n",
       "        fig.buttons[key].disabled = !msg[key];\n",
       "        fig.buttons[key].setAttribute('aria-disabled', !msg[key]);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_navigate_mode = function (fig, msg) {\n",
       "    if (msg['mode'] === 'PAN') {\n",
       "        fig.buttons['Pan'].classList.add('active');\n",
       "        fig.buttons['Zoom'].classList.remove('active');\n",
       "    } else if (msg['mode'] === 'ZOOM') {\n",
       "        fig.buttons['Pan'].classList.remove('active');\n",
       "        fig.buttons['Zoom'].classList.add('active');\n",
       "    } else {\n",
       "        fig.buttons['Pan'].classList.remove('active');\n",
       "        fig.buttons['Zoom'].classList.remove('active');\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function () {\n",
       "    // Called whenever the canvas gets updated.\n",
       "    this.send_message('ack', {});\n",
       "};\n",
       "\n",
       "// A function to construct a web socket function for onmessage handling.\n",
       "// Called in the figure constructor.\n",
       "mpl.figure.prototype._make_on_message_function = function (fig) {\n",
       "    return function socket_on_message(evt) {\n",
       "        if (evt.data instanceof Blob) {\n",
       "            var img = evt.data;\n",
       "            if (img.type !== 'image/png') {\n",
       "                /* FIXME: We get \"Resource interpreted as Image but\n",
       "                 * transferred with MIME type text/plain:\" errors on\n",
       "                 * Chrome.  But how to set the MIME type?  It doesn't seem\n",
       "                 * to be part of the websocket stream */\n",
       "                img.type = 'image/png';\n",
       "            }\n",
       "\n",
       "            /* Free the memory for the previous frames */\n",
       "            if (fig.imageObj.src) {\n",
       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
       "                    fig.imageObj.src\n",
       "                );\n",
       "            }\n",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                img\n",
       "            );\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        } else if (\n",
       "            typeof evt.data === 'string' &&\n",
       "            evt.data.slice(0, 21) === 'data:image/png;base64'\n",
       "        ) {\n",
       "            fig.imageObj.src = evt.data;\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        var msg = JSON.parse(evt.data);\n",
       "        var msg_type = msg['type'];\n",
       "\n",
       "        // Call the  \"handle_{type}\" callback, which takes\n",
       "        // the figure and JSON message as its only arguments.\n",
       "        try {\n",
       "            var callback = fig['handle_' + msg_type];\n",
       "        } catch (e) {\n",
       "            console.log(\n",
       "                \"No handler for the '\" + msg_type + \"' message type: \",\n",
       "                msg\n",
       "            );\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        if (callback) {\n",
       "            try {\n",
       "                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n",
       "                callback(fig, msg);\n",
       "            } catch (e) {\n",
       "                console.log(\n",
       "                    \"Exception inside the 'handler_\" + msg_type + \"' callback:\",\n",
       "                    e,\n",
       "                    e.stack,\n",
       "                    msg\n",
       "                );\n",
       "            }\n",
       "        }\n",
       "    };\n",
       "};\n",
       "\n",
       "// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n",
       "mpl.findpos = function (e) {\n",
       "    //this section is from http://www.quirksmode.org/js/events_properties.html\n",
       "    var targ;\n",
       "    if (!e) {\n",
       "        e = window.event;\n",
       "    }\n",
       "    if (e.target) {\n",
       "        targ = e.target;\n",
       "    } else if (e.srcElement) {\n",
       "        targ = e.srcElement;\n",
       "    }\n",
       "    if (targ.nodeType === 3) {\n",
       "        // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "    }\n",
       "\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    var boundingRect = targ.getBoundingClientRect();\n",
       "    var x = e.pageX - (boundingRect.left + document.body.scrollLeft);\n",
       "    var y = e.pageY - (boundingRect.top + document.body.scrollTop);\n",
       "\n",
       "    return { x: x, y: y };\n",
       "};\n",
       "\n",
       "/*\n",
       " * return a copy of an object with only non-object keys\n",
       " * we need this to avoid circular references\n",
       " * http://stackoverflow.com/a/24161582/3208463\n",
       " */\n",
       "function simpleKeys(original) {\n",
       "    return Object.keys(original).reduce(function (obj, key) {\n",
       "        if (typeof original[key] !== 'object') {\n",
       "            obj[key] = original[key];\n",
       "        }\n",
       "        return obj;\n",
       "    }, {});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.mouse_event = function (event, name) {\n",
       "    var canvas_pos = mpl.findpos(event);\n",
       "\n",
       "    if (name === 'button_press') {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x * this.ratio;\n",
       "    var y = canvas_pos.y * this.ratio;\n",
       "\n",
       "    this.send_message(name, {\n",
       "        x: x,\n",
       "        y: y,\n",
       "        button: event.button,\n",
       "        step: event.step,\n",
       "        guiEvent: simpleKeys(event),\n",
       "    });\n",
       "\n",
       "    /* This prevents the web browser from automatically changing to\n",
       "     * the text insertion cursor when the button is pressed.  We want\n",
       "     * to control all of the cursor setting manually through the\n",
       "     * 'cursor' event from matplotlib */\n",
       "    event.preventDefault();\n",
       "    return false;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function (_event, _name) {\n",
       "    // Handle any extra behaviour associated with a key event\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.key_event = function (event, name) {\n",
       "    // Prevent repeat events\n",
       "    if (name === 'key_press') {\n",
       "        if (event.key === this._key) {\n",
       "            return;\n",
       "        } else {\n",
       "            this._key = event.key;\n",
       "        }\n",
       "    }\n",
       "    if (name === 'key_release') {\n",
       "        this._key = null;\n",
       "    }\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.key !== 'Control') {\n",
       "        value += 'ctrl+';\n",
       "    }\n",
       "    else if (event.altKey && event.key !== 'Alt') {\n",
       "        value += 'alt+';\n",
       "    }\n",
       "    else if (event.shiftKey && event.key !== 'Shift') {\n",
       "        value += 'shift+';\n",
       "    }\n",
       "\n",
       "    value += 'k' + event.key;\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, { key: value, guiEvent: simpleKeys(event) });\n",
       "    return false;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onclick = function (name) {\n",
       "    if (name === 'download') {\n",
       "        this.handle_save(this, null);\n",
       "    } else {\n",
       "        this.send_message('toolbar_button', { name: name });\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onmouseover = function (tooltip) {\n",
       "    this.message.textContent = tooltip;\n",
       "};\n",
       "\n",
       "///////////////// REMAINING CONTENT GENERATED BY embed_js.py /////////////////\n",
       "// prettier-ignore\n",
       "var _JSXTOOLS_RESIZE_OBSERVER=function(A){var t,i=new WeakMap,n=new WeakMap,a=new WeakMap,r=new WeakMap,o=new Set;function s(e){if(!(this instanceof s))throw new TypeError(\"Constructor requires 'new' operator\");i.set(this,e)}function h(){throw new TypeError(\"Function is not a constructor\")}function c(e,t,i,n){e=0 in arguments?Number(arguments[0]):0,t=1 in arguments?Number(arguments[1]):0,i=2 in arguments?Number(arguments[2]):0,n=3 in arguments?Number(arguments[3]):0,this.right=(this.x=this.left=e)+(this.width=i),this.bottom=(this.y=this.top=t)+(this.height=n),Object.freeze(this)}function d(){t=requestAnimationFrame(d);var s=new WeakMap,p=new Set;o.forEach((function(t){r.get(t).forEach((function(i){var r=t instanceof window.SVGElement,o=a.get(t),d=r?0:parseFloat(o.paddingTop),f=r?0:parseFloat(o.paddingRight),l=r?0:parseFloat(o.paddingBottom),u=r?0:parseFloat(o.paddingLeft),g=r?0:parseFloat(o.borderTopWidth),m=r?0:parseFloat(o.borderRightWidth),w=r?0:parseFloat(o.borderBottomWidth),b=u+f,F=d+l,v=(r?0:parseFloat(o.borderLeftWidth))+m,W=g+w,y=r?0:t.offsetHeight-W-t.clientHeight,E=r?0:t.offsetWidth-v-t.clientWidth,R=b+v,z=F+W,M=r?t.width:parseFloat(o.width)-R-E,O=r?t.height:parseFloat(o.height)-z-y;if(n.has(t)){var k=n.get(t);if(k[0]===M&&k[1]===O)return}n.set(t,[M,O]);var S=Object.create(h.prototype);S.target=t,S.contentRect=new c(u,d,M,O),s.has(i)||(s.set(i,[]),p.add(i)),s.get(i).push(S)}))})),p.forEach((function(e){i.get(e).call(e,s.get(e),e)}))}return s.prototype.observe=function(i){if(i instanceof window.Element){r.has(i)||(r.set(i,new Set),o.add(i),a.set(i,window.getComputedStyle(i)));var n=r.get(i);n.has(this)||n.add(this),cancelAnimationFrame(t),t=requestAnimationFrame(d)}},s.prototype.unobserve=function(i){if(i instanceof window.Element&&r.has(i)){var n=r.get(i);n.has(this)&&(n.delete(this),n.size||(r.delete(i),o.delete(i))),n.size||r.delete(i),o.size||cancelAnimationFrame(t)}},A.DOMRectReadOnly=c,A.ResizeObserver=s,A.ResizeObserverEntry=h,A}; // eslint-disable-line\n",
       "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Left button pans, Right button zooms\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"jpeg\", \"pgf\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n",
       "\n",
       "mpl.default_extension = \"png\";/* global mpl */\n",
       "\n",
       "var comm_websocket_adapter = function (comm) {\n",
       "    // Create a \"websocket\"-like object which calls the given IPython comm\n",
       "    // object with the appropriate methods. Currently this is a non binary\n",
       "    // socket, so there is still some room for performance tuning.\n",
       "    var ws = {};\n",
       "\n",
       "    ws.binaryType = comm.kernel.ws.binaryType;\n",
       "    ws.readyState = comm.kernel.ws.readyState;\n",
       "    function updateReadyState(_event) {\n",
       "        if (comm.kernel.ws) {\n",
       "            ws.readyState = comm.kernel.ws.readyState;\n",
       "        } else {\n",
       "            ws.readyState = 3; // Closed state.\n",
       "        }\n",
       "    }\n",
       "    comm.kernel.ws.addEventListener('open', updateReadyState);\n",
       "    comm.kernel.ws.addEventListener('close', updateReadyState);\n",
       "    comm.kernel.ws.addEventListener('error', updateReadyState);\n",
       "\n",
       "    ws.close = function () {\n",
       "        comm.close();\n",
       "    };\n",
       "    ws.send = function (m) {\n",
       "        //console.log('sending', m);\n",
       "        comm.send(m);\n",
       "    };\n",
       "    // Register the callback with on_msg.\n",
       "    comm.on_msg(function (msg) {\n",
       "        //console.log('receiving', msg['content']['data'], msg);\n",
       "        var data = msg['content']['data'];\n",
       "        if (data['blob'] !== undefined) {\n",
       "            data = {\n",
       "                data: new Blob(msg['buffers'], { type: data['blob'] }),\n",
       "            };\n",
       "        }\n",
       "        // Pass the mpl event to the overridden (by mpl) onmessage function.\n",
       "        ws.onmessage(data);\n",
       "    });\n",
       "    return ws;\n",
       "};\n",
       "\n",
       "mpl.mpl_figure_comm = function (comm, msg) {\n",
       "    // This is the function which gets called when the mpl process\n",
       "    // starts-up an IPython Comm through the \"matplotlib\" channel.\n",
       "\n",
       "    var id = msg.content.data.id;\n",
       "    // Get hold of the div created by the display call when the Comm\n",
       "    // socket was opened in Python.\n",
       "    var element = document.getElementById(id);\n",
       "    var ws_proxy = comm_websocket_adapter(comm);\n",
       "\n",
       "    function ondownload(figure, _format) {\n",
       "        window.open(figure.canvas.toDataURL());\n",
       "    }\n",
       "\n",
       "    var fig = new mpl.figure(id, ws_proxy, ondownload, element);\n",
       "\n",
       "    // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n",
       "    // web socket which is closed, not our websocket->open comm proxy.\n",
       "    ws_proxy.onopen();\n",
       "\n",
       "    fig.parent_element = element;\n",
       "    fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n",
       "    if (!fig.cell_info) {\n",
       "        console.error('Failed to find cell for figure', id, fig);\n",
       "        return;\n",
       "    }\n",
       "    fig.cell_info[0].output_area.element.on(\n",
       "        'cleared',\n",
       "        { fig: fig },\n",
       "        fig._remove_fig_handler\n",
       "    );\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_close = function (fig, msg) {\n",
       "    var width = fig.canvas.width / fig.ratio;\n",
       "    fig.cell_info[0].output_area.element.off(\n",
       "        'cleared',\n",
       "        fig._remove_fig_handler\n",
       "    );\n",
       "    fig.resizeObserverInstance.unobserve(fig.canvas_div);\n",
       "\n",
       "    // Update the output cell to use the data from the current canvas.\n",
       "    fig.push_to_output();\n",
       "    var dataURL = fig.canvas.toDataURL();\n",
       "    // Re-enable the keyboard manager in IPython - without this line, in FF,\n",
       "    // the notebook keyboard shortcuts fail.\n",
       "    IPython.keyboard_manager.enable();\n",
       "    fig.parent_element.innerHTML =\n",
       "        '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "    fig.close_ws(fig, msg);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.close_ws = function (fig, msg) {\n",
       "    fig.send_message('closing', msg);\n",
       "    // fig.ws.close()\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.push_to_output = function (_remove_interactive) {\n",
       "    // Turn the data on the canvas into data in the output cell.\n",
       "    var width = this.canvas.width / this.ratio;\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] =\n",
       "        '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function () {\n",
       "    // Tell IPython that the notebook contents must change.\n",
       "    IPython.notebook.set_dirty(true);\n",
       "    this.send_message('ack', {});\n",
       "    var fig = this;\n",
       "    // Wait a second, then push the new image to the DOM so\n",
       "    // that it is saved nicely (might be nice to debounce this).\n",
       "    setTimeout(function () {\n",
       "        fig.push_to_output();\n",
       "    }, 1000);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var toolbar = document.createElement('div');\n",
       "    toolbar.classList = 'btn-toolbar';\n",
       "    this.root.appendChild(toolbar);\n",
       "\n",
       "    function on_click_closure(name) {\n",
       "        return function (_event) {\n",
       "            return fig.toolbar_button_onclick(name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    function on_mouseover_closure(tooltip) {\n",
       "        return function (event) {\n",
       "            if (!event.currentTarget.disabled) {\n",
       "                return fig.toolbar_button_onmouseover(tooltip);\n",
       "            }\n",
       "        };\n",
       "    }\n",
       "\n",
       "    fig.buttons = {};\n",
       "    var buttonGroup = document.createElement('div');\n",
       "    buttonGroup.classList = 'btn-group';\n",
       "    var button;\n",
       "    for (var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            /* Instead of a spacer, we start a new button group. */\n",
       "            if (buttonGroup.hasChildNodes()) {\n",
       "                toolbar.appendChild(buttonGroup);\n",
       "            }\n",
       "            buttonGroup = document.createElement('div');\n",
       "            buttonGroup.classList = 'btn-group';\n",
       "            continue;\n",
       "        }\n",
       "\n",
       "        button = fig.buttons[name] = document.createElement('button');\n",
       "        button.classList = 'btn btn-default';\n",
       "        button.href = '#';\n",
       "        button.title = name;\n",
       "        button.innerHTML = '<i class=\"fa ' + image + ' fa-lg\"></i>';\n",
       "        button.addEventListener('click', on_click_closure(method_name));\n",
       "        button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n",
       "        buttonGroup.appendChild(button);\n",
       "    }\n",
       "\n",
       "    if (buttonGroup.hasChildNodes()) {\n",
       "        toolbar.appendChild(buttonGroup);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = document.createElement('span');\n",
       "    status_bar.classList = 'mpl-message pull-right';\n",
       "    toolbar.appendChild(status_bar);\n",
       "    this.message = status_bar;\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = document.createElement('div');\n",
       "    buttongrp.classList = 'btn-group inline pull-right';\n",
       "    button = document.createElement('button');\n",
       "    button.classList = 'btn btn-mini btn-primary';\n",
       "    button.href = '#';\n",
       "    button.title = 'Stop Interaction';\n",
       "    button.innerHTML = '<i class=\"fa fa-power-off icon-remove icon-large\"></i>';\n",
       "    button.addEventListener('click', function (_evt) {\n",
       "        fig.handle_close(fig, {});\n",
       "    });\n",
       "    button.addEventListener(\n",
       "        'mouseover',\n",
       "        on_mouseover_closure('Stop Interaction')\n",
       "    );\n",
       "    buttongrp.appendChild(button);\n",
       "    var titlebar = this.root.querySelector('.ui-dialog-titlebar');\n",
       "    titlebar.insertBefore(buttongrp, titlebar.firstChild);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._remove_fig_handler = function (event) {\n",
       "    var fig = event.data.fig;\n",
       "    if (event.target !== this) {\n",
       "        // Ignore bubbled events from children.\n",
       "        return;\n",
       "    }\n",
       "    fig.close_ws(fig, {});\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function (el) {\n",
       "    el.style.boxSizing = 'content-box'; // override notebook setting of border-box.\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function (el) {\n",
       "    // this is important to make the div 'focusable\n",
       "    el.setAttribute('tabindex', 0);\n",
       "    // reach out to IPython and tell the keyboard manager to turn it's self\n",
       "    // off when our div gets focus\n",
       "\n",
       "    // location in version 3\n",
       "    if (IPython.notebook.keyboard_manager) {\n",
       "        IPython.notebook.keyboard_manager.register_events(el);\n",
       "    } else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function (event, _name) {\n",
       "    var manager = IPython.notebook.keyboard_manager;\n",
       "    if (!manager) {\n",
       "        manager = IPython.keyboard_manager;\n",
       "    }\n",
       "\n",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which === 13) {\n",
       "        this.canvas_div.blur();\n",
       "        // select the cell after this one\n",
       "        var index = IPython.notebook.find_cell_index(this.cell_info[0]);\n",
       "        IPython.notebook.select(index + 1);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_save = function (fig, _msg) {\n",
       "    fig.ondownload(fig, null);\n",
       "};\n",
       "\n",
       "mpl.find_output_cell = function (html_output) {\n",
       "    // Return the cell and output element which can be found *uniquely* in the notebook.\n",
       "    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n",
       "    // IPython event is triggered only after the cells have been serialised, which for\n",
       "    // our purposes (turning an active figure into a static one), is too late.\n",
       "    var cells = IPython.notebook.get_cells();\n",
       "    var ncells = cells.length;\n",
       "    for (var i = 0; i < ncells; i++) {\n",
       "        var cell = cells[i];\n",
       "        if (cell.cell_type === 'code') {\n",
       "            for (var j = 0; j < cell.output_area.outputs.length; j++) {\n",
       "                var data = cell.output_area.outputs[j];\n",
       "                if (data.data) {\n",
       "                    // IPython >= 3 moved mimebundle to data attribute of output\n",
       "                    data = data.data;\n",
       "                }\n",
       "                if (data['text/html'] === html_output) {\n",
       "                    return [cell, data, j];\n",
       "                }\n",
       "            }\n",
       "        }\n",
       "    }\n",
       "};\n",
       "\n",
       "// Register the function which deals with the matplotlib target/channel.\n",
       "// The kernel may be null if the page has been refreshed.\n",
       "if (IPython.notebook.kernel !== null) {\n",
       "    IPython.notebook.kernel.comm_manager.register_target(\n",
       "        'matplotlib',\n",
       "        mpl.mpl_figure_comm\n",
       "    );\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAAAXNSR0IArs4c6QAAIABJREFUeF7t3Qm8TeX6wPEHmQsZ/yWHIrdSqAxJA7d7U1IpFaGEBmSMkAhFkpJM0cS/JFPh1r2VRoWupKQSJRkyRTiRkOH/ed579/4fxznOPvvZe6+91/q9n49Pg/2u4buec55nr/W+78pz5MiRI0JDAAEEEEAAAQQQCIxAHgrAwFxrThQBBBBAAAEEEHACFIAEAgIIIIAAAgggEDABCsCAXXBOFwEEEEAAAQQQoAAkBhBAAAEEEEAAgYAJUAAG7IJzuggggAACCCCAAAUgMYAAAggggAACCARMgAIwYBec00UAAQQQQAABBCgAiQEEEEAAAQQQQCBgAhSAAbvgnC4CCCCAAAIIIEABSAwggAACCCCAAAIBE6AADNgF53QRQAABBBBAAAEKQGIAAQQQQAABBBAImAAFYMAuOKeLAAIIIIAAAghQABIDCCCAAAIIIIBAwAQoAAN2wTldBBBAAAEEEECAApAYQAABBBBAAAEEAiZAARiwC87pIoAAAggggAACFIDEAAIIIIAAAgggEDABCsCAXXBOFwEEEEAAAQQQoAAkBhBAAAEEEEAAgYAJUAAG7IJzuggggAACCCCAAAUgMYAAAggggAACCARMgAIwYBec00UAAQQQQAABBCgAiQEEEEAAAQQQQCBgAhSAAbvgnC4CCCCAAAIIIEABSAwggAACCCCAAAIBE6AADNgF53QRQAABBBBAAAEKQGIAAQQQQAABBBAImAAFYMAuOKeLAAIIIIAAAghQABIDCCCAAAIIIIBAwAQoAAN2wTldBBBAAAEEEECAApAYQAABBBBAAAEEAiZAARiwC87pIoAAAggggAACFIDEAAIIIIAAAgggEDABCsCAXXBOFwEEEEAAAQQQoAAkBhBAwHOBRYsWybx586R79+5SokSJhB/Prl27pHfv3jJ79mzZu3ev1KlTR5588km54IILcjyWzz77TCZPniyLFy+W5cuXy8GDB+XIkSPZ9nvhhRfkiSeekJ9++kkqVKggXbt2lS5duhzz+Y0bN0qPHj2cy+HDh6Vhw4by1FNPyRlnnJHjMfEBBBBAICcBCsCchPh7BBCIu4AWRPfff78riipVqhT3/WXcgRZXl156qXz11VfuGEqXLi3jx4+XDRs2yNKlS+XMM8887vEMGjRIHn30Ualevbrs3r1bvv/++2wLwIkTJ0qHDh2kWbNm0qhRI/nkk0/k5Zdflscee0z69OkT3s+ePXtc8Zmeni49e/aU/Pnzu+JPC8tly5ZJqVKlEmrEzhBAwH8CFID+u6acEQIpJ+BlAThjxgxp3ry5zJw5U2666SZnt23bNqlatapcffXVMnXq1ON6bt26VYoVKyaFCxeWzp07y7hx47IsAP/44w93x++iiy6SN998M7zN1q1by5w5c1zBefLJJ7v///jjj7uCUO8u1q5d2/2/lStXyrnnnuvuVGrBSUMAAQQsAhSAFj36IoCAWUDvoA0ePPiY7STqbuAtt9wiH3/8sWzatEny5s0bPo577rlHpkyZIjt27JCCBQtGdJ7HKwD/9a9/yTXXXCP//Oc/pXHjxuHtffrpp3LxxRe7O4FaDGrTR9DatADM2PSu4Y8//iirV6+O6Hj4EAIIIJCdAAUgsYEAAp4K6Lg5fQT66quvusec+ghW2w033CBFixbN8th0nJ7+yanly5cvfFctu8/qI179owVaxqZj9e688043ru+8887LaVfu749XAA4dOlT69+8vesewbNmy4e0dOHDA3T3U8Y867lAfSRcpUkTatWvnHkVnbAMGDJAhQ4bIb7/9JieddFJEx8SHEEAAgawEKACJCwQQ8Fwgt4+As7trmPlEKlasKGvXrj3u+Z144onuEbAWfBlb6I7d22+/7cbrRdKOVwDq302YMMFNEsnctCC84oorXBG8fft2KVOmjDz88MOiBV/GpgXhvffe6x4H/+Uvf4nkkPgMAgggkKUABSCBgQACngvktgBcs2aN6J+cmt5Zq1+//nE/pncJ9XFv5rttH3zwgSvKdGZw06ZNc9qV+/vjFYDt27d3BV5Wdy7T0tLcpI/QWED97+HDh7vxfhnbiy++KLqdL7/8UmrWrBnRMfEhBBBAICsBCkDiAgEEPBfIbQEYywPmDmAsNdkWAgikigAFYKpcKY4TAR8L5LYA1GVS9E9OTe/u6ePU4zXGAOakyN8jgIAfBSgA/XhVOScEUkxAJz/06tUr4nUAYzkG8Oabb3br8WWeBXz33XfLK6+8ErNZwDr7t0mTJsfMAtZFsPUx9UsvvSS33Xabu3K69EuePHmOmQV85ZVXulnA+oeGAAIIWAQoAC169EUAgZgI6OSIjh07Rjy2LZZjAKdPny4tWrQ4ah1AnYihdwZ18se0adPC5xgqvCpXrpzleee0DuBpp53mlnx54403wv216Hv99dfdOoAlS5Z0/1/H//Xt21eWLFkitWrVcv9v1apVUq1aNVco66xpGgIIIGARoAC06NEXAQRiIqCFjq59p+vjaTGmb7649tprs10GJiY7/e9GDh06JJdccol88803R70JZP369a4AyzjbNvSWkowzi9etW+fW8NOmCzzrK+EeeeQR9986Czl0V0//OzSLVxecDr0JRO/86RIx/fr1C5+WvlHk/PPPd28W0YJPPUaOHCl6rPomkJwea8fSh20hgIA/BSgA/XldOSsEUk5A17fTO4GbN292a+ElaiFohdq5c6cr/nQWrr6xQx/B6rjE0N23EGZWBeBHH33k3tObVbv88stF/z5je+6559x6f6F3Aetdw27durlHvhnbzz//fNS7gBs0aODWSaxSpUrKXVsOGAEEkk+AAjD5rglHhAACCCCAAAIIxFWAAjCuvGwcAQQQQAABBBBIPgEKwOS7JhwRAggggAACCCAQVwEKwLjysnEEEEAAAQQQQCD5BHxTAH788ccyYsQIWbp0qRtEHsnrm3Rw9n333SfffvutVKhQwb2o/Y477kiKq7R//34ZNmyYPPDAA1KwYMGkOCY/HgTOibmqOOOcGIHE7IV4xjkxAvHdi28KwLfeeksWLlwoF154odx44405FoA6A+/cc8+VDh06yJ133invv/++dO/e3S3SGumL3+N5aX777TcpXry4pKenS7FixeK5q0BvG+fEXH6ccU6MQGL2QjzjnBiB+O7FNwVgRiZdTiGnO4B9+vRxxZ6u/RVquv7Yrl275O23346vegRb5xdMBEgx+AjOMUCMYBM4R4AUg4/gHAPECDaBcwRIMfgIzjFAPM4mAlsAXnbZZXLBBRfIqFGjwjyTJk1ydwH1rpvXjcBPzBXAGefECCRmL8QzzokRSMxeiOf4Oge2AKxataq0bdvWjbELtX/9619yzTXXyN69e6Vw4cLHyOu4D/0TarpY7Y4dO6RUqVLHLOJqvWwa+DouUV8PxSNgq2b2/XGOn23GLeOMc2IEErMX4jn1nY8cOSK//vqre/1i3rx5wyekY+6DMu6eAjAXBWCkL6BPzI8Ge0EAAQQQQACBWAoMHDhQNNcHoQW2AIzmEXDmO4D6qDgtLY27dEH4SeEcEUAAAQR8IxC6i6vv/NYJl6HGHcAUv8SRTgLRR75ff/11+GxbtmzpHulGOgmE8QkpHigcPgIIIIBAIAXI3yK+uQO4Z88eWb16tQvk888/X0aOHOle0K7P9/UunY7127hxo7z00kvuM6FlYO69915p166dfPDBB9K1a9dcLQNDAAXy9wYnjQACCCCQ4gLkbx8VgLqosxZ8mVubNm1k8uTJboHntWvXin4u1PTfe/ToIStWrJDTTjtNBgwYkKuFoAmgFP8NwOEjgAACCARSgPztowLQiwgmgLxQZ58IIIAAAgjYBMjfFICmCCKATHx0RgABBBBAwBMB8jcFoCnwCCATH50RQAABBBDwRID8TQFoCjwCyMRHZwQQQAABBDwRIH9TAJoCjwAy8dEZAQQQQAABTwTI3xSApsAjgEx8dEYAAQQQQMATAfI3BaAp8AggEx+dEUAAAQQQ8ESA/E0BaAo8AsjER2cEEEAAAQQ8ESB/UwCaAo8AMvHRGQEEEEAAAU8EyN8UgKbAI4BMfHRGAAEEEEDAEwHyNwWgKfAIIBMfnRFAAAEEEPBEgPxNAWgKPALIxEdnBBBAAAEEPBEgf1MAmgKPADLx0RkBBBBAAAFPBMjfFICmwCOATHx0RgABBBBAwBMB8jcFoCnwCCATH50RQAABBBDwRID8TQFoCjwCyMRHZwQQQAABBDwRIH9TAJoCjwAy8dEZAQQQQAABTwTI3xSApsAjgEx8dEYAAQQQQMATAfI3BaAp8AggEx+dEUAAAQQQ8ESA/E0BaAo8AsjER2cEEEAAAQQ8ESB/UwCaAo8AMvHRGQEEEEAAAU8EyN8UgKbAI4BMfHRGAAEEEEDAEwHyNwWgKfAIIBMfnRFAAAEEEPBEgPxNAWgKPALIxEdnBBBAAAEEPBEgf1MAmgKPADLx0RkBBBBAAAFPBMjfFICmwCOATHx0RgABBBBAwBMB8jcFoCnwCCATH50RQAABBBDwRID8TQFoCjwCyMRHZwQQQAABBDwRIH9TAJoCjwAy8dEZAQQQQAABTwTI3xSApsAjgEx8dEYAAQQQQMATAfI3BaAp8AggEx+dEUAAAQQQ8ESA/E0BaAo8AsjER2cEEEAAAQQ8ESB/UwCaAo8AMvHRGQEEEEAAAU8EyN8UgKbAI4BMfHRGAAEEEEDAEwHyNwWgKfAIIBMfnRFAAAEEEPBEgPxNAWgKPALIxEdnBBBAAAEEPBEgf1MAmgKPADLx0RkBBBBAAAFPBMjfFICmwCOATHx0RgABBBBAwBMB8jcFoCnwCCATH50RQAABBBDwRID8TQFoCjwCyMRHZwQQQAABBDwRIH9TAJoCjwAy8dEZAQQQQAABTwTI3xSApsAjgEx8dEYAAQQQQMATAfI3BaAp8AggEx+dEUAAAQQQ8ESA/E0BaAo8AsjER2cEEEAAAQQ8ESB/UwCaAo8AMvHRGQEEEEAAAU8EyN8UgKbAI4BMfHRGAAEEEEDAEwHyNwWgKfAIIBMfnRFAAAEEEPBEgPxNAWgKPALIxEdnBBBAAAEEPBEgf1MAmgKPADLx0RkBBBBAAAFPBMjfFICmwCOATHx0RgABBBBAwBMB8jcFoCnwCCATH50RQAABBBDwRID8TQFoCjwCyMRHZwQQQAABBDwRIH9TAJoCjwAy8dEZAQQQQAABTwTI3xSApsAjgEx8dEYAAQQQQMATAfI3BaAp8AggEx+dEUAAAQQQ8ESA/E0BaAo8AsjER2cEEEAAAQQ8ESB/UwCaAo8AMvHRGQEEEEAAAU8EyN8UgKbAI4BMfHRGAAEEEEDAEwHyNwWgKfAIIBMfnRFAAAEEEPBEgPxNAWgKPALIxEdnBBBAAAEEPBEgf1MAmgKPADLx0RkBBBBAAAFPBMjfFICmwCOATHx0RgABBBBAwBMB8jcFoCnwCCATH50RQAABBBDwRID8TQFoCjwCyMRHZwQQQAABBDwRIH/7rAAcN26cjBgxQrZs2SI1atSQMWPGSJ06dbINrlGjRskzzzwj69evl9KlS8tNN90kw4YNk0KFCkUUkARQREx8CAEEEEAAgaQSIH/7qACcPn263H777TJhwgSpW7euaHE3c+ZMWbVqlZQtW/aYwJs6daq0a9dOXnzxRbn44ovl+++/lzvuuENatGghI0eOjChQCaCImPgQAggggAACSSVA/vZRAahFX+3atWXs2LEuyA4fPiwVKlSQLl26SN++fY8JvM6dO8t3330n77//fvjvevbsKYsXL5YFCxZEFKgEUERMfAgBBBBAAIGkEiB/+6QAPHDggBQpUkRmzZolTZs2DQdZmzZtZNeuXTJ37tws7wB26tRJ5s2b5x4Tr1mzRq655hq57bbbpF+/fhEFKgEUERMfQgABBBBAIKkEyN8+KQA3bdok5cuXl0WLFkm9evXCQda7d2+ZP3++u6uXVRs9erT06tVLjhw5IgcPHpQOHTq4MYHZtf3794v+CTUNIL3LmJ6eLsWKFUuq4OZgEEAAAQQQQCBrgVABuGHDhqPyd8GCBUX/BKHlOaLVT4q3aArAjz76yI33GzJkiBszuHr1aunWrZvcddddMmDAgCxFBg0aJIMHDz7m7ygAUzyAOHwEEEAAgUAJhArAzCc9cOBA0VwfhOaLAjCaR8CXXnqpXHTRRW7WcKhNmTJF7r77btmzZ4/kzZv3mOvPHcAg/EhwjggggAACfhfgDqBPHgFroOpdPB3Lp0u/aNNJIGlpaaKTPbKaBHLhhRfK3/72Nxk+fHg4zl999VVp37697N69W/Lly5dj/DOGIEciPoAAAggggEDSCZC/fVQA6jIwOulj4sSJrhDUZWBmzJghK1eulHLlyrklYnScoK7zp01v8epyL88++2z4EXDHjh1FC0PdViSNAIpEic8ggAACCCCQXALkbx8VgBpaugRMaCHomjVrik7y0DuD2ho0aCCVKlWSyZMnu//WSR9Dhw6Vl19+WTZu3ChlypSRa6+91v2/EiVKRBSpBFBETHwIAQQQQACBpBIgf/usAEx0dBFAiRZnfwgggAACCNgFyN8UgKYoIoBMfHRGAAEEEEDAEwHyNwWgKfAIIBMfnRFAAAEEEPBEgPxNAWgKPALIxEdnBBBAAAEEPBEgf1MAmgKPADLx0RkBBBBAAAFPBMjfFICmwCOATHx0RgABBBBAwBMB8jcFoCnwCCATH50RQAABBBDwRID8TQFoCjwCyMRHZwQQQAABBDwRIH9TAJoCjwAy8dEZAQQQQAABTwTI3xSApsAjgEx8dEYAAQQQQMATAfI3BaAp8AggEx+dEUAAAQQQ8ESA/E0BaAo8AsjER2cEEEAAAQQ8ESB/UwCaAo8AMvHRGQEEEEAAAU8EyN8UgKbAI4BMfHRGAAEEEEDAEwHyNwWgKfAIIBMfnRFAAAEEEPBEgPxNAWgKPALIxEdnBBBAAAEEPBEgf1MAmgKPADLx0RkBBBBAAAFPBMjfFICmwCOATHx0RgABBBBAwBMB8jcFoCnwCCATH50RQAABBBDwRID8TQFoCjwCyMRHZwQQQAABBDwRIH9TAJoCjwAy8dEZAQQQQAABTwTI3xSApsAjgEx8dEYAAQQQQMATAfI3BaAp8AggEx+dEUAAAQQQ8ESA/E0BaAo8AsjER2cEEEAAAQQ8ESB/UwCaAo8AMvHRGQEEEEAAAU8EyN8UgKbAI4BMfHRGAAEEEEDAEwHyNwWgKfAIIBMfnRFAAAEEEPBEgPxNAWgKPALIxEdnBBBAAAEEPBEgf1MAmgKPADLx0RkBBBBAAAFPBMjfFICmwCOATHx0RgABBBBAwBMB8jcFoCnwCCATH50RQAABBBDwRID8TQFoCjwCyMRHZwQQQAABBDwRIH9TAJoCjwAy8dEZAQQQQAABTwTI3xSApsAjgEx8dEYAAQQQQMATAfI3BaAp8AggEx+dEUAAAQQQ8ESA/E0BaAo8AsjER2cEEEAAAQQ8ESB/UwCaAo8AMvHRGQEEEEAAAU8EyN8UgKbAI4BMfHRGAAEEEEDAEwHyNwWgKfAIIBMfnRFAAAEEEPBEgPxNAWgKPALIxEdnBBBAAAEEPBEgf1MAmgKPADLx0RkBBBBAAAFPBMjfFICmwCOATHx0RgABBBBAwBMB8jcFoCnwCCATH50RQAABBBDwRID8TQFoCjwCyMRHZwQQQAABBDwRIH9TAJoCjwAy8dEZAQQQQAABTwTI3xSApsAjgEx8dEYAAQQQQMATAfI3BaAp8AggEx+dEUAAAQQQ8ESA/E0BaAo8AsjER2cEEEAAAQQ8ESB/UwCaAo8AMvHRGQEEEEAAAU8EyN8UgKbAI4BMfHRGAAEEEEDAEwHyNwWgKfAIIBMfnRFAAAEEEPBEgPxNAWgKPALIxEdnBBBAAAEEPBEgf1MAmgKPADLx0RkBBBBAAAFPBMjfFICmwCOATHx0RgABBBBAwBMB8jcFoCnwCCATH50RQAABBBDwRID8TQFoCjwCyMRHZwQQQAABBDwRIH9TAJoCjwAy8dEZAQQQQAABTwTI3xSApsAjgEx8dEYAAQQQQMATAfI3BaAp8AggEx+dEUAAAQQQ8ESA/E0BaAo8AsjER2cEEEAAAQQ8ESB/UwCaAo8AMvHRGQEEEEAAAU8EyN8UgKbAI4BMfHRGAAEEEEDAEwHyNwWgKfAIIBMfnRFAAAEEEPBEgPxNAWgKPALIxEdnBBBAAAEEPBEgf/usABw3bpyMGDFCtmzZIjVq1JAxY8ZInTp1sg2uXbt2yYMPPiivv/667NixQypWrCijRo2Sxo0bRxSQBFBETHwIAQQQQACBpBIgf/uoAJw+fbrcfvvtMmHCBKlbt64r5GbOnCmrVq2SsmXLHhN4Bw4ckPr167u/69evn5QvX17WrVsnJUqUcMVjJI0AikSJzyCAAAIIIJBcAuRvHxWAWvTVrl1bxo4d66Ls8OHDUqFCBenSpYv07dv3mMjTQlHvFq5cuVLy588fVWQSQFGx0QkBBBBAAAFPBcjfPikA9W5ekSJFZNasWdK0adNwULVp00b0Me/cuXOPCTR9zFuyZEnXT/++TJky0rJlS+nTp4/ky5cvy8Dcv3+/6J9Q0wDSIjM9PV2KFSvmaTCzcwQQQAABBBCITCBUAG7YsOGo/F2wYEHRP0FoeY4cOXIk1U9006ZN7hHuokWLpF69euHT6d27t8yfP18WL158zCmeddZZsnbtWmnVqpV06tRJVq9e7f7ZtWtXGThwYJYkgwYNksGDBx/zdxSAqR5BHD8CCCCAQJAEQgVg5nPW/K+5PggtsAVg1apVZd++ffLTTz+F7/iNHDnSPRbevHkzdwCDEP2cIwIIIIBAIAW4AxjgR8CXX365G/v33nvvhYP/rbfecjOA9TFvgQIFcvyhYAxBjkR8AAEEEEAAgaQTIH/7pADUyNJJILrkiy79ok0ngaSlpUnnzp2znASiM3+nTp0qa9askbx587o+Tz/9tAwfPlz0kXIkjQCKRInPIIAAAgggkFwC5G8fFYC6DIxO+pg4caIrBHUZmBkzZrhZvuXKlXNLxOg4wWHDhrko1IGf1apVc310pvAPP/wg7dq1c2MAdW3ASBoBFIkSn0EAAQQQQCC5BMjfPioANbR0CZjQQtA1a9aU0aNHuzuD2ho0aCCVKlWSyZMnh6Pw008/lR49esiyZctccdi+ffvjzgLOHL4EUHL9QHM0CCCAAAIIRCJA/vZZARjJRY/lZwigWGqyLQQQQAABBBIjQP6mADRFGgFk4qMzAggggAACngiQvykATYFHAJn46IwAAggggIAnAuRvCkBT4BFAJj46I4AAAggg4IkA+ZsC0BR4BJCJj84IIIAAAgh4IkD+pgA0BR4BZOKjMwIIIIAAAp4IkL8pAE2BRwCZ+OiMAAIIIICAJwLkbwpAU+ARQCY+OiOAAAIIIOCJAPmbAtAUeASQiY/OCCCAAAIIeCJA/qYANAUeAWTiozMCCCCAAAKeCJC/KQBNgUcAmfjojAACCCCAgCcC5G8KQFPgEUAmPjojgAACCCDgiQD5mwLQFHgEkImPzggggAACCHgiQP6mADQFHgFk4qMzAggggAACngiQvykATYFHAJn46IwAAggggIAnAuRvCkBT4BFAJj46I4AAAggg4IkA+ZsC0BR4BJCJj84IIIAAAgh4IkD+pgA0BR4BZOKjMwIIIIAAAp4IkL8pAE2BRwCZ+OiMAAIIIICAJwLkbwpAU+ARQCY+OiOAAAIIIOCJAPmbAtAUeASQiY/OCCCAAAIIeCJA/qYANAUeAWTiozMCCCCAAAKeCJC/KQBNgUcAmfjojAACCCCAgCcC5G8KQFPgEUAmPjojgAACCCDgiQD5mwLQFHgEkImPzggggAACCHgiQP6mADQFHgFk4qMzAggggAACngiQvykATYFHAJn46IwAAggggIAnAuRvCkBT4BFAJj46I4AAAggg4IkA+ZsC0BR4BJCJj84IIIAAAgh4IkD+pgA0BR4BZOKjMwIIIIAAAp4IkL8pAE2BRwCZ+OiMAAIIIICAJwLkbwpAU+ARQCY+OiOAAAIIIOCJAPmbAtAUeASQiY/OCCCAAAIIeCJA/qYANAUeAWTiozMCCCCAAAKeCJC/KQBNgUcAmfjojAACCCCAgCcC5G8KQFPgEUAmPjojgAACCCDgiQD5mwLQFHgEkImPzggggAACCHgiQP6mADQFHgFk4qMzAggggAACngiQvykATYFHAJn46IwAAggggIAnAuRvCkBT4BFAJj46I4AAAggg4IkA+ZsC0BR4BJCJj84IIIAAAgh4IkD+pgA0BR4BZOKjMwIIIIAAAp4IkL8pAE2BRwCZ+OiMAAIIIICAJwLkbwpAU+ARQCY+OiOAAAIIIOCJAPmbAtAUeASQiY/OCCCAAAIIeCJA/qYANAUeAWTiozMCCCCAAAKeCJC/KQBNgUcAmfjojAACCCCAgCcC5G8KQFPgEUAmPjojgAACCCDgiQD5mwLQFHgEkImPzggggAACCHgiQP6mADQFHgFk4qMzAggggAACngiQvykATYFHAJn46IwAAggggIAnAuRvCkBT4BFAJj46I4AAAggg4IkA+ZsC0BR4BJCJj84IIIAAAgh4IkD+pgA0BR4BZOKjMwIIIIAAAp4IkL8pAE2BRwCZ+OiMAAIIIICAJwLkbwpAU+ARQCY+OiOAAAIIIOCJAPmbAtAUeASQiY/OCCCAAAIIeCJA/qYANAUeAWTiozMCCCCAAAKeCJC/KQBNgUcAmfjojAACCCCAgCcC5G8KQFPgEUAmPjojgAACCCDgiQD5mwLQFHgEkImPzggggAACCHgiQP6mADQFHgFk4qMzAggggAACngiQvykATYFHAJn46IwAAggggIAnAuRvCkCLpaOxAAAgAElEQVRT4BFAJj46I4AAAggg4IkA+dtnBeC4ceNkxIgRsmXLFqlRo4aMGTNG6tSpk2NwTZs2TW699Va5/vrrZc6cOTl+PvQBAihiKj6IAAIIIIBA0giQv31UAE6fPl1uv/12mTBhgtStW1dGjRolM2fOlFWrVknZsmWzDbq1a9fKJZdcImeccYaULFmSAjBpfjw5EAQQQAABBOIjQAHoowJQi77atWvL2LFjXbQcPnxYKlSoIF26dJG+fftmGUGHDh2Syy67TNq1ayeffPKJ7Nq1iwIwPj9rbBUBBBBAAIGkEaAA9EkBeODAASlSpIjMmjVLmjZtGg6wNm3auKJu7ty5WQbdwIEDZfny5TJ79my54447ciwA9+/fL/on4yNgLTLT09OlWLFiSRPYHAgCCCCAAAIIZC8QKgA3bNhwVP4uWLCg6J8gtDxHjhw5kuonumnTJilfvrwsWrRI6tWrFz6d3r17y/z582Xx4sXHnOKCBQukRYsWsmzZMildunREBeCgQYNk8ODBx2yLAjDVI4jjRwABBBAIkkCoAMx8znpjSHN9EFogC8Ddu3dL9erVZfz48XL11Ve768wdwCCEO+eIAAIIIICACHcAA/oIWO/6nX/++ZIvX77wz4GOGdSWN29eN3GkcuXKOf6MMIYgRyI+gAACCCCAQNIJkL99UgBqZOkkEF3yRZd+0aYFXVpamnTu3PmYSSD79u2T1atXHxWQ/fv3F70z+PTTT0vVqlWlQIECOQYsAZQjER9AAAEEEEAg6QTI3z4qAHUZGJ30MXHiRFcI6jIwM2bMkJUrV0q5cuXcEjE6TnDYsGFZBmIkj4AzdySAku5nmgNCAAEEEEAgRwHyt48KQL3augRMaCHomjVryujRo92dQW0NGjSQSpUqyeTJkykAc/zR4AMIIIAAAgj4V4AC0GcFYKJDlQBKtDj7QwABBBBAwC5A/qYANEURAWTiozMCCCCAAAKeCJC/KQBNgUcAmfjojAACCCCAgCcC5G8KQFPgEUAmPjojgAACCCDgiQD5mwLQFHgEkImPzggggAACCHgiQP6mADQFHgFk4qMzAggggAACngiQvykATYFHAJn46IwAAggggIAnAuRvCkBT4BFAJj46I4AAAggg4IkA+ZsC0BR4BJCJj84IIIAAAgh4IkD+pgA0BR4BZOKjMwIIIIAAAp4IkL8pAE2BRwCZ+OiMAAIIIICAJwLkbwpAU+ARQCY+OiOAAAIIIOCJAPmbAtAUeASQiY/OCCCAAAIIeCJA/qYANAUeAWTiozMCCCCAAAKeCJC/KQBNgUcAmfjojAACCCCAgCcC5G8KQFPgEUAmPjojgAACCCDgiQD5mwLQFHgEkImPzggggAACCHgiQP6mADQFHgFk4qMzAggggAACngiQvykATYFHAJn46IwAAggggIAnAuRvCkBT4BFAJj46I4AAAggg4IkA+ZsC0BR4BJCJj84IIIAAAgh4IkD+pgA0BR4BZOKjMwIIIIAAAp4IkL8pAE2BRwCZ+OiMAAIIIICAJwLkbwpAU+ARQCY+OiOAAAIIIOCJAPmbAtAUeASQiY/OCCCAAAIIeCJA/qYANAUeAWTiozMCCCCAAAKeCJC/KQBNgUcAmfjojAACCCCAgCcC5G8KQFPgEUAmPjojgAACCCDgiQD5mwLQFHgEkImPzggkhcCBAwfk559/lnXr1snatWuP+qceYMWKFaVSpUpH/bNChQqSP3/+pDh+DgIBBHIvQP6mAMx91GToQQCZ+OiMQEIE/vjjD1m/fv1RhV3GYm/jxo1y5MiRXB1Lnjx5pHz58kcVhRkLxbS0NClcuHCutsmHEUAgcQLkbwpAU7QRQCY+OiMQE4Hdu3e74i70J/NdvK1bt+a4n0KFCrliLmMRp/+uLfN29b/37duX4zbLlSsXvnOY1V3EE088Mcdt8AEEEIiPAPmbAtAUWQSQiY/OCOQooHfmdu3aleXj2dBdvB07duS4HS22Mj/GzViUlS1bVvSuXiRNj+mXX37J9o6iFqB79uzJcVMlS5Y87jGVKFEi4mPKcWd8AAEEjhIgf1MAmn4kCCATH50RcI9et23bdtxiSu/w5dROPvnk495t07+PtMDLaV85/b2e086dO49btOrf59SKFSuW5V3JUCFbpkyZhJ1TTsfK3yOQagLkbwpAU8wSQCY+OgdEQAuib775xv3J/HhW7+LpGL2cmt6hCxU+WT1OPemkk3LaRFL9vf7uON5ja73DmFPTMYaZH1ur0XnnnSfVqlWjOMwJkL8PtAD5mwLQ9ANAAJn46OxzgZ9++kmmTp0qr7zyinz33XfZnq3emTv11FOznG2rBY5OqChSpIjPtY4+vb179x534sqmTZuOO3HlnHPOkVatWknLli1d4UxDAIGjBcjfFICmnwkCyMRHZx8KbN++XWbOnOmKvoULF4bPsGDBglK3bl05/fTTj7lrpUuqFChQwIca8TslXbpmw4YNxzw616J78eLFsn///vDO69ev74rBW265RUqVKhW/g2LLCKSQAPmbAtAUrgSQiY/OPhHQu1X/+Mc/XNH39ttvy8GDB92Z6Z29v/71r674uPHGG6V48eI+OePkPo309HR5/fXXZcqUKfLhhx+G7xSecMIJctVVV0nr1q3l2muvDdxd1eS+ahxdogXI3xSAppgjgEx8dE5hAS3yPvjgA1dkzJ49+6hZr+eff74rMlq0aOEe7dK8E9A1DqdNm+aK8y+//DJ8IDorWotyLc61SNfikIZAkATI3xSApngngEx8dE4xAZ3M8fnnn7tiQouKjOvr6aNdHW+mBcXZZ5+dYmcWjMPVcZh67fSPTsYJNV2vUIt1vXa1atVi8kgwwiHwZ0n+pgA0/RAQQCY+OqeIwOrVq13RoBM6vv/++/BR63iy5s2bu8KhXr16FA4pcj21kP/000/dNZ0+fbr8+uuv4SOvWrWqu576p3LlyilyRhwmArkXIH9TAOY+ajL0IIBMfHROYgFdhkSLAy0SdFJBqOnSI9dff717xHvllVfyPtwkvoaRHNqff/4p77zzjrvOc+fOPWpJHp20o4WgFvm6DA8NAT8JkL8pAE3xTACZ+OicZAL69oo5c+a4YuDdd9+VQ4cOuSPMmzev/P3vf3fFQNOmTSXV1txLMuakPRxdcDvj9T98+LA71nz58h11/XmFXdJeQg4sFwLkbwrAXITLsR8lgEx8dE4CAb0DNG/evPAdIJ3RG2p16tQJ3wHScWK04Ahs2bIlfAd4yZIl4RPX9Rj1S4B+GdAvBfnz5w8OCmfqKwHyNwWgKaAJIBMfnT0S0DFg//73v8NjwHTtvlCrUqVKeAzYmWee6dERsttkEtBxn6EFvXU8aKiVLl06PAb0oosuYgxoMl00jiVHAfI3BWCOQXK8DxBAJj46J1hg5cqV4ckca9asCe9dx3eFZoHWrl2bRJ7g65Iqu9MvDno3UJf+0fGhGV9Xd8YZZ4RngZ911lmpckocZ4AFyN8UgKbwJ4BMfHROgMDmzZvdki2atL/44ovwHosWLSo33HCDm8xxxRVXsA5cAq6Fn3ah60C+99577guFrgP5+++/h0/vwgsvdHeR9UvFKaec4qfT5lx8JED+pgA0hTMBZOKjc5wENC71TRCanHWx5tBgfl3st1GjRi45X3fddaJFIA0Bq4AWfxnfBJNx8lDGN8EUK1bMuiv6IxAzAfI3BaApmAggEx+dYyig74Z96623XNH3xhtvyL59+8Jb1zX69E7fzTffLGXKlInhXtkUAkcLbNu2TWbMmOHiUNcaDLVChQq5189pHOrr6Hj3M5HjtQD5mwLQFIMEkImPzkYBvbO3cOFC93h35syZsnPnzvAWdRyW3unTt3Po+CwaAokW0HGmockjOv401EqWLOm+jGh81q9f3y0zREMg0QLkbwpAU8wRQCY+OkcpoK/0eumll1xyXb9+fXgrOt7q1ltvdYlV38ebJ0+eKPdANwRiJ6CTR/Q9xPpFRcej6rjUUEtLS3NfUm6//XZeIRg7crYUgQD5mwIwgjDJ/iMEkImPzrkU0GLvwQcfdIk01HRR5mbNmrmir2HDhm7RXhoCySqg4wM//PBD94j4tddeE118OtT08fDQoUNFi0IaAvEWIH9TAJpijAAy8dE5QoH09HQZNmyYjBo1Svbv3+966Xiq2267TZo0aSL6ejYaAqkm8Mcff8ibb74pL7/8shu3qq1gwYLSo0cP6du3rxQvXjzVTonjTSEB8jcFoClcCSATH51zENC3dEycOFEGDx4socWaGzRoIE888YToUhs0BPwi8Pnnn0uvXr1k/vz57pR0kemBAwfKPffcw9tG/HKRk+w8yN8UgKaQJIBMfHTORkDHTM2dO1d69+4tP/zwg/uUTup4/PHH3R0/xvYROn4U0LjXO4L333+/rFq1yp1i1apVZfjw4XL99dcT93686B6eE/mbAtAUfgSQiY/OWQh89tln0rNnT1mwYIH7W1225eGHH5Y777yTxZqJmEAI6J3v559/3t0B1GVltF166aXuzre+n5qGQCwEyN8UgKY4IoBMfHTOILB27Vp54IEH3CxJbbpumhaCeheQBXQJlSAK6O9Xvfs3cuTI8LqWOsv90UcflUqVKgWRhHOOoQD5mwLQFE4EkImPziJu7T5NaKNHjxZdzFkf7+qSGEOGDJHTTjsNIwQCL7BhwwYZMGCAW/pIHxPrItLdunWTfv36SYkSJQLvA0B0AuRvCsDoIue/vQggE1+gO2ux98wzz7jHuzt27HAW+k5efcxVs2bNQNtw8ghkJaBrCepEEX29oTZdUFofE3fo0IE3ixAyuRYgf1MA5jpoMnYggEx8geysdzD0Pb19+vSRH3/80Rmcc845rvDTV2QxwSOQYcFJRyigPz/6ykOdKLJixQrXq0qVKvLYY4/JjTfeyM9PhI58TIT8TQFo+jkggEx8gev873//243rW7RokTv3cuXKySOPPCJt27ZlgkfgooETtggcPHhQXnzxRXnooYdk69atblMXX3yxPPnkk3LRRRdZNk3fgAiQvykATaFOAJn4AtNZ7/TpBA99X6+2IkWKuEdZehfjxBNPDIwDJ4pArAX0TSJ693zEiBGiC0tru+WWW9zC6bwDO9ba/toe+ZsC0BTRBJCJz/eddWyfTuYYO3as6NIW+ni3Xbt2btzfqaee6vvz5wQRSJTAxo0b3d3ASZMmuYki+fPnly5durhXJ+pYQRoCmQXI3xSApp8KAsjE59vO+rq2cePGuce7u3btcufZqFEjt5Bz9erVfXvenBgCXgssX77c3VmfN2+eO5STTz7ZzSDu1KmTe80cDYGQAPmbAtD000AAmfh811nvPOhjXn2P6U8//eTO77zzznOPp7QApCGAQGIE3nnnHTfM4ptvvnE71MfBOlHkpptuYqJIYi5B0u+F/O2zAlDvumiy3bJli9SoUUPGjBmT7crxzz33nFtXKvQLQt+tquux5WaleQIo6X/GE3aACxcudAlHJ3poO+WUU9zj3zZt2ki+fPkSdhzsCAEE/iNw6NAhmTx5srsDuHnzZvf/dIKIThTRCSO0YAuQv31UAE6fPt0toDthwgSpW7eujBo1yt2N0XdKli1b9phIb9WqldSvX9/9ItC3LuiK87Nnz5Zvv/1WypcvH9FPBgEUEZOvP7R69Wp3x++1115z51m0aFH39g6d7av/TkMAAW8Ffv/9d1f06RAM/XdtzZo1c3cEdQkZWjAFyN8+KgC16Ktdu7YbcK/t8OHDUqFCBTcQWBN0Tk2/Lep4Ee2vhWQkjQCKRMmfn9m+fbsb4zd+/HjRJSny5s0r7du3l8GDB7u7fzQEEEguAb0LqAtHv/DCCy4/6EQRHRuodwhLlSqVXAfL0cRdgPztkwJQ36qgS2vMmjVLmjZtGg4cffymg/Dnzp2bYzDpcgJ6p1DvGjZp0iTLz+vgfv0TahpAWmSmp6fzvtYchf3xgX379rmhBUOHDnXXXVvjxo3d3YVq1ar54yQ5CwR8LKDDfvQuvS4ora148eLSv39/6dy5s3saRAuGQKgA1FcNZnzfuk4WCsqEoTxHdOR6irdNmza5x7a6wG69evXCZ6M/5PPnz5fFixfneIb6TVAHDusj4Ox+CQwaNMjd4cncKABz5E35D+gdAx1moOv5rVu3zp2PjjPVNcj+9re/pfz5cQIIBE3g3XffdTOGv/rqK3fqFStWdOsHNm/e3N3Rp/lbIFQAZj5LvUusuT4IjQJQxI0F0Ts4H3300XGX6eAOYBB+JI49x48//thN8FiyZIn7S/2yoXcAW7duzQSPYIYEZ+0TAR36M2XKFLdeoK4lqE2HEumYwUsvvdQnZ8lpZCXAHUAeAbs7ODpb87333pNatWrl6ieFMQS54kq5D3///ffunb1z5sxxx65v7dDxpD169HBDDmgIIOAPgb1798pTTz3lbgbs2bPHnZQOJ9LJgVWrVvXHSXIWRwmQv31SAOpV1UkguoSLjs/Spo/s0tLS3LiO7CaB6F0/vZOjj36jeX8kAeTP3yjbtm1zb+vQGeU6wUOXcbn77rvdAHJ9fy8NAQT8KaDvFdbHf7pMmN4dPOGEE6RDhw7uLSNlypTx50kH9KzI3z4qAHV8lk76mDhxoisEdRmYGTNmyMqVK13S1pm9+uhOx3ho0292+kM9depUtxxMqOldnkjfz0oA+es3h75L9Omnn3YxotdW27XXXuti5eyzz/bXyXI2CCCQrcCKFSvc3f8333zTfUYnCfTr10+6du0qhQsXRs4HAuRvHxWAGo+6hEtoIeiaNWvK6NGj3Z1BbQ0aNJBKlSq5hUG16b+HBvNnjOXcDAAlgHzwW+C/d4v1i4D+gtcZYdouuOACNzygYcOG/jhJzgIBBHIt8OGHH7rxv1988YXrq0+V9IUBt956KxNFcq2ZXB3I3z4rABMdXgRQosVjvz+d+KOLNod+weuyPvoLvmXLlvyCjz03W0Qg5QR0OFFWXxB1oojeWKClpgD5mwLQFLkEkInP08460Pu+++5zY3208YjH08vBzhFIeoHQEBH9gqjrxmrTscFaCEY6bCjpTzJAB0j+pgA0hTsBZOLzrLO+r/e2224TfY1bnjx5pGPHjm7gN4O8Pbsk7BiBlBHQSWI6VEgniekyuvo6uZdffjmqiYQpc9I+PFDyNwWgKawJIBNfwjv/+eef7vVtOvM79KrAl156icc4Cb8S7BCB1BfQ8YE68VDHDetKAbqWoL5RRF8xR0t+AfI3BaApSgkgE19CO69atcrd9Qst5tyqVSs3aahEiRIJPQ52hgAC/hHQV43qUmOvvPKKOyldRFoXlmbtwOS/xuRvCkBTlBJAJr6EdNZHNPqoRid66Biek08+WZ555hn3uicaAgggEAsBXYZM1wvUglCXidFxgfrfOsSElpwC5G8KQFNkEkAmvrh33rx5s7Rv3z780nd9Z68uA6TrQdIQQACBWAr8/PPPcscdd8j777/vNnv11VfLiy++KP/zP/8Ty92wrRgJkL8pAE2hRACZ+OLaefbs2XLXXXfJr7/+KgULFnSLOXfp0oWlXeKqzsYRCLaAji3Wt1HpItL67vhSpUq5lQZuuOGGYMMk4dmTvykATWFJAJn44tJZr0n37t1l0qRJbvu6ILiOyalWrVpc9sdGEUAAgcwC3377rbRu3VqWLVvm/qpt27buLUMnnXQSWEkiQP6mADSFIgFk4ot55wULFriJHmvXrnVjb/Rb+ODBg6VAgQIx3xcbRAABBI4noHcAdbkYfee8jkU+/fTTRVcduOSSS4BLAgHyNwWgKQwJIBNfzDofOHDA/aLVx7z6i1Zf86e/aC+99NKY7YMNIYAAAtEIfPLJJ+6Lqb56NG/evO6Lqa47yhfTaDRj14f8TQFoiiYCyMQXk8760nZ91PLll1+67ekgbH3Uom/2oCGAAALJIKC5omvXrvK///u/7nD0XeM6NOXss89OhsML5DGQvykATYFPAJn4TJ11sLWu46ffpvft2+cGW0+cOFGaNWtm2i6dEUAAgXgJvPbaa+71cTt27JBChQq5x8P33nsvk9PiBX6c7ZK/KQBNYUcAmfii7rxx40Y3qPrdd99127jqqqvccgunnHJK1NukIwIIIJAIAV2eql27dvL222+73V155ZXu9xfLUyVC///3Qf6mADRFHAFk4ouq84wZM9wCqzt37nQLrj7xxBPuXb4suBoVJ50QQMADAR2rPH78eOnVq5d7gqEL1OsTjJtvvtmDownmLsnfFICmyCeATHy56pyenu5euaTjZrTVqlXLvYD9rLPOytV2+DACCCCQLAIrV650Y5iXLl3qDkn/XYe2FC9ePFkO0bfHQf6mADQFNwFk4ou48/z58+X222+X9evXu7Ey/fr1k4ceeoiXrkcsyAcRQCBZBf788095+OGH5dFHHxUd25yWluZWMbj88suT9ZB9cVzkbwpAUyATQCa+HDvrOlr9+/d379XURyaVK1d2d/3q1auXY18+gAACCKSSwKJFi9xyMWvWrHFDWvTx8COPPOLeZESLvQD5mwLQFFUEkInvuJ2//vpr9zhk+fLl7nN33nmnPPXUU3LiiSfGb6dsGQEEEPBQYPfu3XLffffJ888/746ievXqbtjLeeed5+FR+XPX5G8KQFNkE0Amviw76yOQUaNGyQMPPCC6wHOZMmXcL8Prrrsu9jtjiwgggEASCsydO9d96d2+fbtbMHrYsGHuFZc6BIYWGwHyNwWgKZIIIBPfMZ11jJ8u5Pzhhx+6v2vSpIkr/sqVKxfbHbE1BBBAIMkFtm7dKu3bt5d//vOf7kgbNmzoFpKuUKFCkh95ahwe+ZsC0BSpBJCJ76jOU6dOlU6dOonO9i1SpIh73HvXXXexvEvsiNkSAgikmICOfX7uueekR48esnfvXjc7WJePadmyZYqdSfIdLvmbAtAUlQSQic911vX8tPCbNm2a+++6deu6iR5nnnmmfeNsAQEEEPCBwA8//OAmiCxevNidTYsWLVwhqOsH0qITIH9TAEYXOf/tRQCZ+OT999+XNm3aiL7ZI1++fG5pF13i5YQTTrBtmN4IIICAzwQOHjwoQ4cOdTODDx06JKeddppMnjxZrrjiCp+daWJOh/xNAWiKNAIoOr4//vjDFXo62UOb3u3TmW516tSJboP0QgABBAIioHcB9W6g3hXUpo+HdQ1BfbcwLXIB8jcFYOTRksUnCaDc8y1btkxatWolK1ascJ31NW4jRoyQokWL5n5j9EAAAQQCKPD777+7dQInTJjgzr5atWruS3TNmjUDqBHdKZO/KQCji5z/9iKAIufTRxb63t4BAwaIrnyvM3v1BeiNGzeOfCN8EgEEEEAgLKAzhHWmsM4Yzp8/vwwZMkR69uzphtTQji9A/qYANP2MEECR8a1du9a9yu2TTz5xHZo2bSrPPvusW+OPhgACCCAQvcC2bdvcigm6dqC2yy67zC0XU6lSpeg3GoCe5G8KQFOYE0DH59MlDPSdll26dBFd4V7f4jF69Gi31p++6oiGAAIIIGAX0N+1kyZNkm7dusmePXvkpJNOkrFjx7qxgvyuzdqX/E0BaPrJI4Cy5/v111/lnnvukddee819qH79+q4YPOOMM0zmdEYAAQQQyFpA3yOsRZ++V1hbs2bNZOLEiVKqVCnIMgmQvykATT8UBFDWfG+//ba0a9dONm/e7JZ0GTx4sPTp04dxKaZoozMCCCCQs4AuFzN8+HAZNGiQ6L+fcsop7u5go0aNcu4coE+QvykATeFOAB3NpyvV9+7dW8aNG+f+4uyzz3Yz0y644AKTM50RQAABBHInsHTpUmndurWsXLnSdezcubMrDPVNSzQR8jcFoOnngAD6f77ly5dL8+bNw79sdNyf/rIpXLiwyZjOCCCAAALRCeiXcn36ouMBtZ111lkyffp0qV69enQb9FEv8jcFoCmcCaD/8OnYvg4dOogu8MzjBlNI0RkBBBCIucA777wjbdu2dcNy9Eu5jgvUsYJBbuRvCkBT/Ac9gPbv3+9mnekvE21XXXWVe49v6dKlTa50RgABBBCIrcD27dvdI2EtBrXpl3Z9G1PBggVju6MU2VrQ87depjxHdP44LSqBIAfQunXr5KabbpLPP//cLTOgA4779+8vefPmjcqSTggggAAC8RXQBfl1sWidmKepv3bt2jJz5kypWLFifHechFsPcv4OXQ4KQENgBjWAdJavvs5tx44dUrJkSXnllVfc3T8aAggggEDyC2T+HT516tTAzRIOav7OGJ0UgIaf1aAF0OHDh+WRRx4Jf3usVauWzJo1K5DfHg1hQ1cEEEDAcwF9Q5M+xdHZwvoUZ+DAge5VnUF5ihO0/J1VwFEAGn4MgxRAurCzjh/Rb47agj5+xBA2dEUAAQSSQmDfvn3SvXv3o8Zx69JdQVg4Okj5O7tgowA0/BgGJYCWLFnivimuX7/ezSCbMGGCe7cvDQEEEEAg9QX03cH6pV4LQh0PqE929AmPn1tQ8vfxriEFoCHC/R5AOkj42Wefla5du8qBAwekSpUq7tVurCFlCBq6IoAAAkko8NVXX7lXx/34449SoEABGTNmjNx1112+fZew3/N3JCFGARiJUjaf8XMA6QKiHTt2dGv8abv++utl8uTJUqJECYMYXRFAAAEEklVg165d0qZNG/nHP/7hDlH/ffz48b58e4if83ek8UUBGKlUFp/zawCtXr3afRPUt3vogOBhw4bJ/fff79tvgoYQoCsCCCDgKwGd7DdixAjp16+f6L/rEx998qNPgPzU/Jq/c3ONKABzo5Xps34MoLlz57rxfXpuZcuWda8NatCggUGJrggggAACqSbw4YcfSosWLeSXX36R4sWLi44T1CdBfml+zN+5vTYUgLkVy/B5PwXQwYMH3ULO+v5ebfXr186APdQAABFPSURBVJcZM2bIqaeeahCiKwIIIIBAqgps3LhRbrnlFlm0aJE7hb59+7qlwE444YRUPaXwcfspf0d7MSgAo5UTcXfJ9JtRenq6FCtWzLAlb7tu3brVfdP76KOP3IHosgCPP/645M+f39sDY+8IIIAAAp4K/Pnnn9K7d2/32jhtDRs2lFdffVXKlSvn6XFZd+6X/G1xoAA06PkhgBYuXCg333yze0n4iSeeKC+88IL7xkdDAAEEEEAgJKBPhNq1aye///67ezKkr5C7+OKLUxbID/nbik8BaBBM5QDSJV5Gjx4tvXr1En38e/bZZ7uBvvpPGgIIIIAAApkFvvvuOzdBUP+pj4GffPJJ6dKlS0pOEEzl/B2ryKQANEimagDt3r1b7rzzTjfGT5s+/n3uuefcHUAaAggggAAC2Qns2bPH5Q+dIKitefPm8vzzz6dc/kjV/B3LyKQANGimYgCtWLHCfYNbuXKl+wY3cuRI6dy5c0p+gzNcOroigAACCEQpoE+QdKHonj17puwTpFTM31Fermy7UQAaRFMtgKZNm+a+uekYjvLly7s7gKk8hsNw6eiKAAIIIGAU0NnBOoZ806ZNUrRoUTeGXO8IpkJLtfwdD1MKQINqqgSQvsZNF3LWMX/a/vrXv7pZXLrOHw0BBBBAAIFoBXSdQB1GpOsGauvWrZtbRUJfJ5fMLVXydzwNKQANuqkQQJnXcdLV3R9++GHJly+f4czpigACCCCAwH8EdCLhQw895N4apU2fLOkTJn3SlKwtFfJ3vO0oAA3CyR5AH3zwgftmtm3bNrde4csvvyzXXnut4YzpigACCCCAQNYC+g5hfZOUro2rT5j0SZM+cUrGluz5OxFmFIAG5WQNIH1/o96Cf/DBB927HGvUqOGWeKlcubLhbOmKAAIIIIDA8QV+/PFHN9Hwq6++cu+SHzp0qFtIWv89mVqy5u9EGlEAGrSTMYB27dolbdq0Ef0mpu2OO+6Q8ePHS+HChQ1nSlcEEEAAAQQiE9i7d6906tTJvT9Y23XXXef+vUSJEpFtIAGfSsb8nYDTPmoXFIAG8WQLoGXLlslNN90k+g2sYMGCMnbsWGnfvj1LvBiuMV0RQAABBHIvoEvF6PqAusyYTkTUJ1D6JEqfSCVDS7b87YUJBaBBPZkCaPLkydKxY0fZt2+fVKpUSWbNmiUXXnih4ezoigACCCCAgE3g888/dzcm1q1bJ4UKFZIJEya4p1Ret2TK315ZUAAa5JMhgLTg69q1q3uTh7bGjRu7yR4lS5Y0nBldEUAAAQQQiI3Ajh07pHXr1vLWW2+5Dd59993y9NNPu4LQq5YM+durcw/tlwLQcAW8DqC1a9e6b1ZLly51j3kHDx7sJn4k22BbAzFdEUAAAQR8IKATEocMGSKDBg0SfTysT6j0SZU+sfKieZ2/vTjnzPukADRcBS8DSL9JtWrVSnbu3CmlSpWSqVOnypVXXmk4G7oigAACCCAQX4F33nlHWrZsKXpXUJ9UvfLKK3LVVVfFd6dZbN3L/J3wk81mhxSAhivhRQAdOnTILeT8yCOPuG9RderUkZkzZ0paWprhTOiKAAIIIIBAYgR0PKC+Qm7JkiXu6ZUuIj1gwICEvqDAi/ydGN3I90IBGLnVMZ9MdABt377d3fWbN2+eOxadZj9y5Eg345eGAAIIIIBAqgjs379fevToIc8884w75EaNGsmUKVOkdOnSCTmFROfvhJxULndCAZhLsIwfT2QAffbZZ26834YNG9yafs8++6wbVEtDAAEEEEAgVQV00uI999wjf/zxh3uSpeMCa9euHffTSWT+jvvJRLkDXxWA48aNkxEjRsiWLVvcWkNjxoxxj0iza/roVG8762SKM888U4YPH+5m0UbaEhFA+ph34sSJ7gXbupaSHqeupXTeeedFeph8DgEEEEAAgaQV+Prrr+XGG2+U1atXS4ECBWT06NFuprA+Ho5XS0T+jtexx2q7vikAp0+f7t5BqGsM1a1bV0aNGuXGxq1atcq9kzBzW7RokVx22WXu5dVNmjRxkyi0APziiy/k3HPPjcg33gGkq6l36NDBLeui7YYbbpBJkya59/rSEEAAAQQQ8IuAvj+4bdu2Mnv2bHdKms/18XCRIkXicorxzt9xOegYb9Q3BaAWfXrbWN9+oU2nnFeoUEG6dOkiffv2PYatefPm8vvvv8ubb74Z/ruLLrpIatas6YrISFo8A+iHH35w71PUb0b58uWTxx57THr27BnXb0SRnDOfQQABBBBAIB4C+sTrySefdDlbJzzqky594qVPvmLd4pm/Y32s8dqeLwpAfTSq3xJ07EDTpk3DVrrauL4bd+7cucf46ViD++67T7p37x7+u4EDB8qcOXPcS6yzajpoVf+Emn5j0e3ouLxixYrF7Bq98cYb7q0eu3fvljJlyoi+5eOSSy6J2fbZEAIIIIAAAskqsGDBAvce+23btslJJ53kbsrok7pYNi0A9SbR+vXrj3qqppMqgzKx0hcF4KZNm6R8+fKij3Xr1asXjpHevXvL/PnzZfHixcfEjY4z0JdT33rrreG/Gz9+vFtMeevWrVnGmS5gqX9PQwABBBBAAAH/CeiNIM31QWgUgLkoADPfAdTHzLqYpS7EHOvBqqFvJ7G+uxiEoM7NOeKcG63oP4tz9Ha56YlzbrSi/yzO0dvlpmc8nfVx86+//uoWo8749izuAObmCiXBZxP1CDiRp8r4hMRo44xzYgQSsxfiGefECCRmL8RzfJ19cQdQiXQSiC75oku/aNO7czo+r3PnztlOAtFZtjreLtQuvvhiqV69esSTQOJ5aQj8eOr+/7ZxxjkxAonZC/GMc2IEErMX4jm+zr4pAHUZGJ30oWvmaSGoy8DMmDFDVq5cKeXKlXNTynWcoC77ok3HC15++eVudu0111wj06ZNk0cffTRXy8DE89IQ+PHUpQBMjC7OOCdaIDH74/czzokRiO9efFMAKpMuARNaCFqXc9HFJPXOoLYGDRpIpUqV3IzaUNN1Avv37x9eCPrxxx/P1ULQ8bw0Ot5Qi9UHHnggMDOS4umZ3bZxTow6zjgnRiAxeyGecU6MQHz34qsCML5UbB0BBBBAAAEEEPCHAAWgP64jZ4EAAggggAACCEQsQAEYMRUfRAABBBBAAAEE/CFAAeiP68hZIIAAAggggAACEQtQAEZMxQcRQAABBBBAAAF/CFAAJuF1HDduXHg2c40aNdzahrq0DS02Ajq7+vXXX3dLBBUuXFh0/cfhw4fLX/7yl9jsgK1kKaBLLums9m7durllmmixFdi4caP06dNH3nrrLdE1TqtUqSKTJk2SWrVqxXZHAd7aoUOH3GvCpkyZIlu2bJFTTz3VvbNWV5OI9duggsT88ccfu5y3dOlS2bx5s8yePVuaNm0aJtC3dugr2p577jnZtWuX1K9fX5555hk588wzg8QU83OlAIw5qW2Dup6hrlmoL7/WJWw0UepyNatWrZKyZcvaNk5vJ3DVVVdJixYtpHbt2nLw4EHp16+ffPPNN7JixQopWrQoSnEQWLJkidxyyy1SrFgxadiwIQVgjI137twp559/vrPt2LGjlClTRn744QepXLmy+0OLjYCuFTty5Ej3Hvlq1arJ559/Lm3btpWhQ4dK165dY7OTAG5Fv7QsXLhQLrzwQrnxxhuPKQD1C7p+cVf3008/XQYMGCBff/21+51dqFChAIrF5pQpAGPjGLOtaNGnhYmuaahN32hSoUIF6dKlS5ZvNInZjgO8oW3btrniev78+XLZZZcFWCI+p75nzx654IILZPz48TJkyBDRNTq5Axhb6759+7oE+sknn8R2w2ztKIEmTZq4Fwu88MIL4f/frFkz9yRB7wrS7AJ6JzXjHUC9+6d3Wnv27Cm9evVyO0hPT3fXQdf11S/ztOgEKACjc4tLr2jeaRyXAwnYRlevXu0eJeg3ynPPPTdgZx//09U39OgL15966im3IDsFYOzNzznnHGnUqJH8/PPP7ouMvvWoU6dOctddd8V+ZwHeot4BfPbZZ2XevHlStWpV+eqrr+TKK690dwVbtWoVYJnYnXrmAnDNmjXuLvaXX37pfneEmr7JS//76aefjt3OA7YlCsAkuuCbNm1yv7j1NXX16tULH1nv3r3dL/XFixcn0dH641D0Dut1113nxpUsWLDAHyeVRGehr1jUx2P6CFgf1VAAxufihB6D3XfffXLzzTc7bx1rqUNJtACnxUZAf1/okBF9a1S+fPlExwRqfOvYVlpsBDIXgJoPdcyf5sdTTjklvBMdUqKf1WFTtOgEKACjc4tLLwrAuLAed6M6XkrHn2jxd9pppyX+AHy8xw0bNrgJCO+++65Ur17dnSkFYHwueIECBZy1JstQ0zFpWgh++umn8dlpALeqX2juv/9+N2FBxwAuW7ZMunfv7u4AUmjHJiAoAGPjGMlWKAAjUUrQZ3gEnCDo/+6mc+fOMnfuXNEZaDqwmBZbgTlz5sgNN9zg7pSEmt4x0V/wefPmFX2fasa/i+3eg7W1ihUryt///nd5/vnnwyeusyR1zKXODqbFRkDHY+t4y3vvvTe8QTXW8X+6qgDNLsAjYLthpFugAIxUKkGf00kguuSLLv2iTR85pKWliRYr+ouHZhfQQcU6qUYHGn/00UcsJWAnzXILu3fvlnXr1h31dzpj8qyzznLLlTDeMnbwLVu2FL3jmnESSI8ePdywkYx3BWO3x2BuqVSpUq6o1icHoaazU3W5ne+//z6YKDE+6+wmgegEEJ0Iou23335zE/eYBGLDpwC0+cW8t45n0EcJEydOdIWgzpacMWOG+3aps55odgEdHD916lR39y/j2n/Fixd3s/lo8RPgEXB8bPVRr65nOXjwYLfczmeffeYmgOiEBSYnxM5c1/x777333O9nfQSsExPuvvtuadeunVtLlBadgK4UoJPxtOlyRvpIXZc00sljegNEbXUd0YzLwCxfvpxlYKLjDveiADQCxqO7LgGjY0x0oVGd5TR69Gi3JiAtNgLZLdiq3+L1FzwtfgIUgPGzffPNN91kBF3/T4c06IQQZgHH1lvvausadPr04JdffnHLk9x6663y0EMPiY7DpEUnoE9itODL3PRmiN7lCy0ErV9odMLeJZdc4paV0pnYtOgFKACjt6MnAggggAACCCCQkgIUgCl52ThoBBBAAAEEEEAgegEKwOjt6IkAAggggAACCKSkAAVgSl42DhoBBBBAAAEEEIhegAIwejt6IoAAAggggAACKSlAAZiSl42DRgABBBBAAAEEohegAIzejp4IIIAAAggggEBKClAApuRl46ARQAABBBBAAIHoBSgAo7ejJwIIIIAAAgggkJICFIApedk4aAQQQAABBBBAIHoBCsDo7eiJAAIIIIAAAgikpAAFYEpeNg4aAQQQQAABBBCIXoACMHo7eiKAAAIIIIAAAikpQAGYkpeNg0YAAQQQQAABBKIXoACM3o6eCCCAAAIIIIBASgpQAKbkZeOgEUAAAQQQQACB6AUoAKO3oycCCCCAAAIIIJCSAhSAKXnZOGgEEEAAAQQQQCB6AQrA6O3oiQACCCCAAAIIpKQABWBKXjYOGgEEEEAAAQQQiF6AAjB6O3oigAACCCCAAAIpKUABmJKXjYNGAAEEEEAAAQSiF6AAjN6OnggggAACCCCAQEoKUACm5GXjoBFAAAEEEEAAgegFKACjt6MnAggggAACCCCQkgIUgCl52ThoBBBAAAEEEEAgegEKwOjt6IkAAggggAACCKSkAAVgSl42DhoBBBBAAAEEEIhegAIwejt6IoAAAggggAACKSlAAZiSl42DRgABBBBAAAEEohegAIzejp4IIIAAAggggEBKClAApuRl46ARQAABBBBAAIHoBf4PCsFoE3fF5IkAAAAASUVORK5CYII=\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "<matplotlib.animation.FuncAnimation at 0x7fb8ac5a4be0>"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "animate_line(u, t, sampling=5)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "4eb57be9",
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "/* global mpl */\n",
       "window.mpl = {};\n",
       "\n",
       "mpl.get_websocket_type = function () {\n",
       "    if (typeof WebSocket !== 'undefined') {\n",
       "        return WebSocket;\n",
       "    } else if (typeof MozWebSocket !== 'undefined') {\n",
       "        return MozWebSocket;\n",
       "    } else {\n",
       "        alert(\n",
       "            'Your browser does not have WebSocket support. ' +\n",
       "                'Please try Chrome, Safari or Firefox ≥ 6. ' +\n",
       "                'Firefox 4 and 5 are also supported but you ' +\n",
       "                'have to enable WebSockets in about:config.'\n",
       "        );\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure = function (figure_id, websocket, ondownload, parent_element) {\n",
       "    this.id = figure_id;\n",
       "\n",
       "    this.ws = websocket;\n",
       "\n",
       "    this.supports_binary = this.ws.binaryType !== undefined;\n",
       "\n",
       "    if (!this.supports_binary) {\n",
       "        var warnings = document.getElementById('mpl-warnings');\n",
       "        if (warnings) {\n",
       "            warnings.style.display = 'block';\n",
       "            warnings.textContent =\n",
       "                'This browser does not support binary websocket messages. ' +\n",
       "                'Performance may be slow.';\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.imageObj = new Image();\n",
       "\n",
       "    this.context = undefined;\n",
       "    this.message = undefined;\n",
       "    this.canvas = undefined;\n",
       "    this.rubberband_canvas = undefined;\n",
       "    this.rubberband_context = undefined;\n",
       "    this.format_dropdown = undefined;\n",
       "\n",
       "    this.image_mode = 'full';\n",
       "\n",
       "    this.root = document.createElement('div');\n",
       "    this.root.setAttribute('style', 'display: inline-block');\n",
       "    this._root_extra_style(this.root);\n",
       "\n",
       "    parent_element.appendChild(this.root);\n",
       "\n",
       "    this._init_header(this);\n",
       "    this._init_canvas(this);\n",
       "    this._init_toolbar(this);\n",
       "\n",
       "    var fig = this;\n",
       "\n",
       "    this.waiting = false;\n",
       "\n",
       "    this.ws.onopen = function () {\n",
       "        fig.send_message('supports_binary', { value: fig.supports_binary });\n",
       "        fig.send_message('send_image_mode', {});\n",
       "        if (fig.ratio !== 1) {\n",
       "            fig.send_message('set_dpi_ratio', { dpi_ratio: fig.ratio });\n",
       "        }\n",
       "        fig.send_message('refresh', {});\n",
       "    };\n",
       "\n",
       "    this.imageObj.onload = function () {\n",
       "        if (fig.image_mode === 'full') {\n",
       "            // Full images could contain transparency (where diff images\n",
       "            // almost always do), so we need to clear the canvas so that\n",
       "            // there is no ghosting.\n",
       "            fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
       "        }\n",
       "        fig.context.drawImage(fig.imageObj, 0, 0);\n",
       "    };\n",
       "\n",
       "    this.imageObj.onunload = function () {\n",
       "        fig.ws.close();\n",
       "    };\n",
       "\n",
       "    this.ws.onmessage = this._make_on_message_function(this);\n",
       "\n",
       "    this.ondownload = ondownload;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._init_header = function () {\n",
       "    var titlebar = document.createElement('div');\n",
       "    titlebar.classList =\n",
       "        'ui-dialog-titlebar ui-widget-header ui-corner-all ui-helper-clearfix';\n",
       "    var titletext = document.createElement('div');\n",
       "    titletext.classList = 'ui-dialog-title';\n",
       "    titletext.setAttribute(\n",
       "        'style',\n",
       "        'width: 100%; text-align: center; padding: 3px;'\n",
       "    );\n",
       "    titlebar.appendChild(titletext);\n",
       "    this.root.appendChild(titlebar);\n",
       "    this.header = titletext;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function (_canvas_div) {};\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function (_canvas_div) {};\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = (this.canvas_div = document.createElement('div'));\n",
       "    canvas_div.setAttribute(\n",
       "        'style',\n",
       "        'border: 1px solid #ddd;' +\n",
       "            'box-sizing: content-box;' +\n",
       "            'clear: both;' +\n",
       "            'min-height: 1px;' +\n",
       "            'min-width: 1px;' +\n",
       "            'outline: 0;' +\n",
       "            'overflow: hidden;' +\n",
       "            'position: relative;' +\n",
       "            'resize: both;'\n",
       "    );\n",
       "\n",
       "    function on_keyboard_event_closure(name) {\n",
       "        return function (event) {\n",
       "            return fig.key_event(event, name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    canvas_div.addEventListener(\n",
       "        'keydown',\n",
       "        on_keyboard_event_closure('key_press')\n",
       "    );\n",
       "    canvas_div.addEventListener(\n",
       "        'keyup',\n",
       "        on_keyboard_event_closure('key_release')\n",
       "    );\n",
       "\n",
       "    this._canvas_extra_style(canvas_div);\n",
       "    this.root.appendChild(canvas_div);\n",
       "\n",
       "    var canvas = (this.canvas = document.createElement('canvas'));\n",
       "    canvas.classList.add('mpl-canvas');\n",
       "    canvas.setAttribute('style', 'box-sizing: content-box;');\n",
       "\n",
       "    this.context = canvas.getContext('2d');\n",
       "\n",
       "    var backingStore =\n",
       "        this.context.backingStorePixelRatio ||\n",
       "        this.context.webkitBackingStorePixelRatio ||\n",
       "        this.context.mozBackingStorePixelRatio ||\n",
       "        this.context.msBackingStorePixelRatio ||\n",
       "        this.context.oBackingStorePixelRatio ||\n",
       "        this.context.backingStorePixelRatio ||\n",
       "        1;\n",
       "\n",
       "    this.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
       "\n",
       "    var rubberband_canvas = (this.rubberband_canvas = document.createElement(\n",
       "        'canvas'\n",
       "    ));\n",
       "    rubberband_canvas.setAttribute(\n",
       "        'style',\n",
       "        'box-sizing: content-box; position: absolute; left: 0; top: 0; z-index: 1;'\n",
       "    );\n",
       "\n",
       "    // Apply a ponyfill if ResizeObserver is not implemented by browser.\n",
       "    if (this.ResizeObserver === undefined) {\n",
       "        if (window.ResizeObserver !== undefined) {\n",
       "            this.ResizeObserver = window.ResizeObserver;\n",
       "        } else {\n",
       "            var obs = _JSXTOOLS_RESIZE_OBSERVER({});\n",
       "            this.ResizeObserver = obs.ResizeObserver;\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.resizeObserverInstance = new this.ResizeObserver(function (entries) {\n",
       "        var nentries = entries.length;\n",
       "        for (var i = 0; i < nentries; i++) {\n",
       "            var entry = entries[i];\n",
       "            var width, height;\n",
       "            if (entry.contentBoxSize) {\n",
       "                if (entry.contentBoxSize instanceof Array) {\n",
       "                    // Chrome 84 implements new version of spec.\n",
       "                    width = entry.contentBoxSize[0].inlineSize;\n",
       "                    height = entry.contentBoxSize[0].blockSize;\n",
       "                } else {\n",
       "                    // Firefox implements old version of spec.\n",
       "                    width = entry.contentBoxSize.inlineSize;\n",
       "                    height = entry.contentBoxSize.blockSize;\n",
       "                }\n",
       "            } else {\n",
       "                // Chrome <84 implements even older version of spec.\n",
       "                width = entry.contentRect.width;\n",
       "                height = entry.contentRect.height;\n",
       "            }\n",
       "\n",
       "            // Keep the size of the canvas and rubber band canvas in sync with\n",
       "            // the canvas container.\n",
       "            if (entry.devicePixelContentBoxSize) {\n",
       "                // Chrome 84 implements new version of spec.\n",
       "                canvas.setAttribute(\n",
       "                    'width',\n",
       "                    entry.devicePixelContentBoxSize[0].inlineSize\n",
       "                );\n",
       "                canvas.setAttribute(\n",
       "                    'height',\n",
       "                    entry.devicePixelContentBoxSize[0].blockSize\n",
       "                );\n",
       "            } else {\n",
       "                canvas.setAttribute('width', width * fig.ratio);\n",
       "                canvas.setAttribute('height', height * fig.ratio);\n",
       "            }\n",
       "            canvas.setAttribute(\n",
       "                'style',\n",
       "                'width: ' + width + 'px; height: ' + height + 'px;'\n",
       "            );\n",
       "\n",
       "            rubberband_canvas.setAttribute('width', width);\n",
       "            rubberband_canvas.setAttribute('height', height);\n",
       "\n",
       "            // And update the size in Python. We ignore the initial 0/0 size\n",
       "            // that occurs as the element is placed into the DOM, which should\n",
       "            // otherwise not happen due to the minimum size styling.\n",
       "            if (fig.ws.readyState == 1 && width != 0 && height != 0) {\n",
       "                fig.request_resize(width, height);\n",
       "            }\n",
       "        }\n",
       "    });\n",
       "    this.resizeObserverInstance.observe(canvas_div);\n",
       "\n",
       "    function on_mouse_event_closure(name) {\n",
       "        return function (event) {\n",
       "            return fig.mouse_event(event, name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mousedown',\n",
       "        on_mouse_event_closure('button_press')\n",
       "    );\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseup',\n",
       "        on_mouse_event_closure('button_release')\n",
       "    );\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'dblclick',\n",
       "        on_mouse_event_closure('dblclick')\n",
       "    );\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mousemove',\n",
       "        on_mouse_event_closure('motion_notify')\n",
       "    );\n",
       "\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseenter',\n",
       "        on_mouse_event_closure('figure_enter')\n",
       "    );\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseleave',\n",
       "        on_mouse_event_closure('figure_leave')\n",
       "    );\n",
       "\n",
       "    canvas_div.addEventListener('wheel', function (event) {\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        on_mouse_event_closure('scroll')(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.appendChild(canvas);\n",
       "    canvas_div.appendChild(rubberband_canvas);\n",
       "\n",
       "    this.rubberband_context = rubberband_canvas.getContext('2d');\n",
       "    this.rubberband_context.strokeStyle = '#000000';\n",
       "\n",
       "    this._resize_canvas = function (width, height, forward) {\n",
       "        if (forward) {\n",
       "            canvas_div.style.width = width + 'px';\n",
       "            canvas_div.style.height = height + 'px';\n",
       "        }\n",
       "    };\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    this.rubberband_canvas.addEventListener('contextmenu', function (_e) {\n",
       "        event.preventDefault();\n",
       "        return false;\n",
       "    });\n",
       "\n",
       "    function set_focus() {\n",
       "        canvas.focus();\n",
       "        canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    window.setTimeout(set_focus, 100);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var toolbar = document.createElement('div');\n",
       "    toolbar.classList = 'mpl-toolbar';\n",
       "    this.root.appendChild(toolbar);\n",
       "\n",
       "    function on_click_closure(name) {\n",
       "        return function (_event) {\n",
       "            return fig.toolbar_button_onclick(name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    function on_mouseover_closure(tooltip) {\n",
       "        return function (event) {\n",
       "            if (!event.currentTarget.disabled) {\n",
       "                return fig.toolbar_button_onmouseover(tooltip);\n",
       "            }\n",
       "        };\n",
       "    }\n",
       "\n",
       "    fig.buttons = {};\n",
       "    var buttonGroup = document.createElement('div');\n",
       "    buttonGroup.classList = 'mpl-button-group';\n",
       "    for (var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            /* Instead of a spacer, we start a new button group. */\n",
       "            if (buttonGroup.hasChildNodes()) {\n",
       "                toolbar.appendChild(buttonGroup);\n",
       "            }\n",
       "            buttonGroup = document.createElement('div');\n",
       "            buttonGroup.classList = 'mpl-button-group';\n",
       "            continue;\n",
       "        }\n",
       "\n",
       "        var button = (fig.buttons[name] = document.createElement('button'));\n",
       "        button.classList = 'mpl-widget';\n",
       "        button.setAttribute('role', 'button');\n",
       "        button.setAttribute('aria-disabled', 'false');\n",
       "        button.addEventListener('click', on_click_closure(method_name));\n",
       "        button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n",
       "\n",
       "        var icon_img = document.createElement('img');\n",
       "        icon_img.src = '_images/' + image + '.png';\n",
       "        icon_img.srcset = '_images/' + image + '_large.png 2x';\n",
       "        icon_img.alt = tooltip;\n",
       "        button.appendChild(icon_img);\n",
       "\n",
       "        buttonGroup.appendChild(button);\n",
       "    }\n",
       "\n",
       "    if (buttonGroup.hasChildNodes()) {\n",
       "        toolbar.appendChild(buttonGroup);\n",
       "    }\n",
       "\n",
       "    var fmt_picker = document.createElement('select');\n",
       "    fmt_picker.classList = 'mpl-widget';\n",
       "    toolbar.appendChild(fmt_picker);\n",
       "    this.format_dropdown = fmt_picker;\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = document.createElement('option');\n",
       "        option.selected = fmt === mpl.default_extension;\n",
       "        option.innerHTML = fmt;\n",
       "        fmt_picker.appendChild(option);\n",
       "    }\n",
       "\n",
       "    var status_bar = document.createElement('span');\n",
       "    status_bar.classList = 'mpl-message';\n",
       "    toolbar.appendChild(status_bar);\n",
       "    this.message = status_bar;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.request_resize = function (x_pixels, y_pixels) {\n",
       "    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n",
       "    // which will in turn request a refresh of the image.\n",
       "    this.send_message('resize', { width: x_pixels, height: y_pixels });\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.send_message = function (type, properties) {\n",
       "    properties['type'] = type;\n",
       "    properties['figure_id'] = this.id;\n",
       "    this.ws.send(JSON.stringify(properties));\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.send_draw_message = function () {\n",
       "    if (!this.waiting) {\n",
       "        this.waiting = true;\n",
       "        this.ws.send(JSON.stringify({ type: 'draw', figure_id: this.id }));\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_save = function (fig, _msg) {\n",
       "    var format_dropdown = fig.format_dropdown;\n",
       "    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n",
       "    fig.ondownload(fig, format);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_resize = function (fig, msg) {\n",
       "    var size = msg['size'];\n",
       "    if (size[0] !== fig.canvas.width || size[1] !== fig.canvas.height) {\n",
       "        fig._resize_canvas(size[0], size[1], msg['forward']);\n",
       "        fig.send_message('refresh', {});\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function (fig, msg) {\n",
       "    var x0 = msg['x0'] / fig.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / fig.ratio;\n",
       "    var x1 = msg['x1'] / fig.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / fig.ratio;\n",
       "    x0 = Math.floor(x0) + 0.5;\n",
       "    y0 = Math.floor(y0) + 0.5;\n",
       "    x1 = Math.floor(x1) + 0.5;\n",
       "    y1 = Math.floor(y1) + 0.5;\n",
       "    var min_x = Math.min(x0, x1);\n",
       "    var min_y = Math.min(y0, y1);\n",
       "    var width = Math.abs(x1 - x0);\n",
       "    var height = Math.abs(y1 - y0);\n",
       "\n",
       "    fig.rubberband_context.clearRect(\n",
       "        0,\n",
       "        0,\n",
       "        fig.canvas.width / fig.ratio,\n",
       "        fig.canvas.height / fig.ratio\n",
       "    );\n",
       "\n",
       "    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_figure_label = function (fig, msg) {\n",
       "    // Updates the figure title.\n",
       "    fig.header.textContent = msg['label'];\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_cursor = function (fig, msg) {\n",
       "    var cursor = msg['cursor'];\n",
       "    switch (cursor) {\n",
       "        case 0:\n",
       "            cursor = 'pointer';\n",
       "            break;\n",
       "        case 1:\n",
       "            cursor = 'default';\n",
       "            break;\n",
       "        case 2:\n",
       "            cursor = 'crosshair';\n",
       "            break;\n",
       "        case 3:\n",
       "            cursor = 'move';\n",
       "            break;\n",
       "    }\n",
       "    fig.rubberband_canvas.style.cursor = cursor;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_message = function (fig, msg) {\n",
       "    fig.message.textContent = msg['message'];\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_draw = function (fig, _msg) {\n",
       "    // Request the server to send over a new figure.\n",
       "    fig.send_draw_message();\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_image_mode = function (fig, msg) {\n",
       "    fig.image_mode = msg['mode'];\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_history_buttons = function (fig, msg) {\n",
       "    for (var key in msg) {\n",
       "        if (!(key in fig.buttons)) {\n",
       "            continue;\n",
       "        }\n",
       "        fig.buttons[key].disabled = !msg[key];\n",
       "        fig.buttons[key].setAttribute('aria-disabled', !msg[key]);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_navigate_mode = function (fig, msg) {\n",
       "    if (msg['mode'] === 'PAN') {\n",
       "        fig.buttons['Pan'].classList.add('active');\n",
       "        fig.buttons['Zoom'].classList.remove('active');\n",
       "    } else if (msg['mode'] === 'ZOOM') {\n",
       "        fig.buttons['Pan'].classList.remove('active');\n",
       "        fig.buttons['Zoom'].classList.add('active');\n",
       "    } else {\n",
       "        fig.buttons['Pan'].classList.remove('active');\n",
       "        fig.buttons['Zoom'].classList.remove('active');\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function () {\n",
       "    // Called whenever the canvas gets updated.\n",
       "    this.send_message('ack', {});\n",
       "};\n",
       "\n",
       "// A function to construct a web socket function for onmessage handling.\n",
       "// Called in the figure constructor.\n",
       "mpl.figure.prototype._make_on_message_function = function (fig) {\n",
       "    return function socket_on_message(evt) {\n",
       "        if (evt.data instanceof Blob) {\n",
       "            var img = evt.data;\n",
       "            if (img.type !== 'image/png') {\n",
       "                /* FIXME: We get \"Resource interpreted as Image but\n",
       "                 * transferred with MIME type text/plain:\" errors on\n",
       "                 * Chrome.  But how to set the MIME type?  It doesn't seem\n",
       "                 * to be part of the websocket stream */\n",
       "                img.type = 'image/png';\n",
       "            }\n",
       "\n",
       "            /* Free the memory for the previous frames */\n",
       "            if (fig.imageObj.src) {\n",
       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
       "                    fig.imageObj.src\n",
       "                );\n",
       "            }\n",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                img\n",
       "            );\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        } else if (\n",
       "            typeof evt.data === 'string' &&\n",
       "            evt.data.slice(0, 21) === 'data:image/png;base64'\n",
       "        ) {\n",
       "            fig.imageObj.src = evt.data;\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        var msg = JSON.parse(evt.data);\n",
       "        var msg_type = msg['type'];\n",
       "\n",
       "        // Call the  \"handle_{type}\" callback, which takes\n",
       "        // the figure and JSON message as its only arguments.\n",
       "        try {\n",
       "            var callback = fig['handle_' + msg_type];\n",
       "        } catch (e) {\n",
       "            console.log(\n",
       "                \"No handler for the '\" + msg_type + \"' message type: \",\n",
       "                msg\n",
       "            );\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        if (callback) {\n",
       "            try {\n",
       "                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n",
       "                callback(fig, msg);\n",
       "            } catch (e) {\n",
       "                console.log(\n",
       "                    \"Exception inside the 'handler_\" + msg_type + \"' callback:\",\n",
       "                    e,\n",
       "                    e.stack,\n",
       "                    msg\n",
       "                );\n",
       "            }\n",
       "        }\n",
       "    };\n",
       "};\n",
       "\n",
       "// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n",
       "mpl.findpos = function (e) {\n",
       "    //this section is from http://www.quirksmode.org/js/events_properties.html\n",
       "    var targ;\n",
       "    if (!e) {\n",
       "        e = window.event;\n",
       "    }\n",
       "    if (e.target) {\n",
       "        targ = e.target;\n",
       "    } else if (e.srcElement) {\n",
       "        targ = e.srcElement;\n",
       "    }\n",
       "    if (targ.nodeType === 3) {\n",
       "        // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "    }\n",
       "\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    var boundingRect = targ.getBoundingClientRect();\n",
       "    var x = e.pageX - (boundingRect.left + document.body.scrollLeft);\n",
       "    var y = e.pageY - (boundingRect.top + document.body.scrollTop);\n",
       "\n",
       "    return { x: x, y: y };\n",
       "};\n",
       "\n",
       "/*\n",
       " * return a copy of an object with only non-object keys\n",
       " * we need this to avoid circular references\n",
       " * http://stackoverflow.com/a/24161582/3208463\n",
       " */\n",
       "function simpleKeys(original) {\n",
       "    return Object.keys(original).reduce(function (obj, key) {\n",
       "        if (typeof original[key] !== 'object') {\n",
       "            obj[key] = original[key];\n",
       "        }\n",
       "        return obj;\n",
       "    }, {});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.mouse_event = function (event, name) {\n",
       "    var canvas_pos = mpl.findpos(event);\n",
       "\n",
       "    if (name === 'button_press') {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x * this.ratio;\n",
       "    var y = canvas_pos.y * this.ratio;\n",
       "\n",
       "    this.send_message(name, {\n",
       "        x: x,\n",
       "        y: y,\n",
       "        button: event.button,\n",
       "        step: event.step,\n",
       "        guiEvent: simpleKeys(event),\n",
       "    });\n",
       "\n",
       "    /* This prevents the web browser from automatically changing to\n",
       "     * the text insertion cursor when the button is pressed.  We want\n",
       "     * to control all of the cursor setting manually through the\n",
       "     * 'cursor' event from matplotlib */\n",
       "    event.preventDefault();\n",
       "    return false;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function (_event, _name) {\n",
       "    // Handle any extra behaviour associated with a key event\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.key_event = function (event, name) {\n",
       "    // Prevent repeat events\n",
       "    if (name === 'key_press') {\n",
       "        if (event.key === this._key) {\n",
       "            return;\n",
       "        } else {\n",
       "            this._key = event.key;\n",
       "        }\n",
       "    }\n",
       "    if (name === 'key_release') {\n",
       "        this._key = null;\n",
       "    }\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.key !== 'Control') {\n",
       "        value += 'ctrl+';\n",
       "    }\n",
       "    else if (event.altKey && event.key !== 'Alt') {\n",
       "        value += 'alt+';\n",
       "    }\n",
       "    else if (event.shiftKey && event.key !== 'Shift') {\n",
       "        value += 'shift+';\n",
       "    }\n",
       "\n",
       "    value += 'k' + event.key;\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, { key: value, guiEvent: simpleKeys(event) });\n",
       "    return false;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onclick = function (name) {\n",
       "    if (name === 'download') {\n",
       "        this.handle_save(this, null);\n",
       "    } else {\n",
       "        this.send_message('toolbar_button', { name: name });\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onmouseover = function (tooltip) {\n",
       "    this.message.textContent = tooltip;\n",
       "};\n",
       "\n",
       "///////////////// REMAINING CONTENT GENERATED BY embed_js.py /////////////////\n",
       "// prettier-ignore\n",
       "var _JSXTOOLS_RESIZE_OBSERVER=function(A){var t,i=new WeakMap,n=new WeakMap,a=new WeakMap,r=new WeakMap,o=new Set;function s(e){if(!(this instanceof s))throw new TypeError(\"Constructor requires 'new' operator\");i.set(this,e)}function h(){throw new TypeError(\"Function is not a constructor\")}function c(e,t,i,n){e=0 in arguments?Number(arguments[0]):0,t=1 in arguments?Number(arguments[1]):0,i=2 in arguments?Number(arguments[2]):0,n=3 in arguments?Number(arguments[3]):0,this.right=(this.x=this.left=e)+(this.width=i),this.bottom=(this.y=this.top=t)+(this.height=n),Object.freeze(this)}function d(){t=requestAnimationFrame(d);var s=new WeakMap,p=new Set;o.forEach((function(t){r.get(t).forEach((function(i){var r=t instanceof window.SVGElement,o=a.get(t),d=r?0:parseFloat(o.paddingTop),f=r?0:parseFloat(o.paddingRight),l=r?0:parseFloat(o.paddingBottom),u=r?0:parseFloat(o.paddingLeft),g=r?0:parseFloat(o.borderTopWidth),m=r?0:parseFloat(o.borderRightWidth),w=r?0:parseFloat(o.borderBottomWidth),b=u+f,F=d+l,v=(r?0:parseFloat(o.borderLeftWidth))+m,W=g+w,y=r?0:t.offsetHeight-W-t.clientHeight,E=r?0:t.offsetWidth-v-t.clientWidth,R=b+v,z=F+W,M=r?t.width:parseFloat(o.width)-R-E,O=r?t.height:parseFloat(o.height)-z-y;if(n.has(t)){var k=n.get(t);if(k[0]===M&&k[1]===O)return}n.set(t,[M,O]);var S=Object.create(h.prototype);S.target=t,S.contentRect=new c(u,d,M,O),s.has(i)||(s.set(i,[]),p.add(i)),s.get(i).push(S)}))})),p.forEach((function(e){i.get(e).call(e,s.get(e),e)}))}return s.prototype.observe=function(i){if(i instanceof window.Element){r.has(i)||(r.set(i,new Set),o.add(i),a.set(i,window.getComputedStyle(i)));var n=r.get(i);n.has(this)||n.add(this),cancelAnimationFrame(t),t=requestAnimationFrame(d)}},s.prototype.unobserve=function(i){if(i instanceof window.Element&&r.has(i)){var n=r.get(i);n.has(this)&&(n.delete(this),n.size||(r.delete(i),o.delete(i))),n.size||r.delete(i),o.size||cancelAnimationFrame(t)}},A.DOMRectReadOnly=c,A.ResizeObserver=s,A.ResizeObserverEntry=h,A}; // eslint-disable-line\n",
       "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Left button pans, Right button zooms\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"jpeg\", \"pgf\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n",
       "\n",
       "mpl.default_extension = \"png\";/* global mpl */\n",
       "\n",
       "var comm_websocket_adapter = function (comm) {\n",
       "    // Create a \"websocket\"-like object which calls the given IPython comm\n",
       "    // object with the appropriate methods. Currently this is a non binary\n",
       "    // socket, so there is still some room for performance tuning.\n",
       "    var ws = {};\n",
       "\n",
       "    ws.binaryType = comm.kernel.ws.binaryType;\n",
       "    ws.readyState = comm.kernel.ws.readyState;\n",
       "    function updateReadyState(_event) {\n",
       "        if (comm.kernel.ws) {\n",
       "            ws.readyState = comm.kernel.ws.readyState;\n",
       "        } else {\n",
       "            ws.readyState = 3; // Closed state.\n",
       "        }\n",
       "    }\n",
       "    comm.kernel.ws.addEventListener('open', updateReadyState);\n",
       "    comm.kernel.ws.addEventListener('close', updateReadyState);\n",
       "    comm.kernel.ws.addEventListener('error', updateReadyState);\n",
       "\n",
       "    ws.close = function () {\n",
       "        comm.close();\n",
       "    };\n",
       "    ws.send = function (m) {\n",
       "        //console.log('sending', m);\n",
       "        comm.send(m);\n",
       "    };\n",
       "    // Register the callback with on_msg.\n",
       "    comm.on_msg(function (msg) {\n",
       "        //console.log('receiving', msg['content']['data'], msg);\n",
       "        var data = msg['content']['data'];\n",
       "        if (data['blob'] !== undefined) {\n",
       "            data = {\n",
       "                data: new Blob(msg['buffers'], { type: data['blob'] }),\n",
       "            };\n",
       "        }\n",
       "        // Pass the mpl event to the overridden (by mpl) onmessage function.\n",
       "        ws.onmessage(data);\n",
       "    });\n",
       "    return ws;\n",
       "};\n",
       "\n",
       "mpl.mpl_figure_comm = function (comm, msg) {\n",
       "    // This is the function which gets called when the mpl process\n",
       "    // starts-up an IPython Comm through the \"matplotlib\" channel.\n",
       "\n",
       "    var id = msg.content.data.id;\n",
       "    // Get hold of the div created by the display call when the Comm\n",
       "    // socket was opened in Python.\n",
       "    var element = document.getElementById(id);\n",
       "    var ws_proxy = comm_websocket_adapter(comm);\n",
       "\n",
       "    function ondownload(figure, _format) {\n",
       "        window.open(figure.canvas.toDataURL());\n",
       "    }\n",
       "\n",
       "    var fig = new mpl.figure(id, ws_proxy, ondownload, element);\n",
       "\n",
       "    // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n",
       "    // web socket which is closed, not our websocket->open comm proxy.\n",
       "    ws_proxy.onopen();\n",
       "\n",
       "    fig.parent_element = element;\n",
       "    fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n",
       "    if (!fig.cell_info) {\n",
       "        console.error('Failed to find cell for figure', id, fig);\n",
       "        return;\n",
       "    }\n",
       "    fig.cell_info[0].output_area.element.on(\n",
       "        'cleared',\n",
       "        { fig: fig },\n",
       "        fig._remove_fig_handler\n",
       "    );\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_close = function (fig, msg) {\n",
       "    var width = fig.canvas.width / fig.ratio;\n",
       "    fig.cell_info[0].output_area.element.off(\n",
       "        'cleared',\n",
       "        fig._remove_fig_handler\n",
       "    );\n",
       "    fig.resizeObserverInstance.unobserve(fig.canvas_div);\n",
       "\n",
       "    // Update the output cell to use the data from the current canvas.\n",
       "    fig.push_to_output();\n",
       "    var dataURL = fig.canvas.toDataURL();\n",
       "    // Re-enable the keyboard manager in IPython - without this line, in FF,\n",
       "    // the notebook keyboard shortcuts fail.\n",
       "    IPython.keyboard_manager.enable();\n",
       "    fig.parent_element.innerHTML =\n",
       "        '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "    fig.close_ws(fig, msg);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.close_ws = function (fig, msg) {\n",
       "    fig.send_message('closing', msg);\n",
       "    // fig.ws.close()\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.push_to_output = function (_remove_interactive) {\n",
       "    // Turn the data on the canvas into data in the output cell.\n",
       "    var width = this.canvas.width / this.ratio;\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] =\n",
       "        '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function () {\n",
       "    // Tell IPython that the notebook contents must change.\n",
       "    IPython.notebook.set_dirty(true);\n",
       "    this.send_message('ack', {});\n",
       "    var fig = this;\n",
       "    // Wait a second, then push the new image to the DOM so\n",
       "    // that it is saved nicely (might be nice to debounce this).\n",
       "    setTimeout(function () {\n",
       "        fig.push_to_output();\n",
       "    }, 1000);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var toolbar = document.createElement('div');\n",
       "    toolbar.classList = 'btn-toolbar';\n",
       "    this.root.appendChild(toolbar);\n",
       "\n",
       "    function on_click_closure(name) {\n",
       "        return function (_event) {\n",
       "            return fig.toolbar_button_onclick(name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    function on_mouseover_closure(tooltip) {\n",
       "        return function (event) {\n",
       "            if (!event.currentTarget.disabled) {\n",
       "                return fig.toolbar_button_onmouseover(tooltip);\n",
       "            }\n",
       "        };\n",
       "    }\n",
       "\n",
       "    fig.buttons = {};\n",
       "    var buttonGroup = document.createElement('div');\n",
       "    buttonGroup.classList = 'btn-group';\n",
       "    var button;\n",
       "    for (var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            /* Instead of a spacer, we start a new button group. */\n",
       "            if (buttonGroup.hasChildNodes()) {\n",
       "                toolbar.appendChild(buttonGroup);\n",
       "            }\n",
       "            buttonGroup = document.createElement('div');\n",
       "            buttonGroup.classList = 'btn-group';\n",
       "            continue;\n",
       "        }\n",
       "\n",
       "        button = fig.buttons[name] = document.createElement('button');\n",
       "        button.classList = 'btn btn-default';\n",
       "        button.href = '#';\n",
       "        button.title = name;\n",
       "        button.innerHTML = '<i class=\"fa ' + image + ' fa-lg\"></i>';\n",
       "        button.addEventListener('click', on_click_closure(method_name));\n",
       "        button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n",
       "        buttonGroup.appendChild(button);\n",
       "    }\n",
       "\n",
       "    if (buttonGroup.hasChildNodes()) {\n",
       "        toolbar.appendChild(buttonGroup);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = document.createElement('span');\n",
       "    status_bar.classList = 'mpl-message pull-right';\n",
       "    toolbar.appendChild(status_bar);\n",
       "    this.message = status_bar;\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = document.createElement('div');\n",
       "    buttongrp.classList = 'btn-group inline pull-right';\n",
       "    button = document.createElement('button');\n",
       "    button.classList = 'btn btn-mini btn-primary';\n",
       "    button.href = '#';\n",
       "    button.title = 'Stop Interaction';\n",
       "    button.innerHTML = '<i class=\"fa fa-power-off icon-remove icon-large\"></i>';\n",
       "    button.addEventListener('click', function (_evt) {\n",
       "        fig.handle_close(fig, {});\n",
       "    });\n",
       "    button.addEventListener(\n",
       "        'mouseover',\n",
       "        on_mouseover_closure('Stop Interaction')\n",
       "    );\n",
       "    buttongrp.appendChild(button);\n",
       "    var titlebar = this.root.querySelector('.ui-dialog-titlebar');\n",
       "    titlebar.insertBefore(buttongrp, titlebar.firstChild);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._remove_fig_handler = function (event) {\n",
       "    var fig = event.data.fig;\n",
       "    if (event.target !== this) {\n",
       "        // Ignore bubbled events from children.\n",
       "        return;\n",
       "    }\n",
       "    fig.close_ws(fig, {});\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function (el) {\n",
       "    el.style.boxSizing = 'content-box'; // override notebook setting of border-box.\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function (el) {\n",
       "    // this is important to make the div 'focusable\n",
       "    el.setAttribute('tabindex', 0);\n",
       "    // reach out to IPython and tell the keyboard manager to turn it's self\n",
       "    // off when our div gets focus\n",
       "\n",
       "    // location in version 3\n",
       "    if (IPython.notebook.keyboard_manager) {\n",
       "        IPython.notebook.keyboard_manager.register_events(el);\n",
       "    } else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function (event, _name) {\n",
       "    var manager = IPython.notebook.keyboard_manager;\n",
       "    if (!manager) {\n",
       "        manager = IPython.keyboard_manager;\n",
       "    }\n",
       "\n",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which === 13) {\n",
       "        this.canvas_div.blur();\n",
       "        // select the cell after this one\n",
       "        var index = IPython.notebook.find_cell_index(this.cell_info[0]);\n",
       "        IPython.notebook.select(index + 1);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_save = function (fig, _msg) {\n",
       "    fig.ondownload(fig, null);\n",
       "};\n",
       "\n",
       "mpl.find_output_cell = function (html_output) {\n",
       "    // Return the cell and output element which can be found *uniquely* in the notebook.\n",
       "    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n",
       "    // IPython event is triggered only after the cells have been serialised, which for\n",
       "    // our purposes (turning an active figure into a static one), is too late.\n",
       "    var cells = IPython.notebook.get_cells();\n",
       "    var ncells = cells.length;\n",
       "    for (var i = 0; i < ncells; i++) {\n",
       "        var cell = cells[i];\n",
       "        if (cell.cell_type === 'code') {\n",
       "            for (var j = 0; j < cell.output_area.outputs.length; j++) {\n",
       "                var data = cell.output_area.outputs[j];\n",
       "                if (data.data) {\n",
       "                    // IPython >= 3 moved mimebundle to data attribute of output\n",
       "                    data = data.data;\n",
       "                }\n",
       "                if (data['text/html'] === html_output) {\n",
       "                    return [cell, data, j];\n",
       "                }\n",
       "            }\n",
       "        }\n",
       "    }\n",
       "};\n",
       "\n",
       "// Register the function which deals with the matplotlib target/channel.\n",
       "// The kernel may be null if the page has been refreshed.\n",
       "if (IPython.notebook.kernel !== null) {\n",
       "    IPython.notebook.kernel.comm_manager.register_target(\n",
       "        'matplotlib',\n",
       "        mpl.mpl_figure_comm\n",
       "    );\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAyAAAAEsCAYAAAA7Ldc6AAAAAXNSR0IArs4c6QAAIABJREFUeF7tnXnUrXPhvj/HSTKXscE8r5SZhEKUTKVkaJIo8zyz8KUopDQRiaQyttCigVKiWMgyZBVC5iHKGBHOb32e3/ec73FwjvfuPdv17M/1rOWP8t77vZ/r3t53X+/ez95jxo0bN654SEACEpCABCQgAQlIQAISGACBMQrIACj7LSQgAQlIQAISkIAEJCCBjoAC4h1BAhKQgAQkIAEJSEACEhgYAQVkYKj9RhKQgAQkIAEJSEACEpCAAuJ9QAISkIAEJCABCUhAAhIYGAEFZGCo/UYSkIAEJCABCUhAAhKQgALifUACEpCABCQgAQlIQAISGBgBBWRgqP1GEpCABCQgAQlIQAISkIAC4n1AAhKQgAQkIAEJSEACEhgYAQVkYKj9RhKQgAQkIAEJSEACEpCAAuJ9QAISkIAEJCABCUhAAhIYGAEFZGCo/UYSkIAEJCABCUhAAhKQgALifUACEpCABCQgAQlIQAISGBgBBWRgqP1GEpCABCQgAQlIQAISkIAC4n1AAhKQgAQkIAEJSEACEhgYAQVkYKj9RhKQgASmLoHLL7+8XHTRRWW33XYrb3zjG6fuN3uZW3/00UfLPvvsU84999zy1FNPlZVWWql89atfLcstt9wUu1x11VXllFNOKVdeeWW54YYbynPPPVfGjRv3irmTTjqpHH300eVvf/tbmXfeecsuu+xSdt5555d8/b333lt23333jssLL7xQ1lxzzXLMMceUhRZaaIqd/AIJSEACEpg6BBSQqcPVW5WABCQwcAL1Afnee+/dPShfYIEFBvr964P797znPeX666/vOswxxxzluOOOK3fffXe55ppryqKLLjrZPoccckj50pe+VJZaaqnyxBNPlFtuueUVBeSEE04o2223Xdl4443LOuusUy677LLywx/+sBxxxBFl3333nfB9nnzyyU5+HnvssbLnnnuWaaedtpOPKjbXXXddmX322QfKyG8mAQlIQAL/n4AC4j1BAhKQwJAQeC0F5KyzziqbbbZZOfvss8vHPvaxjuhDDz1UFltssbLuuuuW0047bbKUH3zwwTLLLLOU6aefvuy0007l2GOPfVkBefrpp7tnPFZeeeVywQUXTLjNT33qU+W8887rhOdNb3pT9/8fddRRnZDUZ1dWXHHF7v+76aabyjve8Y7umZoqPB4SkIAEJDB4AgrI4Jn7HSUgAQmMOoH6DMKhhx76ktsd1LMhm266abn00kvLfffdV6aZZpoJPbbddtvyox/9qPzzn/8s00033as678kJyM9//vOy/vrrl5/97GdlvfXWm3B7V1xxRVlllVW6Z0KqjNSjvgSsHlVAJj7qsya33XZbufXWW19VH79IAhKQgARGl4ACMro8vTUJSEACrwmBet1EfQnS6aef3r3MqL4Eqh4f+chHyowzzviynep1GvWfKR1jx46d8KzCK31tfYlV/acKwsRHvVbjc5/7XHddxzvf+c4pfavu309OQA4//PBy4IEHlvqMyVxzzTXh9p599tnu2ZN6/Uu97qS+JGyGGWYoW221VfdSsImPgw46qBx22GHl8ccfLzPPPPOr6uQXSUACEpDA6BFQQEaPpbckAQlI4DUlMNKXYL3SsyaTnsT8889f7rjjjsme20wzzdS9BKsKx8TH+GcsfvnLX3bXa7yaY3ICUv/d8ccf312kPulRhWSttdbqJOzhhx8uc845Z/nCF75QqnBMfFQh2XHHHbuXYy2++OKvppJfIwEJSEACo0hAARlFmN6UBCQggdeSwEgF5Pbbby/1nykd9ZmFVVdddbJfVp8lqS+3mvTZht/85jedFNR3xtpoo42m9K26fz85Adl66607wXi5Z27mm2++7qLz8deC1P995JFHdtd7THycfPLJpd7OtddeW5ZZZplX1ckvkoAEJCCB0SOggIweS29JAhKQwGtKYKQCMpplfQZkNGl6WxKQgASGm4ACMtz7enYSkEBDBEYqIPVtaus/Uzrqsxv15UyTO7wGZEoU/fcSkIAEJDCegALifUECEpDAkBCoF1/vtdder/pzQEbzGpBNNtmk+zyOSd8Fa5tttik//vGPR+1dsOq7X22wwQYveRes+iGM9WVip556avn0pz/dLVrfenfMmDEveResD3zgA927YNV/PCQgAQlIYPAEFJDBM/c7SkACEpgqBOrF2dtvv/2rvrZhNK8BOfPMM8vmm2/+os8BqReC12dG6sXnZ5xxxoRzHv/Af+GFF35ZDlP6HJB55pmne8vd888/f0K+Ssc555zTfQ7IbLPN1v3/9fqP/fbbr1x99dVlhRVW6P6/m2++uSy55JKdqNV3DfOQgAQkIIHBE1BABs/c7ygBCUhgqhCoD7TrZ1/Uz8eoMlA/+XvDDTd8xbfhHc0Szz//fFlttdXKjTfe+KJPQr/rrrs6AZj43abGf0r7xO+sdeedd3af4VGP+gGDV155ZfniF7/Y/e/6Llzjn9Wo/3v8u1jVDzwc/0no9ZmP+ha9BxxwwITTqp+ovuyyy3afrF6Fo/L42te+VmrX+knoU3pZ2Wjy8bYkIAEJSOD/CCgg3hskIAEJDBGB+vkW9ZmQ+++/v/ssjEF9EGFF+Mgjj3TyUd+Fqn5ieX0JVL0uZfyzD+Mxv5yAXHLJJWXNNdd82SVWX331Uv/9xMeJJ57Yfd5HPb/6yej1WZNdd921e8nVxMc999xTdt9993LRRRd1PNZYY43uc1IWWWSRIVrdU5GABCTQLwIKSL/2sq0EJCABCUhAAhKQgAR6TUAB6fV8lpeABCQgAQlIQAISkEC/CCgg/drLthKQgAQkIAEJSEACEug1AQWk1/NZXgISkIAEJCABCUhAAv0ioID0ay/bSkACEpCABCQgAQlIoNcEFJBez2d5CUhAAhKQgAQkIAEJ9IuAAtKvvWwrAQlIQAISkIAEJCCBXhNQQHo9n+UlIAEJSEACEpCABCTQLwIKCGyv173+bbBGo1tnwVnfPLo3CLy1+d8wB7DV6FWad+xMo3djwFt6U5kW2Gp0K01XphndG4Td2jPlBVij0a/zSPnP6N8o6Bbvfv5JUJvRr3Lnvx8e/RuF3eLfHnsA1mh06zz37L2je4ON3ZoCAhtcAYENEtRRQAJooIgCAhojrKKAhOBAMQUENEZYRQEJwTUSU0BgQysgsEGCOgpIAA0UUUBAY4RVFJAQHCimgIDGCKsoICG4RmIKCGxoBQQ2SFBHAQmggSIKCGiMsIoCEoIDxRQQ0BhhFQUkBNdITAGBDa2AwAYJ6iggATRQRAEBjRFWUUBCcKCYAgIaI6yigITgGokpILChFRDYIEEdBSSABoooIKAxwioKSAgOFFNAQGOEVRSQEFwjMQUENrQCAhskqKOABNBAEQUENEZYRQEJwYFiCghojLCKAhKCaySmgMCGVkBggwR1FJAAGiiigIDGCKsoICE4UEwBAY0RVlFAQnCNxBQQ2NAKCGyQoI4CEkADRRQQ0BhhFQUkBAeKKSCgMcIqCkgIrpGYAgIbWgGBDRLUUUACaKCIAgIaI6yigITgQDEFBDRGWEUBCcE1ElNAYEMrILBBgjoKSAANFFFAQGOEVRSQEBwopoCAxgirKCAhuEZiCghsaAUENkhQRwEJoIEiCghojLCKAhKCA8UUENAYYRUFJATXSEwBgQ2tgMAGCeooIAE0UEQBAY0RVlFAQnCgmAICGiOsooCE4BqJKSCwoRUQ2CBBHQUkgAaKKCCgMcIqCkgIDhRTQEBjhFUUkBBcIzEFBDa0AgIbJKijgATQQBEFBDRGWEUBCcGBYgoIaIywigISgmskpoDAhlZAYIMEdRSQABooooCAxgirKCAhOFBMAQGNEVZRQEJwjcQUENjQCghskKCOAhJAA0UUENAYYRUFJAQHiikgoDHCKgpICK6RmAICG1oBgQ0S1FFAAmigiAICGiOsooCE4EAxBQQ0RlhFAQnBNRJTQGBDKyCwQYI6CkgADRRRQEBjhFUUkBAcKKaAgMYIqyggIbhGYgoIbGgFBDZIUEcBCaCBIgoIaIywigISggPFFBDQGGEVBSQE10hMAYENrYDABgnqKCABNFBEAQGNEVZRQEJwoJgCAhojrKKAhOAaiSkgsKEVENggQR0FJIAGiiggoDHCKgpICA4UU0BAY4RVFJAQXCMxBQQ2tAICGySoo4AE0EARBQQ0RlhFAQnBgWIKCGiMsIoCEoJrJKaAwIZWQGCDBHUUkAAaKKKAgMYIqyggIThQTAEBjRFWUUBCcI3EFBDY0AoIbJCgjgISQANFFBDQGGEVBSQEB4opIKAxwioKSAiukZgCAhtaAYENEtRRQAJooIgCAhojrKKAhOBAMQUENEZYRQEJwTUSU0BgQysgsEGCOgpIAA0UUUBAY4RVFJAQHCimgIDGCKsoICG4RmIKCGxoBQQ2SFBHAQmggSIKCGiMsIoCEoIDxRQQ0BhhFQUkBNdITAGBDa2AwAYJ6iggATRQRAEBjRFWUUBCcKCYAgIaI6yigITgGokpILChFRDYIEEdBSSABoooIKAxwioKSAgOFFNAQGOEVRSQEFwjMQUENrQCAhskqKOABNBAEQUENEZYRQEJwYFiCghojLCKAhKCaySmgMCGVkBggwR1FJAAGiiigIDGCKsoICE4UEwBAY0RVlFAQnCNxBQQ2NAKCGyQoI4CEkADRRQQ0BhhFQUkBAeKKSCgMcIqCkgIrpGYAgIbWgGBDRLUUUACaKCIAgIaI6yigITgQDEFBDRGWEUBCcE1ElNAYEMrILBBgjoKSAANFFFAQGOEVRSQEBwopoCAxgirKCAhuEZiCghsaAUENkhQRwEJoIEiCghojLCKAhKCA8UUENAYYRUFJATXSEwBgQ2tgMAGCeooIAE0UEQBAY0RVlFAQnCgmAICGiOsooCE4BqJKSCwoRUQ2CBBHQUkgAaKKCCgMcIqCkgIDhRTQEBjhFUUkBBcIzEFBDa0AgIbJKijgATQQBEFBDRGWEUBCcGBYgoIaIywigISgmskpoDAhlZAYIMEdRSQABooooCAxgirKCAhOFBMAQGNEVZRQEJwjcQUENjQCghskKCOAhJAA0UUENAYYRUFJAQHiikgoDHCKgpICK6RmAICG1oBgQ0S1FFAAmigiAICGiOsooCE4EAxBQQ0RlhFAQnBNRJTQGBDKyCwQYI6CkgADRRRQEBjhFUUkBAcKKaAgMYIqyggIbhGYgoIbGgFBDZIUEcBCaCBIgoIaIywigISggPFFBDQGGEVBSQE10hMAYENrYDABgnqKCABNFBEAQGNEVZRQEJwoJgCAhojrKKAhOAaiSkgsKEVENggQR0FJIAGiiggoDHCKgpICA4UU0BAY4RVFJAQXCMxBQQ2tAICGySoo4AE0EARBQQ0RlhFAQnBgWIKCGiMsIoCEoJrJKaAwIZWQGCDBHUUkAAaKKKAgMYIqyggIThQTAEBjRFWUUBCcI3EFBDY0AoIbJCgjgISQANFFBDQGGEVBSQEB4opIKAxwioKSAiukZgCAhtaAYENEtRRQAJooIgCAhojrKKAhOBAMQUENEZYRQEJwTUSU0BgQysgsEGCOgpIAA0UUUBAY4RVFJAQHCimgIDGCKsoICG4RmIKCGxoBQQ2SFBHAQmggSIKCGiMsIoCEoIDxRQQ0BhhFQUkBNdITAGBDa2AwAYJ6iggATRQRAEBjRFWUUBCcKCYAgIaI6yigITgGokpILChFRDYIEEdBSSABoooIKAxwioKSAgOFFNAQGOEVRSQEFwjMQUENrQCAhskqKOABNBAEQUENEZYRQEJwYFiCghojLCKAhKCaySmgMCGVkBggwR1FJAAGiiigIDGCKsoICE4UEwBAY0RVlFAQnCNxBQQ2NAKCGyQoI4CEkADRRQQ0BhhFQUkBAeKKSCgMcIqCkgIrpGYAgIbWgGBDRLUUUACaKCIAgIaI6yigITgQDEFBDRGWEUBCcE1ElNAYEMrILBBgjoKSAANFFFAQGOEVRSQEBwopoCAxgirKCAhuEZiCghsaAUENkhQRwEJoIEiCghojLCKAhKCA8UUENAYYRUFJATXSEwBgQ2tgMAGCeooIAE0UEQBAY0RVlFAQnCgmAICGiOsooCE4BqJKSCwoRUQ2CBBHQUkgAaKKCCgMcIqCkgIDhRTQEBjhFUUkBBcIzEFBDa0AgIbJKijgATQQBEFBDRGWEUBCcGBYgoIaIywigISgmskpoDAhlZAYIMEdRSQABooooCAxgirKCAhOFBMAQGNEVZRQEJwjcQUENjQCghskKCOAhJAA0UUENAYYRUFJAQHiikgoDHCKgpICK6RmAICG1oBgQ0S1FFAAmigiAICGiOsooCE4EAxBQQ0RlhFAQnBNRJTQGBDKyCwQYI6CkgADRRRQEBjhFUUkBAcKKaAgMYIqyggIbhGYgoIbGgFBDZIUEcBCaCBIgoIaIywigISggPFFBDQGGEVBSQE10hMAYENrYDABgnqKCABNFBEAQGNEVZRQEJwoJgCAhojrKKAhOAaiSkgsKEVENggQR0FJIAGiiggoDHCKgpICA4UU0BAY4RVFJAQXCMxBQQ2tAICGySoo4AE0EARBQQ0RlhFAQnBgWIKCGiMsIoCEoJrJKaAwIZWQGCDBHUUkAAaKKKAgMYIqyggIThQTAEBjRFWUUBCcI3EFBDY0AoIbJCgjgISQANFFBDQGGEVBSQEB4opIKAxwioKSAiukZgCAhtaAYENEtRRQAJooIgCAhojrKKAhOBAMQUENEZYRQEJwTUSU0BgQysgsEGCOgpIAA0UUUBAY4RVFJAQHCimgIDGCKsoICG4RmIKCGxoBQQ2SFBHAQmggSIKCGiMsIoCEoIDxRQQ0BhhFQUkBNdITAGBDa2AwAYJ6iggATRQRAEBjRFWUUBCcKCYAgIaI6yigITgGokpILChFRDYIEEdBSSABoooIKAxwioKSAgOFFNAQGOEVRSQEFwjMQUENrQCAhskqKOABNBAEQUENEZYRQEJwYFiCghojLCKAhKCaySmgMCGVkBggwR1FJAAGiiigIDGCKsoICE4UEwBAY0RVlFAQnCNxBQQ2NAKCGyQoI4CEkADRRQQ0BhhFQUkBAeKKSCgMcIqCkgIrpGYAgIbWgGBDRLUUUACaKCIAgIaI6yigITgQDEFBDRGWEUBCcE1ElNAYEMrILBBgjoKSAANFFFAQGOEVRSQEBwopoCAxgirKCAhuEZiCghsaAUENkhQRwEJoIEiCghojLCKAhKCA8UUENAYYRUFJATXSEwBgQ2tgMAGCeooIAE0UEQBAY0RVlFAQnCgmAICGiOsooCE4BqJKSCwoRUQ2CBBHQUkgAaKKCCgMcIqCkgIDhRTQEBjhFUUkBBcIzEFBDa0AgIbJKijgATQQBEFBDRGWEUBCcGBYgoIaIywigISgmskpoDAhlZAYIMEdRSQABooooCAxgirKCAhOFBMAQGNEVZRQEJwjcQUENjQCghskKCOAhJAA0UUENAYYRUFJAQHiikgoDHCKgpICK6RmAICG1oBgQ0S1FFAAmigiAICGiOsooCE4EAxBQQ0RlhFAQnBNRJTQGBDKyCwQYI6CkgADRRRQEBjhFUUkBAcKKaAgMYIqyggIbhGYgoIbGgFBDZIUEcBCaCBIgoIaIywigISggPFFBDQGGEVBSQE10hMAYENrYDABgnqKCABNFBEAQGNEVZRQEJwoJgCAhojrKKAhOAaiSkgsKEVENggQR0FJIAGiiggoDHCKgpICA4UU0BAY4RVFJAQXCMxBQQ2tAICGySoo4AE0EARBQQ0RlhFAQnBgWIKCGiMsIoCEoJrJKaAwIZWQGCDBHUUkAAaKKKAgMYIqyggIThQTAEBjRFWUUBCcI3EFBDY0AoIbJCgjgISQANFFBDQGGEVBSQEB4opIKAxwioKSAiukZgCAhtaAYENEtRRQAJooIgCAhojrKKAhOBAMQUENEZYRQEJwTUSU0BgQysgsEGCOgpIAA0UUUBAY4RVFJAQHCimgIDGCKsoICG4RmIKCGxoBQQ2SFBHAQmggSIKCGiMsIoCEoIDxRQQ0BhhFQUkBNdITAGBDa2AwAYJ6iggATRQRAEBjRFWUUBCcKCYAgIaI6yigITgGokpILChFRDYIEEdBSSABoooIKAxwioKSAgOFFNAQGOEVRSQEFwjMQUENrQCAhskqKOABNBAEQUENEZYRQEJwYFiCghojLCKAhKCaySmgMCGVkBggwR1FJAAGiiigIDGCKsoICE4UEwBAY0RVlFAQnCNxBQQ2NAKCGyQoI4CEkADRRQQ0BhhFQUkBAeKKSCgMcIqCkgIrpGYAgIbWgGBDRLUUUACaKCIAgIaI6yigITgQDEFBDRGWEUBCcE1ElNAYEMrILBBgjoKSAANFFFAQGOEVRSQEBwopoCAxgirKCAhuEZiCghsaAUENkhQRwEJoIEiCghojLCKAhKCA8UUENAYYRUFJATXSEwBgQ2tgMAGCeooIAE0UEQBAY0RVlFAQnCgmAICGiOsooCE4BqJKSCwoRUQ2CBBHQUkgAaKKCCgMcIqCkgIDhRTQEBjhFUUkBBcIzEFBDa0AgIbJKijgATQQBEFBDRGWEUBCcGBYgoIaIywigISgmskpoDAhlZAYIMEdRSQABooooCAxgirKCAhOFBMAQGNEVZRQEJwjcQUENjQCghskKCOAhJAA0UUENAYYRUFJAQHiikgoDHCKgpICK6RmAISDv3www+Xk08+uVxxxRXlgQce6G7lzW9+c1lllVXKlltuWeacc87olhWQCBsqpICg5hhxGQVkxMhwAQUEN8mICykgI0aGCygguElQhRSQYI6rr766rLPOOmWGGWYoa6+9dpl77rm7W3nwwQfLxRdfXJ566qly4YUXlhVWWGHEt66AjBgZLqCA4CYZUSEFZES4kF+sgCBnGVEpBWREuJBfrIAgZ8GUUkCCKVZeeeWy9NJLl+OPP76MGTPmRbcwbty4st1225Ubbrihe3ZkcsczzzxT6j8TH2+afYmX3GZQERtZcNY3Y7uNVjEFZLRIvja3o4C8NtxH87sqIKNJ87W5LQXkteE+mt9VARlNmsN3WwpIsOn0009frr322rLEEku8bPqmm24qyy67bHn66acne+uHHHJIOfTQQ1/0NWOmmalMM3aWoFU/IgpIP3aaXMt5x87U/5OYzBkoIP2fVwHp/4YKSP83VED6v+HUPAMFJKC74IILduKwxRZbvGz61FNPLQcffHC54447JnvrPgMSwO9BxGdAejCSAtLvkabQXgHp/7wKSP83VED6v+HUPAMFJKB77LHHlj333LNsu+22Za211nrJNSAnnnhiOfroo8sOO+ww4lv3GpARI8MFFBDcJCMq5DMgI8KF/GIFBDnLiEopICPChfxiBQQ5C6aUAhJOceaZZ5ZjjjmmXHPNNeX555/vbmXs2LFl+eWXL3vssUfZdNNNo1tWQCJsqJACgppjxGUUkBEjwwUUENwkIy6kgIwYGS6ggOAmQRVSQP7LOf7zn/+U+pa89ZhjjjnKtNNO+1/dogLyX+FDhBUQxAxxCQUkRocJKiCYKeIiCkiMDhNUQDBTIIsoILBZFBDYIEEdBSSABoooIKAxwioKSAgOFFNAQGOEVRSQEFwjMQUENrQCAhskqKOABNBAEQUENEZYRQEJwYFiCghojLCKAhKCaySmgMCGVkBggwR1FJAAGiiigIDGCKsoICE4UEwBAY0RVlFAQnCNxBQQ2NAKCGyQoI4CEkADRRQQ0BhhFQUkBAeKKSCgMcIqCkgIrpGYAgIbWgGBDRLUUUACaKCIAgIaI6yigITgQDEFBDRGWEUBCcE1ElNAYEMrILBBgjoKSAANFFFAQGOEVRSQEBwopoCAxgirKCAhuEZiCghsaAUENkhQRwEJoIEiCghojLCKAhKCA8UUENAYYRUFJATXSEwBgQ2tgMAGCeooIAE0UEQBAY0RVlFAQnCgmAICGiOsooCE4BqJKSCwoRUQ2CBBHQUkgAaKKCCgMcIqCkgIDhRTQEBjhFUUkBBcIzEFBDa0AgIbJKijgATQQBEFBDRGWEUBCcGBYgoIaIywigISgmskpoDAhlZAYIMEdRSQABooooCAxgirKCAhOFBMAQGNEVZRQEJwjcQUENjQCghskKCOAhJAA0UUENAYYRUFJAQHiikgoDHCKgpICK6RmAICG1oBgQ0S1FFAAmigiAICGiOsooCE4EAxBQQ0RlhFAQnBNRJTQGBDKyCwQYI6CkgADRRRQEBjhFUUkBAcKKaAgMYIqyggIbhGYgoIbGgFBDZIUEcBCaCBIgoIaIywigISggPFFBDQGGEVBSQE10hMAYENrYDABgnqKCABNFBEAQGNEVZRQEJwoJgCAhojrKKAhOAaiSkgsKEVENggQR0FJIAGiiggoDHCKgpICA4UU0BAY4RVFJAQXCMxBQQ2tAICGySoo4AE0EARBQQ0RlhFAQnBgWIKCGiMsIoCEoJrJKaAwIZWQGCDBHUUkAAaKKKAgMYIqyggIThQTAEBjRFWUUBCcI3EFBDY0AoIbJCgjgISQANFFBDQGGEVBSQEB4opIKAxwioKSAiukZgCAhtaAYENEtRRQAJooIgCAhojrKKAhOBAMQUENEZYRQEJwTUSU0BgQysgsEGCOgpIAA0UUUBAY4RVFJAQHCimgIDGCKsoICG4RmIKCGxoBQQ2SFBHAQmggSIKCGiMsIoCEoIDxRQQ0BhhFQUkBNdITAGBDa2AwAYJ6iggATRQRAEBjRFWUUBCcKCYAgIaI6yigITgGokpILChFRDYIEEdBSSABoooIKAxwioKSAgOFFNAQGOEVRSQEFwjMQUENrQCAhskqKOABNBAEQUENEZYRQEJwYFiCghojLCKAhKCaySmgMCGVkBggwR1FJAAGiiigIDGCKsoICH28jO7AAAYGklEQVQ4UEwBAY0RVlFAQnCNxBQQ2NAKCGyQoI4CEkADRRQQ0BhhFQUkBAeKKSCgMcIqCkgIrpGYAgIbWgGBDRLUUUACaKCIAgIaI6yigITgQDEFBDRGWEUBCcE1ElNAYEMrILBBgjoKSAANFFFAQGOEVRSQEBwopoCAxgirKCAhuEZiCghsaAUENkhQRwEJoIEiCghojLCKAhKCA8UUENAYYRUFJATXSEwBgQ2tgMAGCeooIAE0UEQBAY0RVlFAQnCgmAICGiOsooCE4BqJKSCwoRUQ2CBBHQUkgAaKKCCgMcIqCkgIDhRTQEBjhFUUkBBcIzEFBDa0AgIbJKijgATQQBEFBDRGWEUBCcGBYgoIaIywigISgmskpoDAhlZAYIMEdRSQABooooCAxgirKCAhOFBMAQGNEVZRQEJwjcQUENjQCghskKCOAhJAA0UUENAYYRUFJAQHiikgoDHCKgpICK6RmAICG1oBgQ0S1FFAAmigiAICGiOsooCE4EAxBQQ0RlhFAQnBNRJTQGBDKyCwQYI6CkgADRRRQEBjhFUUkBAcKKaAgMYIqyggIbhGYgoIbGgFBDZIUEcBCaCBIgoIaIywigISggPFFBDQGGEVBSQE10hMAYENrYDABgnqKCABNFBEAQGNEVZRQEJwoJgCAhojrKKAhOAaiSkgsKEVENggQR0FJIAGiiggoDHCKgpICA4UU0BAY4RVFJAQXCMxBQQ2tAICGySoo4AE0EARBQQ0RlhFAQnBgWIKCGiMsIoCEoJrJKaAwIZWQGCDBHUUkAAaKKKAgMYIqyggIThQTAEBjRFWUUBCcI3EFBDY0AoIbJCgjgISQANFFBDQGGEVBSQEB4opIKAxwioKSAiukZgCAhtaAYENEtRRQAJooIgCAhojrKKAhOBAMQUENEZYRQEJwTUSU0BgQysgsEGCOgpIAA0UUUBAY4RVFJAQHCimgIDGCKsoICG4RmIKCGxoBQQ2SFBHAQmggSIKCGiMsIoCEoIDxRQQ0BhhFQUkBNdITAGBDa2AwAYJ6iggATRQRAEBjRFWUUBCcKCYAgIaI6yigITgGokpILChFRDYIEEdBSSABoooIKAxwioKSAgOFFNAQGOEVRSQEFwjMQUENrQCAhskqKOABNBAEQUENEZYRQEJwYFiCghojLCKAhKCaySmgMCGVkBggwR1FJAAGiiigIDGCKsoICE4UEwBAY0RVlFAQnCNxBQQ2NAKCGyQoI4CEkADRRQQ0BhhFQUkBAeKKSCgMcIqCkgIrpGYAgIbWgGBDRLUUUACaKCIAgIaI6yigITgQDEFBDRGWEUBCcE1ElNAYEMrILBBgjoKSAANFFFAQGOEVRSQEBwopoCAxgirKCAhuEZiCghsaAUENkhQRwEJoIEiCghojLCKAhKCA8UUENAYYRUFJATXSEwBgQ2tgMAGCeooIAE0UEQBAY0RVlFAQnCgmAICGiOsooCE4BqJKSCwoRUQ2CBBHQUkgAaKKCCgMcIqCkgIDhRTQEBjhFUUkBBcIzEFBDa0AgIbJKijgATQQBEFBDRGWEUBCcGBYgoIaIywigISgmskpoDAhlZAYIMEdRSQABooooCAxgirKCAhOFBMAQGNEVZRQEJwjcQUENjQCghskKCOAhJAA0UUENAYYRUFJAQHiikgoDHCKgpICK6RmAICG1oBgQ0S1FFAAmigiAICGiOsooCE4EAxBQQ0RlhFAQnBNRJTQGBDKyCwQYI6CkgADRRRQEBjhFUUkBAcKKaAgMYIqyggIbhGYgoIbGgFBDZIUEcBCaCBIgoIaIywigISggPFFBDQGGEVBSQE10hMAYENrYDABgnqKCABNFBEAQGNEVZRQEJwoJgCAhojrKKAhOAaiSkgsKEVENggQR0FJIAGiiggoDHCKgpICA4UU0BAY4RVFJAQXCMxBQQ2tAICGySoo4AE0EARBQQ0RlhFAQnBgWIKCGiMsIoCEoJrJKaAwIZWQGCDBHUUkAAaKKKAgMYIqyggIThQTAEBjRFWUUBCcI3EFBDY0AoIbJCgjgISQANFFBDQGGEVBSQEB4opIKAxwioKSAiukZgCAhtaAYENEtRRQAJooIgCAhojrKKAhOBAMQUENEZYRQEJwTUSU0BgQysgsEGCOgpIAA0UUUBAY4RVFJAQHCimgIDGCKsoICG4RmIKCGxoBQQ2SFBHAQmggSIKCGiMsIoCEoIDxRQQ0BhhFQUkBNdITAGBDa2AwAYJ6iggATRQRAEBjRFWUUBCcKCYAgIaI6yigITgGokpILChFRDYIEEdBSSABoooIKAxwioKSAgOFFNAQGOEVRSQEFwjMQUENrQCAhskqKOABNBAEQUENEZYRQEJwYFiCghojLCKAhKCaySmgMCGVkBggwR1FJAAGiiigIDGCKsoICE4UEwBAY0RVlFAQnCNxBQQ2NAKCGyQoI4CEkADRRQQ0BhhFQUkBAeKKSCgMcIqCkgIrpGYAgIbWgGBDRLUUUACaKCIAgIaI6yigITgQDEFBDRGWEUBCcE1ElNAYEMrILBBgjoKSAANFFFAQGOEVRSQEBwopoCAxgirKCAhuEZiCghsaAUENkhQRwEJoIEiCghojLCKAhKCA8UUENAYYRUFJATXSEwBgQ2tgMAGCeooIAE0UEQBAY0RVlFAQnCgmAICGiOsooCE4BqJKSCwoRUQ2CBBHQUkgAaKKCCgMcIqCkgIDhRTQEBjhFUUkBBcIzEFBDa0AgIbJKijgATQQBEFBDRGWEUBCcGBYgoIaIywigISgmskpoDAhlZAYIMEdRSQABooooCAxgirKCAhOFBMAQGNEVZRQEJwjcQUENjQCghskKCOAhJAA0UUENAYYRUFJAQHiikgoDHCKgpICK6RmAICG1oBgQ0S1FFAAmigiAICGiOsooCE4EAxBQQ0RlhFAQnBNRJTQGBDKyCwQYI6CkgADRRRQEBjhFUUkBAcKKaAgMYIqyggIbhGYgoIbGgFBDZIUEcBCaCBIgoIaIywigISggPFFBDQGGEVBSQE10hMAYENrYDABgnqKCABNFBEAQGNEVZRQEJwoJgCAhojrKKAhOAaiSkgsKEVENggQR0FJIAGiiggoDHCKgpICA4UU0BAY4RVFJAQXCMxBQQ2tAICGySoo4AE0EARBQQ0RlhFAQnBgWIKCGiMsIoCEoJrJKaAwIZWQGCDBHUUkAAaKKKAgMYIqyggIThQTAEBjRFWUUBCcI3EFBDY0AoIbJCgjgISQANFFBDQGGEVBSQEB4opIKAxwioKSAiukZgCAhtaAYENEtRRQAJooIgCAhojrKKAhOBAMQUENEZYRQEJwTUSU0BgQysgsEGCOgpIAA0UUUBAY4RVFJAQHCimgIDGCKsoICG4RmIKCGxoBQQ2SFBHAQmggSIKCGiMsIoCEoIDxRQQ0BhhFQUkBNdITAGBDa2AwAYJ6iggATRQRAEBjRFWUUBCcKCYAgIaI6yigITgGokpILChFRDYIEEdBSSABoooIKAxwioKSAgOFFNAQGOEVRSQEFwjMQUENrQCAhskqKOABNBAEQUENEZYRQEJwYFiCghojLCKAhKCaySmgMCGVkBggwR1FJAAGiiigIDGCKsoICE4UEwBAY0RVlFAQnCNxBQQ2NAKCGyQoI4CEkADRRQQ0BhhFQUkBAeKKSCgMcIqCkgIrpGYAgIbWgGBDRLUUUACaKCIAgIaI6yigITgQDEFBDRGWEUBCcE1ElNAYEMrILBBgjoKSAANFFFAQGOEVRSQEBwopoCAxgirKCAhuEZiCghsaAUENkhQRwEJoIEiCghojLCKAhKCA8UUENAYYRUFJATXSEwBgQ2tgMAGCeooIAE0UEQBAY0RVlFAQnCgmAICGiOsooCE4BqJKSCwoRUQ2CBBHQUkgAaKKCCgMcIqCkgIDhRTQEBjhFUUkBBcIzEFBDa0AgIbJKijgATQQBEFBDRGWEUBCcGBYgoIaIywigISgmskpoDAhlZAYIMEdRSQABooooCAxgirKCAhOFBMAQGNEVZRQEJwjcQUENjQCghskKCOAhJAA0UUENAYYRUFJAQHiikgoDHCKgpICK6RmAICG1oBgQ0S1FFAAmigiAICGiOsooCE4EAxBQQ0RlhFAQnBNRJTQGBDKyCwQYI6CkgADRRRQEBjhFUUkBAcKKaAgMYIqyggIbhGYgpII0O/3Gk+88wz5ctf/nLZf//9y3TTTTeUJIb9HD2//t9t3bDfGw77fnWdYT9Hz6/f/w22cB/t/0IvPQMFZBhXfZXn9Pjjj5dZZ521PPbYY2WWWWZ5lal+fdmwn6Pn16/748u1dcN+bzjs+9V1hv0cPb9+/zfYwn20/wspIMO4YXxOw/5Dt4UfSsO+4bCfn/fR+McXJuh9FDNFXGTYNxz282vh52h85wYHfQYEPM7UruYPpalNeOrf/rBvOOzn18IvzmHfcNjPz/vo1P85PrW/g/fRqU3Y208IKCAJtSHJ+EOp/0MO+4bDfn4+uPO/wT4QGPb/Dj2/PtwLJ99x2Dfs/0K+BGsYN4zPadgvvKtghv0cPb/47o8JuiFmiqjIsO/nz9HoboEKeR9FzWGZ/yXgMyDeFSQgAQlIQAISkIAEJCCBgRFQQAaG2m8kAQlIQAISkIAEJCABCSgg3gckIAEJSEACEpCABCQggYERUEAGhtpvJAEJSEACEpCABCQgAQkoIN4HJCABCUhAAhKQgAQkIIGBEVBABoaa942OPfbY8pWvfKU88MADZemlly7f+ta3ykorrcQrGjS69NJLu3O75ppryv3331/OPffcstFGGwW3xIx8+ctfLuecc0656aabyvTTT19WWWWVcuSRR5bFF1+cWXiErb7zne+U+s8dd9zRJZdccsly8MEHl3XXXXeEt9SPLz/iiCPK/vvvX3bdddfy9a9/vR+lp9DykEMOKYceeuiLvqreP+t9dliOe++9t+y7777lF7/4RXnqqafKIossUr7//e+XFVZYofenuMACC5Q777zzJeexww47lPq7YxiO559/vtT76Y9+9KPu9+Bb3/rWsuWWW5YDDzywjBkzZhhOsTzxxBPloIMO6n4H/v3vfy/LLrts+cY3vlFWXHHFXp7flH63jxs3rvzP//xPOfHEE8ujjz5aVl111e53yaKLLtrL8x3m0grIMK87mXM788wzyxZbbFGOP/748q53vat70HP22WeXm2++ucw111y9p1IfEPzhD38oyy+/fPnoRz86dALywQ9+sGy++ebdL5HnnnuuHHDAAeXGG28sf/7zn8uMM87Y+/3OP//8Mnbs2O6XRv2F8oMf/KATymuvvbaTkWE6rr766rLpppuWWWaZpay55ppDJSA/+clPyq9//esJc73uda8rc8wxx1DM98gjj3QP5upm22+/fZlzzjnLX//617Lwwgt3//T9eOihh0p9gD7+qD9f3v/+95ff/va3ZY011uj76XX9v/SlL5Wvfe1r3c+X+nPlj3/8Y/nsZz9bDj/88LLLLrsMxTluttlm3e+G+iC8ClaVrWOOOab7XfG2t72td+c4pd/t9Q9x9Q90ddMFF1ywk68//elP3fm+4Q1v6N35DnNhBWSY153MuVXpqA9ev/3tb3df9cILL5R555237LzzzmW//fYbKir1L1nD9gzIpAPVBwtVHH/3u9+V9773vUO13/iTmW222ToJ2XrrrYfm/J588smy3HLLleOOO64cdthhZZlllhkqATnvvPPKddddNzR7TXwi9edk/SPHZZddNpTnN+lJ7bbbbuWCCy7oJGtYnh3YYIMNytxzz11OOumkCae78cYbd88q1wfqfT+efvrpMvPMM5ef/vSnZf31159wOvUPc/XZ5Pozp8/HpL/b6x+rqmTtueeeZa+99upO7bHHHus2PuWUU7o/2nlwCCggnC0G1uTZZ58tM8wwQ6l/nZz4ZUmf+cxnuqcs6w+rYTpaEJBbb721e7ag/qXnHe94xzDN1/0Vtj47V++f9RmQt7/97UNzfvWcqljVv0jWvyoPm4BUYZx11lm7vzy++93v7v4yOd988w3FfvV+uM4665R77rmnE//61+T68qTPf/7zQ3F+E59E/Z1RH9jtscce3bOtw3LUZ0C++93vlosuuqgstthi5frrry8f+MAHumdFPvnJT/b+NOvLr+ozq/VZyLXWWmvC+ay22mqlPht5ySWX9PocJ/3dfvvtt3fPPtbfE/Vn6fhj9dVX7/53femZB4eAAsLZYmBN7rvvvu6X5eWXX949KBh/7LPPPt0v0iuvvHJgXQbxjYZdQOqzVx/60Ic6efz9738/CKQD+R5Vpur989///neZaaaZymmnnVbWW2+9gXzvQXyTM844o3upR30JVn2APmwCUl8qUZ/hqdd91Ouw6vUg9ZqJ+nKQ+lfZvh/jX85RH5Rvsskm3Y71Gp76stYqlsN0nHXWWeUTn/hEueuuuzoRGZaj/uysQnXUUUd1L/msf+yo/03W67GG5ajXB77+9a/vfn7WZwJOP/307v5Zr1eqL7nu8zHp7/b6mKZe81Ef47zlLW+ZcGr1Ja71a+tLzz04BBQQzhYDa6KADAz1QL5Rff15fbBX5WOeeeYZyPccxDepf3WtD3jqU+j12brvfe97nSAPwzMgd999d3eh8q9+9auy1FJLdTiHTUAmvY9UQZ5//vm7vy4Pw8vo6oO6umF90DP+qNcNVBG54oorBvGfyMC+R32mp55vvTZrmI76R4C99967e2lnvQakvlywvtSs3keHRSJvu+22stVWW5V68XaVrPqSz/psT32Dlr/85S+9nlMB6fV8RQHp935Re1+CFWFDhnbaaafuJXP1l0u94G6Yj7XXXrt7ev2EE07o/WnWayM+8pGPdA8Ixh/1r6/1F+o000xTnnnmmRf9u96f8P+eQL3urO5YX4rV96PKVL0ou4rx+KNe6FtfV1+f6RmWo74T1kILLdS9696HP/zhYTmt7jzqdY/1Wp4dd9xxwnnV/er1H8P0bm315P71r3+Vxx9/vHtmoF6YXp+d/NnPftbrPX0JVq/nU0D6PV/evl6EXt9yt771bj3qU9H1tdn1Aa0XoedcB5WsF9vVNwyoF9fX1/G28BaD73vf+7r7aL2YsO9HfW32pG9xWt99Z4klluje1nXYruOpe9UHPHW/+ranw/AOQ/UlSfWZrIkvQt999927l7BO/KxI3++rda8q/fVc63UDw3TMPvvsnTDWZ5HHH1WO61sp33LLLcN0qhPOpb57W/1jVX3Z2TbbbNPrc3yli9DrBej1QvR6VOmqb9DiRei8qX0GhLfJQBrV10LWp5jrL5YqIvVteOvrfOtfferrRPt+1Ac79cLsetS3yqxPqde3y6wX/A7DRbD1Ytf6mt767MfEn/1RL/it7+DS96O+Bru+S0vdqj5Yr+da317xwgsv7P7qPIzHsL0Eqz4I2HDDDbuXXdWXfdb35q8vcalvh1nfsrbvR32pVX19fb22pb7G/KqrruouQK8XNQ/DBczj/zBVH6x+/OMfL/WzaobtqJ/5US/Qrr8H60uw6sXL9UF5fclS/XkzDEf9mVn/YFV/T9TfifUlZ/X6pSrO0047be9OcUq/2+tu9b468dvw3nDDDb4NL3BpBQQ4yqAq1bfgHf9BhPUdIr75zW92nwkyDEd9VqAKx6RHla5h+Av6K70NZv3LXf2l2vejXiNw8cUXdxcvV6mq10nUZwaGVT7qXsMmIPUtL+tLA//xj390wlHfeade4DsMn5Ex/r+v+ra0VZbrW9PWB+r1gvRhehes+u5Q9fqPerFyvW5g2I5JP6SvXmBfZat+6Gm95mUYjvqHxXofre/WVv8AV99muP53WH+u9vGY0u/28R9EWP8QUK87qz936tucD+P9t4/7TdxZAen7gvaXgAQkIAEJSEACEpBAjwgoID0ay6oSkIAEJCABCUhAAhLoOwEFpO8L2l8CEpCABCQgAQlIQAI9IqCA9Ggsq0pAAhKQgAQkIAEJSKDvBBSQvi9ofwlIQAISkIAEJCABCfSIgALSo7GsKgEJSEACEpCABCQggb4TUED6vqD9JSABCUhAAhKQgAQk0CMCCkiPxrKqBCQgAQlIQAISkIAE+k5AAen7gvaXgAQkIAEJSEACEpBAjwgoID0ay6oSkIAEJCABCUhAAhLoOwEFpO8L2l8CEpCABCQgAQlIQAI9IqCA9Ggsq0pAAhKQgAQkIAEJSKDvBP4fgwuaElLQZqIAAAAASUVORK5CYII=\" width=\"800\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "<matplotlib.animation.FuncAnimation at 0x7fb8ac51aa30>"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "animate_heatmap(u, t, sampling=2)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.9.9"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}