{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Interactive visualization of propagation of light in Jupyter notebook\n",
    "\n",
    "Matplotlib present useful visualization tools for showing the intensity profiles at a distance z in an interactive fashion in jupyter notebooks. Here, we show how to use it.\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In the first place, we well determine the intensity distribution XZ for an object. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2020-05-02T18:16:04.562760Z",
     "start_time": "2020-05-02T18:16:03.313741Z"
    }
   },
   "outputs": [],
   "source": [
    "from diffractio import degrees, mm, plt, um, np\n",
    "from diffractio.scalar_sources_X import Scalar_source_X\n",
    "from diffractio.scalar_masks_XZ import Scalar_mask_XZ\n",
    "\n",
    "from ipywidgets import interactive"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2020-05-02T18:16:04.570012Z",
     "start_time": "2020-05-02T18:16:04.565488Z"
    }
   },
   "outputs": [],
   "source": [
    "# Initial parameters\n",
    "x0 = np.linspace(-1050 * um, 1050 * um, 1024 * 2)\n",
    "z0 = np.linspace(-0.125 * um, 4 * mm, 1024 * 4)\n",
    "wavelength = 0.6238 * um * 10"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2020-05-02T18:16:04.577405Z",
     "start_time": "2020-05-02T18:16:04.572878Z"
    }
   },
   "outputs": [],
   "source": [
    "# Definition of source\n",
    "u0 = Scalar_source_X(x=x0, wavelength=wavelength)\n",
    "u0.gauss_beam(A=1, x0=0 * um, z0=0 * um, w0=2000 * um, theta=0 * degrees)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2020-05-02T18:16:35.160747Z",
     "start_time": "2020-05-02T18:16:33.920954Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAADjCAYAAAB+QFQfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/H5lhTAAAACXBIWXMAAAsTAAALEwEAmpwYAAAnoUlEQVR4nO3deXxU9b3/8dcnkw1CCAlBQDbZVMACsokFqbtotWgXi9aK1Stt1dvW/rwutde2Vltrr73aWvfSShepFRXqjgt4XUACIqtAgAgEkCVsgZBl5vP743wDk5CFwMycM5PP8/EYcuY7Z85558zMh5PvOfM9oqoYY4xJLWl+BzDGGBN7VtyNMSYFWXE3xpgUZMXdGGNSkBV3Y4xJQVbcjTEmBaX7HSCaiEwBLga2quoprq0A+CdwAlACXK6qO0VEgIeAi4D9wDWqutA9ZxLwU7fYe1T16ebWXVhYqL16nRDT38cYY1pi4cIF21W1UyyWJUE6z11ExgHlwNSo4n4/UKaq94nI7UC+qt4mIhcB/4lX3E8DHlLV09x/BkXACECBBcBwVd3Z1LqHDx+h788ritvvZowxzWmTIQtUdUQslhWoPXdVfVdETqjXPAE4000/DcwGbnPtU9X732muiHQQka5u3lmqWgYgIrOA8cAz8c7f2kUiSlU4QmV1hAPVYSqqwuyvDFNeVUNElbIDVUTczkReVgaZaSFyMkNkZYRomxmijbtlhIT0kPUYGnMsAlXcG9FZVTe76S1AZzfdDdgQNd9G19ZYu4mRSESpqAqzbW8lJWX7mVe6i6KSXawp2cnWTWXs27UH9u6Ayn1QXQnhalAFjRxaiKSBCKSFIJQBWTnQNo/MDvkUdC6gZ88OnNqngNE923Nyx/Z0zssiNzvdir4xRygZivtBqqoiErN+JBGZDEwG6NGzZ6wWm3LCEWX3/mqKt5bz0qptvLmwlNXLN1CzaS3s2wWRcMsXqu454RrvP4AD5bD7c6o2w5YV3v/iHwGPi3iFv7AXXfr1ZPSpx/P1IV0Y3CWPznnZZKZbsTemIclQ3D8Xka6qutl1u2x17aVAj6j5uru2Ug5149S2z25owar6BPAEeH3usY2d3Cqrw2zYUcHMlVuYNqeE1R+vhi2roaYqsUFUvcK/cRlbNi7jxdnwYloI2nei4ORBnD+uL/8xogf9OrejfZt0vOPsxphkKO4zgUnAfe7njKj2m0RkGt4B1d3uP4DXgV+JSL6b73zgjgRnTkpVNRFKtu3jzx+X8uyrKyhb+jGUl/kd63CRMOzaQtncLUyb+xbTMrJI63UKZ54ziJu/1Ich3fPIbZPhd0pjfBWo4i4iz+DtdReKyEbgZ3hF/VkRuQ74DLjczf4K3pkyxXinQn4HQFXLROSXwHw33921B1fN4VSVbXurmLliEw88t5wtRfOCWdCbUl1JpHgBbxcv4O0pWWT2G8pXLxnCj8b0pm/nHOunN61SoE6F9FNrOxWyJhxh9ZZyfjFrFa+/OA82fep1gaSSdgX0/dIZ3HP5KZzRpxM52YHalzHmMCl7KqSJv8rqMEs27OEH0z5mxazZybeX3hLlZax5eQZXvP4q7Yeczh1XD+eKoT3Ia2tdNib12Z67k+p77tU1ERZ8tpPv/bmIz+a84x2kbG3SQmQNGMkt14zm+lG9rMibwInlnrsVdydVi3s4oizfuIfv/X0hy199Ayr2+h3Jf2kh2gw6jd/cMIZvDO5OdmbI70TGAFbc4yIVi3tpWQX/9e/lvPr0TO9LRaau9Ey6fek8pkw+jeEn5BNKs9Mojb+sz9006UBVmKcXrOf2370JJZ/4HSe4aqoofetlLpj3IedcdTEPf20wXTpk+53KmJiw4p5CVJXVW8qZ+OiHrJv1WuK/cJSsyst467GpDHpzOI/ffh4TBh1Phn3z1SQ5eweniAPVYR76v7Wcdt3jrHt1phX2oxApXsD1NzzEJY99yOe7D/gdx5hjYnvuKaC0rIJvPjGXZTP+bUX9WFVVMO/pZxgyfx3/+ukFjOlbSJr1xZskZHvuSUxV+bB4B0P/81mWTZ9uhT2GKpfP5SvffZj/fm0lFVVHMTCaMT6z4p6kqmoiPDBnDRd9/1FqVs1v/gmm5fbu4JG7H+O8373LVuumMUnGumWS0J6Kaib9bSGz//ysN1yuiZ9wDcumT2fUZzt4464LOLFrrt+JjDkitueeZLbsOsBZ973D7Cf/ZoU9gXYXzea0H0xj7pod2HdDTDKw4p5ESrbt44s/fYW1r8xMvUG+kkHJJ1x40594YUkpkYhtfxNsVtyTxKrNezn91hfZOe9tv6O0btvXc92PnmTK/BLCVuBNgFlxTwKfbtrLmFue58DSD/2OYgD27uC/bn2KRz5YawXeBJYV94BbUbqHsbdMtzNigmb/bu76yVM8/uE666IxgWTFPcCKt5Qz7rYXCK9OrQHNUkbFXu786V+Ytmi9HWQ1gWPFPaA27azg7J+9Qs1K22MPtPIybrx1Ki8t22wF3gSKFfcA2rmvinN//Q57F77rdxRzJHZ/ztW3PcPHJbv8TmLMQVbcA6aiKsw3npzH5tmv+h3FtMTWdZz3kxdZt3Wf30mMAay4B0o4otz+8qcsmDbdzmNPQpHiBYy/723Kym2MH+M/K+4B8tKyTUx98B8QrvE7ijlKW997g6umFlFZbYONGX9ZcQ+I4i3lXHPXTNi/2+8o5lio8uHfpvPY3BI7wGp8ZcU9APZUVHPJA3Ng4zK/o5hYqK7k5/c8y5INe/xOYloxK+4+q+1n3zLnNb+jmFgqK2XCb95k5z7rfzf+sOLus4/WlvHMI8/ZAdQUtGv+bH48Y5kNUWB8YcXdRzv3VXH5/W/C3h1+RzHxoMqLT75A0bqdficxrVDSFHcRKRGRJSKySESKXFuBiMwSkdXuZ75rFxH5vYgUi8hiERnmb/rDRSLKf7++ivJF7/kdxcRTeRlXPjiHPRXVficxrUzSFHfnLFUdqqoj3P3bgbdUtT/wlrsPcCHQ390mA48mPGkzlpfu4e+PvmjdMa1A2Uez+fU7a+zsGZNQyVbc65sAPO2mnwYujWqfqp65QAcR6epDvgYdqA4z6al5sGuL31FMIkTCPPbwS6zbtt/vJKYVSabirsAbIrJARCa7ts6qutlNbwE6u+luwIao5250bYHw6orNrJ31ut8xTCJtK2HyPxZSE474ncS0EslU3Meq6jC8LpcbRWRc9IPq/c3bor97RWSyiBSJSNG27dtiGLVxeyqq+cHD79v1T1uhBc+/wgIbXMwkSNIUd1UtdT+3Ai8Ao4DPa7tb3M+tbvZSoEfU07u7tvrLfEJVR6jqiE6FneIZ/6A/F62n/JMPErIuEzAVe/nun+fb0AQmIZKiuItIjojk1k4D5wNLgZnAJDfbJGCGm54JXO3OmhkN7I7qvvHNzn1V3PvE+xCxD3dr9dk7b/L+Gjv11cRfUhR3vL7090TkE+Aj4GVVfQ24DzhPRFYD57r7AK8Aa4Fi4EnghsRHPtyf5q+n2i6+0bpVVfDDpxfY3ruJu3S/AxwJVV0LDGmgfQdwTgPtCtyYgGhHbPf+av7nL3Pt1EfDxv97h4WfjeL0fh39jmJSWLLsuSe9mcs3UbnC9toNULmfm6ctsjNnTFxZcU+AA1Vh7v7bIutrNwetfGsOq7eU+x3DpDAr7gmweONuthe973cMEyTlZdz71mr71qqJGyvucRaJKD979VOotG8nmrpenrmAHXZJPhMnVtzjbOueSua+9pHfMUwQbVjGrOLP/U5hUpQV9zibsWIzfL7G7xgmiCJh7n9+hR1YNXFhxT2OasIR/vjvlXb6o2lUydz5bCyr8DuGSUFW3ONoY1kFGz6y0x9NE3Zt4fnlvn952qQgK+5xNOPTLbDb+lRN055+o9i6ZkzMWXGPk5pwhKlvWl+7ad76oo/ZsttGCTWxZcU9TsrKq1i7YJnfMUwy2LmJD9Zv9zuFSTFW3ONkyebdsP0zv2OYZKDK0x9utC80mZiy4h4nf12wCcI1fscwSaLow2L2V9rwFCZ2rLjHQXVNhPfnr/c7hkkiVZ99yuZdB/yOYVKIFfc42FFexfaVK/2OYZLJ/t3MLbWLeJjYseIeByU79sHOTX7HMEnmxUV22qyJnRYXd3fJu1A8wqSKt0t2WH+7abHFS7fY+e4mZpot7iKSJiJXisjLIrIV+BTYLCLLReS3ItIv/jGTh6oyZ/k2v2OYJLRt3Qb2VNhOgYmNI9lzfwfoC9wBdFHVHqp6HDAWmAv8RkSuimPGpFIdVopXbfU7hklGOzawfa99mcnExpFcQ/VcVa2u36iqZcB0YLqIZMQ8WZLaX1lD2cYtfscwyahyH8Vl5ZzYNdfvJCYFNLvnHl3YGyviDRX/1qpsXzXsLPU7hklGqizastfvFCZFHMmeOwAi8hTwZRGpATYBi4HFqvqHeIVLRtv2VsIBuzamOTqfrN/tdwSTIo64uANnAN1VNSwi3YAhwOD4xEpeJXvK7UwZc9RKNu4hElHS0sTvKCbJtaS4zwM6AltVtRQoBV6JS6oktmq7fcvQHL2tW3YRtuJuYqAl57k/DswRkVtE5AwRyYtXqGS2crP1mZqjt2/PPqrsXHcTAy0p7n8DpuLt7d8AfCAiNmB5PXsr7NiyOXrVe/dSXWPF3Ry7lnTLbFTVX0c3iEhWjPMktUhE2VC6x+8YJpkdKKfSiruJgZbsuS8SkR9GN6iqfeOinhr7YJpjcWAvFVU29K85di0p7p2B74nIJhF5SUTuFZFvxCvYsRKR8SKyUkSKReR2v/MYY0wiHXG3jKpeDge7YgYBXwBOA/4Vn2hHzw1s9kfgPGAjMF9EZqrqcn+TGWNMYrSkzx042BWz0N2CahRQrKprAURkGjABsOJujGkVjrhbRkSWiMjfReQ2EblQRLqLyJ3xDHcMugEbou5vdG11iMhkESkSkaJt220kR2NM6mhJn/uXgCeBCmAisBS4KB6hEkVVn1DVEao6olNhJ7/jGGNMzLSkz70MmO1uiEh/4KdxSXXsSoEeUfe7uzZjjGkVWtItc2L0fVVdTXDHlpkP9BeR3iKSifeXxsxErDgUsq+Nm2OQnUubTLvQmTl2LTmg+riI9MXbA14MZANLRSRHVffFJd1RUtUaEbkJeB0IAVNUdVm815uWJvTsnkdJvFdkUld2O7LS7dLG5tg1W9xFRNRzlrvfE29EyKHu50IRQVVPimvSFlLVV/BhYLOcbLtuiTl6Gbm5ZFhxNzFwJHvu74jIdGCGqq5X1fXAehF5HW8Y4ElAUTxDJpOTurbjVb9DmKTVNrctmSEr7ubYHUlxHw9cCzwjIr2BXXhdMiHgDeB3qrooXgGTzYmd2vgdwSSxzl3zCdlwvyYGmi3uqnoAeAR4xF1mrxCoUNVdcc6WlHq3z4FQul2wwxyVnt3a21juJiZa9Pefqlar6mYr7I3r1C4LsnL8jmGS1JCedpkEExvWuRdjBe0yIf94v2OYZCTCqV1z/U5hUoQV9xjLyUonv3tXv2OYZJSVQ7/8dn6nMCmiJV9iGthA25mxDJMKMkJC/5M6+x3DJKOCbnRqb9e/MbHRkj33Z92gYSIibUTkD8Cvm31WKyMifGmgjVNjWq6wd09ys1s8UKsxDWpJcT8Nb7yWD/C+3r8JGBOPUMnurF4dvTNmjGmBwad0tS8wmZhpyTupGm9EyDZ457mvU1W7plwDehfm2EFV02KXDrXuPBM7LSnu8/GK+0i8b6ZeISKBuwpTEBS0y6TjiSc2P6MxtdrmMbp7R79TmBTSkuJ+nareFXWu+wQSNNJisslMT+OLI3v6HcMkkYyeJ3F8h2y/Y5gUcsTFXVUPGz9GVf8a2zip49vDu0GaDd1qjsyI0/vRNsveLyZ27OhNnAw+Pg8Kbe/dHAERJp3eHREbdsDEjhX3OOnYLpPew7/gdwyTDPKP54s9C/1OYVKMFfc4SQ+lcfW5ff2OYZJAj+Gn0tX6202MWXGPo0sHdIE8O73NNO3q8/uRbmO4mxizd1QcdctvQ7cRI/yOYYIsrzNfG2hjEZnYs+IeRxnpadxwyUlgB8pMI3qNHkWPjnaBFxN7Vtzj7LJBx8NxffyOYYJIhFsuO9m6ZExc2Lsqzjq3z+K08aP8jmGCqMcgLujfxe8UJkVZcY+ztDThFxcOgEz709vUNf6SERTmZvodw6QoK+4JMKRHHh1H2ACaJkq7Au4670T74pKJGyvuCZCdGeKnVw214QjMQf3PGkf/LnbVJRM/VtwT5LJB3cg6eaTfMUwQZLXld1cMtQOpJq7s3ZUgeW0z+H/XjLbTIg3HjzmTkSfk+x3DpDgr7gl03ciepPe3LzW1apltePCa4WRlWBedia/AF3cR+bmIlIrIIne7KOqxO0SkWERWisgFUe3jXVuxiNzuT/LDFbTL5Lbrx1jfeyvW40vncEY/GyTMxF/gi7vzv6o61N1eARCRgcBEYBAwHnhEREIiEgL+CFwIDMS7YtRAv4LXd/2oXuQMPt3vGMYPbXJ5/NqRZNteu0mAZCnuDZkATFPVSlVdBxQDo9ytWFXXqmoVMM3NGwh5bTN48KaxkJHldxSTYKdeeiEje1tfu0mMZCnuN4nIYhGZIiK1n45uwIaoeTa6tsbaDyMik0WkSESKtm3fFo/cDbp4QFdOOOf8hK3PBEBhT5781jA7Q8YkTCDeaSLypogsbeA2AXgU6AsMBTYDD8Rqvar6hKqOUNURnQo7xWqxzcrODDH1+tE2HHBrIcL1N15Cn+Ny/E5iWpF0vwMAqOq5RzKfiDwJvOTulgI9oh7u7tpooj0wBnVvzxXfu5Rn7n8CVP2OY+Kow8gzufPsfvZtVJNQgdhzb4qIRA92fRmw1E3PBCaKSJaI9Ab6Ax8B84H+ItJbRDLxDrrOTGTmI5GWJtwz/iRyBtuwBCktpwP/uPlM8tpm+J3EtDKBL+7A/SKyREQWA2cBNwOo6jLgWWA58Bpwo6qGVbUGuAl4HVgBPOvmDZyCdplMu+1caFfgdxQTDyJcfN1XGdXHXl+TeKLWJQDA8OEj9P15RQlfbzii3PDcEp79nyeteybF5I04k4W//QoF7WzkR3Nk2mTIAlWNyTcdk2HPPaWF0oT7Lx7AcWPt7JmUkn88L952rhV24xsr7gGQ1zaDl249C7oN8DuKiYX0TO78yeUM6ZXndxLTillxD4j+Xdrxp7snQJtcv6OYY3Tat77Gf47pbWfHGF9ZcQ+QCad048ofXmljzySxjl88l79fM9IGBjO+s+IeIKE04bcXD+DUb37NhgZOQtJ3GK//5Fw6Wj+7CQAr7gHTNiud5yaPpvO4C5qf2QRHpxN47VeX0rezXV3JBIMV9wAqaJfJWz85h5yhY/2OYo5E+05MuW+iDQpmAsWKe0B1K2jD23d/mZBd3CPYcjrw+99cw4RTutkBVBMoVtwD7MSuucy5/6tI32F+RzENaZPL3fd8hyuH9SQtzQq7CRYr7gE3qHt73vufr5PWb7jfUUy0Nrncefd1fP/03oSssJsAsuKeBAZ2b88HD3ydrIGj/Y5iANoV8Ov7rufmcX1tfHYTWPbOTBInHZ/LB7+9jLyRZ/odpXUr6Mbj/3s9k0fbHrsJNivuSaTPcTl8cO+X6XXBJXYevB96nsK/H76ebwzpbn3sJvCsuCeZ4/PbMOcnZzPm2ish3b4skyi5w8bx4R++xdj+hXZWjEkKVtyTUF7bDKZfN4rb7vk+tE/c5QFbpbQQJ1/6VT76zVc4+Xgb98ckDyvuSSorI8StZ/VjxqM32Jk08dKugGvv/C5v3jKOLh2y/U5jTItYcU9iaWnCuBM7sejhb3LShMsgFIhL4qaEjJNH8fxjN3H/xQPIybLtapKPFfcU0KNjW2bfeia3/+pGKOzpd5zklpHFiKsm8skfLuesk46zM2JM0rLiniKyM71umg+e+h49zrvY9uKPRp9TeeSPP+TlG79IV+uGMUnOinsKEREGdGvPvJ+fxy8f+AH0PMXvSMkhpwPj/uMqlj12JVec2pPMdPtYmORnu3cpqE1miJvG9GHCyVfz4xnLePOv/4byMr9jBU8onS5nnMeU745mZO98+7apSSlW3FNYj45tmXbNCJae05/JUxew6o1ZcKDc71j+EyFrwGncc8MYrhzag7Z2wNSkIHtXp7hQmjCkVwfeu+Ms5k0cwvenzGfje3NaZ5FPC5F50ghunjSa747uRX6OfQnMpC4r7q1ERnoaY/sXUnT3+Xy8fhQ//MciVr31Duzb5Xe0+Aulk/OF0fzX1SO5ZnhP8tpm+J3ImLgTVfU7QyAMHz5C359X5HeMhKkJR1i5uZy7XvuUt2fOhS2rIdXeCzkd6DX2DH458Quc1e842mXbvowJtjYZskBVY3KFHnu3t1LpoTQGdW/Pc9eNZOs3BjN96SYemr6MrQvnJffefHom6X2GcOklQ/nxGb3p17kdGXb2i2mFbM/daW177g2prA6zbtt+nizawPTXlrN72cewf7ffsZqXkQU9BnLG2YO4+cw+DOuRb10vJinFcs89EMVdRL4B/BwYAIxS1aKox+4ArgPCwA9U9XXXPh54CAgBT6nqfa69NzAN6AgsAL6tqlXNZbDiXteBqjDrd+znueWbeW5OCesWfQrbSqCm2U0ZfyKQW0jeSYM4+4x+XD+yBwO65pLXNsNGbDRJLRWL+wAgAjwO3FJb3EVkIPAMMAo4HngTONE9bRVwHrARmA9coarLReRZ4HlVnSYijwGfqOqjzWWw4t64mnCEXfurWfn5XmZ8uo3ZC0tZs2IDkc1roGIvRMLxDZDVFgq606nvCYwY2o3Lh3Xh1K75dM7LIjsjFN91G5NAKdfnrqorgIb2uiYA01S1ElgnIsV4hR6gWFXXuudNAyaIyArgbOBKN8/TeH8RNFvcTePSQ2kU5mZRmJvFmH6FRC46mX2VNWzbW8WaHeXM3bibBet2su6znWzdtIMDO3fBvp1Quc/b0w/XeAuK/k9ABCQN0kLeUAmZbaBtB9I7FJDfKZ+evfIZ3Lsjp/fKZVBhHp07ZNM+O936z405QoEo7k3oBsyNur/RtQFsqNd+Gl5XzC5VrWlgfhMjaWlCbpsMcttk0Oe4HM4b0BmAcESprolwoDrM/qowFVVh9leG2VNVTUSVHQeqiLi/FPMyM8gKhcjNTCc7I0TbrBBtM0NkZ4bISk8jlCbWxWLMMUhYcReRN4EuDTx0p6rOSFSOaCIyGZgM0KOnjaZ4rEJpQsgV6A45fqcxpnVLWHFX1XOP4mmlQI+o+91dG4207wA6iEi623uPnr+hTE8AT4DX534U+YwxJpCC3oE5E5goIlnuLJj+wEd4B1D7i0hvEckEJgIz1Ts6/A7wdff8SYAvfxUYY4yfgnK2zGXAH4BOwC5gkape4B67E7gWqAF+pKqvuvaLgAfxToWcoqr3uvY+eKdCFgAfA1e5A7LNZdgLrIzpL3ZsCoHtfoeIYnmaF7RMlqdpQcsDcJKqxuRivYEo7kEgIkWxOgUpFixP04KWB4KXyfI0LWh5ILaZgt4tY4wx5ihYcTfGmBRkxf2QJ/wOUI/laVrQ8kDwMlmepgUtD8Qwk/W5G2NMCrI9d2OMSUFW3I0xJgW1+uIuIuNFZKWIFIvI7Qlcb4mILBGRRSJSOwpmgYjMEpHV7me+axcR+b3LuFhEhsUowxQR2SoiS6PaWpxBRCa5+VeLyKQY5/m5iJS67bTIfb+h9rE7XJ6VInJBVHtMXlMR6SEi74jIchFZJiI/dO2+bKMm8viyjUQkW0Q+EpFPXJ5fuPbeIjLPLfuf7ouGuC8j/tO1zxORE5rLGcNMfxGRdVHbaKhrj/v72i0rJCIfi8hL7n78t5Gqttob3heg1gB9gEzgE2BggtZdAhTWa7sfuN1N3w78xk1fBLwKCDAamBejDOOAYcDSo82A92Wxte5nvpvOj2Gen+MNA11/3oHu9coCervXMRTL1xToCgxz07l4w0wP9GsbNZHHl23kfs92bjoDmOd+72eBia79MeD7bvoG4DE3PRH4Z1M5j/I1ayzTX4CvNzB/3N/Xbnk/Bv4BvOTux30btfY991G4oYPVu6DHNLxhhv0yAW+YYtzPS6Pap6pnLt74OV2PdWWq+i5QdowZLgBmqWqZqu4EZgHjY5inMQeHg1bVdUDtcNAxe01VdbOqLnTTe4EVeKOM+rKNmsjTmLhuI/d7lru7Ge6meMNuP+fa62+f2u32HHCOiEgTOVusiUyNifv7WkS6A18GnnL3hQRso9Ze3Ltx+NDBiRoiWIE3RGSBeKNTAnRW1c1uegvQ2U0nMmdLMyQi203uT+YptV0gic7j/jw+FW9P0PdtVC8P+LSNXHfDImArXgFcQ+PDbh9cr3t8N94w3THdPvUzqWrtNrrXbaP/FZGs+pnqrTuWmR4EbsW7IBE0PTR5zLZRay/ufhqrqsOAC4EbRWRc9IPq/S3m63mqQciAd6GVvsBQYDPwQKIDiEg7YDre2EZ7oh/zYxs1kMe3baSqYVUdijcC6yjg5EStuzH1M4nIKcAdeNlG4nW13JaILCJyMbBVVRckYn3RWntxb2pI4bhS1VL3cyvwAt4H4/Pa7hb3c6sPOVuaIa7ZVPVz92GNAE9y6E/RhOQRkQy8Qvp3VX3eNfu2jRrK4/c2chl24Y3Iejpu2O0Gln1wve7xPLxhuuPyHorKNN51aal6gwj+mcRtozHAV0SkBK/762y8az/Hfxsd7QGCVLjhjWe/Fu8ARe2BpUEJWG8OkBs1/QFef95vqXug7n43/WXqHvT5KIZZTqDuAcwWZcDbC1qHd9Ap300XxDBP16jpm/H6HQEGUfcA01q8A4Uxe03d7zoVeLBeuy/bqIk8vmwjvFFcO7jpNsD/ARcD/6LuwcIb3PSN1D1Y+GxTOY/yNWssU9eobfggcF8i39dumWdy6IBq3LdRzAtXst3wjpavwusrvDNB6+zjXqhPgGW168XrW3sLWI13MfCCqDfkH13GJcCIGOV4Bu/P+Gq8PrzrjiYD3pDMxe72nRjn+atb32K88f2jC9mdLs9K4MJYv6bAWLwul8XAIne7yK9t1EQeX7YRMBhvWO3FwFLgrqj390fud/0XkOXas939Yvd4n+ZyxjDT224bLQX+xqEzauL+vo5a3pkcKu5x30Y2/IAxxqSg1t7nbowxKcmKuzHGpCAr7sYYk4KsuBtjTAqy4m6MMSnIirsxxqQgK+7GGJOCrLgb0wIi0kZE5ohIKAbLyhSRd6O+hm5MzFhxN6ZlrgWeV9XwsS5IveF23wK+ecypjKnHirsxjoh8L+pKPetE5J0GZvsWMCPqObNF5GQ33VHcVaRE5F8i8rCIvCcin4nIWBH5q4isEpE/RS3vRbdMY2LKirsxjqo+pt5QsSPxxrb5XfTj7lJofVS1JKq5H94YLeCNa7LETX8BWKuqY4HHgT/hjek9EPhy1HjiS936jIkpK+7GHO4h4G1V/Xe99kJgV+0dEekFlKo31C54xX2xiGQDHfBGHwRvsK8/qTfsbA0QBqrAG3scqBKR3Pj8Kqa1suJuTBQRuQboBfyigYcr8EbtqzUEb/TBWsPd/UHAwqiiPwR3xSR3ybVNWnfEvizgQCzyG1PLirsxjogMB24BrooqzAepdy3NkNszB+/KR9nuuf3xrnO5BK9L5pOopw7m0H8Cdf5DEJGOwHZVrY7pL2NaPTsFy5hDbsK7SMM73jWJKVLV/6g3zxt446q/iVeoD4jIJ3gFezkwCa9L5iMA9x9BG/cfA9Qt9ABnAS/H45cxrZuN525MC4jIMOBmVf22iKwGhqnq3mNY3vN4V3Va1ezMxrSAFfeAOf+C8bp9+/YGH9OD/0Tdb+DazA29pPWfW3cZUVPa+NWe9bCJY12G1plP680Y/fvVX542sh2il3/YshtZ5uHLauD5ethv2cTz6qZVrTePa9QG5qm7rPpLqj9f3ZCHv8ba8HvmiJdx+PO9DA2/kfSwBTW8jobbtO42P2wd2kB4beDxqHnqzN/441qx7XVVHX/4L5XcrFsmYHZs3877c4sOvvkjeuiDoHjvZ0W9n9HTcPADFtFD0xr9PFcstKFlNvJYRA9fd1O5msscqbeeyME2Pfi82t8tgnqPuelDv4sSPrgcJeJ+qtJAu7f8cMQtP2pZkYOPufmjskeUus85mCO6Peo5bn2RyKHMtZkiqu52qO3QcxueX936D22PQ88PR+rOf/B5kUiD7apK5AifUzuf9zN6XvfaRiJ17te5tFtE69zn4Lojh8/bxPyqCpGwe5NEon5GGmmvnT/iXoxI3cebec6BRX8sTOynPDHsgKoxxqQgK+7GGJOCrLgbY0wKsuJujDEpyIq7McakICvuxhiTgqy4G2NMCrLibowxKciKuzHGpCAbfiBg3JV8gjD8ayHQ8DgIrSsDBCNHEDJAMHLEOsN2G37AJMIBVR3hdwgRKfI7RxAyBCVHEDIEJUcQMiQD65YxxpgUZMXdGGNSkBX34HnC7wBOEHIEIQMEI0cQMkAwcgQhQ+DZAVVjjElBtudujDEpyIp7gojIFBHZ6k51bOhxEZHfi0ixiCx2l3OrfWySiKx2t0k+5giLyCJ3mxnHDCeLyIciUikit9R7bLyIrHT5bj/aDDHIUSIiS9y2KIpjhm+512GJiHwgIkOiHkvktmgqR6K2xQSXYZGIFInI2KjHYvYZSRmHXR3FbnG5AeOAYcDSRh6/CHgVEGA0MM+1FwBr3c98N52f6BzusfIEbYvjgJHAvcAtUe0hYA3QB8gEPgEGJjqHe6wEKEzAtvhi7esNXBj1vkj0tmgwR4K3RTsOdSUPBj510zH9jKTKzfbcE0RV3wXKmphlAjBVPXOBDiLSFbgAmKWqZaq6E5gFHPUXLo4hR8w0l0FVt6rqfKC63kOjgGJVXauqVcA0lzfROWLmCDJ84F53gLlAdzed6G3RWI6YOYIM5eqqOZDDoauixvQzkiqsuAdHN2BD1P2Nrq2x9kTnAMh2fw7PFZFL45ihMYneFk1R4A0RWSAikxO0zuvw/qoCf7dFdA5I4LYQkctE5FPgZeBa1xyk90Vg2DdUTUv0UtVSEekDvC0iS1R1jd+hfDLWbYvjgFki8qnb84wLETkLr6iObW7eeGokR8K2haq+ALwgIuOAXwLnxmM9qcD23IOjFOgRdb+7a2usPdE5UNXan2uB2cCpcczRomyJFrUttgIv4HWTxIWIDAaeAiao6g7XnPBt0UiOhG6LqHW+C/QRkUIC9L4IEivuwTETuNqdrTIa2K2qm4HXgfNFJF9E8oHzXVtCc7j1ZwG4D9QYYHkcczRkPtBfRHqLSCYw0eVNKBHJEZHc2mm816TBMzxisK6ewPPAt1V1VdRDCd0WjeVI8LboJyLipocBWcAOEv8ZSQrWLZMgIvIMcCZQKCIbgZ8BGQCq+hjwCt6ZKsXAfuA77rEyEfkl3ocZ4G5VbeqAaFxyAAOAx0UkgrdTcJ+qHlVxby6DiHQBioD2QEREfoR3JsgeEbkJ74MbAqao6rKjyXAsOfBGJXzB1Zl04B+q+lo8MgB3AR2BR9z6alR1hKrWJHJbNJYD6EzitsXX8HY8qoEK4JvuAGtMPyOpwr6haowxKci6ZYwxJgVZcTfGmBRkxd0YY1KQFXdjjElBVtyNMSYFWXE3xpgUZMXdGGNSkBV3Y4xJQf8fmBd9ujXcgskAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# insert cylinder\n",
    "u1 = Scalar_mask_XZ(x=x0, z=z0, wavelength=wavelength, n_background=1)\n",
    "\n",
    "u1.cylinder(r0=(0, 1.25 * mm), radius=(1 * mm, 1 * mm), refractive_index=1.33, angle=0)\n",
    "u1.draw_refractive_index(scale=\"scaled\", colorbar_kind=\"horizontal\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2020-05-02T18:16:35.166046Z",
     "start_time": "2020-05-02T18:16:35.163302Z"
    }
   },
   "outputs": [],
   "source": [
    "# loading the incident field in the simulation\n",
    "u1.incident_field(u0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2020-05-02T18:16:38.689621Z",
     "start_time": "2020-05-02T18:16:35.169234Z"
    }
   },
   "outputs": [],
   "source": [
    "# WPM propagation\n",
    "u1.WPM(verbose=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2020-05-02T18:16:39.962992Z",
     "start_time": "2020-05-02T18:16:38.692400Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaQAAADhCAYAAACHrsFeAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/H5lhTAAAACXBIWXMAAAsTAAALEwEAmpwYAAB0NElEQVR4nO29e5wdRZn//66ZySQhAUIMEAhKuAQCRB0xgiIoLisXRVGQXZRFxAuuIvLFGELMxpBfyCsEDGxERQRBZBFW5A6L3ARBkUuQAOFqhAAJhBBCgIQkk5mp3x/VlfOcZ6q6+5w5M3Nmpj+vV7/OOdVV1U/3qa5PPZeqMtZaChQoUKBAgd5GQ28LUKBAgQIFCkBBSAUKFChQoE5QEFKBAgUKFKgLFIRUoECBAgXqAgUhFShQoECBukBBSAUKFChQoC7Q1NsCSBhjLgEOB1ZYayckaSOB/wXGAkuAf7PWvmmMMcB84DPAu8DXrLV/T8ocD/xXUu2Z1trLcly7iH8vUKBAgcqx0lq7dS0qqjcN6TfAoSrtdOAua+044K7kN8BhwLjkOBG4ADYR2AxgX2AfYIYxZqtul7xAgQIFBiZerFVFdUVI1tp7gVUq+QjAaziXAV8Q6b+1Dg8AI4wx2wGHAHdYa1dZa98E7qAzyRUoUKBAgTpDXZnsItjWWvtq8n05sG3yfQzwssi3NEmLpRfoA2hpaWHHHXekubmZ5uZmWltbN53TvwHWrFnD008/zZIlS3pY0gIFCtQafYGQNsFaa2vp6zHGnIgz9xXoIRxwwAGMHj06SC4777wzZ555Js49mB+tra1MnjyZFStWBEnsr3/9a0FYBQr0AfQFQnrNGLOdtfbVxCS3IklfBrxX5NshSVsGHKjS7wlVbK39FfArgAZjilX9agxPPuCIYezYscycOTOz3O23384VV1wBkKohjRw5kjPPPJOtttqK+fPnR+uThOXr+fOf/8zSpUuruq8CBQp0D0y9La5qjBkL3Cyi7M4B3rDWnmWMOR0Yaa09zRjzWeB7uCi7fYGfWmv3SYIaHgH2Tqr8O/Bha632TZWhyRjbAdTX0+hbGD16NAcccAAAe+65J9OnT++Ux5PD6tWrO2lIzc3NPPnkkzzzzDOsX78+9zXHjRtXRnyy3pgc69at47TTTmPlypW0t7dz7bXX5r7PAgUKlOERa+3EWlRUV4RkjLkSp92MAl7DRctdD/weeB8umuPfrLWrkrDvn+ECFt4FTrDWLkjq+Trwo6Ta2dbaS7Ou3WyMHQRsAApiyo+JEyey++67M3z48E3aikfIlPbwww/3uPlMamoAEyZMYNq0aWV55s2bx4IFCwD429/+VmhPBfocpKHbBtK6C7a/ElJvYogxdiSwDmgDNuL+1PbkfEfyWTwtR0I77bQTY8eOZc6cOWXnZs2axVNPPVX3vpsDDzyQrbfemqFDh/L973+flpaWTeek9nTrrbeyZs2a3hO0QAEcsaSRTChcOo2MZD+m82WRWQOl/hCgrSCk2mOoMXYHHBF5QmqnREpea9KfECerDvKhXv6BWAP08h1xxBHstttuZSRkreVHP/oRS5YsYdWqVdx1113dLmetMXz4cA477DCgs/a0cOFCfvrTn7J48WLuv//+3hKxwACFfyctcSLSmpEhfT6P75dCZTsCZY3Ir/s6A2woCKn2GGqMHUvpT5EkFCIi+cdocgp9hzBxadIK/Rt5iK1W/6IciYHz0XzqoIOYMGECP5w8eVP61KlTWbJkCU8//TSLFi2q0dXrAwcffDBbbrklLS0tTJkyZVP69OnTefbZZ7nxxhtpa2vrRQkL9Hdkmdo8mXjyiJGMhCQcXTamFTWIdNnvyWu9WxBS7THUGLsLca0ni0xC5/PWlVZ37Hp5tLG0fzaL5D7Y0sKH992XM848kxGJX+jhhx/mZz/7Gc8+++wmf0t/x+GHH864ceOYO3cuDQ3u9bzmmmu45ZZbuPHGG1m9enXvClig36ESMpLvcWNKGUlGMl+sH5DEpvspTXoFIXUDPCFpVNup532qWUQSIj2dp1rSCsm4V0sLHz/oIH581lmQzAeaO2cOjz/+ODffdBPr168fkP60CRMmMGHCBL74xS9y1FFHAbBkyRJmz57NlVdemTsqsECBNOQhI08oebQiKNd+tEYVu4Yko7TrWGBdQUi1x1Bj7K6B9Hp4OnqkkoW8Zj+Zb8zYsRx01FFMEUT04ylTeHzhQv50112b/Gkh31klsvV1NDU18fnPf56PfOQjTE5MmAsWLOAXv/gFv/vd7wpTXoGqUEsTndec5DvZSNjsFrqGhzbxheouTHbdhBgh1SNCDsu8ixLqxvue0aP5xFFH8cOzz6Z58GAAzpo2jT/ddhuPLVxIO3Qio5g/DQYOMQEcdthhfOITn9hETNdddx2XXXYZN998cy9LVqAvoRqtyKdnkVEerShERFq+tHrXFIRUe/QlQspCHsIywL8cfTTTL7iA4SNGAHDujBnce+edLHzwwTItqB0XeRiKNtRpHgOpVR155JF8+ctf5otf/CIAc+bM4ZprrmHhwoW9K1iBuketTHT+/W4XaZX6iiAc+CBN9KFgiMJk1w3YzBg7rreFqAJ5Q8slxo4fz6H//u/8Z7KCwS/mzOHe22/nofvuAzqr9r6FtAeOLK1poLSupqYmvvKVrzB16lTGjXMtadq0acyfP7/wLxUIIouMGukcuKAJJBYp1xWtKEQ6odBvj8Jk1w3oq4RUKb586qlMOussGhobWfP22/zXySdz/RVXlIV3SoQCI7xm1E76fC0C3/s7RowYwTe+8Q3OOussjDE8+uijnH766X1yflaB7kMaGcWi6PKQUbVakZQrRDohgvPa21sFIdUe/Z2Qth87ln87+WS+esopAJw9bRrXX3YZy5cvDzY+32B1QIWO6AtpTWk+J11Hf0VLSwvHHnssp556KgAzZsxg7ty5RdBDgVxklJYmze/eRBcLeEgjtK4SkU9/pyCk2qM/E9JnvvpVfnzBBTQPHsyt11/PJeeey6Ni1QFJNrEGHJtIp31N/mgjbtbz+WX5/opvf/vbnHvuuQwePJjrr7+eM844o99NJC6QD3lMdBrVaEUxEglF68bMe3mIyNdTBDV0AzYzxu7e20LUGIOHDOFrU6fyrWQpnPNmzODCZJQeW4oEOhOUNOfpZUkksehgh9gx0Mx5Y8eO5eSTT+aUU06ho6ODb37zm/z2t7/tbbEK9CC6SkZZWpFEiIhioeCh/LqeWICDRy0JqS/sh1SgCuw4fjzTL7uMlg9/mDdWrGD6iSdyZxKOHFLNJWRDbaBENvLwDbRJ5JWfjZS0ppDPqV3U09+1piVLljBp0iRWrlzJrFmzuOSSS9hxxx2ZM2dOYcIbAKjWRGfV73bCWpGEXrnBv+sy+i7mn5J1+HezkbBW1l0oNKQE/UlD+tSRRzL78stpHjyYS887j4vmzmXlypWZoyk5Igoh1CB9447VLY+QttRGOXFpmfpb6zz88MO54oorGDZsGNdffz0//OEP63ZF9AJdRzVkJMs2Um556IpW1BUigs6mPo/Ch9QN6C+E9O+nnMLkefMAmD9jBhfMnl1xHWnROVoT0iGhoflOvoVJ8sky6fXneU077LADp556Kqeccgrr1q3jM5/5DPclIfcF+g9828+zSrdP8+WklpTXRJdFRDIvdNZ8fL6QKU/WBeX3VRBSN2AzY+z43haii/jWjBmcOH06bW1tnP61r3HLVVd1yhP7t0MvTZ6WoQkKSsSU5qMKkdNG9VuGklciU1/BjBkzmJ78X8cddxxXX311b4tUoEZIIyNtbgtpRVBu1k6zWvg6QgEHebQinxYiLn0PegAKtQ37LnxI/QBNTU385+zZHD9pEm+89honHXkkTzz4YDBvbJ5CyBSXFfQgTQq60etgCA+f5n1MkpyaKZGSjtKToeQhefsiZs6cycqVK5k/fz5XXnkl22+/PfPnz+9tsQrUACEyyjLR6fXm8pjo8pjndF4bKS/r9dDySjKS1o9aoSCkPo7BQ4bw49/+loOPPJKnFi7k/33pSyxLfBJZ69vp82mBDqFRlYcnOb1CcKzRo9IlOQ1KvvuNEkOBEP1pXtPPf/5zVqxYweWXX868efNoampiXmJyLdA3UY2JLqQVZQUuhIIW0hZgleY/HcGXZ8UXKa+sq5bI6rMK1DE2HzGC6ZddxsFHHsml8+ZxwoEHbiIjk3Jo+PRGcTSIAzoTjD4v0xtxxCIdozrk21/XfzbiRkeDk2NYcmwGDE++D03ONSd5ff1p2lxfwNVXX82nP/1p1q5dy9y5c5k0aVJvi1SgSqSRkYyak5YCb2XI4yvylgdJRp7AtFbk8/kIvdAAMQ8ZaRO8rysWzNQVFD6kBH3Nh/Se0aM59+ab2bOlhUvOOYf/njq17LxuZNopWQl0eLasJ9R6QiMvPTEWSi+Lz+vL6sm2lnJTntSeYoEQfbFV77vvvtx+++0MGzaMadOmMXfu3N4WqUAFiJGRPB8LXIj5eDxCWlHMT6QnweY1z4Vk1+95yBS/ughqqD2GGWP36G0hcmL4iBH89M472aulhUvPOYd5U6cGyUK/DNB56Z8QQr6f2FyhNIILOUB9HbFVIbTGIyfc6nDxjeIztGyRlLOvoKWlhZtuuontttuOKVOmFOa7PoKs8O20cO5YGU0mOmhBk0tPEFGoTyiCGgYwGpuaOO2Xv2SvlhYunjOHn06fHm1AoYannaghxFZqCJVJI7iYedC/jNJ3JK/n5ZSmRHkN72MKEVPWkkX1joULF/LpT3+aO+64Y9Pad0WgQ30jDxlpE1qMjDShhIIW5DuRphVpP5HMH0IaEcl3sDvfpcKH1Mfwn7Nnc8iXvsQlc+fys4SMZEPyjVf6WGLw2oeG9Dc1qt9ZeX3j1f4oKDfB+TJN4pD1+7zyk6SuZpw/yfuXvI/Jf98sOT8kqdf7s/qKr+mZZ57hc5/7HGvXrmXevHmcdNJJvS1SgQBC7Slkleig5COFzr5UTUbSR6OJR5rsQiQny2o/UWziu65DkqU0q+v+ojvep8Jkl6AvmOz+c9YsTpg6lduuvZYffeUrbBTLzhj1CeGgBKmNIM7pRhqqTyLNZJenReny2hQQCkP3n5KE9aTbjepIM+nVe8s/4IADuOOOO2hqauL444/niiuu6G2RCiRI8xfpczJwoRqtSOav1k8U0pBiJrjQPcTe9wZgVeFDqj3qnZCOnzqV786axf9dfTUzTziBDWrTN9mJxzQZHRWj5w35euSnLC8/s1TrtFneBM7p62jHrc4j70WbNCQ5tSZHWjBESI56wWGHHcaNN96ItZbDDz+c22+/vbdFGvDIoxX59Ni8P00GPsIuFD2niUmbzzSJhZCHiCrRduT1Xi8IqfYYZozds7eFiODz3/kOPzr/fN5+802O/tCHeG3p0k3nQpPndEeu7cIhcpLQduasTjumlflzaX4tLUNIW/PpscANHRDhZfXmvjZKxNRKWHOqZ63pq1/9KpdccgkbNmzg4x//eLE1ei9Ctq9QsIKHNpfHCCHk75H5ZNm0cmkDxFB5X4eWOwuh6xWE1A2oV0L6yEEHcf5tt2E7OvjPgw/m4XvuCeZLCxvNO88gD2JkFavXf6ZpR5VG/cSITc5L0p2Bn2CryUlrTnKSYT29GZMmTWLu3LksXryYj370o6xevbq3RRpQ0B25HBCFTHSxgZbUeNoj5bOISMuUNbCsNRHJoAxLbU12RVBDHWPk6NGckeybM+Pb346SEZQHB0g0qHMN6gDXsNoombo8TODQAQtNxOuXdWsfDqK+QZQCD9KCIdqTMoNwE2QHiWtCOel4cvGBE36y7XBgS2AEMDL53JJSMMQQXNCEl6FegiDmzZvHrFmz2HXXXTn55JN7W5wBBd0GYmTUgGuTejCkyci35SwyksEG0JmMYgPJUMCDLy8DivLA35MPkoLy97gIaugm1JuG1NjUxPTLL+fQo4/mFzNmcHGyanfe8M0YYk5MWaccWcUctaHrhGRKCwvXdYfk1ia/0Eup69Z5JLnoIA/vZ2oDNlDud/LkVi/h46NGjeKxxx5j22235aSTTuLCCy/sRWkGBrSJzkNqCjLCVCIUGZdXK4qFWYfMc1JjkbLJ8pUSh7Q2+Dpk+5fa0orCZFd71BshfXnSJE6ZO5dbr76amccd16WN3EIjLJ+uHaaovDpKJ82hmxdp5BGqVwc66HIyT8wM6F8eT0wyzNyP+KQpbwOdiakeAiEOOOAA/vSnP2Gt5YADDuDByCK6BbqOEBnpTjqNjGRQg1wWSCLNPOev4ctVSkRZCEXxxXyxNnDOoyCkbsAwY+yE3hYiwT4HH8x//9//0dbWxr+3tPDSM88AXe8A0xprHqdoHvLIcy2NSu4rpk2FgiHkiE6/tCHTIJTMfpKY1tPZ59SbgRAnnXQS8+fPZ+HChRx44IGsWbOmhyXo/4hZBKTGrrVuCGtFIZ9NjIjylOkKEYUsCmlEpPMgznvUMqih8CHVGYaPGMEPzz8fgNnf+x4vPvNMJ9tvtUcI2oQV8iFpf5FerSEUZJAXfhQpzWp57kFfMzR69JNiB+H8QoMovcw+uMGb6by/aRDOl7QlsBXwHmBrYBTO5+T9TUNFnV7b6il/089//nPmzp1LS0sLp556ag9ddWAg9j+GTHQhMpJ+JD0JXNbl80Op7RMpIyNJJRlpH5GWV373h54g699r7R+S9Uv59PSS2Gov1aLQkBLUi4b0H1Om8N3Zs7l41iwunDkzV5mujCpis7fzmurSNIRamPcqhTbd+ZdYdxx6Rrp/4aUD15fxAR9SY1qffPemvZ7WmkaNGsXjjz/ONttsw+c+9zluvfXWbr5i/0fMX+Qntfo8vhMOEYvUcPKa52JlQlpJLERc1h+CHrzpYIeYxhSzjshrLR+IJjtjzBLgHZKpJdbaicaYkcD/AmOBJcC/WWvfNMYYYD7wGeBd4GvW2r+n1V8PhLTf4Yfzk+uvp3XDBr601168mmwlkRd5Ovy0UXwsHDRtBneewIZY3lpoFDH/mD+n02P2fk0iUiuUkVPe1+QJaV1ySLNeT81tOvbYY7nssssK010NILUfD91OQqvT5/EVhbSoNCKCctJJe/9CckroCbUxs5yUyZfLCqbwqCUh9TWT3aestS3i5k8H7rLWjgPuSn4DHAaMS44TgQt6XNIKseWoUUxNoqbO/M53eG3JkjJVO8+RhSwfUsinAp3Vdglt6tP1SNk0WYUCE6o1Q2q/kAmkG8Ih7j703JvgZOfiSWdD8rsJZ7IbCWwDbJcco3GmvZHAFrgQc1lfd5jzrrjiik2mu+985zvdcIWBgzQy0iY6afrybctH0KWZ52R+n0Ygv15XzufR+ULvvc8j30m5P5mvV5rl5N5l/j60RqifR5YroFr0NULSOAK4LPl+GfAFkf5b6/AAMMIYs11aRV31z3T1OOz443nPttvyq5kzuS2Ze1QrGT3SRuohYgk1eG1XDuULmctihCZfOvldy5CHjLXZQc7Lkr9lun4BQ/OiLKVAh/XJd3CEMxznX9oeGJMc2+HIaQTl/iZvCqzlS3zppZfS1tbGnDlzmDixJoPUAQmpTevON6Q9+I7c//Z5PTRxSd9SyIpQKRFp2X0e7f/xc+okEUE5Sfl70LvVSo2quwhIoy8RkgVuN8Y8Yow5MUnb1lr7avJ9ObBt8n0M8LIouzRJq0vsue++fH/uXDra27nj6quj+SydO+0YZGcdIpBQvjRNSzZI/cJkNVrppNWTZ2P55YoMmmhCh56Yq59F6F4k8fjr6pGldPp67UpG3tnk3GBcwMPWOHLaAXhf8n0bXIDEcMon9NaCnBYvXszkyZMBOProo7tY28CFDGv20IOpkFakoQdIofI6f2jyq3wv0rShmMYjAxU0EYUmnfv6QiQUk6E7TNF9aT+k/a21y4wx2wB3GGOekSettdYYU9EzSojtRHAjid7CRw8+GIBzTzuNl555pkudVNoIQ9YrG2ssnyX9RUGd09fQZhBts5afelKuRuzll5AvjzaLpNUZs9fryMYGkd5BKaDBdwDNyTEEF6nnTX7rcY7MtcnxLuXBEF31NV1wwQUccMABTJo0iX/84x9cfPHFVdY0sOA7c01EUE5G0HlCa1p7DAU5xObFST+Sv66uT8vlod+ntFUdZISplkXef9r1Q/ORiig7wBhzBrAG+BZwoLX21cQkd4+1dndjzIXJ9yuT/M/6fLE6h/dSUMPEgw5i/m23sfLVVzn2Qx/irZUra1JvXmLqKtLqivmOYvlCMuctG5KlmmeQZuL0L6L2C8jOZBDlyw/5TsCb+zwxrUm+y4AIHwwRunYWDjzwQO68807efPNN3v/+97N8+fIKaxhYyCIjSSaQ3vHK8mlRcNDzRBQildB15GBO1pU2X8k/wxcHWlCDMWaYMWZz/x04GFgE3Agcn2Q7Hrgh+X4j8FXj8FHgrTQy6i00NjXxuROd9fHnP/5xlIyyTGx5zW6Q3hF3xT8Vkk/az/PIGpoHFVrfLnafuhMPjQRDcsfS9bW1jDIYwvuI2nGE8xaOdDYkeTej5GsaC+wC7ETJrLc1zuTngyEqNefdc889zJgxg6222oovf/nLFZQcmIiRkf9PffuVWpJGSBMPReLJukMBDdosJiG1EGma8/4hPwdOzlUMBSGF3kUC+bUJ0MvlrQBy7Uhp7qsV+oSGZIzZGbgu+dkE/M5aO9sY8x7g97j3+kVc2PeqJOz7Z8ChuIHoCdbaBWnX6A0N6V+OOYZZ//M/3HHNNZxx7LHYGiwPVCtk1ZdnJFPtaCc0igu95JVoTx4m8j2tnJYhpjnpkaUMlGjG+Y985B24l9lrTWuAt5PPNTjNya8M0a7qi2HixIk88MADdHR0sPfee7No0aKMEgMLvrP137WGoLWiLE0HwhqRtgyk1SPza1l9uiQHTTQhX65MD0Hnlfm1zFr7CoWK11JD6hM+JGvt88AHA+lvAAcF0i1Q8b7Pte7U09A8ZAhHJas233b55bQLMqq0I88rd3f5pmL5Q3bzELRcOnpJwncA1ciUl4Bi56RPKaTZyZe4SZTxxLOe0krlg3GmvS2Toy05vwY32c5rWJKcsvxNCxYsYPLkyZxzzjl8/OMfLwhJITTfB7LJKJQemqTaVSKSRKEDFUJEpAnDp8cg65cmw9CUDVln2rSOWqNPEFJ/xCeOOoqWffflN3PmcP/NN2c23log9kLmQcz2nKczD70waXJI00KaLCEbeN6XM62+LM1KygidTTqhKD0ohdV635Enp6E4chqKi8bzc5/eoUROq0U5PflWYsECZwg499xzue2221hS4eTq/oyQlq0j4EJEJNNlhGmjyAeVEVFIW/Npfo6cN5VprSdGRNKvI+9ZkpCXM0R2mqhkPdJUV4mVolL0CR9Sf0PzkCF8IZnI+NSCkiXRBo5KUYnPJ68/Ks1vEwrN1vKE5NPI4wPzZeULJF9C/dx0fl0mhpAvKySrh1Xp/vCjXE8cPjzc77sEjmRWAStxpLMhObc1zsc0Hng/8AFgjyRtNI64Qv6m++67jxkzZjB48GAOOeSQlLscWIiRkaEzkUgfqCwT6ui16S7N9yd9TtrfAyXzLJSvkSgJMKTBQPk74H/LuUb+fuRcI70WoyRDL29orzN/v7UmIyg0pF7BPoccQstHP8pv5szhrzfcUBPNqNI68o5EunPEklW3f0li95bmD5LltJmv0hFezO+grxvq2HQ9bZQ6i+bkuyctHxLejNOWhiWffvLtBpwJz2tMb1Iy7a2j1KHde++9AJx55pnccMMNAz7iTms82vwlNfmYeQ6VxyMUhRfTmHTgAZQGLTqcGpU37f0OaVpSBnlebrsiy8rry3dE3n8o4q7WKDSkXsAeyYz6fzz+eKo2lKXh5BntQ7aWk1YuC5VoZHnqluWyTIz6mcn7G0S++9WRdKHJw/J60sGsIcvpvPqZ+G0u2pNr+rlLzUnaO8BruNneb+LIZhguIm93nMb0YWBvYAKwM05z2hJ45L77mJlE3B144IEpd9//oQkkRkZSi9D/Y6guacrT0Z/6nK/HE4IMatETViVC2pCWTdcr5035uv05HZEn26ZcGV/KoaPttHmv1gRSaEg9jI8edhjHT53Kwgce4P6bbso98kF8zxrhd7WR1LKR5a0rRLah9DTIFyfkZwpBP+NK712Wj/kQYpqTNJP4dO9zaqDUqbyNIyhPWsOBzXFmu+1w2pHXmt5IjlceewyALxxxBH/4wx+6tMFjX4XUeELmtAaVx6fF6vFIC2ZoEOe0RuSvs1HUFSIhf00pp7y2JCHtF/Jl9cBK+oWMyiMHgDJYQ8uXZh2oFQpC6mGM3nlnAK674AJa16+P5svycXh0t+mtK42ukmvqkahOT7PLp0F3NFn3U2n+GGJE5zUnmS47vHZxvglHQt6stxHnb3orSR+eHCNwi7p2UJoDtfqPf+Suq6/mS0cfzY1XX82V115b9UoQfRHSXBuKCJPaS0fgvIecmyPTfB06L2QTUahj14MYSWp6hQQ9v0ib9uR9yXlC2mclZdV5YnLK9IKQ+jhGjx3LKeecA8DKl16qyR+aJ7S6Eo0B8je0niK5asgodq5SjShv/iwSlZ1eKN2TkOyY/Krk3rnsAxg8Oa3EaUNDcSuMb4Ejp1FAW1sbr99wAxx9NAfssgu34chsIJCSfL5aC5BpocFHyF+Y5kOS19MBBf4anhSyiEhrK+0qPc0vFCIpeY2Qz0ebxNM0IF8+1M5riYKQBLrboTbugx9kUHMzF0ybxmP33ZdbprzzeWKoFRnVksy6Wgaq0xT1tSp11KZpUNW0n7zmPU9OhtLqENLk9xqOnPwir1sC6xcswLa18a05c1hy993874IFvIILG++PiJnjQh1+aIAgtZBQ3VJr0eX87xARheoJ1evrkoESMU0pi4RCRKwRWoFCP6sYEXbXwKa7++ACAuM+8AEAlj3/fK78DeJTHlllQgdkh2hDekPLWz5PY9VOfo3QPYRCUEMhqbEOJXQt/YxCQSYyjzb/5AlGSQvXDc350OY9TVIbKW0M6Pdo2hxHRm3ACuAfwEOLF3NNshL4x3baiaOAliRvd5hbehOVkJHUGKSvJ2Ta02auUJScJwgZrCBl0j4bKPkK/TwjGYTgBxx+fpSel+TPQfn76MO8Y3twxfoRL5venkWb/kKaYq01mkJD6iF86MAD+fqMGaxfu5Yn//a3YJ5KRgeVagSh0VAlqES2rM6uKxqRJJ3YCDULecxuHiGfTyx/SGMzpE+iTYPWnKQz3XcIGyn5Jnwk1eDkdzuwIgn5Nrvuyh7Ae4G/An/DRfD1h1AHTUa6M46FeEO6ViT9RzKf76A9Efk0GbUmZZN1S3LzfqUGOpOeNteFJl9LcgmRj5Qh5KMNaUGyfgmZpyElX1dRaEg9hM223BKAcydN4vWlS4N5OsinheTJm6fTj2lTMQ0rdI3QkZU3j0wxyHsOrVRQLWLaU1YoOKQ/Lxm5JLWh2GTakK1eyqJ9DH6UOihJ83OV1iX5lv3tb7S++y6fnTWL1fvtRxNuF8vjgIk4n1NXByu9iRAZpf2WWhGETVshDUiGU0NnDUJDD5xCPkBL+d5FXkvygRhaU5F1+3Qdgh2bAK61ILlhZMiCoi0QchqF19rkRN5aoSCkHsLue+8NwLtvvQVUTwIhhPJKW2+s066EACtBHpJKu+es55Jn5YW0sjFNSCP0bELXS8uXRWh6ZK072JA5Ueb1ncNGSh2Z7zxagVeWLuXKKVMA2HLMGJ4E/okz3X0TR07jcOHkfc2Mp7UgGTKtSUb/7yHfjTarZRGRhv6P/eRnPe/IpzepdE8AnrhkXSGy0W069p7Jvbo8Acby6d2SPdopkVAH5fXWEoXJrgfwkYMP5mvTpvHGq6/y2N13p460awXtfM0bPh0LgIh14nk691CetFFlVhBHtaHgukNKy5tWNs+1vIwhc1/a84jl0b9j0X+yc/Mj83eTbU2Gjh1LI87HtAo3ofaLwF7AfcBCXHCE3wm3nqHJJkQiEDdnpZnspGnO/y9NoiyE3xPp55Hak/Qzecj5PlLLkTLq+9Syh94XqQXK6+kBk88XMtWFVmaA0sTaBvVZSxSE1ANoHjoUgF9Mm5a6AV9W1Filf341ppi8jSyP+S1Ud3fmz0LIjl4JAXukyZVWLkY6MQ03jcxifg+drx1Y8uijdLS1ccycOSy++25eW7CAJuAR3Jp5Lbj9mf6O8y89jSMsb1aqN2gfiL/nNN+RzA+dTXG+Xk1E+ho63ROL13p8+VAYv25/cmKrli/UHjX5SVlClhV9PkQiMqpPElGojCQ7f5+1bh+Fya4HsH0yGbZ13bqMnCWEzGlWHWmoZOSfxwekkSZHNSbIvPmz5Mwii7z3mLYsUazONNNkyMzoETMJSrOMlktH4BHI5/O+sXgxVyTRdlvvtBMdOB/TYNxE2ruBZcD+wLeBLwP74MjKh5fXG3QnKcnIP28JHaAQ0qp8Jyvr0Ok+TWo+0iwnzVmhzfP8OWnG1deS96RlkHJIU5wmN19+kMinVw6XC/9KjdtH8g0Wh/Q5+akG/qglCg2pm7HThAl87+yzAVj+wguZ+WNqeAgxUvK+hVjdsTrydjxZ+dLmPVUyAsorT546qyEw+WxiZsSs5xeLvItdX9fXQeeORuf319AjWynDO0m0XWNzMxtxL35rknc48DzwMm7TscOBDwEP47So53Hr6dWDKS80YteahX6+shMOaUQeMopR/y++I/bmUG+Wk89eB53IwJXQZ0zTkyZ2KZOe3iARqkvK0qE+fRkdFKFl1ME4sj3H/LFdQUFICWImkK5izC67YIzh/MmTeVZsNZFHY6jFn511naxr1NosV42JLE+9tagnpuFUcr1KTa2ScGL5ZB7pt0vLL+UxwPo1awDY9zvf4bFrrqF9/XracB3A25RG2g/jJtXuDRyVfD6M8y+9gFsvz+/H1NPQJCSDGCCsFYV+x4hokMqvtQltltNh+FrT8JAdfRZxyM88JBQyp0n/pSYR/9ykuVIPgrVVRpb18uhnXysUhNTNaGp2cSgrVKh3luMeOgclVPPn5w0Q8AiNwEPpsfqqzZMnECINefNrjSMmQ5aPKa1+jxCJ6HOh68f8RqH8IZ9D6D9+9u67WfiXv9Cy//7sfsghLLrhhk2jfoPrvN7CrShugbtwK4vvDezcAPt0OB/TozhiWoULMe8OP0IIaWSkB29pvtPQOU1EOrLOm8c8Ogh3yDESkudkmLkmMZ9Ph6LLdqFNlNKcqKdC+HwyGEETsfz00KQQ0z6lbLVCQUjdDE9IUL0GlkVGXdWCuiJDFuFVqpHEECKSWmg0lZRNK59VVj+/ENHpurPMfb5enTck08b163nwooto2X9/Gpqby3YMHURJA9iA20p9BG5vpv8Dxna47S52bYZ9W522tBAXPv46pRUjuouYdEesfT8eIXOdNsNJDSJmvpNEFDP/xdpe6P/0BCrPec1FEkWjOi+j83Qknh8IyDBsAvlC2o/MG8oj78cTZyhfd1iUCkLqRjQPGcLnv/1tADYkJhONWpBFV7WgrshRrYkO0rWvarWqroSMp8laCRHl0YbS6k0jqVDZWP3a1OfNdu2trWUaoF8twJummnDBDuBWEX8TuAXYqTXZe2kY7LsWnsAR0z9w4eJrqO1kST0yl5pDmlbUoNIlEcX8SDIyTgdI+Hyx5x0aLMRIKCSXPC81QOnX8Wa00OrdIZkI5NPtKhS4E9LQpOxpFoZaoCCkbsTEQw6hZb/9WHj//fz97rtzBxp4dHU0H0NI28kyzcWQx/QYyifzxialZtVbDYlWo+FUWiYUmBArE3rumkj0+ZgmpUfFOq/X1kfttlunuTFeTt/5tSbn3krSRuIWb70JGLvWbav++ZGw7yp4EngcFy6+LCnj/UzVak1au5AmsphWpJ+tjKzzv1HfdeRayDylfU76vEzTHb88FzLD6XK6rA7vJpDHy0jgfFq7kdpZ6B5D5thQKH0tUTEhGWOGAeuttb3h1+yTuPb882lbvz73yDcPKhmpVHKtvP4Oje6QpVIzWAiaDGLaYp6oRom8BFkJ0cXyZg0e0joq+YxeTuYjfXnOHJ67+26WLliwKWxXbm3t25YnlEZKxLQ1TmO6CRizykXlHbId7PMaLO5wWtOTOD/TSpzZr9L5TDEyipnQqiEif5+hegz5SMg/m9hKCfKZyoCFmBlOQxJUljYs88i8MVOvvEaIJOW5tDp6nJCMMQ3AMcCxwEdwZubBxpiVOE3+Qmvt4hrL1S8g/Uce/mWv1IlfyR8fq6taEszTWeclyEruO8tWnecFyXvPabKHiC8kWywoIq9MWSHesTSp2UiZdN4VixfzP5Mn89XzzmObnXZi2YIFZX6BBjprZ96kB25TwNXJ921xvqY/AiNedZNr9xkDLevgpVVOW3oCZ857lVLYeJbW5GUOLQXkoc12IeLxvzsC50ITP/11YpNp/bPxJCRNe9rUJmFFXm2Gi0GbJEOk5tNjWo4sE4q6DGl4/jzqXOz/ktpbrZBHQ7obuBOYCiyy1nYAGGNGAp8C5hpjrrPW/k+NZeuXyGvm0ajEkRhrQDGzXCVElSZHGnFVO7oKlUvr0NJekDSZski3kijHGHn5uvPO/Qrly2NajZk7DaX5SHJOiuxwvWbkv0sfhZ+3NBgXKt4BbJOk/QVoXgZ7AruNgt2HwQHLYHEbPIUjqBdwyxatSerSJiHZqctOs0Gdl2k2ksfDTwbV/ht/TyFtSJOQLBcy5UkzXSjEuoGwtqTvScohw651lF3IxyOhy8jtKHTgRAixyEl9vWpWgslCHkL6V2vtRp1orV0FXANcY4zRkZMFgFHbbw9AW2trcAQag3/J0sKQJarVCuT1QuXS5KxUVQ81cLk0Tqz+0DPIY76Q+bKQ1/yWlj9WrhpiTrP7h8qkldfn2ltbN33q0XpohO1HwY3i06cPBt7BkcvmuB1rnwUWroT3rYQJTbD/rjBxHSxb5s49kxxLcT6ptZRHimlznZYrpJVA5zahJ6nqZxLSUkLEkoeEfH5pYkwjDRmGHYPWerLafIiUtYlTT4xNu3alWlStkElIkoyMMYMi5NQpbaBj15YWvj9vHljLihdfLHtR8qLaEUjWCDxvQ6pEHa+0U4+VqaR8Grl31VRXDUmHkGcAkYeYPXw+WSYtAERvn92YmJEbmps3kY0cnWvzkx65+0GVn8O0ETdZtA1nmhsKjMFpQTe1wZaLE63pvTBuGHxiJby0Ep7DaU3/xAVCvInzBYTWgJOyaPOcJprQ6gHap5KlDXmtQpJQjER8mtyeQV9Tm9Y0QsQqoZ9FSFPU14RyLUv7y2Q+TbYSsT7Aa2G9YbIDwBhzMfBZY0wb8AousOZxa+35NZapX2D7ceMA+MXpp7N44cJN6XnNNRD/syvRiPI4k6vVjLIaY5aZK8u8liW7idSR1yxWjfkzRB6ynlpMJM6bLzTyT6sjq534tc18xyTvVa6n5j/9MjrgCMl/fyM5PwZHVouAR152mwPu2QQTxsMHGmHVq/DiKkdOz+FMeq/hAijW05lMJVnmISKP0NycUGc+SHxPIyFN1ForSzMBor7H5v745x0KQU8zyYbklNqWhn6fQz4neS5t3lstUEmU3QHADtbadmPMGFyAzQe6R6y+j7bENLJ0cTzeI9bhpnUyWfV0dUQP1TW2an03leQJ5UtDXrNb1vPOI281+fNqr1lEEysbeqYdSbscOXZsWWca6rBCHajvsKT27onJUlq9YQjOz7QBt9fSe3Ha1B1tYJ6BXYE9R8JHWmCideT00gpYjAuEWILb0fYtShNvpVzSRBbqsEOaSYiENHGESEjWoZ9XbCKq9oP5ZyfrhHL/jr8fMsrI5Xuy6tfndJ/TEUgPaUn6OXUXKiGkB4H3ACustctwmvb/dYtU/QA+ws5/VmtvrTRQAMobVx5NoVoSi3XqldYTQky7yhO11hUzXuh5Z/ny9PXyzKuSCGlY8lp6fkXMxJdHjqVPPknHxo18ec4c/vnXv7Lk/vs3aZohcoqZv6wqYyjvZP1oejCuk1mZ5N8aR1CrgBtWwZBVsBuw+xjYZx+Y2AFvvgbLXnYa02LgJZxJ5i1cGPn6yP16hDSDELki8ugVIOQ5XU+I6NJ8RrJePblXP7uYZhUz+3WITxmskWdQ2UGcANLad3dpR1AZIV0I/NkY82scOT1urX2re8Tqf8ijDUG+TrzScOhKyLDakVAsaCFLjrRrxMxfWS+Evm6l2kysnrS6uqIVSdNX3jpDz6ZdpYXayfLFi7n8tNM4/rzzGDFmTKc9lwyd92GSHaP28ejNAP2nLLcu+T4ERybv4EjrfTjCehV4ehkMW+Y0p113hA99HD7cAG+vglf+CS+vd/6mF3Ga00qcBraR0m65IXnl89DEkEZCejM6n0fPHwrVrYlEnw+RZZo2p9uBNKsNiuTxkFpjWnpsOkbIRNcd0XUelRDS/wC/Ssp8F/iAMWaItXaXbpGsj8Ob7PxnqMFIm6yHHv2GtIFQg4jVLz+z8kN5FE6oQ84yb+UhGImsa2QhT0ee5vSvJHghrwkwljdNo+xOk54us/qllwBnvktbJFR2jmnwmpL3I2mSsriO0wdR+K2yvcYzGLeQ61Cc2eWJF2HQi7AjsPs2MO5D8IFhsH6dM+0tfd5tlfFScrxKybTnTYb+PvScHUkoHjESkmY0aR7M0njkM5RErSPwZL5YWY9YqHjWOxgbvMWIKOucT88TrVcNKiGkpdbaOTLBGDO4xvL0G2iTXWx2fRaqUY/zdI5pGlsW6aShK2HpWZ1yd2gnleTvyn3kKZd2vUo1qFgeAwxO2uTg5uZOG7ulyRJru7LTDnXovj7p9/GrN2yG64TewQVDDMb5m5pxkXp3rABWuIm4OzfDjnvCRw+FjwLvrIJVr8GrL5bI6TWcOXB1Uj5m3pME5clC/pb3IufwhLSfkPlOk5gmfq1RyWccIydtrtPoUJ8eoWAEPR8pZPqT53oioAEqI6SFxphTrLXzfYK1dkM3yNSvoDv4NO1Dl9OopENOy1MpyVTTWcfyhfJWShRZ5Sp9TpU+j7z15CWePHLkmb+WB9q/qMOBYx2ulC/UltJk8+d8RydX0/aTcNtxJr1BOFOczzMat4FgK/BEKzyy0JHVjsDOY2HMLrBrC9gOWPu206CWLYZlHU5zeh1HUG9S8j9J02aIgELajL9nTTSa1BooJ6eYaS6kjennpckpRDQyXVtEYpqO11RDvmYZaSfbhp8e4Feq6I6Qb6iMkLYF/tUYMwW3NcpjwEJr7dXdIFeXYYw5FJiPe34XW2vP6snra5PdJrnE92oJJi8ZpCHPtSshjiw5KjHDSaR1xDFTZJ4XRdaZx2wYMvPpvHmvGyuvfUC6XFb93lwWqtufB2hN2uTG1tZOnY/MJ79rktIdqex0/VYWUCId34nJdNmp+Ym368U1N8OZ4d5Jzg3DEZTXqB5YAquWOHPfSGDcXrDdzrBLCzAINq6B1Svg9aWw/GWnPa3AkdNbOOJbR1izkBqTXm1Bkpb/DOXTmlUsoEE+a/3fh7QS+dxi5CPz63BuSVyagIzK00HnOrUstURuQrLW/htsMtPthVvwd1+g7gjJGNMI/Bz4NG5S+MPGmButtU/1lAzaZAdxJ2UMefOl+UgqNV2l1ZuFmNbXFc0uVm8eEk7TRCr1Xemw4zxlQ/nSSCetXIzEZP5QEIMst0kzEZNjQwEgIfOyjqrTEXYh+X25BpHf+0h9FGWoczM4LehdSoS2WZL3DZy2NBinUe2ZnF8DPPgkLH/SBUsMB7YbAzvtBXt+AvYcnFzgHXjrDXhzBbyxDFauLmlRfkmjVhwxymAJTzqDxGcjnTUpSUKh9BCpS4QIJu17Wlp7IE2Xkf9ho8rn033dUm5NYLVCxat9J2a6vydHvWIfYLG19nkAY8xVwBG4ZbV6BG1qFAqVmelCI+lQ3lCDkC95V5dkr9acJlENoWSZtfS+MGnQ16qEjEJy5g3hTpMxTfvME3oeIhp/Lkvb2Si0d21GCsnqy+r0kEkHyrUgSUY+6MGnaS3Kf5f+ilYcuWzABT2sxW2x7rfHWJqUaU7q/HATdHTAqg54aRn8Y5kLlNgd2GU8NA+BEaNg7D4wdotEmFZgDbSuciHnq1fA22/ChrWwvq2kpb3lsrExueYgHDEOxZHj4ESOIYlsjQ3Q2ASmARoCkUgd7dDR5uS1HSXzpSSUDsqJRJo4035L05pPg3SCsqKc/90m/guppetBUK1QyUoNT5CszpAcTwDHW2tn11imWmAMLhDHYylOm+sxhDQkjdjoCNKJJG8j6IoW4cunkWLWCCmtk88KfohdNy8ZZJXRebNCzGOy5cnvy8RIJ5Tf151GTvraWcTs6/JtsjHRkHw5vxVFyBwnrxnzM+nRv/aj+HRPRv6Z+A7Rj9KhtCSRl7tDpPmOsw3HJdvhiKIZeLvNEdC2wJbNsKrVaVFbbwcrXoYX1rp193bZEZoHu6CI1nXOD7XtHsmxeTM0DnFssfFdx0jrKC17PlgcsmduoxSHviERboM7Nm5012nbCK3rYeMG933jBmhrhfb2zt/b26C9w1XrNz5sVJ/++Wk/jyQpnyYXzQ1pu9qk5wcpWmPy9ev1/GqBSjSkT+JWZvgAbjuKK3Fbn9QjIeWCMeZE4ERwbas3UYktNmt0kreRdFXdrrW63hVU+mKk+aZ6Gl0xq3YVWZpbQyDNwz9DWYcmXp0ee+6V/BfaZ+dNYj7Nv0uNg6C5NbluO7RvLGlaAM1DYfgIaB8Ow7bAdQKNJGrLGmjrcCrZGmAtdKxLrjeUkmokwxQlc7ayiZza2x252A6nDXW0J5+JZgSdv3ukhV6noRa+ne4I685CJT6kVcA9yYExZhzwX90iVdexDBc96rFDklYGa+2vcHOr2NyYmj7/WFBD2fXF9zzBDqG8tei40vwTIVTSeVTb6ef1ycTyVnr9Skx2lV4v615qUR90JoMQuUCpTbYLk11oboz2G8TIy4oDyk1wMo82M8fkg9JqD1LDasT5lYZTMpn5VSCGJflHNMEWHbCiA5avdXKsAMavgHEtsOeWsMVIaNwaRyodOO3nbXjrH/DmX2HVqx28s6rkR9qQXPdtShsOmuT63kTXTMmvNUgcjSQmu+TmjLhJK0hJm+ukhhMyzaWZ7dL8SASuo012HSKfTNdlEZ+1QiUmu92stc/539bafxhjPlBjeWqFh4FxxpidcER0DPCVnhTAm0UGNTd3egG7GtQQM3nJ86HvWdfWefM0tmoIp1oSrZXm15Ugi64ENOS9TtY18mg1adceJMzJ/j/WBOJNdp4MQis66Hwyv0zzeWVn482AXmvyFgrvw/CmoUE4YhiWpG2G8+UMT36PxRHEO4mMf2+DnYDth8DosfDe3WHQyCTzBmh9G1Yth5V/T4Ia2lygxAoc2bRTsry10pkUPdH4qLomnBKlQ7/Lghs63KFNmoi6PbTZLOT7IfI9z2eIgEL1aTmkLLHftUBFSwcZY3bBdfCP4wYDi4wxw6y1a7tBtqphrW0zxnwPuA3XJi6x1j7ZK7JQ+lNDEVbVdOh5fCJpeWPlslayJnA+S4ZaaHBd9YV1ta6ukEtWPaHylVy/krr1gKiBUgeQRmqaiDQhaUifRCjc24hPSWByLbYGHEl5DWA9LphhyyTPWzgSWo4zhezQDDt/ALZ5r3P/sNHNS1q+BFb8GV5dXQr7Xo2zwq2htOWFhJTLk0toXlJWHn2voTlPoevqZ+kRIqm03/qcJB1NRG0iPS3SLyRPLZFnC3NjHT6V/H4fbqXvluTz78YYrLW7d4N8VcNa+3/04uKv0iyibex5ghfStKNQQ7CB72lO7yw/lIceOcfq83mr1SCqJa5amPYqHRTkrSeWr1KNKm/9ea/fIdqmL6fJRmtKqDR9vZBW5M83qk+vFUkzoSctv7dSO05DGYlz1byLM82twmlHuwA77uW0oEHNsGGdW7nh6Yfg1WUuomkZpTXv3iK8tJAkD5mmiUU+iwZxrkH9lhNmZR4ofz5pqz7oZyqh331NRjGS0vlD5CTTO+isNfm8IeKrFXJtYW6MuQa4wVr7krX2JeAlY8xtuC0pjgcWdINsfRoykinmo6mkI8xrs5V1ap+ClMOPiNJIS9aZt/HlDfGuRNPKo1HKPJWQaKX55YBAn8sTOQdhx3+WdpqnreiAgZic0pwsO0yvjRjiC9Rqbcl3qB2UIrJCg50GXJvzHb2/njd9teHIohG3pQA4DeglHEm9F/jEKBgzDjYfCW0bYPVKeOpv8MoqF0b7Ek5jWs2mOIRNBCQ7f/9bk4lOCxGITPfPQz4TWZc248m8+pnKQYE8h8gj8+nnKxHSYtJIJMsnpMkrRna1QB5COhT4OnBl4pNZjTPXNQK3A+daaxd2g2z9CjHTRuhPDTVa3SGk1euR1mCyOriQFpfHlJSXZGOy+dFzDF1dgTutTF5TWSWaVlqdWdfrEHmyVk5vUPnSSL1dfLZF8oQ6P0/GUc0r+dR+F6kN+Tw+pHstzj/0niTPatyWE1sAOwOfHAvb7AiNjSUT3IN/Ky2uugxHQm/jyGdjyjPwJOGDLmSUniQQTbbSJBcyX0qS0gO9kEYVWklcyqv/h9CANqQ9he45hCxCSTsvSb07kGcL8/XAL4BfGGMGAaOAddba1d0kU7+AjrJL64AlshpUtRpINT6WSk1ZaajWXFatKa8Wdcc0npjGG8ofM5loNFDZdh1a+0mbHCvLQWnpoNbW1mCwQqjtxDpLfQ3vK5Gdr58nMyT57ieWbpXkX4Ujl1GUJrCO2NqFSr/5Giz8s9t24gVKu8quprSygjQvSQKQq45vunc6ay1+1N8m0rS/aKOqW5+T9y8/Y8/PP6eYmTNkzvOQWgqBPDHNOA1pg1dNUF15J7NQ0UoN1tqNuDULC2Sg2g36sqLn0kbYErqBpY1sZbocNVYaZVftiC2GvCSW97qVvEh5y3eF9NLKppXLI0eadrftbrsB0Jys9u3TpdkozUQnf/v9eGRwhPQN+c7bR62NTMq8iSOXLYHxwG67wpZbO1/QyqXw5DPwPE5behlHQm/jfEl+9QApp4cMXW5Q+ST5ys5VmxGhnMy0v0fWE1ucVSNEVpIAs6wiUhZ5zod2d6h0CW3K9Wlavtg5fT2P3jLZFagBKvnzQo2zkjplA692NJM3kq/aoIhKNcM0hJ5DXhNX6HwewstCtSTZVbJJy2eAHSZO5Ctz5kBHB8uffrpTB6e1Pd8Rdqg8vpO06rwv5zWjdTjz27Y4MnkZF6SwB7DHe50m1LoBVrwEf19c2sL8RVxEnI6EkyajWCctlyEKEZPPI4mhTdTlCUiaNkNmPP1cQtpNaMKuhyaokH8qpo1rDcujI5LX30tWUISG9i3pe+21eUgFKkPWxNiuaEyVdJih8qH8eQmiUuKo5j5jpjGPrAVK8zyfWACArD8GWV+lUYixvKH8sZc+jcRk3aGoxxFjxgDw20mTeGXRoqAMstOVHX/MrOQ7Eb/oqNdgRibpb+B8RbsA+20J277PLY3z+jJ47GV4FkdEL+I0IW+KCz0v30lLs2RM+5fl/SHJRZKOJCdJaJqcfF5NHpIcO0SZjZQTjNSgGiiZ+xBpIYKSAwB/b22BciFoEtWQ7Swr+KEt5XwtUMnE2D31atnGmAOttffUXKp+AG+q23XCBO6/4Yayc3k76bydf5r2lYd8QsiK4AnJVo2WU63ZKq1sJSRYqcmtK/LGyueV10S+p0Hfo2+X7yxfXlZXKAxZawVQGvEPEjIMSs6vx/mJtsKZ114EtsEtIjl2LDQ1u7XjHnjCkdDTOLPcclxYtp6I6mXzaXKxUCmfJCYZeCCJS0bWeWhy2qh+e3KSBOFJQZKBnG+k/Sz6GcoOX2ulUsP0ZIbKFzINajOi1phiPidN9ATStcYl8+l6a4FKNKTfG2MuB87GtbuzgYnAx2osU7/Acw8/zMYNG/j6zJksvPdeFt53X+bIXyOPkxsq7+RqEWBRC7KEdCd/VgRiSKuS18xLvJVqnLF7D43U89SdlifPf16N+VFqDf63vpZ8vnJ+je8APRFtwGlDI3Da0BpgN+Bft4SttoY1q+HZJY6EnsJpQ0txfqR1lC/6KTt/38GHCEf6ifx9SeLRxIQ6L7UZKN/DyQde+Pv2bTFETr6sHKj5c02UP0Mdfaghn4OvS5rtvOkt7X2PDSrks8yzA4AeFGjoAUqtUAkh7QvMBe7HLZh7BfDxGsvTb7B8yRLO+8EPOO3nP2f4yJEVaQ+VEE8MlWhNseuljX6ytIusjjjvvXTFtKkRkynv/KFQfp1Xk1JW3SYjT4xc5Yg47b+QHZihtPXExtbWssg0eW3Z8YZ2Q/VLwDXgtKF23ATUzXEj1J2cVZAVy+CRtxwJPYOLpHsdZ76TpjEP6Yfyz0V2/P45SW1DyiqfkzS7+U49NMlTk6wsF0vT15Wdvyc9KEXmyahDVJ1aE5XPQvudQtcLtYtQGa0FQ7yP8PK3qXR5re6KuKuEkDbiBjR++48XrLV5+oEBizWrVgHpW1B4VDuar7Serjag2PWz5NKNP8sPlLfuajSSSupPqyMvqVaiwYbMKHnKheoI5R0UmRTry+mOS3aCfuuHdpxGtB54Bbdj6yHANmPg3bfhyWWOhBbhtKHluEmu0i8UgjYB+dBtrxlJrUhqJJog5HdvhpMj+fZAPk+CmpwaVJoOfmgnTE5QIhQ/IViW07JLaG1HPh/pO5Lms6ZAOSvKxEx2Me1HE6QkdI/u6PwrIaSHgRuAj+CmDPzSGHOUtfbobpCrX6BWgQ0SaT4cyDbH5TUXVqvRVFN3CHmuF6onK2w+hCwtJo8pLiZPV8x7oVGoHv3qcr5zlSZPWX7LHXYAXJuUHXzMryCJCBwRrcVFwO0EHNjkVk54fQX8eZnbJO0pnDa0ipJJLm+7k057T0T+t4xY0/fqy7arMiGNSWtc0pRmCf83ek5TiJykduTl01F2uox+5qFrQvi/tuKzlXKZZb2aeKSZUtavNUhZrpqBXjWohJC+Ya31SwS9ChxhjDmuG2TqN/Ca0bCRI1N9JSHk1Ziy6vGQnVPeOUmxa6Wdr4S40vLmMaOl3UdauUrINu15xAhEls1jjgvl1fnkaDgkiy+j5ZL53ztxIseecw50dLDin//cNLLW+XxHNEikbYXTiF4HxgF7Nrv9hJa9BfescKstP0u5NhQKda4EIW3Jd+b+vNaWoJx8oLP2BOUddihar0Glyev485KcdGSe/z91GV+vDzWHsPaEyCehtVoP/f/5wYahXAOLtX05IJHPJRZ23l2oZD+kTuvVWWsvr604/Qvr16wB4Os/+hF/u+461qxe3SlPtX+w7ujyklOeTrwSdMVs5VGJZqFRC1t2NT4cWS5r9FjtvYfOZWlIUi5d3zY77QTAFVOmsHzRIgzl24r7DlHOb9kCF7TwJm4VhfHNYBpgyXp4rNUR0T9xROX3CqrFhEk9Wpf+JGkC02Qs7zdETP53yP+kTXgevoPPMrP6zl9PqpUmPEmI/jlLn42/vl5fr0OU1wQjSQXSTYBZZjdtptWWlTwBEV1BMQ+pG/HQrbfym1mz+Nr06bz/k5/kr0n4t+xAYlEzIVTry6iks640hLwSLS1PfaF60+qvlMyq0aigMuKrxF8WIp+0ew+Z7/TCrLqD8uW86fiN55/flO7XdtNmpBE4IlqFi5jbM6ns2VZ4FHgMF9a9Cqc5VWKWqwTSrCjnFcW0pRhJSZ9JiJj8b9RvXZ8kF9T3kHbh328Zpejvyxvy/X3IslqTk4u7enjfUGjredkGQtF1IVNgrN3pOrvyzmehIKRuxvNPualbMrAhpN10pROt1LlYDYFVS4ZpefKGn6chRBRZzyNN5hh5SBOORqXm2NC5PJqPz6d9BTpvrNzwkSMBt2SQXo/Nd3RDk+9v4BY3fX9y/pkOeIQSEb1FdpBCraBJSQY5SG3E/w/a14QopzUiRFrabxkiDuUdtYc00Xly0kEFkpykz0oGKkjtyderZZEajA8v99eKReZpc59+PhqhZ6LfC21S7CoKQuol5OnIa2Va64qfqlIyDF0rj7bT1XuthemuGtLNY57Lq1nqfFkj0RAJh0xXBhg+YgSfO/10ADasWdPJDNWEI6M1wGhgP1x493O4aKbHcUT0NqUJrD2JUISdNq+FNBz/W/p7ZLCEr1ua6/TvDtw9S7Om1poM5c8kRE6hdGnS87LIIAXt19Hk1aA+m1Q9coKtlFeSYkjzCYWd+2vGViuvBQpCStAds46hNO/jfbvuWpZeybW6agqL+TryOCwr1WKyOtFKtKW0DjpWbx5504IHQtAmsdj15Cg9TcaQec1fR5vZdJ6QyS4km07b5ZOfZJtdduGaWbN48tZbGUSpcx2GM7uB2+BsBG7h04dw5rkluECFWvmHqoUmJfnpTXHy3kO/9Xf9/HXEXYjcfCSf/r81cXj4zl1rsTpddvKSUKT5TV8j5E+S9+XXFPRy+Pya1PTE5yZRxhOTDjvXGlctUBBSN+O5Rx5h/dq1fGvWLBb97W88es89FdcRM+mljU66SwOTL0usjrRrV2LeykNSut5KX5C0cjHiyCOPLpNH7hCRhAYR2mTXEMirI8b8/KOVTz21KXrObxH+LjABGIuLkrsbR0Yv4ExztTbLdBWajDRBhHw9mphCgQ1S2wqZ93QdPtCiUZSVvi3/roSi+CTktbTJWP6nMV9RHoLSdXm5fJmQaS5k5oMSSWm/VC1Qa4IroPD60qWcn5hKtth6606hoCFom68cQckjDXmuYyNHVhkpS9q1Q0eWvFl5s+SstpxE6N5icqU9O/3i+kPn1R2FPOQ19f/fqMrqfDJ9q7FjN/1uxJnn1uM2xvtX3DYQdwC/BK4GFuL8SPVGRnK0Lzt72an65621FNnpyvzSzyOJqV2l6+t4UmulFKAgIwHlNhmSIH3dHSLNm9Z8utRItGxejjZx7Y2E26bf9qONcg3Xtx8ZyCLL+Lo34OaS+TpkYIlfTLeWKDSkHsCqV14BYItRo1JHAFmjg0rMJX50FgsBzUKWj6cj5VxIllC5tOv5vF3xf+nyebW5rOecZpaMaXVaDt/hhMw+sTL6GiGzU8g/scPEiXx5zhw6Nm5kxaOPQnJuXxwRPQbch5vU+galzrVeEXtuoSAHrXXGzHgx31FMYwr5h7SfyeeTZWXHL2XR/jAtp0wjIouODtTajSRTrTnJukI+JV1Otq9aoiCkHsDbr78OwA/OO48Fd9zBssWLox1wtajGP5RVR1fqyzrfFWLOgh5NeoTIIm3SYEimLL+WzpdGeNokmHZ97asKkbzeH8d3hlsn849+f9ppvL54Me/DbQXxOvBHnHnuVZzG1Js+okogzV+ys5R+JklaaWa8PEENIT+irMvXB07D8LviyiAIHSmnw861RUTXrUk4ZjbX4d3+elJGqZ01qLy+Prl6uc8vzYZya41aoTDZ9QCeuP9+Lpo+ncamJsZ96ENAulkrj7lNIs3/kWZaC5kF8147j2kvzUyWds95no02f6WZy7Kea+wZhWzooeul5ZMml9C1tQ8g5AwPmeJ8Pm868SNLGVU1BLfg6XuS5YI2LF/OB3Frz90P/Bq4BecrWkffISMPKa82wXXQ+VnK3zqPJCapFemghtCisPIakuy8OU+uHC7NfG0q3dfvz8n2Emrz8l1Iew/aKTfb6X2UPPx5bzqU9ylXMJfao97PqasoNKQegp+PNHL06GievKODPPnydC5p2kGs88x7DV2fHuHHIF/oELpjBBWTLfR88mq2WuPJq2HF8upRvTYrNeDs+YMpdW4bga0mTODfzz4bOjoY9vTTvATcBCzArdBdbz6iSiE1JUncMjRbT55FlAmZPaWGpIMffL1yTTppavPQmowmH1+P3uZCbxMhtSlZtzYZxjR1iVDARCg4QsofKhPT4mqBgpB6CCuXLQPg1HnzePTuu3khsFtnnkYVypdWRjp+s/J5pAUrxMpLZJnrYtFGWYiZMmJpMbnS8mo7fVp92mQSI52slzdWTjugpXbk/9dBlDQk31GuoxTOvecee4Ax/Om00/jrokXch5tb1Bc1ohg0KTVQ7oMJTZ6VnbE20cl65TJFeYkppHX5uttVHYNEGblBoCcgL0ObuEbMotFA5/uUkJo4hCPrdFvuCJSJTR6vBQqTXQ/h2QULuGDqVExDA+/bY49cnUE1pry0/CHTkkYl5r1KymvkiZbTh38xQ2m6bF7o56SXeInJQiQfdDZXxsyLsgPR9ciRsfSNWFwnNpTSqgobcfOENuA0pDG4BVB3BbZLouv+vnQp1+G2hHiX/kNGHiHznCSHmAlPRsHFBm+6E46Z8qS24/NJ+aTp1lIe+SZJSJrYoLS8k9eopElPk5M0uWWZsbVsberQJjttMq7UtZAHhYYkUOuHq/Hic88BsMuee3IP1WsKlWhIulzWCCQmU9boXpbPQpo5L0/+2LWzENIC88ofyh+qL80MqjUhbZ7UkViSJP1o2qf7AASvCW0ODMeRVAOOnAZNnMi/zJlD28aN/OHRR1mWcn/9AVIz0QMxmabNbPL5p2lLOl2m+f8stCWFDLQIRUjKgAnpZ9IDJW1SSzPpSU1MDnwk0iZ863OyfMiXWSsUhNSDWPjnP7N08WK+Pn06Tz70EA/cemvwD62EpCrplPPmrSSirhpyrMT35PPHSDKvWbMSH041+UNlpdy6c9Pko89LJ7K/7rrksxHYDEdAw3HBCwY3Mn4dtzL3G8AHd98dgB+ddhoLFi/OIXn/gHy+XlPwHbYJpGmy8HVoP5JO1yYwXz+E9z2S0XSxunTEHSJdmh6lVuQ1mhA5yflvmpy0Fufr89DXzzJN1wKFya4HsWb1ai6eNQuAESnBDVkRanqklNekFzMBZpXX8qSZ2vT18iL0guQZiWXJnxYlF7qXvPmznkPMnCf9Df5+panITzb0I/h1lLb8HgxsDbwXFyk3Iim7ErcFxELgb8ADwKIhQzj8u98F4OkBREYeUivxGoOET5N5Ym1f+0yk1qPLeLLz/6005fnlmnRbleY8X17nkcELvk5PULEoPUmoUvvzZj95X/I9a6dcJmlmlCZoHYlXCxQaUg9j9fLlAJx69tk8dOutvJH8hs4jjpi2kOVUjGkHeREbJYVkyRolVWrq0457E8jj8+nzlUTAyXok8rxgcoQau75VeUPkBOV+In99GZE1GGeSG5IcDbgOZRXOb/Rm8n01bpmftThz3RGf+xx777svc+bM4dZbb81xV/0PkpSgvINGnGtXadD5v9Qahq9Pm/70b0lMOpBBay6ha4UCFBoD5WUbkgQcM9lJmaB8kCQ1wnZVJut96yoKDamHseCuu7h01iyGb7UVex900KbRSYx80kbjeTWWSrShSlGN9uBlSqsPwiYFmS9kTpMvatr8HwJ59YhQlvP/kbblh65PIK8MavCj2mZKHeJ6nDbUlqSPALbBaUSbJ/neAl7G7cy6ELcdxELcluEvAK8ledYBe33gAwA8/vjjEWkHBrQmo0f1ss3IIAH5P8q65DwlqX3J6zSoc74eHcjQSPn8MbngqZQ1ZlbTvqvQnCdfj59/JMvKduzl0ssQycg++b7oZ1ErFBpSL+DJhx4C4Ivf+Q5/vuYaWtevj+btDjttJaikwdVaE7PqMxYMETofGsl1qE8NLUNIC9J+K60JyXQ9StWmE7+NgyenwbigBK8JdeC2g1hDSRt6KznexZGYN814kwrAYYcdxtSpU3nggQe46aabInc7cOAHL9JXJMO5Eec02fjvcnAkNS/dtuT8J9Q5qaVILUjPPfL16ECZkLaNOOcHPG0qr9T+Qv4kXd7LYlQ+7SqQn7VCoSH1Ah649VYumzOHlo9+lE8cddSmdBM4KkEef4fWXPL6lWKaRqixaplC8mlU4svSYdTQOQLIBvJ3kC6DvqfQPUhZPYxK11qSnKi6jtJilY24oISRybFFkrYWp+38E6f9PIEL134WtyfRCkqakPcFSBn3228/AC644ALWpwx2BhL885H/m39uklRkPq0ZaG1J+3hknSF0qPO+bepJsz5ARfqBsrQmq/LJurXG49+PjTjzrtecDKWVGHwe6bOSGlbWe18tCg1JIK2TqjWeWrAAgCO/+13uCWhJlfzReYirq5pW3mi2vOW1PyiUV9YtTXR5tKkQYs9A1xfS3GKaEJRHRsmRqX/p5Yi6CacBDcZpRZ6svMnOa0Ork+/vUk48aR0elLSjxYsXc8stt6TkHHjQpOTNpR2EQ8FDpjn534Z8SLr+UJsL+alCPh9JYDI4IZQ/LfhHak26Lj84k5Nytc9JDuxkhF53EFLda0jGmDOMMcuMMQuT4zPi3FRjzGJjzLPGmENE+qFJ2mJjzOm9I3k6/nLDDfxmzhxa9t2XTyZaUrWjjqwRv8yjI2by+n6q8T3lvR99bT25L0ZGlWpzMRJMK6ft9VLjkv4IbY7ZQCnayeBMcVsCW+G0okHJ+dXAUuB5nAb0LLAY5yt6DRew8K6oK+t/bmlpAWDWrFmsXr06I/fAREi78NpITFtCfJf/uc4X0qzSNCbd3qRvpgHXTgwlE68kKN0etZYUum+pNSHya3+T14h8nX4whahD+qxqhb6iIZ1nrf2JTDDG7AkcA+yFi4K90xizW3L658Cnce/6w8aYG621T/WkwHnwxP33A3Dij3/M3265hTU5O5A85rEYYo21Ug0qjx+mLZCWByHfjUdaFFwoyCEWUZemoYVGvFpr0qNfGRAhgxb84UfgG3EmuXcp+YfWUIqOi/mF8uDAAw9k1qxZrF27lj//+c8VlBx4kFqrh/TfQLmZS0Oe09qD1qx0gECaj0nKp309chAkfZKh/FDu99LXCPmaZISdz6+v5WWRWlMtUfcaUgqOAK6y1m6w1r6AG1jukxyLrbXPW2tbgauSvHWHB269ld/MmcN7d9mFQ447Lne5anw9HrEOrhL/U15NSo4U88ga8t2kRcFl3VvITBiSPZYWmmshZZRzOvxI0Ydpb4HThJqTfGtwGs8yXETcC7itwV/G7dK6inS/UB4cdNBBAJx++uksXbq0wtIDC9JHIuH/U9l+0/4LqS2h8sbqlvmtyp+mMcnIzQ46azEhDd4j9C4SyC9NenJAppcsgnLNqlboK4T0PWPM48aYS4wxWyVpY3Dvs8fSJC2WXpe494YbADcvaftdd81lfktDmolMq/O1IKDYOfkiyTQtp1TRvQNVHnmJVt5bXiLUPpmQSRM6a0Ly5QRnVhlGaeWE5qSOd3FEsxzXCF/CBSUsBV7BrarQVRLyOOaYYzZF1l166aVV1jJw4DveECnJTho6m+Y0dHCDNutpyLq1L0fXJ2WS5ussYpLkJAlERwbq9i7fC0T+0CCsK+01hrogJGPMncaYRYHjCOAC3H5iLbh9xObV8LonGmMWGGMW1NoWmhfPLljAz6dOpWnQIP5FRNxVSgweWdqRfNny+oKgcwRPXnKKdey+4cu9V9ojdTREjti1vQ0869nJF1OayeQ9+/Myj03KDsYt4TOUErGuA97GLd3zCiUiehnXeD0Jeb9QnkCFLIwaNYpzzz0XgLPOOquIrMsJORAJDZ46RLr0sUC449SWAd2x67zyvI4OlfXFyClGTKE6Y1Fy8l4lMcm0UBn/LPqlD8la+6958hljLgJuTn4uw62i4rFDkkZKur7ur4BfAQw3ptZknxt/uekmvjN7Nt+bPZtH7rmHpx58MJgvr58n1IBjdYRIKS2KTRNTDPo6Vn3G8stysZc+hJAdPpYuOxt5LR1Bp4nCv/x+CRjtCN5AaamftTh/kPcJySCNSn1DWTjhhBPYZpttmD17NjfffHN2gQKdILUl3VF7ctEmPO1rQvzWWkUsv2yfoXfLBPIRyBPy/YRWX5BmOa0t+fM6gEj6xbobdaEhpcEYs534+UXclAyAG4FjjDGDjTE74Vbcfwh4GBhnjNnJGNOMC3y4Mes61ZivanW8+Mwz/OTUUwH45Oc/X3MZPdIITWsesdFZzJ8T0iy0RqYdoNK8EJIvj69Mj+BkOTki9L/lSNHfrw6n9eQiZ737ze+GUFp1ewMuRPt1nDb0Mm7ksxynIXlNyJvkpIZVK+ybLA+0ceNGLr/88hrWPLDgO2norCnJgQSU2pn0Nfm8HiEzns6Pyi/bRswaEXrfZB69OogcCIWi6WREnzYfysCaPNaZWqAuNKQMnG2MacE9iyXAtwGstU8aY36PmzvYBpxkrW0HMMZ8D7gN94wvsdY+2QtyV4SbLryQlgMO4GtTprBq5UquPO+8iuvIE14dIyU5atd5JJHE5gzJemLII1+l0FFG+pxM185aD/3i+Xx+8ztZpjU5NuDIxmtFPsxbmk2k47g70NTUxKRJkwD44Q9/yOIBuIhqrSDn6EBn4pCaA8S1JdnWZD1aW/Iain7vtMYE6WZvnxeRV6ZJ06AkHa1RSZl8OUlMUmvqTlIy1vaapaquMMwYO6GXZdh94kQufeABsJbjPvYxnkkmz2ahK2puHlNXiMRk+GqoBcXKdDdCJg9titNOXC+bD7LQJCTX+VpPKTxbhml3p0kuhpNPPpnzzjuPefPmMWXKlB64Yv+GbDsx07EMk5adN8R3jdX16Hcrq4yUq4P0dzIGOVDUVgEbOB8iJw/9Pi2HR6y1E1Munxt1b7IbSHh2wQJnujOGYydNorGpKZeaXIlpSx8xE51vtBpeHhmMEDLt5ZWzq/cRk092HJbSRD5vivNmO7mOnDfHGUpbgb+FW0NuJW7JnteTYxU9Y5KLYfz48ZxzzjkAXH311T1wxf6PmIWgQ6TLAYcezGhNRNbp65Eak9deZP2+jBVltOYO4XYWei+0jD5NL00kV52XWpUkW22NyPveV4K+YLIbULjhgguYsN9+HHr00TzxwAP87/z5QG06udCoKqYhyWuGGn9W8EFeSDNbmiz6ulrj0Wm+Xvli+/PSbKFt7u2UyEtqRP63jgqU1+gpNDU1ccYZZ9DU1MSUKVNYkFOTLpCNNFKSmokPGpDtSk41aBTndZ26vXoSM6qMJyGdX2pVENaYQqQk65Bkogdw2pwn8+hpI7VGYbJLMMwYu2dvC5Fg15YWLn/4YbCWbxx4II8nKzpUi1hUnG7YErIjh/KXQaIaMoqZ+tLMJBpprTZ0P/Jl0xvk+YCHVsqJqJXyJVK0c7u33pzJkyczZ84c5s6dy7Rp03pJiv6PkAlPm9BifkkZKADlwQN5zHgemsz0e6E1tbyRuD6/flfku+Gv7wnSEL7Gihqa7ApCSlBPhARwxHe+w9Tzz+flxYs57qMf5e0K1iWLvRyQPsLRJjv5mWWzTiO1WJ7YyymhX7qYLTuN5OQqyjK/n3zr/UP68CaNnghQyIujjz6aK6+8kra2NiZMmFAEMnQz0vxKMhAiZOLO61vy9ckymhSkPPo9ySKm2DXTEPMzybr8dWpJSIUPqU5xwwUXcOmcObx31105fvLkzPwybDvkD/Lw4Z7ehqxXRfCOfdkgYw1ZT2htV4fs9LUtW6+NhZLb16G1E0kw8n5CNm+/jlwTpZe2Fef3eYfSqtryeAu3zI/3C+lFXnsT48eP3xTafdxxxxVk1ANIM+HpqDwd5al9SxD3AYX8Szowx8sTWwFF+6Z0+UrasPYzedn1nl76nruKwodUx7h45kzGjBvHCVOm8OaqVfzPvPgiFW3ie0MkPQStCUmEGm/MNh2r2x8xs1toprfU6Hw98oXcELlOKAJKkqXXhnxgg1wCSBOpv1a9oKmpiRkzZmzyGxWBDD2H0Nwg79uRbVP6j3T79Xm970l26lKD6lBl/HdJKoiyIY0ppm2FIlCzEPIzSR9Trd+RwmSXoN5Mdh7vGz+eKxcupLGpiSn/8R/cdtVVwXxppquYLTymHodIJ62VaPt2qN4O9emRFZwgz8XCYjUZyWAMr6FJM1woOKEnw7WrwcyZM5k2bRpXX301xx13HG1tWUONAt0B/S7FfJ2alHRgQZ5Q7iwzni6X9T7HTHla9jToAR/A64UPqfYYZozdo7eFiOCTRx7J2b//PR3t7ZzwqU+x8P77gw0rNvLRpjEPTRK6Hl2H/Eyz9WrS0fZ3fU5fR76IMVm05iXNHNIM6cknpAnJCbH1jFNPPZVzzjmHxYsXs//++7Ny5creFmlAIy2wRw+WpF8pFO0mo/ZkGR1IIMunEVOeQWYtiQngjYKQao96JiRwQQ4/Ov98Xn/1Vb756U/zwjPPlJ0P+XpC/hltRw5pJZp8NLIacJ4WpclTvyCxCDxNRv6F1EEKWhPSewz1BSICt07dRRddxNq1a/n4xz/OokWLsgsV6HakBfnoc973EpqAKvPL8zoQIi3UOxTQoLUzTWq6Dl+Ph76HNDPf6iKoYeDhhgsu4L+nTGHr7bbj13feyZixY8t8NFDekXsbsw5h1p2wDl4IEYM+iOSRpBc6L81j+lqeUPSqxIi8eo8Wf29y6++3Kd/0bh2dAxR6avJqV3H00Udz0UUX0dHRwZFHHlmQUR0hzb+q3xMd1BAbJPo2Hirj84cCF+S7EZrUGvJPoeqA8ndSvyOS/LozwKcgpD6EK+fN45czZjBq9Gh+8oc/sPmIEWUNDUomK90oQ4hF0IUIJNYIZd5QlJx+EUILOMoVDrxMcqSnX7ZWSqtqS/LxO7GGVk+ol0i5vDjssMO44oorADj22GO56667elmiAhqh9hSLftODsZB/yefX7V1qOPJ3HmLS0XtaRggTkw4K8ulZloyuoiCkPoZLZ8/m4lmz2KulhR9deCHNQ4aUmaukLTp0eEjS0IeHDPP0kKsZyA30QpAvpC7vo9z00imSgDxxyWV/1lEiHUlIPlTba4N9TRuSOOigg7juuutoaGjg+OOPLyLq6hxZpCR/+3YptSFfRmtLkhRCYeMhYiJQVpv0YtDkFqonZN2oJYqw7wS64dQzLpw5k0FDhnD85Mk0NTRwxre+xTti4myIJNICCzRCZWNIC4DwztmQ30ZqQVBuT5cmQEmU0qQn5x/JuvuKbyiGww47jOuuu46mpia+9a1vbdKSCtQ39PulCUAGH/hBVsi3JPPI8tr349u5Dv3WyxDJumNBEhpadl0m5iurBQpC6qP42dSpNDQ1cdyppzJmp5347uGH8/ry5ZvOhyLtIJ1AiOTJM7oKlfe7tsYcrqEXOA8RtVM+mqvHeUPV4Oijj+aKK66goaGBk046qdiKvI8hNOiTEXCSlKB8TpIMWNDz8KBcC5LmuRAx6fetkc7z6yohJllXKOKvlihMdn0YP508mUvnzGHPlhZ+fvPN7DR+/KZzeQISdF7t94k5QkOHPydNctKsoBcy1TJ64vFmOW+ak8EJ3hznj77oG4rhhBNO4Morr6ShoYFvfetbXHjhhb0tUoEqEPMrSbKS70wo4CFmrdHmOa8NSb9szL+kgynkux6aEqIRMglK8qsVCkLq4/jF9OmbfEp/eOghWvbbb9O5GHHo83k6dm3j1pF0WUEMEL6+XzlhA+VBCT4wwZ/zAQp+183+QEIeJ510EhdddBEAxx9/fKEZ9QPESEmn+bxa84F4UJImFq2FpfmXNDHF/Etp/YD0MxcaUoFOuHDmTH45YwZDNtuMS++4g4OOOCI6ctGrE+jQzljYdyg4IC3AQc9x8NAk5HdelSHa68UhAxT6GxEBTJs2jfnz59PW1sZRRx1V+Iz6CeSKIR46ms2neWhtSWo9GiFtSZfTEXohYoppSzETv/+dFszUFRQTYxNsZozdvbeF6CK+fOqpTEo2bZszeTKXVbENOmSPyuSnNEWEVH89eTUW3RcLUPDob610xIgRTJ8+nVNOOYV3332Xz372s9x33329LVaBGiLN6R8asOkgAkksMl1Dvid6fy9dTvuFQmnaFwXlwUj+t//+TrFSQ+3RHwgJ4JBjj2XWJZfQ0NjIJfPm8ZNp0zate5bWmEMNNZRPt5bYskS+jCQhSUahBU3l6Ct0rf6Cfffdl8suu4xdd92V66+/nh/+8IcsWbKkt8Uq0A2olpR8WW1ey7M8kCaV0PyiGDHJazSqvNqa4lFLQipMdv0Mt11xBV/9+Md5cuFCvj5pEvOvuoo9WlqAcoLImkektRnpQ4KSrdtvVSHrlfsLbaBkfltHuVkutPFdfzTLSXzzm9/k3nvvZdddd+Wcc87hmGOOKcioHyOtHaf5lfTvkClOI+RbgpIpPmTGk9eQJOStFiEzHmQHSlWLQkNK0F80JI/Nhg/nWzNnctwpp2Ct5f/7wQ+4/PzzN51Pm5ckG5psvHqynD9vRTkZhqrNcaEw7f6uDXmMGjWKU089lSlTpgAwadIk5ifb0xfo/8ias6M394POJJBXW0ozwWktR5vmYkQV0tY81hQmu9pjM2PsuN4WohvwH5Mm8YO5cwG46JxzuOLCC1mqRuSaUDzkC6Ebb4iE/MgsZo4bSGY5iYMPPphf//rXbLfddrz++uscf/zx3H777b0tVoFeQDUmvDRS8nWGJq7G3mVSyoauo+vQ5FYQUjegvxISwAf3249TfvITWvbZh42trUw76SSuufTSTuTgoYkIOqvneh6FPqQzVgcpDJQWN2rUKL7xjW8we/ZsAObOnct5551XbB8xwFErUoJ4wIMnpjykpMuHNDBdj/QvFUEN3YD+TEgAQ4YM4WtTpnDi9OkA/GLOHG665hqeWrgwunoCVE5CA10b8jjmmGP4yU9+wujRo+no6OC73/0uF198cW+LVaBOkIeUsgglpC3F9l8K1e+RZsbLIiZDQUjdgr5ISNXMA/jUEUfwg7PO4r3j3N2ePW0av5k/n/Xr12/Ko+c/+Bai/UIhAhqo2pDH+PHjOfbYY5k6dSrgtKLLL7+cZ9T+VQUK5PErZZESdDahxVYT12XSSEnXq9Mk3i0IqfYYaozdtbeFqAFCTk6dtsWIEXz+G9/g/511FhjDwkcf5Tfz5nHjVVelkpCcgCcn+Q2EeUNZGDt2LEcccQRnn302jY2NvPDCC0yfPp2rIlvOFygA2aQUCr9OCxevhbbUQVhLixFTQUjdgN4kpFh8fwixOP1q4vfHtbRw2LHH8tVTTwXg4nnz+PuDD3LjtddGFzGVwQtQENHo0aM56qijOPvssxk8eDDgVl+48MILWS1WYC9QIIY873weEx6ka0tpFhVJXjrCNouYiqCGbkCIkGr5ZCohna7W7xFqgKF8nz7ySA455hgOO/JIAH4xbx733nMPt916a9nchZhvKFZvf8aoUaM46qijOPPMM9lqq60AmDVrFn/84x958MEHe1m6An0NefoEvQ4dxJfogvJ3MmspolBdUlvyMoYIr9CQugFDjbG7BNLzdvRp+fOUT5sgF/LHxNbDSqsn7V6ampo48itf4YTvf58PJhNpfzJnDo8//jj/d+utvLNmTVSWgYTx48ez//77c8YZZzB69GjAEdFf/vKXYlfXAl1GrUx4EpVqS6G17DoIX8cC6wpCqj08IYUc8qFOvSPlXCy6LO9vmZY3XZ+TqCT4Yfjw4Rx59NGcetpp7JIEPjy2cCE//elPeeqpp3h4wYIKaus/OPzww9lll102+YjABSz86U9/KoioQE1RiQkP8mtLnlS0byk2SV4TU+w6hYbUDRhqjB1LeBkb7UNBfA9FlWURj0wjcC6WR6M7/7kRI0bwr5/+NB+eOJEfTJq0KX3q1KksWbKE6667btMaef0VO+ywAx/72MdoaWnZtMICOB/Rk08+yc0339yL0hXoz+hOUoLwNhSx68qFXkM+rIKQugFDjbHbU1qLratzbHqTTLqCkEp+5JFHsuOOOzJ37lwaGlwzvvbaa7n22mtZt24dt9xyS78hpx122IFPfvKTbLHFFpx11lkMGzYMgBUrVvDjH/+YJ554ovARFegxhFbXlgj5fdIWYIXyvkkHM+h17rQcMq9P31AQUu0x2Bg7ErfwZ9raawMZEyZM4IMf/CAnn3wyEyeW2p8np9bWVm6//XbWJP6mvoLx48fz/ve/n1GjRpWREMDMmTN55plnuPvuu4sVFgrUDdJWCdeBCB6x/ittMnyWppZYkPoXIRljjgbOAPYA9rHWLhDnpgLfwPHD9621tyXphwLzcSR/sbX2rCR9J+Aq4D3AI8Bx1trWLBmajbFNuBWo9eKBBcoxZMgQPvOZzzB06FC+973v8ZGPfGTTuYWJv6m9vZ1169bVJUFNmDCBPfbYA4Dtt9+ec845Z5PmBzB79mwWLVrE66+/zj333NNLUhYoUDn0xoB5CCVPvow6+h0h7YHjgQuBH3pCMsbsCVwJ7ANsD9wJ7JYUew74NLAUeBj4srX2KWPM74FrrbVXGWN+CTxmrb0gS4YmY2xBRJVjyJAhfOELXwDgiCOO4Oijjy47LwmqtdWNCx5++OEe23LhwAMPZOutt970e4cddigzPXrMmDGD5557riChAgUqR/8iJA9jzD2UE9JUAGvtnOT3bThNCuAMa+0hMh9wFvA6MNpa22aM+ZjMl4YGY+roSfRNNDU1ceihhzJ06FCGDh3KYYcd1omgAFpbW5k6dSqvvPJKWXpzczOtra289NJLLFu2jKVLl6Zeb8KECWy//fZsueWWneoAZ4qbMWNGsOz06dNZvHgxAMuXLy92ay1QoHrUjJCaalFJN2IM8ID4vTRJA3hZpe+LM9Ottta2BfKnoiCjrqOtra0s8ux3v/sdV1xxBUOHDt1EFLvtthszZ85k3rx5qXWtXbuW//qv/2L58uWb0iTZxDSdkEyTJ08uq2fZsmXcf//91dxigQIFuhE9RkjGmDuB0YFT06y1N/SUHBLGmBOBE3vj2gMBmqA87r33Xt73vvdtIheP5uZmRo0axdlnn82wYcM477zzMq+xdu1aTj/99E0BB5K0AB599NFNmlCBAgXqGz1GSNbaf62i2DLgveL3DkkakfQ3gBHGmKZES5L5QzL9CvgVgDGmUJJ6CFnmsdtuu43Ro0czZsyYMnLRZPPCCy+wfPnyTNNegQIF+gbq3WR3I/A7Y8y5uKCGccBDuKCQcUlE3TLgGOAr1lprjLkb+BIu0u54oFe0rwLV45lnnim2ayhQYACimkWiaw5jzBeNMUuBjwG3JMELWGufBH4PPAX8ETjJWtueaD/fA24DngZ+n+QFmAL8wBizGOdT+nXP3k2BAgUKFKgGdRVl15swxrwDPNvbcuTAKKDeZ2j2BRmhkLPWKOSsLfqKnLtbazevRUX1brLrSTxbq9DF7oQxZkG9y9kXZIRCzlqjkLO26Ety1qquujDZFShQoECBAgUhFShQoECBukBBSCX8qrcFyIm+IGdfkBEKOWuNQs7aYsDJWQQ1FChQoECBukChIRUoUKBAgbrAgCckY8yhxphnjTGLjTGn14E8S4wxTxhjFvroFWPMSGPMHcaYfySfWyXpxhjz00T2x40xe3ejXJcYY1YYYxaJtIrlMsYcn+T/hzHm+B6S8wxjzLLkmS40xnxGnJuayPmsMeYQkd5t7cIY815jzN3GmKeMMU8aY05J0uvqeabIWW/Pc4gx5iFjzGOJnDOT9J2MMQ8m1/xfY0xzkj44+b04OT82S/5ulvM3xpgXxPNsSdJ78z1qNMY8aoy5OfndM8/SWjtgD9xeSv8EdgaagceAPXtZpiXAKJV2NnB68v10YG7y/TPArbiVKz4KPNiNcn0C2BtYVK1cwEjg+eRzq+T7Vj0g5xm4VeR13j2T/3wwsFPSFhq7u10A2wF7J983x22lsme9Pc8UOevteRpgePJ9EPBg8px+DxyTpP8S+E7y/bvAL5PvxwD/myZ/D8j5G+BLgfy9+R79APgdcHPyu0ee5UDXkPYBFltrn7duE7+rgCN6WaYQjgAuS75fBnxBpP/WOjyAW8dvu+4QwFp7L7Cqi3IdAtxhrV1lrX0TuAM4tAfkjOEI4Cpr7QZr7QvAYlyb6NZ2Ya191Vr79+T7O7jVRsZQZ88zRc4Yeut5Wmut3wVyUHJY4F+APyTp+nn65/wH4CBjjEmRv7vljKFX/ndjzA7AZ4GLk9+GHnqWA52QxtB5G4tc21V0IyxwuzHmEeNWIwfY1lr7avJ9ObBt8r235a9Urt6U93uJ2eMSbwpLkafH5ExMHB/CjZbr9nkqOaHOnmdiYloIrMB10P8kvhXNJnmS82/hlhnrcTmttf55zk6e53nGmMFaTiVPd8v538BpuE1TIX1bn5o+y4FOSPWI/a21ewOHAScZYz4hT1qnD9ddaGS9ypXgAmAXoAV4FUjfjKmHYIwZDlwD/D9r7dvyXD09z4Ccdfc8rVvjsgW3wv8+wPjelSgMLacxZgIwFSfvR3BmuCm9JZ8x5nBghbX2kd64/kAnpLTtLXoF1tplyecK4Drcy/WaN8UlnyuS7L0tf6Vy9Yq81trXko6gA7iIkumg1+Q0xgzCdfJXWGuvTZLr7nmG5KzH5+lhrV0N3I1bqHmEMcYvjyavuUme5PyWuK1rekPOQxPTqLXWbgAupXef58eBzxtjluBMq/8CzKennmUtHGB99cCt5fc8zunmna179aI8w4DNxff7cbbhcyh3dp+dfP8s5U7Ph7pZvrGUBwtUJBdu9PcCzhG7VfJ9ZA/IuZ34firOtg2wF+WO1+dxDvhubRfJc/kt8N8qva6eZ4qc9fY8twZGJN+HAvcBhwNXU+6I/27y/STKHfG/T5O/B+TcTjzv/wbOqpP36EBKQQ098ixregN98cBFsjyHszlP62VZdk7+xMeAJ708OJvsXcA/gDt940sa6s8T2Z8AJnajbFfizDMbcfbgb1QjF/B1nINzMXBCD8l5eSLH47g9tmSHOi2R81ngsJ5oF8D+OHPc48DC5PhMvT3PFDnr7Xl+AHg0kWcR8GPxPj2UPJurgcFJ+pDk9+Lk/M5Z8neznH9Knuci4H8oReL12nuUXONASoTUI8+yWKmhQIECBQrUBQa6D6lAgQIFCtQJCkIqUKBAgQJ1gYKQChQoUKBAXaAgpAIFChQoUBcoCKlAgQIFCtQFCkIqUKBAgQJ1gYKQChQoUKBAXaAgpAIF6gDGmKHGmD8bYxprUFezMeZesdRLgQJ9AgUhFShQH/g6cK21tr2rFVm3xcNdwL93WaoCBXoQBSEVKNDNMMb8p9gN9AVjzN2BbMcCN4gy9xhjxiff32OSHXCNMVcbY35mjPmLMeZFY8z+xpjLjTHPGWN+Leq7PqmzQIE+g4KQChToZlhrf2ndlgMfwa2vd648n2wHvbO1dolI3hW39hu4NdCeSL6/H3jeWrs/cCHwa9zeNXsCnxV76SxKrlegQJ9BQUgFCvQc5gN/stbepNJHAav9D2PMjsAy67Z3AEdIjxtjhgAjcCtCg1v49NfWbV/QBrQDreD23QFajTGbd8+tFChQexSEVKBAD8AY8zVgR2Bm4PQ63KrJHh/ErQjt8eHk917A3wVRfZBkB9dk2+lXbPlqyYOB9bWQv0CBnkBBSAUKdDOMMR8Gfgj8hyCTTbDWvgk0JhoQuJ1YhyRlxwFH4Ex278dtTeLxAUrEVUZixpj3ACuttRtrejMFCnQjirDQAgW6H9/Dbap2tzEGYIG19psqz+24/YfuxJHLemPMYziSeQo4HmeuewggIa+hCZlBOTkBfAq4pTtupkCB7kKxH1KBAnUAY8zewKnW2uOMMf8A9rbWvtOF+q7F7T77XGbmAgXqBIXJrkCBOoC19u84DWpL97NLZNQMXF+QUYG+hkJDKlCgQIECdYFCQypQoECBAnWBgpAKFChQoEBdoCCkAgUKFChQFygIqUCBAgUK1AUKQipQoECBAnWBgpAKFChQoEBdoCCkAgUKFChQFygIqUCBAgUK1AX+f480NzEAfe7hAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# simple drawing\n",
    "u1.draw(logarithm=True, scale=\"scaled\", draw_borders=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY8AAADiCAYAAAC2h0vXAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/H5lhTAAAACXBIWXMAAAsTAAALEwEAmpwYAAAUbUlEQVR4nO3dfbBd9V3v8ff3JBCwpUmAY0xCayhFFFEoHHvp9GGwlBZqbagiAz40VcZM74jTWp0aLzN3vDPemVZHEEfHToTaqLUFsQyxjrZpJD6MFntCAwTSNoELmgeSQwqB2k4g7O/9Y/1OujnZ5+F3cs7aO+X9mtmz1/qttfb65rfPOZ+s58hMJEmqMdTvAiRJJx7DQ5JUzfCQJFUzPCRJ1QwPSVI1w0OSVG1hvwuYrTPPPDNXrVrV7zIkqe+2bt36VGYOt7nOEzY8Vq1axejoaL/LkKS+i4gn2l6nu60kSdUMD0lSNcNDklTN8JAkVTM8JEnVDA9JUjXDQ5JUzfCQJFVrNTwi4ryI2Nb1ejYiPhQRp0fEpojYWd6XtlmXJKlOq+GRmV/LzIsy8yLgEuBbwN3AOmBzZp4LbC7jkqQB1c/dVpcDj2bmE8BqYENp3wBc3a+iJEnT62d4XAd8ugwvy8x9ZfhJYFl/SpIkzURfwiMiTgbeA/z1xGmZmUBOstzaiBiNiNGxsbF5rlKSNJl+bXlcBdyfmfvL+P6IWA5Q3g/0Wigz12fmSGaODA+3evdhSVKXfoXH9XxnlxXARmBNGV4D3NN6RZKkGWs9PCLiFcAVwGe7mj8KXBERO4G3l3FJ0oBq/WFQmfnfwBkT2g7SnH0lSToBeIW5JKma4SFJqmZ4SJKqGR6SpGqGhySpmuEhSapmeEiSqhkekqRqhockqZrhIUmqZnhIkqoZHpKkaoaHJKma4SFJqmZ4SJKq9eNhUEsi4q6I+GpE7IiIN0bE6RGxKSJ2lvelbdclSZq5fmx53Ar8Q2b+IHAhsANYB2zOzHOBzWVckjSgWg2PiFgMvBW4HSAzn8/MZ4DVwIYy2wbg6jbrkiTVaXvL42xgDPiziPhKRNxWnmm+LDP3lXmeBJa1XJckqULb4bEQuBj4k8x8PfDfTNhFlZkJZK+FI2JtRIxGxOjY2Ni8FytJ6q3t8NgN7M7M+8r4XTRhsj8ilgOU9wO9Fs7M9Zk5kpkjw8PDrRQsSTpWq+GRmU8C/xUR55Wmy4FHgI3AmtK2BrinzbokSXUW9mGdvwp8KiJOBh4DfpEmxO6MiBuAJ4Br+1CXJGmGWg+PzNwGjPSYdHnLpUiSZskrzCVJ1QwPSVI1w0OSVM3wkCRVMzwkSdUMD0lSNcNDklTN8JAkVTM8JEnVDA9JUjXDQ5JUzfCQJFUzPCRJ1QwPSVI1w0OSVK3153lExOPAc8CLwJHMHImI04E7gFXA48C1mfl027VJkmamX1seP56ZF2Xm+EOh1gGbM/NcYHMZlyQNqEHZbbUa2FCGNwBX968USdJ0+hEeCXwhIrZGxNrStiwz95XhJ4FlfahLkjRDrR/zAN6cmXsi4nuBTRHx1e6JmZkRkb0WLGGzFuA1r3nN/FcqSeqp9S2PzNxT3g8AdwNvAPZHxHKA8n5gkmXXZ+ZIZo4MDw+3VbIkaYJWwyMiXhERp40PA+8AtgMbgTVltjXAPW3WJUmq0/Zuq2XA3RExvu6/ysx/iIgvA3dGxA3AE8C1LdclSarQanhk5mPAhT3aDwKXt1mLJGn2BuVUXUnSCcTwkCRVMzwkSdUMD0lSNcNDklTN8JAkVTM8JEnVDA9JUjXDQ5JUzfCQJFUzPCRJ1QwPSVI1w0OSVM3wkCRVMzwkSdX6Eh4RsSAivhIRnyvjZ0fEfRGxKyLuiIiT+1GXJGlm+rXl8UFgR9f4x4BbMvN1wNPADX2pSpI0I62HR0ScBfwEcFsZD+BtwF1llg3A1W3XJUmauX5sefwB8BGgU8bPAJ7JzCNlfDewsteCEbE2IkYjYnRsbGzeC5Uk9dZqeETEu4EDmbl1Nstn5vrMHMnMkeHh4TmuTpI0UwtbXt+bgPdExLuAU4BXAbcCSyJiYdn6OAvY03JdkqQKrW55ZOZvZeZZmbkKuA74x8z8OeBe4Joy2xrgnjbrkiTVGZTrPH4T+HBE7KI5BnJ7n+uRJE2h7d1WR2XmFmBLGX4MeEO/apEk1RmULQ9J0gnE8JAkVTM8JEnVDA9JUjXDQ5JUzfCQJFUzPCRJ1QwPSVI1w0OSVM3wkCRVMzwkSdUMD0lSNcNDklTN8JAkVTM8JEnV2n6G+SkR8R8R8UBEPBwR/6e0nx0R90XEroi4IyJObrMuSVKdtrc8DgNvy8wLgYuAKyPiUuBjwC2Z+TrgaeCGluuSJFVo+xnmmZnfLKMnlVcCbwPuKu0bgKvbrEuSVKf1Yx4RsSAitgEHgE3Ao8AzmXmkzLIbWDnJsmsjYjQiRsfGxlqpV5J0rNbDIzNfzMyLgLNonlv+gxXLrs/MkcwcGR4enq8SJUnT6NvZVpn5DHAv8EZgSUQsLJPOAvb0qy5J0vTaPttqOCKWlOFTgSuAHTQhck2ZbQ1wT5t1SZLqLJx+ljm1HNgQEQtoguvOzPxcRDwCfCYifgf4CnB7y3VJkiq0Gh6Z+SDw+h7tj9Ec/5AknQC8wlySVM3wkCRVMzwkSdUMD0lSNcNDklTN8JAkVTM8JEnVDA9JUjXDQ5JUzfCQJFUzPCRJ1QwPSVI1w0OSVM3wkCRVMzwkSdXafpLgqyPi3oh4JCIejogPlvbTI2JTROws70vbrEuSVKftLY8jwK9n5vnApcCvRMT5wDpgc2aeC2wu45KkAdVqeGTmvsy8vww/R/P88pXAamBDmW0DcHWbdUmS6vTtmEdErKJ5JO19wLLM3FcmPQksm2SZtRExGhGjY2Nj7RQqSTpGX8IjIl4J/A3wocx8tntaZiaQvZbLzPWZOZKZI8PDwy1UKknqpfXwiIiTaILjU5n52dK8PyKWl+nLgQNt1yVJmrm2z7YK4HZgR2be3DVpI7CmDK8B7mmzLklSnYUtr+9NwC8AD0XEttL2v4CPAndGxA3AE8C1LdclSarQanhk5r8CMcnky9usRZI0e15hLkmqZnhIkqoZHpKkaoaHJKma4SFJqmZ4SJKqGR6SpGqGhySpmuEhSapmeEiSqhkekqRqhockqZrhIUmqZnhIkqoZHpKkam0/SfATEXEgIrZ3tZ0eEZsiYmd5X9pmTZKkem1veXwSuHJC2zpgc2aeC2wu45KkAdZqeGTmPwPfmNC8GthQhjcAV7dZkySp3iAc81iWmfvK8JPAsslmjIi1ETEaEaNjY2PtVCdJOsYghMdRmZlATjF9fWaOZObI8PBwi5VJkroNQnjsj4jlAOX9QJ/rkSRNYxDCYyOwpgyvAe7pYy2SpBlo+1TdTwP/DpwXEbsj4gbgo8AVEbETeHsZlyQNsIVtriwzr59k0uVt1iFJOj6DsNtKknSCMTwkSdUMD0lSNcNDklTN8JAkVTM8JEnVDA9JUjXDQ5JUzfCQJFUzPCRJ1QwPSVI1w0OSVM3wkCRVMzwkSdVavSW7NNcys+drqmlTzTtRRMx6+HiXn83nDg0NERHHzC/NtYEJj4i4ErgVWADclplTPhTqyJEjPPXUU3Q6nZ6vF198cdJp3w3TB6GGidNh5n+wj+cP/GR/6PVSEXE0TLpDZWLbbOZp6/OGhoZYsGDB0eGJ41NNa3u8n+vuy8/XIPwSRsQC4OvAFcBu4MvA9Zn5yBTL9L9wDYzxPzzdr8nap5t3XPfvRu3w8S4/2881WF+2tmbmSJsrHJQtjzcAuzLzMYCI+AywGpg0PBYsWMCSJUteksQTXxOT+kSbPgg1zGR69/8Wy/dX/Qf7eObXsTKTTqfT832qabOZd67mmTjvdFu6/RwftFoOHTrU+s/YoITHSuC/usZ3A/9j4kwRsRZYW0a/efDgwa9N8nlnAk/NaYVzx9rqDWpdMLi1DWpdMLi1DWpdMH1t399WIeMGJTxmJDPXA+unmy8iRtvehJspa6s3qHXB4NY2qHXB4NY2qHXBYNY2KKfq7gFe3TV+VmmTJA2gQQmPLwPnRsTZEXEycB2wsc81SZImMRC7rTLzSETcCHye5lTdT2Tmw8fxkdPu2uoja6s3qHXB4NY2qHXB4NY2qHXBANY2EKfqSpJOLIOy20qSdAIxPCRJ1U7I8IiIn4mIhyOiExEjE6b9VkTsioivRcQ7u9qvLG27ImJdV/vZEXFfab+jHLCfixrviIht5fV4RGwr7asi4ttd0z7etcwlEfFQqeUPY56ugIuI346IPV01vKtrWlX/zUNtvxcRX42IByPi7ohYUtr73m8T6mylP6ZY/6sj4t6IeKT8LnywtFd/t/NQ2+Pl+9gWEaOl7fSI2BQRO8v70tIe5TvbVb7zi+exrvO6+mVbRDwbER/qV59FxCci4kBEbO9qq+6niFhT5t8ZEWvmssYp1dyLaFBewA8B5wFbgJGu9vOBB4BFwNnAozQH4BeU4dcCJ5d5zi/L3AlcV4Y/DvzPeaj394H/XYZXAdsnme8/gEuBAP4euGqe+u+3gd/o0V7df/NQ2zuAhWX4Y8DHBqXfutbXWn9MUcNy4OIyfBrN7X3Or/1u56m2x4EzJ7T9LrCuDK/r+l7fVb6zKN/hfS313wLgSZqL6/rSZ8BbgYu7f65r+wk4HXisvC8tw0vb6MMTcssjM3dkZq+ry1cDn8nMw5n5/4BdNLc+OXr7k8x8HvgMsLr8D/VtwF1l+Q3A1XNZa1nHtcCnp5lvOfCqzPxSNj8Vfz7XtcxAVf/NRwGZ+YXMPFJGv0Rzzc+k+tRvrfXHZDJzX2beX4afA3bQ3KlhMpN9t21ZTfP7BS/9PVsN/Hk2vgQsKd/pfLsceDQzn5hinnnts8z8Z+AbPdZZ00/vBDZl5jcy82lgE3DlXNU4lRMyPKbQ6zYnK6doPwN4puuP1Xj7XHoLsD8zd3a1nR0RX4mIf4qIt3TVvrtHjfPlxrL5+4nxTWPq+2++/RLN/7bGDUK/ja+zH/3RU0SsAl4P3Feaar7b+ZDAFyJiazS3FAJYlpn7yvCTwLI+1NXtOl76H7p+99m42n7q28/iwIZHRHwxIrb3eLX6P7ypzLDG63npD+k+4DWZ+Xrgw8BfRcSrWq7tT4BzgItKPb8/1+s/jtrG57kJOAJ8qjS10m8nmoh4JfA3wIcy81n6/N0Wb87Mi4GrgF+JiLd2TyxbiH27RiCa45rvAf66NA1Cnx2j3/00nYG4SLCXzHz7LBab6jYnvdoP0mz+LSxbH1W3RZmuxohYCPwUcEnXMoeBw2V4a0Q8CvxAWW/3LprjukXLTPsvIv4U+FwZre2/eaktIt4PvBu4vPwCtdZvMzQQt9OJiJNoguNTmflZgMzc3zV9pt/tnMrMPeX9QETcTbOrZ39ELM/MfWV3y4G26+pyFXD/eF8NQp91qe2nPcBlE9q3zHONwABveczSRuC6iFgUEWcD59IcTO15+5Pyh+le4Jqy/Brgnjms5+3AVzPz6G6ViBiO5vklRMRrS42PlU3VZyPi0nKc5H1zXMtRE/YpvxcYP9ujqv/mqbYrgY8A78nMb3W1973fuvT9djrl33o7sCMzb+5qr/1u57quV0TEaePDNCdAbC/rHz8TqPv3bCPwvnI20aXAoa7dNvPlJXsD+t1nE9T20+eBd0TE0rK77R2lbf61cVR+rl80X/Bumv+J7gc+3zXtJpqzIr5G11k3NGcrfL1Mu6mr/bU0PxC7aDZjF81hnZ8EPjCh7aeBh4FtwP3AT3ZNG6H5wX0U+CPKHQDmof/+AngIeJDmh3L5bPtvHmrbRbMPd1t5fXxQ+m1Cna30xxTrfzPNLo0Hu/rqXbP5bue4rtfSnKH0QPm+birtZwCbgZ3AF4HTS3sAf1zqeoiusyfnqb5X0OxxWNzV1pc+owmwfcALNH/PbphNP9EcG9xVXr/Y1s+gtyeRJFX7btttJUlqgeEhSapmeEiSqg3sqbq9bNmy5ZyhoaEPAz/f6XROGxoaeg74y06nc/Nll132aL/rk6SXizk5YB4RdwEfyczH5uCzvgj8TDaX2h+1ZcuWq4aGhu5auXLlSStWrDhp0aJFHD58mL17976wZ8+eFzqdzjWXXXbZ30/ysZKkSZRT5G+luefXbZn50WmXOd7wiIgfBn4nM997XB/0nc9bA5yVmf93vK1scTx44YUXfs/ixYuPWebQoUM88MAD3+p0Oj/qFogkzVy5furrwBU0pwx/Gbg+Mx+Zarkpd1tFxAeAD5TRxcDjmfnjE2b7ObouyoqIb2bmK8vwNcC7M/P9EfFJ4Ns09+D5Xppzk98HvJHmDpHvLx+xEfgX4Gh4DA0NfXjlypUn9QoOgMWLF7NixYrvOXjw4K7t27f3nEeSvltdcMEFx/MYgqM3+gSIiPEbfU4ZHlMeMM/Mj2fmRcCP0STSzT1mexOwdYZFLqUJi1+jCYlbgB8GfiQiLirrfBpYFBFndC338ytWrDhpqg9euXIlzz///AzLkCQVs7q54kwPmN8K/GNm/m2PacuBsRl+zt9mZkbEQzR3mn0IICIepnlew7Yy3wFgBc2VoHQ6ndMWLVo05QcvWrSIF198cYZlSNLLRzmW/H09Jt0028+cNjzKTeq+H7hxklm+DZzSNd59EOWUCfMeLu+druHx8e5aTimfC8DQ0NBzhw8fftWpp546aZ2HDx9mwYIFk06XpJernOJmpBHxRmZxA8jpjnlcAvwG8JbM7Ewy2w7gdTRPD4PmrpA/RHMvmPcCz01XxIR1Bk1CPt7V/Jd79+795XPOOWfSXVd79+59ITPXX3DBBZOFnCTpWEdv9EkTGtcBPzvdQtNdJHgjzeMN743m2b639Zjn73jpLYHX0dzS+N9obvpV6xLgS/mdBzTR6XRu3rNnzwuHDh3qucChQ4cop+veMov1SdLLVvlbeyPN3Xh3AHdm5sPTLTcXp+qeSnNb8zdl5nEfdIiIW2lul765u93rPCRpcMzVRYLvpHmuwH/OwWf9cmb+aa9p5XqPXwN+odPpvHJoaOibwF90Op1bvL5DktrjLdklSdW8MaIkqZrhIUmqZnhIkqoZHpKkaoaHJKma4SFJqvb/AThgZRXMKbWvAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "u1.draw_profiles_interactive()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Directly with interactive"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib widget"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2020-05-02T18:16:40.107619Z",
     "start_time": "2020-05-02T18:16:39.965223Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "75.72461367055705\n"
     ]
    }
   ],
   "source": [
    "Imax = u1.intensity().max()\n",
    "print(Imax)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2020-05-02T18:16:40.114014Z",
     "start_time": "2020-05-02T18:16:40.109768Z"
    }
   },
   "outputs": [],
   "source": [
    "def __interact2__(i):\n",
    "    plt.plot(u1.x, u1.intensity().transpose()[i], \"k\")\n",
    "    plt.title(\"z = {:2.3f} mm\".format(u1.z[i] / 1000))\n",
    "    plt.ylim(0, Imax)\n",
    "    plt.xlim(-100 * um, 100 * um)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2020-05-02T18:17:12.846865Z",
     "start_time": "2020-05-02T18:17:12.626323Z"
    }
   },
   "outputs": [
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "e5e95dc3328d48228cb0585fbcb3f7e9",
       "version_major": 2,
       "version_minor": 0
      },
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAeAAAAFoCAYAAACPNyggAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/H5lhTAAAACXBIWXMAAAuJAAALiQE3ycutAAAS+UlEQVR4nO3dUWjd9fn48UeTCxmdoxXEtDltSUJjT2rXQhCkZW4FWbeCdaGVylbLyHGDBa0XyoTeDcGL9kYl200qG2SMmVW2Ijgo1dhJQVJpKTS0MzpNj0qgmpqLykzT7//ih+Ef0s6T1uY5ia8XnIsv5yE8fGjOm+854fSWoiiKAADm1a3ZCwDAt5EAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJKjLAO/duzeam5ujsbHxmjODg4PR0dERbW1tUalUYmpqah43BIAbU5cB3rlzZ5w4ceKaz1+5ciUqlUoMDAzEyMhITExMRH9//zxuCAA3pi4DvHnz5rjrrruu+fzQ0FAsX748yuVyRER0d3fHoUOH5ms9ALhhdRngr1OtVqNUKk1fr1y5Ms6fP5+4EQDMzbU/ZK1jRVHUPNvb2xu9vb3T1++//360tLTcjLUAmEfVajUmJiay17huCzLApVJpxh3v6OhoNDc3X3W2p6cnenp6pq/L5XIMDw/f9B0BuLm++hhyoVqQb0F3dnZGtVqdDunBgwejq6sreSsAqF1dBvjXv/51NDc3x9TUVDQ3N0dPT0+cOHEifvrTn0ZERENDQ/T19cWOHTuitbU1lixZErt3707eGgBqd0sxlw9UFwFvQQMsDgv99bwu74ABYLETYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJ6jbAg4OD0dHREW1tbVGpVGJqamrWzIEDB6KjoyPWr18fW7dujbGxsYRNAWDu6jLAV65ciUqlEgMDAzEyMhITExPR398/Y+bdd9+NP/zhD3HixIk4ffp0bNiwIQ4cOJC0MQDMTV0GeGhoKJYvXx7lcjkiIrq7u+PQoUMzZoqiiMnJyfjiiy+iKIqYmJiIpqamjHUBYM4asxe4mmq1GqVSafp65cqVcf78+Rkza9asiccffzxKpVJ897vfjfb29njxxRfne1UAuC51eQdcFMXXznz66afxyiuvxMjISHz00UfR0tIS+/fvnzXX29sb5XJ5+jE+Pn4zVgaAOanLAJdKpRl3vKOjo9Hc3Dxj5vXXX4+WlpZoamqKhoaG2LlzZxw/fnzWz+rp6Ynh4eHpx9KlS2/6/gDwdeoywJ2dnVGtVmN4eDgiIg4ePBhdXV0zZlatWhVvv/12TExMRETEkSNHpj8zBoB6V5cBbmhoiL6+vtixY0e0trbGkiVLYvfu3XH48OGoVCoREXHvvffGnj174t5774177rknPvzww3jmmWeSNweA2txS1PKB6yJSLpen76wBWLgW+ut5Xd4BA8BiJ8AAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJ6jbAg4OD0dHREW1tbVGpVGJqamrWzNjYWGzfvj3uvvvuaG9vj1dffTVhUwCYu7oM8JUrV6JSqcTAwECMjIzExMRE9Pf3z5rbs2dP/PznP4+zZ8/GmTNn4r777kvYFgDmri4DPDQ0FMuXL49yuRwREd3d3XHo0KEZM+fOnYuxsbF4+OGHIyKisbEx7rjjjnnfFQCuR10GuFqtRqlUmr5euXJlnD9/fsbM2bNn484774xHHnkkNm7cGI8++mh89tln870qAFyXugxwURRfO3P58uU4duxY7Nu3L06ePBmrV6+Op59+etZcb29vlMvl6cf4+PjNWBkA5qQuA1wqlWbc8Y6OjkZzc/OsmXK5HOvWrYuIiF27dsU777wz62f19PTE8PDw9GPp0qU3d3kAqEFdBrizszOq1WoMDw9HRMTBgwejq6tr1syXX345HeojR45ER0fHvO8KANejMXuBq2loaIi+vr7YsWNH/Pe//437778/du/eHYcPH47Dhw9HX19f3HrrrfH73/8+tm/fHpcvX44VK1bESy+9lL06ANTklqKWD1wXkXK5PH1nDcDCtdBfz+vyLWgAWOwEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0CCug3w4OBgdHR0RFtbW1QqlZiamrrm7LZt26KtrW0etwOAG1OXAb5y5UpUKpUYGBiIkZGRmJiYiP7+/qvO/vnPf45ly5bN84YAcGPqMsBDQ0OxfPnyKJfLERHR3d0dhw4dmjV34cKF6O3tjX379s33igBwQ+oywNVqNUql0vT1ypUr4/z587PmnnzyyXj22Wfjtttum8/1AOCG1WWAi6L42pnXXnstGhoaYsuWLf9zrre3N8rl8vRjfHz8m1oTAK5bY/YCV1MqlWbc8Y6OjkZzc/OMmWPHjsXRo0dj9erVcfny5RgbG4v169fH6dOnZ8z19PRET0/P9PVXb2sDQKa6vAPu7OyMarUaw8PDERFx8ODB6OrqmjHz3HPPRbVajQ8++CDeeuutWLVq1az4AkC9qssANzQ0RF9fX+zYsSNaW1tjyZIlsXv37jh8+HBUKpXs9QDght1S1PKB6yJSLpen76wBWLgW+ut5Xd4BA8BiJ8AAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJ6jbAg4OD0dHREW1tbVGpVGJqamrG86dOnYpNmzZFR0dHrFu3Ll544YWkTQFg7uoywFeuXIlKpRIDAwMxMjISExMT0d/fP2PmO9/5Trz00ktx5syZOH78eLz44otx6tSpnIUBYI7qMsBDQ0OxfPnyKJfLERHR3d0dhw4dmjGzZs2aaG9vj4iI22+/PdauXRvnz5+f910B4HrUZYCr1WqUSqXp65UrV/7PuL733ntx4sSJ2LRp03ysBwA3rDF7gaspiqLm2YsXL8ZDDz0Uzz//fCxbtmzW8729vdHb2zt9PT4+/o3sCAA3oi7vgEul0ow73tHR0Whubp41d+nSpdi2bVs89thjsXPnzqv+rJ6enhgeHp5+LF269KbtDQC1qssAd3Z2RrVajeHh4YiIOHjwYHR1dc2YmZycjK6urnjggQfiiSeeyFgTAK5bXQa4oaEh+vr6YseOHdHa2hpLliyJ3bt3x+HDh6NSqURExMsvvxxHjhyJv//977Fhw4bYsGFD/O1vf0veHABqc0sxlw9cF4FyuTx9Zw3AwrXQX8/r8g4YABY7AQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEhQtwEeHByMjo6OaGtri0qlElNTU7NmXn755VizZk20trbGvn37ErYEgOtTlwG+cuVKVCqVGBgYiJGRkZiYmIj+/v4ZM59//nk89dRT8eabb8a5c+fijTfeiDfffDNpYwCYm7oM8NDQUCxfvjzK5XJERHR3d8ehQ4dmzPzzn/+MH/7wh9HU1BSNjY2xZ8+eWTMAUK/qMsDVajVKpdL09cqVK+P8+fNzngGAetWYvcDVFEXxjcxERPT29kZvb+/09bvvvjt9Z821jY+Px9KlS7PXWBCcVW2cU+2cVW3ef//97BVuSF0GuFQqzbibHR0djebm5lkzp0+f/p8zERE9PT3R09MzfV0ul2N4ePgmbL24OKfaOavaOKfaOavaLPSbqbp8C7qzszOq1er0P8CDBw9GV1fXjJmtW7fGG2+8EZ988klcvnw5/vSnP82aAYB6VZcBbmhoiL6+vtixY0e0trbGkiVLYvfu3XH48OGoVCoREfG9730v9u/fHz/4wQ+ivb097r///vjRj36UvDkA1KYu34KOiNiyZcust2AefPDBePDBB6evd+3aFbt27ZrTz/3/347m2pxT7ZxVbZxT7ZxVbRb6Od1S1PrXTADAN6Yu34IGgMVOgAEgwaIMsO+Rrt3XndWpU6di06ZN0dHREevWrYsXXnghadNctfyb+sq2bduira1tHrerL7Wc1djYWGzfvj3uvvvuaG9vj1dffTVh01y1nNOBAweio6Mj1q9fH1u3bo2xsbGETXPt3bs3mpubo7Hx2n+yNJffz7pSLDJTU1NFa2trcebMmaIoimLnzp3FH//4xxkzFy9eLEqlUvHxxx8Xk5OTxX333VcMDg5mrJuqlrM6d+5ccfbs2aIoiuLzzz8v2traipMnT873qqlqOaev9Pf3F7/4xS+K1tbW+VyxbtR6Vj/+8Y+Lv/71r0VRFMXk5GRx4cKFed0zWy3n9O9//7toaWkpLl26VBRFUfz2t78tnnrqqXnfNdu//vWv4pNPPikaGhqu+vxcfj/rzaK7A/Y90rWr5azWrFkT7e3tERFx++23x9q1a791X/lZyzlFRFy4cCF6e3u/1e+o1HJW586di7GxsXj44YcjIqKxsTHuuOOOed81Uy3nVBRFTE5OxhdffBFFUcTExEQ0NTVlrJtq8+bNcdddd13z+Vp/P+vRoguw75Gu3VzP4b333osTJ07Epk2b5mO9ulHrOT355JPx7LPPxm233Taf69WVWs7q7Nmzceedd8YjjzwSGzdujEcffTQ+++yz+V41VS3ntGbNmnj88cejVCpFU1NTnDlzJvbu3Tvfq9a9hfx6vugCXHyD3yO92M3lHC5evBgPPfRQPP/887Fs2bKbuFX9qeWcXnvttWhoaIgtW7bMw0b1q5azunz5chw7diz27dsXJ0+ejNWrV8fTTz89D9vVj1rO6dNPP41XXnklRkZG4qOPPoqWlpbYv3//PGy3sCzk1/NFF+Bav0f662a+DWo9h0uXLsW2bdvisccei507d87ninWhlnM6duxYHD16NFavXh2bN2+ODz/8MNavXz/fq6ar9fevXC7HunXrIuL/vlDnnXfemdc9s9VyTq+//nq0tLREU1NTNDQ0xM6dO+P48ePzvWrdW8iv54suwL5Huna1nNXk5GR0dXXFAw88EE888UTGmulqOafnnnsuqtVqfPDBB/HWW2/FqlWrZvxnId8WtZxVZ2dnfPnll9MvmkeOHImOjo553zVTLee0atWqePvtt2NiYiIi/u+cFvp/PnAz1HKWdSvtz79uoqNHjxZr164tWlpail/+8pfF5ORk8Y9//KPo7u6envnLX/5StLW1FS0tLcUzzzyTuG2urzur/v7+4tZbby2+//3vTz8GBgaSt55/tfyb+sp//vOfb+1fQRdFbWd17NixYuPGjcU999xTbN26tfj4448TN85Ryzn97ne/K9rb24t169YVP/vZz4rx8fG8hZP86le/KlasWFFERLFixYriN7/5TTE0NFT85Cc/mZ652lkuBL6KEgASLLq3oAFgIRBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAn+H9SFNNSVrCuaAAAAAElFTkSuQmCC",
      "text/html": [
       "\n",
       "            <div style=\"display: inline-block;\">\n",
       "                <div class=\"jupyter-widgets widget-label\" style=\"text-align: center;\">\n",
       "                    Figure\n",
       "                </div>\n",
       "                <img src='data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAeAAAAFoCAYAAACPNyggAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/H5lhTAAAACXBIWXMAAAuJAAALiQE3ycutAAAS+UlEQVR4nO3dUWjd9fn48UeTCxmdoxXEtDltSUJjT2rXQhCkZW4FWbeCdaGVylbLyHGDBa0XyoTeDcGL9kYl200qG2SMmVW2Ijgo1dhJQVJpKTS0MzpNj0qgmpqLykzT7//ih+Ef0s6T1uY5ia8XnIsv5yE8fGjOm+854fSWoiiKAADm1a3ZCwDAt5EAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJKjLAO/duzeam5ujsbHxmjODg4PR0dERbW1tUalUYmpqah43BIAbU5cB3rlzZ5w4ceKaz1+5ciUqlUoMDAzEyMhITExMRH9//zxuCAA3pi4DvHnz5rjrrruu+fzQ0FAsX748yuVyRER0d3fHoUOH5ms9ALhhdRngr1OtVqNUKk1fr1y5Ms6fP5+4EQDMzbU/ZK1jRVHUPNvb2xu9vb3T1++//360tLTcjLUAmEfVajUmJiay17huCzLApVJpxh3v6OhoNDc3X3W2p6cnenp6pq/L5XIMDw/f9B0BuLm++hhyoVqQb0F3dnZGtVqdDunBgwejq6sreSsAqF1dBvjXv/51NDc3x9TUVDQ3N0dPT0+cOHEifvrTn0ZERENDQ/T19cWOHTuitbU1lixZErt3707eGgBqd0sxlw9UFwFvQQMsDgv99bwu74ABYLETYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJ6jbAg4OD0dHREW1tbVGpVGJqamrWzIEDB6KjoyPWr18fW7dujbGxsYRNAWDu6jLAV65ciUqlEgMDAzEyMhITExPR398/Y+bdd9+NP/zhD3HixIk4ffp0bNiwIQ4cOJC0MQDMTV0GeGhoKJYvXx7lcjkiIrq7u+PQoUMzZoqiiMnJyfjiiy+iKIqYmJiIpqamjHUBYM4asxe4mmq1GqVSafp65cqVcf78+Rkza9asiccffzxKpVJ897vfjfb29njxxRfne1UAuC51eQdcFMXXznz66afxyiuvxMjISHz00UfR0tIS+/fvnzXX29sb5XJ5+jE+Pn4zVgaAOanLAJdKpRl3vKOjo9Hc3Dxj5vXXX4+WlpZoamqKhoaG2LlzZxw/fnzWz+rp6Ynh4eHpx9KlS2/6/gDwdeoywJ2dnVGtVmN4eDgiIg4ePBhdXV0zZlatWhVvv/12TExMRETEkSNHpj8zBoB6V5cBbmhoiL6+vtixY0e0trbGkiVLYvfu3XH48OGoVCoREXHvvffGnj174t5774177rknPvzww3jmmWeSNweA2txS1PKB6yJSLpen76wBWLgW+ut5Xd4BA8BiJ8AAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJ6jbAg4OD0dHREW1tbVGpVGJqamrWzNjYWGzfvj3uvvvuaG9vj1dffTVhUwCYu7oM8JUrV6JSqcTAwECMjIzExMRE9Pf3z5rbs2dP/PznP4+zZ8/GmTNn4r777kvYFgDmri4DPDQ0FMuXL49yuRwREd3d3XHo0KEZM+fOnYuxsbF4+OGHIyKisbEx7rjjjnnfFQCuR10GuFqtRqlUmr5euXJlnD9/fsbM2bNn484774xHHnkkNm7cGI8++mh89tln870qAFyXugxwURRfO3P58uU4duxY7Nu3L06ePBmrV6+Op59+etZcb29vlMvl6cf4+PjNWBkA5qQuA1wqlWbc8Y6OjkZzc/OsmXK5HOvWrYuIiF27dsU777wz62f19PTE8PDw9GPp0qU3d3kAqEFdBrizszOq1WoMDw9HRMTBgwejq6tr1syXX345HeojR45ER0fHvO8KANejMXuBq2loaIi+vr7YsWNH/Pe//437778/du/eHYcPH47Dhw9HX19f3HrrrfH73/8+tm/fHpcvX44VK1bESy+9lL06ANTklqKWD1wXkXK5PH1nDcDCtdBfz+vyLWgAWOwEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0CCug3w4OBgdHR0RFtbW1QqlZiamrrm7LZt26KtrW0etwOAG1OXAb5y5UpUKpUYGBiIkZGRmJiYiP7+/qvO/vnPf45ly5bN84YAcGPqMsBDQ0OxfPnyKJfLERHR3d0dhw4dmjV34cKF6O3tjX379s33igBwQ+oywNVqNUql0vT1ypUr4/z587PmnnzyyXj22Wfjtttum8/1AOCG1WWAi6L42pnXXnstGhoaYsuWLf9zrre3N8rl8vRjfHz8m1oTAK5bY/YCV1MqlWbc8Y6OjkZzc/OMmWPHjsXRo0dj9erVcfny5RgbG4v169fH6dOnZ8z19PRET0/P9PVXb2sDQKa6vAPu7OyMarUaw8PDERFx8ODB6OrqmjHz3HPPRbVajQ8++CDeeuutWLVq1az4AkC9qssANzQ0RF9fX+zYsSNaW1tjyZIlsXv37jh8+HBUKpXs9QDght1S1PKB6yJSLpen76wBWLgW+ut5Xd4BA8BiJ8AAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJBBgAEggwACQQYABIIMAAkECAASCBAANAAgEGgAQCDAAJ6jbAg4OD0dHREW1tbVGpVGJqamrG86dOnYpNmzZFR0dHrFu3Ll544YWkTQFg7uoywFeuXIlKpRIDAwMxMjISExMT0d/fP2PmO9/5Trz00ktx5syZOH78eLz44otx6tSpnIUBYI7qMsBDQ0OxfPnyKJfLERHR3d0dhw4dmjGzZs2aaG9vj4iI22+/PdauXRvnz5+f910B4HrUZYCr1WqUSqXp65UrV/7PuL733ntx4sSJ2LRp03ysBwA3rDF7gaspiqLm2YsXL8ZDDz0Uzz//fCxbtmzW8729vdHb2zt9PT4+/o3sCAA3oi7vgEul0ow73tHR0Whubp41d+nSpdi2bVs89thjsXPnzqv+rJ6enhgeHp5+LF269KbtDQC1qssAd3Z2RrVajeHh4YiIOHjwYHR1dc2YmZycjK6urnjggQfiiSeeyFgTAK5bXQa4oaEh+vr6YseOHdHa2hpLliyJ3bt3x+HDh6NSqURExMsvvxxHjhyJv//977Fhw4bYsGFD/O1vf0veHABqc0sxlw9cF4FyuTx9Zw3AwrXQX8/r8g4YABY7AQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEhQtwEeHByMjo6OaGtri0qlElNTU7NmXn755VizZk20trbGvn37ErYEgOtTlwG+cuVKVCqVGBgYiJGRkZiYmIj+/v4ZM59//nk89dRT8eabb8a5c+fijTfeiDfffDNpYwCYm7oM8NDQUCxfvjzK5XJERHR3d8ehQ4dmzPzzn/+MH/7wh9HU1BSNjY2xZ8+eWTMAUK/qMsDVajVKpdL09cqVK+P8+fNzngGAetWYvcDVFEXxjcxERPT29kZvb+/09bvvvjt9Z821jY+Px9KlS7PXWBCcVW2cU+2cVW3ef//97BVuSF0GuFQqzbibHR0djebm5lkzp0+f/p8zERE9PT3R09MzfV0ul2N4ePgmbL24OKfaOavaOKfaOavaLPSbqbp8C7qzszOq1er0P8CDBw9GV1fXjJmtW7fGG2+8EZ988klcvnw5/vSnP82aAYB6VZcBbmhoiL6+vtixY0e0trbGkiVLYvfu3XH48OGoVCoREfG9730v9u/fHz/4wQ+ivb097r///vjRj36UvDkA1KYu34KOiNiyZcust2AefPDBePDBB6evd+3aFbt27ZrTz/3/347m2pxT7ZxVbZxT7ZxVbRb6Od1S1PrXTADAN6Yu34IGgMVOgAEgwaIMsO+Rrt3XndWpU6di06ZN0dHREevWrYsXXnghadNctfyb+sq2bduira1tHrerL7Wc1djYWGzfvj3uvvvuaG9vj1dffTVh01y1nNOBAweio6Mj1q9fH1u3bo2xsbGETXPt3bs3mpubo7Hx2n+yNJffz7pSLDJTU1NFa2trcebMmaIoimLnzp3FH//4xxkzFy9eLEqlUvHxxx8Xk5OTxX333VcMDg5mrJuqlrM6d+5ccfbs2aIoiuLzzz8v2traipMnT873qqlqOaev9Pf3F7/4xS+K1tbW+VyxbtR6Vj/+8Y+Lv/71r0VRFMXk5GRx4cKFed0zWy3n9O9//7toaWkpLl26VBRFUfz2t78tnnrqqXnfNdu//vWv4pNPPikaGhqu+vxcfj/rzaK7A/Y90rWr5azWrFkT7e3tERFx++23x9q1a791X/lZyzlFRFy4cCF6e3u/1e+o1HJW586di7GxsXj44YcjIqKxsTHuuOOOed81Uy3nVBRFTE5OxhdffBFFUcTExEQ0NTVlrJtq8+bNcdddd13z+Vp/P+vRoguw75Gu3VzP4b333osTJ07Epk2b5mO9ulHrOT355JPx7LPPxm233Taf69WVWs7q7Nmzceedd8YjjzwSGzdujEcffTQ+++yz+V41VS3ntGbNmnj88cejVCpFU1NTnDlzJvbu3Tvfq9a9hfx6vugCXHyD3yO92M3lHC5evBgPPfRQPP/887Fs2bKbuFX9qeWcXnvttWhoaIgtW7bMw0b1q5azunz5chw7diz27dsXJ0+ejNWrV8fTTz89D9vVj1rO6dNPP41XXnklRkZG4qOPPoqWlpbYv3//PGy3sCzk1/NFF+Bav0f662a+DWo9h0uXLsW2bdvisccei507d87ninWhlnM6duxYHD16NFavXh2bN2+ODz/8MNavXz/fq6ar9fevXC7HunXrIuL/vlDnnXfemdc9s9VyTq+//nq0tLREU1NTNDQ0xM6dO+P48ePzvWrdW8iv54suwL5Huna1nNXk5GR0dXXFAw88EE888UTGmulqOafnnnsuqtVqfPDBB/HWW2/FqlWrZvxnId8WtZxVZ2dnfPnll9MvmkeOHImOjo553zVTLee0atWqePvtt2NiYiIi/u+cFvp/PnAz1HKWdSvtz79uoqNHjxZr164tWlpail/+8pfF5ORk8Y9//KPo7u6envnLX/5StLW1FS0tLcUzzzyTuG2urzur/v7+4tZbby2+//3vTz8GBgaSt55/tfyb+sp//vOfb+1fQRdFbWd17NixYuPGjcU999xTbN26tfj4448TN85Ryzn97ne/K9rb24t169YVP/vZz4rx8fG8hZP86le/KlasWFFERLFixYriN7/5TTE0NFT85Cc/mZ652lkuBL6KEgASLLq3oAFgIRBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAkEGAASCDAAJBBgAEggwACQQIABIIEAA0ACAQaABAIMAAn+H9SFNNSVrCuaAAAAAElFTkSuQmCC' width=480.0/>\n",
       "            </div>\n",
       "        "
      ],
      "text/plain": [
       "Canvas(toolbar=Toolbar(toolitems=[('Home', 'Reset original view', 'home', 'home'), ('Back', 'Back to previous …"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "interactive(__interact2__, i=(0, len(u1.z) - 1, 1))\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "ExecuteTime": {
     "end_time": "2020-05-02T18:59:08.658238Z",
     "start_time": "2020-05-02T18:59:08.539538Z"
    }
   },
   "source": [
    "<img src=\"interactive.png\" width=\"400\">"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.12"
  },
  "toc": {
   "base_numbering": 1,
   "nav_menu": {
    "height": "2px",
    "width": "252px"
   },
   "number_sections": true,
   "sideBar": true,
   "skip_h1_title": false,
   "title_cell": "Table of Contents",
   "title_sidebar": "Contents",
   "toc_cell": false,
   "toc_position": {},
   "toc_section_display": "block",
   "toc_window_display": false
  },
  "varInspector": {
   "cols": {
    "lenName": 16,
    "lenType": 16,
    "lenVar": 40
   },
   "kernels_config": {
    "python": {
     "delete_cmd_postfix": "",
     "delete_cmd_prefix": "del ",
     "library": "var_list.py",
     "varRefreshCmd": "print(var_dic_list())"
    },
    "r": {
     "delete_cmd_postfix": ") ",
     "delete_cmd_prefix": "rm(",
     "library": "var_list.r",
     "varRefreshCmd": "cat(var_dic_list()) "
    }
   },
   "types_to_exclude": [
    "module",
    "function",
    "builtin_function_or_method",
    "instance",
    "_Feature"
   ],
   "window_display": true
  },
  "widgets": {
   "application/vnd.jupyter.widget-state+json": {
    "state": {
     "006a137aeb83431093680f3a24a7670c": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "1.4.0",
      "model_name": "VBoxModel",
      "state": {
       "_dom_classes": [
        "widget-interact"
       ],
       "children": [
        "IPY_MODEL_d65ee933ec8a4bce83e400934ce1785a",
        "IPY_MODEL_a6470499fb564faea0661e190de0b7f8"
       ],
       "layout": "IPY_MODEL_1f6e1e518c244ec2ad6c511a1a245570"
      }
     },
     "1801440298a248e8a9689c2fa4c1c063": {
      "model_module": "@jupyter-widgets/base",
      "model_module_version": "1.1.0",
      "model_name": "LayoutModel",
      "state": {}
     },
     "1f6e1e518c244ec2ad6c511a1a245570": {
      "model_module": "@jupyter-widgets/base",
      "model_module_version": "1.1.0",
      "model_name": "LayoutModel",
      "state": {}
     },
     "2261c03730b04048a095e0633aa83099": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "1.4.0",
      "model_name": "SliderStyleModel",
      "state": {
       "description_width": ""
      }
     },
     "51d75a73448542c79218c2bf3464fe61": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "1.4.0",
      "model_name": "VBoxModel",
      "state": {
       "_dom_classes": [
        "widget-interact"
       ],
       "children": [
        "IPY_MODEL_f507918388474e55a593e4ce89070464",
        "IPY_MODEL_c175ae3db1084cefbca99f3f7c4edfc9"
       ],
       "layout": "IPY_MODEL_dbd5e35759af4870a87a88b5a4736540"
      }
     },
     "72fb652c4f5e416885c871c85a593063": {
      "model_module": "@jupyter-widgets/base",
      "model_module_version": "1.1.0",
      "model_name": "LayoutModel",
      "state": {}
     },
     "7dc29c00dbc146d6b06bdd1c6d0ff0dc": {
      "model_module": "@jupyter-widgets/base",
      "model_module_version": "1.1.0",
      "model_name": "LayoutModel",
      "state": {}
     },
     "a6470499fb564faea0661e190de0b7f8": {
      "model_module": "@jupyter-widgets/output",
      "model_module_version": "1.0.0",
      "model_name": "OutputModel",
      "state": {
       "layout": "IPY_MODEL_7dc29c00dbc146d6b06bdd1c6d0ff0dc",
       "outputs": [
        {
         "data": {
          "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAEICAYAAABRSj9aAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nO3de3zU9Z3v8dcn9xCSAEmIgXBVLgZRkIsEUMJFuVkBtRW3R2j16LGubd3teRS23V33tOtjtd2tx+q2lqOutLq14g1KAQOUhDsKgsj9JnIRMASSSUjIbT7nj5mJAXKZJDNJ5jef5+Mxj8z8rh9+Gd7zzXd+v+9PVBVjjDHOFdHeBRhjjAkuC3pjjHE4C3pjjHE4C3pjjHE4C3pjjHE4C3pjjHE4C3pjjHE4C3oT0kTkDRE5IyIuETkkIv+zkWVvEpEPReS8iDR4AYmIDBCRyyLyxlXTvy8in3v3tV1Exgfy32JMsFjQm1D3b0BfVU0C7gH+VURGNLBsFfA28EgT2/xP4OO6E0TkNuBZ4H4gGXgVeF9EIltRuzFtwoLeBI2IPCAipXUeFSKSF8h9qOpeVa3wvfQ+rm9g2YOq+iqwt5Ga5wJFwNqrZvUF9qrqDvVcTv57IBXo3sB28kTkX0Vks/ff/mcRSRGRN71/EXwsIn3rLK8i8oSIHBaREhH5uYhc713fJSJvi0iMP8fEmKtZ0JugUdU/qWpnVe0M9ACOAX+sb1kR+Y2IFDXw2N3YfrzrlgEHgDPAipbUKyJJwM+Av69n9kogUkRu87biHwZ2AWcb2eRc4CGgJ54Pny3AfwHdgP3A01ctPxUYAYwBfgwsAv4H0Au4CXiwJf8uYyzoTdCJSATw30Ceqv6uvmVU9QlV7dLA4+bGtq+qTwCJwO3Ae0BFY8s34ufAq6p6qp55JcC7wEbv9p8GHtPGB4v6L1U9qqrFeD4ojqrqGlWtBpYAw69a/heq6lLVvcAeIFdVj9VZ/+rljfGLBb1pC8/gCeIfBGsHqlqjqhuBTOB7zV1fRIYBU4DnG1jkEeC7wBAgBk9Le7mI9Ghks+fqPC+v53XnVi5vjF+i2rsA42zePu8HgVGqWtXIci/jCc/6fKGqQ/zcZRQN9NE3IQdPP/wJEQFPqEaKSJaq3goMA5ar6iHv8qtE5AwwFninBfszps1Yi94EjYgMB14EZqtqQWPLqurjvv78eh71hryIdBeRuSLSWUQiRWQqng+Vq79I9S0vIhKHp0WOiMSJSKx39iI8HxDDvI+Xgb/g6TcHz1k4M0Wkv3c7dwID8XSxGNOhWYveBNMsoCuw0dtKBtigqtMDtH3F003zMp5GyxfAU6q6DEBEegP7gCxVPQH0AT6vs365d52+qloGlPlmiEgpcLnOB9Tv8XwQ5Hn/TaeA/6WqBwL0bzEmaMRuPGKMMc5mXTfGGONwFvTGGONwFvTGGONwFvTGGONwbXrWTWpqqvbt27ctd2mMMSFvx44d51U1raXrt2nQ9+3bl+3bt7flLo0xJuSJyBetWd+6bowxxuEs6I0xxuEs6E1Y27t3L5MmTWL8+PH8/Oc/b+9yjAkKC3oT1pYuXcq6des4efIkL7zwQnuXY0xQWNCbsHbw4EEyMzN56qmnKCws5Pz58+1dkjEBZ0FvwtqBAwcYNGgQgwYNAjzBb4zTWNCbsKWqHDhwgMGDBzN48GDAE/zGOI0FvQlbZ8+exeVyMXjwYPr06UNsbKwFvXEkC3oTtnyhPnjwYCIjIxk4cKAFvXEkC3oTtnz98b7++UGDBlkfvXEkC3oTtg4cOEBCQgI9e/YEPC37Y8eOUVFR0c6VGRNYFvQmbPnOuImI8Pw3GDx4MDU1NRw9erSdKzMmsCzoTdjynXHjY2feGKeyoDdhqaysjBMnTtT2zwMMHDgQsKA3ztNk0IvIIBHZVefhEpGnRKSbiKwWkcPen13bomBjAuH06dOoKnXvj5CYmEhqaionT55sv8KMCYImg15VD6rqMFUdBowAyoD3gYXAWlUdAKz1vjYmJBQUFADQvXv3K6anpaXVzjPGKZrbdTMZOKqqXwCzgMXe6YuB2YEszJhg8oV5amrqFdNTU1Mt6I3jNDfo5wJ/9D5PV9Uz3udngfSAVWVMkPkGL0tLu/LubGlpaTawmXEcv4NeRGKAe4AlV89TVQW0gfUeE5HtIrLdWkqmo/C9F+sLenufGqdpTot+OvCJqp7zvj4nIhkA3p9f1beSqi5S1ZGqOvLq/1TGtJeCggI6depEp06drpielpZGYWEhbre7nSozJvCaE/QP8nW3DcAyYL73+XxgaaCKMibYCgoKrumfB08fvdvt5sKFC+1QlTHB4VfQi0gCcCfwXp3JzwJ3ishhYIr3tTEh4fz589d028DXXTnWT2+cJMqfhVT1EpBy1bRCPGfhGBNyCgoKGg36goKCK66aNSaU2ZWxJiz5E/TGOIUFvQlLjfXR++Yb4xQW9CbslJWVUVZWZn30JmxY0Juw09DFUgCxsbEkJiZai944igW9CTsNXSzlYxdNGaexoDdhx4LehBsLehN2fF039X0Z65tuffTGSSzoTdixFr0JNxb0JuwUFBQQFRVFly5d6p3vC3rPWH3GhD4LehN2fOfQi0i989PS0qioqKC0tLSNKzMmOCzoTdg5f/58g/3z8HXfvfXTG6ewoDdhp6HhD3xsGATjNBb0JuwUFhaSkpLS4Hxfi76wsLCtSjImqCzoTdgpLi5u8ItYgOTk5NrljHECC3oTdoqLi2vDvD4W9MZpLOhNWKmqqqKsrMyC3oQVC3oTVlwuF0CjQd+pUyciIyMt6I1jWNCbsFJUVAQ0HvQiQnJycu2yxoQ6C3oTVnyt9MaC3jffWvTGKSzoTVixoDfhyK+gF5EuIvKOiBwQkf0iki0i3URktYgc9v7sGuxijWktC3oTjvxt0b8ArFLVwcAtwH5gIbBWVQcAa72vjenQ/A36Ll26WNAbx2gy6EUkGbgDeBVAVStVtQiYBSz2LrYYmB2sIo0JFGvRm3DkT4u+H1AA/JeI7BSRV0QkAUhX1TPeZc4C6fWtLCKPich2EdluY4eY9mZBb8KRP0EfBdwK/FZVhwOXuKqbRj0Dd9c7eLeqLlLVkao6srGBpIxpC8XFxcTHxxMdHd3ocsnJybhcLhuT3jiCP0F/Cjilqtu8r9/BE/znRCQDwPvzq+CUaEzgNDX8gU9ycjJut9vGpDeO0GTQq+pZ4KSIDPJOmgzsA5YB873T5gNLg1KhMQHU1IBmPjYMgnGSKD+X+z7wpojEAMeA7+L5kHhbRB4BvgC+FZwSjQmc5rTofctnZmYGuyxjgsqvoFfVXcDIemZNDmw5xgRXUVFRs1r0NgyCcQK7MtaElZa06I0JdRb0JqxY0JtwZEFvwoq/Qe/r3rGgN05gQW/CRlVVFeXl5daiN2HHgt6EDX+vigWIj48nKirKgt44ggW9CRvNCXrfzUcs6I0TWNCbsNGcoPctZ0FvnMCC3oQNC3oTrizoTdiwoDfhyoLehA1/bgxel90g3DiFBb0JG9aiN+HKgt6EjeYGvd1O0DiFBb0JGy6Xy6+bjvgkJSXZzUeMI1jQm7BRUlJCYmKi38snJiaiqpSVlQWxKmOCz4LehI2WBL1vPWNCmQW9CRsW9CZcWdCbsGFBb8KVBb0JGxb0JlxZ0JuwYUFvwpVf94wVkeNACVADVKvqSBHpBvwJ6AscB76lqheDU6YxrVdSUkJSUpLfy1vQG6doTot+oqoOU1XfTcIXAmtVdQCw1vvamA7LWvQmXLWm62YWsNj7fDEwu/XlGBMcbrebS5cuWdCbsORv0CuQKyI7ROQx77R0VT3jfX4WSK9vRRF5TES2i8j2goKCVpZrTMuUlpYCNCvoO3fuDFjQm9DnVx89MF5VT4tId2C1iByoO1NVVUTqvU5cVRcBiwBGjhxp15KbduEL6+YEfUREBAkJCRb0JuT51aJX1dPen18B7wOjgXMikgHg/flVsIo0prVaEvS+5S3oTahrMuhFJEFEEn3PgbuAPcAyYL53sfnA0mAVaUxrWdCbcOZP10068L6I+Jb/b1VdJSIfA2+LyCPAF8C3glemMa1jQW/CWZNBr6rHgFvqmV4ITA5GUcYEmgW9CWd2ZawJCxb0JpxZ0JuwYEFvwpkFvQkLFvQmnFnQm7DgcrmIiIigU6dOzVovMTERl8sVpKqMaRsW9CYslJSU0LlzZ7xnj/ktMTGR8vJyqqurg1SZMcFnQW/CQnMHNPPxreMbQsGYUGRBb8JCa4Pe+ulNKLOgN2HBgt6EMwt6ExYs6E04s6A3YcGC3oQzC3oTFizoTTizoDdhwYLehDMLehMWLOhNOLOgN45XWVlJZWWlBb0JWxb0xvFaOs4NQHx8PBERERb0JqRZ0BvHa03Qi4gNbGZCngW9cbzWBL1vPQt6E8os6I3jWdCbcGdBbxzPgt6EO7+DXkQiRWSniCz3vu4nIttE5IiI/ElEYoJXpjEtZ0Fvwl1zWvQ/BPbXef0c8Lyq3gBcBB4JZGHGBIoFvQl3fgW9iGQCM4FXvK8FmAS8411kMTA7GAUa01oW9Cbc+dui/7/AjwG393UKUKSqvtvunAJ61reiiDwmIttFZHtBQUGrijWmJVob9ElJSRb0JqQ1GfQicjfwlaruaMkOVHWRqo5U1ZFpaWkt2YQxrVJSUkJ0dDSxsbEtWt9a9CbURfmxzDjgHhGZAcQBScALQBcRifK26jOB08Er05iWa+k4Nz6JiYlUVVVRUVHR4g8LY9pTky16Vf0HVc1U1b7AXOCvqvptYB1wv3ex+cDSoFVpTCsEIuh92zEmFLXmPPoFwN+LyBE8ffavBqYkYwLLgt6EO3+6bmqpah6Q531+DBgd+JKMCSwLehPu7MpY43gW9CbcWdAbx7OgN+HOgt44ngW9CXcW9MbxLOhNuLOgN46mqhb0JuxZ0BtHKy8vx+12W9CbsGZBbxyttePcALXDJ1jQm1BlQW8cLRBB71vfgt6EKgt642gW9MZY0BuHs6A3xoLeOJwFvTEW9MbhLOiNsaA3DmdBb4wFvXE4C3pjLOiNw/nCuXPnzq3ajgW9CWUW9MbRSkpKiI+PJyqqWbdeuEZiYiKlpaWoaoAqM6btWNAbR2vtODc+iYmJuN1uysrKAlCVMW3Lgt44WiCD3rc9Y0KNBb1xtJKSEpKSklq9HQt6E8os6I2jWYveGD+CXkTiROQjEflURPaKyP/xTu8nIttE5IiI/ElEYoJfrjHNY0FvjH8t+gpgkqreAgwDponIGOA54HlVvQG4CDwSvDKNaRkLemP8CHr1KPW+jPY+FJgEvOOdvhiYHZQKjWkFC3pj/OyjF5FIEdkFfAWsBo4CRapa7V3kFNCzgXUfE5HtIrK9oKAgEDUb4zcLemP8DHpVrVHVYUAmMBoY7O8OVHWRqo5U1ZFpaWktLNOY5qupqeHSpUsBOevGtw2Xy9XqbRnT1pp11o2qFgHrgGygi4j4LjfMBE4HuDZjWsXX+g5E0CckJCAiFvQmJPlz1k2aiHTxPo8H7gT24wn8+72LzQeWBqtIY1qiuLgYgOTk5FZvKyIigqSkpNptGhNK/BkAJANYLCKReD4Y3lbV5SKyD3hLRP4V2Am8GsQ6jWm2QAa9bzsW9CYUNRn0qrobGF7P9GN4+uuN6ZB83SyB6Lrxbce6bkwositjjWNZi94YDwt641gW9MZ4WNAbx/KFcqC6bizoTaiyoDeO5etPD1SL3vroTaiyoDeOVVxcTFRUFPHx8QHZnrXoTaiyoDeOVVxcTFJSEiISkO0lJydTUVFBRUVFQLZnTFuxoDeO5XK5AtZtAzYMggldFvTGsYqLiwMa9L5tWfeNCTUW9MaxfF03gWJBb0KVBb1xLGvRG+NhQW8cy/rojfGwoDeOZV03xnhY0BtHUlXrujHGy4LeOFJ5eTk1NTXWdWMMFvTGoQI9zg1AdHQ08fHx1qI3IceC3jhSoEeu9LFhEEwosqA3jhSsoLfbCZpQZEFvHCnQd5fySU5Otj56E3Is6I0jWdeNMV9rMuhFpJeIrBORfSKyV0R+6J3eTURWi8hh78+uwS/XGP9Y0BvzNX9a9NXAj1Q1CxgD/K2IZAELgbWqOgBY631tTIcQrK4bu/mICUVNBr2qnlHVT7zPS4D9QE9gFrDYu9hiYHawijSmuYJxeiVYi96Epmb10YtIX2A4sA1IV9Uz3llngfQG1nlMRLaLyPaCgoJWlGqM/4qLi+ncuTORkZEB3W5ycjIlJSXU1NQEdLvGBJPfQS8inYF3gadU9Yq/XVVVAa1vPVVdpKojVXVkWlpaq4o1xl+BHv7Ax/cXQklJScC3bUyw+BX0IhKNJ+TfVNX3vJPPiUiGd34G8FVwSjSm+VwuV8C7bcDGuzGhyZ+zbgR4Fdivqr+qM2sZMN/7fD6wNPDlGdMyFy5coFu3bgHfbteunpPLLl68GPBtGxMs/rToxwEPAZNEZJf3MQN4FrhTRA4DU7yvjekQCgsLSUlJCfh2U1NTa7dvTKiIamoBVd0ISAOzJwe2HGMC4/z584wYMSLg2/V9eJw/fz7g2zYmWOzKWOM4qhq0Fr1vm9aiN6HEgt44TllZGRUVFUEJel+/vwW9CSUW9MZxfCHs608PpJiYGBITEy3oTUixoDeO4wvhYLTowfMBYkFvQokFvXGcYAd9SkqKBb0JKRb0xnF8Z8QEM+jtrBsTSizojeMEs48erEVvQo8FvXEcXwgH48pYsKA3oceC3jhOYWEhSUlJREdHB2X7KSkpFBcXU11dHZTtGxNoFvTGcYJ1sZSPr0vowoULQduHMYFkQW8cJ9hBb1fHmlBjQW8c5/z580H7IhZsvBsTeizojeNYi96YK1nQmw6lurqasWPHEhcXR58+fdi6dWuzt9FRg97tdnP//feTlJREfHw8q1evDkZ5xlzDgt50KL/73e/YsmULDz74IBUVFSxcuLBZ61dVVeFyudrky9jmBv3y5ct59913mTFjBtdddx1PPvkklZWVwSjRmCtY0JsOo7CwkH/6p39i0qRJvPbaayxcuJD8/Hw2bdrk9zZ8Z8IEs4++U6dOxMbGNivoVZVnnnmGfv368cYbb/DSSy9x6NAhXnzxxaDVaYyPBb3pMJ577jlcLhcvvPACIsKjjz5KamoqzzzzjN/bCPY4NwAi0uyLptauXctHH33EggULiIqKYubMmcyYMYOf/exnlJaWBq1WY8CC3nQQbrebN998k5kzZ3LTTTcBkJCQwA9/+ENWrlzJiRMn/NpOsMe58WnueDe/+93vSE9P5zvf+U7ttAULFuByuVi+fHkQKjTmaxb0pkPYuHEjX375JQ8++OAV0+fMmQPAmjVr/NpOW7Tofdv3t0VfU1PD2rVrmTlzJrGxsbXTx48fT48ePXjrrbeCVaYxgB9BLyKvichXIrKnzrRuIrJaRA57f3YNbpnG6d566y06derEN77xjSumZ2Vlcd1113W4oE9NTfW7Rb9z504uXrzIlClTrpgeERHBAw88wMqVKykqKgpGmcYA/rXoXwemXTVtIbBWVQcAa72vjWmR6upq3nnnHe6++24SEhKumCciTJkyhTVr1uB2u5vc1unTpxERrrvuumCVC0CPHj348ssv/VrW9yE1adKka+bNnTuXyspKPvjgg4DWZ0xdTQa9qq4Hrh7UYxaw2Pt8MTA7wHWZMJKfn09BQQEPPPBAvfOnTJlCQUEBe/bsqXd+XSdOnCAjI4OYmJhAl3mF3r1743K5KC4ubnLZNWvWcPPNN5Oenn7NvFGjRtGvXz/efvvtYJRpDNDyPvp0VT3jfX4WuPYd7CUij4nIdhHZXlBQ0MLdGSdbtWoVMTExTJ06td75kydPBvzrpz9x4gS9evUKaH318e2jqS+Jy8vL2bhx4zXdNj4iwqxZs1i3bh3l5eUBr9MYCMCXsaqqgDYyf5GqjlTVkWlpaa3dnXGgDz/8kPHjx1/TbeOTmZnJ4MGD/Q763r17B7rEa/j20VTQb9q0iYqKigaDHmDatGlcvnyZ9evXB7RGY3xaGvTnRCQDwPvzq8CVZMLJmTNn+Oyzz7jrrrsaXW7ChAls3ry50X56VeXkyZNtGvQnT55sdLmNGzcSERHB+PHjG1zmjjvuIC4ujpUrVwa0RmN8Whr0y4D53ufzgaWBKceEm9zcXIAGu218srOzKS4uZv/+/Q0uU1hYSHl5eZt03aSnpxMVFdVki37Lli3cdNNNJCYmNrhMfHw8OTk5rFq1KtBlGgP4d3rlH4EtwCAROSUijwDPAneKyGFgive1Mc2Wm5tL9+7dufnmmxtdLjs7G/AEZ0N8odsWLfrIyEgyMzMbDXq32822bdtqa2/MtGnTOHjwIJ9//nkgyzQG8O+smwdVNUNVo1U1U1VfVdVCVZ2sqgNUdYqq2q12TLO53W5Wr17NnXfeSURE42/FAQMGkJKS0mGC3refxoJ+//79FBcX+xX006dPBzzfVxgTaHZlrGk3u3btoqCgoMluG/CcnTJmzJhGg97XX96WQd9YH72vVn+CfsCAAfTr18+6b0xQWNCbduPrn2/sjJS6srOz2b9/PxcvXqx3/okTJ4iLiwvqyJV19erVi1OnTlFTU1Pv/C1btpCSksKAAQOa3JaIMG3aNNauXWtDF5uAs6A37SY3N5ebb76ZjIwMv5b3tYwbuhmJ7xx6EQlYjY3p3bs31dXVnD17tt75mzdvZsyYMX7XM23aNEpLS5s1LLMx/rCgN+2itLSUjRs3NnlaZV2jR48mIiKCzZs31zu/rc6h92nsXPoLFy5w4MABv7ptfCZOnEh0dLR135iAs6A37SI/P5+qqiq/+ud9OnfuzIgRI8jPz693fludQ+/T2Ln0voufJkyY4Pf2EhMTGT9+vAW9CTgLetMucnNziYuLa/RCovrk5OSwbds2ysrKrpheVVXFl19+2Sbn0Ps0NgxCXl4e8fHxjBo1qlnbnDZtGrt37/Z7wDRj/GFBb9qcqrJixQpycnKIi4tr1roTJ06ksrLymrNvjh8/jqrSp0+fQJbaqOTkZLp06cKRI0eumbdu3TrGjh17xfjz/rDTLE0wWNCbNnfo0CGOHDnC3Xff3ex1x40bR2RkJHl5eVdM932BedtttwWiRL+NHj36mu8MCgsL2b17Nzk5Oc3e3k033USPHj1sOAQTUBb0ps35bp3XkqBPSkpixIgR1wT9+vXrSUlJ4cYbbwxEiX674447+Oyzz2pvSu6rBTx/fTSX7zTL1atXU11dHbA6TXizoDdtbvny5QwdOrTF3SwTJ068pp9+/fr13H777U1eYRtod9xxB+AZvMynpf3zPtOmTaOoqIiPPvooIDUaY0Fv2lRRUREbNmxoUWveZ9KkSVRVVbFixQrAc1epo0eP1oZuWxo1ahSxsbG1rfiamhqWLVvG7bff3uKbn0yZMoWIiAg7+8YEjAW9aVOrVq2ipqamVUE/efJk+vbty69//WsANmzYANAuQR8XF8dtt91WG/R//vOfOX78OI8++miLt9m1a1fGjBlT+0FmTGtZ0JsWUVWqqqqavd7vf/97evbs2aovTSMjI/n+97/Phg0b2LlzJ+vXrycxMZFbbrmlxdtsjQkTJvDJJ59QUlLCCy+8QO/evZk9u3V315wzZw47duxg3759zVrP7XZb3765hgW98VtVVRWvv/46WVlZREdHEx8fz9/93d/5fQu8EydOsGrVKh5++GEiIyNbVcvDDz9MQkICf/M3f8PixYsZP348UVFRrdpmS+Xk5FBTU8PEiRPJy8vjySefbHUt8+bNIzo6mldeecXvdf7617/Sv39/oqOjSU5O5vHHH+f48eOtqsM4hKq22WPEiBFqQofb7daKigo9dOiQvvzyy9q3b18FdPjw4fqTn/xE58+fr4DedNNNeu7cuSa39/TTT6uI6PHjxwNS349+9CONjo7Whx56SD///POAbLMl3G63vvjii9q/f3/t1q2bXrhwISDbvf/++zUlJUUvX77c5LLPPvusAjpw4EB9+umndf78+RoTE6NRUVH68MMP65YtW/TChQtaU1MTkNpM2wK2ayuy14LeqKpqdXW17tmzR5955hmdNGmSpqam+u4FXPsYPXq0/vnPf1a321273sqVKzU+Pl5HjBihLper0e336tVLp06dGtCaS0tLA7a91gp0PR9++KEC+tZbbzW63Msvv6yAzp07Vy9dulQ7/eTJk/qDH/xA4+Lirvg9xsTE6M0336yPPPKILl26tNHfm+kYWhv04tlG2xg5cqRu3769zfbX0ZSWlrJ9+3YOHjxIeno6t956K5mZmUREROByuXC5XKgqmZmZDY54WF5ezsGDB7l8+TKpqan069ePiIgITp8+zZo1azhw4ACXLl1i6NChzJkzh6tvyH769GneeOMNVq1axdmzZykrKyMmJobTp0/XdsHceuutDB8+nJ49exIdHU1mZiZZWVmMGjWq3rr+8pe/MGvWLCZOnMhf/vKXes82efHFF/nBD37A+++/3+r+63DhdrsZOHAgUVFR7Ny5k/j4+GuWee+99/jmN7/JtGnT+OCDD4iOjr5mmXPnzrFlyxaOHj1KaWkppaWl7Nu3j02bNlFcXAx4hnOIiopCVUlNTWXIkCHMmzePnJycK05ZrampYcOGDaxYsYLq6mrS0tKYMGECI0eOJCYmhosXL3LkyBFEhOuuu47MzMwG/31FRUW4XC5iY2NJTU0lMjKSkpIS9uzZw+7du0lJSWH48OFcf/31ATiaoU1EdqjqyBav7/Sgd7vdHDp0iK1bt3Ls2DGys7OZMGECnTp18nsbx44dY/HixWzbto1jx44xYsQIZs+ezb333lvvf6yrFRUV8etf/5rnn3+eoqKiK+b5+rpdLlfttF69enHvvfdy3333MXbsWESE9evX84c//IF33nnnimVjYmKoqampHRM9OjqauLg4SkpKiIyMZNKkSUyePBkRYd26deTm5uJ2uxkxYgT9+/enU6dOVFRUkJGRwdChQ7nrrrvo2bOn38fGZ/HixXznO99h7ty5vPnmm1eEw7Fjxxg6dCh33HEHK1asaLNhhJ0gNzeXqVOn8n/5yn8AAAxdSURBVOMf/5jnnnvuinn5+flMnTqV4cOHs2bNGhISEpq17crKSvLz8/n44485cOAAviw4f/48mzdvxuVy0bt3b+bOnUv37t35/PPPeeeddzh37hwxMTFER0dz6dKl2u3Fx8df831NdnY2Dz30EN/61rdISUnhyy+/5P333+fdd98lPz+/9mbvUVFRJCQk1H7w1DV9+nSefvppv7/A37dvH2+++SYrV64kIiKCoUOHcu+99zJ9+nS/vztRVfbv38/q1aspKSlhyJAhjBs3ju7du/u1fqCFfNBXVlayc+dOtmzZUvs4f/48GRkZ3HDDDWRlZTFkyBCGDBlCVlYWycnJje6jqKiIbdu2sXXrVrZu3cq2bduuuVFF165d+clPfsITTzxxTeBfunSJ3bt3c/DgQfbs2UNeXh47duwgIiKCW265hd69e7Nlyxa++uorevXqxXe/+13GjRtX23KuG2KnT59m06ZN/OY3v6G4uJh77rmHxx9/nKysLM6cOcOuXbv44osvuHTpEr169aJr165UVlby4Ycf8uGHH1JRUQF4wryyspLOnTtz3333MWPGDBITEzl79iwHDhwgOjqa7t27k5OTw9ChQwH49NNPWbJkCW+//XbtWCy9evVi/vz5zJs3z6+bYTTXL37xCxYsWMDUqVN5+eWX6du3L/v37+ehhx7i8OHD7Nmzp00HHXOKRx99lNdee42XXnqJxx57DICXXnqJn/70p/Tu3ZsNGzaQkpIS0H2Wl5ezdOlSXn/9dVavXo3b7SY+Pp6ZM2fyzW9+k5kzZ5KQkMCFCxdYu3Zt7W0Tr7vuOgYNGoSIsG/fPt544w327NlDREQEERERtWcEDR48mPvuu4/+/ftTXl7Ol19+icvlIjMzk0GDBjFs2DAuXLjAqlWr+NWvfkVhYSF33XUX3/72txk0aBBxcXFf9z/jOVHg008/ZdmyZSxfvpzIyEjGjRtHbGwsn3zyCYWFhaSlpTF58mSys7NrMyU9Pf2ahsfWrVtZuHDhNaOkigijRo1i7NixjBo1ilGjRnHDDTc02nBRVQ4dOsSOHTs4duwYR48erf1ZWFhIv379uOWWW5g4cSITJ05scHshF/TLly9ny5YtbN68mS1btrBjxw4uX74MQJ8+fcjOziYjI4MzZ85w6NAh9u/ff0UrYeDAgUydOpWpU6cyYcIEEhISOH78OGvWrGHJkiWsXbsWt9uNiDBkyBCys7MZM2YMY8aMoU+fPmzatInnn3+eVatW0blzZ+6++25SU1O5ePEin3zyCQcPHqxtZURHR5Odnc3UqVOZN29e7Z+hbrebVatW8ctf/pL8/HyaOoazZ8/mn//5nxk+fLjfx6qkpIQVK1awb98+ysrKuPXWW5k1a1az/hIBzxutpKSEiIgIEhISgtqaVlV++9vfsmDBAsrLy8nIyODs2bN07tyZ1157jTlz5gRt305WXFzM7NmzycvLIz09nbKyMkpKSpg+fTqvvPIKPXr0COr+L1++TGVlJbGxsc0epE1V2b17N++99x5VVVWkpKQwffp0srKy/N5GaWkpv/nNb/j3f/93CgoKGl02PT2d733vezzxxBO1jS/fxXVLlixh3bp1V4wM6useGjBgACJCfn4+e/fupXv37ixcuJB7772X1NRU9uzZQ25uLrm5uezYsaM2k3r37s0DDzzAzJkzGTNmDFFRURw9epS8vDzWrVtHXl7eFTem6dGjB9dffz39+/cnJSWFY8eO8dFHH9XW1LNnT8aNG8fYsWMZO3Ysw4YNIzo6un2DXkSmAS8AkcArqvpsY8vHxsaq7zZpMTExjBgxgrFjx5KdnU12dna9b1i3283x48fZu3cve/fuZf369eTl5dV7Sl///v2ZO3cukyZNYtSoUSQlJTVYy8aNG3n99ddZsWIFFRUVJCQkMGzYsNr+6aysLPr27dtk14zL5eLjjz+mtLS0NvB9LY3u3buTlZVFt27dGt2G05w8eZJFixZx8uRJUlNTWbBgwTXfFZjmUVWWLFnCBx98QGpqKjk5OcyZMyesusGqqqo4evQohw8frv3LQERqHzfeeKNfLexz587V5slnn33Gjh07+OKLL1BVsrKyePDBB5k3bx6JiYn1bqO6upp9+/axbds2PvjgA3Jzc+u9diEjI4OJEyeSk5NDdnY2119/fb3fs/ha/b4Phs2bN9fe4yA+Pp7Ro0eTn5/fPkEvIpHAIeBO4BTwMfCgqjZ4hUe3bt30H//xHxk7dizDhw9vduvA5/Lly2zYsIGPP/64tpUwadIkbrzxxrB64xtj2l9xcTF5eXns2rWr9kvonJyc2r8SWuLUqVNs3ryZzZs3s2nTJrZv395uQZ8N/IuqTvW+/gcAVf23htYJ97NujDGmJVrbddOay/d6AnXvoXYKuOZrcRF5DHjM+7JCRPa0Yp9tJRU4395F+CEU6gyFGsHqDDSrM7AGtWbloF8zrqqLgEUAIrK9NZ9KbcXqDJxQqBGszkCzOgNLRFrVFdKasW5OA3XPlcv0TjPGGNOBtCboPwYGiEg/EYkB5gLLAlOWMcaYQGlx142qVovIk8CHeE6vfE1V9zax2qKW7q+NWZ2BEwo1gtUZaFZnYLWqzja9YMoYY0zbs/HojTHG4SzojTHG4YIS9CLyTRHZKyJuERl51bx/EJEjInJQRKbWmT7NO+2IiCwMRl1N1PwnEdnlfRwXkV3e6X1FpLzOvJfburar6vwXETldp54ZdebVe2zbqc5fisgBEdktIu+LSBfv9A51PL01tet7rz4i0ktE1onIPu//pR96pzf4+2/HWo+LyGfeerZ7p3UTkdUictj7s2s71ziozjHbJSIuEXmqIxxPEXlNRL6qe41RQ8dPPH7tfa/uFpFb/dpJawazb+gB3IjnBP88YGSd6VnAp0As0A84iueL3Ejv8/5AjHeZrGDU5mf9/wH8s/d5X2BPe9VST23/AvzveqbXe2zbsc67gCjv8+eA5zro8exQ7706dWUAt3qfJ+IZbiSrod9/O9d6HEi9atovgIXe5wt9v/+O8PD+zs8CfTrC8QTuAG6t+/+ioeMHzABWAgKMAbb5s4+gtOhVdb+qHqxn1izgLVWtUNXPgSPAaO/jiKoeU9VK4C3vsm1OPINTfAv4Y3vsvxUaOrbtQlVzVdU30tNWPNdZdEQd5r1Xl6qeUdVPvM9LgP14rkYPFbOAxd7ni4GOdLeZycBRVf2ivQsBUNX1wIWrJjd0/GYBv1ePrUAXEcloah9t3Udf37AJPRuZ3h5uB86p6uE60/qJyE4RyReR29uprrqe9P7Z9lqdP4k70jG82sN4WiE+Hel4duTjBni6u4DhwDbvpPp+/+1JgVwR2SGeIU8A0lX1jPf5WSC9fUqr11yubMh1tOMJDR+/Fr1fWxz0IrJGRPbU82j31lBD/Kz5Qa58E5wBeqvqcODvgf8WkYbHPw5+nb8FrgeGeWv7j2DW0oo6fcv8FKgG3vROavPjGcpEpDPwLvCUqrroQL//Osar6q3AdOBvReSOujPV0+fQIc7jFs/FnfcAS7yTOuLxvEIgjl9rLpia0oLVGhs2IejDKTRVs4hEAfcCI+qsUwFUeJ/vEJGjwEAgaMNw+ntsReT/Acu9L9t8SAo/jud3gLuByd43a7sczyZ02KE8RCQaT8i/qarvAajquTrz6/7+242qnvb+/EpE3sfTHXZORDJU9Yy3a+Grdi3ya9OBT3zHsSMeT6+Gjl+L3q9t3XWzDJgrIrEi0g8YAHxExxlOYQpwQFVP+SaISJp4xt5HRPp7az7WDrX56qnbHzcH8H1T39CxbRfiuSnNj4F7VLWszvQOdTzpOO+9K3i/K3oV2K+qv6ozvaHff7sQkQQRSfQ9x/Ml/B48x3C+d7H5wNL2qfAaV/zF3tGOZx0NHb9lwDzv2TdjgOI6XTwNC9K3yHPw9B1VAOeAD+vM+ymesxwOAtPrTJ+B58yCo8BP2+nb79eBx6+adh+wF9gFfAJ8oz1qq1PPH4DPgN3eX3pGU8e2neo8gqcvcZf38XJHPJ7emtr9vVdPTePx/Lm+u84xnNHY77+d6uyP50ylT72/1596p6cAa4HDwBqgWwc4pglAIZBcZ1q7H088HzxngCpvbj7S0PHDc7bNf3rfq59R56zGxh42BIIxxjicXRlrjDEOZ0FvjDEOZ0FvjDEOZ0FvjDEOZ0FvjDEOZ0FvjDEOZ0FvjDEO9/8Ba8g6jWfZ9b4AAAAASUVORK5CYII=\n",
          "text/plain": "<Figure size 432x288 with 1 Axes>"
         },
         "metadata": {
          "needs_background": "light"
         },
         "output_type": "display_data"
        }
       ]
      }
     },
     "c175ae3db1084cefbca99f3f7c4edfc9": {
      "model_module": "@jupyter-widgets/output",
      "model_module_version": "1.0.0",
      "model_name": "OutputModel",
      "state": {
       "layout": "IPY_MODEL_e9ffcbaefa8b4588b38d900caa2865c8",
       "outputs": [
        {
         "data": {
          "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAEICAYAAABRSj9aAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAcCUlEQVR4nO3de3Bc533e8e+DOwkQNxIEAdIkSPGiS1VRNiQ5tmSnlh3bamLJqePKTROm1YzG47i1m3RcJU7SpIlnbHeS1O00cRhLNZPxNb5EqmfsWFGt2rGrC2RRF5IiwStICAQBgrjfCOLXP/Zgu4QAAgQWBHT6fGbO7Dnvvgfnh3fPPnv24OxCEYGZmaVXwXIXYGZmS8tBb2aWcg56M7OUc9CbmaWcg97MLOUc9GZmKeegNzNLOQe9ve5J+qikFkljkr44j/7/TtJZSf2SHpFUmrRvljQ4bQpJv5mz7r+QdErSkKS/lVS7hL+aWV446C0NXgX+CHhkro6S3g08BNwNbAG2AX8AEBFtEVExNQE3A5PAN5N1bwL+AvgVoB4YBv4s77+NWZ456G1JSfrn046QxyQ9mc9tRMS3IuJvgfPz6L4HeDgiDkTEBeAPgV+bpe+vAj+MiJPJ8i8D/zMifhgRg8DvAr8oac1MKyfvBj4iqVXSgKQ/lHSdpJ8k7ya+Lqkk6fuzks5I+oSkc5I6JN0n6R5JRyT1SPrteQ+KWQ4HvS2piPhazhFyI3Ac+MpMfSX9maTeWaYX81TSTcALOcsvAPWS1k6rRWSCft9s60bEMWAc2HmF7b0beBPwZuATwF7gXwJvAP4R8KGcvhuAMmAj8HvAXyZ93wTcBfyupK3z/D3Nshz0dk1IKgC+DDwZEX8xU5+I+EhEVM8y/eM8lVIB9OUsT81PPyq/k8zpmW9cYd2p9Wc8ok98NiL6I+IA8DLw/Yg4HhF9wHeBW3P6XgQ+FREXga8C64DPRcRAsv5B4Ja5fkGz6Rz0dq18ikwg/ttlrmMQqMxZnpofmNZvD/DN5BTNbOtOrT993VydOfMjMyxX5Cyfj4hLOffNtH5uf7N5cdDbkpN0P5lTFB9IjlZn6/f5Ga56mZoO5KmcA1x+VHwL0BkR2fP7klYBv8Tlp21es66kbUApcCRPtZktCQe9LSlJtwL/DbgvIrqu1DciPpx71cu06aYrbKNIUhlQCBRKKpNUNEv3vwIekHSjpGrgd4AvTuvzfuAC8INp7V8CfkHSXZLKgf8EfCsirnREb7bsHPS21O4FaoB/yDk6/26et/E7ZE5rPETmj5cjSVvutfGbASLie8BnyYR4G3AK+I/Tft4e4K9j2j9rSM6Tf5hM4J8jcyrqI3n+XczyTv7HI2Zm6eYjejOzlHPQm5mlnIPezCzlHPRmZik32yVoS2LdunXR1NR0LTdpZva699xzz3VHRN1C17+mQd/U1ERLS8u13KSZ2euepFOLWd+nbszMUs5Bb2aWcg56M7OUc9CbmaWcg97MLOUc9GZmKeegNzNLOQe9mVnKOejNzFLOQW9mlnIOejOzlHPQm5mlnIPezCzlHPRmZik3Z9BL2iVpf87UL+njkmolPS6pNbmtuRYFm5nZ1Zkz6CPicETsjojdwJuAYeDbwEPAExGxA3giWTYzsxXmak/d3A0ci4hTwL3AvqR9H3BfPgszM7P8uNqgvx/4SjJfHxEdyfxZoD5vVZmZWd7MO+gllQDvA/5m+n0REUDMst6DkloktXR1dS24UDMzW5irOaJ/L/DTiOhMljslNQAkt+dmWiki9kZEc0Q019Ut+H/bmpnZAl1N0H+I/3faBuAxYE8yvwd4NF9FmZlZ/swr6CWVA+8CvpXT/GngXZJagXcmy2ZmtsIUzadTRAwBa6e1nSdzFY6Zma1g/mSsmVnKOejNzFLOQW9mlnIOejOzlHPQm5mlnIPezCzlHPRmZinnoDczSzkHvZlZyjnozcxSzkFvZpZyDnozs5Rz0JuZpZyD3sws5Rz0ZmYp56A3M0s5B72ZWco56M3MUs5Bb2aWcg56M7OUm1fQS6qW9A1Jr0g6JOlnJNVKelxSa3Jbs9TFmpnZ1ZvvEf3ngO9FxPXALcAh4CHgiYjYATyRLJuZ2QozZ9BLqgLeBjwMEBHjEdEL3AvsS7rtA+5bqiLNzGzh5nNEvxXoAv6HpOclfUFSOVAfER1Jn7NA/UwrS3pQUouklq6urvxUbWZm8zafoC8C3gj8eUTcCgwx7TRNRAQQM60cEXsjojkimuvq6hZbr5mZXaX5BP0Z4ExEPJ0sf4NM8HdKagBIbs8tTYlmZrYYcwZ9RJwFTkvalTTdDRwEHgP2JG17gEeXpEIzM1uUonn2+zfAlySVAMeBf0XmReLrkh4ATgEfXJoSzcxsMeYV9BGxH2ie4a6781uOmZnlmz8Za2aWcg56M7OUc9CbmaWcg97MLOUc9GZmKeegNzNLOQe9mVnKOejNzFLOQW9mlnIOejOzlHPQm5mlnIPezCzlHPRmZinnoDczSzkHvZlZyjnozcxSzkFvZpZyDnozs5Rz0JuZpdy8/mespJPAAHAJmIiIZkm1wNeAJuAk8MGIuLA0ZZqZ2UJdzRH9P4mI3REx9U/CHwKeiIgdwBPJspmZrTCLOXVzL7Avmd8H3Lf4cszMLN/mG/QBfF/Sc5IeTNrqI6IjmT8L1M+0oqQHJbVIaunq6lpkuWZmdrXmdY4euDMi2iWtBx6X9ErunRERkmKmFSNiL7AXoLm5ecY+Zma2dOZ1RB8R7cntOeDbwO1Ap6QGgOT23FIVaWZmCzdn0Esql7Rmah74OeBl4DFgT9JtD/DoUhVpZmYLN59TN/XAtyVN9f9yRHxP0rPA1yU9AJwCPrh0ZZqZ2ULNGfQRcRy4ZYb288DdS1GUmZnljz8Za2aWcg56M7OUc9CbmaWcg97MLOUc9GZmKeegNzNLOQe9mVnKOejNzFLOQW9mlnIOejOzlHPQm5mlnIPezCzlHPRmZinnoDczSzkHvZlZyjnozcxSzkFvZpZyDnozs5Rz0JuZpdy8g15SoaTnJX0nWd4q6WlJRyV9TVLJ0pVpZmYLdTVH9B8DDuUsfwb404jYDlwAHshnYWZmlh/zCnpJm4B/CnwhWRbwDuAbSZd9wH1LUaCZmS3OfI/o/wvwCWAyWV4L9EbERLJ8Btg404qSHpTUIqmlq6trUcWamdnVmzPoJf08cC4inlvIBiJib0Q0R0RzXV3dQn6EmZktQtE8+rwVeJ+ke4AyoBL4HFAtqSg5qt8EtC9dmWZmtlBzHtFHxG9FxKaIaALuB/5XRPwy8APgA0m3PcCjS1almZkt2GKuo/8PwG9IOkrmnP3D+SnJzMzyaT6nbrIi4kngyWT+OHB7/ksyM7N88idjzcxSzkFvZpZyDnozs5Rz0JuZpZyD3sws5Rz0ZmYp56A3M0s5B72ZWco56M3MUs5Bb2aWcg56M7OUc9CbmaWcg97MLOUc9GZmKeegNzNLOQe9mVnKOejNzFLOQW9mlnIOejOzlHPQm5ml3JxBL6lM0jOSXpB0QNIfJO1bJT0t6aikr0kqWfpyzczsas3niH4MeEdE3ALsBt4j6c3AZ4A/jYjtwAXggaUr08zMFmrOoI+MwWSxOJkCeAfwjaR9H3DfklRoZmaLMq9z9JIKJe0HzgGPA8eA3oiYSLqcATbOsu6DkloktXR1deWjZjMzuwrzCvqIuBQRu4FNwO3A9fPdQETsjYjmiGiuq6tbYJlmZrZQV3XVTUT0Aj8AfgaollSU3LUJaM9zbWZmlgfzueqmTlJ1Mr8KeBdwiEzgfyDptgd4dKmKNDOzhSuauwsNwD5JhWReGL4eEd+RdBD4qqQ/Ap4HHl7COs3MbIHmDPqIeBG4dYb242TO15uZ2QrmT8aamaWcg97MLOUc9GZmKeegNzNLOQe9mVnKOejNzFLOQW9mlnIOejOzlHPQm5mlnIPezCzlHPRmZinnoDczSzkHvZlZyjnozcxSzkFvZpZyDnozs5Rz0JuZpZyD3sws5Rz0ZmYpN2fQS3qDpB9IOijpgKSPJe21kh6X1Jrc1ix9uWZmdrXmc0Q/AfxmRNwIvBn4dUk3Ag8BT0TEDuCJZNnMzFaYOYM+Ijoi4qfJ/ABwCNgI3AvsS7rtA+5bqiLNzGzhruocvaQm4FbgaaA+IjqSu84C9bOs86CkFkktXV1diyjVzMwWYt5BL6kC+Cbw8Yjoz70vIgKImdaLiL0R0RwRzXV1dYsq1szMrt68gl5SMZmQ/1JEfCtp7pTUkNzfAJxbmhLNzGwx5nPVjYCHgUMR8Sc5dz0G7Enm9wCP5r88MzNbrKJ59Hkr8CvAS5L2J22/DXwa+LqkB4BTwAeXpkQzM1uMOYM+Iv4B0Cx3353fcszMLN/8yVgzs5Rz0JuZpZyD3sws5Rz0ZmYp56A3M0s5B72ZWco56M3MUs5Bb2aWcg56M7OUc9CbmaWcg97MLOUc9GZmKeegNzNLOQe9mVnKOejNzFLOQW9mlnIOejOzlHPQm5mlnIPezCzl5gx6SY9IOifp5Zy2WkmPS2pNbmuWtkwzM1uo+RzRfxF4z7S2h4AnImIH8ESybGZmK9CcQR8RPwR6pjXfC+xL5vcB9+W5LjMzy5OFnqOvj4iOZP4sUD9bR0kPSmqR1NLV1bXAzZmZ2UIt+o+xERFAXOH+vRHRHBHNdXV1i92cmZldpYUGfaekBoDk9lz+SjIzs3xaaNA/BuxJ5vcAj+anHDMzy7f5XF75FeD/ALsknZH0APBp4F2SWoF3JstmZrYCFc3VISI+NMtdd+e5FjMzWwL+ZKyZWco56M3MUs5Bb2aWcg56M7OUc9CbmaWcg97MLOUc9GZmKeegNzNLOQe9mVnKOejNzFLOQW9mlnIOejOzlHPQm5mlnIPezCzlHPRmZinnoDczSzkHvZlZys35H6Zs5ZicnKSg4LWvzRMTExw6dIjDhw8zMjLCxo0buf7662loaEDSMlRqaTYyMsLBgwdpa2ujr6+PrVu3cvPNN1NbW/uavhEB4P1wmV3ToG9vb+dTn/oUBQUFRER2J5h+K4ni4mKKioooLi6muLiYkpISKioqKCgoYGxsjPHxccbHxxkbG6OsrIytW7dSX19PeXk558+f59VXX6Wjo4OOjg46OztpbGzkxhtvZOPGjQAcPnyYw4cPc+LECYqLi6mvr+euu+7ijjvuoKGhgdHRUV555RV+9KMf0drayujoKI2Njdx2223ccMMN1NXV0dnZyaFDh3jhhRfo7+9HEtu3b2fHjh1s2LCBS5cu0d3dTXd3N11dXXR3dzM5OcnatWspKChgaGiI06dPc/bsWYaHh6moqGDz5s1s3ryZdevW0d3dzaFDh3jqqac4ceIEFy5coKGhgeuuu47rrruOoqIi2tvb+clPfkJ/f/9rxnvt2rVs2LCB8vLy7FiNj49TWlrK6tWrWb16NQD9/f309fUxODhIcXExZWVllJWVUVpaetl8aWkpFRUVbNiwgYaGBtatW8fIyAh9fX309vbS19fHwMAA1dXVbNy4kU2bNlFRUUFPTw8nTpygtbWVvr4+ALZv387111/Prl27qKmpoaCggIGBAfr6+rK/S3l5eXZavXr1a+aLi4uBTPBcuHCB3t5eLly4wODgIEVFRdmpuLiYiKCjo4Ouri7Gx8epqalh586d7Nixg6qqKiDzgnny5ElaW1s5deoUABUVFWzatIm1a9eyevVqxsfHGRwcZHBwkIigpqaG6upqampqqKyszO7bIyMjDA8PMzQ0xPDwcHYaGhpicHCQoaGhy8a6rKyM4uJiTp8+zeHDh3nllVc4f/485eXlbNq0iZ07d1JfX09JSQkdHR2cOXOGV199lYmJCUpKSqitrWX9+vXU1dVx6dIlenp6LptGRkayz62IyD5mvb29SKKqqorKykpWr17NxMQE4+PjXLx4kfHxcSYmJli1ahUAPT09tLW1MTk5edm+JolbbrmFG264gaqqKs6cOcOxY8c4ceIEk5OTNDY2csstt3D77bezceNGJNHW1sapU6fo6OiguLiYqqoqtmzZwrp16yguLqa/v5/e3l4igvLy8uxzo6Ghgb6+Pk6fPs3Bgwfp7e2luLiY7du3s3v3bjZv3kxBQQFHjx7l+eef56WXXso+Vm95y1u47bbbWL9+PSMjIxw4cIAnn3ySgwcPMjg4SE1NDbt27WLnzp00NjZy7tw5jh07xsGDBxkZGaG4uJimpqbsvlNaWkpPTw+nTp3i1VdfZWRkBElUVlZSWVlJeXk5Y2NjjI6OMjIywsjICKOjo1y8eJGJiYnsdOnSJSYmJogISktLWbVqFatWreLSpUszPrevlqbC9VqQdO02ligoKGDt2rV0d3cz/XedetAmJyezDxLAqlWrsvMA1dXVlJWV0dnZ+ZqfAZmdfNWqVdknyNVYs2YNjY2NlJeX09/fT1tb22U/o6KigjvuuIPrr7+empoa2tvbOXbsGMeOHWNycpL6+npuu+023v72t3PTTTdRXl7O6dOnOXDgAC+//DLd3d0MDQ1RWlpKSUkJJSUljI+PZ0MnIrJP8oqKCiYmJhgdHc3unLnTxYsX6evr4+zZs4yOjl72exQVFVFVVcWaNWu4cOFCNtCnFBYW0tTURG1tLZcuXeLIkSMMDg5e1VhNV1xcnH3hX4yqqiqKioro6+tjYmJiwT9HEmVlZZftOwtRUFBAU1MT69evZ2hoKHvknKuoqIiGhgZKS0sZHR3l/Pnzr9luSUkJa9eupba2NhvUEZEN9urqaqqrq4mI7Iv98PBw9sBq6rawsJDR0VEmJyepqalh69at7N69m23btrFmzRqOHz/OM888w5NPPsnJkyfp7e1l48aN2QOSwsJC2traePbZZzl69OhlNdbX19PY2Jh9cWpvb7/sOVZUVIQkLl68OOt4TR3IzNanqqqKkpISLly4kH18pz/HGxsbqayszB6YzfQzKisrGRsb49y5c7PWcjWmDmYLCwuzt5IYGxtjZGQkW2t5eTlDQ0PPRUTzQre1qKCX9B7gc0Ah8IWI+PSV+jc3N8ePf/zj3PWzb+lybycnJ5mYmODixYvZV77R0dFsMJWUlGSDq7S0lKGhIU6cOEFXVxcDAwOsW7eOhoYGGhoaWL9+PYWFhQwPD3P06FE6OjqYnJxk586dbNmyhaKizJua8fFxnnnmGV544QWOHTtGTU0N27dv561vfSubN28GYHBwkP3799Pa2kp3dzcbNmxg27Zt7N69m/LyciYnJ2lra+P48eN0dnZSVFTEunXrstPUkXxPTw8AZWVlrFmz5rK3tZOTk5w7d47z589n1yssLFzwY7QUpoKhu7ub1atXU1VVxapVqy77PQYHB2lvb2d4eJiamhoaGxspKSm57Gd0dHRw5MgRBgYGmJiYYM2aNVRVVWWPsIeGhl4zTb1ATU1T4ZM7lZeXZ4+QpqbJycns/lBaWkpXVxetra20trbS1tZGRFBZWcmOHTvYsWMHW7dupaCggP7+fs6cOUNPTw/Dw8OUlZVRUVFBeXk5QPZdxNTt8PBw9t3S1LuP6fNT60/t11PT2NhYNiDLysouG6upAJp6Z1lXV/ea03hDQ0N0dXVRVFSUDfeVdspkYGAge8C0adOm7AvQlPHxcfr7+7l48WL2HYYkRkZGOHPmDG1tbXR0dGTfNe7cuTP73Dtx4gQvvvgi7e3tTExMsG3bNm6++WaampqyP+Opp57ipZde4uTJk9TW1rJr1y7uvPNOGhoasjX09PRw5MgROjs7Wb9+PVu2bLnsNOjAwABHjx6ltbWVS5cuUV1dzebNm9m0aVO2loGBAfr7+xkcHMweoZeVlWVvp3LnSiYmJigoKKCgoABJyxP0kgqBI8C7gDPAs8CHIuLgbOs0NzdHS0vLgrZnZvb/q8UG/WKuurkdOBoRxyNiHPgqcO8ifp6ZmS2BxfwxdiNwOmf5DHDH9E6SHgQeTBbHJL28iG1eK+uA156oW3leD3W+HmoE15lvrjO/di1m5SW/6iYi9gJ7ASS1LObtx7XiOvPn9VAjuM58c535JWlR57wXc+qmHXhDzvKmpM3MzFaQxQT9s8AOSVsllQD3A4/lpywzM8uXBZ+6iYgJSR8F/o7M5ZWPRMSBOVbbu9DtXWOuM39eDzWC68w315lfi6rzmn5gyszMrj1/qZmZWco56M3MUm5Jgl7SL0k6IGlSUvO0+35L0lFJhyW9O6f9PUnbUUkPLUVdc9T8NUn7k+mkpP1Je5OkkZz7Pn+ta5tW5+9Las+p556c+2Yc22Wq8z9LekXSi5K+Lak6aV9R45nUtKz73kwkvUHSDyQdTJ5LH0vaZ338l7HWk5JeSuppSdpqJT0uqTW5rVnmGnfljNl+Sf2SPr4SxlPSI5LO5X7GaLbxU8Z/TfbVFyW9cV4bmfoWyXxOwA1kLvB/EmjOab8ReAEoBbYCx8j8Ibcwmd8GlCR9blyK2uZZ/x8Dv5fMNwEvL1ctM9T2+8C/n6F9xrFdxjp/DihK5j8DfGaFjueK2vdy6moA3pjMryHzdSM3zvb4L3OtJ4F109o+CzyUzD809fivhCl5zM8CW1bCeAJvA96Y+7yYbfyAe4DvAgLeDDw9n20syRF9RByKiMMz3HUv8NWIGIuIE8BRMl+lsGK+TkGZby76IPCV5dj+Isw2tssiIr4fEVNfBfkUmc9ZrEQrZt/LFREdEfHTZH4AOETm0+ivF/cC+5L5fcB9y1jLdHcDxyLi1HIXAhARPwR6pjXPNn73An8VGU8B1ZIamMO1Pkc/09cmbLxC+3K4C+iMiNactq2Snpf0vyXdtUx15fpo8rbtkZy3xCtpDKf712SOQqaspPFcyeMGZE53AbcCTydNMz3+yymA70t6TpmvPAGoj4iOZP4sUL88pc3ofi4/kFtp4wmzj9+C9tcFB72kv5f08gzTsh8NzWaeNX+Iy3eCDmBzRNwK/AbwZUmVy1jnnwPXAbuT2v54KWtZRJ1TfT4JTABfSpqu+Xi+nkmqAL4JfDwi+llBj3+OOyPijcB7gV+X9LbcOyNzzmFFXMetzIc73wf8TdK0EsfzMvkYv8V8YOqdC1jtSl+bsORfpzBXzZKKgF8E3pSzzhgwlsw/J+kYsBNYsu9bnu/YSvpL4DvJ4jX/Sop5jOevAT8P3J3srMsynnNYsV/lIamYTMh/KSK+BRARnTn35z7+yyYi2pPbc5K+TeZ0WKekhojoSE4t5Oe/dSzee4GfTo3jShzPxGzjt6D99VqfunkMuF9SqaStwA7gGVbO1ym8E3glIs5MNUiqU+a795G0Lan5+DLUNlVP7vm49wNTf6mfbWyXhTL/lOYTwPsiYjinfUWNJytn37tM8reih4FDEfEnOe2zPf7LQlK5pDVT82T+CP8ymTHck3TbAzy6PBW+xmXv2FfaeOaYbfweA341ufrmzUBfzime2S3RX5HfT+bc0RjQCfxdzn2fJHOVw2HgvTnt95C5suAY8Mll+uv3F4EPT2v7Z8ABYD/wU+AXlqO2nHr+GngJeDF50BvmGttlqvMomXOJ+5Pp8ytxPJOaln3fm6GmO8m8XX8xZwzvudLjv0x1biNzpdILyeP6yaR9LfAE0Ar8PVC7Asa0HDgPVOW0Lft4knnh6QAuJrn5wGzjR+Zqm/+e7KsvkXNV45UmfwWCmVnK+ZOxZmYp56A3M0s5B72ZWco56M3MUs5Bb2aWcg56M7OUc9CbmaXc/wVJEQweXZp7IAAAAABJRU5ErkJggg==\n",
          "text/plain": "<Figure size 432x288 with 1 Axes>"
         },
         "metadata": {
          "needs_background": "light"
         },
         "output_type": "display_data"
        }
       ]
      }
     },
     "d65ee933ec8a4bce83e400934ce1785a": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "1.4.0",
      "model_name": "IntSliderModel",
      "state": {
       "description": "i",
       "layout": "IPY_MODEL_72fb652c4f5e416885c871c85a593063",
       "max": 4095,
       "style": "IPY_MODEL_d7ebb3d9d7a942fb9f83e331af2bbaf9",
       "value": 3223
      }
     },
     "d7ebb3d9d7a942fb9f83e331af2bbaf9": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "1.4.0",
      "model_name": "SliderStyleModel",
      "state": {
       "description_width": ""
      }
     },
     "dbd5e35759af4870a87a88b5a4736540": {
      "model_module": "@jupyter-widgets/base",
      "model_module_version": "1.1.0",
      "model_name": "LayoutModel",
      "state": {}
     },
     "e9ffcbaefa8b4588b38d900caa2865c8": {
      "model_module": "@jupyter-widgets/base",
      "model_module_version": "1.1.0",
      "model_name": "LayoutModel",
      "state": {}
     },
     "f507918388474e55a593e4ce89070464": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "1.4.0",
      "model_name": "IntSliderModel",
      "state": {
       "description": "i",
       "layout": "IPY_MODEL_1801440298a248e8a9689c2fa4c1c063",
       "max": 4095,
       "style": "IPY_MODEL_2261c03730b04048a095e0633aa83099",
       "value": 1095
      }
     }
    },
    "version_major": 2,
    "version_minor": 0
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}