{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "57b20b67-4e25-4965-a52b-646444440043",
   "metadata": {
    "editable": true,
    "slideshow": {
     "slide_type": ""
    },
    "tags": []
   },
   "source": [
    "# Align two tracks \n",
    "\n",
    "Align two tracks with different coordinate systems with a geometric affine transformation.\n",
    "\n",
    "<div class=\"alert alert-block alert-warning\" style=\"padding:1em\">The function to map a track on a track is: <b><i>mapOn</i></b>.</div>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "00056bc1-5ccf-49fd-b634-13603ba2664d",
   "metadata": {},
   "source": [
    "## Let's start by importing tracklib library\n",
    "\n",
    "The first task is only useful for the online notebook and import the local tracklib code source. It’s not necessary if tracklib is installed from PyPI."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "a1d8611a-73a6-418e-b5fd-16331bfb9201",
   "metadata": {
    "editable": true,
    "slideshow": {
     "slide_type": ""
    },
    "tags": []
   },
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import os\n",
    "import sys\n",
    "\n",
    "#-------------------------------------------------------\n",
    "# 1. [if tracklib is not installed using pip] add tracklib's local path in the python path\n",
    "module_path = os.path.abspath(os.path.join('../../..'))\n",
    "if module_path not in sys.path:\n",
    "    sys.path.append(module_path)\n",
    "\n",
    "# 2. Import\n",
    "import tracklib as tkl"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "3c8d6092-5a12-4138-b938-feb8a6960dcf",
   "metadata": {},
   "source": [
    "The following two imports are necessary for the tutorial:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "e5cf671b-e7d4-484d-91a1-5ee1758f6e04",
   "metadata": {
    "editable": true,
    "slideshow": {
     "slide_type": ""
    },
    "tags": []
   },
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "cfe6964f-fdff-40a0-9c05-47f4469432f2",
   "metadata": {
    "editable": true,
    "slideshow": {
     "slide_type": ""
    },
    "tags": []
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Generated track from 05/02/2042 00:44:56 to 05/02/2042 01:44:56 [360 pts, 159.24m]\n",
      "360 360\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAD4CAYAAAAJmJb0AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAuO0lEQVR4nO3debzV4/bA8c9qTprnQSpkiAydokxRaECKSGaubpSKa869rrGUWblEkSH9EuGKVJKQcEJp0ogiFaLQdM5Zvz/WPveccuqc2nufZw/r/Xrt1977u4fv+rZPaz/7+T7PekRVcc45l5pKhA7AOedc/HiSd865FOZJ3jnnUpgneeecS2Ge5J1zLoWVCh1AfjVq1NBGjRqFDsM555LK7Nmzf1LVmgU9llBJvlGjRmRmZoYOwznnkoqIfLuzx7y7xjnnUljcW/Ii8g2wEcgGslQ1I977dM45Z4qru+YkVf2pmPblnHMuwrtrnHMuhRVHkldgsojMFpFeOz4oIr1EJFNEMtetW1cM4TjnXPoojiR/rKoeBXQE+ojICfkfVNURqpqhqhk1axY4Asg559weinuSV9UfItdrgQlAq3jv0znnnInriVcRqQCUUNWNkdunAnfGc5/OJStVWLcOvvnGrn/6CX7+GbZsgexsu5QtCxUrQqVKUK0aNGwI++5r950rSLxH19QGJohI7r7GqOqkOO/TuYS3aRN8+SV89pldFiyApUthw4Y9e79q1aB5czjySLscdxw0bhzTkF2SimuSV9XlwOHx3IdzySA7G2bPhilT7DJzJmzbZo/Vq2cJ+thjYf/9LTnXqgU1akD16tZ6L1UKSpaErVvti2DDBmvpf/uttfyXLYM5c+A//4HNm+19GzWCk0+GTp2gY0fYa69QR+9CkkRaGSojI0O9rIFLFdnZ8MEH8PLL8Oqr8OOPtv2II+CUU6BNG2jZEurXj90+s7Jg4UJ4/32YNg3eew9+/dUSfKdO0LMnnH46lC4du3268ERk9s4mmnqSdy7GliyBkSNh9GhL7OXLQ+fO0LUrtG9vrfTikpVlXzTjx+d90dSrB3/7m1322af4YnHx40neuTjLzoYJE2DYMGtFlyxpLeYLLrAWdIUKoSO0GN96C554At5+G0SgWzcYONB+Xbjktask7zNenYvC77/DY4/BAQdA9+6wciUMGmTXr71m2xIhwYN98ZxxBkycCMuXww03wOTJdqL2jDPA21epyZO8c3vgzz9h6FA7udmvH9Sta90hixfDzTfb/UTWqBEMHmwnbu+6y04Et2xpvzy+3WnRWpeMPMk7txs2b4ZHH4UmTeDGG6FVK0uQH31kfe4lS4aOcPdUqQK33QYrVsCtt9oX1YEH2rZNm0JH52LBk7xzRaAKL7xg3TL9+8PBB8OHH1ofd+vWoaOLXqVKcM899kvknHPsdvPmNkLHJTdP8s4VIjPTxrBfdBHUqQNTp1ryO/bY0JHF3j772JfZ1Kl2v107uPRS+OWXoGG5KHiSd24n1q61YYatWtmJylGj4JNPLPHZJO7U1a4dzJ1rXTgvvgiHH26jhlzy8STv3A5U4bnn4KCDbKz7P/5h3RiXXQYl0uh/TPny1m0za5bdPukkG26ZO1PXJYc0+pN1rnDffWcTly65xPrd5861UTTpXACsRQv4/HO4/HK49144/ngbIuqSgyd557DW+xNPQLNmMGOGjaCZMcMSvYO994ann4Zx46yYWosW3n2TLDzJu7S3Zo213q+6Co45BubNg2uuSb7hkMWhe3f49FOretmunU0ES6BJ864AnuRdWps4EQ47zAp5DRtmM0AbNQodVWI76CBL9J0720Sw/v2tZIJLTJ7kXVratMla66efbsMiMzOhT5/UHzUTK5UqWa2e666z1ny3bvDHH6GjcgXxJO/SzvLlVuZ32DBrhX76qfXFu91TogQ88IAl+TffhLZtrevLJRZP8i6tvPWWnTT85hv473/h4YehXLnQUSW3vn2tGNuCBXDCCT7yJtF4kndpITsbbr/d+pEbNbJVmk4/PXRUqeOMM+x8xo8/2hDLZctCR+RyeZJ3Ke/XXy2h33mnTdGfOdMKjLnYOvZYK/fw+++W6BcsCB2RA0/yLsUtXWoFxN5918bBjxplszddfLRoAdOn27DKtm1tKUIXlid5l7KmT4ejj4Z162zx7L//3UfPFIdDD7WJUiVK2HKH3nUTlid5l5KeesoWy65Tx0bPnHhi6IjSS9OmVslyyxabNPXdd6EjSl+e5F1KycmBm26CXr2sFen97+EceqidjP31V0v0q1eHjig9eZJ3KWPrVrj4YhgyBK6+2oZIVq4cOqr0dtRRMGmSjbpp3x7Wrw8dUfqJe5IXkQ4i8rWILBWRm+O9P5eeNm60ETQvvmjlcYcNg1KlQkflwOoB/fe/dhL8rLNsCUVXfOKa5EWkJDAc6AgcApwvIofEc58u/fz4o/W5T5sGzzxjC134CdbE0rat1eafMcPKOOfkhI4ofcS7rdMKWKqqywFEZCzQBfARtC4mFi+GDh1sFac337TbLjH16AGrVsENN9gyg/ffHzqi9BDvJF8fyD/JeRVwdP4niEgvoBdAw4YN4xyOSyVz5tgIGrAqki1bho3HFe4f/7CRNg88YIm+f//QEaW+ePfJF/Sjebvq06o6QlUzVDWjZs2acQ7HpYrMTFuOrmxZ+PBDT/DJQgQeegi6doVrr7VaQi6+4p3kVwH75LvfAPghzvt0Ke6jj2xIXpUq1sfbtGnoiNzuKFkSnn/eFgc//3xYtCh0RKkt3kn+M+AAEWksImWAHsAbcd6nS2HTpsGpp9okpxkzoHHj0BG5PVGhArz+uv0S69LFh1bGU1yTvKpmAX2Bd4CFwDhVnR/PfbrU9dZb0KkT7LefJfgGDUJH5KLRsCG8+iqsWGEt+qys0BGlpriPk1fVt1S1qarup6r3xHt/LjW9+aaNsW7WzE6y1q4dOiIXC8cdB48/Du+8Azf7LJq48OkiLuFNmgRnn219uFOmWF+8Sx1/+5uNlHrgAasYevbZoSNKLV7WwCW0qVPzWvCTJ3uCT1UPPGAVQy+7zOY+uNjxJO8S1vTpcOaZcOCB1oKvWjV0RC5eypSBcePs+pxz4M8/Q0eUOjzJu4T0wQe2VF/jxtaar149dEQu3ho2tNpD8+ZZgTnVwl/jCudJ3iWcjz+2UTT77GMrOvkcufRx2mnwr39ZnZuRI0NHkxo8ybuEkplp9Wfq1rUx8XXqhI7IFbd//tPKEvfr58sHxoIneZcwFi2Cjh2hWjVL8PXqhY7IhVCyJDz3nE2YOv98L00cLU/yLiF8950VGytZ0k6y+kSn9Fa3rpWNnjPHx89Hy5O8C27dOitVsHGjTYrZf//QEblEcPrpcM018MgjXsgsGp7kXVAbN1oXzbff2upBhx8eOiKXSIYMgebN4dJLfY3YPeVJ3gWzebNNdPryS3j5ZTj++NARuURTrhy89BL8/rsleh9Wufs8ybsgsrKgZ087wfrss/bT3LmCHHKIzYidPBmeeCJ0NMnHk7wrdqpw1VUwYYL1t154YeiIXKLr3dvG0F9/PSxZEjqa5OJJ3hW7u++Gp5+2Bbf79QsdjUsGIjY5qkwZWwg8Ozt0RMnDk7wrVqNH24zGiy6yZO9cUdWvD8OH24zooUNDR5M8PMm7YjN1qpWVPflka8lLQSsAO7cL558P3btbQ2HOnNDRJAdP8q5YzJ0L3brBQQfZakBlyoSOyCUjEVtkpFo1uPhi2Lo1dESJz5O8i7tVq6zgWKVKNqmlcmVsTNzXX8PataHDc0mmRg0YMcIaDoMHh44m8XmSd3H122+W4DdsgIkTrbIka9fCiSdas752bfjkk9BhuiRz5pnQo4ed15k3L3Q0ic2TvIubbdus/3Thghxe+cdMDv/qBRscX78+fP65TW+tUMGaZcuXF/1/67Zt8OGHdgbuq692vsKEKsyf/9cZNBs22LfPrkyaBA89BDk5RYvJFbtHH7VfhVdc4aNtdklVE+bSokULdakhJ0e1d29VUB1VLnIDVGvUUD36aLvdrp3q5Zerli6tWqKEavnyqsuWFf7mDz6Y935gr3355b8+77HH7PHRo7cPrE0b1dq1VZcvL/j9n3/e3hNUL7pIdevWPftHcHE3Zox9TA88EDqSsIBM3UleDZ7Y8188yaeORx6xv64be660G0OHqn7xhWp2tur69arHHaf60Ue2bf/9Va+9VrViRdWTT7ZEvCutW6seeqjq22+rjh1r7/+vf6nedJPqfvupvv66Pe+ii+yxs87Ke+377+d9MRxxxF/39cknqiKqJ51k7wmqXbsWHpMLIidH9YwzrH2wZEnoaMLxJO+K1VtvWQ7t0mmLZl/5d/szW7eu8Bc+8YQ9d8CA7bfn5Ki+9JLqwoWqq1fbc+66K+/x8uVVzz9ftVIl+1UAqi++qNqxo92uW1d1yxZ77r332rZHH7XradPsiycryx6/9VbVUqVUf/3V7t9xhz1v9uzo/2FcXKxaZR9927b2UaYjT/Ku2MybZw3yw2v9oBtLVLI/sU6divbi7GzdrhsmI0N1yhTVIUPsvkjeY5Mm5b2uXDnbVrGi6qefWndQw4aql1yS9/w2bVSbNFGtWVO1QgXVP/9UrV5d9dhjVdu3V23USPXJJ1VbtrTn5vrlF9WyZVWvuSam/04utp56yj7mJ54IHUkYQZI88G/ge+DLyKVTYa/xJJ/c1q5VbdxYtU6tbP2uVGNL7h9/vHvNq/fft6Sam9Rr1LDrbt1U7747L2m/917ea3K3TZtm959+Wv/XJdOypbXGS5dW7dBBtWpV1SOPtOeNHm2t9vxfLKA6cOD2MbVpY3F4l03Cysmxnr6KFVW/+y50NMUvZJK/fnde40k+eW3ebI3icuVUP+n/ov1pffXVnr1ZTo7qzz+rtmhh79Osmervv9tjrVrZtpkz856f28LPtWaNarVqqj16WP+/quq2bXb9xx95XTGqqrNmWfIfONC+oUB1wYLt47nwQtt+3312f/Vq+wJr0ED1nXf27BhdzC1bprrXXtZHn27fx57kXWzl5OQlTVXN+WW9Xtx0poKdB9VmzVSPOSb6/fTqZX+ic+bkbZs0yVrjGzbkbfvhh7xknmtPOmezs62Dt6Dt3bvbL4MHH7Qvg732stjOOGP39+PiZuhQ+1hefTV0JMUrZJL/BpgLjAKq7uR5vYBMILNhw4Zx/8dwUZo4UbVePRsdk52tet99OqjOwwqq/+Zfea3eIUOi39eGDaqLF0f/PrGwcaP13YMl+FmzVK+/3rp7VqywLqLNm0NHmfa2blU9/HDV+vVVf/stdDTFJ25JHpgKzCvg0gWoDZTEJlzdA4wq7P28JZ/gfvxRt+u7fvppfYPTVcjWHozRnPyP5Q5jTCVbtqjedpudDFa1rpr8x/z442Hjc6pq378iqv36hY6k+OwqyYs9Hl8i0gh4U1UP3dXzMjIyNDMzM+7xuN305JOwYIFNMQS4/XaYMoVFC5VW6ydxAEv48EOhfJsjrXpU376wbBk0aRI27njbvBnKl7fbdetCo0Ywc2bQkJzp2xf+8x+rmJGRETqa+BOR2apa4JHGLcmLSF1VXR25fS1wtKr22NVrPMkngJwcKLFDtYsdawJ/8w2/LVlLq1MqsZ6qZLbuR8OZY/Me//lnqF49/rEmgi++gFKlrPLazTfD0qWw336ho0p7v/0GBx9s372ffGIfUSrbVZKPZ+2aISLylYjMBU4Cro3jvly0Vq2yJXcqV4YPPij4OXPmwJw5ZDfYlwseymA5TRjPOTR88/Htn5cuCR7gyCPhsMPgggvsy/D++0NH5LA/40cesRJJw4aFjiawnfXjhLh4n3wg779vs0bLlFGtVUt1333zzlqtX2/9zYMH/+/pAwfapuHtXlEdMSJIyAlpwAD7h5k8OXQkTm0QWKdONvct1cfOs4s+ea9Cme62bYM+faBOHavv/tprsHIl9O9vj3/4oV0feSQA48fDPffYCk9XTekGV14ZJu5ENGiQddXccINXr0wAIrZcYE5Oeq8l7Ek+nT32mC3RNG+eLZrZqBG0bg033QTPPguLF8OYMVC1KrRty9y51qPTurX9BPbl+3ZQrpydlJ4zB0qWhI8+Ch1R2mvUCO64w9our70WOJhAPMmns7GRk6XXX29r8+XKHY4wcyb83//BpZfy88YynHWW9XW+8gqULVvs0SaHCy6Adu3s9rnnho3FATBgADRvDtdcAxs3ho6m+HmST1e//mrDDgYOtFZ8/mZ57sKZ118PNWuSdeu/6NEDvv8eJkywEQtuJ0qUgOeft9sXXxw2FgdA6dI2Cvj77+2HVrrxJJ9MtmyxSzRycuCNN+Dvf7fldE477a/Pad8eLrrInvvoo9w0qApTp8ITT8DRR0e3+7SQOyy5USMby5eVFTQcB8ccA7175424SSee5JPFDTdYn2+5ctGd1LvnHujSxc6gXnEFtGnz1+fUqAHPPQe//MLYnHN58EE7N3vZZXu+27Tyxx92/e230LgxnHJK3q8jF8y990KtWtCrV3otF+hJPlnMmZN3e8yYPXuPr76Cu+6yvuLff4enn7YThDuxYIGNomnTBh58cM92mZYqVbLrQYPsC3n6dPuWdEFVqWIt+dmzbdRNuvAknyyaNbPOxSZNLPOeeio8/LAl/O+/L/g1S5ZYhp40yboMLrvM/tKHD8+bjr8TGzbYudgKFWDcOBuE44qodm2480446yz7cr7hBvtC/eyz0JGlve7doWNHOxW1alXoaIrJzgbQh7j4ZKid2LZNtU4dW2t0+XLV/v23L4wFqhMmqH79dd5rVq9Wbdo07/GLL7brceMK3V1OjurZZ6uWLLn92hxuD23YYCtSlSiheu65tui4Lw4ezPLlNveva9fQkcQOvvxfkpsyxT6q8ePzth13nG277rrtk/0JJ1gFyEMPtZK4Tz2lesAB9tjZZxdpd/ffr/9be9vFyMyZthxh7dr2j9uokdXBd0Hcd599DK+9FjqS2PAkn+z69rWmxx9/5G1btcqWr9u4UXX4cNVnn7XsXK+efaxlyqhOnWrP/f571ZtusvX5CjF9urXgzz47/VbXKRZZWfZZgeqYMaGjSVtbt6oedpgt7pV//ZlktaskXyylhovKq1AWQNVGaDRvbkMfC7N5M7zwgr0md1JOEf3wAxx1lHXbf/pp3vlDF2O//mqziB94AK67LnQ0aevjj+HYY22yVLIPLAhVhdLFwrx5NhTvjDOK9vxy5ezE7G4m+G3bbNDNxo02o9UTfBxVrmyf0w8/hI4krbVubdNFUn3svCf5RDd1ql136hTX3dxwg5VaGTnSBvK4OBKBevU8ySeAQYOgZs3UHjvvST7RzZoFDRtC/fpx28XYsdaa6dcPeuxyWRcXM3XrwurVoaNIe+kwdt6TfKL75BObkx0n+Sc8DR0at924HXlLPmGcey506JC6Y+c9ySeyH3+0/vg4FYzJnfC0997w8ss+4alY1a4Na9eGjsJhvWePP27dNalYd96TfCL65huYPDnvbFDLljHfhapNgF261KoJ16sX8124XSlXzuvZJJDGja1C5YQJ8PrroaOJLU/yiWb0aPuLO+20vJK1cTgT+tBD8OqrMHgwnHhizN/eFaZMGasomkBDmNPdddfZcr19+ljx0FThST6RfPEFXHpp3v2xY61iZLVqMd3NzJm2+FPXrvCPf8T0rV1RVali/QMrV4aOxEWULm0lhlavhhtvDB1N7JQKHYDLZ/Bg+0sbM8YGrv/5Z8zr+/70E5x3ng3YGTXKl/ALpnZtuz77bC9clkBatbKGz9ChdkJ2N6ebJCSf8ZooliyBgw6yJsSgQXHZRU4OdO4M06bZbL+jjorLblxRbNxoM87KlrVZyi5hbNoEhx9u7ayvvrKBCYnOZ7wmg8GDrZ92wIC47eK++6zq8MMPe4IPrmJFK2uwZYv9vHIJo3x5+5X77bdw662ho4meJ/lEMGsWPPOMrU+W+zM+xt5/H267zSY79e4dl1243dWihV2PHBk2DvcXxx0HffvCY4/BjBmho4lOVEleRLqLyHwRyRGRjB0eu0VElorI1yJSwEKijpwcuOQSK6JRv74tNBEHa9bA+efD/vvDiBHeD58wTjjB+uT/+U876e4Syr332ho9F1+c3KNtom3JzwO6Adt914nIIUAPoBnQAXhcRHa+zly6euYZW0sVYNgw+wkfY9nZcMEFsH69TXiKwy7cnhKBJ5+0kTap0C+QYvbe2wq6rlqV3JOkokryqrpQVb8u4KEuwFhV3aKqK4ClQKto9pWS/vMfq46UO1QyDu66C959175DmjePyy5cNKpXt3V2q1QJHYkrQOvW1s353HO2DGYyileffH0g/wDgVZFtfyEivUQkU0Qy161bF6dwEsymTZZxZ8+2n+rnnReX3Uydaj1AF18Ml18el124WNh/f1i+3LpscquOuoQxcKBVFundOzlr2xSa5EVkqojMK+Cyq6ZnQb2+BY7VVNURqpqhqhk1a9YsatzJ7ZVXbGxW1apw4YVx2cUPP0DPnnDwwVaXw/vhE9gBB9gqLUcdBaecYudqXMIoXdq6bbZutQZTspUkLjTJq2p7VT20gMuuKjysAvbJd78B4CX3co0YYa23n36yRB9jWVl2ovWPP6wfvkKFmO/CxdKOX/QFTY7avBmuvtoKrGzaVDxxuf/Zf3/r8nzvPbjjjtDR7J54dde8AfQQkbIi0hg4APg0TvtKHllZ1kXzwQdw5ZVQIj7//P/6lw37evJJOOSQuOzCxdLJJ9tQjjp17Bv5mGPgqqtgyBBYuNCeM3asncO580677YrdpZdat+ddd8Hbb4eOZjfsbPHXolyArlirfQuwBngn32MDgWXA10DHorxfSi/k/cEHqvvvbws4ly6tumZNXHYzcaLt4sor4/L2Lt4ef9w+wNxLnTq20nSbNnnbTjwxdJRp688/VZs3V61WTXXFitDR5MEX8g7sxx9twG3uz+y2be13X4ytXAlHHAH77GNlC8qXj/kuXLypws03Wyv+tNOs5HSzZrbW74EH2oDtWrVgzpzQkaatJUus+neDBrZkZuXKoSPysgbh3X+/TV/PLR18ww0x38W2bTZIZ9s2G+rlCT5JicA999hadK+8YpfVq63kxYwZcM45Nt8+gRpn6eaAA+xj+fpr+zi2bQsd0a55ko+3n36yvtQLLrATbL/9FpdFuW+5xVrvTz8NTZvG/O1dcSpVyk6yVqhg9aAXLbJzObVq2Vi+336zk/cumHbt7COYOtVOn0T7nfv99/EbteNJPt5Gj7aSwbkFqitVivkuXn/dal316WPlUV2KqVEDDj3UbvfsaX1yub8KXTCXXWZj6EeOtP/ee5rot22znrk4TZfxevJxpWpf923a5P0njYMGDaB7d0v0LsWVKGErhy1eHDoSh420Wb/eemRFrNLr7s5JeeghmD8/bhXGPcnHRU4OvPiinTRbvDjudUlatEjeKdduD1SunNwVs1KIiI2fB1toZM0aa9eVLVu013/xhQ157tIFzjgjPjF6ko+HJ5+0PtWqVe3ntfehuFiqVAk2bAgdhYvITfS1a9tcteXLrY3XsOGuX/fdd3bKpWZNeOqp+MXnffLxMGGCDXf76Sf7tH2oi4ulypVtZanhw212ztatoSNKeyLWIn/pJfjySytNNXz4zkfefPaZVZr+9VdLF/Gs6OJJPtY2bIDp0+23V5xmtLo016SJne/p29fKVa9fHzoiF9Gjh01hOOoo+3j2289Ozk6aZOWJXnnFBtm1bm29ulOnQkaBo9tjx7NQrE2ebF/f8epgc65nTyumkmv58nCxuL9o0sTKe7/5ppUVGTQIOna00a/nnANvvGH16efMiX+CB++Tj73cUQ/F8em59FSmDIwZY2vTPf+8ZYvWrUNH5fIRgc6d7fLbb9aF8/vvNhr2iCOKfmI2JrF4WYMYGz7cfqetWWOTV5yLl5wcW3AEfAZsmttVWQNvycfKQw/BtGl2JgWsn9STvIunr77Ku7106fZdOM5FeJ98rIwcaSdclyyx32S+nJuLt7p1824fd5wN7XBuB57kY0UVOnSwipPr1tmgWefiqVYt+7v7/HMblN2zJ0ycGDoql2A8ycdKTo6vsefCOPJIq3lbt65NxHMuH0/ysZKd7ePiXTilS9vSRRMn5q0m5Rye5GNj+HDri2/UKHQkLp1dey1UrGjXPtrGRXiSj4WxY23w6513ho7EpbOaNa14yjvvwJQpoaNxCcKTfCz88gvUq2eTVJwLqU8f+1u8//7QkbgE4Uk+Ftq2tSIU69aFjsSluzJlLNFPmWIlNpYvh2++CR2VC8iTfCxs3mzXOTlh43AO4O9/t+vTTrMKWR07wtq1MGuW99WnIU/y0Zo3zyoB9u3rY+NdYqhe3WZgg61KtmiR/W22bg2nnmqFzF3a8CQfrVdesetbbgkbh3P5DRhgrfann95++4cfwt13BwnJheFJPloffwzNmlkpA+cSzUEHQatWdvv//s9a85mZ1jjxxUbSQlRJXkS6i8h8EckRkYx82xuJyCYR+TJyeSL6UBPQli3Wz+llXl2iEoFPPrFW/bnn2oLyX3xhhc2ffTZ0dK4YRNuSnwd0A2YU8NgyVT0icukd5X4S02uvWbHobt1CR+Jc0XTqBPXrWwG9994LHY0rBlEleVVdqKpfxyqYpDNiBOy7r53Mci4ZdOgAq1ZZsp8+3UfbpIF49sk3FpEvROR9ETl+Z08SkV4ikikimeuSZZy5KowaZfXjr7zSa9a45HPSSVYx9ev0baOli0Kzk4hMFZF5BVy67OJlq4GGqnokcB0wRkQqFfREVR2hqhmqmlEznkuWx9K4cXDFFdCyJfROzZ4ol+LatrXrKVPy5nm4lFRoklfV9qp6aAGX13fxmi2q+nPk9mxgGdA0dmEHlJMD99wDBx9sJ12rVw8dkXO7b7/9rN5Sv35QvjyMHx86IhcncelnEJGaIlIycrsJcACQGkvKv/mmLbt2yy3eTeOSlwhMmpT3N/zzz2HjcXET7RDKriKyCmgNTBSRdyIPnQDMFZE5wHigt6r+El2oCUDVWvGNG8P554eOxrno1K4Nc+bY7cqVw8bi4iaqhbxVdQIwoYDtrwCvRPPeCWn+fPj0U3jsMSjla6C7FJDbH1+hQtg4XNx4f8Pu+Owzu27fPmwczsXK77/b9V57WcmDlSvDxuNizpujRfXuu/DAA7D33tA0Nc4hO8fGjXad23CpWxd++CFcPC7mvCVfFGvWwFlnWavnqaf8hKtLHTsudLN6tY2fdynDs1VR/POfsGmTLcLQo0foaJyLnVNOsQlRQ4fmrVE8ZIjPhE0hnuQLM3q0td4HDPBuGpd6SpSwv+vrr4cVK6z+/EMPQf/+oSNzMeJJflfuvRcuvdRmBw4aFDoa5+Jv/Hg47zwYPhwWLw4djYsBT/I7M3UqDBxo4+HffhtKlw4dkXPxV7cuPPIIlCsHd90VOhoXA57kC7JhA1x2GRx4IIwcaX/wzqWL2rVtMfAxY2DBgtDRuCh5ki/Iww9bOdZnnrG6Hs6lmxtugEqV4Jpr/CRskvMkv6NffrHx8F27+opPLn3VrAl33GHltD//PHQ0Lgqe5Hd0//02QeSOO0JH4lxYucOFfQWppOZJPivLSgZv2gTr18Ojj9of92GHhY7MubCqVrXrzz7zLpsklr5JXtWGSNarZ90yuaNo/vjDxwg7l9+4cbZUoEtK6Zvkb7/dhki2bAlXXw2vvw5XXQW1akFGRujonAuvVCmbHAXwzju7fq5LWKIJ9DMsIyNDMzMz47+jnBwbNXPqqfDGG7atZ09YuNCGTLZoEf8YnEsGS5bkzfROoFzhticis1W1wNZpelahLFHCJn1UrGgr5AC89FLYmJxLRAccYNd9+oSNw+2x9OyuUbXKkvXrh47EucRXuzZs2xY6CreH0jPJr1ljK+LkVt1zzu3cgQfaoIStW0NH4vZAeib5b76xa0/yzhXu0kttxai5c0NH4vZAeib5FSvsunHjsHE4lwwWLbLru+8OG4fbI+mZ5FevtmtfvNi5wt10k13Png1btoSNxe229EzyzZvb9axZYeNwLhlUq2aroq1a5RVZk1B6JvkZM2wY5Yknho7EueSQu9A32KAFlzTSK8l//LGVLHj0UUvwdeqEjsi55CBi5Q0APvggbCxut0SV5EVkqIgsEpG5IjJBRKrke+wWEVkqIl+LyGlRRxqtCRNsivaTT8JJJ8Fjj4WOyLnk0rkzlC1rwyld0oi2JT8FOFRVmwOLgVsAROQQoAfQDOgAPC4iJaPc157LzrY6NQcfDOvWWcJv1ixYOM4lpb32sl/AnuSTSlRJXlUnq2pW5O4soEHkdhdgrKpuUdUVwFKgVTT7ispLL1ldmjvvtFIGzrk9c+yxNqTS++WTRiz75C8Hcr/i6wMr8z22KrLtL0Skl4hkikjmunXrYhhOPqNGwUEHQbdu8Xl/59LFl1/a9Y03Bg3DFV2hSV5EporIvAIuXfI9ZyCQBbyYu6mAtyqwhJ2qjlDVDFXNqFmz5p4cw8799JMt5/fRR9afWCK9zjM7F3PZ2XY9Y0bYOFyRFVqFUlXb7+pxEbkEOB1op3l1i1cB++R7WgPghz0Nco9kZ8Mhh9giIFu3bj8EzDm3Z6680spzeznupBHt6JoOwE3Amar6Z76H3gB6iEhZEWkMHAB8Gs2+dtvy5XaStX59K7B0/PHFunvnUtKpp9p1kyZh43BFFm09+WFAWWCKWF32WaraW1Xni8g4YAHWjdNHVbOj3NfumT/frl94AVqFO+frXEopU8Zmvf72W+hIXBFFleRVdf9dPHYPcE807x+VefPs+pBDgoXgXEqqXduWA7ztNqhUKXQ0rhCpeSZSFaZNs1LCe+8dOhrnUsvjj8OCBdYVmtuYcgkrNZP8U0/Be+9Bv36hI3Eu9XTqZEMof/8dDjssdDSuEKm3xmv37jB+vJUu6N8/dDTOpab8J15zcnx4cgJLvU/mjTfs+okn/A/PuXi5+GK48EK77bNfE1pqZcGcHNhnHzjhBGjaNHQ0zqWu0qXh6KPt9saNYWNxu5Q6SX7NGpuosWyZXTvn4it3hnrnzmHjcLuUOkm+b194/nno2RPOOSd0NM6lvu7d4frrbVnAb78NHY3bidRI8vPm2cnWG2+EF1/0JcqcKw4lSsAFF9htX0gkYaVGkq9WzVry114bOhLn0sthh0H16vDaa6EjcTuRGkm+Xj1b6al69dCROJdeSpaEK66wJP/dd6GjcQVIjSTvnAvn6qut6uu+++bVm3cJw5O8cy46++4Lw4fb7T59wsbi/sKTvHMueldfDeeeC0uX2nyVRYtCR+QiPMk752KjdWtYuxaGDYODD4bXXw8dkcOTvHMuVipXtuvcmlHvvhsuFvc/nuSdc7HRoQOcfTaccordX7s2bDwOSMUqlM65MOrWtUmJYFVg33rL6tpUrBg2rjTnLXnnXOz9+9+W4CtVgq1bQ0eT1jzJO+di78QT824/9VS4OJwneedcnD37bOgI0poneedcfMybZyNuMjNh/vzQ0aQtT/LOufho1gzef99uL1wYm/e85BJfEGg3+ega51z81Kpl1+vWRf9eWVnw3HPRv0+a8Za8cy5+ataEMmVgxYro3mf8eFvaE/LWlnVFElWSF5GhIrJIROaKyAQRqRLZ3khENonIl5HLEzGJ1jmXXEqVgowMePJJG2WjuvvvoQo33wxVqsCrr8LIkTEPM5VF25KfAhyqqs2BxcAt+R5bpqpHRC69o9yPcy5ZPf88HHUU9OoFb79d9NepWqGzs86ytZuvuQa6drVfBq7IokryqjpZVbMid2cBDaIPyTmXUpo0gUmTrFXfuTMsWZL32MSJMHdu3v2cHHjkEfjiCzjhBCt09sYb9ph30+yRWPbJXw7k/5puLCJfiMj7InL8zl4kIr1EJFNEMtfF4uSMcy7xlC0Ld99tif6882DLFpg+HU4/HTp2hM2b4f77LekPGGAt/w8/zHv9ypU2e9btNtFC+shEZCpQp4CHBqrq65HnDAQygG6qqiJSFthbVX8WkRbAa0AzVd2wq31lZGRoZmbmHhyGcy4pvPEGdOlifezPPgs//mjb+/WDRx/96/NV7SJSrGEmGxGZraoZBT1W6BBKVW1fyJtfApwOtNPIN4aqbgG2RG7PFpFlQFPAM7hz6ezMM+Ggg2DwYKhQIS+5F5Tgc6tZeoKPSrSjazoANwFnquqf+bbXFJGSkdtNgAOA5dHsyzmXIm6/HS67DObMgRtvzNs+d6612v/2N7v/8sth4ksxhXbX7PLFIkuBssDPkU2zVLW3iJwN3AlkAdnA7ar638Lez7trnEtDjz9u68R27mz3N22yyVMNG4aNK4nsqrsmqiQfa57knXNu9+0qyfuMV+ecS2Ge5J1zLoV5knfOuRTmSd4551KYJ3nnnEthnuSdcy6FeZJ3zrkU5kneOedSWEJNhhKRdcC3kbs1gJ8ChhMLqXAMkBrH4ceQGPwY4mNfVa1Z0AMJleTzE5HMnc3gShapcAyQGsfhx5AY/BiKn3fXOOdcCvMk75xzKSyRk/yI0AHEQCocA6TGcfgxJAY/hmKWsH3yzjnnopfILXnnnHNR8iTvnHMpLOGSvIgcISKzRORLEckUkVb5HrtFRJaKyNciclrIOAsjItdE4pwvIkPybU+aYwAQketFREWkRr5tSXEMIjJURBaJyFwRmSAiVfI9lhTHALbMZiTOpSJyc+h4ikJE9hGR90RkYeT/QP/I9moiMkVElkSuq4aOtTAiUlJEvhCRNyP3k+sYVDWhLsBkoGPkdidgeuT2IcAcbLnBxsAyoGToeHdyDCcBU4Gykfu1ku0YIvHuA7yDTVCrkWzHAJwKlIrcvg+4LwmPoWQkviZAmUjch4SOqwhx1wWOityuCCyO/LsPAW6ObL859zNJ5AtwHTAGeDNyP6mOIeFa8oAClSK3KwM/RG53Acaq6hZVXQEsBVoV8PpEcBUwWFW3AKjq2sj2ZDoGgIeAG7HPJFfSHIOqTlbVrMjdWUCDyO2kOQYsrqWqulxVtwJjsfgTmqquVtXPI7c3AguB+ljsoyNPGw2cFSTAIhKRBkBn4Ol8m5PqGBIxyQ8AhorISuB+4JbI9vrAynzPWxXZloiaAseLyCci8r6ItIxsT5pjEJEzge9Vdc4ODyXNMezgcuDtyO1kOoZkirVAItIIOBL4BKitqqvBvgiAWgFDK4qHsYZOTr5tSXUMpULsVESmAnUKeGgg0A64VlVfEZFzgZFAe0AKeH6w8Z+FHEMpoCpwDNASGCciTUiuY7gV6+74y8sK2JaQx6Cqr0eeMxDIAl7MfVkBz0/UscTJFOtfiMjewCvAAFXdIFLQ4SQmETkdWKuqs0WkbeBw9liQJK+q7Xf2mIg8B/SP3H2ZvJ9Jq7A+4lwNyOvKKXaFHMNVwKtqnXafikgOVtQoKY5BRA7D+qrnRP5TNgA+j5wET4pjyCUilwCnA+0inwck2DEUIpli3Y6IlMYS/Iuq+mpk8xoRqauqq0WkLrB25+8Q3LHAmSLSCSgHVBKRF0iuY0jIE68LgbaR2+2A2ZHbzdj+ZNlyEvdkWW/gzsjtptjPbUmmY9jheL4h78Rr0hwD0AFYANTcYXsyHUOpSHyNyTvx2ix0XEWIW4DngId32D6U7U9aDgkdaxGPpy15J16T6hiCtOQLcSXwiIiUAjYDvQBUdb6IjMP+02YBfVQ1O1yYuzQKGCUi84CtwCVqfxHJdAwFSrLPYRiWyKdEfpHMUtXeyXQMqpolIn2xUU4lgVGqOj9wWEVxLHAR8JWIfBnZdiswGOu+vAL4DugeJryoJNUxeFkD55xLYYk4usY551yMeJJ3zrkU5kneOedSmCd555xLYZ7knXMuhXmSd865FOZJ3jnnUtj/A87fQbAIqAxEAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "tkl.seed(123)\n",
    "\n",
    "reference = tkl.generate(0.2, dt=10)\n",
    "\n",
    "noised = tkl.noise(reference, 3, tkl.ExponentialKernel(100))\n",
    "noised = tkl.noise(noised, 20, tkl.GaussianKernel(1000))\n",
    "original = noised.copy()\n",
    "\n",
    "print (reference.size(), noised.size())\n",
    "\n",
    "plt.plot(noised['x'], noised['y'], 'r-')\n",
    "plt.plot(reference['x'], reference['y'], 'b-')\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "b9a61fc1-e981-49fa-a997-2126282710b6",
   "metadata": {
    "editable": true,
    "slideshow": {
     "slide_type": ""
    },
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAD4CAYAAADhNOGaAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABBp0lEQVR4nO3dd1xW1R/A8c8XBBy4xT1wb9PEUZZ7Z5qZZuXI8meWVrbVrEzLUhuWmmlqaqZlaY7c2zQXztxbQRzgRmWf3x8XcQGK8HDheb7v1+t5+XDu4HtV+D73nHPPV4wxKKWUcl1udgeglFLKXpoIlFLKxWkiUEopF6eJQCmlXJwmAqWUcnEZ7A7gQeTJk8f4+vraHYZSSqUrW7ZsCTHG+NzZni4Tga+vL/7+/naHoZRS6YqIHI+vXbuGlFLKxWkiUEopF6eJQCmlXJwmAqWUcnGaCJRSysVpIlBKKReniUAppVycJgKl0qFZe2ex6+wuu8NQTkITgVLpzKy9s2g3ox2PTHiENcfX2B2OcgIOTwQickxE/hOR7SJy1+PAYvleRA6JyE4RedjRMSmVFm0I3MCxi8cS3efQ+UN0m9ON6gWqUzhbYZpPbc6iQ4tSJ0DltFLrjqCBMaaqMcYvnm0tgNKxrx7AmFSKSalU5x/kz8RtE+9qvxJ+hYaTG1JnYh1OXj6Z4PGd/+qMu7gzs8NM1ry4hrJ5ytLmtzacuHTCkWErJ5cWuobaAFOMZQOQQ0QK2B2Ucm6BlwPZGLjxnp/A79f1yOsEXg7k/PXzJFb+tcZPNXh57st39e/P3jeb61HXCb4azDN/PBPvOc5fP8+GwA28++i7FMtRDJ8sPszsMJOI6Ah+2/VbilyHck2pkQgMsEREtohIj3i2FwICbvk6MLbtNiLSQ0T8RcQ/ODjYQaEqVzBh6wSKfFuE2hNqU3t8bSKjIx/4XC/PeZkcX+Yg85DMFPm2CLmH5ebluS/Hu290THTc+2VHlt22bda+WRTNXpTvmn/HhsANbDq56a7jt53aBkCNgjXi2krkLEHNQjX5fffvD3wNSqVGIqhjjHkYqwuol4jUvWO7xHPMXR+HjDHjjDF+xhg/H5+7VlFV6r4sOrSIcVvHAfDh4x9y5uoZVhxdgTGGy+GXk3SuK+FXmLxjMlXzV2VIwyGMbTWWCj4V2BuyF4CGkxtSfVx11p5YC8DZq2fjjt1+evtt5wq6EkT5POV5ocoLZPbIzMRtEzHGMGHrBE5dOQXA1lNbAahWoNptxz5b8Vm2ntqaYnc3yvU4fBlqY0xQ7J9nReQvoCZw61SHQKDILV8XBoIcHZdyLccvHmfwmsFM2DYhru2NWm8watMoxviPYeSmkSw9spT1L6/n4QL3N19hfeB6ok00/R/vT9OSTQFYfHgxB84dIOBSACuPrcTDzYN6k+oxp+McaheuHXfsymMruRJ+hWuR1wi4HMClsEsUy16MbF7ZaF+hPdN2TaOebz26z+sOQNX8VfFw86Bo9qLkyZzntjgye2QG4HToaXxz+Ma1B14OZNyWcRhjGNxw8AP9vSnX4NA7AhHJIiJZb7wHmgJ3Tn6eC3SJnT1UG7hkjDnlyLiUa1h9bDXPz3yeJ6Y9QamRpZi8YzLPV34+brtPZh+eLv80c/bPYdmRZWT2yEyvBb2IMTH3df41x9fgLu48WuTRuDZvT29CI0IZtWkUAGtfWku1/NXo/FdnvD298XL3wq+gH4GXA2kwuQElvy9JjZ9qsP/cfrJ5ZQOg/+P9yZUpFy/MeiHuvNtPb2dz0GaqF6h+VxwVfCqQ1TMr7Wa049y1cwCEXAuh7s91GbxmMJ/98xn7QvYl/S9QuQxHdw3lA9aKyA5gEzDfGLNIRHqKSM/YfRYAR4BDwE/Aaw6OSbmIr9Z/xex9szl+8Ti9avTiyBtHmNp2Klk9s/JokUcREd555B3aV2iPfw//uP75SdsnJXjOIxeOxPX1bwzcSGWfhwi77M3Jk3D8OLhFeXPs4jGG/TuM7tW6U7NQTT6u9zEXwy4yZccU3N3cqZa/GkMbD2XLqS00LdmUH1r+QJ0idWhUvBEAZXKXYdsr22hTtg3FcxTnz/Z/Mqj+IAC6PNTlrpjqFqvL6hdXE3w1mM5/dSY8KpznZz5P0JUg/mz/JwAz98xM4b9d5UwksRkOaZWfn5/RCmXqXgp/U5gGxRvwS9tfbmu/FnkNN3EjY4aMt7UbY6g7qS5rT6wl5L0QcmfOzYULsGsX7N8P8w/PZnbGtniGloKLvkQUXgbXc8DQCzdP0vgDeGwYnC9Brul7KFzAi3xFrrDMLxdGosjukYeFz67ikZIVORN6hnze+RK9BmMMItYw2tELRymes3iC+47bMo5X/n6F0rlKc/D8QX5o+QOv1niVRyY8QnhUOJv+t4kv/vmCDhU7UDZP2aT9ZSqnICJb4pvGny5LVSp1L2evnuXklZNUy1/trm03+tTvJCL0KN+ftSdakmd4HjKfbsy1q8C+p2BHV+j9Op4RZfGWfLjlvUgIQKaLjBwJnp7g7g4TQ87zbxjUy9WJ8u28CAqCkyezkuFEYyLz+HNp8goe/bAivr7w2GP5aNgQGjaEYsXiv44bSQBINAkA9KjeA093T16e+zKPFnmUV/xeAaBJiSYMXjOYV+a9wsTtE9kTsofp7abf+y9RuQxNBMop+QdZd4zxJYJbRUfDv//C/PmwcCHs/K8pDMgA7lFcy7+MTOTCq9xOOrx/nHG7A1n50joeLfIo0THRZBhs/fj07n3zfEdX5OPff2DMax0pf8vktoth0zlzBo41ysG2bbBlCyxeDFOnWtsrV4b27a1XuXIPft0vVn0Rv4J+FMpaCDexen79ClofACdun0h2r+z8feBvrkdeJ5NHpgf/RsqpaNeQcipnQs8wbN0wxviPwcPdgxN9TpA9Y/bb9jEGtm2DadPgt9/g5EnIkAEeewxatoQWLeBStn9ZdWwlGdwy0Hd5XzK4ZaBzlc5MbHPzqWD51Pq0bj65+TMUFhXGkQtHqOBT4Z6xGgO7d8PSpTBzJqxbZ7VXqgTdukHXrpA7d/L/TiKjIykzqgyFsxXmgzof8OT0J5n97GzalGuT/JOrdCWhriGMMenuVb16daPUnSKiIkyRb4oYt0/dTJe/upgDIQdu237pkjE//GDMQw8ZA8Z4eBjTurUx06db2+Iza88sw0BMxs8ymsBLgbdtyzAog3l70dspFn9goDHffWdM7dpWfF5exnTqZMy2bck/94XrF0xkdKSJiIowOb/MaTrN6pT8k6p0B/A38fxOtf2X+oO8NBGo+Kw8utIwEDP9v+m3tR89akyvXsZkyWL9j69a1UoI587d+5x7g/caBmLeW/KeY4JOwM6dVsxZs1oxt2xpzLp1KXPubrO7maxDspqQqyEpc0KVbiSUCNLCWkNKpYh5++fh6e5JqzKtAPjvP+jUCUqVgnHjrP73jRth61Z49VXIleve5yyXpxyLOy1mUINBDo7+dpUrw6hRcOIEfP45bNoEdepA/fqwYUPyzv1W7be4HnWdt5e8nSKxqvRPE4FyGvMOzKNh8YYc3e9N69ZQpQrMng19+sDRo/Dzz1CzJkh8i5okomnJpndNNU0tOXJA//5w7BiMGGFNY33kEXjuOavtQVTOV5m+dfoyZccU/jn+T8oFq9ItTQTKKewP2c/B8we5sOFJHnoI/vkHBg2yPlF/9RUUumsZw/QlSxZ48004eBA+/hjmzLFmF/XtC1evJv18/R/vj09mH75Y+0XKB6vSHU0EKl0KuBTAU789Rci1EC5dgte/XgXAthktePddOHwYPvro/rp/0hNvb/j0UzhwAJ59FoYOtWYZLUpibZpMHpl4o9YbLDy0kLcWvcWSw0scE7BKF3T6qErzomKiiIqJuq17pskvTVh2ZBmZ3LLhtqkPV0v8iuQ8xpFuEfgWc53PN//8Az16wL598MIL1rhCjhz3d+yV8CvUHF8zbh2iW6fBKueU0PRR1/mJUenS1Yir1PypJnUm1iE6Jpr5B+Yjn0rcev7XYy5ztfpgyHUYI9EulQQAHn8ctm+HgQPh99+tcZFVq+7v2KxeWfnr2b8cGJ1KL1zrp0alKxfDLuL9hTfbTm9j66mtTN05lXn7F8Rtdw98nAG5dnK5b6iNUdrPyws++cR6QjpTJmvJivffh/Dwex9bLk85CmYtSJuy+nCZK9MlJlSate6E9ajtk2WeJOhKEH2XfIz3/JnQ8AcAAgavpEB+dwB+ffrXBNcQchU1alhTY99+G4YPt+4MZs6EIkUSPy6LRxZdbsLF6R2BSpPeWfwOraZbzwP82GIClU5/yenrJwgo9iUAz1R4Ji4JADxf+XmeKveUHaGmKVmywNixVgLYtw+qV4eVKxM/JiI6Ak93T85ePUvX2V3vqqesnJ8mApWmfLbmM3xH+DJi44i4ttaNfZj8SWNyX69JeMmZeHt683Obn+0LMh14+mnrIbTcuaFxY2sKbULzQiJjIvFw86D3gt5M2TGFNr+14fz186kbsLKVoyuUFRGRlSKyV0R2i8ib8exTX0Quicj22NfHjoxJpW0frfyI45eOE2NiaJvrYzz/+YyjR61PuM1rlgasGr3ent42R5r2lStnJYO2beG99+CllyAi4u79sntl5489f/DHnj/oVKUTAZcC6PLX3QVwlPNy9B1BFPCOMaY8UBureH18yzL+Y4ypGvtK3Wf5VZpyo0pXrQvf8tcbn1Iz/EO2b7c+4RbMWhCATlU62Rhh+pI1K/zxhzWraNIka3XVS5du3+f1mq9zOfwyr9d8nclPTWZIoyHMPzifNcfXxHdK5YziW4DIUS9gDtDkjrb6wN9JOY8uOpe+7QveZ16c/aI5cv7IXdv8Rtc1mXs9bkSM+fBDYyIjb267En7FLDiwIBUjdS6TJhmTIYMxFSsac/z4zfaYmBiz8/ROExMTY4wx5mrEVZNveD7jN87PTNk+xRw+fzhum0rfSGDRuVR7oExEfIE1QCVjzOVb2usDM4FAIAh41xizO57jewA9AIoWLVr9+PHjjg9apbgFBxfw3MznuBx+mXrF6rGi64q4Aipz58JTq4vhdaoec7tNoUkTm4N1QsuXW3dX3t6wbBmULx//fj9v+5nu87oTY2IAKJ+nPD2q96BXjV54uHukYsQqJdn6QJmIeGP9su9zaxKItRUoZox5CBgJzI7vHMaYccYYP2OMn4+PT3y7qDRu4raJtJrWipI5S/Jp/U9ZfXw147eOJybG6rpo0/4KJmsgr3YsoUnAQRo1sgrgxMRA3brWw2jx6VatG9c/vM6OnjsY2WIk3p7evLX4LVYeu8cUJJUuOTwRiIgHVhL41Rgz687txpjLxpjQ2PcLAA8RyePouFTqCo8Kp9/yfjxW9DHWvrSWj+p+RAPfBry35D2adzjJp59Co5fWgMTwRKXH7Q7XqVWqBGvWWA+fNWiQ8LLWnu6eVMlXhd41ezOj/QwAAi8HpmKkKrU4etaQABOAvcaYbxLYJ3/sfohIzdiYzjkyLpX6ZuyewdmrZ/mo7kdk9siMiPBxtZ+4ci2SZaYv338Pldssx8vdi0eLPGp3uE6vdGlrnaLcuaFJE1i9OvH983vnB+DUlVOpEJ1KbY6+I6gDdAYa3jI9tKWI9BSRnrH7PAPsEpEdwPdAR5NaAxfK4S6GXaTLX13oMrsL5fOUp3GJxoBVIObZJiVxD3qECo8d4X+vhjF91zQal2isT7mmkmLFrDuDwoXhiSesJSoSkjFDRnJmzMmpUE0EzsihS0wYY9YCiZYBMcaMAkY5Mg5ln4UHF/LLzl8AGFB3ACLCrFnWSpkFC0Kx6tFkypyBX3b8wpmrZ3j7Ea2alZoKFoQVK6BePWjRwhpArlEj/n0LZC2gicBJ6ZPFyqE2ndxEpgyZuNz3Ms9Vep6vv4ZnnoGqVa2+aa9M0Zy8fJJh/w7j4QIP08C3gd0hu5wCBaxkkDs3NG2a8AByAe8CBF0JStXYVOrQRKAcalPQJqoXrE4Wj6y8/Ta8+y60a2f94vHxgaYlmhJ4OZBD5w/xQZ0PkKTWkVQponBh698ka1ZrSYp9++7ep0DWAjpG4KQ0ESiHiYyOZOuprfjlr0m3blbN3TfesNbNzxQ7DPBRvY8IeT+Ezf/bTPsK7W2N19X5+lrJwN0dmjeHoDs+/BfwtrqGdAjP+WgiUA4RER3BD5t/ICwqjDW/1WTKFKuG8IgR4HbH/zpvT2/8Cvrp3UAaUKoULFgA585ZYwa3LkdRwLsAEdERXAi7YF+AyiE0EagUN2ffHPJ9lY8+i/uQ43wjtk5vzejRVg1h/V2f9lWvDrNmwZ498NRTNwvc3JhCejr0tH3BKYfQRKBSVER0BG8seoN8mQtSauPfXBmzmGmTM/Haa3ZHppKiSRNrkbpVq6BbN2sJ6xsrvl6LvGZrbCrlaYUylaImbZ/EiUsnKL5uIUFrmjPnL2uOukp/XngBAgKgXz+oUAFqP28N7FyPvG5zZCqlaSJwYTEmhonbJnLu2jk+eOyDZJ8vIjqCwas+J9O5WgT904y5c63piCr9+uADq4voo49gcNGMAIRFhdkclUppmghcWOmRpTly4QgAjxR5hLrF6ibp+HPXzvFvwL+0KtMKEeG71ZMIDD2B14qxLJgvNGzoiKhVahKBcePg4EH4/OOs0A0uhV+694EqXdExAhfWrny7uPefrPokycd/sOwDWv/WmjcXvcnxwAj6L/oct6BaLB7dTJOAE8mYEWbPhlxeeQH4Y8ff9gakUpwmAhdWLk85ALpX686qY6toNrUZA1YMYPa+2UzYOoE1x9cQHhUe77FrT6xl0vZJAIzcNJIKXzQhKssJhrUcSL16OjXI2eTLB7OnW4sCzzgwmdVH1tsckUpJmghc2JLDSyjgXYDvW3xP7xq9OXXlFJ//8zltf29L93ndqTepHu8tfY8zoWcIuBQQd9zZq2d5ZsYzlMhZgsM9zpHzRGeu5V1D+ay1eLt1MxuvSDlSjYc9aJHrdQA+Gf2fzdGolKSJwEVFx0Sz9MhSmpRsQiaPTIxsOZIdPXfEbZ/YeiJZPLIwbss48n+dn2IjitF6emtWH1vNq/Nf5ULYBX554i+eeyoXob/+TPdiw/mz80R9KMzJ/d17BBlMJlbv3sfUqXZHo1KKDha7qG2nt3H++nmalrg5rUdEaFqyKUsOL6FjpY4sPbKUi2EXqe9bn9CIUMb4j6H+5PoADK47lPderMjWrTBzpjutW79rz4WoVOUmbpTPV5qTpffTowc89BBUrmx3VCq5NBG4qJVHrZKDjUo0uq19ZoeZnAk9QyaPTExrN+22bf0e68cvO3/hQPARVg15h3XrYNo0aN061cJWaUDlfJUIuLSQrPmCadfOB39/yJbN7qhUcmjXkItaF7COUrlKxS0bcIO3pzclc5WM95hMHpl46aEeHPjhS5Yvc2fiRHj22dSIVqUl/R7rx9XIUCq//yZHjtx88lilX6lRs7i5iOwXkUMi0jee7SIi38du3ykiDzs6JldnjOHfgH+pU6ROEo+DV16BefNg9Gjo2tVBAao0rVLeSgyoO4DlZ6fT+bO5zJoFY8bYHZVKDkfXLHYHRgMtgArAcyJS4Y7dWgClY189AP0v5WAHzx8k+FpwkhPBgAEwcSJ8/DG6dpCL6/tYX0rlKsUBn6E0bw7vvAO7dtkdlXpQjr4jqAkcMsYcMcZEAL8Bbe7Ypw0wxVg2ADlEpICD43Jp606sA6BO0ftPBCNHwpAh0KMHDBzooMBUuuHp7sm1yGv45vRl0iTInh2eew6u6zJE6ZKjE0EhIOCWrwNj25K6DyLSQ0T8RcQ/ODg4xQN1JRsCN5AzY864B8ruZcYMePNNaNPG6hLSGaIKrGQQER1BvnwwebJ1R/Dee3ZHpR6EoxNBfL8y7hxWup99MMaMM8b4GWP8fHx8UiQ4V7UreBdV8lXBTe79z79iBXTuDHXqwPTpkEHnmalYtQvXZkPgBs5fP893IS3p9u5BRo+GuXPtjkwllaMTQSBQ5JavCwN3Vr++n31UCjHGsCd4DxV87hyqudv27VZhktKlrR/uG+UllQKoVagWgZcDeWPhGyw8tBDz6FCqVYOXXoJT8ZQ23hi4kT92/8GF61rhLK1xdCLYDJQWkeIi4gl0BO78vDAX6BI7e6g2cMkYoxWyHeR06Gkuhl28ZyI4edKqI5A9OyxaBDlzplKAKt2onNd6kuzX/34FIOhqANOnw9Wr1ljSrVNK94fsp/aE2nT4swOjNo2yI1yVCIcmAmNMFNAbWAzsBWYYY3aLSE8R6Rm72wLgCHAI+AnQ+SgOtCd4D0CiiSA0FJ58Ei5fhvnzoXDh1IpOpSeNSjSKSwZgPaTolucgQ4YY/v4bJk26mQnm7r/5+W/g6oFEREekaqwqcQ5/jsAYs8AYU8YYU9IY83ls24/GmB9j3xtjTK/Y7ZWNMf6OjsmV7Q3ZC0D5POXj3R4dDc8/Dzt2wO+/Q5UqqRmdSm9WdF1Bx0od2fDyBiJjIikzqgxvX3bDbYA3L+/Li/9+6+Z+3oF5ccfEmBiCrmjvb1qiTxa7mH0h+8julf2uJ4pveOcd64GxkSOhZctUDk6lO3ky52F6u+nUKlyLT+rdrGmRP3tujNdF2o4YSMjVc6wLWEfZ3GXj9gm8HGhXyCoemghcQER0RNzaQscuHqN4zuLxrhI6ejR89x306aMPjKmk+6TeJ4x5wnoetHze0jTwfo3A/OOoNqIRMSaGKW2n0KFiB4DbljVX9tNE4AK+/vdrGk5pyNLDSzl4/iDFshe7a58FC+CNN6wF5L76yoYgVbonInR5qAtPlnmSr5p+xaK3v6Lk8UEEhu8mb6YC+BX0o0g2a4Lg4QuHbY5W3UpMOlwtys/Pz/j761DCvURGR/LMH8/EDdR5e3oTGhHK6Jajea3GzY/8u3dD7drWNNE1a8Db266IlbMJCIDyj++nfMVINs6rhJsbPDrhUQ6dP8SeXnvIkzmP3SG6FBHZYozxu7Nd7wic2OLDi+OSQMPiDQmNCGVQ/UG86vdq3D7nz1tPDHt7W2MDmgRUSipSBL7/uCz+CyoxYYLVNrrlaIKvBfP7rt/tDU7F0edEndjUnVPJnSk3u17bRd4sedkXsu+2aaNRUdYy0gEBsGoVFLprYQ+lkq9bN/jlF2v5iVatoHI+a8ppyLUQmyNTN+gdgZO6FHaJOfvn0LFSR/J758dN3O56duD992HZMvjxR3jkEZsCVU5PBMaNg7Awa82qDG4ZyOqZlYthF+0OTcXSROCkZu6dSVhUGJ2rdI53++TJ8O231gBxt26pHJxyOaVLW8uX//GH1QWZPWN2LoVfsjssFUsTgZOaunMqpXOVpmahmndt27jRKjDTsKHOEFKp5913oVIl6NULsnnk0DuCNEQTgRM6euEoq46tolOVTnc9LxAUBG3bQsGC1vLSHh42Balcjqen1UUUGAgXT2siSEs0ETih4f8OJ4NbBl6q9tJt7WFh8PTT1hpCc+ZA7tw2Bahc1iOPWA8rBh3JQdD5i2w7tY3a42tz8NxBu0NzaZoInEzQlSAmbptI14e6UjjbzdXijLFuyTduhClToHLlRE6ilAMNGQKZJBeHggMYtekHNp7cyE9bf7I7LJemicCJDFo9iELfFCI8OpwPHvvgtm3jx1v1hgcMsO4KlLJLtmzwdpPnifYKYd7O1YBVR1vZRxOBkzh+8TifrLIW9KqavyqlcpWK2+bvD717Q9OmWm9YpQ2DujQlx6XHCY6xEkBkdKTNEbk2TQROYvqu6QA0L9WcWR1mxbWfOwfPPAP588O0aeDubleESt3k5ib82WkSbgvHkDu8ug4c20yfLHYSfx/4m4cLPMzCFxbGtUVHQ6dOVtnAtWt1cFilLY0eLkHfRj0ZsmMH67x+tDscl+awOwIRGS4i+0Rkp4j8JSI5EtjvmIj8JyLbRURXknsAGwM3si5gHa1Kt7qtffBgq8zk999DjRo2BadUIvr3B2pYSWDJweX2BuPCHNk1tBSoZIypAhwA+iWybwNjTNX4VsVTiVt7Yi21J9QGoFWZm4lg4UIYNAi6drXqxyqVFmXJcvN9r98/ty8QF+ewRGCMWRJbsxhgA6CVbx3gj91/ADCyxUj8Clp59NgxeOEFq8zkDz9Ya70olVZ5uFlPNR6KXsnT057j5OWTNkfkelJrsPglYGEC2wywRES2iEiCn11FpIeI+IuIf3BwsEOCTG+MMcw9MJdWZVrRu2ZvRITISOjY0Rof+PNPyJzZ7iiVSlxo/1D8215FVn/C3AN/UW50Of4785/dYbmUZCUCEVkmIrviebW5ZZ8PgSjg1wROU8cY8zDQAuglInXj28kYM84Y42eM8fPx8UlO2E5j19ldHLt4jDZl4/66+fBD66Gx8eOhVKlEDlYqjfB096R6lcz0qjiQmFG7iYqO5kd/HTxOTclKBMaYxsaYSvG85gCISFegFfCCSaAUmjEmKPbPs8BfwN2rpKl43Sg6c2NsYOFCGD4cevaE9u3tjEyppBs4EHLElCTb6dbM2DNDny1IRY6cNdQc+ABobYy5lsA+WUQk6433QFNgl6NichZXI67y+67f+WvfX9QqVIv83vk5eRK6dLHGBb75xu4IlUq63Lnh00/h7NIXCLkWwux9s+0OyWU4coxgFJAVWBo7NfRHABEpKCILYvfJB6wVkR3AJmC+MWaRA2NK9y5cv0DTqU3pOLMjW05toXXZ1kRHW4PD167B779Dpkx2R6nUg+nZE8plaIHHxfJ8tPJjomKi7n2QSjaHPVBmjIm3hzq2K6hl7PsjwEOOisEZtf29Lf5BNx+3aFO2DYMHw+rVVrGZcuVsDE6pZPLwgBHfZKB5n8/Z3/FppuyYctcquirl6RIT6cius7tYfXw1XzT6Iq7tzK4KDBpkdQt16WJjcEqlkGbNoGXJp3A7VZOPV3yqdwWpQBNBOjJp+yQ83Dzo8lAXlnRawp9PLqdTJ6FMGRg92u7olEo533wt8E8/ToaeiJsUoRxHE0E6ERkdydSdU2lVphV5MuehUfEm/NS/IefPW5XGvL3tjlCplFO2LLze7Em44MuXK0bZHY7T00SQTiw/upwzV8/Q9aGugFVrePFiGDHCmimklLMZ+LE7mfZ1Z3PISr769ytGbhxJs6nNCI0ItTs0p6Orj6YTc/bNIYtHFpqVasb69dZiXe3bW0XolXJGOXLA+8268OnlAby39L249lXHVrH55GaCrgTxYd0P8c3ha1uMzkLvCNKBG0tJNCvVjOtXMtKxIxQtCj/9pOsIKefWv1cRsu7vCcAXDYfi5e7F3P1zGbRmEOO3jafMyDJM2DrB5ijTP00E6cDWU1sJuhLEk2Va07MnBAXBb79B9ux2R6aUY3l6wrjWY2CgoeDR96lVuFZcfeNeNXpRNk9Zxm4Za3OU6Z8mgjTu3LVzjNg4AjdxI3RbS2bMsJaXrqkLcSgX0aED+PlZ9bbL56oc196oeCPqF6vPvpB9xJgYrkdetzHK9E3HCNIoYwxvLX6LH/1/JDw6nA6lXqb/qz7UrQvvv293dEqlHjc3GDYMGjaE3Hv7MaRhIZqVaka1/NU4HXqaKxFXqPpjVTzdPfHvobWtHoQmgjRqT/Aevtv4He0rtKdfnY/o9Uxl3NxgyhStO6xcT4MG8MQTMPqLQhw+3C+u7Gp5n/IA/HfWWrb6dOhp8nvntyvMdEu7htKodQHrAPi84efMG1+Z9eutIjPFitkcmFI2+fJLuHIFhgy52VYtfzUq5a1En1p9AFh3Yp09waVzmgjSqHUB6/DJ7EPIgVIMGmQtKvf883ZHpZR9KlWCbt1g1Cg4etRqy54xO/+9+h9DmwwlY4aMcR+gVNJoIkij1p1YR60CdejcWShcWJeQUAqsZard3a2B41t5untSo2ANTQQPSBNBGnQm9AyHLxzmjH8djh6FX37RqaJKARQqBG+/DdOmwbZtt2+rU6QOW09t5VpkvOVPVCI0EaRBm05uAmDzX7Xp2xcef9zmgJRKQ957D3LmtMqy3qpO0TpExUSx+eRmewJLxxxZoWygiJyMLUqzXURaJrBfcxHZLyKHRKSvo+JJT86dt/4sXzojAwfaGopSaU727NC3r1Wa9Z9/brY/WuRRPNw86LO4D1V/rMrHKz+2L8h0xtF3BN8aY6rGvhbcuVFE3IHRWIXrKwDPiUgFB8eUphkDP36bE4B3BpzHw8PmgJRKg3r3hgIFrDW3blRDz5UpF9PbTWf76e3sOLODSdsn2RpjemJ311BN4JAx5ogxJgL4DWhjc0y2Gj8eNq7KBYB3ngs2R6NU2pQ5M3z0EaxdC4tuKW7brkI7xjwxBoDahWvbFF364+hE0FtEdorIRBHJGc/2QkDALV8HxrbdRUR6iIi/iPgHBwc7IlZbGWPoN+8rXvv7dfI9Y02UPn/9vM1RKZV2vfwylChhjRXExNxs7+nXk5I5S+Lp7mlfcOlMshKBiCwTkV3xvNoAY4CSQFXgFPB1fKeIp83E972MMeOMMX7GGD8fH5/khJ0mnb4SzJdb3yOq8kSiiy3FN4cvFXxcupdMqUR5elrTSbdtgz//vH1bRHSEJoIkSNYSE8aYxvezn4j8BPwdz6ZAoMgtXxcGgpITU3o1NnYBxc75hjPl9dfsDUapdOK552DoUKub6OmnIUPsbzRNBEnjyFlDBW75si2wK57dNgOlRaS4iHgCHQGXK1C6ezcM+dz6p6hRM+YeeyulbnB3h88+gwMHYPLkm+03EoExhjGbx3D26ln7gkwHHDlGMExE/hORnUAD4C0AESkoIgsAjDFRQG9gMbAXmGGM2e3AmNKciAjo3BmyZbX+KQyaCJRKitatoVYtq5soLMxqC40IZeSmkYzePJrXFrzGsHXD7A0yjXNYIjDGdDbGVDbGVDHGtDbGnIptDzLGtLxlvwXGmDLGmJLGmM8dFU9a9dlnVh/niG+t4ZIYo4lAqaQQsRaiCwiAH3+02iJjIgF4feHrAMzYPUN/thJh9/RRl7Zpk/UfuGtXeLKV9U+h/1mVSrqGDaFRI+vn6coVyO+dn3J5ygHgLu4EXA5gQ+AGm6NMuzQR2OT6dejSBQoWhO++g1l7ZwGQ3UsXFVLqQQwZAsHB1s9T4FuB7O21l4/qfsTYVmPxcvfih80/6AetBIgx8c7WTNP8/PyMv3/6rkT01lswYgQsWwY1HruMz3AfHi/6OAtfWIiHuz5OrNSDaNsWVqyAI0eIK14D0Pb3tszeN5vahWuz5sU1LvszJiJbjDF+d7brHYEN/vnH+tTSu7d1O3v26lkioiPo+lBXl/0PqlRKGDzY6hoadsfY8PfNv+edR95hQ+AGpu+abk9waZgmglR27Rq89BIUL25VXAK4FHYJgGxe2WyMTKn0r1Il6NQJvv8egm55IqlI9iIMbzKcSnkrMWzdMNJjT4gjaSJIZQMGwKFDMGECZMliteXzzocgbDm1xd7glHICAwdCVBR8fsccRBHh/UffZ3fwbhYeWmhLbGmVJoJUtG6dNS7QqxfUr3+zvXC2wrQs3ZLxW8fbFZpSTqNECejeHX76CY4du31bx0odKZKtCEPXDbUltrRKE0EquX7dqrdarNjNLqFbhVwLIXtGnTGkVEoYMADc3KyHzG7l4e7Bu4++y5rja1hyeIk9waVBmghSyUcfwcGDVpeQt/ft2zaf3MzGkxt5zU/XGFIqJRQqBK+9BlOmwL59t297pforlMhZgneWvKNjBbE0EaSCf/+Fb76Bnj2tB1/u9IP/D3h7etO1atfUD04pJ9W3L2TKBJ98cnu7VwYv+tbpy66zu9h1Nr4l0FyPJgIHu37dmiVUpMjdU9puWHF0BU+UfkJnDSmVgvLmhT59YMYM2L799m0tSrcAYPnR5akeV1qkicDBPvkE9u+3uoSyZr17+5XwK5y4dIIq+aqkfnBKObl334UcOeDjO8oXF8pq1b+6cF2rAIImAofasAG+/hp69IDGCVRu2BO8B4CKPhVTMTKlXEOOHPDeezBvnvXzeIOI4CZurDi2gpBrIbbFl1ZoInCQsDBrllChQjB8eML73eijrJhXE4FSjvDGG1Y30YABt7fHmBjWnljLF/98YU9gaYgmAgcZONCarTB+PGRLpOt/+dHl+GT2oUTOEqkWm1KuxNsb+vWD5cth5cqb7aVylQLgYthFewJLQzQROMCmTdZdQPfu0LRpwvtFRkey4OACWpVphZvoP4VSjtKzJxQubBW6vzFjdF8va17pxO0TbYwsbXBkqcrfRWR77OuYiGxPYL9jsZXMtotI+l5SFKvi2EsvQYEC8NVXie+75vgaLoVfok3ZNqkTnFIuKmNG61me9ethwQKrzd3N3d6g0hBHVih71hhT1RhTFZgJzEpk9wax+961PGp68+WXVg3iH3+E7Pd4UHju/rlkzJCRxiUSGElWSqWYbt2gZElrrCAmtixBfd/65MqUy97A0gCH90eIiAAdAKdf+3XPHqv0ZMeO0KpV4vvuD9nPjD0zaFyiMVk8s6ROgEq5MA8Pa+xu+3aYOdNqK5e7HO6idwap0TH9OHDGGHMwge0GWCIiW0SkR0InEZEeIuIvIv7BwcEOCTQ5oqOtMYGsWa1aA4k5HXqamuNrEhUTxYDHByS+s1IqxTz3HFSoYD1XEB1tLf1+KfySyy81kaxEICLLRGRXPK9bO72fI/G7gTrGmIeBFkAvEakb307GmHHGGD9jjJ+Pj09ywnaIH36w+h9HjLCmqiXm152/cjn8Miu7rqRW4VqpEp9SCtzdreI1+/bB1KlQwacCEdERLDuyzO7QbOXQUpUikgE4CVQ3xgTex/4DgVBjTKLDrGmtVOXx41CxIjz2GCxcCCKJ7//Qjw+RKUMmNnTXYtpKpTZjwM8Pzp+HnbvDKTG6MC1KtWBK2yl2h+ZwdpWqbAzsSygJiEgWEcl64z3QFEhXq0AZY01NAxg79t5JYMfpHew8s5POVTo7Pjil1F1ErLG8Y8fg1ylehFwL4Zedv9gdlq0cnQg6cke3kIgUFJHYCVzkA9aKyA5gEzDfGLPIwTGlqF9/hUWLYMgQq9bAvUzZMQUPNw+erfSs44NTSsWreXOoU8fqJrrBldcdcmgiMMa8aIz58Y62IGNMy9j3R4wxD8W+KhpjPo//TGlTcLC1umHt2lbVsXuJioli2q5pPFHmCfJkzuPw+JRS8ROxSlkGBUFjT2vCxoRtE2yOyj76OGsyvPkmXL5sLSPhfh8z0FYdW8Xp0NPaLaRUGlCvHjRpAtu/GUy9Ig0Zum4ol8Mv2x2WLTQRPKD582H6dOuR9Yr3uV7ctlPbAGjg28CBkSml7tdnn0FICJQ4OcAqF/ula5aL1UTwAC5ftgaIK1a0FrO6XwfOHcAnsw85M+V0XHBKqftWsya0aQMzv64f1xZjYuwLyCaaCB5Av35w8qTVJeTpeX/HGGPYE7KHsnnKOjY4pVSSDB4MVy4LTSNHAbA3eK/NEaU+TQRJtHat9fDYG29Yg8T341LYJVpOa8m/Af/yWJHHHBugUipJKle2loVZ83MzwFoM0tVoIkiCsDBrGYlixay+xfs1ctNIFh1axIhmI/isYRIOVEqlik8/hcizJfGOLsbiw4vtDifVaSJIgs8+s+oPjxtnFbu4H1ExUYzdMpYmJZrwZu03delbpdKg0qXh1Z7C1e3NWXp4ORHREXaHlKo0EdynnTth6FDo0iXxYjN3mrV3FoGXA+lV4z4eNFBK2ebjjyHjyaZciwqNm+HnKjQR3IeoKHj5ZciZE7755v6PuxZ5jfeXvk+lvJVoVeYe61IrpWzl4wMvPe0LwM8rV9kaS2rTRHAfvv8e/P1h5EjInfv+j+u/vD/HLx1ndMvR2iWkVDrwRnfrif+xR/rGFa9xBZoI7uHYMavEXatW0KHD/R83fN1wvtv4HW/UfIO6xeJdWVsplcYUznVz6ZcZM2wMJJVpIkiEMfDaa9a6JKNH33tl0RuOXjhK/xX9eabCM3zb/FvHBqmUSjGZPTLHve/XD8LDbQwmFWkiSMSMGVZ9gc8/h6JF7/+4QWsGkcEtA981/w430b9ipdKjY8esD4CuQH9LJeDCBWtROT8/6N37/o87cO4AU3ZM4VW/VymYtaDjAlRKOcTL1V6mgHcBmjWzpoxfcIHVqTURJOCDD6zFqMaNu7+VRW/4+t+v8XL3ou9jfR0XnFLKYbJ7ZedS+CWGDYOLF60eAWeniSAe//wDP/0Eb70F1ard/3GhEaFM2zWNDhU7kDfLPQoXK6XSpBwZc3At8hrlK0by4ovWbMGjR+2OyrGSW7y+vYjsFpEYEfG7Y1s/ETkkIvtFpFkCx+cSkaUicjD2T9uX5QwPh1desZaRGDgwacfO2D2D0IhQ/vfw/xwSm1LK8XJntuaIbzy5kcGDIUMGePddm4NysOTeEewCngZuW6VJRCpglamsCDQHfhCR+DpY+gLLjTGlgeWxX9tq6FDYuxfGjIEsWZJ27Pit4ymfpzyPFnnUMcEppRyuQ8UOlMhZgra/t+V6pkN8+CHMmgVLltgdmeMkKxEYY/YaY/bHs6kN8JsxJtwYcxQ4BNRMYL/Jse8nA08lJ57k2r/f6g/s2BFatEjasbvP7mZ94Hq6P9wdud95pkqpNCdP5jwsemERxhhKjyxNtoY/Uro0vP66804nddQYQSEg4JavA2Pb7pTPGHMKIPbPBDvWRaSHiPiLiH9wcHCKBmt9f6tLKHNmGDEi6ceP3zoeDzcPLUOplBMonbs0LUpbnwZfX/wq338PBw7At076WNA9E4GILBORXfG82iR2WDxt5sHDBGPMOGOMnzHGz8fHJzmnitfPP8Pq1TB8OOTLd+/9I6MjWR+wPq6a0fyD82leqjk+WVI+NqVU6gu5FhL3vlGTSJ56yipiExCQ8DHp1T0TgTGmsTGmUjyvOYkcFggUueXrwkBQPPudEZECALF/nk1K8Cnl7FlrMOjxx+GllxLf1xjDrL2zqDSmEo9OfJSha4dy/vp5Dp4/qGMDSjmRMU+M4eECDwNw/vp5vv0WYmKcc+DYUV1Dc4GOIuIlIsWB0sCmBPbrGvu+K5BYcnGYt96C0FAYOxbc7vE38s36b2g3ox3u4k4D3wZ8suoTXl/4OgCPFH4kFaJVSqUG3xy+9KnVB7Cmhvv6Qv/+1ooDi52sdk1yp4+2FZFA4BFgvogsBjDG7AZmAHuARUAvY0x07DHjb5lq+iXQREQOAk1iv05VixfDtGnWP3D58onvu+3UNvot78dT5Z5i56s7+bPDn+TJnIdp/02jZ/WeuricUk7G29OqQBUaEQrA++9DuXLQowdcuWJnZClLjElW170t/Pz8jL+/f7LPc+0aVKpkFaDfsQO8vBLfv+HkhuwN2cuuV3fFzTXedmobm4M287+H/6ezhZRyMksPL6Xp1KZ0q9qNsa3G4uHuwfr1UKeOtSDlqFF2R5g0IrLFGON3Z7tLP1n86afWE4Njx947CYC1MmGOjDnikgBAtQLV6FG9hyYBpZxQzUI1aVO2DT9v/5l5B+YB8Mgj8MYb1oJ0//xjc4ApxGUTwY4d8PXXVuWxevXu75ii2Yty9qot49lKKRtkz5id35/5HU93TzYGboxr//xz8PW1fn9cv25ffCnFJRNBdDT8739WtbFhw+7/uLNXz5Ivy33MLVVKOQ2vDF5Uy1+NDSc3xLVlyQLjx8PBg9DX9vUQks8lE8EPP8DmzdaDY7ly3f9xZ66eIZ+3JgKlXE3twrXZfHIzYVFhcW2NGllL1X//PSxYYGNwKcDlEkFAgDVDqFkzaymJpDgTekbvCJRyQa3KtOJ61HVm75t9W/uXX0KVKtCtG5w5Y09sKcHlEsHrr1tdQ2PG3H/pyRtOh54mv3d+xwSmlEqzGhZviG8OX8ZvHX9be8aMMH06XL4MnTtbv1tSWlRUyp/zTi6VCP76C+bMsWYLFS+etGOvR17nSsQVvSNQygW5iRsvV3uZ5UeX8+T0J7kWeS1uW4UKVs2CpUutOscpadMmKFECtm1L2fPeyaUSgb8/VK0Kffok/dgzV637Ph0jUMo1dX+4OwB/H/ibj1Z8dPu27tZzBcOHWw+oppRPP7WedypVKuXOGR+XSgSffw7//gseHkk/9sa00VyZkjC6rJRyGvm98xPaz3rC+PCFw3dtHzEC6ta1ppRu2HDX5iRbv94ahH7nHciaNfnnS4xLJQKATJke7LiSOUuSxSMLc/fPTdmAlFLpRhbPLPSs3pNlR5YRHhVO8NVgAi5Zy5F6eMAff0DhwtCyJeza9eDfJzoaevWCggWhd+8UCj4RLpcIHlSOjDnwdPfEJG81baVUOlenaB2uRl7l2MVjlBlVhqIjisZty5vXqmSWMSM0bWpVO3wQo0ZZ4wLffuv4uwHQRHDfdpzZwYWwCzQq3sjuUJRSNsqRMQcAO8/s5GLYRQAOn7/ZVVS8uJUMYmKspe03b07a+bdssRa3a9kS2rdPoaDvQRPBfVp+ZDlgTSNTSrmunBlzAtB1dte4tiMXjty2T6VKsG4dZMtmLWEzadL9nfvwYXjySas41uTJSZ/i/qA0EdxDyLUQ3l78NgNWDqBa/moUzFrQ7pCUUjbyK+hHv8f60aJ0Cya1mQRA8LW7y+eWLGkN+NaubT1w9uyzcPx4wudduxYeewwiImDhQsiTx0EXEA9NBIlYfGgxJb8vyXcbv+OFyi8w77l5doeklLKZVwYvhjQawswOM2lVphUAa46viXfffPms5wsGD4a5c6FsWWuq6aJFcPIknD5tbe/UyZpx5O1tlcytWDE1ryiZ9QhEpD0wECgP1DTG+Me23ygy4wlEAO8ZY1bEc/xA4H/AjXTa3xhzz1U7Uqoewb3U/Kkm56+fZ+5zc6ngU8Hh308plb4YY2g1vRXLjyzn4OsHKZK9SIL7njhhJYTp0+Hq1du3eXvDq6/CgAFWd5KjJFSPILmJoDwQA4wF3r0lEVQDzhhjgkSkErDYGFMonuMHAqHGmK+S8n1TIxHsOL2DqmOrMqLZCN6s/aZDv5dSKv06dvEYxb+zlio4/c7pez50evUqbNwIBw5Yy0eULGmNI2TO7PhYE0oEGZJzUmPM3tiT39l+6wPRu4GMIuJljAlPzvdLTeO3jsfT3ZNOVTrZHYpSKg3zzeEb9/7r9V8zrEnia9tnyQING1qvpDLGOKQIVmqMEbQDtiWSBHqLyE4RmSgiORM6iYj0EBF/EfEPDr57YCYlhUWFMfW/qbQr3+62amRKKZWYfSH7HHbujYEbcRvkluB4RHLcMxGIyDIR2RXPq819HFsRGAq8ksAuY4CSQFXgFPB1QucyxowzxvgZY/x8fHzu9a2TZdPJTVwMu0jHSklcp1op5ZK+bfYtAPMOzHNYMth0chMAh84fSvFz3zMRGGMaG2MqxfOak9hxIlIY+AvoYoy5e2EO69xnjDHRxpgY4Ceg5oNcRErbErQFsOqVKqXUvfSp3YdN3Tfh6e7JmM1jUvz8MSYmbuHLOkXqpPj5kzVGkBARyQHMB/oZY9Ylsl8BY8yp2C/bAslYnSPl+J/yp1DWQlp7QCl132oUqkGJnCUICg1K0fPuC9lHtznd2BC4gfq+9SmZq2SKnh+SmQhEpC0wEvAB5ovIdmNMM6A3UAr4SERurNfa1BhzVkTGAz/GzjAaJiJVAQMcI+EupFTx7J/Pcu7aOfYE79G7AaVUkvlk9iH4asqOYfZZ1Id9Ifv4uc3PdHmoC26S8kO7yZ019BdW98+d7Z8BnyVwTPdb3ndOzvdPSVExUczaOyvufe3CtW2OSCmV3uTNkpfdwbtT7HzRMdH8G/Avnap04sWqL6bYee/kkK6h9CjoShBRMVGMeWIM1fJXo0q+KnaHpJRKZ3xz+PL3gb+Jiokig1vyf71+t/E7rkRccci4wK10iYlYRy8cBay6A7UK1yKTxwMWLlBKuazKeSsTHh3OyI0jiTExD3yeTSc30X95f95Z8g4ALUu3TKkQ46V3BLGOXTwGQPGcSSxmrJRSsZ4u/zRTdk7h7SVvM2f/HBZ1WkTGDBmTdI41x9dQb1K9uK//bP8nOTMl+IhVitA7gljHLh5DEIpkS3itEKWUSkxWr6ws67yMoY2Hsvr4atYHrL9rn7CosLvuFk6HnqbU96X4ZccvNJx885Hj+r71aVehncPj1kQQ69CFQxTKVgivDF52h6KUSsdEJG5gt+GUhuw6e3NW/LXIa+QZloc3F1rrl4VFhREaEcqCgws4fOEwXWZ3IdpEs7PnTswnhhVd7lqr0yE0EWANFM/cM5MGvg3sDkUp5QTyZsnLa36vAdBtTjcioyMB+Hb9t1yNvMqozaMwxtDi1xZk/SIr7y5597bjK+erDNy9jpujaCIABq0eRFRMFAPrD7Q7FKWUkxj9xGj+aP8H/kH+fLbmM4KuBPHF2i/itm87vY1Vx1YBcCHsgk1RWlw6EcSYGL7b8B3jt47nleqvUCJnCbtDUko5kWcqPMMTpZ9g4vaJvDr/VSJjIuPWJao1vlbcfl0e6sL/Hv4fAAMeH5Dqcbp0Ipi1dxZ9Fvch2kTzaYNP7Q5HKeWE2pRtQ+DlQObun8vQxkNpXbY1YD242uWhLoR9GMbkpyZTNX9VANpXTKWK9bdw6URwo8xcvWL1yJUpl83RKKWc0VPlnqJh8YbMeGYGfWr3oUTOEjxe9HEyuGVgeJPhcRNUevr1JOCtAFseZnXp5wjOhFqr+T1Z5kmbI1FKOSufLD4s77L8trbRLUdz6Pwh8mbJG9fmJm4UzlY4tcMDXDwRzDtgFaN/sqwmAqVU6qmcr3LczKC0wGW7hnad3cW4LeMom7ssZXKXsTscpZSyjUveERw8d5Dq46rj5e7FuCfH2R2OUkrZyiUTwcRtE4mOiWbX67somr2o3eEopZStXK5rKDommik7p9C8VHNNAkopRTITgYi0F5HdIhIjIn63tPuKyHUR2R77+jGB43OJyFIRORj7p2OX2AMuh1+mSYkmvFLd1mJoSimVZiS3a2gX8DQwNp5th40xVe9xfF9guTHmSxHpG/v1B8mMKVE5M+Vk0lOTHPktlFIqXUnWHYExZq8xZn8yTtEGmBz7fjLwVHLiUUoplXSOHCMoLiLbRGS1iDyewD75jDGnAGL/zJvAfohIDxHxFxH/4OCULQ6tlFKu7J5dQyKyDMgfz6YPjTFzEjjsFFDUGHNORKoDs0WkojHm8oMGaowZB4wD8PPzMw96HqWUUre7ZyIwxjRO6kmNMeFAeOz7LSJyGCgD+N+x6xkRKWCMOSUiBYCzSf1eSimlkschXUMi4iMi7rHvSwClgSPx7DoX6Br7viuQ0B2GUkopB0nu9NG2IhIIPALMF5HFsZvqAjtFZAfwJ9DTGHM+9pjxt0w1/RJoIiIHgSaxXyullEpFYkz662738/Mz/v539jIppZRKjIhsMcb43dnuck8WK6WUul26vCMQkWDguN1x3Kc8QIjdQaQyV7xmcM3rdsVrhvR73cWMMT53NqbLRJCeiIh/fLdizswVrxlc87pd8ZrB+a5bu4aUUsrFaSJQSikXp4nA8Vyx8o0rXjO45nW74jWDk123jhEopZSL0zsCpZRycZoIlFLKxWkicDAReVdEjIjkuaWtn4gcEpH9ItLMzvhSkogMF5F9IrJTRP4SkRy3bHPKawYQkeax13UotsCSUxKRIiKyUkT2xlYmfDO2PdUrDaY2EXGPXVb/79ivneqaNRE4kIgUwVpD6cQtbRWAjkBFoDnww40F+pzAUqCSMaYKcADoB859zbHXMRpoAVQAnou9XmcUBbxjjCkP1AZ6xV7rjUqDpYHlsV87mzeBvbd87VTXrInAsb4F3gduHZFvA/xmjAk3xhwFDgE17QgupRljlhhjomK/3AAUjn3vtNeMdR2HjDFHjDERwG9Y1+t0jDGnjDFbY99fwfrFWAgnrzQoIoWBJ4DxtzQ71TVrInAQEWkNnDTG7LhjUyEg4JavA2PbnM1LwMLY9858zc58bQkSEV+gGrCRJFQaTKdGYH2gi7mlzamuObnF611aYtXbgP5A0/gOi6ct3czhvZ+KdSLyIVY3wq83Dotn/3RzzffgzNcWLxHxBmYCfYwxl0Xi+ytwDiLSCjgbW2Crvs3hOIwmgmRIqHqbiFQGigM7Yn9ICgNbRaQm1ifGIrfsXhgIcnCoKeZeFetEpCvQCmhkbj6kkq6v+R6c+druIiIeWEngV2PMrNhmZ640WAdoLSItgYxANhGZipNds3YNOYAx5j9jTF5jjK8xxhfrl8XDxpjTWFXZOoqIl4gUx6retsnGcFOMiDQHPgBaG2Ou3bLJaa8Z2AyUFpHiIuKJNSg+1+aYHEKsTzUTgL3GmG9u2eS0lQaNMf2MMYVjf447AiuMMZ1wsmvWO4JUZozZLSIzgD1Y3Se9jDHRNoeVUkYBXsDS2DuhDcaYns58zcaYKBHpDSwG3IGJxpjdNoflKHWAzsB/IrI9tq0/VmXBGSLyMtYMufb2hJeqnOqadYkJpZRycdo1pJRSLk4TgVJKuThNBEop5eI0ESillIvTRKCUUi5OE4FSSrk4TQRKKeXi/g94ojzU+eJvMQAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "tkl.mapping.mapOn(noised, reference, TP1=range(len(noised)), verbose=False)\n",
    "\n",
    "plt.plot(reference['x'], reference['y'], 'b-')\n",
    "plt.plot(noised['x'], noised['y'], 'g-')\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "d7e2de41-995b-4a37-a077-b0b24e40e37f",
   "metadata": {},
   "source": [
    "## For gallery"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "0af91c36-f81e-48bd-af05-7821b8a44d52",
   "metadata": {
    "editable": true,
    "slideshow": {
     "slide_type": ""
    },
    "tags": [
     "nbsphinx-thumbnail"
    ]
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiAAAAEdCAYAAAA1hdzyAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABhyElEQVR4nO3dd3iTVRvH8e/pZu+y995DQBBQNoJQAUEBZaiIKDhwC+5XRUTEAYqoICAqOBCQDSIiG2RP2auUsinQlrbn/eOkC9JBm+TJuD/X1Stp8iT5MZreOc8591Faa4QQQgghXMnP6gBCCCGE8D1SgAghhBDC5aQAEUIIIYTLSQEihBBCCJeTAkQIIYQQLicFiBBCCCFcLsDqACkVLlxYlytXzuoYQgibTZs2ndFaF7E6x62Q9xEh3Eta7yNuVYCUK1eOjRs3Wh1DCGGjlDpidYZbJe8jQriXtN5H5BSMEEIIIVzO6SMgSqnDwGUgHojTWjd09msKIYQQwr256hRMK631GRe9lhBCCCHcnJyCEUIIIYTLuaIA0cBipdQmpdQgF7yeEEIIIdycK07BNNNan1RKhQJLlFJ7tNZ/J95pK0oGAZQpU8YFcYQQQghhNaePgGitT9ouTwOzgMY33D9Ra91Qa92wSBGPajcghBBCiCxyagGilMqllMqTeB1oD+xw5msKIYQQwv05+xRMUWCWUirxtX7QWi908msKL5KQACdOwMWLcP06BAVBaCgUKgR+MoVaCJEOrTXHLh2jTD45ve+OnFqAaK0PAnWd+RrCu5w9C/Pnw4oVsHo1HDgAsbE3HxcSAvXrQ+PGcM890LIlBAa6PK4Qwo39sP0HHpr1EHN7z6Vzlc5WxxE3cKtW7MI3xcfD3Lnw7bewcCHExUH+/NCsGXTpAhUqQMGCpsCIiYGICDh0CDZuhIkT4dNPzf19+8LTT5vjhWfIqFGhMsOnnwKdgKvAAK31v67OKTyP1prRq0cDMHT+UFqVa0WuoFwWpxIpSQEiLBMXB1OmwIcfwr59ULIkPPss9OplRjcyc4rl2jVYtAhmzIDx4+Hzz83j338fypZ1+h9BOEZ6jQo7ApVtX7cDX9ouhQBg5MqR5AjMwbNNnk11+/LDy9kasZVH6j3CpC2TeGfFO4xqN8qakMIuOYsuLLFwIdStCwMHQu7cpoA4fBhGj4bbbsv8/I4cOaBrV/jxR/P4F16A336DqlXhzTftn74RHuVeYKo21gL5lVLFrQ4lrLHz9E66zejG+WvnATPKMfzP4QxbNIyTl0+mOnbS5kkUzlmY8feMp0/tPny+/nPiEuKsiC3SIAWIcKnISOjZEzp2NKdTfvvNnEq5/34IyOZ4XMmSMGqUGU257z545x1o2hR273ZMduEUGTUqLAkcS/H9cdttwgOcijrF+hPriU+Iv6XHXY65zLXr1266/cUlL/L7nt9ZdGARAIcuHEq677N1n6U69r9z/1G/WH1CAkLoWKkj1+KusTtS3gzciRQgwmV+/x1q1oQ5c+C992DXLujWDcwiKccpXRqmTzfFzdGj0KABTJ7s2NcQDtNMa90Ac6pliFLqzhvut/e/Q994g1JqkFJqo1JqY2RkpDNyils0e89sio8pzu3f3M7MnTMz9ZhP1n5CyY9LkveDvOR8Pyefr/s81f35QvIBsPfMXgD+DTfTgXIF5uLPQ3+mOvbYxWOUzlsagNuK3wbApvBNWf8DCYeTAkQ4XVwcvPSSKTbKlIF//4Xhw82SWmfq1g22b4fmzeGRR+D5582EV+E+MmpUiBnxKJ3i+1LAyRuOkYaGbkZrzYL9C5K+X31sdaYeN2XrFHIE5GBkm5HkCcrDjtM76DerX9LoxtXrVwE4fuk4QNKpmG7Vu7H51GYuRl/kVNQpYuNjORV1itL5zH+dKoWqkDsoN5tOSgHiTqQAEU51+jS0b2/mdjz5pFlaW7Om616/WDFYsACeego+/hh69DCnfoT1MtmocA7QTxlNgIta63AXRxW36K2/3uKrTV8BUKlgJdYcX8PgPwaz4L8FaT5Ga83+c/vpVLkTrzR/hQI5ChB5NZJp26bxzMJnmLR5UtIcjqOXjgIQHRcNQOtyrYlLiKPdtHYUH1Ocx+Y+hkYn9f/w9/OnYI6CXIy5mPR6CTohaSRFWEMKEOE0+/aZPh1r1pjVLuPHO3/Uw56AAPjsM/jkE3MaKCwMrl51fQ5xk6LAP0qprcB6YJ7WeqFSarBSarDtmPnAQWA/8DXwpDVRRWZcib3CyJUjeefvd5Juu6P0HWwK38RXm77i9eWvp/nY01dOExUbReWClQEI9g9OOsUC8Nqfr5EjIAcAfx/5m+HLhvP0wqcBaFOhDTkDc7Lh5AYApm6dCpB0CgYgPiGeVcdWcT3+OnEJcTT4qgHVxldj2cFlDvrTi1slBYhwivXrTR+Pq1fh77+hXz+rE8Ezz5heI0uXmiJERkKspbU+qLWua/uqqbV+z3b7BK31BNt1rbUeorWuqLWurbXeaG1qkZaRK0eSe2Ruhv85nN61eifd3qlSJwBqFKnBpvBN7Du7z+7j95/bD5gRE4DggGCOXDxCgF8A4zqOIzwqnCMXjwBm5GPkPyOTHlsiTwm+6vzVTc9Zo0iNpOsftf+Ig+cPMmvPLGbunMnWiK0ArDm+Jjt/bJEN0gdEONzixWb+RdGipkdH5cpWJ0r2yCNmRKR/f3jwQbP819/f6lRCeL61J9ZSKm8pPunwCd2rd+eeyvew/fR2ulfvzuKHFlO9SHXKjC3Dj9t/5M2WbyY97tjFY1yIvsChsycAOH+kNCtOwqWo6wA8UP1B2ldsD8D2iO2UyFOCAXUHANCndh+2nNpCgF8AD9V5iFNRpyiVtxQxcTHM3TeXknmTF0z1rNGT4cuG89Hqj7iecJ2aRWoSHRfNllNbXPMXJG6itL5pQrllGjZsqDdulA84nmzxYjO6UK2a6fVRrJjViewbOxaee870IZk40fErcbyFUmrTjd1J3Z28j1ijwVcNKJa7GPMfnJ/mMS2/a2kmi75ykWPHzAeUl/bfxvkcKZrbfrYPzlWGt2w/lN/Pp/i1dkQ8mpsEvxg+aDaBl9s+nqWM07ZOo9/vZjj2s7s/Y+XRlWwK38TQRkMpm78s3at3z9LzivSl9T4ip2CEwyxfDvfea4qPP/903+IDYNgwGDECvvnGTJAVQmTP8UvHU825sCd/cEEuxVxCva0oM3Qgj/1vZeriA5j8TRDLUkzLGDGgMa3uCiDgbF0Ir88r7Qdy111m9DLuFvuKPVTnIdqUb0OQfxC9a/emWuFqHDx/kOcWP8dHqz+6tScT2SYFiHCIf/6Bzp2hYkVYssTszeLu/vc/0xTtlVfMBnhCiKyJjosm8mpk0rLXm+6PhnHj4J+xKeYQN/iWwk/eD8DeoXspktMsn27fNpDWrZMPe3dEIaZPh6MfzGHdkGX87x1/Tp40Wy5UqGAmuCckZC6nUoqfe/7MqkdWUThnYYrnTm6quzVi6y03TBPZIwWIyLZdu8ymcaVKmQmentKGQSnToKxOHejTB/bKijwhsuSnHT8BUDZf6g2YtDZNAStWNEvha4S0ZXydf2lYwozGn4k+Rd2idalSqArKdh7UX5lJWVsHb2XNo8kTRIvmLkrj2gV47TXzszp3LhQvDgMGwO23myX+mVEgR4Gk17+j9B1Jt1+9fpX/zv2XpT+/yBopQES2nDxp2qqHhJj5H+582sWeXLlg9myz027PnuaTmhAi88asHsPDsx/mzrJ30q16t6Tb9+2DVq3goYfMNgnLl5sVcU92q8+GxzZQtVBVAO6vaUZBEle/BPiZtRF1itahSakmdl/Tz8+MuK5ZA99/D+HhpuHgCy/c2s9wnaJ1+KLTFyx+aDEAm8M33/KfX2SdFCAiyy5fhnvugXPnzCkMT919tmxZmDrVdE198UWr0wjhObTWvL3ibTpU7MCSvkvIHZQbrU3Pn3r1YNs2M8l77Vpo2TL1Y0NzhQJmdQrA7F6z+fX+XymUs1CmX9/Pz6xm27sXBg+GMWNM76GDBzP3eKUUTzR6gpblWhLsH8zmU1KAuJIUICJL4uPhgQfML+1ffoH69a1OlD0dO5qJqePGmb1qhBAZO3LxCJdjL9O9eneC/IO4dMnsTj10KNx5J+zYAY89Zn9369qhtbmr7F1ULmTW6RfOWTjLq1By5YIvvjAfhI4fh0aNzET4zAr0D6RO0Tr8feTvLL2+yBopQESWvP66aXE+bhx06GB1GscYOdIUUg8/DKdOWZ1GCPe3PWI7YIqJvXvNXIx588wy9wULoESJtB87rtM4lvVzbBfSjh1NE8RixcwWEF9/nfnH9qzRk3Un1kl7dheSAkTcsp9/Nr+sBw0yw57eIjgYfvgBrlwx+9a4UYscIdzS9tOmADm1vRaNG8PZs7BsGTz7bMa9dZRS+Ps5vgtgpUpmbki7duY96pNPMve4vnX74q/8mbZtmsMzCfukABG3ZNs2M+u8aVOzv4q3qVYN3n4bZs0yhZYQIm3bT2+nSGA57r83D+XLw8aNcNddVqeCvHnN5PLu3c2p1fffz/gxxXIX447Sd7Bw/0LnBxSAFCDiFly+DPfdB/nywa+/mhEDb/T88+Yc8pAhEBlpdRoh3MulmEusP7EegL937yRyRy3uuANWrIAyZSwOl0JQkGlW9uCDpungR5noM9auQjv+Df83zf1qhGNJASIyRWt44gkzu3zGDLP+3lsFBMCkSXDxoqyKEeJGrae05vZvbmfU6OucvBhB2YKlWLjQfDBxNwEBplHZAw+Yn+VJk9I/vmu1rmg0VcdVZcXhFa4J6cOkABGZMmWKaSj01lvQooXVaZyvVi0zEjJlSuYbHAnhjeIS4pi+bTqx8bEcv3ScTeGbAHglKgfkPk2PLnnJkcPikOnw9zfL7Nu3NytyFi1K+9jaRWsnXd8VucsF6XybFCAiQ3v2mNMRLVvC8OFWp3Gd114z3V2HDDHLjoXwRZM3T+ahWQ8xbes0So9NbrUeFF8AgIsx56yKlmlBQea0ce3aZjRkz560j32n5TsAVCtczUXpfJcUICJd0dHmBzZnTjMC4ktb1+fKBR9/DFu2wIQJVqcRwvUOnT/EoD8GAfD73t+Tbi974nl+e+gHAOK1Z1TnuXObianBwWbriHNp1E31i5umRrmDcrswnW+SAkSka/hws/JlypT01/R7qx49oE0bM4nt9Gmr0wjhWltObUm6vnj/kqTrK0e+Sseqbfi4/ceMajvKgmRZU7asWeF29KjZ/8neJnaJG9I5Y4mwSE0KEJGmlSvNGvonn4ROnaxOYw2l4PPPISrK7J4rhK84c/UM3WeazqTv3/4NsQkxSfeVLlQIP+XHsKbDKJLLQ3aftLnjDvj0UzMXZPTom+9PHNE5FXWKGTtmuDidb5ECRNgVFWX6fZQvD6M85wOOU1SvDgMHmtMw+/dbnUYI59txegc9ZvZI+v6b53pCXAgAL93xklWxHObxx83mkyNG3DzJPC4hDoBev/Si16+92HBigwUJfYMUIMKul1+GQ4fgu+/MuVNf99Zb5tyxL03CFb6rzdQ2rDiSvAz1xMG8EGC2mU25hb2nUsq0aS9TBnr3hvPnk+9L0Oa8zOXYywCMWuXjn8CcSAoQcZNly8zGTs8+6xtLbjOjWDGzLPfnn81eEyL7lFKllVLLlVK7lVI7lVLP2DmmpVLqolJqi+3rDSuy+pocAWZdbdn4NvBROF9+mXxf09JNLUrlWPnywU8/wYkTpltqolqhtZKud63Wld92/8aJSycsSOj9nF6AKKXuVkrtVUrtV0q94uzXE9lz+TI88ghUrQrvvWd1GvfywgsQGgovvST7xDhIHPC81ro60AQYopSqYee4lVrreravd1wb0TfdVc70Uz/yyRQGP1SMhx+Ge6veC0BorlArozlU48bw6qtmkv3cuea2mkVq0rRUU95p+Q7vtHwHjZb27E7i1AJEKeUPjAc6AjWA3mm8wQg38frrcOwYTJ6MWzcXskKePPDGG6bl9K1s9S3s01qHa63/tV2/DOwGSlqbyvdcib3CgN8HJLVXB4g8F4M6V5kmNUsmbeb26/2/Ej0i2pqQTvT661Cnjtm47tw5s0ne6kdX8/pdr1MrtBYl85Tkzb/eZFvENqujeh1nj4A0BvZrrQ9qrWOBn4B7nfyaIos2bjQrPp54wmw2J242cCCULGnmhMgoiOMopcoB9YF1du5uqpTaqpRaoJSqmcbjBymlNiqlNkbKBj63ZNiiYUzZOoU3/3oTML1/Vm48i39sIX75JXnPJ38/f4IDvG8DqKAgM9ftzJnUp2LAFCMD6g3gxOUT1J1Ql84/dObk5ZOW5PRGzi5ASgLHUnx/nBs+4cgbh3uIizNtiosWzdzOkb4qcSLqP//IKIijKKVyA78Cz2qtL91w979AWa11XeBz4Hd7z6G1nqi1bqi1blikiGctC7XSb7t/4+t/v6Z03tIsPrCYE5dO8PLLEKUjqV+1ECV9ZDyqfn1zanXqVDPCmdK7rd/lzItnGNlmJH8e+pPBfwy2JqQXcnYBouzclupzo7xxuIdPPzUdPz//3D03lXInjz5qRkHefltGQbJLKRWIKT6ma61/u/F+rfUlrXWU7fp8IFApVdjFMb3W+yvfp1ZoLRY+tJAEncCImdP47DNNYNEDNKlawep4LjViBJQrZ0aAY2NT31coZyFeaf4KbSu05ejFo5bk80bOLkCOA6VTfF8KkPErN3P4sJnb0KULdO9udRr3FxxsJq6tXAnLl1udxnMppRTwLbBba/1xGscUsx2HUqox5j3rrOtSeq8L0RfYfGoz91W/jxpFatC4WHO+3zGZGref5LqK8rm9UHLmhHHjYPduGDvW/jH5QvJxKebGQTqRVc4uQDYAlZVS5ZVSQUAvYI6TX1PcoqeeMuvix40zlyJjiaMgMhckW5oBfYHWKZbZdlJKDVZKJY5z9wB2KKW2Ap8BvbSWv3FHWHlkJQk6gZblWpKQAJdW9SK+wD56vb4AgKqFqlqc0PXuuQe6dTOjm8eO3Xx/3qC8XIy56PpgXirAmU+utY5TSg0FFgH+wCSt9U5nvqa4NfPnwx9/wIcfmqY8InNCQuCVV0zxtnIl3Hmn1Yk8j9b6H+yfpk15zDhgnGsS+Y55++YR9lMYAE1KNeGLL2DP1rxQEXZH/4m/8qdhiYYWp7TG2LHmfXH4cJg2LfV9iSMgWmuUfFrLNqf3AdFaz9daV9FaV9RaS2cJNxITY5qNVa0Kz9zUAkpk5NFHoUgR+/tJCOHOFh1YBMCwJsM4eTSEl1+G2nVMB9Blh5bRuGRj8oX45mSwsmXhuefg++/NysCU8gXnIy4hjmtx16wJ52WkE6oP++QT+O8/cxkUZHUaz5MjhxkB+eMP2CnjesKDHL14lFqhtfio3cc8+ij4+0P/AaYAOX3lNO0qtLM4obVeecU0HXzuudSnWPMG5wXgYrSchnEEKUB81MmTZnfXsDC4+26r03iuJ580k9c++sjqJEJk3uELhymXvxxffQV//QUffwyNK1VKur9thbbWhXMDefPCO++Y06uzZiXfnjgqJBNRHUMKEB/10kum90das71F5hQqZJqTTZ8Ox49bnUaIzDly8QgF/cry4ovQrp05ndiibAsOPXOIH7r/QPMyza2OaLlHH4WaNc17ZeKy3KQREJmI6hBSgPigVavML8wXXoAKvrXU3ymGDYOEBNNLRQh3dyH6AheiL7BhSTm0NrvCJs6nLJe/HL1r95YJlkBAgJnfdeAATJpkbssXLCMgjiQFiI+JjzfzFkqVMr0sRPaVKwcPPABffQUXLlidRoi0xcbHMnLlSAB2r6jF22+bSZfCvrvvhjvugHffNS3qE0/ByBwQx5ACxMd8/z1s3myW3ebKZXUa7/HCC2Yn4cmTrU4iRNoGzhnIh6s/JOfeh6mds4OsfsuAUqb4OHECJk6EPEF5ALgce9niZN5BChAfcvWqaTfcqBH06mV1Gu9Svz40a2aaucXHW51GiJtFx0Xz6+5fqRU7kKs/TuKrCYrAQKtTub9WrczX+++Dvm62CL92XZbhOoIUID5k7FhTyX/0kXQ8dYann4aDB2HBAquTCHGzvw7/xdXrV9n1S3cGDZIdr2/F//4HEREwfUoIYIo5kX1SgPiIiAj44APo2lW6djpLt26mPfvnn1udRHi71cdW3/IvwT/2zcMvPgcFL7Xkgw+cFMxLNWsGHTrApx/ZRkCkEZlDSAHiI95+20yiGjXK6iTeKzAQBg+GxYthzx6r0whv1fWnrjSb1IxXlr6S6cdorZm5eR4JB1ozZlQOChRwYkAv9dZbcPZ0EAolIyAOIgWID9i920ygevxxqFLF6jTebdAg01V2nOxeIpxk9t7ZAIzfMD7DY19Y/ALtp7Vn4+E9RMYdonzcPfTt6+yE3qlJE2jZUkFcCFdipABxBClAfMCrr5punW++aXUS7xcaaib4TpkCF2WlnnCC1+98HYC4hDhOXj6Z5nFxCXGMWTOGJQeX8ODk4QB8OuQemf+VDa++CjrgGh+vG02CTrA6jseTAsTLrVsHs2fDiy+ajdOE8z39NERFwXffWZ1EeKM9Z5LP7zX5pgnNJzVnxLIRbDixIdVxM3bMSLr+n//v5I+tRZc7Zcvr7GiXYoucQ+eOWhfES0gB4uVGjDCFx7PPWp3Ed9x2GzRubE57pdzISghH2H1mN52rdObnnj/jp/xYdWwV7//zPi0mt+Ds1bPEJcShtebD1R9SK7QW9Q9PBaBPo04WJ/d8SkFYoRcBmDkvwuI0nk8KEC+2bJn5GjEC8uSxOo1vGTQIdu2C1autTiK8SXxCPPvO7qN64er0qNGD+sXrA6ZBVkx8DIVHF+a2ibcxfft0tkVso32+p9g85UG6Bkzg9bbDLE7vHYbfez8A3/wUIR8wskkKEC+lNQwfDqVLm8mnwrUeeMAUfRMnWp1EeJNDFw4RGx9L9cLVAahaqCoAbSq0oVz+cjx7+7Mcv3ScvrP6kiswF8s/602Z0n78+PzjFMtdzMroXqNE3qIAHIyI4O+/LQ7j4QKsDiCcY/ZsWL8evv0WQkKsTuN7cueGBx8080A++QRZ9igcInH+R7XC1QB4q+Vb1ChSg751+iZtIDes6TAem/sYeS405dd1eZgyRd4DHCk0VygAOYqe4NNP4a67LA7kwWQExAvFx8Nrr0HVqtCvn9VpfNfjj5veK99/b3US4S0OnDsAQKWClQAICQihX91+qXavLZOvDHPvX8TmsW9Rt64phIXjBAcEc1vx28jbaDa//266H4uskQLEC/34I+zcadoHB8gYl2Xq1TP77shk1LQppe5WSu1VSu1XSt3UWUsZn9nu36aUamBFTndx4PwB8gTloXDOwukeN2GC+cU4ahT4+7sonA95pP4jRPhtwa/EZun5kw1SgHiZ2Fh44w2zOdp991mdRgwaBDt2wNq1VidxP0opf2A80BGoAfRWStW44bCOQGXb1yDgS5eGdDMHzx+kQoEKqUY8bnTxIrzzDrRpA+3buzCcD+ldqzfB/sFU7DGZb781O2GLWycFiJeZNAkOHTI7N/rJv67levUy80G++srqJG6pMbBfa31Qax0L/ATce8Mx9wJTtbEWyK+UKu7qoO7iwPkDVCxYMd1jPvwQzp41l9J0zDkK5ChA6/KtiS21lEuXpOdPVsmvKC8SEwPvvZe8cZKwXuJk1Bkz4Px5q9O4nZLAsRTfH7fddqvH+IQEncCh84eokL9CmsecOGF2ve7TBxr49Mkq5zt77SwVQ0vQpAl89hkkSGPUWyYFiBeZPBmOHzct19Wli1C+vOlCtmKF1dF8WuJk1OnTrU7idux9Pr9xtkxmjkEpNUgptVEptTEyMtIh4dxN+OVwYuJjqFAg7QLk7bfNJPR333VhMB+ltcbfz59nn4X9+2H+fKsTeR4pQLxETIw57XLHHdD2tvNw991w+DCcOQO//ur8ALGx5kvcpH598yXDtDc5DpRO8X0p4MbNTTJzDFrriVrrhlrrhkW8dM+B45eOA1A6X2m79x84YE7BPv64+ewhnKtc/nIcvnCYck03EVo+kgkTrE7keaQA8RLffQfHjsGbkUNR5culnvU4axZ06gSnTtl/8G+/QefO0LevaRyScsnGd9/B+BS7br79tmmvmlJ8PLRsaaqf69cd8wfyMgMGwKZNsH271UncygagslKqvFIqCOgFzLnhmDlAP9tqmCbARa11uKuDuoNTUebnt3hu+1Ng3n0XAgPNhmnC+crlL8eBcwdoMrkhul8r5s+HI0esTuVZpADxArGxZvSjSb5dtPtvPHTpknwC+I47zHmZBQvg66/tP8Fnn5nTNN9/DwMHmuIlcRj74Ydh6FBTZJw6BW+9BW3bph7tmDAB1qwxv2FHj3bqn9VT9eljfjlMmWJ1EvehtY4DhgKLgN3ATK31TqXUYKXUYNth84GDwH7ga+BJS8K6gfAoU3fZ62j6338wdSo8+SQU99kpuq6VIyAH8ToegEi1E6Xgm2+S75+9Zza9funF9gj51JEWKUC8wHffwdGj8FbIKFSfPqaQ+Osv+PJL+Pln+OADaN7cHHjjTCmtYetWM1Pyt9/MbZ06mX3lN29OPi48HP78M/n7uXOTr48aZdoB3nefWf93Mu0twn1V4cJmkOn772WQKCWt9XytdRWtdUWt9Xu22yZorSfYrmut9RDb/bW11hutTWydU1GnUKikTpwpvfOO6Xb60ksWBPNRN87F6djRFCCJP99dZ3Rlxs4ZzNl746CeSCQFiIdLHP24/bbrtI+YBlWqmDvy5IHBg6FECXj5ZdOQ4uDB5PkgWpvrmzfDhQtQty4EB5v7Llwwl337Jr/QwYPmHA+YtX1//WWuX75sbu/Qwex6FxMjk17TMGAARETAokVWJxGe6FTUKQrnLEygf2Cq23fvhh9+MAOVRYtaFM4H9a/XP9X39z1yjFOnzGez2PjkEeLXlr/m6mgeQwoQDzdlijnv+NauB8ya/9at7R94993m8v774d57zehIjx5m73gwH9GDgpKPv/9+00410eOPwyu2RpVt2sAXX0COHPDoo+a2qlWhdm3ImVO6bqWhY0ezKEkmo4qsCI8Kt3v65Z13zI/diy9aEMrHzeszjycaPgHAI9vLwOvBDF3WP2nPHpE+pxUgSqm3lFInlFJbbF+dnPVavur6dTP60Tj/XjoUWG9mOLZoYf/gIkWSd02aMweGDEl9/7lzyQVIYKAZNUlpT4ofqP79zamc6GhzigfMyEtAANSsaVp/ipsEBsJDD5m//rNnrU4jPM2pqFM3FSA7dpgeM08/bT5DCNfqVLkTX9zzRfIN/rGEh07ll41/Jd+k/NGyF4Ndzh4BGau1rmf7klXSDjZ1qllp+2b8m6gO7c0v//QsXw5LliR/P3q0mWQK5rfjtWvmelgY1EjREfvsWdNe9fhxWL3azBfp0iX1c1eunPw8mzZBVJRZAvzJJ7BvX3b+mF5lwABTOP74o9VJhKc5e/UsRXKlXmL89tum2d3zz1sUSgCwadCmVN+P3TCKQL9A3m/9PvE6nsux0qvdHjkF46Hirmvee0/TqPY1Ol6eYdqfZkQps4Il0cCBUK6cuV6iBLRqZdbwff21mdG2dKmZUFqwoDmuZElo2tQ8z+zZpshYtsx02EqcP/Luu2Yzim++Mad6hg0zu+IJAOrUkZ4g4tZ8svYTXvvzNc5dO0fBkIJJt2/bBr/8Yn7EChZM5wmE0zUoblYdlshTguLn7yNKnaRq4WpJE4YvRF+wMJ37cvZeqUOVUv2AjcDzWmtpRp0dR4/ClStQvTo/tZzAoUNP8Mkza1DbMYVBZr32mjklkz8/DB8ODRsmzxF5//3k49q0Sfs5lIJcuW6ec9KqlTkNNGxY8m3//pv5bD6gf3949llzxqx2bavTCHcWnxDPsEXmZ0mhUu2CO3KkmWv+7LMWhROpHBt2jJyBOVm6MAcP/HaZmsVvT/r3OnLhCGXylbE4ofvJ1giIUmqpUmqHna97MbtWVgTqAeHAmDSew+tbKDvEf/9B2bJQowYJz7/IyNV3UovtdP7UNqJRIe32zDf53//MSWMw8zY6OXh6TmJLwDx5TF/43bvhuedkcqpNnz7mr116goiMVB1XNem6v58/PWr0AMzbwcyZpu9HgQJWpRMplcpbioI5CtKtcw5CFy/i+qJ3aFmuJSEBIfy4Q8652pOtERCtdduMjwKl1NfAH2k8x0RgIkDDhg1lpk5aWrVKujr3433soibT6YNf4rYYISEWBbOjRg0zV6R4cfMxX2uzQ9aZM9CkidXpLFekSHJPkJEjzbQZIW4UFRvFgfMHkr5/oekL1Aw187w+/NDMGU850CjcQ+Jk888/h+tR+WhdvjUrj660OpZbcuYqmJT9+LoBsjQiK44eNXMpTpwAQOfIyfsMpwIHuP/oGKhVy+KAaWja1MwbSbkl5/r1lsVxN9ITRGTk33Bz6vKP3n9w9qWzvN/GnB49dsyMnj36qPT9cFeJk82nT4f8Ifm5ev2q1ZHckjMnoX6olNqulNoGtAKkVs+KgQPNuk2AoCCWN3qJ9dzOS3xIQOnisGGDmfTprkqWNMt1nngC9u5NbnLm4zp1Mssm5TSMSMv6E6Zgb1SyEQVzFEQpszHwmDFmUFH6friv2rVNi6XJk03L9mvXr1kdyS05rQDRWve1tU6uo7UO89UNpG5ZRITpKHrmjPk+f35z2bo1rFvHyL3dKUY4/T+0jXyEhEDevJZEzbS+faFbN3NdRkEAM0zbp4+pLc/L1Gxhx/oT6ymXv1yq1uuRkTBxohniL1vWwnAiQw8/bHa5uHIhJ9fipACxR5bhupMTJ0yzsPffN3MmIHmMdelSNlyvx9KI2jzHx4R0ztT0G/fRtKlpWPDTT1YncRv9+5tW+jNmWJ1EuKP1J9bTuGTjVLd9+qnp/3djn0Dhfnr3NvN09u2SEZC0SAHiTp5/3jT7qlcPvv3W/HYKD4fq1UEpRo6E/Pk1gzcMNLd5kty5TQOzH380/UME9eub3nFTp1qdRLibiKgIjlw8QuMSyQXIxYswbpzZ87FaNQvDiUwpWBC6doVd23IQEx9Dgk7I8DG+RgoQd3HqlGnu1bcvvPeeORUzezbs2gWVKrFrF8yaBU89pcjTsGrGz+eOevUyH9+WLrU6iVtQyoyCrFkjzWJFamuOrwFINQLy5ZemCHn1VatSiVvVty9EX8oBIKMgdkgB4g4uXDBLVqOjzQZxHTqYUy/Tp5s9WG67jVGjzIZTie07PFKzZqZ5Wbdu0pzM5sEHwc8Ppk2zOolwJ4sPLCZXYC5uL3U7AFevwscfm36BKReWCffWvj3kDMwJQJ/f+jB61WiLE7kXKUDcQeLmbTVqQMuW4O9vmnj98QdozeGqHZg+HQYN8vANpwID4Y47zPXEPWh8XIkS0K6dOQ2TICO0AtBas2D/AtpUaEOQv9kgctIkMwF1+HCLw4lbEhQEjevnAmDO3jm8tPQlixO5FylA3MHu3eZy3jxTfIBZRB4fD8WL89HKxvj5ecmGU6NtnwAefNDaHG6kXz/T7uXvv61OItxB1xldOXzhMB0qdgAgLs4svb3jjrQ3uxbuq3+7RlZHcFtSgLiDbdvM+ZUyKfYKuOceACK7P863k/zo1w9KlbIonyMldmz1ij+MY3Ttaga8pCeIOHzhMHP2mr4/d1cy+zPNmmV2vX7hBQuDiSzr274WxAVbHcMtSQHiSHFx8NFHsHFj5o7X2rTMGz/eLFP1S/HPMW4crFnDuLzDiY72ojef2Fhz+ddfpkuPIGdOuP9+s7PplStWpxFWSiw+ACoUqIDWZvSjYkUIC7MwmMgyf39F64QPkr6Pjou2MI17kQLEkUaMMO0JG2VyyG3VKvOx96GHzPLUlJTiap0mjJ8YSFiYFy27SyxAvv4aHnkEDh60No+b6NfPrE6eNcvqJMJKW09tJX9Ifq4MN5Xo6tWwbp3Z8yXx7KzwPJ/0ehZmfQfA8UvHLc3iTqQAcaSaNZOv2/ZuSde335r+GF9+aXYou8HkyXD2rJe1XC5fPvX3o0ZZk8PNNG9u/mp85TSMUmq0UmqPUmqbUmqWUip/Gscdtm3psEUplcmhRc+1I3IHDYo3SFo5MWaM2e12wABrc4nsqV0byhYuBsDRi0ctTuM+pABxpHPnkq+3bGlmjR44kPqYLVugShXYtMnsp92rl1maeoP4eLPsrkkTs3rVaxQsaLaJ/P13M/Lzyy/m1JWP8/MzPQOWLTO96HzAEqCW1roOsA9Ir7tFK611Pa11Q9dEs8aV2CusP7Ge2qG1AfjvP/Nj8uSTdt8ihIcp3XIxAC8teMPiJO5DChBH2rDBbL7266+wf7+pICpXTt146/nnzTtL48Zmcf+jj9p9qlmzzNmJF180Dau8ytChZoffbt1M0dazp9ni08f162emBX3/vdVJnE9rvVhrnVh5rgV8flby3dPNpNNaoWafp08+MSvXhw61MJRwGB26FYBNkassTuI+pABxpA0bzPyP7t2Tb9ParGjp2tV82v/zT1NRJCSYvh+3337T02htVqtWqmR+T3utzp1hyBCYPx9q1fL5DqkVK5rRrqlTzf8BH/IIsCCN+zSwWCm1SSk1yIWZXEprnbT7bc8aPTl71pyCffBBKFbM4nDCIcoWLpp0fU/kXguTuA8pQBwlKsqMbCS2Kfz5Z9MzuW9feOAB0+OjZ08oV8586m/RwtxvZ3hj5Uqzaexzz3n5xLOgILPaZ88eyJdP5oNgWrPv3p35hVTuTCm1VCm1w87XvSmOGQHEAdPTeJpmWusGQEdgiFLqzjRea5BSaqNSamNkZKTD/yzOFh4VTmx8LJ93/Jx8IfmYMAGuXTPvAcI7fNX5Kwbl+w2uFOHOSa05dP6Q1ZEsJwWIo+zcaS5rm/O39OhhdrWdOtV8LV5shjQ+/hjy5zddpx56yO5TjR5tOp76zMSz8uXNCIjsS8/990NwsHdsUKe1bqu1rmXnazaAUqo/0Bl4UGv7Yz5a65O2y9PALKBxGsdN1Fo31Fo3LGJnQre72xaxDYA6ResQHW2mSXXoYH4shHfIHZSb//Xphpr6J+evneOTtZ9YHclyUoA4yvbt5rJOHfv3t2plRki6dUv3aXbvNh3Yhw6FHDkcnNGd5ctndtrycfnymbN1P/6YvGLZGyml7gZeBsK01lfTOCaXUipP4nWgPbDDdSldJ7EAqR1amx9+MHtRekXnY5FKaCi0qlmLwMiGbDi5weo4lpMCxFG2bzdT1cuVy9bTjBljmoUOGeKYWB5DCpAk/fub5dfz5lmdxKnGAXmAJbYlthMAlFIllFLzbccUBf5RSm0F1gPztNYLrYnrXNsitlEqbynyhxTg44/N55i2ba1OJZyhZ0+4dqAR/57czPX461bHsZQUII6yd6/pFuaX9b/S8HCzK+rDD3v4pnNZIQVIknbtzMRDbzgNkxatdSWtdWnb8tp6WuvBtttPaq072a4f1FrXtX3V1Fq/Z21q51hzbA3Tt0+nTtE6/PmnOZs7bJgXrn4TgFmjoE42JiYhms2nNlsdx1JSgDjK4cNQoUK2nuLzz80edD458SxfPnPOISbG6iSWCwgwqx/mzYMzZ6xOI5zpyIUj3DHJ7BBdJ7QOn31mPnz06mVxMOE0oaHQrHhbSPBj9p45GT/Ai0kB4ghaw5Ej2Tr9cvmyaYjavbuZq+pz8uUzlzIKApieINevw08/WZ1EONOzi55Nul7Mrw5z58KgQcl7Ngrv9GC3wnC0BTO3zbY6iqWkAHGEU6cgOjpbBcikSXDhgpe1Xb8VefOaSylAADMHoF4932nN7ouuXr/KH/v+SPp+88I6+PnBE09YGEq4RPfuoPbey/5LOzh43nf3w5ICxBEOHzaXWSxA4uPh009NEyo7fcl8Q2Kv6WvXrM3hRvr3N/1Adu2yOolwhn/D/yUuIY4uVbrgr/z5/dsqdO8OpXy+J6z3Cw2FxgU7ArDs4J8Wp7GOFCDZ9f77pnEHZPmdY/ZsOHTIR+d+JAoONpfRslV1ot69TSM6b56M6svWHl8LwDdh3zC+aBwXzwXy9NMWhxIu0/+eqnClMH9s+8fqKJaRAiQ7PvoIRoxI3kM9iz2TP/7Y9OLy6rbrGUk86S0FSJKiRaFjR7M3THy81WmEo609vpYKBSpQJGcon39uTrl51caTIl333afgWHNWHpECRNyquLjUEzZat4YsdGDcsAFWrYJnnvHytusZSRwBOXvW2hxupl8/OHHCbCEkvMv6E+tpXLIxy5ebpbdPPy1Lb31JaChUCmzBeXWAh2c/Qr0J9Th/zbe6QUsBklXrzcZRvPYajB1r9s3OwrvH2LFm/uUjjzg2nsdJHAEZPNjaHG6mSxfTuV8mo3qX89fOc+zSMeoVrcfnn5ult717W51KuNp9jZoD8N2WyWyN2MoH/3zAy0teJjbei9sgpxBgdQCPtWiRaTo2bBgULJilpzh2DGbONKMfefI4OJ+nSdz04rpvdwa8UUiI2ctw6lSzVNvn/594iR2nTUf5otRhzhx4+WVZeuuLnuxen1HfJH//4eoPAdh8ajNze88lOCDYomSuISMgWbVoETRqlOXiA8xGsFojE8/AvPs++qi8C9vRv79ZHPTLL1YnEY6SuPfLurm1UUqW3vqqMiUDabIogXqzU0/yWnJwCX8e8v7zrlKAZEVEhJm80aFDlp8iKgomToT77oOyZR2YzZMVK2b+br15F7YsaNIEKleW1TDeZPvp7RQIKcBPX5ekWzcoXdrqRMIq3bsptmy++VfxycsnLUjjWtkqQJRSPZVSO5VSCUqphjfc96pSar9Saq9SKuu/qd3RV19BQkK2+iV/951pPDZsmMNSeb477jB/r7N9uzvgjZQyk1H/+iu55YzwXLsid7Fg/wKKJNThwnnFU09ZnUhYKXGD9PL+ZgnUpLBJACw9tJSdp3daFcslsjsCsgPoDvyd8kalVA2gF1ATuBv4QinlHWs8Pv4Y3nwTSpaE6tWz9BQJCabxWJMm0LSpg/N5sg4dTDO3L7+0Oonb6dvXXE6bZm0OkT2rjq6iwVcNuHb9Gix/h5o1oUULq1MJK1WqZKbAlfxzMedfPs/D9R8m0C+Qn3b8RJ0JdayO51TZKkC01ru11nvt3HUv8JPWOkZrfQjYDzTOzmu5jeefN5f335/lp/jjD9i/X0Y/buLvbyY8LF8uu7DdoGxZs9L7m2+kJ4gnm7tvLgk6galNt7NvyZ088YQsvRVmFGT1ipxcv5wfgOsJZjJ+gk6wMJXzOWsOSEngWIrvj9tu82wXLpjLxo1h1KgsP83YsVCmjNkPQNwgNtYUIvnzW53E7QwZAkePmgJWeKaD5w9SLn85Zk4qSq5cySNbwrd162ZGxufONd8/1yS5LfblmMsWpXK+DAsQpdRSpdQOO1/p9e20V9PrNJ5/kFJqo1JqY2RkZGZzWyOxV/oHH0BgYJaeYvNmcy7/6afNtuviBoUKmY/4J0/KR/0bhIWZyYqff251EpFVB84foHTuCvz0E/Tpk7wHo/Bt9eqZUc7EptpjOoxhZo+ZAF69WV2GBYjWuq3Wupadr/RmCh4HUs7rLgXYndKrtZ6otW6otW5YJAudRF3m4kWYPNlcb9Qoy08zdizkzg0DBzool7epWNFcli0LL71kbRY3ExBglmsuWwa7d1udRmTFwfMHiT1VkWvXZOmtSKYUdO0KS5aYfj8AFQua98ID5w9YF8zJnHUKZg7QSykVrJQqD1QG1jvptVxjrdk4imnTTAWRBSdPwk8/ma6n+fI5MJs3KVEi+frEidblcFMDB5qu9ePGWZ1E3Krz185zIfoCe9dV4PbboX59qxMJd9KtG8TEwMKF5vsKBSoAPj4Ckh6lVDel1HGgKTBPKbUIQGu9E5gJ7AIWAkO01p49nr5qlZmbkI0d48aPN1vISOOxdKQcXcqd23RqS9SgAdSp49MTVIsUMau/p0wxg3LCcyR+ko3cW1FGP8RNmjc3LfkTT8PkD8lPwRwFOXBORkDs0lrP0lqX0loHa62Laq07pLjvPa11Ra11Va31guxHtdCFC6YNZd26We6FffUqTJhghtkSzzIIO5Qy45Bly8KpU6bd/YgREB5uJtBs327OQfiwoUPhyhX4+murk4hbcfzScQDy6NLZWUQnvJS/v5nnNW9eci/GCgUqcPCCjID4rogIKFDAnHRv3jzLTzNtGpw7J0tvM6Vt2+TN/gDefx9efz35+0WLXJ/JjTRsCG3awJgxEB1tdZqsUUq9pZQ6oZTaYvvqlMZxd9uaGe5XSr3i6pyOFHXRTFy/p0scOXJYHEa4pW7d4NKl5N2va4fWZvGBxai3Ff8c/cfacE4gBUhGRoxIvt6sWZaeIiHBTD697bZs1TC+JTQ0da+Vb79Nvq7tLqjyKcOHmwGixHnRHmqs1rqe7Wv+jXfamheOBzoCNYDetiaHHmn54lwAdAy7anES4a7atjVnnhNPw7x+Z/IHr/UnPHsapT1SgKRHa1iQ4uxRFtuWLlwIe/eaVbzSdOgWTJgAr72WvAlK/vxmDsi5c5bGcgetWplOuqNGefUGwo2B/Vrrg1rrWOAnTJNDjxMfD/N+zwlAoaJSgAj7QkKgY0ezG0VCApQvUJ5u1Uyv9uK5i1uczvGkAEnPvn1m6crYsWbzuSzuGDV2rOnc3rOng/N5uwIF4H//M92afvzRfCwoWhT++8/nR0GUMqMgR47A9OlWp8myoUqpbUqpSUqpAnbu95qGhnPmxRBxJgaAq9elABFpu/fe5P1OAd5p9Q4AAX7e1zhKCpD0LF9uLu+5x5x4z4Jt22DpUjNxMIu9ywSYpR8tW0KPHmY+zpYtVieyXOfOZmHQG2/AtWuZe8zataaduyvmjmTQxPBLoCJQDwgHxth7Cju3eVxDw+OXjtNjYwF45E4A4j18QaBwro4dzYTUOXPM9/EJ5v+Lv593bKeWkhQg6Vm+3AxdVKqU5acYPx5y5oRBgxyYy5fVrWsuw8OtzeEGlDITUY8dM6NsmfHpp/Dyy645FZheE0OtdYTWOl5rnQB8jf29oryioeHm/yJI8L9Gg4RBTOs2jbCqYVZHEm6sYEGzQWFiARKXEAfICIhv0doUIK1aZevd+sMPzX+kggUdmM2XJS4fyOxHfi/XsqUZsh05Ek6cSP/Yixfh99+hd2/TzMxKSqmUJ7S7YXbWvtEGoLJSqrxSKgizw/YcV+RzpN9nmfePJ9t14qE6D5EzMKfFiYS7CwuDHTvg0KHkETMpQHzFxYtmXDsy0hQg2ZAvn1kyKRwkJMRceur6UycYM8ZMchw8OP2pMT/+aP7a+vVzXbZ0fKiU2q6U2ga0AoYBKKVKKKXmA2it44ChwCJgNzDT1uTQY8TFwdw5pgApVNi35y2JzAuzDZLNnZs8AnI9/jpfbPgi6Xtv4H0llSO89JJpA96+Pdx3n9VpREoyAnKTihXhvffMKqtvvoHHHrv5mLg4GD3aNJrNxlZGDqO1trsPrNb6JNApxffzgZuW6HqK+fMhMtIUINrHJ06LzKtYEWrUMKPn9e81IyBdZ3QFIGdgTgbUG2BdOAeSERB71q+HDh1MwyvZtMW9JI6A/PsvHD5saRR38vTT0K4dDBkC/9jpVzR2LBw8aPq5yVJw1/nmGyhUyFaA2J8/K4RdYWGwYgVcuJR6xOPQ+UMWJXI8KUBuFBdnVlnUrm11EmFPLtPMiS+/hPLlrc3iRvz9YcYMKFcO7r7btHNOtGSJaacSFmZWzgjXOHHC/Dt0vde8zcoIiLgVYWHm19E/61KP9u6ItDddyjNJAXKj//4zWxJKAeKecuY0LWXFTQoUMJ+YKlUyhcZdd5kJqnffDZUrm66pMvrhOpMnm2ZSXbslADICIm5N48amIfSRP9vz2/2/Jd2+6eQmC1M5lhQgN9phqy6lAHFffVNMH0hIsC6HGypeHNasMdvnXLpkBvOeeQZWr5aVWK6UkGB2D2jdGlae/wGFok7ROlbHEh7E3998kFg4P4DOlbpx5sUzvNbiNY5cPMKWU1usjucQUoDcKHF1Rd681uYQaXvqKXjFti/ZqVPWZnFDOXLAq6+azYP37YOPP5b/zq62bJmZovTYYzBz10y6VO1ClUJVrI4lPExYmFmUuXIlFMpZiHuq3ANAqymtuB7v+XswSAFyo5y2NfpXpV2y2/LzS/44v8l7hiOF9/j6a/NftFs3uBJ7hZJ5PLKDvLBY27amZ09iU7ImpZowvtN4LkRf8IrN6aQASWnCBJg501yXAsS9jRtnLj//3NocQtwgMtI0fOvXz/zyuBZ3jZCAEKtjCQ+UK5cpQubMSe7x07tWb/yUH0sOLrE2nANIAZLSsGHJBUhO6Vbo1u6/31zu3Gl2bhLCTUydanYoTuzHEhIQwoXoC5ZmEp4rLMx0RN1pa8FXIEcBGpZoyOIDi60N5gBSgKSUkGAaKhw7JpNQ3d2HH5pW+efPm92bZDKqcANam9Mvd9xhGkkBNCvdjJVHV1obTHisxKXzc1JsQtC+QnvWn1jPxeiL1oRyEClAUtLajHyUKmV1EpERpcxGKGPGmNmWu3dbnUgI/vkH9u6FgQOTb6tcsDKHzh/yqhbawnVKlDDdi+fOTb6tXcV2xOt4lh1aZl0wB5ACJCWtpVGCp2nb1lyuWmVtDiEwnU/z5k0+Qwiw/uR66hev75WbiQnXCAuDdeuSF/01LdWUAiEF+H3P75bmyi4pQBJpbYbxpQDxLJUqQbFiZt2jEBa6cAF+/hn69Elu2BsdF8264+u4q+xdlmYTni0szPyK+uMP832gfyBdq3Vl9t7ZHr0cVwqQRNOnmwKkcmWrk4hboRR06QILFpgOtkJYZPp0s0diys0A1x1fR0x8jBQgIltq14ayZVOfhmlZriWXYi5x4PwB64JlkxQgif76C4oUcZu9ysUt6NoVLl+GP/+0OonwUYmTT+vXhwYNkm/ffcbMTapfvL5FyYQ3UMqMgixZktwholrhagDsObPHwmTZIwVIoqgoyJ/fNLkSnqV1a8id2zRfEMIC//4LW7emHv0AOHn5JApFsdzFrAkmvEaXLmaELfFsc6EchQC4FHPJwlTZI79tExUtCsePyzC+JwoJMTuuzZ2b3K1HCBeaPNn8N+zdO/Xt4ZfDCc0VKhNQRbbddRfkyZN8GsZPef4uy1KAJLrtNlNe7vHc4Syf1q4dhIebn9DPPrM6jfAh0dFm/ke3bmYQNaXwqHCK5yluSS7hXYKCTMujuXPNdMXEAmTA7AHEJ8RbnC5rpABJtGKFWT9XrZrVSURW3HmnubxyxWz/OmwYnDhhbSbhE2bPNitgHn745vtOXD5B8dxSgAjH6NLFLMXduBFUihWbZ6+dtTBV1kkBAhAfb95FOnc2mzcIz1OtGowdm1xAfvKJaSgnbdqFk02eDKVLm6lIKcUnxLPnzB6qF65uTTDhdTp1An9/MwqSsrC9et0z9y6TAgTMJ+WzZ5M/RQvP9Oyz8O67qW+bPduSKMI3HD8OixdD//7mF0NK+8/tJzoumjpF61gTTnidggWheXPTlj3QP5CZPczeZVdir1icLGukAAGz0w9AhQrW5hDZ17Rp8nWlzFjltWvW5RF2KaVmKKW22L4OK6W2pHHcYaXUdttxG10cM0NTp5p5zwMG3HzftohtAFKACIcKC4Nt2+DwYcgZaDZNvXLdBwsQpVRPpdROpVSCUqphitvLKaWupXiDmZD9qE508KC5LF/e2hwi+0qUML8REr++/hpeftnqVOIGWusHtNb1tNb1gF+B39I5vJXt2IbpHONyWpvTL3feCRUr3nz/toht+Ct/qheRUzDCcbp0MZdz50KuINNy11dHQHYA3YG/7dx3IPENRms9OJuv41ybN5tN6MqWtTqJcKQ8eczlihXW5hBpUmYm3f3Aj1ZnuVWrVsH+/fYnnwJsO72NqoWrEhIQ4tpgwqtVrmymus2dC7kCTQHik3NAtNa7tdZ7HRXGMv/8A02aQGCg1UmEI61da5ZXR0VZnUSkrQUQobX+L437NbBYKbVJKTXIhbkyNHmy2fOlR4+b7/s3/F8W7V9E01JNb75TiGwKCzPNuxNibCMgvngKJgPllVKblVIrlFItnPg62XPpkmlh2Ly51UmEo9WoAffcY06WSoM5l1NKLVVK7bDzdW+Kw3qT/uhHM611A6AjMEQpZXemuFJqkFJqo1JqY2RkpAP/FPZduQIzZ5pdb3Pnvvn+91e+T76QfIxsM9LpWYTvCQuD69dh4xpTgHhqN9QMC5BMvoncKBwoo7WuDzwH/KCUypvG87v0jeMm335ruro0a+b61xbOV6WK+fdNnOcjXEZr3VZrXcvO12wApVQA5hTujHSe46Tt8jQwC2icxnETtdYNtdYNixQp4vg/zA1++cUMrD3yiP371x5fS9sKbSmSy/lZhO9p0gQKF4bVC0tRIKQAq4+ttjpSlmRYgGT0JpLGY2K01mdt1zcBB4AqaRzr0jeOVM6dgxdfNF00b1zEL7xDFdt/u337rM0h7GkL7NFaH7d3p1Iql1IqT+J1oD1m3pnlJk825+LtfW45cekEJy6f4PaSt7s+mPAJ/v5mcHfBPH+alW7O+hPrrY6UJU45BaOUKqKU8rddrwBUBtzvI+jSpaYJ2VtvQYDs1eCVqlc3PYy7doXXXrM6jUitFzecflFKlVBKzbd9WxT4Rym1FVgPzNNaL3RxxpscPGjmNQ8YYFZ632jdiXUAUoAIpwoLg/Pn4VDEGXZG7rQ6TpZkdxluN6XUcaApME8ptch2153ANtsbxy/AYK31uexFdYJFiyBfPmhsd1RXeIPcuWH4cHN90iRrs4hUtNYDtNYTbrjtpNa6k+36Qa11XdtXTa31e9YkTe2778ym2f362b9/3fF1BPoFUrdYXZfmEr6lfXvz2WrnpTWAZy7Fze4qmFla61Ja62CtdVGtdQfb7b/a3jDqaq0baK3nOiauA2kNCxea0y8y+uHdXn/dXLZta20O4fHi400B0q6d6fRvz/qT66lXrJ4svxVOlTu3mTkQfMZ8gN540u369GXIdzuh7toFJ0+abdyFd/PzMw3KgoKsTiI83J9/wrFjaff+iE+IZ+PJjXL6RbhEWBjEfLMIf+XPkoNLrI5zy3y3ANm+3VzeLm8UPiEmRjYaFNk2eTLkzw/3prEG8OD5g0TFRnFbidtcmkv4pi5dgOj8lFK3897K99BaWx3plvhuAXL0qLmU7qe+ITYWLl6EcePMslwhbtGFCzBrFvTpAyFpnF05FXUKgJJ5SroumPBZpUpB/foQecn0Afnf3/+zONGt8c0CJPGdpECB5HbdwrtdvgzTp8NTT8ESzxuqFNb76SeIjk779AvA6SunASiau6iLUglfFxYG15aY/a4OnD9gcZpb45sFyMcfw/r1MGaM1UmEFXa4RSsJ4WEmT4ZatUx3/7REXIkAIDRXqItSCV8XFgZ6Wx8CCKJgSEGr49wS3yxAfvvNbGGZ3kcZ4V02b4a8tma8a9dam0V4nF27zGeWhx+23/sDIDoump93/UxIQAiFcxZ2bUDhs+rXh5Il/Mh1rQZ7z3rW1my+V4CsXw87d9rfQUp4r3r1IDISevY0fbS3bbM6kfAgkyeb1foPPZT2MZ+t+4y/Dv/FhHsmEOAnS/uFayhlRkGijlRhzxkpQNyX1vDGG2Yae1pdhIT3Cgoyc0By5DDdb4XIhLg4mDbNtL4OTefMyrJDy6gdWpv+9fq7LpwQQPfuEB9RlSMXDhMT5zkbb/pWAfLtt6b76TvvyORTX9WiBfTta9rwx8ZanUZ4gEWLICIi/TO2cQlxrDq6ijvL2t2sVwinuusuyHmtKgkkeNREVN8pQLSGzz6DBg1g6FCr0wgrdexoVsWsWWN1EuEBpkwxO4927Jj2MetPrOfK9StSgAhLBAbCXXXLAbD/zBFrw9wC3ylA1q0zzccefzztWWTCN7RubU7oL1hgdRLh5s6fhzlzoHfv9Bvpztw5k2D/YDpU7OC6cEKk0KmNmWS/YN1/FifJPN8pQGbMMJ0we/WyOomwWt680Ly5FCAiQz//bJro9k9nWsf1+OvM2DmDe6rcQ76QfK4LJ0QKbVvkBmDCkWcsTpJ5vlGAaA1z55pPvolLMYVv69jRrIQ5edLqJMKNTZkCNWqYM7dpmbhpIqeiTjGw/kDXBRPiBoXzJs9rjI+3MMgt8I0CZM8eOHAAOne2OolwF4mbEC5caG0O4bb274fVq82CubTO2iboBN5b+R53lb2LuyvJxpbCOrmDciddX7nSwiC3wDcKkD/+MJdSgIhENWuayxMnrM0h3NbUqabwSK/3x/oT6wmPCmfQbYNQMrdMWCjYP3mzzV9/tTDILfDeAkRr03Bq4UJYvBhq14YyZaxOJdzNl19anUC4oYQE0/ujbVsomc6+cnP3zsVf+dOxUjpLZIRwAaUUAX4BFIpqwW+/ecaem97Zri8hwYx2LFhg+n3Ex8MDD1idSrij8HCzS24+mTwokv3zDxw+DP/LYHPRP/77g+ZlmlMgRwGX5BIiPY1LNuZicBA7T5qm302aWJ0ofd45AvLjj6b4GDjQ9Hu4etU0zBcikb8/VK9url+5Ym0W4XamTIHcuaFbt7SPiYiKYFvENjpV7uS6YEKkI0dADnLlv0ZgoGechvHOAmTqVKhYEb76yoyEdOgAgwZZnUq4m+HDzeXVq9bmEG7l6lWz/LZHD8iVK+3jVh1bBUCLMi1clEyI9OUIzMF1fY22bU0BorXVidLnnQVIkSJmAwc/P9NFaOFC0wNEiJRy5jSXMgLiNEqpnkqpnUqpBKVUwxvue1UptV8ptVcpZbeDl1KqoFJqiVLqP9ul0891/P67GThNr/cHwOpjqwn2D6ZB8XTW6ArhQjkCcnA59jLdu8OhQ7Bli9WJ0uedBci1a8m/XGRmukhL4sdbKUCcaQfQHfg75Y1KqRpAL6AmcDfwhVLK387jXwGWaa0rA8ts3zvV1KlmvvqdGXRVX3VsFY1KNiI4QD7cCPdwe8nb2X9uP8Ubr8bPz6zDcGfeWYCcPQuFClmdQri7xALk/Hlrc3gxrfVurbW9PcLvBX7SWsdorQ8B+4HGaRw3xXZ9CtDVKUFtTp6EJUvMfoV+6bw7Xrt+jU0nN3FHqTucGUeIWzK44WBCc4XyxtqnaNUumh9/dO/TMFKACN9Vt67Z4GPFCquT+KKSwLEU3x+33XajolrrcADbZai9J1NKDVJKbVRKbYyMjMxyqOnTzSK6fv3SP27jyY1cT7hOszLNsvxaQjharqBcDGsyjH/D/2Vni3ocOgRr11qdKm1SgAjflScPxMbCDz9YncSjKaWWKqV22Pm6N72H2bkty5/VtNYTtdYNtdYNixQpksXnMKtfmjSBKlXSP3b1sdUANC3VNEuvJYSzBPoFAnAqbi8hIe799uZ9fUC0lgJE3BrphpotWuu2WXjYcaB0iu9LAfY25olQShXXWocrpYoDp7OSMTM2b4adO+GLL9I+JiIqgvwh+dl1Zhel8paiSK6sFTtCOMuV68lz2rp0Mfuwjh1rNgB3N943AnLtmvlUK5vOicxo3txcnjplbQ7fMwfopZQKVkqVByoD69M4LnE9Sn9gtrMCTZ1qzsil1bNw7t65lP2kLE8teIrDFw5TPn95Z0URIsuea/ocCkXRXEXp0wciI2HZMqtT2ed9BUiOHKZ38vbtVicRnmDcOHM5YoS1ObyUUqqbUuo40BSYp5RaBKC13gnMBHYBC4EhWut422O+SbFk9wOgnVLqP6Cd7XuHu37dDFV36QIFC958/8ojK+k2oxsJOoFp26axK3IXZfOXdUYUIbIld1BuhjQaQmx8LB07Qv78Zm6TO/K+AiQ+Hs6dgxIlrE4iPEGdOuZy0iRrc3gprfUsrXUprXWw1rqo1rpDivve01pX1FpX1VovSHH7QK31Rtv1s1rrNlrryrbLc87IuXCh+aSYVu+PD1d/SOGchVnWbxnRcdGcuXqGHtV7OCOKENkWHBBMTHwMwcGmod6sWe7ZbcD7CpDdu81pmNtuszqJ8ARKma654Bm7NwmnmDoVCheGu++++b5D5w8xb988HmvwGC3KtqBthbbcV/0+wqqGuT6oEJkQ7B/M1eumw3PfvhAV5Z49QbJVgCilRiul9iiltimlZiml8qe4L8Muhw7366/mHSQgAO6Q9fkik15+2VweOWJtDmGJc+dMw+Q+fSAw8Ob7/zn6DxpNz5o9AVjSdwk/9/wZJU0OhZsKzWVWq3+69lNatIDKleHbby0OZUd2R0CWALW01nWAfcCrcEtdDh3nwgXo1QtCQ01fhwoVnPpywovUrGkud+60NoewxMyZZt56WqdfSuY17UkiryT3F5HiQ7izQbcNomy+svy862eUMvuyrlwJe+21BLRQtgoQrfVirXWc7du1mKV0kPkuh46zaJHZ/2XcOBn9ELemRg1zKQWIT5o61dSgaW2YnSvQdMyNjot2YSohsi5HYA5alW/FoQuHANNYLyDA/UZBHDkH5BEgcSJZZrscOs7cueYkbpMmTn0Z4YXy5zdt2ffvtzqJcLF9+2DNGvMGndagxqWYSwDkC8nnwmRCZE/5/OU5efkk0XHRFCtmVnhNmWJG+9xFhgVIZrocKqVGAHFA4mKfTHc5dFQLZZYuNfM//J17pkd4qStX4JtvrE4hXGzaNLPny0MPpX1MYgGSJyiPi1IJkX3l8pcD4OjFowA8+iicPm0+q7uLDAsQrXVbrXUtO1+zAZRS/YHOwINaJ217k9kuhw5pocz58xARYfb2EEKITEhIMKdf2rZNf9X+5djLAOQNluaGwnNUK1wNgG0R2wDo0AFKl06/06+rZXcVzN3Ay0CY1vpqirsy2+Uw+44cgcceM9erVXPKSwghvM/ff8PRo2lPPk2UNAISLCMgwnPUK1aPkICQpH2LAgJg6FD480/Yts3icDbZnQMyDsgDLFFKbVFKTYD0uxw63PjxZvktQNWqTnkJ4QPyyfl9XzN1KuTODV27pn/c5RgzAiKnYIQnCfIPolGJRqw6tirptoEDIWdO+PTTrD3nuXOObZeU3VUwlbTWpbXW9Wxfg1PcZ7fLocOdOWMuX345uaGUELfqvvtMC3/hE65cgZ9/hp49zRtyei7FXCLYP5jggGDXhBPCQZqVbsb6E+sZs3oMYLYZ6N/ftGaPiLj15xs4EBo7cD2r53dCDQ+Hhg3hgw/MbDIhsiJnTrh6NePjhFc4dQrq1cv49AuYOSBy+kV4ogfrPAjAC0teSLrt2WfN3kejR9/ac8XGmrUejmwy7vm/sU+elH1fRPblzQuXLsHZs1YnES5QsaJpzHTXXRkfeynmkpx+ER6pVmgt+tbpC0BMXAwAVaqYVV/jx5vP75m1YgVcvgydOjkun2cXIFqbAqR4cauTCE/Xp4/5//Tee1YnEW6mcM7ChEeFJ+2tIYQn6VTZVAx7zya3QX39dTMKMmpU5p9n5kwzZ6qDAzdW8ewCZN06MwdENp4T2VWzJrRrB2PHmn3ZhbCpWaQm0XHRRERl4aS5EBarXrg6AHP2zqHldy05f+08lSqZ048TJmRuC6zYWPjtN7j3XggJcVw2zy1A4uJM+ZYnD/TubXUa4Q0OHjSXUoCIFPaf20+gXyBl8pWxOooQtyyxf83ry19nxZEVLDqwCIA33zTTJp97LuPn+OUXswImvYZ9WeG5BcjEifD779CqlRkXEiK7xo83lwULWptDuJW9Z/dSqWAl/P2ky7LwPDkCc6T6PnEkr0wZeO01M7KRXndUrc2y3SpVoH17x2bz3AIkceJpZqaxC5EZ7dqZfavj4jI+VviMvWf3UrWw9BgSnil3UOoP6ImN9QCef940EB8wwDTls2f+fFi/HoYNc/xCU88tQA4cMJey861wpKAgiImxOoVwE6NXjWbPmT00LN7Q6ihCZEnuoNzM6zOPt+56Cz/ll7S1AEBwsOmHc/26aYV08WLqx166BE89ZUY/Hn3U8dk8swC5csV0UqlQAYoVszqN8CYFCpjxyFOnrE4iLPbm8jd5aelLPFDzAV5s9qLVcYTIsk6VO/FmyzcplKMQp6+cTnVf5crw/fewZQu0aZO8Kfi5c9Ctm5mkOnkyBAY6PpdnFiCvvw6bN5u2bEI40quvmo8Dy5dbncQrKKV6KqV2KqUSlFINU9zeTim1SSm13XbZOo3Hv6WUOmHb6mGLUsqBXQjSFh0XzQerPqBnjZ5M7z6dIP8gV7ysEE5Vv3h9FuxfQIJO3U89LAxmzYL//jNbqjVoAOXLm94f333nvBMNnleAaA2zZ0PnzuaXhRCO1KYN+PvDrl1WJ/EWO4DuwN833H4G6KK1rg30B6al8xxjU2z3MN9JOVPZdHITsfGxPFj7QZl8KrzGgLoDOH3lNF9u+PKm+zp3ht27za4mRYuabQo2bYK+fZ2XJ8B5T+0kp06ZJvYvypCocILgYKhUCXbutDqJV9Ba7wZQSt14++YU3+4EQpRSwVprt5iAk7iBV9PSTS1OIoTjtK9olrEMXTCUIY2H3HR/iRKu7cXoeQVI8eKmXXa8czbXFYKaNWHHDqtT+JL7gM3pFB9DlVL9gI3A81rr884OtOrYKioXrExorlBnv5QQLlMwR3KLAa31TR8MMuti9EUuRF+gbP6y2crjeadgwHxKzWgLSyGyqkED2LcveadlkS6l1FKl1A47X/dm4rE1gVHA42kc8iVQEagHhANj0nieQUqpjUqpjZGRkVn7g9horVl9bDXNyjTL1vMI4W5SFhyx8bFZfp7q46tT7tNy2c7jmQWIEM7Upo25LCOdLzNDa91Wa13Lztfs9B6nlCoFzAL6aa0PpPHcEVrreK11AvA1YHczcK31RK11Q611wyJFimTrz3Pu2jnOXD1DndA62XoeIdxRvuB8AHy7+dssP0e8dswZCClAhLjR7beby2vXICrK2ixeSimVH5gHvKq1XpXOcSl3muyGmdTqVEcvmo5M2R1eFsIdnXv5HCXzlGTef/Oy/ByNSzambtG62c4iBYgQN1LKrLQC2L7d2iweTinVTSl1HGgKzFNKLbLdNRSoBLyeYoltqO0x36RYsvuhbanuNqAVMMzZmRMLENn7RXgjP+VHrdBanLmatVPM32/7nnn75lGvWL1sZ/G8SahCuEKNGubyzTehdm0YY3fqgciA1noW5jTLjbe/C7ybxmMGprjuxEWAN1t3fB1LDi4BpAAR3it3UG6OXTp2y4/bHrGd/r/3p1X5VnxxzxfZziEFiBD25MplLpcsMV8jR5o27cKr9fy5J8cuHSMkIIQiObM3l0QId5UnOA+XYy5nfOANNpzcQIJOYGLnieQMzP5CEDkFI4Q9efOahmSJpDW719NacyrqFCXzlOSBmg9keYmiEO4uT1CeVJvSZdbB8wfxV/4Omx8lBYgQ9uTKBcuWQYsW5nuZjOr1omKjuJ5wnWduf4bvun5ndRwhnKZorqJcjLnItevXMv2Y+IR43lv5HqXzlSbAzzEnT6QAESItd91l9h0C2LrV2izC6c5eOwtAoZyFLE4ihHOVzFsSgPUn1mfq+KjYKBp93QgwI4WOIgWIEOlp2dJc9ulj9iESXuvsVVsBkkMKEOHdOlfpTLn85ejxcw/2n9uf4fEjlo1g8ymze8KihxZlcHTmSQEiRHpS7kE9d651OYTTJY6ApGxXLYQ3KpyzMIsfWsyF6AtM2jwpzePm/zefEctG8Nn6z5Juq1q4qsNySAEiRGbtz/iTgvBcSSMgcgpG+IDKhSpTPHdxpm2blnRa5e8jfzNh44SkY+754R7e/+d9AFY9sopDzxxyaAYpQITIyMiR5nLxYmtzCKc6cvEIYCboCeELyhcoz/FLx/lxx48A3PXdXTwx7wmOXzrOxE0TUx3btFRTyuUv59DXlwJEiIy88gp07w7rMzdhK5Vz5xyfRzjFH/v+oH6x+jICInzGjB4zCPAL4KtNX3HsYnJjsrAfw3j8j9T7QzpjWboUIEJkRu3acP48xMVl/jGrVkGhQvDLL87LJRwiIiqC1cdWc2/VDDfwFcJrFMtdjHdbvcvfR/5m1KpRSbcnTjh1NilAhMiMQrZPxRcuZP4xhw+by4sXHZ1GONiEjRPQaLpW62p1FCFcqnOVzgCM3zD+pvveuPMNAGoWqemU15YCRIjMKGhbGZHZUyo7d8J774GfH7Ru7bxcItt2nt7JWyveAqBO0TrWhhHCxSoXqpx0/d1WydszrRiwgm7VuwHQrVo3p7x2tgoQpdRopdQepdQ2pdQs2xbbKKXKKaWupdjlckIGTyWEe0scAdmcyaHJ//0PwsNh3jwoX955uUS2Jb4Bt6vQTtqvC58T5B9E01JNaV6mOcNbDE+6vXmZ5tQrVo/1A9fzdqu3nfLa2e2nugR4VWsdp5QaBbwKvGy774DWul42n18I93DnnVCnDgwcCLfdBpUqpX2s1rB0qemkevfdrssoskShCPQLpH6x+lZHEcISKx9eiVIKpRSj243m6vWr+CkzPtGoZCOnvW62RkC01ou11omz8tYCpbIfSQg3lDMnzJxp9oT57be0jztwwJx2OXsWGjRwXT6RZfvO7uN6wnVqF61tdRQhLOHv559UcLxwxwu8cdcbLnldR84BeQRYkOL78kqpzUqpFUqpFmk9SCk1SCm1USm1MTIy0oFxhHCwSpUgOBhefhmuXIGFC+GZZ8x9hw6ZUZKUIyNDhliTU9yS7ae3A1A7VAoQIVwpw1MwSqmlQDE7d43QWs+2HTMCiAOm2+4LB8porc8qpW4DfldK1dRa37T/r9Z6IjARoGHDhrLZhnBf/v7w8cemsJg6FZ580txeqxa89RacPJl87LhxyfNGhNuKS4jjl12/EOAXQLXC1ayOI4RPybAA0Vq3Te9+pVR/oDPQRtv6uWqtY4AY2/VNSqkDQBVgY7YTC2GlJ56ACROSiw+AQYNSH3P1KuTI4dpcIku+3PAlv+7+lTfufIPggGCr4wjhU7K7CuZuzKTTMK311RS3F1FK+duuVwAqAwez81pCuAWloHHj9I+R4sNjLD+8nIoFKjptlr8QIm3ZnQMyDsgDLLlhue2dwDal1FbgF2Cw1lp6Ugvv8MADULUqLFuWfNvp09blEVmitWbN8TU0Ld3U6ihC+KRsLcPVWttdi6i1/hX4NTvPLYTbatcO9uwx1xOX5xYpAnPnSs+PGyilegJvAdWBxlrrjbbbywG7gb22Q9dqrQfbeXxBYAZQDjgM3K+1Pu+IbKevnOZC9AWalGziiKcTQtyi7PYBEcK3rViRfL1zZ+tyuK8dQHfgKzv3ZaZX0CvAMq31B0qpV2zfv5zBYzKlaO6iXHzlInEJt7C/jxDCYaQVuxDCabTWu7XWezM+Mk33AlNs16cAXbMdKoUg/yByBuZ05FMKITJJChAhhFUy0yuoqNY6HMB2Geq6eEIIZ5JTMEKIbMlMryA7Mt0rKJMZBgGDAMqUKZOVpxBCuJgUIEKIbMmoV1Aaj8lsr6AIpVRxrXW4Uqo4YHe5kTQ0FMLzyCkYIYTL3UKvoDlAf9v1/kBaIypCCA8jBYgQwmmUUt2UUseBpsA8pdQi211p9gpSSn2jlGpoO+4DoJ1S6j+gne17IYQXkFMwQgin0VrPAmbZuT3NXkFa64Eprp8F2jgtoBDCMjICIoQQQgiXkwJECCGEEC4nBYgQQgghXE5p7T4r1pRSkcARoDBwxuI44B453CEDuEcOd8gA7pHDVRnKaq2LuOB1HCbF+4irucP/i8zwhJyekBEkZ2bZfR9xqwIkkVJqo9a6YcZHen8Od8jgLjncIYO75HCHDCI1T/k38YScnpARJGd2ySkYIYQQQricFCBCCCGEcDl3LUAmWh3Axh1yuEMGcI8c7pAB3COHO2QQqXnKv4kn5PSEjCA5s8Ut54AIIYQQwru56wiIEEIIIbyYWxUgSql6Sqm1SqktSqmNSqnGKe57VSm1Xym1VynVwck5nrK9zk6l1IdWZEjxmi8opbRSqrCrcyilRiul9iiltimlZiml8rs6Q4rXu9v2WvuVUq84+/Vsr1laKbVcKbXb9n/hGdvtBZVSS5RS/9kuC7ggi79SarNS6g+rMoi0Wflzmsl8bvOznBErftYzw53eDzKR1TPeL7TWbvMFLAY62q53Av6yXa8BbAWCgfLAAcDfSRlaAUuBYNv3oa7OkCJLaWARtt4oFvxdtAcCbNdHAaOs+LsA/G2vUQEIsr12DRf8fywONLBdzwPss/3ZPwResd3+SuLfi5OzPAf8APxh+97lGeQrzX8bS39OM5nRLX6WM5HTkp/1TGZzm/eDTGT1iPcLtxoBATSQ13Y9H3DSdv1e4CetdYzW+hCwH2hs5/GO8ATwgdY6BkBrfdqCDInGAi9h/l4SuSyH1nqx1jrO9u1aoJSrM9g0BvZrrQ9qrWOBn2wZnEprHa61/td2/TKwGyhpe+0ptsOmAF2dmUMpVQq4B/gmxc0uzSDSZenPaWa40c9yRiz5Wc8Md3k/yIgnvV+4WwHyLDBaKXUM+Ah41XZ7SeBYiuOO225zhipAC6XUOqXUCqVUIwsyoJQKA05orbfecJdLc6TwCLDAogxW/ZmTKKXKAfWBdUBRrXU4mDclINTJL/8J5hdcQorbXJ1B2OGGP6eZYeXPckbcLY9dFr8fZOQTPOT9IsDVL6iUWgoUs3PXCMy228O01r8qpe4HvgXaAsrO8VlevpNBhgCgANAEaATMVEpVcHSGTOQYjhk2velhjsyRXgat9WzbMSOAOGC6MzJkgqtfL/WLK5Ubs3X8s1rrS0rZi+O01+4MnNZab1JKtXTZC4sk7vBzmhke8rOcEXfLcxMr3w8y4mnvFy4vQLTWbdO6Tyk1FXjG9u3PJA8hHcecZ01UiuTTM47O8ATwmzYny9YrpRIwffQdmiG9HEqp2pjzsVtt/7lLAf/aJuW67O/ClqU/0BloY/s7wdEZMsHVr5dEKRWIebOZrrX+zXZzhFKquNY6XClVHDid9jNkWzMgTCnVCQgB8iqlvndxBp/mDj+n2cmZyE1+ljPibnlScYP3g4x41vuF1ZNQUn5hzqm1tF1vA2yyXa9J6olSB3HexMvBwDu261Uww4HKlRnsZDpM8uQ2V/5d3A3sAorccLtL/y4whfJB22slTkyr6YK/dwVMBT654fbRpJ7Q9aGL/h+0JHlSmSUZ5Cvdfx9Lfk4zmc0tfpYzkdOSn/VMZnOr94NM5HX79wuXj4Bk4DHgU6VUABANDALQWu9USs3E/ADFAUO01vFOyjAJmKSU2gHEAv21+VdzZYY0ufjvYhzmjWmJ7RPeWq31YBdnQGsdp5Qaillp4A9M0lrvdNbrpdAM6AtsV0ptsd02HPgAc2ruUeAo0NMFWW7kDhlEGlz9M5IJbvGznBELf9Yzw53fDzLilhmlE6oQQgghXM7dVsEIIYQQwgdIASKEEEIIl5MCRAghhBAuJwWIEEIIIVxOChAhhBBCuJwUIEIIIYRwOSlAhBBCCOFyUoAIIYQQwuX+DyM1TmjGivs1AAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 648x216 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(figsize=(9, 3))\n",
    "plt.subplots_adjust(top=1.3, wspace=0.2, hspace=0.5)\n",
    "\n",
    "#\n",
    "ax1 = plt.subplot2grid((1, 2), (0, 0))\n",
    "ax1.plot(original['x'], original['y'], 'r-')\n",
    "ax1.plot(reference['x'], reference['y'], 'b-')\n",
    "\n",
    "# \n",
    "ax2 = plt.subplot2grid((1, 2), (0, 1))\n",
    "ax2.plot(reference['x'], reference['y'], 'b-')\n",
    "ax2.plot(noised['x'], noised['y'], 'g-')\n",
    "plt.show()"
   ]
  }
 ],
 "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.10.12"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
