Skip to content

Instantly share code, notes, and snippets.

@fonnesbeck
Created August 12, 2019 18:21
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save fonnesbeck/21f2377b26f421c6355f798bf4a4b5dd to your computer and use it in GitHub Desktop.
Save fonnesbeck/21f2377b26f421c6355f798bf4a4b5dd to your computer and use it in GitHub Desktop.
StudentT sampling
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"%load_ext autoreload\n",
"%autoreload 2"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"%matplotlib inline\n",
"import pymc3 as pm\n",
"import arviz as az\n",
"import numpy as np\n",
"import seaborn as sns"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD5CAYAAADcDXXiAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAOqUlEQVR4nO3df6ykVX3H8fenrGjFENS9Wtxlu5hsq4TWaG4AS2KNa40/KPCHNrjWbpFmY0IVq1ZA0+A/TTA2KomNyQa0awSUoA3E+AsRY/sH295FrYXVQtEsV1b3GsUfNSlZ/faPeba9Xe7euzPPzJ2ds+9XQmae85yZ5zsJ+7nnnjnPuakqJElt+Y1pFyBJGj/DXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQRvW6pDko8BFwKGqOrdrez/wx8DjwH8Cl1fVY925a4ErgF8Bb62qL651jY0bN9bWrVtH/QySdFLat2/fj6pqbqVzWWude5KXAL8APr4s3F8BfKWqDid5H0BVXZ3kHOBW4DzgOcCXgd+pql+tdo35+flaWFgY8mNJ0sktyb6qml/p3JrTMlX1NeDHR7V9qaoOd4f3Apu755cAn6yq/66q7wIPMQh6SdI6Gsec+5uAz3fPNwGPLDu32LVJktZRr3BP8h7gMHDzkaYVuq0475NkV5KFJAtLS0t9ypAkHWXkcE+yk8EXrW+o/5u4XwTOWtZtM/DoSq+vqt1VNV9V83NzK34fIEka0UjhnuSVwNXAxVX1y2Wn7gQuS/LkJGcD24B/6V+mJGkYx7MU8lbgpcDGJIvAdcC1wJOBu5IA3FtVb66q+5PcBjzAYLrmyrVWykiSxm/NpZDrwaWQkjS8XkshJUmzx3CXpAatOecunehu2XvgmOd2nL9lHSuRThyO3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ3asFaHJB8FLgIOVdW5XdszgE8BW4HvAX9SVT9JEuAG4NXAL4E/r6r7JlO6WnPL3gPHPLfj/C1jfc9R30+aFWuGO/APwIeBjy9ruwa4u6quT3JNd3w18CpgW/ff+cBHukepl9WCX9ITrTktU1VfA358VPMlwJ7u+R7g0mXtH6+Be4Ezkpw5rmIlScdn1Dn3Z1fVQYDu8Vld+ybgkWX9Fru2J0iyK8lCkoWlpaURy5AkrWTcX6hmhbZaqWNV7a6q+aqan5ubG3MZknRyGzXcf3hkuqV7PNS1LwJnLeu3GXh09PIkSaMYNdzvBHZ2z3cCdyxr/7MMXAD89Mj0jSRp/RzPUshbgZcCG5MsAtcB1wO3JbkCOAC8ruv+OQbLIB9isBTy8gnULElaw5rhXlWvP8ap7Sv0LeDKvkVJkvrxDlVJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUG9wj3JXyW5P8m/J7k1yVOSnJ1kb5IHk3wqyanjKlaSdHxGDvckm4C3AvNVdS5wCnAZ8D7gg1W1DfgJcMU4CpUkHb++0zIbgN9MsgF4KnAQeBlwe3d+D3Bpz2tIkoY0crhX1feBvwMOMAj1nwL7gMeq6nDXbRHYtNLrk+xKspBkYWlpadQyJEkr6DMt83TgEuBs4DnAacCrVuhaK72+qnZX1XxVzc/NzY1ahiRpBRt6vPblwHeragkgyWeAPwDOSLKhG71vBh7tX6ZacsveA9MuQWpenzn3A8AFSZ6aJMB24AHgHuC1XZ+dwB39SpQkDavPnPteBl+c3gd8q3uv3cDVwNuTPAQ8E7hpDHVKkobQZ1qGqroOuO6o5oeB8/q8rySpH+9QlaQGGe6S1CDDXZIaZLhLUoMMd0lqUK/VMtKsWu1Gqh3nb1nHSqTJcOQuSQ0y3CWpQYa7JDXIcJekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QG9Qr3JGckuT3Jt5PsT/LiJM9IcleSB7vHp4+rWEnS8ek7cr8B+EJVPQ94AbAfuAa4u6q2AXd3x5KkdTRyuCc5HXgJcBNAVT1eVY8BlwB7um57gEv7FilJGk6fkftzgSXgY0m+nuTGJKcBz66qgwDd47PGUKckaQh9wn0D8CLgI1X1QuC/GGIKJsmuJAtJFpaWlnqUIUk62oYer10EFqtqb3d8O4Nw/2GSM6vqYJIzgUMrvbiqdgO7Aebn56tHHToB3bL3wLRLkE5qI4/cq+oHwCNJfrdr2g48ANwJ7OzadgJ39KpQkjS0PiN3gLcANyc5FXgYuJzBD4zbklwBHABe1/MakqQh9Qr3qvoGML/Cqe193leS1I93qEpSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1KC+NzFJzVlt64Qd529Zx0qk0Tlyl6QGGe6S1CDDXZIaZLhLUoMMd0lqkOEuSQ1yKaQ0BJdJalY4cpekBhnuktQgw12SGmS4S1KDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSg3qHe5JTknw9yWe747OT7E3yYJJPJTm1f5mSpGGMY+R+FbB/2fH7gA9W1TbgJ8AVY7iGJGkIvbb8TbIZeA3wt8DbkwR4GbCj67IHeC/wkT7X0Ylpte1vJU1X35H7h4B3Ab/ujp8JPFZVh7vjRWDTSi9MsivJQpKFpaWlnmVIkpYbOdyTXAQcqqp9y5tX6Forvb6qdlfVfFXNz83NjVqGJGkFfaZlLgQuTvJq4CnA6QxG8mck2dCN3jcDj/YvU5I0jJFH7lV1bVVtrqqtwGXAV6rqDcA9wGu7bjuBO3pXKUkayiTWuV/N4MvVhxjMwd80gWtIklYxlj+QXVVfBb7aPX8YOG8c7ytJGo13qEpSgwx3SWqQ4S5JDTLcJalBhrskNchwl6QGGe6S1KCxrHOXtPoumTvO37KOlUiO3CWpSYa7JDXIaRmtyj/IIc0mR+6S1CDDXZIaZLhLUoMMd0lqkOEuSQ1ytYy0Do616sibmzQpjtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDTLcJalBhrskNWjkcE9yVpJ7kuxPcn+Sq7r2ZyS5K8mD3ePTx1euJOl49Bm5HwbeUVXPBy4ArkxyDnANcHdVbQPu7o4lSeto5HCvqoNVdV/3/OfAfmATcAmwp+u2B7i0b5GSpOGMZVfIJFuBFwJ7gWdX1UEY/ABI8qxjvGYXsAtgyxZ3xtPJabW/UeuOkeqj9xeqSZ4GfBp4W1X97HhfV1W7q2q+qubn5ub6liFJWqZXuCd5EoNgv7mqPtM1/zDJmd35M4FD/UqUJA2rz2qZADcB+6vqA8tO3Qns7J7vBO4YvTxJ0ij6zLlfCLwR+FaSb3Rt7wauB25LcgVwAHhdvxIlScMaOdyr6p+BHOP09lHfV5Oz2pd3ktriHaqS1CDDXZIaZLhLUoMMd0lqkOEuSQ0y3CWpQYa7JDXIcJekBo1lV0hJ4+eOkerDkbskNciRe2PcYkASOHKXpCYZ7pLUIMNdkhpkuEtSgwx3SWqQ4S5JDXIppHSS8Kaok4sjd0lqkCP3GeSNSlrNKP9/OKpvjyN3SWqQI/cpc8QkaRIcuUtSgxy5r4NR58idW5c0KsNdmkH+4NdanJaRpAZNbOSe5JXADcApwI1Vdf2krnUicCQlTdax/o258GBlExm5JzkF+HvgVcA5wOuTnDOJa0mSnmhSI/fzgIeq6mGAJJ8ELgEeGPeFHDFLkzXuf2OzMNIe928J01jyPKk5903AI8uOF7s2SdI6mNTIPSu01f/rkOwCdnWHv0jynQnVsh42Aj+adhE9zHr94Gc4ERxX/W8Y80XH/H6rfoZx1z6G9/ztY52YVLgvAmctO94MPLq8Q1XtBnZP6PrrKslCVc1Pu45RzXr94Gc4Ecx6/dDGZzhiUtMy/wpsS3J2klOBy4A7J3QtSdJRJjJyr6rDSf4S+CKDpZAfrar7J3EtSdITTWyde1V9DvjcpN7/BDPr00uzXj/4GU4Es14/tPEZAEhVrd1LkjRT3H5AkhpkuI9ZkncmqSQbp13LMJK8P8m3k/xbkn9Mcsa0azoeSV6Z5DtJHkpyzbTrGVaSs5Lck2R/kvuTXDXtmkaR5JQkX0/y2WnXMookZyS5vfs3sD/Ji6ddU1+G+xglOQv4I2AWb5u9Czi3qn4f+A/g2inXs6ZGtrk4DLyjqp4PXABcOYOfAeAqYP+0i+jhBuALVfU84AXM9mcBDPdx+yDwLo66YWsWVNWXqupwd3gvg3sTTnT/u81FVT0OHNnmYmZU1cGquq97/nMGoTJTd3Mn2Qy8Brhx2rWMIsnpwEuAmwCq6vGqemy6VfVnuI9JkouB71fVN6ddyxi8Cfj8tIs4Dk1tc5FkK/BCYO90KxnahxgMan497UJG9FxgCfhYN7V0Y5LTpl1UX/6xjiEk+TLwWyuceg/wbuAV61vRcFarv6ru6Pq8h8FUwc3rWduI1tzmYlYkeRrwaeBtVfWzaddzvJJcBByqqn1JXjrteka0AXgR8Jaq2pvkBuAa4G+mW1Y/hvsQqurlK7Un+T3gbOCbSWAwpXFfkvOq6gfrWOKqjlX/EUl2AhcB22s21siuuc3FLEjyJAbBfnNVfWba9QzpQuDiJK8GngKcnuQTVfWnU65rGIvAYlUd+Y3pdgbhPtNc5z4BSb4HzFfVzGwC1f1xlQ8Af1hVS9Ou53gk2cDgy9/twPcZbHuxY5buhs5gNLAH+HFVvW3a9fTRjdzfWVUXTbuWYSX5J+Avquo7Sd4LnFZVfz3lsnpx5K4jPgw8Gbir++3j3qp683RLWl0j21xcCLwR+FaSb3Rt7+7u8Nb6eQtwc7cX1sPA5VOupzdH7pLUIFfLSFKDDHdJapDhLkkNMtwlqUGGuyQ1yHCXpAYZ7pLUIMNdkhr0P31gtWksVdPAAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"t_data = pm.StudentT.dist(nu=5).random(size=1000)\n",
"sns.distplot(t_data, kde=False);"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Using the Student-T directly"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"Auto-assigning NUTS sampler...\n",
"Initializing NUTS using jitter+adapt_diag...\n",
"Multiprocess sampling (2 chains in 2 jobs)\n",
"NUTS: [ν]\n",
"Sampling 2 chains, 0 divergences: 100%|██████████| 6000/6000 [00:03<00:00, 1603.59draws/s]\n"
]
}
],
"source": [
"with pm.Model() as m0:\n",
" ν = pm.HalfNormal('ν', sigma=5)\n",
" t = pm.StudentT('t', nu=ν, observed=t_data)\n",
" \n",
" trace0 = pm.sample(1000, tune=2000, cores=2)"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAA2gAAACYCAYAAACRbMq5AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOydd7gc1X33P2d29/amjiSqKDZgmimmmthxA+MSijsY44J7S2K/iRPbsZPYca8BNzAucQdjsOkGLFQQkpAE6u1e3d737t2+M3PeP6bPzq6upNsknc/zgO7uzs6cqfv7nl8TUkoUCoVCoVAoFAqFQjHzaDM9AIVCoVAoFAqFQqFQWCiBplAoFAqFQqFQKBSzBCXQFAqFQqFQKBQKhWKWoASaQqFQKBQKhUKhUMwSlEBTKBQKhUKhUCgUilmCEmgKhUKhUCgUCoVCMUtQAk2hUCgUCoVCoVAoZglKoCkUCoVCoVAoFArFLEEJNIVCoVAoFAqFQqGYJSiBplBMIkKIK4UQUgjxOt97pwshCkKIf5jJsSkUCoVCMRmo3zqFYmpRAk2hmFyeAkYB90dLSrkV+D/g8zM0JoVCoVAoJhP1W6dQTCFKoCkUk4iU0gAeAq4RQgjfRw8BZwshTpmZkSkUCoVCMTmo3zqFYmpRAk2hmHzuBxYD5/veS9v/HjP9w1EoFAqFYtJRv3UKxRShBJpCMfk8ABj4Qj+AlwEmsG1GRqRQKBQKxeSifusUiilCSClnegwKxRGHEGI50CilfLEQ4lhgC3CflPLtMzw0hUKhUCgmBfVbp1BMDcqDplBMDfcD5wkhjgd+AiSBj8/skBQKhUKhmFTUb51CMQUoD5pCMQUIIc4Engc2AycDV0gp187sqBQKhUKhmDzUb51CMTUogaZQTBFCiN3A8cA/SCnvn+nxKBQKhUIx2ajfOoVi8lECTaFQKBQKhUKhUChmCSoHTaFQKBQKhUKhUChmCUqgKRQKhUKhUCgUCsUsQQk0hUKhUCgUCoVCoZglKIGmUCgUCoVCoVAoFLMEJdAUCoVCoVAoFAqFYpYQP5CF58+fL0888cQpGopCoVAoFBbr1q0bklIumO7tqt85hUKhUEwH1X7nDkignXjiiaxdq/oPKhQKhWJqEUJ0zMR21e+cQqFQKKaDar9zKsRRoVAoFAqFQqFQKGYJB+RBUygASA9Cug9Keahvg9bjIFE306NSKBQKxSRS0A1SOZ0FzbUzPRSFQqE4qlACTbF/pIR9q2DD/8GeJ2CsM/i5FodFZ8JpV8FZ18P8U2dkmAqFQqGYPFbuHiaVK/H6c5YghJjp4SgUCsVRgxJoiursXQ6Pfh6610JNE5zyCnjJ+6HtOIjXQ24UBrdZAu7J/4Envwynvhr+7tOw9PyZHr1CoVAoDpJUrjTTQ1AoFIqjEiXQFNHkU/DQv8KzP4eWpXDNN+HsN0NNY+XvjPfB+p/B0z+AH70cznsH/P3noWnaC7EpFAqF4hARQiClREpQDjSFQqGYPpRAU5TTvwV+/TZIdsDln4ArPw2J+v1/r/kYuPJTcPEH4MmvwOrbYMt98PJ/gwvfA5qqSaNQKBSHG3KmB6BQKBRHGcpiVgTZ/iD8+BVQysLNf4FXfH5i4sxPbTO86ovwwVVw7PnwwD/Dz98IY91TMWKFQqFQTCFSKommUCgU04kSaAqPjb+2PGfzT4X3PQknXHJo65t/Krzjbnjdt6FrLdx2CTz3+8kZq0KhUCimFCeqUckzhUKhmF5UiKPC4rnfwz3vhxMvh7f+yvKCHSRSSjpHcqxpH2F7X4q9Q2cSb/42H01+lTP+8G4e//OvePjEf+KkJQu5eNk8zlzSSkxTCQ4KhUIxGzkcHWiZgk4iplETV/PQCoXi8EMJNAXsfATuuRVOuBTe/rsDD2nEEmXrOka5f1Mvj2zppzuZA6AmrrFsfiMLmo/ntrnf59VDP+Xq5C85adsW3r/xo/y3PJ75TTW84dylvOXC4zh10cELQ4VCoVAoAB7d2k9dIsarzzxmpoeiUCgUB4wSaEc7+1bDb26EhWdYnrMDFGf5ksHv1nZy16oOdg2kqY1rXHHqAt5/5TIuOmkepyxsCnnHLoQ9b+LEu9/LA7HPselFn+b29JX8bFU7P3lqL1e96Bg+8crTOE0JNYVCoZhRhAAkyMM0yDFfMmZ6CAqFQnFQKIF2NDOwDf7vTdC61MoVq2ud8FfTBZ1fru7gR8v3MpQucM6xrXzlurN57dmLaazdz2W17Ep4/wrEH9/PORu/wG2nv57RT36DO9aNcOeKdh7e0s87LzmRj7/yVFrqEoe4kwqFQqE4FA7HEEeF4ojCuQlVv4ujBiXQjlZySfj1WyFWCzfeM+FeZUXd5Ger2vnuX3cxlitx+Snz+dDLzuPiZXMRB/LgaFoAb/sdrPoePPYfzOnZwD9efwfvuuxlfO3h7dy5ci8PPt/L1244h0tPmX9w+6hQKBSKI5KBVJ5Ve4Z5xemL9j8pqFAc7ux6DPJJOOv6mR6JYppQ2bNHI6YBf3g3JDvhzT+HtuMn9LXHtw3wmm/9jf/881bOPa6Nez54Kb94z0u45OR5BybOHDQNLvso3PKQVS7sztcw99n/5b/fcCZ3f+BS6hIx3vbjp/nKg9swTDWFq1AoFNOJsOs4zkYPWudoFoDRbHGGR6JQTAP55EyPQDHNKIF2NPLYF2DXo3D1V+H4i/e7eOdIlpvvXMO7fvoMCLjz5gu565aLOO/4OZMznmMvgFuXwwtfC49+Dn55PefN1fnzR6/gLRcex/8+sZub71zDWLY0OdtTKBQKxYSZjTlophvxpUK+jgbSBZ17N3QzOF6Y6aEcXZgmdK+DUm6mR3LUoQTa0cbme2DFt+CCW+CCd1Vd1DQlP12xl1d982+sbR/l3157Og9+7KW87IULJ39c9W1ww11wzTehYwXcfhn1Xcv58nVn86Vrz2L1nmHe9INV9I3lJ3/bCoVCoSjD0T6z0YPmpuTM7DAU04QjzHqSSihMK+O9MLIXejbM9EiOOpRAO5oY3g33fgSOewm85n+qLrp3KMObf7iKz9+3hQtPmstDn3gp77li2dT2lBHCEo7v/atVsORnb4THvshbz1/CT991Ed3JHNfdtpJdA+mpG4NCoVAoAsxCfYa0FZp2FHvQ8iUD8ygJ/zdME0D1TJ12jo7razaiBNrRQikPv3snxOJw/R0Qr4lcTErJb5/p5OpvL2d73zhfvf5s7nrXhSxtO/DeaAfNojPhfU/AeW+H5V+Du67hsgV5fv2+iynoJtffvpL1+0anbzwKhUKhmFWYqqgdD23uY91R8luo2yc8EVNm67SiqkfOGOpKP1p4+DPQ9xy88XZoPTZykXRB5+O/2cCn/rCJc49r4+FPXMkNFxw3MzH+NY3whu/DtT+2xn3bZbxofAV3f+BSWusTvP1HT7N6z/D0j0uhUChmMXuHMty7oXtSPCvOk1/OwhhHJy/uaLUbnXNytIT86Ya1v4ezB21n/zh/2tgz08M4QGbfvX+0oATa0cDmP8IzP4ZLPwIveE3kIs93j3HNd5Zz38Ye/vGVp/GL97yEY1rrpnmgEZx9A9z6N5hzAvz6rRy/5j/4/XvO59g59dzy02d4pn1kpkeoUCgUs4ZtvSnA8zhMBrPSRHNz0MoN9tkoKCeboySy0cW5nuNTKdCyI7Dv6SlLutzSmzr8rk13vIevMD5cUQLtSGesC/70UTj2Qvj7z0Uu8sdnu7nutpUUdJNfv+8SPvL3p86uWap5J8O7H4GXfACevp0Fv72GX1+/gGNa6rj5jjUq3FGhUChCTGblxdloU1YLcZyN4/XTk8zROZI9pHUcba1nnBy0Kd3rfathrBOKmancyuHJ0eqqnkGUQDuSMU344wdBGnDtjyCWCHysGyb/9ectfPw3GzjnuDbu/8jlXHTS3Bka7H6I18JVX4a3/AqS+5j3i1dy9xXdzG+u5Z0/WcPGTtUjRKFQKCYT1yabhVrADXGM/Gxy2DUwzp7ByS9K9Uz7yCFPLJqzXYVOMoalz6bWAyUck9jeRikPXWut3rFHLUfXdTabUALtSGbND2Dvk/Dq/4a5JwU+GsuWeNdPn+FHy/dy0yUn8Mv3vIR5TbUzNNAD4IVXw/ufgmPOou2BD/LACb/mmAaDm+5Yw1Y7tEehUChmBHuWf9oY64bMUOCto2Wiu5oDabKM+M09KZ7rHpuUdU02R5k+c5nS3Xb7Stj3cd8mGG23vGpHK26RECUXpht1xI9UBrbBI5+D066CF98U+Kh3LMcNP1jJ6j3DfPnas/jCG150eFVGaj0W3nk/vPSfadjyG/5S91leFOvixp+sYe+QCk1QKBQzwHgfbL7bymOZLvatgj1PTPlmZmOjakeERY1sMkY723OFjFk+vqliWjxo4byrSd7mrLy2ipkKoZ2zcKxHCYeRVa6YMHoR7nkf1DbD678TmFLdNTDOdf+7kp5knrtuuYi3XHT8DA70EIjF4eX/BjfdS6KU4ufy//EPxkO840erj5qqVgqFYhYx3mv9O50CrQqTYwOKSVzX5FJtSJMx3tm4z36OthBHh8ne7aBYCnnQZqM7eqrO+/YHrP8qMguPxRGOEmhHIk/+D/RuhNd9G5oWum+v3zfK9bevomhIfv2+i7n05PkzOMhJYtmV8P4VaCdezmfkD/ls/ivc+uPHGEoXZnpkCoXiqOTIM5xn4x65HrSIwU2Gx2827rOfo6VBdZjJ3O11HSPBsvci7DGbmiTMQ9JYUoJesCbipwM3xHF6NqfwUALtSGPf0/DUN+Dcd8Dp17hvP75tgLf9aDWt9Qn+8IFLeNHS1hkc5CTTtADe/nt45Rd4pbaW28Y/zn/94OeM5UozPTKFQqGYZibfkpqykCwpwdAP+qsH89nE11++kmzx4MY6FXhVLI8Oy9kLaZ28a7FrNBRt44Y4hoqCTHaI46F+e+t9sPVPkzSaEP778SDvzQmR6oXkUZzbNwGUQDuSKKThnlutHK3XfMl9++71XbznZ2s5ZWETv3//pZwwr3EGBzlFaBpc9jG0Wx5kXlMtX0l9ige+93GyuUMrZaxQKBQTY3YZypNhU7oOhUNfVTTtT8GWPx7UV539m6r8OL+nRkpJdzLHI1v6GRjPT8n2DpSjrcy+s7dTWofHueCdqo2z65a2mOrQ1pJtM2WGrXsz5XgYJ/lgdKyAzqcnd51HGEqgHUk8/Bmr4tA//ADqWgD41Zp9fPK3G7l42Vx+9d6LWdB8GFRqPBSOu5D6D69g4PireUvmFwx+4zKKXRtmelQKheJoYTIMqP7N0Lvp0IZxOPRBS/cf9Fer5WBNigfNd/xMCaMZK6QslZsdXjTHozQbNcRUYE6BB62MMg/a1BzdQ/JIyymuFGva13fOzqV17tGJeGqTnWVVZQ9npJQzOhGiBNqRwo6HYN1P4bKPwgmXAvCbZ/bxL3c/x8tesIA7br6Q5rpE9XUcKdS3sfTdv2DFBd+hvjhC7Mcvx3zk86r5pEKhODwY2ApDOw7qq2VpNEcosuwP/2eHtvOmKdnaO+69noUHM8puNE3JvRu66RiO/q3LFHTu3dDNSGaa8pcmE8djOrV19u2NhETQQW60byxPUS8XVIcc4ngA6IZJQT+APm4VBWB1gZYrGsh9q6elqux0sbknxf2bemZMpCmBdiSQGYJ7PwyLXgQv+wwAv1vbyf+7+zmuPG0Bt73jfGrjsRke5PRz2TXv5JGX3cs9xmVoK76J/N6FsPmeI99yUSgUCCEWCyHuEkIMCiHyQogtQogrp3CDU7bqg2EynnJeiYTZ98ycyhy00WwxIHJm40+GU2bff9k5723uie4JOjhuFc/qHDn8Qv/NKRRorkfL8aC5IY4HXyQkXzJ4eu8wz7RPclXXAzwAj24d4MHn+wDYN5zdv9gwD9xDnC8ZPLylj64jrIJ216h1n5SMae5vaaME2uGOlHDfxyCfhGt/CPFa7l7fxaf+sInLT5nPD248n7rE0SfOHN7+d+fS83ff4LrC5+gr1cPvboafXgPtK2Z6aAqFYooQQrQBK7A0xmuB04GPAAMzOa5JZT+G2qQW9piFAmUq+6CFjdhZ6UGrYmhXGu/s24uJ40wSTMVkgVe00QlxdAzyg590cdaZKZQLnkO7nA7sy473rG8sz7Odo2zrixbvLpUEWpUJqHzJ2sb4LCrMtmLXEOs6Rg9pHU4Bnpm6/ZVAO9zZ+CvYdj+8/N9h0Zk8uqWff/rdRi49eR4/uumCo1qcOXzk5afw4stfw+Wjn+PRkz5lhQ799Gq463XQsXKmh6dQKCafTwG9UsqbpJRrpJR7pZSPSSm3Tv2mp+nXfD+5KJM5iinfo4Oo/OCMKbLMfjWLyihB3/NVt3k4NIE+lBy8WebsnRCOHp2KaDN3lWJyQxwnSr5ksKEzObHWCQeZg+YItaiQS6Dce1i+QMV1OxMamjZ7LqyhdMH1gFVDN8yKFb9j9vWgT2llmsoogXY4M9oBf/kUnHAZXPIh1nWM8KH/W89ZS1v54Y1KnDkIIfjXq0/nTRedyHu2nstPzr8HXv3fMLAN7rwKfvwK2PzHKg8mhUJxmPFG4GkhxG+EEANCiA1CiA+LKa1JPs3GSQVDbTJHMX0zyNNYxnxgCwxug2R7xUUOCw+a257KO+PSFTG+8RbSbs+syW6X0DmS5d4N3RR1k/ahDMns1OW2uUVCpuBceMcrXMXx4EMcq3n6wp9t7hmjYzhDb6pyhVAppTXOg9x/53rRKj0CNdtePAgPmrvugxrZNDCyF7b9ufzYGTpr9wzwxPaBSHGsOZp1ZvTZ7D2eiv1gGvDHD1h/v/E2dgxmueWna1nSVs8dN19IY218Zsc3yxBC8J9vPItrzl7MFx/u4Mf6VfCxjXDVVyEzCL97J3znPFh9u/WDplAoDmeWAR8E9gCvBr4NfBn4UNTCQoj3CSHWCiHWDg4OHtqWp8uY358HbVIjHKd4nw6hMl3U2Kruu7OtKlZX+KNZqM9coaKbJg9v7qNzJEsxKldmx4Ow98nAW/smKQdt16D1W5krGmzsSvLkjkO8d6oxSR40KSWPbe2nd8zLl8oWjFBeXkiwFcbhud9Dfj/hgT4ObJzWdqp50Lb0pli/b/Sg7xVHhFYUaMIWaAexflfgzlbXbDEDpVz5vm39E4nt9wHRkzDaRDxoY91T9oBQAu1wZdX3rT4SV3+FHrGQd96xhpq4xs9uuYh5TUd4Kf2DJKYJvvnmc7n6rGP4zz9v5fZVvfCS98FH1sObfg7Nx8CDn4ZvngGPft5qpKhQKA5HNGC9lPJfpJTPSinvBL5DBYEmpfyhlPICKeUFCxYsmNAGMgXdzb2YEfZnSB2szWAargfB9R9MtUA5CKOwugjz/izvW1YhjM1HOMRxNnvQAHIlg/X7RstzbpxGw/kxwDsshiknp5JjSMdUpYqRmy8Zkblaga+75+DQzkXRMEkXdDZ2Jt33ntw5aIkfh7DgSO4L/jsBDiQE1YkMrHaZ5dxnzcHtv9uWodK52p8HbSIhjrNEoJUVQrHbJoxl82zoTHpeWN8zIEobN6Z2U5cbwND16PYByX2wbxUM7wq8PZwusGr3MAXdsAqMZEesvnIHiBJohyN9z8NfvwgvvIbkqddx0x1rSOd17nrXRRw3t2GmRzerScQ0vvOW83j9OUv48gPb+O5jO60H0xmvh3c/DO9+BE66ElZ8G751Ftzzfuh7bqaHrVAoDoxeYEvova3A8ZO1gUe39vPQ5j7vjek2TsICIzcKesELSzxYQ3bzPVY4kH9TB7emCZEvGWQLQbGwdygzofwRqJCD5hvxqt0hw6isEEQ5YQNv9smz6FC/XNEzrsdyJa/psPsd7+/9VfPbNZBmwAm5691Y1cDc76WvF2Dz3TAY3Trioc19PLq1ek88L3xzP9uaIP5j4RWcqSACy/qjHdj694dzz05oImCqQhzd8M6DLxKizZCikDI44VBWddG+1/+2vY+O4QyFyNYH5ce1ObmVecPPIHrWW+0DwtFVJfv+CLVwWtsxysB4ngef77N+I3b/FfY8fsD7pQTa4YZegHtuhbo2cq/+BrfctZZ9w1l+eNMFnLGkZaZHd1gQj2l8883ncu15S/n6Izv4+sPbvR+74y6CN//c8qpdcAts+RPcfjn84rqKPy4KhWLWsQJ4Qei904COqd+090MvpWQ4XZiizYSMjF2Pwa5Hfds+hHUb09cn6/meMf4aMs43dSUPqQJb1X2vItAc0RI2lMOL5kvG/qvhTYBDyaeK7IPme++J7QOeQEvUH/D6N/eMsWrPsHUwh3ZGGpgTngRwxhHyQI3nSzy+fWKFVb0ctIltMkymoFMyzEDOXtk2ymr528u63qUJCDSjBIXxAzq3rgdtIgsfZK68Uxyksgh0FPCBr/+5bstDGyUo+sbyPLljcEpyBx12DqRZvnPQfdaWCTTbeyuqCOyo4bk1T/K2tzUsXiuIVr8IPpQeakqgHW48/l/Q/zzm677Lx+7r5NnOJN9+y7lccvK8mR7ZYUVME3z1hnN48wXH8d2/7uKz924O3khzT4KrvwKf3Ax//1nofAZuuxQe+azKUVMoZj/fBC4WQnxGCHGKEOIG4KPA9w91xQcS1ri9f5yndg0xPJaGng0TM34maiBFeYBKOQ6+rEE5XlE708rBGZiqIpgHU8VREi+mDrzMviPQQksVdIP7N/WwayAd4UELvl7fMcr2vvFDLopxKDZrlDgqM74dg1JzctKDnw+M53l8+0D1cMcq16Pr1TpII3R73zipCZZm97RTxLaKGUhXF3qPbu3nqZ1DvnL9FnOH1jJvcA0AK3c7OXThnKr9h8W67H0SdjxU1dNXHuJ4IB60gxNolSYfytd/cI2qkZKoGkwbOkdJZouRXqsoupM5xvMHVq7fqcKYt7dRKcRRyPJj4FwP1XLQDvTynqxilkqgHU7s/Rus+A68+J38z54TeHhLP//+2jO46qzFMz2yw5KYJvjydWdx65XL+PnqDj7yq/VuKVqX+jlwxT/CR9bB2W+2Qh+/f5Hqo6ZQzGKklM9gVXJ8E/A88F/AvwP/eyjr3TOY5qHNfRUM8/Jf5XTeMpD1gW1WnsLInv1v5FAE2hQhnfyhCgItW9QZOgRPoTgIpdI4vpdFA8vRsuW5IVVn6yuUUi8Z1nc294xFVHEMrqJ0CLPiMmAcWmW+nQbSB0JUSlfZbru5NtHjHUgVSOVKjGSqbH8C15kedTyGdkGys+r3Kh3G6EkQx5CO+GjHg5Z9tB9SEYZ/fb6fuoItzMqOk32t6FYomzT1/XtEckl7nBO/RnrH7PVP5DsH4OHyr88ZT8Xx+5atVHZ+P1tzBU22qPNc1xhSSmrjlvdxohNba9tH+Ou2AWs87Suic78qDN15ApcLtKAHbX/eZ4+JCufg55OVi6cE2uFCdgTuvhXmnczv5n+QH/xtDzdefALvuuzEmR7ZYY0Qgn+56nT+7bWn85fn+rj5jmciH+I0LYA3ft/KUUvUWz3UVn53dpb3UigUSCn/LKU8R0pZJ6U8TUr5HXkIcTb5kuGG8lTsJWRt2P3T/Z12jKoKxq5pSnLF6stU2070xwe+q6Yp6U5mMSoWdIhe5yNb+lmxaz+GVNfayga7b5939I9HLpIf6qB310b3daJkhxgWyyMaJuRBCx1n//HaPZiu+Bl4BptAWPl6/eF0x8qYA9tJFJPuejd0Jlm5e8g7/+HlTRl5LqOMxvL3gjlV4Y+d6zOq+KO3isrjStuFPco8aIYOvRvYsPJBdg+mWb8vGSngo/Zh18A4D23uc9ftLWsPJ7rxXZUdmPiirlld4R58vmuU+zf1HPJ2/N7P/lSerJ07OBHdXxgfmvC97RcqXh+56iGOfWN5ntg+MOHQbGcbQpru9bR6zzB7htKMF3Rq49b9lisZPL59gO5krtKqgpRyMN4L+1bvd9FwAZSyPZQm7cMZhKwc5lntmFb+qFJFzMpjPRCUQDsckBLu+xhkBtl40df4l/v38NLTFvC5150R6VJWHDjvuWIZ33zzOTzTPsINt60Kldz1cdxF8N7H4YVXw8P/Br+9ySrBq1AojmiyPgN64s9dp2BH9TJt2/rGeXhLnzXLPNEQpiloVN01mqN3LE9P0prRD4Q4wqFNSI22Q+fTgCVwV+8ZdsvCC99ot/Z6uV3Pd4+5M+/tzz5G97anyeedqoyVj+nB5KDlqszwlxnO7mvTMiQHygVa12g20lsh+zaxcGCFuxpnQrBSKe/7NvWwqWusfD0TORUVPGgxPQOm7r5dXunSv47ocflL9Zd50HKj7vvPd4/RN5ajfThYSAGiDeUB25uYDQk0Z9FqTcR39o1x30afgBrtgOd+j1kqFxuRq5Gh6zx0n49lC/bHks09Y1W9QhM1+P1CdCLn9LlnV9HhHPut91Ut3hI8Vo4HzXtn71DG2wd72UzBuh4zxXCxEN+6siPWscUvAj2BNp73hLtwvWoGqVypam5pUOhHn4NqiAoer67hcYbSBdeDFnWcq4ekhiZoTEn/uj8he56NXN7vQYuXDj4lRgm0w4FnfwFb/8TwSz7FjQ8UOWl+I99723nEY+r0TSb/cN6x/PRdF9E7luMN31/Bmr0j0QvWtVhl+V/5RWv29K7XWQ8shUJxxOI3xsqML2N/IUHVM8OSOStkMpktHUCIY/S6RHUtCMO7oWtd5EduGFS4SEYVuXcw5dq7RrP0p/L0jTmGs8QwJXrIlbN7MM3mHkucODksjlg0q+QFVS9eUX6ABscL5dUeq6zPeS306Im84XSBdR2jrsc19GUXU8rqnhb7wyhxM7ECHSEPmv3qmL4nSOx7yt32SKZY2bNRyevrG7hfiEopkaUMppSYIlFhPM6y5eud07eKOSMbygzmshA9o2TlRQ55Jc639owExtW9+3mG0wXMiP5lUcdPlL0XFgfW6+HRUXb3JYOl+X3opsnz3cnIz8L4nyvO2MeypbJ7wc+oc8/pBRjYXHE5v+Z3D5v9R75ksKkryeo94es+2tsaYPdfoeuZwPpExPJP7hh0xf9E8hQDzx1XJO/fznWbtld4zO4byQQ+kFKW7WC0oI4e886BNJ29fd6zz3/NJXOBEJJVGFMAACAASURBVPiFA8v3O/5KKAt/tjO8Gx74NKXjL+e6jeeTiGnccfOFtNSFH3yKyeDyU+fzxw9dRlt9grf/eDW/eaZC3xMh4LKPwlv+zwpvufMq1TdNoTiC8XtYHENCmEUY74ct98LgNvtTCaPtjIwMW0ZWNSN825+h73nmME68lCFdKPGXTV0TayS8Xw9ahY32PAujeyM/8gyc4Hc7IgSCw/KdB96cOG7X43aPozR5bGs/f35uIs9QZ2ye0CrL4aqqz8oF2v7yY/yHw8wOu4JcONUJ48Heo45HKSpsMeAjkNHvR30eZkJpcFVy0LRcUMxUCrEMTBj4lq9Uqe5PG3t4dncvUkpMrbqdEuVBqykM05DtLq+kGf6OE9rqu5Y10/GQ2MJ2NMfe4QyG6Z9c8dYZL3rCLV5KU1twwnSrhwHW7n6I+UNrKuZz9Y7lqxa68H+rc8QTxqaUFHSDJ3YMsMkv7qur+Iof+b2yrgfSDApdJ/fSHVWl9VV435ASzchTn6t+77oiqtq6ojxoE8CdMAlty0G4RUJ8OWih0vjV76fgh87zwr0WR/ZAMUNRN1nbHpysF4eQK6wE2mzGKMEf3oOMJfhI7lZ6UkV+eNMFqtfZFLNsQRP3fPAyLl42j0//4Tm+cN+WyrNZL3gNvOP3MNYFd74GRqINH4VCcXgTnOm2/l04sBIzXJwgO8K+TU+y/ZmHGR3oYmn3XxBFOww6bJiUcjC4jZae5Szqf4KSIRHSrB5y5uJXDZWNgJW7h6p6h/w4IUJDmSL5kuG+HptgtcKJ5sbEY9Z6HdtQSLNiiGE44V64FR+t93uTOVbuHgr0TQvnLwVXEPS86UblbTu4OWemTur5h2jqX2utQrePS6ym6ve9FZmVTf+ooh8TGFMUMT1Ha3KLb6XRXhH/OuKVSs/5B+YTOjHf8mGP6/j4mLWtsvC04Osyo9j0wuQClfakl4fn/hQ7x94njv0GuOWxsUxcafgEmu+oLvJ5NxqyXd44KoQ4Osa2JgQ1xdHKRv3EtUWgMJmUkMpZ126gcXdZ8mD5enTDZChdoG8sz70busmXjJAHzRFkZmCIXt5dSKhNEMOQzB98mrbk81WXC1yvE/DK+lsdDPa0s7mjep88e9Hy9QCavU8tqZ3MGdmIaZplRWWinl3uEQk9X931+8/D0A4eeH4CE0wjeyzP7wQiJZRAm8088SXoWc9d8z/Jg50xvnbDOZx/wpyZHtVRQWtDgjtvvpCbLz2RO1bs5Za71jKWrTAjdtJL4aY/QX7MCncc657ewSoUiinHXxjECdeJ6xnW7xtlOF1ASrtoglliYLxAXM9QU7RCoGLpiMICvrBI5/d+NFtkwiXn/d6PQtHN5/I+t/4ZHC9MUPAF7dFMQfe9npjRNtGeP44YMGR5TFJzaifzhp5xX2tCWPvqWku22LAH5+TM5EvW/o/nS9GhhQ6uV8lafvnOoUDeW7WvOAZ6LDvo5nHZnwBWTk+mmjhEBg6l3yaMyq2qnsdU+bM5oxtpSu91c8EqUTu8laVdfwYpPcFlhMYfcB96n/mvlcB5l5L6bC8m1ftOQYTILGU9L0hEeF5gW3ZlRWJ+gWaNT+5+nLVrnkKKmD3s8ntNhq7VoKej0rF1ciYrjN/GchAf+LkzpVd4pbE27v9GxXVRyoJRYv2+JCt2DbHV7s+XypUC11SZQAsV1nAQFcR85BhME900iRlOFcrKwwxeI9HPOP8yz3Z6nqiO9Q8zvvWxyit39ZLwv3TR7O0lSikasl2Yer6siXvUo8vbn7BAs9frV2j2ws1jO4jpVQqhOAWF9P0XYVECbbbS/hQs/wZbjnkDn991Cp94xWm8/pwlMz2qo4p4TOPzrz+TL117Fqt2D3H1d5azobNCXPmx58ONf7RK7P7iWpWTplAcYRR1kxo77zdsmA1nCuwdzrCtL0Um5/zwCkzN8qz4fAHel3w/0M7qBscLEw+J8S336JZeNnVZz6aKlcwmgN9Yk0Ai5vUMq5YT41DZo1AppMneru/zltQO6vJeT6uYJgKzzY5h6RReccLXNEwojJMpTLQMubWecNXeeGmclrFt7uvGdDumXgx9U3JM3xMIR0UUUuhj1jlYuXu4kgMGpElXMhtYj7NsZGU53/dc9jwJPRuqGsOaae+TI6jcYyYD56JudIe7JSGwUiq2/BFp9/qcP7g62KA63KjXxm9YC6mjyRLI/Yd3ledyFt2iErppcu+GbtZ1jAaOjSml5YXotz02cc976RaByAyhDW1F2vlLpq/xugwtS+iTom5af3attaoI+gi3g6gktDQhqraO8OYawuvDraAaE+XGfxRjY6Nktz7sy4fyPgvmBlr/OiGN7nEITcJUa+S9tTcVDOEzLRHoCOFqz5zALlTYH3+Yba/jEbevuZiRrZjHVu6IDYnvkMCSpgmNC9zXMT1Tfi6ldNdc9plzH/gPldBIFMdoGd/J3JFnQRpoRpFSoiV6tBMofqIE2mwkNwp330q66Xiub38Dbzx3CR/9+1NmelRHLW+96Hh+e+slANxw+0rueGpv9IN5ybnw1l9ZYY6/vEE1tFYojiCKhkldwp6RD93+uiHd3kEl32y9Y6BG2iOO0ajFAgbG/rwOHtIei8SZ4S35RFQlmy6VL1XsWRauTul5jiTtw3bFt/F+KzSTcsOlLEQpuc96Dtrv50o6Y9mSV7Cgwgy1H8uD5gsvdcSa00TWNmjrh56DHQ8hM9UbFnvji97mgoGVNI/vRphF6nIDtCU3E69QiEH6hWOv1QKgrJdmaJv+oipr20fdannOsZNS0pN0jq+VG7W0+wE6dm9jJFNEpgfY/Nw6esYmUHnR8TIhYbyPpj0PoJneuXfzgqRJfnzYrUYpC1Z+Tm0hFBprRlccDAo051wb+/eghU+BlG6TXyePr2s06x6buKZZ2+pe716D/vBSNwfNWZ0t0KRujdsJlxVmqWxszr26qTuJbhpW1dFQDy4hDZDSvVqjIotN6VQu9I5DpRzHqEIoznFNjHfC1vttj2Zl6bNzIM2Wdi9nz72OiBZoumlaIaPO5EhZ0Zygh9q/gh3948FiMqaBYUqk0wjduZ7M8pDoiYQ4rgoULJHWd+zJBiFNSsXo3FwpYc7IRkRyb3BX3GEGj79133oLHdP3RESlVlnmtXffsK8zv5g1fIJOSINF/U+xuPcRV7y61+kEw8BBCbTZh5Twp48gx/t4Z/J9nHniYv7n+rNVOf0Z5rzj5/Dnj17Olact4Av3b+EDv1gf3czxpCvg+jugZ71Vgn+/1d0UCsXhQFE3aTDHqSmMlHk7dFO6P6amfc9LhC8nxomr8v2btAsQaYmA0JmwQLO/s61v3PU4BYoChIw6x1jb0T9uVQWMMBQCE8JjnT7PlSSZs4Vd+3LYtyqwTofAccmOQOca6FyNY7hs7knxxI4Bd2wl14NWWaCVe9CCfzjHK5a3GwQX9tNnKRTiGEazw+SElAjbOPQ8MCGPR4RgqWp/hT4c9eX2Oadu50CaZ9pH6BkaQex80C2+0NG+i+1P3Y1hSnIlg9pM5X5cbjXCki3QjBLsW40wCjSlO7zheKU3SG56gJJ97PzGeUBcVCi2EVXcYftg3ndeIw5KqodEts99qRkFKHihprqn3t1jE48JW8SEJgKcdWDfb+54bBFve6stra+zpOdhFvY9ERiO/xo0K+UHSRNXOGBd7/mSERjP5p6UnQ/pvRdu2+NdK+XHxZTQNL6Hhv51lsDW8xUvKv/3HUHrvw4NKYnpORoyncSzXg5XyfD2QQjrHA9ni5a43E+BlOBgdUyJJ0Js5g+uidiv8mukGkKa1mSH71yY2x+KGIMBRp6GbBeJ3vVQzJQ9n2vs4R1v128wDaMslLc7GRZ/0qfPgje3Mynkf14amRESJa/lUlxP2/vheLFD98IEhJoSaLONNT+ErffxLd7GUOuZ/ODGC9xO7IqZpa2hhh/ddAGfufp0Ht3azzXfXe6GFQU4/Rq45luw+zF48F+mf6AKhWLSKeomi/qXs2BwVVkITdEwPY9OyTcpU9YY1f535yMwbJcHr2kMCTTPUEznizy6eh3pyIpw1neyRZ2E3WtHSl8ehgx61Mo8OxECxT8PGO9eQ9PoZtvrEDT+HIvZ7X8lJUgjYKiP53JW+GBhvMwYccOtXIFWWZRqgpDnxjGOzMB3nbL7ur6/giY+8dK+gtbk5mixKg2EbfAbaNaxDR8zv1Ef8mBEbzr4/fpsN5ot/hxx7YR5yZEOZCHtCSqhkSilXKMxoVeO0HA8R/hzYdxQMc/zJqWkLh5z98vwGflOHs3zPWOM5+2y735PZoQ4AJ/XWItDKKQyQMdKGvuepj5r5Wwv6l8O3esiGyo7fzshxqGrKfi3ND3T2t0vr7iHtI+3Vhau6fM2GdHXY0JPo5lF934t6CYPbe5jc48lLE1TuveZv+y843n3053Msa0v2EPV8sFIWse24p9qYSy6wbs/x0yGrj+JxDQlLantzBndREvfavdaKxmmT6AJNvek2DMwbk86h+51H3W5AYT/uJm61Ug9JNBqSuU5oH5BOJECGUKaAZEOYJaK0BcsRlLc9Th1u/7ivbH9AYQ9WRPTczDejzQNFjTV0lwXt/fNcD1zDl2joYkd6ZviCj+/okIcM8PMGd0UuR/2t4LrmkAouxJos4nu9ciHPsOq+EXcKV/LHTdfyNzGCVaHUkwLQgje+9Jl/ObWSzAMyXW3reR/n9hVnhx//jvh0o/AMz+CZ348M4NVKBSThiElMbs8fFljXjwxZPoKR1RsjFpIQU0T1LeBNAIBfn6xMrJ9JW1D6+hr3wKpkMfE9wPvfsfvWCDohSgbs20kDaTykYUtTBM0PWvn04Q+rGuxlzFpTHcwf2g1S7sfDIQJrdg1xI7+cdfr4A3MNzNtG3bVBJqlD72QJG10NyQ7XZHoHmNnP0v7iVqwPTAyNwrjvTSl24kZeeYNrqEh4xnCQpqucWUibAMzZKj5ZuHNiHi3shQ0337G9BxzRzYwd8TqSec08PVJHftLofwZ+1/NZ2Au7HuSpvG9gX20BhUK7ZK+GX37tSa8Y+i0K8gUiswfetpdbnv/ONv7xwNCOVBl0cmnlKZbDdF0z60emY/leObmjmyw96fA2o4RN/zW78VyQxyd6p+hoiQOwt4/LzTXWscmu+CEJjwhE/YwB67BcKEUHwsGny7LG+wYtrwvBd1/rrz1N+38Y1BYSFjbPsLuwaDIln496/cURzREHs4U2BESeNaYfLthSleUSSAurP3Kl4xA2XvrXFiCzi0SgiRT0MmW7GNRyjFv+BnmD/m8Y6aBlGCKeGiPywk8fszye1RKSU1hBM3wwnJdUedMvpimr52JRS7phTSbta3WHwXruMwdXofc+zdMvUQ8JryQRMPwQswrIU3fuQg1MDfLn+uVJmZEKA+0zJNWBSXQZgu5JPJ3NzMs5vDh7Hu5/cYLOHlB00yPSlGB80+Yw58/egWvPGMRX3lwO2/94eqyMAZe8R9w6qvgL5+yErsVCsXhgVHeMNqUnoepWnPmQD4OoRw0vwXWdjyImGUM+n6r/cZsTdoydmv7N0LHyuCGIrxuYYlQMk1qCqM0ZDqR2dCsdtcaGO9j1Z5hHts2EF6lFUqGILL6un1sjNF9tCWfp7YwYn/f9iyMdHqGXMDqhDmjmxCZQeqzPa4RX+7NCO2mT/yI0b3Q+bSvt5HtOXDC2SKM63Ao2bp9o+wc8IxbIXXqCoOBGXAhDW/d0jNbgyv2l2+30IzK1dmcyoGNNXEvNLPS8rbj0vHiuV6xkPjQjCIJPU3r2BbfV4PjzBZ01naMkC0arrCTtkTRNOFeP/12tc/VuwaJGUGPQq5kUPCJX+cyj+kZWvY+QHNqJ0u7H6AlZRUekXYPNBEW6AB2aKD/dZjEWLsXYuqEOGrVPWiasEW1e4zCd4So6LnQAl7aytdQXE+X56DaAwzmmtniwixafw7tYP/4i8Y4b0WPd+9QhmxEbpsrQKXdo8ypbCmhKW5VRNU717k58v7mzuFrfGtfii22d9AZh1OZFoBCynpOaJ4HrWKFyoBy9EcZWO+XsikWDK5i/qA1MSCkSf9Ynj1Daa/1ghn8TpjS3NPsPyxbLGYWrBw5JHFNc3MQTdOoWPCG8T5ID+CX8IF9Mg1auy17zi/yKwo0d/JMhTgenkiJ/NOHMZNdvC/7QT597SVcevL8mR6VYj/Maazh+297Md940zls7U1x1beX89tnOr2bWYvBdT+B+ada+WjDu2d2wAqFYmJsuTdYvQ7rR9oxFCoV2QBfKA9meQ6aP+RLCBAamKWAwdec8gw5z3ayF+jfbPXQcdblDq7ciJNSoqcGWDC4kjmjm4jveTQgdMxUr1UtGL9XwbcfWALNX/BgLFfimY4k6bxlyJuh0CZXnGbCfde8NTdku6jdt5y5I88yp9EqkV7NgyYB+jaF3vD94YhTO2xrOFXeVDtoC1kvUnnPQNMiihokNHjxcS3EhLBDHGW5J9Fn5I1li9QURlnc+yhaKjokzQnPtESvszKnjYBdcCEc2ue+dEqIO4a/te2EHU5mar5m2SHjL2nnS6fypcC+mtKpOBgOf40+H4MpT7Q544wZVshfSyooQByvSkCg5ZNWOxppBrxNMaO8+EP94CZa7Wqabohj3DGwoz1oSBPNNHyd8hwR76/ZHy14Ap7FiBBH/xGtJEKCHjQ4Y3ELMaNkfTfmNe2u1EReSmtfTREvD4uugive8YSdYUoMU1oeRWF1AquPm9QUR5Aje4j1rgcs75TXnz7iGsdqXUHIo7lnMM3uPbstgeY+ByrtWciDFiHQikXrmVqj2/evNEkVdEYzJfe5qzvn8fk/WM9Cgn0SZbwOYgm3gbxEo2Sf77gm3AfqvmGrcFGUqCru+hv71j2AaRdTcfbL/beU83nUKuwfkNA05jTUWNVlHXqejXCRVkYJtNnAmh8itt7Hl0pv4eIrr+JNFxw30yNSTBAhBNe++Fge+PgVnLmkhU/9YRM33bHGa5pa12JVdhQCfvUW68dJoVDMfnJefqnbM6hKCWqHXru6nj9EToaMSGsBzfrPKAWKe2g+Q9F2GHgG6cBW69/x/qAHDV+um68om1H0GdSA9IVKmWa5gRIQB/aQ/GIirxuYWoJNHUOsbR8hXHnfzR8KGR+9Y1nWdpS3HtFs95xjHPtzdTSjYBncqS7IDrsGf9g74vq2TJ2u0SyFoiee46UMS7ofRDoFWZwDE8Kfl+Uwt15jTn0cTVj92hw/XQBfvtu+kSxx3fLKadnoSpJucQGnt5u9O/OH1rC0+8FAzqBjlbpeWKcioTP/53jQbMFlal7frHBZcc3nBtX8JeelLCsJbxnpFbxM9lGwmiLboXNCYETYmo5XpWxdOx8BaQb6CtblB8u+b+KFcTrXqRNiHLy+fB40AGm4yzue3caaGPFSmsW57VCMzt3z51bJCO+KXxRWajfu3ychJfOaar2cNS0R+R0/AQ0Tb7DfLN9WWeis8M69mwtpGohUD5qpY2h1mKakRpMkBJR0E5lP0prc4oofgbSjBJyJD2/9Y7lSYBLINCUj2SJ9yXFfqGSEePaP0b8fEeGFhaL1Xk1cs9dnX/dId+In8LyxK2w6BfSk0DDq51t98fQiSEnCyLlh3vGYcAVPJmfd71Gh6p2jWQbCvSPd8ERJwLPm+154vzVNUJfQgpNn/kl65UE7DOheh/ngv/KocR69Z7ybf3rVC2Z6RIqD4Ng5DfzqvRfzxTecybqOUV71zb/xs1Xt1k07dxm86edW75bf31I1vl2hUMw+TGkZcF44UNDojBZu0jX63KWlLBNoY9m825w2vF6nEWqZudy+3PcD74Wohe1Gv8FvSkkp44kkE89AcULJokMcReADM1aDkAbdyRxmOD/KFiDhJsDtg0GvlvOpoxuc8LKamGeSLO59lKXdD1DfZzWtHm89LfBdp7qh60EzDfJFI2Box4wsQhrIXn/yvje2sVzRXs4TdcfNaWBJaz1nL20BaSKEwLQ9aGVGle7z/EiJ6+WqkP/vep004YkoIagtWMam/3wJRKQ/wvWgyRKt9QlPqMrykEFX2Pguz5hdZt8VeiL4XUNK1/MUsQMA/G3jdor71oAt5qKMcs+rYpbnoEkzUOAiUSpvFC6lfY6l4S6bMPPOKj384a+YaNIoK4gjMGlKt1M/ttuqTuqjrd7K89ekjinilBLNYOrsG/FNKpxwGWZda9QmA5QCH5howpcrGPNEVKUQaWma1A5uQpM6uuaUZfeXC4k22f1FOpxzro3to67naWJGDjNWi0QSkwY1MXucepGm9F5fLqM9YeDXIvj+DhTEsf4V0nQFvabZHkB7mZpCsEl6RQ+aM/lTtN5LxJxJm/JryvAfX1vkOftbqF1grUmLYUqDulwfEkl/yrreE5rmiVF7Qih83Uop3Sf52r2+KAB/eKI03X2pFuIo8Lx74SIqwXVWRgm0mSQ9QP4Xb6PHbOO3S/+Vb7z53MBMl+LwQtMEN15yIg9/4qWcf8IcPnvvZq6/faUVw33SFXD112DXo/DIZ2d6qAqF4gAwpUQzS+5srRcOZb3WKvySajKYQxMsmCFAiEBjZRGR+wKVJlttg0ponqfOZzBIGZwhllJS0v0eMqsogGYUWNLzMAyGQtSkpDbTQ21xJJDTZGo1lkEtTcyQC81pGh0WbqVwHyJ7fY7Z4giEmrjGKQuaaKnzvA2maVWqLCVafXvt4Y7NaVgt/cfTLP+O72A6ldv8HrRFLXUsaasnIayqhZoQduhheW6XsIsReN6BYBhiGKdIiBDe2PyluUuG/5vRBc+dTb1wQQNnLmn1ifPK3/SHgfkxbU0pQsZ3JQ9aSS+xpSdFS89yGjOdCGlV/StFKBavAIx1zQcmMULhZVHbc45pQ7YHKSFRHKNt38NW4/DAgv7zbd1DTtEOd588i9rXG87CuXeFqVs5oQikYbgeFAkMMIe1+zwRWcn3Ea46qGkCzSxaHjefB+257uhImrrkLmrH9gCQr1toD96X51jBC2dGvD+a1d3r0NRqrHBWdOaMbbHzstxvuy45fyEcw79dZOA46z7BYh1aSzrafwIEchgTxSSxjNVSQSCCRULs85wv2B40d5Km/JpwHzexhHceA0JSWjeXaboTEUl7EqauJuY+v51+ec510eD0t0wPuZNGzeO78V82AJlCiUy+REE3GGs9w/ugtqUsxNF/y5ni4KSWEmgzhVFi/OdvR2ZH+Erbv/P1m1+myukfIRw7p4Gf3XIRX7/hHNqHs7zue0/xxfu3kD7rRrjoVlj9fVj/s5kepkKhiMTvRvJyyKyy63DivEavIa7zgx/yoKWaT7XeD83yOjOwgJuD5oYxiri7Xj3eFBhJWZVY8ELkRCxQLVIAMT1L/a77ifU/F/hKPuRRM01J3Mn5SHUHhGDvWJ5cyaQh3RE4JnHba9SaaS/rF2XY0QHSt52iYVol2gNjdw6B0zjYK2TQ1lDjzqKDVbhiS2+KgowWGW6BFHvbcT2DZuS55OR51Lo/qdEmtSNg63JWj6hlC5p9YzRhaCdCWIbZwHihfD1OGXcpbX+X40ELhoPtGUxb5cjt92OhvK8Wu/x3yfAqx23pS5HMBJueC4Q7BE3qLGiuZX6TXYwjaKla264SRmVKy6B1wgIBxlpPZyCVLxNMhp3ftnsgHSiuopkGmi8ELXBo3AIwRTSz6Iau2QfIriDpTHiUG+OBkEIpiRk5YppGfa7PFW9juSJ9o17IohBmBbFn4gQmipBAizn3sNSRQljn0CcgDGn1Ngt4sioVqfB7QO3908ySdc60OA2ZThb2/y3yu9YKsq6hX6iZa2/Ma4psVBRo8bL3+lM59xjGjCxSQuPABmrNbCDXUZMExL3j/S4aQdHuX2anXT1SOC0NHGEEbjSA/zwsHFhBS79XFTTKg5bM5KmJacxrqmFOQ03keXRDwRMNoBcgO1I+GSI0pCz32mrCigYQCHciTHc8s/a1KXc/7jpJGjMdZf0rl+8c5Mkd/RimxNTivvRQEbombC93VQ9a5XvTHfN+l1BMCcl7/onm/jV8re5DfPa9b6W5bv/xyYrDByEE151/LH/9xyt50wXH8ZOn9vKKrz/J3Qs+iFz2Mrj/k9C+YqaHqVAowvh/N02ddEEnVzQsgWY5vWhMt9Oa3IIebwRAJLwCDeNNy6wwKay+SUasPhg3FApxdIxUI17vCS2hBRvyRv2YO94h4c5du0OP6xn0UoG+IS/MSMpgjowpLfHjbPOpPUme7QyGJUlkWRVHZ2Y8gV5WWj7WvwnG+wPCbVNXklwhWFbb8zc5RmLQc1hBitn7ET4W1hiGxnOYWFX2Fvc+xsLGGk6eX+/uu7ttn+FnmFYOVtwuUjG30VdowzYiY5ogU9ApGeVGn/MyW7DbADgDL6Zd0dkxkuW57jH2DGVc49fK+/LGsaTNyjcazwVD33rGPC9EY6YDTVh7K4WG0LMgJYtb7DEHDFoZGJ+RaLTf947snqE0RcMMVHE0YrX0pfKEi4bo9velNAKha0LqxCqIXydncP7QGuaObCjLwzRM6YZe+o/FsvlNJFoWBfwnht1rKyagtjCMKSUjmSI7B9Ls7vdyRR0PWhhp54NK8Jp32wSKTIhYmbFtmpJMUQ/k9VUSvsGcJmkLcd32VGrMGd1keUzL+ulJEsWkXfHRFtfuxny5cbZAC3vMosIBhTTc3MBc/RJ3+iCuaaQLOj12rqzmC2c1fWG847lS8H3fPV1yXUsmpgkxYU9NSMmeIadBs7X8svkR1cj9EztSIkt5Cv07mNNYQ2t9DS89bUHkMyBXNBhM5yk1LbbeCE0qSfs4S9PwXcOCFx7T4i6zuLUOzSwhpXTFcELzQsmF73rwesp5uBNCIu7btnDDcBc0ec8QNyI+SqCpIiGzk9EVd9D2/E/5pfY6br71U8z3nVDFkUVbQw1fYgogGQAAIABJREFUuvYs7v7gpSxoruWTf9jMm4dvJdt4LPz2Rhhtn+khKhSKAD7jzNB5bGs/T+4YtDwYwvKWNaf30JTeC1KSqz8Go67NWl6rJdV2Ovn6RZQSLZgiQabh2KCx5fyq2x40L69Nlgk0r+hG5XHW1dSU9VvTjCIjmaKVd+LM4spgGfBUrsRwpugWmZBaLCQEBaaIE5OlgAE41nqGva81yJAHTUvuI7v98TLPWn0u2MPNn/8EwZ5e7rEp21snfNCjFG9CM0vMSW628nb8IWYdT6EVndy3cqPL+VvWesZbYLtOQ2O8ggUnzrPy07wxAW3Hs2Ng3DpG9qr19BBtyc0In7E/kMqj28a2EDBveK27npqYRm08xlgmRzXDzbBFe6bxeMtQ7FrrazcQzMVyx2fn0EF51U2wxaK7vLW3Ye+F49EV0gyErglpBPLbHBpq4oHS6xC6tOzcJc12H/tFlRCQXvDigAdNSkuEON6N4UyxTAhYmGXi0nrbd2+FrrWYbwbCusaCAk2X2NUQ/ddN+SYA9EDum0TTfOPzXVvhthJ1+QEWDqygJt3jnn3d8Rj7qyfaeWlS7H9CX0gTwzDpXfwKxptPdofg9JJzvH0CbyJISi88OVfyjVEC0gzkiDrfNe28LSGCgbULGhNceOJczlzSQm08FuqD5t9/ib77SWqKo976TZN4xHWVKxl0DGfp05usPpKFdPndIuzCHPZ+xBM1NNX6CugI63rrS+UZtyu5OsVnrKI59qhEzGvc7h5Tr4COGauhMO90OPYCEN414Rf8zp9FKdg3kvFyjaFcpEegBNo0k9r6OE2P/DOr5Fmcd8t3OG5uw0wPSTENvPj4Odz7ocv49lvOpbtQy2uHPkwmVyD3sze5TRUVCsVswPqlLegGewe9XBErxFEEbHjLJBFoToU99wONgUVX0Lv0VRRr53rGpp6nZOgMpvPkDXsNjl6DkEDzthNZGc0pABGLUxsL5j5pskTOFmMvXGx580wkRcNE2t49p9qkE7LorwQHUKxpo1TTaud0edsvtC5jUXMd0iz3oOWKOlv6UuwdDD7TmtLt5eMHtwR1XE/TkPEqLYZts/5Ff+dZO45XoG4RRqyeRGmcZZqV3+Kvhkl6gJrRbf6vWMch5G3QW473XvgFTMkWIra368wlrZx+THMg/BJANh9DuukkqyGzvT9F3SRRtDw7Tu+uwXSBrXbuUSzslhRQl9BI56OLRyxstiZxS/POoPCCN5BtWGqFOyY7iNmVgQWme004YmJrT8r+xL4+IkLhYpov3FJomFqCuB6sdDjW6hUv84spzSwRj6jLPr+phtOXzA29G/QuGdIT6P7CLpoAoWnl50wa7mRDMusP/fTloBGsQJltWEop0eJ50GRwGIZWG8j7N2M1ZSPN2jmi6aYTvNFH3o4SXS9551ZKt4WBhIBBLmRIJNqiV080u9t251Ls7+m1rV4lz/B5jDD2hTQomdLytglBuulEhCi/9mKYrmfXv19CSlrq4tQnbHEljYB3ydmGlFZ1z1LtHEq1c12RddI8K5dT06zsQ7f9iIBwGweZs65hN5fX1ImFjtFpi5o5e6k1EZYtSsaMGihmIjxoMft8G5anOXSshBZDM4p0J3PsxKqYHqu3JiD8zwYjVhcco3vxOMtoFOeeBnNOtMMq7eex7/g6f43lTQbGC+wZtLzWumlGX0QhlECbRsY6tyB++w465CLib/05ZxwbfoApjmQ0TfCGc5fy2D9eyVuvehn/yMdJjOzk2W/dwMaOcN8ghUIxI9g/nPtGsmzv8YdzmfZMsbeoZpaQCLSYZdgnQnnEDTVxtwcRAHqBvmSOjuEsuwYylnfD+VD6Z/81K2cpFLoYHKc9k4sgLrwqjnXZnkDZcqc4Q7ZgMJwu0lQXZ0FTrbvuJml5mfyhkgCGVmPlxUndNdgLtfPQNM0yoE2jTKDl7RDKTEhoOMZOsu1MJN7x8Bf0mDP6HNjePP8xlmiYNU045ooz7mLtXGu2HEjYRqEXYuaUqI/wxEmJvxpezBeeit/rU0gFxmIJMxk0UiUUaue7eYDO+UvmSpgxa73+aoVO1ciwoasBtXGN8WzBLVzi99i01idoqUsQT9RgCg0jVueOSzO9CpQxI0+8lME5jyXTRDela9ibolygaUJw1hJLtEuhkW1YSn2u392+Hm9CajU4njW/90dIwxVoqXnnee8DS+Y0BsLbgoU9LA+aY8xqIYNdiKCXygovNd3lDV8+ovCX1Y+LwPiyDUuRCM9gD+27EzbpjtHxUEmvqMneYev+yDSdyOics63lIoxrKa1S8Y5AEfgEmj+0GRAhD7NVRVQwsuRKzwvrtMGw9yfddKJ7HovSu34LtfPcUOHmOp+nSJroJu49AoKYprnja6lLUJ+IWWGK/lA+37411SasKqb2+MPhxUJ6HjTn/zVxjZa6OC213hgtUSY9Yeb3oA3vcouOxJyxmjqNwrqul7bVs7i1jubaODVxDYFg93CWDf0lirlUIAdNIq39NU07J1QDW3RzzNlw4hUQqyFmP2f0eCP5ukXE49YyZm0bulZPsaYt6OWUzrPFmwCQQrgVIkHYx0EECkY597lhP2+KhsmmriTPdY1R4akeQAm0aSI52Evmzmspmhojb/gFF77wpJkekmKGqEvEeN9LT+arn/44K0/9J87LrWLljz7GO+9Yw8pdQxUTkBUKxfQhTevHOF5MWTP89gyy3+jXzGIgj2x+Ux2NNZ6RdPaxrSxua/SMzcI40vbMjBf0wMzr+Se00VrjJZWbSNdWch1wpsnGrqTVOFYaVr6MXZAdrJ/8lsH11OW9PlzO2PrtqnQtdXEa7OoZMU0w37Qmhxoz+wLeiHTzMqSWICYNSzwB482nWB6OWBxp6vQlLe/byNxzrd0reWFTDsVEq1vG3tBq6Vv8copNx1pjC/Ufi5UimkwLjabaOIEaIQteYI3P3jfH6MzXLgAg27DEHoe9jnARBIEVmhqrDwo0vyci0LNSEo9p0LHSfSfZdiaFJRewsTeHFDHmN9awoCFGbdwKjRKmjjAKiNE9LOl+kNr8IE0l61iHwwIdQRIMLQyW3BcCkgWTVbuHMbVar8CK7oUczhtey6L+JwLrHsqU3BDH05fOIUxMEzT44skSjXMIegqcIyBoHt9Ffc7yVi5qruP0hfVomrVcvGWBb38EaHHaGhK01VtFHwxTsrlnzCo+YYfvjh/zEvJ1i9w2D853hYi5XuM5o5tI9uymuUa4QtzUEr4CI1bLh2xRpy7TwxIGOG2RIzhj7j0mpGlNlvgbRvvuXWu9NQizRMnwG/0isDxEN6ov6CbpbD4gkqwy++UtNMJNvTWzYFVaRJBc+nf0HfNyEJoVYmca6KZp55MJd5ze+LyJlVMXekVuhDQoyeDY5zfWEK+zl7HviZae5dSP77HfsiaF4sUUZ7QWWNhca3k57Rw0v2nSkIghZInU9iehMGY/G+3cQk0L5pmB3dcuQqANbHWr3LoRlL0bOKc5xTnHtrG4tZ6lbQ3u9R7XBBINPd5IoVAIekRNydrOJIPjOetZJmJgiy/q26B5ESJWY/VYtI+JFAKnfo00dPS6uXZlXF/uHz7vvFsZVmMobfdLEwLTdDyA5SGOhgxO3BlSsqUnyciOVVRDCbRpIDmWpPsH1zLPGGLfq37MRS9+8UwPSTELaK5L8NK3/yvF827mA/H7WNZ1D2/78dO89jtP8cdnu4NNSxUKxfTgGA/Cmo1fNLCcecNrIz1oYM2kCttw07RQCKSw3vPPuMdGdgKQypvsG827nyU0QbNtuEmhIU1vpt4xVp1iFb1jOdsAsgSaU6ghUNzArrzndyK1NdSwtLWWeYuXsezMizj9mBY3zBCgIdsNwFjrCynWzkWKWKCAhxTCConTLI9RpuA0K7aMXqfcuj9PyYzVel4eoWHGainVzQkPzfpeKTrc2zLMfY2a7SIJjicsEdcYXHCpKwyzDUuhzQtJo5SD3o2gF90cm5F559O3+OXU1flCmfyGo10kxAo9tXpaYXqly0uJFvTmY+lL5TG1OM11Cc6ZZ3LyQqugRk1pjDkdDzK0fRVCGpyS2UBibC9QXvZeCMH8ploSmnSrOQcKPwhLpLkeMOFNE2glXznziH5io9miKzLaGsrz3a0+d15BhdrGttASdmiafZ3UFC2v8nFzGzi+rdat2NlY64kGARCzRNQpC5uotS3gXMlgU1fS9r5K4rG4bQx7xrwmhN17znctD2xlXkPc6yulxZnbaG1PmAb9qTx7hzKIUpZ4KePrP2UZ4FIaCKmTq5lP4ZSrfcdVC4Y4ajXEzBLtwxnMqDzIyGIPFqs2bSZmZGlwJ2gsb6twPFS+81mft6qG1hSGQUpiRhEjVkMyW6Qzl8CI1wOCjV1J+pJpNnQmLc+zcASaJwKt4ydZ2Pckc0c3ersmTVeYO2dFCBB1LXhn0xpWwi6B7/T6O3ZoBcdoKeIxDQ17qiDQIsQSSkKa1Of7rWbcQrNz2OzQ1ZBX1Mrj02kYb4es1V9uS2+KLT1j7Bm0PPSanZdLqgdNeILcTzxmLaPHG8jrRuA6Gc4UyRbtiru2IBdOQRXnfo7XEHP66aGhaZo7YWIaJQwRQ4p44Jo0TLuxupScubiJUxY0uZ7JgZTl/TSlRAv57J3qvnqEJz/Z8Rx7djxX9r4fJdCmmJFUmp3fu47TS1vZdfnXOe+yV8/0kBSzCSGoueZrsOxlfFb+gLsuH6VomHz8Nxt46Vce54d/200qX9r/ehQKxeSg26JDShrikvlNtdQWhpk7usESXGWGm3B/rIW/6AfYIS9awDgThiNqBAOpnOU5QYBeIF5KU0o0Y8TqAqJONyXdySzddvibJoQlJmwDNGabWx0jWbJFu4R0vIGFzXUB06CuNAqlLJqmMXfxMuoSsUC4nSMK6hKWAWhq8VD4mV2qOha3wrQcARkqCOHv7WVqNbTVJ9zvA251OU1YIYKZRktMOf3IgrkwJomYcI1T59jZBxywjEWrAIbvi7XNrgKUSBjaCem+QIhjS12C1havAbE765/w54ZboVvhwAbHILX+jlmLpXpoSHjGc9JXCa/G56UqC3EUgqbaOK86fQGvPnORs4XACIQI5h7JtuP9h6AiphRomQHrGsslyz6PaXiFYoQGda2US+cgx7TZoYv/v73zDpfjKg/+70zZXm7dW6V7paveZUm25Cobd8dWwIQSIJA4lFDyPSQEAgn5QkggjQTTCeQD4wQDIQkJDgSDjY1xAcvY2LgX2bJkWf1KunXbfH+cmdmZ2dm9e6VbrfN7nn3u3dkp77w7M3ve87ZykXxRXm+JaKBwhceYCco4li9QxkLXhfRi+fLIQNN1XzPoIjqZaEWsshZxwx2dgbT7XlRy2ypFT6QHbbRo8ci+Sn6dDHEMeNCcJsZaxDYChW99nx723Sa9wJZF5qV7aTv0c9/+QFZJtCw4POT3Fi/SD9N+8F7iIy/KXmmByozONf74vkGPEu1JCuE10OR1bxaHiA/vqejRzsFy10MavomoAVieZ4Mlq0diVwgtl/1fvycM23cLBL5TozQCljRmdE3gbyBekSk7+Ct3+Ui+yEih5Jb01zUBZqUID1q1Qaxr8llQ0qOyKJBHqHyx7PaFTIzsAatMPiPzzIjZxYCMSnN3Swh0oVWul1KRMrpPbwDHx4o8dWCI1iO/oPXogzQlIu73c2JMnsOR4XxVH2PnUijYHrSinnC9/M5ERz2UgTaN7D0yxIOfeh1bCjt58sy/YPUlb55tkRRzEd2E196I6FzLBQ++j1uujfGVt2yhvzXJx773OGd//DY+evOj7Dk6MvG+FArFqWGVGE30UCxbJI2SO/sP9iA52eZfXWjuCETgz3uSHjRPBTMhEGNywNWZjZMnApYdGmPPUB9tXg/IfIXKwMhi37ExRgpOo2PHQBOULc3nBXMo6THa0hHf4NhwRyJlNzfDayykhqSHZ3N/C9m4KSuZlfwGmpzkloabnIevDBwBmdfioTmdZEnOCTmTx3ru8Lh77JIeZbBpNQDJkKIegrLdQ0tUPnOKJXg8l4tzabe5b9FIgKjMZ7v7KxV8hm9zMgKRJPs7zq80BgZp3Nl0ZWNENEFryvYQud+Jf/jk7YV3sH0bQZxiIflIM7oGR5vXkzez/pXKnkIMXqNeSL17B+bl3i3QsthdfzTe6dtVSzJCJiZ7NeWT3XYeUPV1ogmBfnSXfRwNwzQ52rymaj3vfnvbbLmtMgXbQPN50ITwDa6rJjWe+4ksz645g2HvucqGxL7cIqET0+VVdrB9KwUz7cn1KtvHsLf3Hk9UrhtpsOhVE55+D5pRmXQQRshkjB+9NEp0/DBmoRIO6za+tvejlWXFwCdeOuYb+McK0osUH91HJH8Ui6Ax4jdA8XjQvJMV8sKwj+XLjyxLA9RVjLwbIobO5r4WaXzalRedSSMsqXfvfjRRqWrpvS+DExZ6KQ9YlMpBryyupeLNjS0Gti+YaayBi8GZhMj2+u7Dijy2LtDYOzjKsKcyYqEkDTRnwkEvj1NKdcPaV7v7Ep7iH5bQKCKLFTUdfZhy0fGg1faUmm6Smfw7ki/ywO6jlOwKkFrgq4FKDpql6RxrWllz31Xn2vCaiknx9P7j3PuZ3+Gi4p28cMYHWHHV78+2SIq5TDQNb/g2ZLrQbnoNF7Yc4aa3beXm95zLxStz3HD3c1zwd7fz7q//gl++UD0LqlAopggjyo8PZxgtlNCtYpW3w+o7l729V1HSYyxuS3HGwmbXYEAEqzzKEs4VA013S8FHDJPBZD9DrWsZT8i+PoVYO4VIFrDYdWi4Zj6qsL01IMAqVvUqA2mgabbHy8FbYQ4j6srofu70ONN1lnemKWumL9TaCXHUNMMTRiR8A894xm/ACt2wjxNQjPfYVXGj/vM29crg9MXBUfLu2E8u0wSkohGG0ovZ13UxJSMJml7Zv7O7Fx9AO/S0u11T3EQIWT2vpHvC/zwDw2Q8yrbFLUT16kqdjph6adynyKJRXZ3ZKW9eNJKMr3ilbLXiuW7kzqWuo7pV1RsO/GXyBUD7cvewQ6lFvnVNTUMTGhaCkda1odcISN26hxcCTUBJr1Nd2qIid3GM5NDzACQ9vVzlGNpb2tx/8GJJ5jppmk7VMFTAkZFA3zxhoGMhNI18tBWgYjDblJ0oQoHvfJwBuwzDC5T+9zSKB7A00zX4LM2UOvN6mMtFurL+6n5Ng4+RO1DpaeqUbE8N7YKHv41mVyOM5Ac9/ebixE7IqqVOuCMBr41zTzkFReT17/qifArz9mi7oC/KluZhzOIJLKH5SuNXVWHE1pkd/lf2/O/gvW4sC1oSEXqbE1XVIJ1iLRZ28ZdySXrRju1xpc0ee4xS2eLYaJ6C5T/ffLSVSLqlYtjXaMqtOc8b+7s8NFwpSCRzB/1yJaOBKo66972GsD1oTn+9sVRvlQfNe2xD908OeXvfCSF8h3f0Pe5MxFjV8tVDGWjTwC93H+Vnn38715Z/wKH172TBNR+abZEU84FUO7zpP+XA6cZfh8PPsKYnyydft5GfvP9CfvfcRdzxxEF2fPYuXvOFe7jlkZfCy28rFIqTx0y4DagjhWOBaB+BZv/AF40kzQmTTLwysNdLY/6QQSHztSxkr6Hj+bJbrCNi6BTKMJLuc71ZhZis7OvMvg+P+71XDu4Pd7mIMXYkNMzNW+nPwfHiYJXdgacz0PL29xJCYGgaljDcXkDuse0Qx+j4IVJDu2SVRXswbuoaLU3+HCZnQBRMoK+ckieELCyUCnwetNFCiZ88fUiu53jQhCAWkQM2p3oinsHSi4OVXD+vDykbNz2Go0e2VA56zoAFZ8pwR09fJdMzQHO8PKPxzoAB5B9adWfj7sy7JXS3/5eLU2lucDeMHOGi0t1syx7xqUnDH+KoCQFG5TsOr9AIo+N5njlaCB0WrurKBMJc5ffufJ/ZePUg2dQ1aLaNwcPP0GTIAXLEE9opBfYXy/DiTFLohlY1GNYQgYbPtiElSu69B/JaXtVV6WE3NF5kzO5D5npOhS7PqTgsQxzjHQTxej6Fx+sn7x/hmywoGXHpdfVuX/YXDHHOJj66T55jQUa+OM2Q85FmtPI4wfQqK3CzLrULnSTsHoIyn86ZJBCMxjoq23lkbNpzO21HHsQsnKCoJ+x7R2oliDykNz/OQpTy/jVFpUiRZcl7sTMTq96f3XcNkNd3uQBPfh9231v1HDo8nOfIaEXmBc0Jzh5ok9eWY0TrRmhIrmaHOIYZUUXbgwaygu7SXIrlnX4vnFNxF2BRe5rlnWlp/FtlnrfaKepJnyHvneBqSXgLzNgTLp6CMpqQ154T4u1sO66nONK8gSMtG6u+53ooA22K+d+H9/HQl9/BG/gexze8jbZf/9hsi6SYTzT3w5u+I2exvnoVHH4GgO6mOB+8ciX3fOgVfPjXVrF3cJS33Xg/r/iHO7jx3ucZzZfq71ehUDSMZc/eRkde8pfVFxWPVNFI2jOmAj1v/yAHvFlynGrncxlxDg9VZnsNe+C05+iom0NWqc4mf/S9xpF3AO4utcoyRDLkR98ZYHo/qyTd23sQGrl0lEVtSToyFUNTGmgyr8vvxZPLheEJZ7PKlLUIuXSM9b1NZJN+D4OTpK9RySNzPAm1QsjKWoSjzWsBwfHMMlm5zTcg8w86BZCO+wfOMsRRMjia5+CJcSzL8oWKZ+KmpzGttz62Di2LoWmh3Lun3HlzIkJ/a9KXg1Y0U1iximFqBU4rEzfRnV51HrncYzrbjg3CM7e5jasB9necL69BIXwhjpoTRujIH2wMbStGs2QobDGk5lQi5i8aYgmNiCHcfKiWRAQzcC69zXGZz2N/d0tyKdb1NPmMzmoPWmX7vJmtFJMRWpVXS9obFofazvTIpaOXiwjdbzCGXT8aeHrV+XMX85FKSOmC5gQdmZgvhNnrtVnSk6MrG3M9YACFSJPPoHML8XjlDxpe9r2WsY3dgplGWGXSsUDOWWA43pVNuv+3JSM4YX3yIIIjrZvY23M50iTwG7QRQ2N/7jyOtG5yDbRaXiGvfeaU0w8+8ypl+C33s4V203bnJZD5ZwCaEYWhA5VeggGGx0vsPlZ5FuqaqLQoca5jPRKyJWRjMvQ67HzydjsGkJVSs/FIVaERb1+0FV1Z2lJxu6hJmRPjJY6PFdE915n3ee43zh0PWqV4itfHKfvDYa8jGE32UIxkqDJs66AMtCnCsiw+/+OnOfDN9/Am7X8Z2fR7ZHb8bXXohkIxER2r4M3fhVIevvprrpEGkIoaXHfuIu74o+18+vUbycQMPvydX3Hu39zGF+94xi0QoFAoTh2NyiN8NNZBcdWrKnkFpjMzKyjl5Sy6bMrq9aDJwcr+ju08EDuTg8crxoFhVAYKcWcUbA86Btrl4KzkrZQgNAbaUmiawWCnzHEaLRTtnJFq2Ut6rGooEIt4PGgAmoGha7T2rXVDsxx0XWBpBp4JYmKmwZJcCs30Vuwr26Fhjpg6S3MpUnZokTAcD1pl+OIMopriJn2tlYGow1DuDEaSC9nbcwUnMkttL5+n4IGbi+MUZxHETIM1PXIAvro7C0L3/fyOF0uMFksMjcuG0ucsaZP7ddep8VstNKkvT9nwTNy0m4n7C3lU5PMbHZpWaQVgCb/HCoBIEtKdgdL+YGkRiqYsdqIJf/U+ZxeOweBU0vSJ7vm/pJkyr8cnmGM8V7aQ4aROlUyqxjBOCX1nVG9oGhHDX65eCGRuNYBm+I0aPcaQnTekaZo/tNTZPzAea+fF7kvcHDVtbLBiyNo43l+vvsv2PnqbE/a2HmPNU54+HTdY0p7yh/159Nva3FzxOHuJVnq75c1M9efJ9uplyII0qzoz7jWbjlZ7PL1oHsNC1+wQ32Dume0hdHKuvLS3toHtCbc3CD3O8fGi26svUGxSyiFkbpeT6+U83yK6RndT3H1BZUJJM/zXYvDI48WSb7JBeNdyBKgR4tiainLRis6aBqdRkEVghAC6N1Z97jPynQk2YdlGrsZIvuifgPLm6Hq+s7DjV2nYXqVQtkLXOtJSLV/I5opTIV8s84F/+wWtt/4Bv2X8kOLW95D4tY8r40xx8jhGWnFMGmkHHvd9bOgaV6/v5jvvOodvvX0bq7ozfPz7j3Pe38jKj8pQUyhODl8FQSEqg0u76IQzKMz2LJUDgNYlRC05U5yJmQTHdJoQFM0ke48XfNXqUp4Z9NWdctDnDFoM+7fjRMQz2LPK6P1nMdq3nQIRjo3meeTF4xTLFiJk8FXSo76foLIWJebMUrsGmv3e8Hu95KBFelK8RshFK3K0pqKYHm9NMCcHBNl4hGw6w2DTajR70Kt5whkdA00IQV9rZcBbCbTy52Zp9kBqOCkrsml23yyf10FoDLSn2LGhhyU5WQbb+13mi2WIt9gHqBRCcHRXMCtyUPRU3BNAfsi/DGkUeMddmt2wmdYBgkM10bXB9dRYnqpxLkbMZxw434fTvsEdwHoGhcHhRbCSJvibau/rvgQWbvWvYg9WvTlbEUPzGYKOrANtqUpIYUgRhYjHE5UfuFTutO9sWHqpz3gzTJMjdt6QTpnhVL9vPz5DV4tQ0hNE8sfQCicgIfPPeptljpxzL47GOtym2OMFWQBntOdsLM10B9IlLcQjE1SiJ/xNjyZsefzraB49n8gMcCxbKfpwpHkDonNt9XHwFiOxv80JxoeGVxa795eXqB1quKQl3JBJ2aHXjp1X8SQK2bS5ZRF74svZ37G9Ev5rQdkq+76vJtsz7Xieq8R2Kp4K4U4o6UbA6A6RzxeSK6hc20U7ZNQMPlcqpOORmoU8Su6zTPhbbdjo3vvENuAF2EVVNFv+yrXi6KIzE/NPYnnvxUBLDHdb+03R8k4uVf4ftfs11kIZaKfIvmOj/PYXf8wVD7+X1xh3YF3wAYzLPqqMM8Wp07Ea3nJ/8upcAAAgAElEQVSzrIb0lcvhhZ9XrSKE4MxFLdx43Vn8++9JQ+1j33uc7X93O9/4+W6KqpeaQnFqeELghCbDBC9akWNjX7scjJsxFi5Zy6quDObqa6pm3f3J9JUf8ow9G5uIGJWyz/a25dbFAAw2VarpCcqIpoWIaIaCFmO8WKake0pSByjpMcTARZBbyf7ceezrvrgy8HJnqe1Bkh6Yzbcsu5x1YBDk5q151SNC17EyPQyn+j3V9iqDxHwk490BAEeb17mLvAaCZ5cUDTkIdwy0QsRbBTEoh+YLET1ejlHK9tlrWm7RDmffI8mFsOAsmQOc6vDth/ET8MxtgSMFPWj2/5kezl9W2f5YdhVapsMdkFu2ob9hQVOlQbMZ83sM7JLgwmNEBkcUbp7ViqvZ13VJqNHk3aYtJQfNvoG+PVj1HicbN93Bp+Vp0xyP6JUeXzXKn8fitjc0YnuXM90QSfgmLTYvrlTLFHYxlnyk0kDb0ATbl1fWsYSGURyS125zPzs29Lge5kpfNJOEnYM4ViiBEaWUcPYh1ykHQuZCR2ie70CP1CiU4hnbLe5sYzjZ6/tMq+HZ0YWgsOQSu7G8f3IgdH2P0nRN2CGsHkNX0zhrcSs9aT10Pys605yxsJmOTIih074MejYxnOqnaFY82GXLolCyPCGi0lvcFDc5MR6Y9M10w8prYNll9rlUPGhCDz43qkXwhuQKPN7BcbuXnzOZEiRnG8Q19HwsuwqwJxZCrtO4N6xXj7gecuEx0Jb3VgodeYpgTozw5zQ612deRPwrNYgy0E6Bnz51iDde/z988MD72K4/DFdfj7jwQ8o4U0wdHavhd34A8Wa44Rp48paaq27qk4bat9+xjd7mOH/8Hw9zxfV3cutj+2tWg1MoFH685b1LVsU75c05SsdMX3luc8FGEhtfA2asqrpZsDeOu5wy5y1t57ylba5Hy/WgJVpkpUijUlziYPtWNNvDcWi0zM8jZ3Igdy5jsVzVT85g02oQOlqyFTpWs2ZgIRcsa6+EuDkn4gzqg+FElhWS9wXO4EIXHsPEPnYh1SVzaO0FXU1x1vRk6WuJV+khH22VzaRtkhGDkeQCNDPEiEAOULcubmUo1c+J9BKGbK/LuF3RzyeIg6b7ik0UCnkOD0uDqCMTJWN7MH1bNS2AlVcHZu+rvz9hH89noDn/aoYvV8USAkP3FhfQERr0tSbZsCDL8o607UHzDCbtUD7DqhSgqDo9Z+BoRquMD5AeLcfwWtyW4uyBVvv41TlXzq6XdmZIRAwuW91JWypKazLqGrK+X5CAAU3vFgBOLLyEvT2XVxWW8RoQhm5g6hpDqX6EHQp3qG0zB9u3kY4aCCF8xUmcSQJZFKU65wugrJmuhzIeCVZqtA00LWigeYwQe7LEW4QlmAN1zpI2lnWkiZgV2ZZ3NXHNhoUszTkecA09pLnyeLRVNniPZcDjQfU5ZAJ5ZN79yHvRZDzaQsFMMx5tq2xbKrC+N8u6Xn/4p6YJFrQkfBMDtXC8WRay6IUROAfvvVuINstw3AVngRGR1+3CbRzt2e6OM4yqW6b6WeILcfR411m4DdJd4DWQvRd/LNhI3TkH+b232bl7Tn5wEMN0rnlhG2hCGmiUsdBY3Z0lk6w2zn35jgu3sbo7WzF+ve0rvIe0/y97FgbzU+uhDLSToFy2+NStT/Gxr3yTf7U+yCrzRcTrvw6b3jLboilejrQskkZa+zK46XWw8//VXX1zfwv//ntn84U3nkGxbHHdDTt5/Zfu5aE9qjy/QjER3rmMUtny/MZr9fNG7BGTETTQPG870p6BfzRDSzJCzNRxegZ5+3o59Kw5l0NtZ8oCBUK4xSNKRpKyHuFw6+aQcYg/jKqvNUlTIgLty+2TDIQ4Bgfclmw26y23Liusyf2ZuWVVpz/SeSb0bnZntoWmM9CeQotJL0lw0F7SK6FIF6/qYE1PlsVtSVsPfnk0UfHsHM8ud4u4+At7BIYzQrNzRgQ9TXGEVWKsKOXPxALFNuoRNrB1Bl7ecZljoQVm7Zd2pNFtD6VpG73OERO960l3DkCizf8dtEgPajFa8SIYul9Ob2PmIL3NCTo8TcoNXVRC6ryeNlG5Tjb3tbDU9kzFTJ3+c16Dlu1hYNOlrOjMEN/w6sp2weulWXomY1EThM540V+0yiejXslXrHjAIuSjLRSSHVV5Q25onqbXnPwuawZC01ndnWFJewoQdGYr19dAewrDCC86AdCTk8ar5j2vwPeeihqs7MqAN2dO00HTyboFarTQ6+lo81rKa66V50OlFYdWx5vinejRNEFZM8hHWzjQcT7jsbbKccoFDN1fTt+LCPkviGMEjxYtCuWyz4MGuM23D7Zv43DHOdB/rv86z/ZQiqRl6wg8+XMLt7rXSjmQI+kt3iLtM3ubTBf0n+MXcOll0CX7Q5KQ98R5S9vtapKSwebV9r4c72QNNIN1vU2s7c3K68n3DBEyxzbkWvGFpEbTLMml3OMbxWFn5zJSAE1mtHkNT//ZNoQy0CbJ0eE8v3PDfTx961f4TvQj5FI62lu+B8uvmG3RFC9nUjl4880wcBHc/F74n/dV9SvxIoTg8jVd3PLe8/mLHat5av8Q13zmLt77zQd56dhYze0UitMdx4w42L6NUrlSgCMVj9b0hnkJetB0zw+7E4a1t+dyf1ihYzCFGAOprqWMx2QumhB+D5+z0Jv/Uxd3/+HGhFceZwZ6JNHL+t4m1nRXKpBpiWZe6txun2AwB8YxBOxjJVpAj5AvViqsSQn8Mg+0p9xBZjkwiKlpQ/n0VR3imIwabB7opDkR4Xh2OQXHcRgsaAE1B7hhB3eW+DxozsJAqGFfa9rVs6Frsuqlo4ZMlyzlb0T8Rk80Bat2MNghc8YsrEr+oHs8+7vwyOd4JZviJmLZpe5y7zXp7aXm5g+5K3oGpmYc+s4mnm2Vpcq910qN62ZFZ4b2dJSurD/01me0tCxmPNnLifSSqntluHOrncNXwW3g7PXypmXfwIqlY8KqHcSXXCC9P0KQihrk0jHZ1kDYRVJsXuq8iNLiV7jvS7lV7O29Ck2TuoqbetW96Mrau7mysKqvmqi6XASy+qapa67n1b1UGhyry/o6Gu2pKD12QQ5XnlJ1gRDf8d1j1DHQNNOdzJD7DnjQnFC9aAujhfCK0QLBicxS9vZcjuZ4oDUTFm9npGmp3ztJiAetnjLMBLQtlQ2nTXn+LcmIG7b7YvclbgEYTcDh1k2MtqwI35fQifgMWnncFZ0Zzlosn7Oax1MrAn/pWOP2SXTuftMumnS071Ly/Re456oJwZGW9b7DOxNswedfGMpAmwQ7nzvCjk/dwfnPfpJPRT6LueAMtLf/BHo3zbZoitOBWAZ+85tw9nvgvi/Bja+EkSN1NzF1jd/a1s/tf7Sdd24f4H8e3seFf3871//oKVWaX6GoQz7S7PP6aLWMmQD1QhxdL0gwX8gesccichBpaIL+1iQr7T5V3n2nAl68K9d20Z6qnVDvwxl0OoaSOyAPhkBX3qfiUUxdk4M2zyCqrEWIGTr57OL6xwKIZWlPR30GWpWxCJCVBkY5oB89ZOBbfbxqA839u/bVMt/GTtb3GgzOd5yYoKqe/1B22JzX3Zq0c55sg3WgLcWitqQ8vj1g622OYxRHqsLwpFBOPmDE3U8sKv8vt60g0rUqVBbvuQw2SS8CRhRiWc86lfV9XkcjWsnpya0MMbZrEPSg2cQjOmcPtFVNGPhsHT3C8faNlANFbCA8HNiV1ytbphtWvwoissT/ucu7pNEYkN+yLLdlgPe8S0YckfCGylU8ievXb2FVf3eVgebq2Qh40Kpk9Z/D8s40G/vb3PtYNlmuNq4B373uxRuO6ejI3bZjTdgmVVh1DLTFuSzJiOEaDS0JvwfJ+0zLpes/a3TdILrgDOmtT+Ug3sRIU7WxVJ3fWsccCaumid8D6+5GCMbineRbltcQ0L52bUPPOW4qapCMRWrK4oqQW+F+f05I5+HcVg63bgY9ip5sw0Ia+aJrHaOeUG6onPeIXfCoHo0/kU5jCqUy1//oKW66/Rd8Mf45NusPwZlvQ1z2scYfaArFVKDpcOlfQm41fPf34QvnwrX/DH3b6m6Wjpm8//IVvP7MhXz8+4/xjz96km/et5sPXLGCa9Z3T1hRSqGYzwghPgT8FfBZy7LeXW9d1ygTgtaUZya1xiAhiDfEUdeEr9KfHkyed7ANl9W9TbSORWhKyFcQTQgWtSXJl8o88ZLsvRbs81OXoIHWtUHmP6U6Q+W5fE0nxksH4OhRZwfuKpesWYDBG3lg7zDkR92crqpjAZhxepsTdPd38d1n7AqMzgDNa+R0bYTcKsq7h4FAmNxEoUG1DDTPR0VL7sWrMdM2JhzPRMiOK/92roODj0NJFkzwRW12nwGscgfwbh6akwuz5GKyL/yc7eu3QpjH0xnse+Te1NfMi4OjJNt7KJTKcGRfPemwNJPjmeWUF8mCK4543ue7b2CsGbIs/IHH/IVRatE6INu+eHXdf96EmzmhfAIBWqX9QdBACQsPLDuhsMGeWpoGpYI0BiNOc3K/V6Rsyfw0IQRW2V8wy3c9eY6r9WwANrjfcS4dZQ9+I3dlZ0a2CqiSV6talooapFrS9mGE78jBzS9b3cmzB4dI2pMFa/s6GRsdqoTUenoeurpqHZDf3ZP/W9lR29Lq86xz++SaM+w+cRiQ4bGRiOlrK+EYhaauyQqpdViaS0vjp0Y1S4eBhb0ceORpj4yTH4OEDVvcSN5aY5pYkwyjtUMlCQtr9dyDbmppiHzOZ2UzxRgVz97KrgwJUUCEVZoUGi92X1oVxh2G8qBNwDMHh7j283dz/x3/xa3JD7NJPAE7PgdX/p0yzhSzx4bXw3W3yNnWr14Jd/yd74FaiwUtCT73hk18821baU5G+D/feJBrP383v3xB5acpXp4IIbYCbwUeamT98UKJlmSEa9Z305mJVQoeNuhB62tN0p6KsnFBM+mY6QtxNFZeyYHcudUbxeVsvmnGWNBSo3occpAohPDnsk0G10CzT8qMQfeG6hlq20CLGrrfqPScSzyiY0ai9LclMHXNk/PjmgWV7bo3Qtd6tHTOHVgucMrre71qmiYHd4FQxKIvFxB/c+GJzhVPc3HnUJ7NMzGTi1d21B54+g6cglXXwOpXAv68OqFpbuhTqBzxJlh2afg6IL0NIAut2MRMncXtUq5gbqND0Ot0IrMEEXWqZFYa6DoUfBU0dUi2wZpr5d+J6N4ow8y8JGpU2wshbtr5ibZO9cAgOuh9Bli7YQurujNuo3AfThXMeHP1Z0gPpwxxFFW/j0IATgXDsMGy/b0taElw9Tr/JGZy0RY6cv5JjYG2lAxnDTM0fOdptz9oWxqosCi/q8XtKXd5dOVlZOMRd48WFgPt0rO+sstzHUVTleIZAxdV8rV8hxbBBRXizViWhWYV7OqHfn04c0Bh34+D00stbJ3gNXq4dQstrW2VdScKcaxBmDiizmdyBSGNWuea8VVqcWYOKs+8sUSnf8cegrXXhJC5npmY4Ybahh3e0syGzlcZaDWwLIsb732eHZ/6MTsOfZmvRz5GUzaLuO4W2PiG2RZPoZA/lm//iQz1+PFfwtd2wNHnGtr0rMWt/Pe7z+Vvr13H7iOj7PjsXfzBt1R+muLlhRAiC/wrcB1wdILVATlL3d0U94SxOftq7JgRQ+PsJW0sbJWGlvf3P5ZIB0rD2/RslgOrOr1/pAxSCCdUMrSRbj2c8DmPERCKd+TRVl0QxEsuHePKtV318+B0053Vv3p9Nzs29LCsyx4gheTSOkdf0Zl2C4R41d9wzh2AEK4hULBLqAVnw5ONhje6/clsOT16ql1spNELJynLlnesDt9LAxegWy4/6EwMhEG6BRtcr90pRFDU6EflJWbq9DTFWZKzPUmuXP71wgbVqXQTCdOohKR56T9f5vAF7xt7PwtaZLNqeb2EVDLuXAv952KlO6s/89xbVaGXrQOw9GLfouZkhFU9svBESzJS7VFGXi/CstycsrZUlM19dQxcIwpLLqZo58tZloyGWdaRDvGwh+eVVnRdrdwWx8vbOkBRj8nqkEJUGWhuKGGdatDlstO4PMRAE4KxWKV1wpq+HClPuXunz+Fk0ZoXuhVdHWp5Z2viK5pT7cUetJtJh+2vIyPPYZFd3CgeDFH17KehSaUAKsQxhBcHR/nQfz7Mricf5r/SX2Qg/wSc8Vtw+V/Lh6hCMVeIZeDaL8PAhfD9P4bPnQ2XfAQ2X1czbttB1wSv2bKAK9Z28tkfP8P/++kuvv/wS7xz+wBvPX9xzXh4hWIe8U/Aty3Luk0I8WeNbJCI6PR5vFjO7HHsJG+HhgYKujEpT4RjmAWrrVWwOGNhiFdB0+WEzkSGnderZcakYVIcq54yDt22QYvWyeWpU+SgJRnh6vXd9m7lfqOGTkTXgRCPihdncGQm3AF2oRxuwEywo8q/toEbNN7r7rNebk2QOpUGvVy4Ihe6XOA3Q0K/LqFzIrOE7LHH/N/zydLgJIG3cIhb+j0gX9i9omsCll8ZXpgk2Qq0Vi+3v7NFbUn6UysQT+1mJBYSwmlEZMn40doFtyaFnYPmNM0OIq+XEA9zPeJNCCsPHAwzMSss3ApHnoVoJvTjsKOdu6RNyqQJBhdeQvHZ+xBitMrodoyuYrm2BM5nYR40XYNj2ZWkhnYB0NOSksfw7S5Ewr6zYXyo5jHFwq0cO7KXZR1pntwvQ74nna5R1bga381cxslbrd40HTPZsUHmmMVMne5gmLSmc+GKHMdHC8RNnZ8+fWhSgZzKg+ahXLb4l3uf59J/vIPeXd/mR4k/ZbG2H17zNbjm08o4U8xNhICNb4R33iNnE7/3PvjaNXDoqYY2T8dM/viKFfzoDy7ggmXtfOKHT/KKT9zBd3/5ouqfppi3CCHeCiwBPtzAum8TQuwUQuw8dvSwrw9QNh5hUWuSgdaTCyt0csRqDdpOhpip0d+a5KzF/sHpuKcse81QyToD6lXdGdnTqapkfSXIaspwvHkhBppr4/lEEFy2upPty9vpaY5PnHsXScjn4cKt7qKxouXu66TwVncTwl9EptY+J2OgNUiYd0bK5H9fyxQYSi2SjcEn8qROE85g1wLW9za5IX01DbRIomYPNP+ObZ+Dx0gRsQysfhWjiW7fqtOSdi206lDMFVe5/1qWBZblVnN0MDV/r7yq3bo5dXXuv2hahjZWl5GUf0T1tkII16CK6Dpm4Zg0xlLtlZUWbnXvtWCrDC/OZ8F2EFDxkLUkIjLXU+g+w0jqI+Q+yXTL9kJ12LGhR7Y/8O6LSXjQ6rRWANmoG6orWwZZ0JLwGKeO0jUyMZPe5sRJPXOUgWbz7MEhXvele/n8d37MTbG/4S+1L2Iu2IT4vbtg1Y7ZFk+hmJimBfCm/4Srr4d9D8HntsGtH4X8SEObL2xN8IU3beKmt24lEzd5z00P8BtfuIf7n69fKVKhmGsIIZYDHwPeYFlW/TrUgGVZ/2RZ1mbLsja3t3sGJz1nANCaiqJxct6GiKFx2epO2XfHpuHBQw2EEKxf0ORr6HvKZHtJ9Kwl27cecoGKgU2yz1WweW84DXoI6njQnFYCwUFNzNSJmTqL2pJcubZrYlGaFvoG9mV7MDY5B5rXg+bXt99Aq7X9zA2znOvKEXk4t4kjzRuqjREhZBW5WcqjX9fTRGsySjZu0t+WdEPtwgb3kyKWkT267PvWRdPckEC3f5y3UXXMoLc5waa+8Fy2CXH0KDT/RH5upS80UwiBwCJi6HivwnWbzmHruRfV3H3Uzt3rzEx+kmjC4jo2q7ozLGmJyZ5u3vs/21vVuDqMUh0PmnNdpmKG9KRqujSMvGlxU2Qxu+0nGr3tRIgHzcPGhc0sOeMVVdVz6+PMMHkKFQXka4TTPsSxWCrz5Z/u4pM/fJw3Gbfyr8mbMCwNrvoH2PTbk/iWFYo5gBCyYfryK+GWD8Odfw8PfQsu+XNY9cqGrudtA63c/J5z+bedL/D3tzzBtZ+/h/OXtfPei5eyMSxsSqGYe2wD2oBfeX4QdeB8IcQ7gKRlWeMT7qVlsfwB33PfKQnjDRe+fE3npA20pbk0+4/XyQ9NtAJHKJgZouNHKGuNhcr58HiaquhYDe0r/L3bThW9toHWmYlxZDjv9o2bKpyKgA20s6uQXQDHX4Rll1V9dHhYXkLZuFl74DWFrpq4qYf2odrS30IqZnDXU4d8y/OpHkYbnKCbVtqW+iI6sgmTc5dWipIsbkuSL5YZaJ8CL3NYPhmwbXErI/kidz9zmFLAEyyEOHnjDGTBiaED8rvWzepCKjbZuMnKzhRtowk5eO9YA0YE0VKjVYVNzNS5fE1n7V59DTCRoaZrgtya7TB8sNKmYewYID18E1G0q2SGyegYbe58hgg0Hbdz8k6FzkyMXCbGMTtctVHDNDTE0UPE0Ih0D8DhByYvlGd/JzMpd1obaHc/fYiPfPdRCgee4L+zX2PZ2EPQf5H0QDQtnG3xFIqTJ5WDV31R5k5+//3w7d+Bruvh4o/IfLUJ0DXB685cyDUburnxnuf5wh3P8MrP3c2Fy9t57yXLWNfbNOE+FIpZ5DvAzsCyrwBPIT1rE3rVXLIL5EClvUbj00kSNSZvdKzqzshKdrVYdAH7Bvspayb5SFNVONcpI0TjxlmjOWhO3l1IEZKlHWn6WpOTKwbSIINNqxFdDVQsdGhaANneqvOJ6BrjRWksrekJKfziMIUetItXdoQGmTq5L7VDHKu/i/OWtlctmza61stwypCCMCAbeNfV4RQQMTQiRoSmhMnBExPPzUyKhdtgaH94IZMAnZko5O0Kf+01enWFcDLPDZhk0YxYRr4ABl4BVsneh0AgGjKgw+SsKqrjhBX6QplPbSLDCfd+YLesBdXwJIzXQIvVuQZbB+RvwWTw3vvC96chTksDbffhEf7qe49y5yPP86HUd/nN2M1oVgJ2fBY2vGGagpMVilmg/xxZ6fHhf4Pb/gpu/HXoOwfO/n1YeumEHrVExODtFwzwhq193HD3c3zpzme55jN3cfHKHG+/YIDNfc2qh5pizmFZ1iDg6x0hhBgGjliW9atJ7UzToGvdFEo3DWgaZTv8MNgYdU4zUDusazqMM4DhVD9abpI6CnnGbRto5fYnDsiPw7aJpCA/FF7c4iQJa+Tsx198I8xW3r4sh6bJ3ONTwkxAYRLeuXqD3xlkS38Lx0YLU3t96aY04hvBMU5qNPueajIxk85MjIV6Co5PYkNNw5sFtamvGRowosOqFSYi8lzdCo+BcYfg5Ko4hlGuU+4/FG+Io6caaEc6RsnbO6974+SFCenF6D3N3ubaLVXgNDPQhsaLfO7HT/PlO5/lGv0efp65iVT+IGx4I1z8fyt9SBSKlxOaDutfJ3v37PwK3P1puOm10LYczn43rH3NhOW9U1GDd124hN/a1sdX73qOf75rFz/6wj1sWNDEW89bzGWrOxqKU1coFNNDOmZwYmyCyoYzwswVFjpnSRtpa3YG/tm4SSZucny0ED5J1X+O9LzOoGHSkYmy+8hI3cFpNjFFeWdLL3U9LPMJU5fl7WeN9hXSKJggrHGqiBia9C4NDtsG2kkYQpFUwxU/wyYR+lsTmLqg7QW/3o9nlhE//ow0WupUc50MThhlwy1IHCMqkI9Zrx9lw4TkoHk5Y2H9SKTTwkAbL5b41n0v8OnbniY39Bjfa/42S0YehNYNcOVNsGDLbIuoUEw/RhS2vgO2XAePfAfuvh7++z2ykMjm35aVICcI7U3HTN7ziqVcd94i/v3+PXz5p7t419d/QU9TnNdtWcBrtiwINN9UKOYGlmVtn20ZppNzlrQxMl7iJ08dnF1B3FLV0z9h05aKwiRDvzb3N97OYCIi9arbRdO1m1JPE+t7m1jemXar7i1ojrPn6Eil39VUohtM5xDynCVtbm+tlxWaDrmpCZeeMZZfPuEqhqa5eWhBhBDSW7THb6YMZZdRECaCPZxqiKODc83ojRacMWIyh7BjzZQc38cEz8CJoo9e1gZaoVTmP36xh0/d+jTRY8/wycx/cXb0TrBa4Nc+KfNzpjD8QKGYF+gmrPsNmci86w64+zNwx9/K18BF8r5YfmXdnjyJiMGbtvXzm2f18cNHX+LGe5/nEz98kk/e+hSvWJHjtVsWcP6y9olLYSsUiikhaugnnacypeRWy4GJU/lxull6aXVp8xCycRNdE7LM9xThPN9qDUxnGk0TbjgZQC4Tc/s0zTdm1cv1cmSaW+ZcvCrHhLfB4gthzBd5znCqj/LCPmiamrzZkhPi2HCZfQ2WvGJKjl2Fx0Bz80EnYYe+LA20E2MF/m3nHr5y9y5iR5/iLzI/4KLY7VCOwfnvl2FdcyQeWqGYNYSAxdvla3A3PPCv8MC/wL+9GRJtsOZaGRa54KyauWq6Jrh8TReXr+li16FhvnHfbr69cw+3PLpfNpld18WOjT1sXNCkctUUihlgQUuClsQ0eEwaxYjIohAzRaxO8RQP25dPfQrDut4shi7IpVXUgGKO4xQwiU5dP0YvDU0OJVvt5uJ+RKZzynLQHDu04Ry06eQUowheVgbarkPD3HjP83xr527WFh7ik6kfsin6c6xSHHHm2+C8P/Q34FMoFJKmhXDhB+GC98Mzt8Evvga/uAF+/kVId8PqX5fGWs/mmsbaorYkH7xiJX94yXJ+8uRB/vPBvXzjvhe44Z7n6c7GuGRVB5eu7uTMRS3Ks6ZQTBNnqFYYM0bM1JW+FfODVA76z5uTtRYaLonfAPX6sU2K1gEonmJeXIiBdlr1QTsynOfmh17kP36xl70vPMdvGHdyW/xOcmIPGK1w9ocQW3431GpXKBQBNB2WXiJf4yfgif+FR/4T7vsy3Ps5SObszy+V4ZAhs9cRQ+PiVR1cvKqD43TIVxsAABBeSURBVGMFbnlkP7c88hLf3CmNtUzM4KIVOS5e1cHZA23Tkx+hUCgUCoWiQrpjtiUIZSqDa8qTDXGsxclUbQzikSEdNRhoT9HfluTupw+F9jMMMu8MNMuyeObgMLc/cYDbnzjInl2PcyH38ZH4g6yNPYJGGbrPgY1/Cqt2QGQKKrEoFKcj0bTMVVv3G7Ia2ZM/gCf/Fx6/GR78V1kmeOE2WHS+bHLbs7nqfsvETF69qZdXb+plNF/izqcOcsuj+7n1sf1858EXAVjZleGcgVbOXtLKlv6WUy//rFAoFArFqZLpbii/UTF3cAy0GahRVBvNBPy99oQQbq+/C1fkKJYmzgkU1iQSBzdv3mzt3Bns/Tm95ItlHtt3nAd2H+Wh3Yc48NyjdA49whbxBOdEnqS3LAd55FbBiqtg3eugbcmMyqhQnFaUivDCz+CpH8BTP4IDjwKWNNi61sOCrbJvVMca2YQ2pNhIsVTml3sGueeZw9z9zGF2Pn+UfLGMrglWdWVY25tlXU+Wtb1ZlnWkVUjkaYgQ4n7LsjbP9HFn43dOoVAoTjcGR/LsPjLCut765eYnw48e3c9wvsgrVnaQis6SD2r8BBzbA7mVE65a73duxg00y7LIl8qMF8uMF8qMF0sMjRcZ2fckY4d3Mzw0xOCJIY4dG6R0/CUio/vp4AgD4kUWay9hImczytEmtL5tstfI8itlvKhCoZh5Ro/CC/fB7ntg972w934o2bNHmiGruTX3QXM/pLtkSdtYk/wbbwIzwXhZ8KuXRrj/hRP8ct8wj+wb4vhYmTIC09BZ0tHEklyaRbk0A7kMi9rTdDUnMA1TNZZ/maIMNIVCoVBMhmcPDvHw3mNctbZrXvRmrfc7N+Pm5fcefol3ff0XVcs/YnyFNxs/rFqeNxPkEzmM9pUY3a+WFmnnOrT2FTWLFSgUihkk3gzLLpUvkB62w0/D/l/B/kfgyLMw+Dy8+CCMHgndRRTYZL9cvIXRDtuvx6q3LaGR1+IUtDhFI0HZTGKZSYimsKJZyrEmrHgLxJsRiRa0RDNWvAWRaEHEmyGWRdc1dCEQmoxd1zWBEN7/lRGoUCgUCsVcZnF7isXt01OpcqaZcQNtRVea91++3O7ZohE1NBIRg55yjj28nUwyRSqZRIskIN1BJJpGlRBQKOYRuiEbceZWyF5rXop5mc82elT2Qxk9CoVRGedfLkKpAOUCWGVZL9cqQ7lkvy8zMl7g8NAYR4dGOTGa58RYnpHRUcrjw2jFEczxYRKMkxTDJDlEhhGaxBBZMVJT3KKlcYwkh6wUx0kyaCU5RpJjnr/HSXKCFCdEihERJy+i5EWMgohQ0KKUtCixaISYqROP6MRN+YrZ/yfsvzFTJ25qJExBQoeEaRHXIWaUiWoWMR1iukVUKxPVLSKiTDSeItLWr4xEhUKhUChOE2bcQBtoT/HO7WE5Yl0zLYpCoZhpjIhsdXGS7S4S9mtBjc/LZYuhfJHB4QKDo3n2F8q8UCxTKIxTHh2EkaNoY0cxxgcx8/IVyR+TfwvHaC8cp6dwnGjheaLFE0SLJxBMEAZelq9i0aAsNCwEWLIxpYWgDAgsdKuMRhlTTFy9ycutpY1cV/gjIoZGzNCImTpRUyNm6EQMDUPXMDXp6TN1DUMXGJrA0Dz/65r9Vy6fDJZlUbbAwv5ryWWWJROync8sC1Z0pnn7BSrcXKFQKBSKU2HeVXFUKBSKWmiaIBMzycRMFhKs4Noz+R2WyzB+3Pb2Dcq/4yegMAbFUd9fozgqvX2AtNAChp2mYwmDAhoFq/LKl3UKFhQsnUJZkLc08mWNvCXIlzWO6a38fnIp48US44UyY4US48XK32LZoliSf0fyRUpli0LJolh2Pqt87qw7WZ0KQBMy7FOIynvNeS9kKqAxF5qDKhQKhUIxz5lUkRAhxEHg+ekTJ5Q24NAMH3MqmK9yw/yVfb7KDfNX9vkqN8xf2U8Xufssyzo5V+spIIQ4ATwx08d9mTJfr9W5htLj1KF0OTUoPU4NNX/nJmWgzQZCiJ2zUcnrVJmvcsP8lX2+yg3zV/b5KjfMX9mV3NPLfJFzPqB0OTUoPU4dSpdTg9Lj9KPKICoUCoVCoVAoFArFHEEZaAqFQqFQKBQKhUIxR5gPBto/zbYAJ8l8lRvmr+zzVW6Yv7LPV7lh/squ5J5e5ouc8wGly6lB6XHqULqcGpQep5k5n4OmUCgUCoVCoVAoFKcL88GDplAoFAqFQqFQKBSnBcpAUygUCoVCoVAoFIo5wqwZaEKIdwkhHhJCHLdf9wghrqqzfr8Qwgp5XT6TcofI9SFbjs9MsN5aIcQdQohRIcReIcSfCSFmratrI3LPFZ0LIf48RIaXJthmTuh7srLPFZ3bsnQJIW4QQhwUQowJIR4VQlwwwTazrvfJyj1XdC6EeK6GHP9TZ5uFQojvCiGGhRCHhBCfEkJE5oHcYeu/YyblriHXO4UQu+zr5n4hxHmzLdNcQgjxQSHEffZv9kH72lsTWEfYz70X7efA7UKI1YF1moUQNwohjtmvG4UQTTN7NnOHsN9jpcfGmOh5r/TYGEIIXQjxUc/zb5cQ4i+FEIZnHaXLGcSYeJVpYw/wAeAppKH4ZuA7QohNlmU9VGe7y4Ffet4fmT4R6yOE2Aq8FagnL0KIDPBD4CfAFmA58FVgGPjE9EoZKk9DcnuYCzp/AtjueV+qteJc0zeTkN3DrOrcfqDeBfwUuAo4CCwGDtTZZtb1fjJye5jt63wLoHvedwH3A98KW1kIoQP/AxwGzgNagRsAAbxnWiX1Mym5PbwVuNnz/tgUyzUphBCvBa4H3om8ft4JfF8IscqyrN2zKdscYjvwOeA+5HX2F8CPbB0598v7gT8E3oJ89v0Z8EMhxHLLsk7Y63wdWAhcAVjAl4Ebgatn5jTmDnV+j5UeJ6DB573SY2N8AHgXciz+MLAO+XsyDnzUXkfpciaxLGvOvJADorfX+Kwf+WVvnm05bXmywDPARcDtwGfqrPt7wHEg7ln2p8Be7EItc1TuOaFz4M+BX01i/bmk78nKPld0/jHgrkluM+t6P0m554TOQ+T6E2AQSNT4/AqgDCzwLHsjMAZk5qrc9joW8OrZ1nFApp8BXwosewr4+GzLNldfQAo54XS1/V4A+4A/8awTB044v+3ASvv7P8ezzrn2suWzfU4zrL/Q32Olx4b1V/d5r/Q4KV3eDNwQWHYDcLPS5ey85kQOmu1afR3yYX/3BKv/hxDigBDiLiHEq2dAvFr8E/Bty7Jua2DdbcCdlmWNepb9AOhGDg5nksnI7TAXdL5YyJC5XUKIbwghFtdZdy7pGyYnu8Ns6/zXgZ8JIb5py/GgEOLdQtQNV5wLej8ZuR1mW+cutrzXAf9iWdZIjdW2AY9ZlvWCZ9kPgCiwaZpFDKVBuR2uFzIs8z4hxDuEELMZch9B6uyWwEe3AGfPvETzhjQyAuao/X4R0IlHj/bz4CdU9LgNGML/W38X0tN+uum61u+x0mNjTPS8V3psnJ8CFwohVgAIIVYhJw6+Z3+udDnDzKqBJmS+yhDShfoF4JWWZT1cY/Uh4H3Aa4ArgVuBbwoh3jgjwnoQQrwVWAJ8uMFNOoH9gWX7PZ/NCCch91zR+c+QLvUrkKEgncDdQojWGuvPCX3bTFb2uaLzxcgQr2eBy5ChX3+NDIGoxVzQ+8nIPVd07uUS5A/il+usE6bvQ0iPxkxf5w6NyA0yNOa1wMXAN5AhsB+aXtHq0oYM0wy7fmdLl/OB64EHgXvs946u6umxEzho2dPrAPb/BziNdD3B77HSY2NM9LxXemycv0GGIj4qhCgAjyA9ap+zP1e6nGFmMwcNZAzrBqAJuBa4QQix3bKsXwVXtCzrEP48lp1CiDZkTOy/zISwAEKI5Ui3+nmWZeUnsWmw4ZyosXxaOBm554rOLcv6vve9EOJe5AP5zcA/1Nos8H5G9e0KMUnZ54rOkZM3Oy3L+qD9/gEhxFLkD1+9gjizrfdJyz2HdO7lrcB9lmU9OMF6tfQ6Ww0uG5LbsqyPet4+aOfT/Qnwl9MpXAOEXb+qWWgIQoh/QIYvnWtZVjCvdiI9hun0tNH1JH6PlR7r0+jzXulxYl4L/Bbwm0jjbAMyymGXZVn/7FlP6XKGmFUPmmVZecuynrYsy7nBHgTeO4ld/AxYOj3S1WQbcrb1V0KIohCiCFwAvNN+Hw3Z5iWqZw9y9t/gbMR0cTJyhzEbOvdhWdYQ8gFSS465oO9QGpA9jNnQ+T7g0cCyx5DJv7WYC3o/GbnDmLXrXAiRA3YAX5pg1TB91/IETTuTkDuMnwEZIUTH1ErVMLU8jzlm+ZkxFxFC/CPweuAiy7Ke9XzkVKitp8eXgJw37Nj+v53TR9d1f4+RhX9A6XEiJnreq+uxcf4O+HvLsr5hWdbDlmXdiJxEdoxfpcsZZk7koHnQkPkTjbIBeYPOJN8B1trHdl47kWE6G4Cw2bB7gPOEEDHPskuAF4HnplNYDycjdxizoXMfth5X1JFjLug7lAZkD2M2dH4Xsgqjl2XA83W2mQt6Pxm5w5jN6/wtyLDvb0yw3j3ASiFEr2fZJfa290+PaHV5C43JHcYGZHGTwakUqFFsL8b9SP15uYSJ86JPK4QQ1yNn2S+yLOvxwMe7kIO0Szzrx5BVRh093oPMN9/m2W4bkOT00fVEv8dPovTYCBM979X12DgJqitMl6jYCUqXM81sVSdBxgmfhywesBb4OLIi2RX25x8HbvWs/2bkj8JK5A35PqRR8d7ZrrRCoBpiiOxZ5IX9DWAN8Cpktbs/nONyzwmdA3+PnF1cBJyFrDZ0HOib6/o+Cdnnis63AAVk2NkS4DeQZdDfVed6mXW9n6Tcc0LntiwCOTj7Ushn7wYe97zXkeWQbwM2IvO59gKfnuNyX40MhVwDDAC/a39H18+03AE5X2t/779rXwvXI/MT+2ZTrrn0Aj5r39MXIWfSnVfKs84H7HVeZX/H30BO0qQ963zfvna3IgdwDwPfne3zm2Xd3o7/91jpcWKdNfK8V3psTJdfRba/ugo5Ln8lsm3BJ5QuZ+k7meWL4XnkjOsB4EfAZYHPn/O8fzPSlT1sXyA7gTfOtgJt2YIPVp/s9rK1yGo3Y8iZ+f/LDJd8n6zcc0XnnodAHjkA/Xdg1XzQ92Rlnys6t2W5CtkXbAw5+P59rw7nqt4nK/cc0/mFyFj9M0M++3PsnGvPsoVIo38EGRb1aSA6l+VG9pt7AFmeeRj5A/5/AGM2dB6Q9Z1Ib6/jhTx/tmWaSy/7Ow57/blnHWF/5/vse/AOYE1gPy3I/M7j9utfgKbZPr9Z1u3t+H+PlR4b09tEz3ulx8b0mAY+iRyXjyJz5T8GxJQuZ+clbIUqFAqFQqFQKBQKhWKWmWs5aAqFQqFQKBQKhUJx2qIMNIVCoVAoFAqFQqGYIygDTaFQKBQKhUKhUCjmCMpAUygUCoVCoVAoFIo5gjLQFAqFQqFQKBQKhWKOoAw0hUKhUCgUCoVCoZgjKANNoVAoFAqFQqFQKOYIykBTKBQKhUKhUCgUijmCMtAUCoVCoVAoFAqFYo7w/wFklsXS6rewWwAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 864x144 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"az.plot_trace(trace0);"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"Auto-assigning NUTS sampler...\n",
"Initializing NUTS using jitter+adapt_diag...\n",
"Sequential sampling (2 chains in 1 job)\n",
"NUTS: [σ2, ν]\n",
"Sampling chain 0, 0 divergences: 100%|██████████| 2000/2000 [00:15<00:00, 126.54it/s]\n",
"Sampling chain 1, 0 divergences: 100%|██████████| 2000/2000 [00:16<00:00, 119.73it/s]\n",
"The number of effective samples is smaller than 25% for some parameters.\n"
]
}
],
"source": [
"with pm.Model() as m1:\n",
" ν = pm.HalfNormal('ν', sigma=5)\n",
" σ2 = pm.InverseGamma('σ2', alpha=ν/2, beta=ν/2, shape=len(t_data))\n",
" t = pm.Normal('t', sigma=pm.math.sqrt(σ2), observed=t_data)\n",
" \n",
" trace1 = pm.sample(1000, tune=1000, cores=2)"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAA2gAAACYCAYAAACRbMq5AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOy9Z5gcV5m/fZ/qNHkkzSjLCrbkiLMcwBgDhrWxCWYJBmMMJqdl+b+bYFmWTSwsyxJMXHL0ggnGOWBkW85BVs5ZmtHkPNOxqs774VTu6gnySDOSzn1dtnq6q6pPxX5+50lCSolGo9FoNBqNRqPRaKYeY6oHoNFoNBqNRqPRaDQahRZoGo1Go9FoNBqNRjNN0AJNo9FoNBqNRqPRaKYJWqBpNBqNRqPRaDQazTRBCzSNRqPRaDQajUajmSZogabRaDQajUaj0Wg00wQt0DQajUaj0Wg0Go1mmqAFmkaj0Wg0Go1Go9FME7RA02g0Go1Go9FoNJppghZoGs0kIIS4QgghhRBvCLx3hhCiIIR481SOTaPRaDSaI4H+7dNojgxaoGk0k8PjQB/g/UhJKbcCtwL/MkVj0mg0Go3mSKJ/+zSaI4AWaBrNJCCltIAHgNcLIUTgoweAc4QQy6dmZBqNRqPRHBn0b59Gc2TQAk2jmTzuBuYDFwbeG3b+nXf0h6PRaDQazRFH//ZpNJOMFmgazeRxH2ARCPUAXgXYwLYpGZFGo9FoNEcW/dun0UwyQko51WPQaI4bhBCPAbVSyguEEIuALcBdUsp3TfHQNBqNRqM5IujfPo1mctEeNI1mcrkbOF8IsRj4EdAPfGpqh6TRaDQazRFF//ZpNJOI9qBpNJOIEOIsYBOwGTgFuFxK+fzUjkqj0Wg0miOH/u3TaCYXLdA0mklGCLEbWAy8WUp591SPR6PRaDSaI43+7dNoJg8t0DQajUaj0Wg0Go1mmqBz0DQajUaj0Wg0Go1mmqAFmkaj0Wg0Go1Go9FME7RA02g0Go1Go9FoNJppghZoGo1Go9FoNBqNRjNN0AJNo9FoNBqNRqPRaKYJyYks3NzcLJcuXXqEhqLRaDSaE5E1a9Z0SylnT/U4KqF/+zQajUYz2Yz22zchgbZ06VKef173HdRoNBrN5CGE2D/VYxgN/dun0Wg0mslmtN8+HeKo0Wg0Go1Go9FoNNOECXnQNJpphZSQ7YWhNrCKYCTBSEAiDXVzoaphqkeo0Wg0Gs2xhVWCkW5omD/VI9FoTli0QNMcGxRHYN8TcOBJaFsP3buUMLNLldfJNMKMk2DmUph/Hiw4DxacD7XNR23YGo1Go9EcU7Q8D4OtcOrVkKmb6tFoNCckWqBppi+lHGy5Ezb/AXavcrxkKZhzBiy+FBoWQP18qJ8LyWqwTfWfVVTirf8gDLRA13bYdg8g1XZnnaJ+eE69Cpa8DBKpKd1NjUaj0WimDcVh9a9tTu04NJoTGC3QNNOP4S548hZ44eeQ74eGRXDRB2HFa2DxSyFVPfFt5gehfQO0vgB7V8NzP4Snvw2ZBjjzTbDyfbDwgsnfF41Go9FojknkVA9Aozlh0QJNM33I9cNT34KnvgNmHs54A1z0AVhyGRgvsp5NVQMsfbn677JPQmEY9j6qPGubfg9rfwHzz4ULb4bzboBkZnL2SaPRaDQajUajmQBaoGmmHtuCZ74Hj35ZeczOejO86rPQvGJcqw/lS+zoGOJgb45s0SJhwMyaNMuaa1nWXEsyESPuMnVw+rXqv6u/CBtugzU/hbs/BY/9D1zxD3DuOyGhbxGNRqPRnEiIqR6ARnPCo61PzdTSsQXu+DgcegFOuRJe83nlyRqDHR1D3L2hjT9v7WBr2yB2hUiMmnSCS5bN4rVnzuPas+fTWBOTb1bVCBd/UHnr9jwCq/4d7vwEPPENuPKflSdP6B8sjUaj0ZxASB3iqNFMFVqgaaYGswiPfxVWf0WFH77lR/CSt4wqhEzL5r5N7fz4ib2sPdCPIWDlkln81atXcM6iRpY01VKXSWLaNj3DRXZ3DbP2QD+rd3bxj7dv5F/u2sw7LzqJD19xCgtmxOSxCQGnvApOfqUKfVz1H3Dbu2H5a+Har6hqkBqNRqPRHM/o+UiNZsrRAk1z9Ok/CL99L7Q+D2e/Da7+0qil7/Mli988d5AfPLaHlr4cy5pr+adrz+C68xfSXBefK7ZoZg3nnjSDv7xgEVJKNrUO8sun9/OrZw5w67MHeNvKk/jUlSuY01BVvrIQcMbrVaXH536ghNp3Xw5v+Dqc/dZJOggajUaj0Wg0R4ZHd3RRnUpw8bJZUz0UzWGgBZrm6LLzIfjDB8Ay4W0/g7Ouq7ioZUtuX9vK/zy4nbaBPBcsnsE/XXsmrz1zLglj/FN8QgjOXtTIf731HP7qyuV879Hd3PZcC3esbeWTV67g5suWkU7G5KklknDpR1We2u8/AL9/P+x5GF73ZUjXHs7eazQajUZzbKBDHI9p+rNF+qd6EJrD5kWWxtNoxoltw6ovwK/eCg0L4cOPjirOVu/o4vXffJy//e16musy3PqBS/jDxy7j6pfMm5A4i7JoZg3/cd3ZPPj/XsFLT2nii/dt46qvr2b1jq7KK81YDO+9Fy7/W1j7K/j+q6Bz22GPQaPRaDTHNrYt6RspTvUwjhDub6wWaJojwHAnHFo71aOY9miBpjnyFLMql2v1l+G8d8EHHoKmU2IX3XxogHf/6Blu+vGzDBdK3PLO87nj45fxsuWVQyAPh6XNtfzwPRfx05svAuCmHz/LJ259gc7BfPwKiSRc+Tm46Y+Q64Mf/YUqKKLRaDSaE47NhwZZvbOL4cJx2MzZzQWX9tSOQ3N8snc19Oye6lFMe3SIo+bIMtQB//cONVty9X/BpR+JXaylL8tXH9zB7etaaaxO8bnXn8mNly4mk0wc0eG98rQ5XHpyE//76B6+/cguHt3exd9edRo3Xrok3lN38ivhg6vg1rfDL98Cr/86XPDuIzpGjUaj0Uwv+nPKe1Y0bZiObTOLI2CkIJk+/G3oEEfNkURKXSF7FLRA0xw5OrfCr94G2R54x61w+jVliwzkSnzn4V385Ml9AHzoFSfzsVcup7E6phz+EaIqleCvX7OCN563gH++YxOfv3Mzv11zkM9ecyYvPaWpfIUZJ8H77leFTu78BPTthVd/Tj9oNBqN5gRDTlcRs/0+SKThzDce/ja0B01zJJE2iCM7CT8d2NDSz3DenHAkmBZomiPDwefgV2+BZBXcfC8sOD/0ccG0+MVT+/nWw7sYyJV48/kL+Zu/OI2FceXvjxLLmmv5+fsu5q4NbXzx3q288wdP85oz5vDp153O8jn14YWrGuGG2+Cev1GNrUe64fVfA+P4f9hoNBrN8cy+7hEk6jehEuJYqEVvvdgcuWkqPjXHB9IGjn+baW/3yGGtpwWaZvLZuxpufQfUzYH33KmKbDhIKblnYxtfum8bLX05Ll/RzKdfdzpnLWicwgH7CCF447kL+Isz5/LjJ/bynYd38xdfW80bz13AJ169guVz6vyFEyl4wzegdjY89hUo5eC676p8NY1Go9Eck6xvUbXvRhNoxwy2dRgTh6PkoB3W9jSaGLSHdlS0JamZXHY8qAqCzFwKN90B9fO8j3Z1DvHPd2zmyd09nDG/gZ+/72xecersqRvrKFSlEnzslcu5fuVJfH/1Hn7+1H7uWH+Ia8+ez3tftpQLl8xECKHCGq/8HKSqVL80M6+abr+YuH+NRqPRTDkjBZPazDFoJtkBw9c2D19QRcM3B1rgwNOw/Eqonnn445uOmAVITsdkwsrs6BhiVm3a6webLZp0DBaOnYmF41igSSlZ3zIQntSfIMfgk0czbdlyJ/zufTD3TLjxdqhV+Vv5ksUtf97J91fvoSad4N/fdBY3XFKhCMc0o6kuw2euOYMPveJkfvT4Xn7x1H7u3tDGaXPrueGSxbz5goU0VKXgFX8HqVp44DPwmxvh+l8ccw97jUZz/NE1VODJ3d1cecZc6o5FsTGFmPYxGuJnl/zXffth9qkTW79SFcehdvVvrv/4Emi5Ptj1Z1h0EcxcMtWjGTdb2wYBeNN5CwF4ek8PQ3mTBTOqjniBtUlhuuZvTgIDuRL7e0boGioc9jZ0mX3N5LD9fvjdzSrX7D13eeJsY8sAb/zW43znkd286byFrPrbV/Luly49JsRZkKa6DH9/9ek889kr+a+3nE0mZfD5OzdzyRf+zN/9dj3P7etFXvpRuParsPMB+O3NYJXG3rBGozmqCCHmCyF+JoToEkLkhRBbhBBXTPW4jhQtfVkAeoeP155dRw57HAJtqk3MohkXhhgo/d++QfXtNA/HUIzsnSvYxHFmOg60qn/zA1M7jokirZDIca8FKadx8Zog09yD1jmYZ8uhwcNa181RDd6fEz0nejpN8+LZvUqFNc47G278HVQ1YtmSb63axTdX7aSpLs1Pbr6IV502Z6pH+qKpSSe5/qLFXH/RYja2DHDrs/u5c90hfrumhWXNtbxt5RXc+Kr/pOHhf4Tff0CFO+qcNI1mWiCEmAE8ATwOXAt0AScDnVM5rqOBnHIpcexhjWJQHc2ivXesa2XhjGpWLp0Ven9HxxBb2wa56qx5VAkLSlmonhEWaAAdm2CkC5ZdPrEvjhrQ7vE43ioWm07/02Ms6mVh6/1kqxcAi5x3/PNyTDh/p7FAK1k2T+3pAeC0efUTdyo4i5uBcGPLliQT49/OcTYNojnq7Hsc/u8GaD4VbvwDVDXSM1zgph8/w9ce2sHrz5nPg5+64rgQZ1HOXtTIF//yHJ77p9fwlbedy+z6DF++fzvn3b+U3zZ9GLb8EfnHj6qkao1GMx34e6BNSnmTlPJZKeVeKeWfpZRbp3pgmumHPY4Z76PlqGjtz1V8r1CyYc/DsOsh9UFcE2BXhIwLN8QxunMy/Plxg7M/UWE7jbEtZfjX5A5574nAaRvPtTvlTGOBtqHF96beu7FtUrY52oRPHHpqX3P4tK6BW69XVRrf/UeomcWa/b18/Fdr6csW+fJbzuHtF5001aM84tSkk7z1wkW89cJF7Ose4TfPH+Q/n02xr9TH3228jf0DJie954cYiWMgJlyjOb65DrhfCPEb4FXAIeCHwLflMRETpDma2OOwH6eNZ7Iw5L/u3VP++eF4h6K3hGdQT84+DxdMBEx9IRbpTKIeQ2kJMthCwWn47MpmW8rp60ELTlhPoUAzLZvHdnVz/kkzmFFTXtStUPLHeVhiN2aVfMmeUG6g9qBpDo+e3fCrt0NNk6rWWDebO9a18o7vP006afCHj73shBBnUZY21/IPV5/OU5+5kpPe9Dl+mbmeJQf+wD3/9S7u39h2bMSFazTHLycDHwP2AFcB3wC+BHw8uqAQ4kNCiOeFEM93dXUd3VEeAY7ZghdTyLhmvKfysDrfHRpnpTEnq8a/XS+EMSrQnL8nKSrkz1s7eGhrx6Rs60Xh7s8xFO1ilfycwpEtD8C+J7y/bSmnb/2NYC7kFA6yN1tkMFeqmGNmvMg6CXGi7pHtE4uk1wJNM3GGu+CXb1GzHzf+AVk/j/99dDd//et1XLB4Jnd94uXTpq/ZVFGVSvCOixdzw99/j10r3s8bivfR+pv/x1u/+yTb2g8v6VSj0bxoDOAFKeVnpJRrpZQ/AW4hRqBJKb8vpVwppVw5e/b0bAcyEUxr6owh07IpWdM3nKkS4ykSMpW61/XeBfNcKAyBkQy1uAEgcRitX6JGZlbl5NB/ADb+Drq2H1Nep4q4HjT7CO9Lrh9KEwk1rUzQg7Z17wEYakM4Hin7CIU4DmRLjBReZBho0PO35+HYRdoGctyxrjW+AM4k4adTxguxF1vIbjKOvxZomolRHIFb367K7d5wG/asU/i3u7fwxfu2ce058/n5+y+msSY11aOcNhgJg+U3/A/2xR/h/cn7eH3X97n2lsf44r1byRaPnXh3jeY4oQ3YEnlvK7B4CsZyVLAcBTGVAum+Te2TlsdxNLGlhE2/h5Y1ZZ/5Pqapd1VYQZW480GVS5WsDi90OL3QSiOB1znfuB5xPMrtG1VftDHY0zVM59DkCJMjgn2UQhx3PQQ77puUTckYoZeQyqaQRyjE8ZEdnS/e4zmOY7yjYxjg8MWgbcG2e2Co8lij9W4KpoUZeEYaAeGWNCYulSbj+GuBphk/tqXKx7etg7f+GLloJf985yZ+8sQ+3nfZMr75jvOPjd4bRxshMF73JVj5Pm6Wf+T7Cx/gf1fv5qqvr+Zpp0qQRqM5KjwBnBZ571Rg/xSM5ajghjZOZYjjVBcseH5fr9duYCJYtlSWXN/eistMqQfNDXGMG0RtM6Rq2NE+yMG+7ATDyRzjNFhsxKrQpmEcpek3tg7w1O7J+617cnc3u7uGJ217vgftKEyaTlYYZefmsrcMqcRP1IP2YlIr9nQNUzAnMfSz0nUUXMTxCL9woK/iMlJK7ljX6vWCC1EYVBMK7Rsrr+9MrLhC7P5N7Tyy3Q9lTxiCTL4Lwyoc1vMrOHHTVHt41UG1QNOMn4c+r3p8XfPfyNNex7/etYVfPn2Aj1xxCp97/RkvOmb3uEYIuOZ/4Px3c2XXz3jikmcwhOAd33+af71rM7nisRP7rtEcw3wNuFQI8VkhxHIhxNuATwLfPmoj2HoX7Hn0sFbtGipwx7rWCT0v3DC9IxkuNF2wbRmaBQclXlr7c6zZX9nYi27Dez2uKo6Tp9AO9mZZf7B/3NsfVaA1LIDTr2GwYNIxmD/8ggxFR9hWEhZjlNxfO4qRfbh0DRXY1FpBGHbvrHx/VdoHN0R0HOJhuiAD5rs7MW44HjRbytDpPtxJhKF8iY2tA6zZF3MOh7sOr7feOESwG449PIoHzb32d3QMlX/ons9Rrk1bli8yUjSRUnL3hkPs7xmhuftZzsytUUVXJngQg7ftzNoUS5tqAciXxv/s1gJNMz7W/wae/CZc9EHkyvfzhXu28tMn9/H+ly/jH64+rWIcryaAYcAbboHzb2Th+lt46Pwnec9Ll/CTJ/Zx7Tcfi58J0mg0k4aU8jlUJce3A5uALwCfA75z1AZhFvwQsQly0PECdQ2N3zBy7YoTQaA9sqOTeyKhlBOd/Q+KsvHYZJPpHHzhQB/7ekZC7422fXessd7RRDTVYCIDDSy7/V7nyw5vEvFA78Q9ly+KtvXx95dVgs23Q0e558nzoE2oFcHUIgPnw50bN5z9sKWkYL3IKoT4+Vkj0QkhKWHvo7B39cQ3Kse+joITDpVCs0fdo3E0VHcnPqKWq2VLz3NelUxQJ9SzdqIRCMFjbghB0dmP6ATMaGiBphmb1hfgzr+CJS+Hq7/It1bt4oeP7+W9L1vKP117hhZnE8Ew4A3fhPNuJPX4l/nXxrv51QcuYShvct23n+DWZw7oSo8azRFESnmPlPJcKWWVlPJUKeUth1NifyBbivdcxNG3//BmmyOkE+onu2hNwIPm7NqkhilNU4by5TPuhdLEhGnwlI5n0uxIP65H23zBEd3jug4nMtC4ZSt5PkYxgo8Eo96qI90ADOSKdA9H7jc396k3Jlw1mIM2nt4K0wDpjFMgvNNlOEVOpIThwL0Qd8g6h/IM5ceXc5ctmhzoCQhtVwCNI7y1jHEcX3M8Am2062BcAs1ZJGK/ep41aSGEMtky+a7xP+sBcv3InC/EjEALhIlEmmmBphmdoQ749bugbi68/Wf8ek0b//OnHfzlBQv5/BvO1OLscDAMeOM34bx3wSNf5LLWH3PfX1/Oxctm8Y+3b+STv1437genRqM5+uRLFo/s6GTdeGZDC8PQ8hwcfOZFf28mqX6yCxPwhrkCbUIetFJOjXsq2PUQ7HiA1Tu6uH9T+4ve3ESOFUzc23AkioRY4wizzBUt77NKxqOUknzVXPevCYxgAgLtKDOqody2DoCdncM8sas7/JljrJcsU+UjDndBtld9FvTqWGNPpGxqHeCJXd1HZjK1lFcVMvv2jbqY60FLGOps5UsWIznlAbSlDIUHxl1DT+3uYdW20cu+B/dv7cFAmGP0WpiQ+B97oig43pJZ4doe9TsmItCi23W3bGPZkoQwqMkeCldKHYtdD5HZu8r70xBw9iJV2by+avw9/7RA01TGLMBt74Z8P7zjVzy03+Ifb9/IK06dzX+95Rwtzl4Mrkg7953w8BdoXnMLP7v5Yv7uqtO4Z8Mh3vDNxyvH2Gs0minFNRL7RsaRs+IaM5NQXjvpetAmIDpcQ6Ro2QwXTAqlCpM/B56G/oP+6x33w/4nK1ddKwzB7lVgTnLeTq4fCkP0ZYuH7fUL5osEtzEegzq0yGjLOz9/k1YkpH0TNSPq+Ae9BpWG4FYBFnaJTOtTscuYtqSneSW2kZmgPhtFoDVG+5vG2wGdQ/lQYZbDLZQQZdTedKN689R6ezoGWbO/j+LOh9X1C8qD5vaJMwuQGz1vbnfXMN3DhVFzpA6bojMxEtdsPIB0zofrm9nXM4LwctCgFGirMdohG8yX2NDSH3tvVLy2g+GuB55RlU7HywRDZYsVPGhBERfy7oEvAkepXOpeR4Lw88ITblJStGwlgIUxMQ8aYS+gEIJMMoEQYkJaVgs0TTxSwj1/o2Z9r/sO68zFfOL/XuAlCxv57rsuIJXQl86LxkjAm74N57wDHv4PjMe/wsdftZxff+il5Es2f/mdJ/nZk/t0yKNGcyzjzuYWBsc0/ECJr7UH+mJDe8bymMQRNGRWbW5ly59+Ch3hTgOFkX5a9+9AuiXTS47BM3gIDj4bv+GOLcoDMTT9yueXArPdQQ/aeA6ba7jVZZKAPebzN/bzkR4lcify7O7axsy+DUDYuKvkoXOXqc51YAxHPI2GmqV3r5Pg/8dHzLLudRztr1Zhovap3T2jFmYJHrdoAYYNLf38aUt5ifSCaY3hVQ2MJVoUxfm+kqkmHMzg59Ly+8R174Bdf4bhsZsKH5EKnq7IHKOoi23blFL12Ok6CqbFcMH0QhzV8VSDa+56BnugpeJ2NhwcYG/3CH3ZmImYSvsX9IINHBx1nBXXbVru/D36Qawc4ui/LpvIGYcHzX0uGkKERL97zwhpIxBKVIlESPCORddw3ssXBv8WEUzskaCtbE08z/0Q1v4CLv9bOk56HR/6+fM012X48XsvojYzfhetZgyMBFz3HTjnelj1H7D6K1y8bBb3/vXlXLa8ic/fuZmP/HINA3EPT41GMyVMyC4LhgPt+vOYi+/uGuZAb5a93SNln7k/7hP5/qARKaSpDJDe3aFlOp6/k7aBPL1ZxxsWbGqcH39S+1gcrUIlQWMqmIM2nvBFd5mEIRDSHvNYb2wd4Ik1ayEfyFdreQ4GWpSX8TAwx+FBCxquQ3mT5/f38vz+XuXVPfUqtR0vocaYYBhazHlyr+Nk1BMWL9Dc8uWzatM01WbKjn3wuox6xfZ2j5AtmmUV74YjOYZlxntALIoKAsc1ekMRa1L6RVWGHWFYysWuH+SITJ66+zDWtqVFIdNMtulsf1Xb7YPmHF9pkyl0Iw6EPaxBQZxMOIVAYryBFe+XuHDX8XrGbEvZPa7HMnKeomJ9Y+vAmMe5KhXxlLnjG02guUJMRFoS+Etw9qJGEkJgi+SEJsV6hp3naEAEut81kZBoLdA05bSugfs/AyuuIn/5p/nQz59nuGDyw/espLlucsIUNAGMBFz3XUek/Ts8/J/Mqknxo/dcxD9eczp/3trJNbc8NmpPEI1Gc/SYkGE2wea3o9ln7vdO5OttKb2IB9do7c0W6QkUUXADIjyDODjmsYpATCDU/ZHtY3slJoODgcqBwdn1cZXNd+zFpGEgpCw711LKMuGQaXtONYh2STlNoksTq2DohqtFZ+sbBraTyYcrE7qlyNMJEdrH3d3D3vdblus7E0xI1sct6hq9QfE+CjNrlDC7fMVsUglRduyDBm8l47cvGw6fNSLXWrl3Jfg9Ee+nc2INQ5AsjXi9tjxcgeYW8zHGnog+IsEt3izM6JMZ0raRIoFI+OIkWGbftqXvUYuMM+hhdsVNLqb8e2WBFjO28T7npA0iUdFTeNeGQ6G/8yWLjsHyvMDg0MrG6YrFUZ5dwUmB4O642xLSJp0wMAy1nYnkoLleezfk1PA8aDrEUfNiyPXBbe+F+vnIN3+Pz9y+mfUtA3zt+vM4fV7DVI/u+MUVaefdCI/+Fzz0LxgCPvSKU/jtR16KEPD27z3F9x7dPeF+HBqNZnIZ1x1olZRXxZ6gQHOM9LiZVvfWn4hAlCGBpgyXXR3DPO4WUbBt73MvHDA05koCbOLPoTgjcDJxjc2gNyAY4jiew+aFPhkAsmydLW2DPLC5ffSQp7TqeUSx3As6Gm4RmN5AbqO0beqHdtHcHQ41dQ3GdDJsxlUlE2XLIATSMYRLll0xbMwnLgfNUmI8mtdTQaBLpHMMic29CRrVlUqPR7060VLn0QISfXmTfY7nWUgZOUfOeQXmdjzib8srgRhpS1AhfynYg3BCBWXGu6wrWEbtr6DEpxQGRkCEuB40W0ps6Vd1tCP3cFwRmjiRHDeCTNKIL/QxziIyOzsG2NdX8I9vwPNWybYR/fvDTdMJPx/LVhtLoJXyWAXlIZUyLNYedgqnZJzhJQxjwjlo7rKGre5jd2LBEDrEUXO42Dbc/lGVU/C2n/L95/q4fW0rf/PaU7nqrHljr695cRgJVThk5fvgia8rL6aUnL94Jvd88nL+4qy5fOm+bdz80+dCs98ajeboMp6evx0bV7F21W+QEyyiMaoHzTFK+rKlcYs0W/pGvIgzrPY/Xt7wOORBG8NDNtZMv1lgQct9ZPLdoy43Gbj7GTwy4ym4EcQ1WJOGQSbfXWakHuxVht2oE2Wu92WC3tOUE24W9BzZXpikUJU1neup5HnQfDMunTBCp8s3PAXP72jhYGcv925s47GdY/Thi734Ip4Pj/jrw5b+ZIMRE9oV/Ir2Qb+ATtAQHi6Er9eoIHp0Z1cobHZz66AXpitkvBB1y5xbUYGWiHjMYq7rnR1DPLjFz4Gb0BRF3DEtDKu+bTGevlGrHUobKVXxirAHze2D5ghkt+x+xNQ3Y7yXpiXLWhNEh9xQ5YjY2BBH/1rf0zVcsT0MKiAAACAASURBVCFzW98I7UPFwIMucH/GrgGptjVwaG3FsZWJp0oNx/c8qv7bdjd1e+/3vjNOaJ+7UDkk1OUSFfvjQzhC0ctB0yGOmsPmqW/Cjvvgqi/w8PBJfOn+bVx79nw+8erlUz2yEwfDgGu/Cpd8FJ75Ltzz/4Ft01id4ts3XMC/X/cSntrTw+u+8RiP7zzyBo9GoylnPDPnbW0HsWxJPjuxPCTX3I0tpOfWGzEtNrSMr8qrCnEUVGXbqBvaW77AcKdnMli2jBEVFQSaO8BRck/6s0X6ujsQ2NQP7a643KjbH4OiaXPHulbaBnKeaLVtSa5osaNjKHSuRq0A6ODaeilrmFl968rMKdfoH/UacC2yCZamd7c4a88dfkNlp6qfmahWlTX3PKz+tpRnNCH8cdRVJcPeBOn/I/K9dDz7OyC+X5xLrmgxko+ZAJS2Emci6kGrsC9S+qFdQpRFxQXPxZKmWu91MJTM9aDt6hzmke2dZR40KSWH+v1cMUHwWrTD5zuSD2TZUp1LVyBEPWgxAi0aajexEMeYhQ8+Dd07VQGh6PdW2Hi+ZIG0HA9agoQR8KBJde9atvL8poU6fts7w55cKyA23Ot5T7dqTRAMc45OLqWThrq+4u55571c0WJj6wDP7O2NHb+QFlIY5E1XoJV7JE+bVx9aJ07UBN8J3ot7uobZ39ETXqqYVaGrI11eE3NX1LnhoGXjdLxfCaFyUSdaxRFAUB6Cqz1omomz/yl46F/hzOvYtfQGPvl/azlzfgP//TZdTv+oIwRc/UW47FPw/I9Vk3DbQgjBuy9dwh8/dhn1VUlu/NEz/PvdWyrOVGk0miPDeARaVVrl62aH4g2VSrhGZdx3hEpL946d31Q0lWGRNAyael+gJufmd5QbuuB4hSba72oUgfboji6e3Tex/XeJ9fbF4PaM3N054okAW8Ka/X1sbRtkKG96v2FSStYf7Oe5UcbkHot00TGaK4Tm+SFyMR7E8XhBYr878G/nVvXaqTroekfcwiOmbZM0BEnDb1icECJkbHoCJfAb7pbyR0pVbTLCg1va2dQSc3ykVAJtnEVCJH4TYEOUX8+V7qGgkHOF5OZDAwzkSqHiKS5ukQvblp7HAlSIoxUT4ujS2p9jfUs/+YLjvUuMLdCiY55QLmrsdRJzDXnXTvnyvSNFHtjcTkvvsLM3AiPoQXO8WCXLVhMzOEVDECFvYlDomhHPkLfcljtI77gr9Fkm6YT6jSPEMc4ukVKqwjsiwWo3xFr6zdbd41uxSnggzDF4LoKvN7YO0Nrpbtt5c/u9sNXfFyml99yQUsaKr3SLKqwihCAhmFgfNAc35zd4G05E5mmBpoGRbvjdzTBzCQOv/Rof/MUaMimD79+0kpq0rtg4JQgBr/kXuOLTsO6XcPuHwVIPwDMXNHD3X13OTS9dwo8e38t1336C7e2HVy1Mo9FMDCklXeMIMU4mnVLnIxPrZ+gaTGPN2I4lEqWU3LdJlcCvSYe9HoZdJFUcCCyr/s2WLMyWF8IbGmuCLmKYFU2bJ3d1h3J1DodKVfhGHYpXREWGCiGkHFeOLVW/qKDXJcoLB1Q+VNoaCW0zSr5kOTljstygdMc+XrEbGDcy7C+QbqEBu8CermE2tWeRbo6RECTd3CpD/R08au4lVJP2xYdbyt/u2KK8cUMdsP0+aN8EQGP/FhJ2oTyUUarqeyUpsE+6NPhBhV2KeNBk9PPwsi6uqGyoSlEwLZ7c7UeKuH3HXn36HO8999gXLdsXsQAy4EHL9sJwF9mi6VcqdbC33q1eRIuCtK4p26cygVa2xCjEXUeepzU87tC/AVxR0TOU8zxowVy5hKnEphI7kBZ+DmKwAmbIqxw5McG/oyNWHjQZPynjeN79EO3y/VXnyEQKQ40dYKidnR1D3L3hkBeuGi0G4w0kEOYYKhISOFSpYj/pYp8SiH17VRhphPbBPCW3BYWMb5cQHMNEc8fcNcsEGhPLW9QC7UTHtuD3H4BsL+Zbfson/rCTlr4s37vxQhbOqJ7q0Z3YCAGv+gxc+c+w8bdKRDshB9XpBP/2ppfw4/eupHu4wBu+9Tg/fnyvLiCi0Rxh9vdkyydEunZAbyR80DWcrIk1qPY9NDa2LVl7oM/LOa344962AfK+4Fq/eQv3rPFnmzNJgzn1VaFVqvJ+RcWgJOg+FG2QW0GgueLD/bc4Arl+WvqydA0X2NmpjlGZJ8y2xmXtBNezbclgPhB6OXjIFzWBddzXKv/Jx23wPZbHw7alVxExMQ7vVzppcNLMqsoG5Xhn3V1DTq0Ufo4HKjT2ZosMWUnHAyIRApJOWJ+BYHjpa2md/xf+Zp39XT43HDIGUMo6hTn2PabOXdc2KOWpG3au46hgcUIc793YxuOdGWhYqN7P9auG5tFdkr6x7ubxhD8PekD8912B0FCtRGXXkD8Z4oqM6mhZdVQ7BSFNbKHWEwTC0navgrZ1DOTKcwK9UxcVpLZVFu4bvXomVCQkVs45X+7kTPVniwx6AjIu7M5mdueTNBx4CFuCLRIEzXhDlhC2iWmWEMUhUrjjF2RL5aGEUF4JM9SDz3nphqB6OZ5mzASVm+8m3e8oX6RgqnMkRYpCZqZ6M9tDizNh4oaQGgJeekoTwi4i7GI4VHXYyZ8MCrRABECyNBLej8FwVUh3HP668ZNhwTvawJ7guXZxn1HOJIpuVK2ZEKu/ombRrvlvvrQ2xWM7u/mP617CyqWzpnpkGpfL/wau+k/YeifcdpNfBhh49elzuf9Tr+Dy5c38291beM9PnqVzcGIGoUajGT+xlQjbN5TNuEsnp2Wiv+uuYWHZKgzrQG+WfT0jsdsqmjaU8qq57t7HcBcq7X2c2Z1PeMsJIcpklh3IuQlu11vypEsgXVd5oK6AcQXajvth10Pex24fNxGdbd98O+x9tPJ2ge7hgifQpJSsa+nn4W2dSjwNdypB4OZoAbXD++kdGvLCqiQyFJrvhcGNci4sW1IMGKsJw5VLlVcaypsk8EOzGOlWnijX+9G/P3YGP0THZnVMHIS0Yj1oLnYi41fpE4K04S8tMrVYIhEIWXX2JSZkrFiK8e4FiisU7ajotD0R05ctwkkX+yLONYJtG1rWQCmnRLIX4ljuQQv/HRQF6nVjdXn0jutBSxj+2IITGsK2sJ02AEJKTFtSNG22tA0yXDBJGuWTDWX3lGXTPphX44h4QMtbLpRtrjKxHjTnvDhC8NEdXTyzp0IBFymp3XUn6WIfSp5DvnpOqEgIQE22ldS2O2hqfZj6gR3e+0HRHzz20edZMOTRFRanzK7lTect9I6fbRYoq3LpTB6PVhWyULIxbBPbSCKNNHbtXCiOeJVHNx9Sk0yGEMypr2LBoT+x4NCfwtfK3kehb39oUsktqlNafxsz+zeExhFHNuDdbxvI8cze8lDf4JyLiKnmWpHt9/vrxYQ4TsTvqgXaicyeR+CRL8K57+T38tX88PG9vPdlS7n+osVTPTJNlJd+HK75iiri8usbQk00m+sy/PA9K/mP617Cc/t6uerrq7l/U/soG9NoNIeLW2kvlqBXIKFydUb7OX5yV3eoZxeEQxxdwZB0CgFEjY6QcRUJq0vY/kROjF3qCaCiaYeq6HkkUlA9Iz7fpH2TEiPgG/XO2ESkrYAhY4TAyOgFjvb1jHg9hKT0G7+alvS9Gvl+JQiyfczo38TM3g3e+kq8+NvzDMsKVpaUkrs3HGLT1s0sbLkHwyqScM7cQHb0SoxCSv8k73lEeaKCVeQGDo66vptr5u6rYVthY9CKHn9VcENKQEAqICQTIixEvZYBZeGKMr5peMBj9EKLysFrG8ip5tfSDlutRgKqZ4bXb1+vwsra1qt8o8BHo+VvxXnQ6qsiOWEogZYwBEIIZtdnQutaUiKkhS1cYac8kZ271pAtmnQO5b0KjpFDoRAG+YWXsqFFeYH7s6UyD1rUIXrYHrRSTnnc3eMZ+p4KXldpByoIOcdPJDAiRVvqB3dStGxyJcsTssKp+uhtapRxm7bNof4ce7qGvSH7jZad66tULO+H51zzvhc78B19+6E4QsG0ELaJdM6Rla6F4og3gVKJwVwp7EHPdof2p3ekyI6OISxLBkSRexOo54hp2ezpHqZvpEi2OHbocVSgjetcSwnFYe+a9CaZAtvUHjTN2AweUqGNs09n/Tmf4zN/3MTLTmnis9eeMdUj01Ti4g+qMvy7/gy3vj3UY0cIwY2XLuHuv7qchTOr+cgv1/APv9tQ1kdGo9G8OIJV08oINCYuMxRi6B4pljXjdfNM4qqLlVWTtqUvzFxREFPavSwED6egQPumsmeENzMtEspDEpdH1bXNf10K53PV7r6HhoHt3t+ukZIpdJMu9JVvK/TlviFkBHo6uXa17RaqALWf2+8htW8VEBaktpQEfYauwI0WRPCXV/8W29W4k+YwhnMcWkfJVztjfkNsCfnQMYlWBxwFico1y5sWrf3qWop60NzqhK4ASjmHY6B2WVkJeW9UkfMvpEmpTPgBln8MpdP0ubU/x+qdXSEPmr+hwHZtO1LEIdj/aXQPWvAWccPZkobg6peUt/dxz+U5i2Y4X+sIVMtGYHueayFtTFsysF8J92gepr+f/r6UZCDvyKDs2o8a6esO9rPlkF+BcVfnMC19WXZ0DLF6R8QTFlx372PK4+7eq4FJDS/30rZgKDDRKu1Aj8SAVybwPJIYoftACKcvnrRC1+hoIkFm+3huXy+92WJIWEBAH1oFSESKxVh+0Y2y/W55DnY/TMFUOWi243mVySqwTcxSOGQyKqQH8yV2dATCyku5somvrYf6aQtMNHnDcI5ta3+O3pFiKH+4sbrCvRl5figvefyiIWxT5QZ65VPVvw1Van8FkftglAJLoAXaiYlVgt/eDMUs3dd8nw/+eitzGzJ8+4YLKlfP0UwPLrgJ3vy/sO9x+OVbVCPcAMvn1PGHj17Gx155CretOci1tzzGCwfGMIo0Gs24Cc/8S+5Y10r7oGOQBz0nXoje6NsKejIsW3phXKbtGwXJkXZVPCKysbKKav0HiWuMnYgYPLaRUQItKLQcWvpy5Eqm8pAIY0wjAtOtgufPqNcO7wecvlyB8c3uKs9VCuEY1aBExKyeF7B7dvu1FGx8o9kqgVnwjlG62E/SqbwYPU6u13NrW/h56RIXjhUsX19JZJ86tx4Du/wcm0GBFi8M4pBSMscJTW0bUGF2MiKkhCNibQlVI61kRtowEzUMNJ4R8KDJ0L+GMFjSVONtY8GhB7EH2soHEAiflyIR9rLZVoxAC/xdDOZlCtwcOVBipyw8kPB95L12wzINQSZZfuzccxltR2E5AsEy/BDHzsG852WWTlxgKdUQ2p4/CSKwZTQPLTp5EcaypZdrubtrmM2HBrzqodGJl9Dabln9aB6nM26PfY/7Oyjt0D67x08EzoGV9M8xKAP/1Hn1COyQMKjkDUoV+0nve5j6wV3O90iQssyDhllUHvZlVwQORiG6l9641Tp56NgKSM+DZhsqL9aOCLQx64Y7Qii0jjS9Y26LJJYtVeVF59h6uaVG+cRNFCFNf1IrmUF4QaVjYBW9Y1udSnDpshlcddY8muqUmFUeNH87Bx//1aib09b4icifPg8Hn6b4+m/w/nuHGC6Y/OCmlcysTY+9rmbqOfd6eOuP1azUL96skrQDpJMGf3/16fz6g5dSsiRv+95TfO1PO2JLFGs0mokRNG5c46+t3xEpVnAm3M+Hit2OYzEF856CCfu2U3a6KtdBbdvT0LVdeViEYFmzStqX0YpqPTvDs+4OVamE50GwRRLbSGHYTrPrGGtoMGcq49tIlHvQooLN9RYFyq9Lx2i88oy5iMj6pm1j2rZnMIWQNpZjoRu2SXWuDQ6t9QzDoMHlGn7B8zG38zHmtj9Kou0FekZ8o6++Zz0LW+4JhUmFc3LC52h211Okh1q8vwtx4YDutlseKT/DoSIK449pii6ZL9nI6PGXNiXTRiKZ0bOWqoTybM1vrPaMT8+D5oadArWZcE7XoYGcV3zGIxiOKQxfoEkZ70ELHrdgrp3jMfMbVful8ONWjfOgxYUjgl/wxd2me3+57QikkWJGdQoIh+7aUmKjKh+6Ho3g+giDUiDvzrYpLxJSQdjkihabWkev1rp6RyfrD0YmS93tW6VAWfrItSZtVYSoMBzIh1Q5UYYQkedLRLQ4uaciMrlTyRuUsArYUpIu9nnLiYAH1zs6tqkmZIKTD26IY/RWCTwvUj1bAFVFFsByJnXsiBfeECK2wE7bgLOcMMr2wQh8jxQJ9vaMsO5gP1im84x17434522QVGlQLb7gfBVGOl4PmlXylmuuy9CQSVAVKGgjhIjkl46+US3QTjQ23w5Pfxt58Yf59LYVrD/Yz1fffh6nz2sYe13N9OGsN8Pbfw5t6+Hnb1RlhCNccnIT933qct507gK+8eedvPV7T3mJ+xqN5vCI+0233ZwQ1zC3bRJF31vTmlpMy0h0HfXj3DNc5E9bOiiYVqiCmmWr/jxJM6tMto5NpPt2U1eVZMksJdAs13B2KeWgfWPZ+GrSCZrqMjRWp5BGyhdoFfbRtB1j3Ej62x/pgZ0Phb2EoD63ba9ghADPkE8Ywit9biZVwZF1B/vZ2jYUnycrbe/4zuxb773tGoaWHRCkFTyUSXOY2pFw3ldm6IAaj+mfBCtksJYfieDkuismy6rkmkVVNnw0O2sclRyLls329sEyg3GkaHr763qTBDaP7/JzcFIJgwuXzuK8k2b4IY6eB83ZFyExYpT43p6RsJE4HKjsKYzAeJQnpaxJddB7awbyGIXhV3HcvYqGXX+kOnsotndVwhBh8eAMOhETlgt+ywTPq+p50JSQuGDZXE6eXeeFCqYMA+E0CB7KmdgiSToQKeQ70ERo0mB39zAycq1XOs87O4fKvNTBfQEYypW80u4eXv6mzSPb1bEvay+x+XZVhGjvo4FqrngCLUh0XSGEV5giLI4rTBoZSQZyJa+fmi1VTpfvQXPWtwIedpeRbsj2htpb9Lq5i5HjUUqpZ4Ep1KSOMRT25hqC8ucMKkzRa6XhPL1Om1fPyqWzQkJSBsKK+4aH2dDS703OlAJhzrYVmPyQkjntj1KVa2d219NqDMJw9jFQJMQqwaF18ZEFtv8sMAzKLhjblnQM5r2847GqbmuBdiLRtQPu+AQsupgf1byPP6xt5f+95tTYOG/NMcDp18I7/w86t8FPX++Xnw3QUJXiq9efxzffeT57uoa59pbH+PWzBybWYFOj0XjEhcNJJGsO9JHNOTO8vXtC5aa3d4zQ0h/2VgQLOWSLJu0Dea+xbiZpOHlGzgy2IwBLrWspmTbGUAs1IwfVNkJNbuPva7cseW06SSHT5Ai0Ii8c6POaAQcxLekYYI5Bnh9U1X7z/ZAtr3iGDJfOd3scGQLm1iWxElVkaxZ6n0e9Z9miyR3rWukazIVmuNUuSapGWpnb9jDJnff5IZzujH1on+ONeteom9f+CHWOJ8kKedDK13G9doXMLPKmWsCKHt9cnzL+46Ruxpn0HEe5/q6hPEMx+cL5kuV5DP3eSuq7hvKmZ8ClkknSScMTNdvahtjQ0u/tt0BW9Ejl4jyZauD++tLGti1kWTuBSteecHIHBWR7EQhm9a4NHT/3ZcL1KhSGyeZyXki+K3hcceD+HfSgGVZRGdk7HiRzUFUxzVTXYAhBdUKNL5kQGIYqhtObLSKNRCiP1B+SYH1LOBrFNAMetP1PMe/A3WVHKVXsp3fXc7HPhfD1MspvrlXyvbRSNVHe1Tkcao5s2xaD3r0qHW+6+muw4TS6my8pc18JAp7GgCc2OlT3+VCVSpIrWZ7YUQLNxDj0PPTsDm8rKtCAvo33s++Rn3ml7oumHbr+3Zbq2ZqTAHhszwC2lGT6d4UH5Ai0cxbOKDtUXtSAWQJpMae+iqbadLgtR2AiITcSbolSwI8UO2mmCrFcOKMawy6RModp6lHVeBNCqEkHIRwPmnPQOjZBzy5VodUZz/6eESW2rCKWDPRyi5wPVyS+cKCPfNEc0yunBdqJQmEYbns3JDM8fM5/84UHdnPN2fP4q1cvn+qRaV4MK14L77oNevdU9KQBvOHcBTzw/17B+Ytn8Ok/bOSDP19THt6i0WjGpMzbEvj7UI9j4AWrpklACHqaLooYXOEf76Jpe5+nk4bTbFYibBNbSkac0tC5kkWq9Vlm9m1QRkHICCj/xZ9dl/EMuZJI0zfzJZ4HDVTvJYBS0i+pb7neEreMeqBBbFwDXdXbzH9fejPugtPmVHPekjkkEn5YmW2omfPn9/UipfSqNLb0jpR5KKWEGT1rSVpZZDELXdv97yRsbJqJsXt3uqH8ozXrBd+wLWSaKZSUAA2uI4SAwqDnVXCLengk06Fxlu1UAPdSkBGD17Skl6/kaSPnOBdMK6DaHNHirN42kGNv94hybDql4yrl9VTqB56wioFj6wi0gMm4p2s47FELCtF8v3vZh8Zll3zBI4sjZPLdGIbTG2rH/XSvvcv7POq1cUux+8IN5rf9iZq9D0JhEFlSHrxESh330+fWIqSlhBzCC/GXIhnyjnoCXwhKqXoKmSbvM7sY8AoOtsYa1LM7n1K942KEePi6GsUaD4gnQ5ps7xiiP1dUlSTdRUJVVdS4hdNXa6hhOYWqZkT0O0I1XAL3Z+T6O2fRDBqqU8yoTqhrxi3QY6s2GUb//nCTaNtS5z5yvbp95jIFVaU16uG3bchkqlgxV01eWMlqp3BP+OpUIY6m13ctiJSAbZHZeRfz2h9RIZzCH7Parp+L19oZrhibF34/yFNn1/DKU+ewcuksalPh7xJCQP08EAaGCNz7TjsB99pv6cux7mA/u7qGVYijm0MpxKiTM6Y9dm81LdBOBKSEu/4aunew64pb+OhdbZx30gy++vbzKs6qaY4hTn4l3PAbVUHrV2+FwlDsYvMbq/nF+y7hn649g9U7u7jq64/x8LbO2GU1Gk08+bI+aL4BcrB7UBkp0gYkhUyzM9ctMFO1rDvY73mIokbSlrZBP/8GyZz9d9G9b6OqeibDuRNeTk80B83yk+RdFsyo9paZueQsEAlsI+WVw5fASO3iUEU+A9QMuWvJjgS887GhPSUoDkON6p8pA98vbItMOsUZCxv9xR3h19qfC+d3yVJsMQl3ZGp/K1fWs5K+8eUWOmiqzbBwhi/catNJkBZiy52qqArxIV/BX0bLqWsf1NQpQ0BhyDtsbQN5FdLlbcBwqgLEHK8dD8CWOz2h5n57vmpOaDEzPwSlnBceCpC0sgGPpVteL5yX5dI5lFfHR/qVMKNUMhFLqXpa+pRHOFPoQ9qm5xkF2Ng6wKP78/51GVR6uT5S+V5foLn5P4Eqp1W7H6C5+xkSQvjhjoGm7lEPmpvL4ze/dgtW5MkWTa/aZiKlroGqhPI21mWSCOH3F5Qi7EHz8rqEASJB9+xLKWSa1XgDlZJtKWM9pUFRtGJOuCG4GRJVzvfHGeWBazpZiu+bF52GkTHntKdpJdmaRd7fhggcp8B1GB3B3IYMrzptDo1VSXW/BSuoykBerXcjmmryJloAx91HN2fUskPPC1vaJGWJeY3+fbolPwMrUhFSCbR4cWM7+yKlul6Es49BD1pwsika9llMB7xy0qaxxu1XGfm++nlOGwkRbqXhnivnGeY+v4qmDWbBE+WG4XjQsr2hitsulmXrEEcN8OwPYNPv6L/073jHnzLMrs/wg5tWhpIXNcc4J18Bb/uJio3+9Q2qeW0MhiH4wOUnc+cnLqO5Ls3NP32Of/rjRnLFscNwNJoTnYJplZVdDxoAhrTY3j4EdgkbAzNZw1DeDHlG3BSIqAcN8IoymEXl3a4f2q22L6VXIv7yhg6/fVK0iqNjIFkJxwCStlrWEW4N1U7VNCPt9ScrWTZSGCQD+Vk2ODPkMb8RroFSPQPmn6deH3hGGSOpGqy6BaF8ENeYSwSMORExhtJta5Q3xcyXNYYO1jGJq5If8mpJ28vVahhUnraXr2gOVQOsSSeUd8gsqkJLgFUoN6A8h5VbkMA2PQHdUJ3i5SuawSqFynHv6R7mYF+WkYJJX9ZUx882GcqXWLO/zzfIisPqnPTuCe2DFMnQ9qy86kVVyMwKNd6e2bcJsH0D3bm+oo2YXY8GyDLxFnf8AEjV0Db/tQzVn+IJkqae55GlfJlAGao/mXwpUBY+sM1UaSBUZh9gd6sv9N0tJYyAQDOE14oh6CkLjrPeCVH1dIdUkxseTuGJ+jRcumwmC2dWYwjhiSVbJAi23XJ3P7hn3bMvoZiegV3w73WvzHsFl6OQNg2R5tpW4IJ1hVxsz8HAsUuX/GIjoUqXAWPetCU2bpEQn2JmJisufFXgO4X3b6gojmVSN7TX23n32nIrl7r3pyUlhlUMV3GUqkBNKATa3W7gm91xhjxoUp3PYH7hQN4KV650167kfXI8aNJbVjgCzX/mlNLh0EgzWeuPwclPO3VOPcGKtzJ6Xt08NiEwglUwI1Vy3euyd6QIZs5vE+GGOO5eBdvvK9sNy7IqdbzzhzDG55pjnT2PwP2fpnTKVbx1o6rq95P3XkxzXWbMVTXHGKdfC9d9F/auht+9D6zynAZv0XkN/PHjl/HBy5fxy6cPcO03H2NDJP5eo9GECZYdf/XpytshQgLBVLPGVtHzEuVKVjh0zRVoMYae2bWLhS33BnK01MywjT8b3zC4K1wgIWY7tiPQElZBGWmDreqDutnq87LeXH5eyEjtYvWWYcQbSW7V2HnnQNoJJcr7oZ22SGAEE/xtC4wkRsBrETxmz+3tpv3Adpq7n8Ewc57B7JaFl4FYuWh+GqiwKSkM8lVzWT4r41Twq0xDVQqQauZb2pQG2ml56vdUZw95y6yYUxcwfFW4KVvv9IyxU+fWq0bKMS0NOgbzbG0fZMuhfkwMfjhKXAAAIABJREFUsG2e399HS1+WobzKtfPCIQ+tpWjZXohnvqo5FH7X1dOLLSWlVAPJgKhIlQacqnUBlQIhERfCCYcLehb8jyLhlgjsRJpSupGBxjP99+3ykv9WsoaB+uUM5Ipl10rCyntS092nzl7nOrFMT3wYhmAgV3L6uglmdz1JVc4vIOPu08mza1kxp55TZqt9EBCf/2ckPM/l7LqUMt6FXxxCGsnQcfJyRR1zeEaNEniWkWHj/nbPY+62v1jYeh8LW+4hXXBzMf1YZrdNUSbfhWEVQiHNGWd5t9DFnu5htc1E2pv0SJgjGHYRM6Gu/f09vscxKKTzJRtkuJANwKzaNLMDtp3Xv0ygCns4pHp30DiwhZpsuJiOe40JRzpYtiQtfUFpCD8nFiMZU9Uz/GewT6MtVdEjIYxIRKPBsqaa0CSKG+IYpW/m2Vgzlqp+k0HHJ76ozFXNJV81O7ReMT2TfNVcepouxP3yVFLAYJsq+mHb5cV83JBsI4mQJonBfWpZd/JbuiHWaiB92SLZ7IgnhD0PWgUsy9IetBOant1w23uwm1fw4ZEPs78vz/+++0KWzyl/SGuOE869Hq75Cmy/B+74+KgVxKpSCT577Znc+oFLyBUt/vI7T/KtVTtj8zE0Gk34dqqvSpE0wn2+hG2qSm2WiRTxQsHrURVzbxYOrgckFyzwcyiEVLO3JUs1qk0mjMohjg5mUoX0zWtfRfXeB71kd6oaefnyZk5f0BRaXooEdctW0rrwavpnvISeJdeqD+ImeZzk+FCOmrdzFmZ1Mwm7QLrQq0r+5/uVBy2QgxY05kd6271qgUYpi+00C84kEvTNPNsRUuqYlWJcaCXLZrDhVKQwMKQZqiIXR11VEiEl3UN5+rJFDrW3U7Jtr7S4QNBYnfbN7kAIXlmFQatUMXQwU+hWIYHSouQIe1dMtA0oI29b+yAbWvq9961EVajaYik3yGDBopSq84pjqDEpo9HPSyvv8RRGbf+kU84sm5wdKVhec3TwhQr4PcXcsdvzzi7b8raOYXZ2DlPo9dsS2FJi2KYnhFxPTsJSHqmRXM6b7CiW1DW2v2fELwRi+sLEvdarUwnOXNDgp2U43uGy0+16fgN5kSNN5/jl+CNen4N96rvc3MdFM9W9I40khrQ42BlTFAeoH9oTGp+QjkCTkubuZ5nd9TRrD/Sra6aYpXFgq7duybLpHSmyq3MIGhZ4Qr86p1IOipmZoWNp2XZIrEkkedN2qlP698Q5C2eEPKXhPD7/QBWLBZWbF3l2BHv/gZqQqrH97xUIhJtf67bhCOB5tZzrsXOwwMFeFbL5woE+ciWLRKQnnhQGCUzmH7ibuqHd7hdVeK7VYjui0PW0C5QYcvelb9a5oXBtUOG6Pc0ryVfP8ybLBCqHlC13wN5HaaqNPMvcZ1u6DrK9pNpeoG3L48r7Dera6t3Lvn17vFXsYi7cJsIqn8BxsexorEA5WqAdr+T64dbrkcLgX+v+mVX78nz5redw6clNY6+rOba5+IPw6s/Bhl/DfX9flpAe5WXLm7n/r1/B686ez1ce3MH1//uUVwZWo9H4uMb/uYuCITThEEcBIQ+aQgSWdgzFGCExUFBGRlMVzG+sUtlrsqSqPA7mSSbCIWNWJITIxTL8HA/DHFGGgmOYNtVlWNTcyKKZgVnrZIYLlzY5xq3wjdjGhdFN+xhGeQiktDGrVR7a7K6nVKNdACMZqgCYCByz5u5nvNfCNrGcbdZkEuRrT2IgMctb2owRaEXLRooEUiQwZAnTMfxn1qS5eNmssuUThqChKsFgwWR31zBZR4MKaTG/sYpzFjVG1vDDw/z8Eucjq1TZawXYKKHgej+zkVDy4UjlRikSZALFClLmMFlRQzE9C3PGyV6RDikSCGkFcoLUPlcqTU+NyqlaNH8BS5tqWdZU63l7Dg3k2N4xBHVzkFKGBJoMXMPSBpksL8Li5kzJwrDXC0/tro3hCDJXVDUM7gCrxCNbD3ki1baUkOoeKQbukoD/0nmzTHw6xXM6hiIhg0bC6d8XKFwT6tGXiM0Dc68xr1KhMEhYOer2PlAxHwoIjbkuk/TESdIcZqRo8vz+vrJedqE2CMkq1ahZSqpzbZSSdaE8KcuWdA4VVMuFANlEfUiHnHfSDBprUqH3DAMQSvK7z5t8yaJjUOUmevdkx2bYvapssqFo2TTmfeGtinE44YVGskwIudezewx6Rgpsbu1Xje8D22isTnmtDiTC89y5eaOGoMyDNrL41RQzTUi3SIfrQXPOsRvi6Oa/9i98pXc/uM+zTNLAPWNCoDxxANkezlnYwFkLAu2m3CiDTJ0XxtvXH+hjZ1uUDjxHc/ezgfdK3jMiaQjo20clLGmPORmuBdrxiGXC725G9u3jhwv+jZ9thc+87nTefP6isdfVHB9c/jfwsk/Ccz+Ah78w5uKNNSm++c7z+cY7zmN7+xBXf301v33+oC7Hr9EEcH9QG50wOiGiIY6Wyi0a7sAKVRQM3Eduyk6M0efO7iZl0SsHLwLLRXOMdnYOxwu0QLEMA6H6owVnuxNp5jVUUes0rzYSYW+fd9+na6GqvNS1WilZHmMlrZjwScBIkEoEQpgqzB13DWVp6S/QO+t8EvVzqE6nKJD08j/MGFFbcqrASWEgpO31mprfWMX8Rv8cnLNwBmc7hUoWNPqeIbf9gZA2cxqqPOFSm04wv7GKGTUZvz6A88IQQlVGLgxWrI5oJarUbL8sF2gidi2BFMlQQ+lUaZAhUUs6mSCz5ELsRJplTbXYjkDzcEt7V7LoFl6gQuDrVFhuU12GU+f6kTSdcy6D6lkc7MuFwtaC57Jw8mtj5/qKGSWCbVt6gn1v1zBV+S5q9zygxhVcoeU5r4KoGrtFNDYudE95YXrlAi0W17MrrcBxCeyTSJZdfaZtY7tePiG4ZFlTyNNmmeWekGgepZCSdNIoK0rRNpDDioaGWq5Ak5CqZihfxLCLCGliJWuRwmBuQ5W3TNSQ75r9Ug7NvCjkLfPyzULFfoQK0UOoyouoohZGtMdY51bI9pKMeNAsW5IMhiaD7/02Eqw72O9N5uZLlucJrx/ag+EUfKkf2uOJcXAKegjBRe7kiTAc8S290Eq14+FjtmTBfPWiMBh6Xzj7Imwn11cIXrFiNpeffQoNNc6EgSPa0kkjVGE2WFAtIS2qU4FJNfeZmK71C9PkA2kgh9aWPY8sS4kuNy8uRKTSrW3aulH1CYeUcO/fwu5VPLjs7/nC5ll85IpT+PAVp0z1yDRHEyHgtf8GF9wEq/8bnvzWuFZ703kLue9Tl/OShY383e828LFfvUDfSHnDSI3mRMQ1tBOBKgPhEMcSDR2q8IQZCA8L4v5ExxWmkCJBQgiMkU6q02qG2ghsP5UI/+hLKbGyfcogPesvQ9vxB4USaMH3HOPD9WwkU+rvxbNUaGW2ZDHieXdijIiaJsjUx3rQbBGz3+k6ajJpmp0S9wg4txmvkqQ3VKm8YbmaBYiTX0nCEJjJGs+QifOgmbaNFEqgGQJm1yujLJM0QtED6aShPIa2FcrncmfeVUn2wPtCsHBGDdWZlGcge8UshIDuHQAVKyF3zL2ckg29Qzlv/NmiycKWe0jEWF5SCGwjSVXSYOGMalbMqcOwiuTsFEXLZklTLWcvnq3aBDhi1PtmN8QxxoN29sJG9XuQCnu/giLRNtIULEnnUB4pDJY01Trv+wJtmOqyAjkAI3VLsBJVmFJ6kwBDBZOElffCGEOpSlYxVNBhTuuDJE213Thz1RMc0Q8rCTQjEOJou9/vX6fSSPg92JxzZ1rSK0BjGDCjJoXraWntz3HvRj8/0SVT6PH6fQEsn6OO2RnzalnQWO1NsEB57qQvuAQkqyiaNgkrz/z6tNPDS1DjFHFz+yEGKaYbkUaqrEiId2gcQWIIR6AJPIFGKUdNTu3P+cl9nDnPrzwZFWgQfuYIIZzG7BKMFPt7Rnih8UqYd3aouI9hF2kc2OYcp+5QdVNXuLi5otITaL4wl1Ce0+ieq8Zl/jL4968hLc971lidIpNMeNeOO/GVTkSeV2ag1ZDlTNQIQ50713scmOwaLpheLzMofx51D+XoHSnG3t9svp05HY/5X2fbsUWPgmiBdrzx6JdhzU9Yv/R9fHjzWbzjopP4h6tPm+pRaaYCIeD1X4czr4MHPwtrfzmu1RbNrOHWD17Kp193Og9t7eCqr69m9Y7yJtgazYmGO/PterKcn39AGblz6pIk8ipnxUzVl60/0HgGVlq9v23TC1Tl/DYXy+fUqWqKCQP69jkhQOHqZMmYX357uFOVhDYMOPlVWAsuDAk0ASqUJ+hBc0K+XIM+5fTsOn/xTJpqMwzmSjy0tSNUFCVEteNVK8tBs4ntuFU9A4RBjesVBGa1rQ6VVQeoybaGCqokDIFpG0gJtkiFGv/aUnKgd8Qzdt3wvzn1VaxcsViVUq/Qsy1oiCYcw8ywiyRyfWWLG4aBLSWD+RJP7+5y3vMFTybWGgNppNnemWVbm79N14OWiHF1SQwVpikE8xurqa9KIbBD57JmxRXKKBeJsKfGrpyD5npiPGqb/X1wVzfS9OZ8g3j5nDpesrAxJNBeONDH5kN+hcEgVqJK9dUUBjQs8Ixw0wu9NPxKelKGPGhS+gU0XCF7/mLfa+t5kqOHLOBhOTVY3t5tD2Gb3vkXgeO9uKnBqbYoqD3rasxkrapk6oavCkHSEKHjHr1OXZLmkCcU3MmH5c01LJhR/f+3995xdlT3wf5zZub2sr03aaVV74giAQJkQMICYzAGdxNskhgSlziucWL7lwTHSdxel1/cweAYE5uXGGxsmulVxhgwvUiAEKorbS/33nn/ODNzZ+bO3b27Wu3elc7z+Vxp78yZuadMOd/zbWhAvF8G4SgmoGW1MIRiZE0wMv1URIQUINEcc+ZM1mRwxC+M2mkVoDYh72V3BMnu6pVyv4YloFla9Cd+wc4X/+SUq0lE6KrJj3F0z58KtLvua0oI0LLDUmA0XAsxtQsKBGjDCFFx4Cnn++76k5w6g/tZJlzXoin7bKSvQAC3NduZyjlQv9jxKRSjQ7DvRYxMvzNm9vmcFA3hECfNr7XyqglXD7rok/f2MR3VLGpM5zVohtdnc2+vfF6MZnMFJst7ewYYyeYK72/7OZvJp1B4aU+vyoN2VPGHq+DOK9je+hbOfeZNnLWskX89b/mYNvKKIxxNh/O/B/M2wq/+VubeKQFdE/z1KfO44fITqYiFeN+PHuYLv/pzQA4oheLowZlwuictdjQvLUwoOyiDf6RbGEh2FBzfl+pkuH0DPUOjxAZ3UrPvEWdfSNcsDVq+fFt1nDlVEdJROfkKuV781YkwmCbZ0ZH8JCJRw5MDVbhf7U6YfeE1cYS8SZwWjnrLW9z85M7AYCaOYBaKeScwZkDy1XBSCgVCBjfJGAnHBEgLiIKIT0AbjDWS1UL0pToZzeZ4YsdBsrkcPYOj7O61V8C9E2pC+TQDBZhZ4q/dm/8Na/IdGd4fGPBD1zWyOXhxT58z1o6PjBCElr+F3mRn4YHAUNabn2mXFWI9UOkmNCuJcn5yWRELUZWMcPxcy3c8HJf9aeXZczQx1m8Iy3zMjeHTutJ5qiwL9KbmS8FQMxjMyHJ1yQjJiCHzkwWZqwaQEyH29Y/QM5yFjvUkI/K4lior2I2A/kQbQ9EGzOyoZ9HBNPN+RPaVE3f54a1qq2RNe5UVfdP9o1mWNKZZ0VJJ1DLVdYQLn4mjbk3uYyGdxS3VxEIGa+dU01DXSE4Ly7D12byQq2vCs1BgZApzk0UMnXn1eX9FzdaN2z5oukZoVJrQZTJeYcPWxO+tO4FnDmq8cmCEyPBeDHKYwsDIDjrjuG1vPz0+X0X7JtU1QXtNnE1LG53okwBDsUa6q6y5n27IQBqWaWDv6896z+XSIoV1jWTEq2VynnWWL5ueG2bbvn7+vDuvFRvJmo5540Bc+q3GtIxM4G1hurV6FtGQLhelXDdEVffjJHY8UCCg2fXImSakW5yLRX/lHnj9j0SG96Jlhz3H5C0EwtQkIx6/O6H7Fpe6X/Z+t699I+rRYNrtfG5XrxNgJt9Iy6TWfctpBoQSrjLWyQ56I2gGoQS0I4Vnb4abPsru+pM448ULOHF+LV9/x6oxojopjhqMCFx0DbSshV9+AF78fcmHLm2u4Ma/PYm/OHEOV96/jXO+eS9P7gheRVUojnTsiZUtKEkfNCsktR5GIKTZSvMqjybJHTkvJ0I8s7twRd6eFLo1G80pg4aEjqiRJuohQ4NlbwOgKhYmNHqQXGbUo8nqG87nXctqLuHJvaorBFTNRdc0MkYSEa/y7HITFNreUziWP5Z4DTKnq+u30k1OWXluQU/1MjRwzKDceDQXmmDESPBGy5mMhOVkeDiTlW10H4NGTg/nEwNbAqijZYm5/OhyWU9XuLUjBRrKRC26pmNawlB45IBTL7KjoIUIReL0VC5mKNrASLiKgxWL6UvOcdri0XSNEW3RREiTVteurvoUq9trPIl90awAIWbWlZspLwT6zRyNIo5pQkBPuovXW88CYMg6RUdtOl9HayzsxM1BLGupcIKJDFlz6pyZIx0xnEAQcgIuTTjNoYOEXMmYTUyiQzKsvj13dTchpGu0VeejmjrkMsQjBmFDKxR4rfxzeQ2akW+TrUnWI0RC0icpk82RswQCIaSg6xHQHFPG/A8tb6lgfq1rYcNpkPzN9uo4yb6XqczsZSTjXdjM5HKMhKvIhJLs7BkhE0phZAbRyWIKneFIjcyjhQzU4aatKi6TrZO/joJy2g7EWxE186FxhRPFMWeaBT5yHjM/8poqG10T8v4xc9LEMTvCSDbHC/vzAtrNT+6U/rDAgcqlZIwkYXxCqWGNYSwv1G7oqmOpFXnSvpTn1SbRycr7y+UbawtxuRwQq2Sg/bSCNgty+cUM8sKgZpk2yuvLEtBChQFvvA0POf+7Uzn0DmcYGs0yGLBYbfet5/4OxR1hs6s+yfJGS+PZ89zYv48S0I4Mtt0L//MXHKhYwsZXL2HVnHq+9961nrwSiqOccALefR3UdMG174bXtpZ8aDSk8/lzlvKTS47j4OAo533nPr51x/My35NCcZQwms3x9M4edE14hCi3Bk3TrBVeI+p5qbsDP+RM0xPEw8aexHqcy60wzUNanKFoPSFNTuJpXIGmQf3u+6QA5dN0OAKaEcMM2VoM3/sgFEUXMBBv9rwr/CZOgYHG3A789kS2uhNa1lir1G5hUHfKaVa63pwRAw0ntH1Q3e0+yeZyMgeXa/sfB+o8JlXJqMFoKOUIBI6Atu95GQjAvQSeG/WYIMkE3YKWSteELZyERWfDnJM9Qlvt3oeliZ6dp8maxK2bV8O+2rXsqV9PX6qTg5VLnTY7Uf20fACJoGXTg5Uy51iheZRv3IQMdW4nMJedln8W25PD5soYq9uqii7SakJ4tJW2tZZmeP0T32jcyN7atYHnACkw2Jq2rDWlzOZMT2AOTQg658wB5CKHW7MCUnsJ+TyDJS0ruwJd6EIQ0jXmWLnz0AwrQZ6tWbQFRU2OWcMymLuBsK7Rk17AKAaZSIVzLvAGKtFzQ8RCOjubTvfWYTTvk6flo//I39I1GlJR0rsfJmMJaPWWf2Q2l39uCOQ1r2WH0IQU6kciVWTnbQxsdn0qQjrmHaNAhIbZvApi1eT0MPS8EehDiE/r5L9edCEgYeV8zGXRc8PktHCAzanEFAY5LUSU/Pj0JzpACHY1bGCo9SRneyys01jhFb6FQN6/2VEwYlTGbG2/S4MGZLDMGV3V7UvO8Sxm2MfYArp8JtsCWuEzOF8BPO17o3Ej3VUrne+v7B8g5hKKnaBRZBkNpbyCVSjm3J8VsTCdVcG+yUEoAW2288qD8NML6Y01c8auy1kyp5kfX3ysZ0KgUABypfu918tIXte8DXY9Nf4xLjYsqON3H93ApqWN/Octz/G2//9+mcdFoTgKsM3TdCHkxPi1rRgjvS4BLZIPf296ZYKwa5KfM02yeqFGwI5GmIgULqwNpTvZV3MM2UXnyA1Ccybyg6NZz0qzwGUyJjTMiOWf4/cXExpV8TBddTEWugIF+Od8WZcwtrdvmKdeP8jQiJzU3f3cHh595SA7Dw7K54uWF8YcsyVnm7DCfguGCXvyfblxa7Q0TZDNSSHRNpECyOoxjwattTrBScvmk0okoHquDGICsOdZeOF2ORi2KebokKeNenbQEzWPplUw7zRpJqnpGIa337TccF5A0wrNttyYQkfPDtHy2q+ppKfAlA+gLzmXHS1vZjDebJ3LdxL/uYUuo92ZphQJ9BC0He/slqHEIR7Waa8J0DwB1Mwv6P2MFanQCMkJpC2oZI0YIaPQ1LEuFeHkrjpCet4UMmNqHByUocY9ipj6xRCrRsuNFmhk4+F8/xaEzB8LW0Cr7kQIwcquudTMO0Zu02T+OVs7pFmTcdk1AuoXQdTKqRavY1/7WeyzZJf8WLqDXmTI5EzSSV9/7sz7c+UPc6XdkHnKHRPHiKEzHKn29MHBwVFMYRDK9KG5fN9C/oAWzu8I4pZJpz+qqx9dCNA0cuE0g/0HnWeYXRcARrxmeoYuGIrkEz27/S01M4uWHfbkx7M5ULlMRq4VgpwWpjLkMmO1hJ1MKIUe8vp0+a9vIZCLKr07wQjTWZdgZWulc184OQltfzPr+B2tW/ILI3bdrX058ho027JBBCySeSviygdoxBhI5KOg61p+IS2ka5y6sF5GpzRNhqIN7G1xaffCca+5Zqb0oGtKQJvNvLYVrrmAvkgdp+/9O+bOmauEM8XYpBrhfTfIyERXnwf7Xx7/GBdViTDfetcavv2uNbyyf4A3/597+e5dL6rk1oojHlvgWt1eJZObdm+jau9WhJklpGnUVaWtl7aQkQVdt4R7DpLLQVbPT1I2dNWxoauO9uo4zRVRGv1BHfQwy1oqmVefpq7SEqTMnDNheWX/ALv7R3lhd5+j1c4LAjrYCbM1vyZGBiRpq4x6zJr8PkwZQs6k6FGzi4HRLA+/JIM6dA+M0DeSlSvz1oRmT98wptBYYeeKc2mAUpYfUXNNZYFQM2oFVdFckxlNCAZGMmRyOYTLj8MUUrOWMRL0J9oh1UIqEYfF50DLMV5h1DZ1s02att9XoKIxhZ7flKz3+NWFfZYoenZIahlG+guEXr/VilvrV5vZ5YxLzsxPNHvSXZ4LpEDY82sqLBPH6u7HZGj7+WdAZZurvnYAiTEm782rECve7tmUGZWCp+2b41bkyaiGXjpqElQnwgghWN/VCMCBoSx3Pru7MFBCRgrFWm6UHd1eIWxBfZKahHfCP65nxu6n88JRskH+X79YCl6Qj+I4OiA1jtZkXNO0AoEgbGjs7RvmuV29VpnCH9dyttCps7f2OKkR8pfx+aCB1EabmIxYAlq2aQ3dVSsZNtJ0V61wytnXiUY+mbdeJPgM5M0Qg9M15HGiI2oa/igezqLM4H7PdkPT6KlYyEC8lT1168m0Hp+/zs1R9NwIOd0nZCGjeb7RJLV+WT1MKqKxqk0+A9zm3oXXt7xnnAi57uvGilQZqmjK56ezU25gBIezd7dFt5Naa/lj7eJCyPtcD+WvIYCOE6GyXUaptbA1lvk2QP9IhnjYYElTGnI5KzKs1PT35GLQaI2vHvH6wmaHCyJUFkMJaLOV1/8IV59PX6iSM/b9PR1z5inhTFEaVXOkkJYdhqvOmbCQBrBlRRO3fOwUTltYx5dufoYLv/sAL+8tDBuuUBwp2IsQlfGQa0VU+nQsbEyxpK3e8iMTZHOmJ1iGx8LONB2NQ20yQlUiTFUiTNjQmFsTKzRxMyLEwjrLWipcIcezHlOkF/YN8+fXD/Kn1w6CgJGw9CPqS83NJxsuMJWzfsfnk2L/hH3+u4bn8/TOHgZGMgzGmuhPtHOwYlE+wqPrPNmcSc/gqBONEMgLhpZp4HFza1nWXueZJ+9sepNjQuQOINE3lP9bD7uTDWuMZk2yeowDVcvR7MlsgHkSIKPChfLaD1ujmbN89Eyh5fvWJ8gaul9AG0Yb6YGhgx7NJUAqaniCWdjakK76JFWJkCOgSaHTij4pfBo6v4BQIFjrrGyw6h3QVnvyPqaAFsBITkiNiyUIu33Z2qoKNXHu/ZFIlJAm6Bt2+cK5qzU6iCYEWT3CgUGpQRiIt9BRE8fQNUeodJpomrDvRXjiF55cVQ67/pz/u6IFFmyGdHN+mx3FcWQAQnE03TYfBr90HjY0T/CraKhwWizMDBl0NA2Go3UcqFpWUEYb2AtDPbDzcVc1LD+y0SwCgZ6sJWvE2du4gdFw3hfLHX1wOCK1v27T2ljIf01awkqRfIJvWtzAosa045sm/NpzXJN/X+RSQxOMhivorl7JSKSKTLLZuc6Fo0ErFNDc2DkgtYD7scDkNkBgdsgMwpJzoePEvDbMNMnlTJ7bM+CYBRdDt+6dfBJyGAlbPrOVHWDE5GfuyfmDYlXQdpynXifO8/pg7rPSBoR1Td5vuQy6lk+dEAvrULcAll9QaLmQGWZ9R2GE3yCUgDYbefVh+Mm59GlJNu3/BG1z5ivhTDEx6hfDe2+QE5crt8iX4QSpS0X4r/ccw9cvWsXzu3o56xt38+P7XparugrFEYYtcLmTLAtMK3+WAD1sCWgGIxmZFSgUTTGvNonpWnF2C2gOmRHpczGwr/CH9YBIermsZwLeY7mR7DgwKIVDPcqO1i0MR+vyAppf8LMFlpB38m0LovUpK1GuHmVwNMtTO3tAaByoWk5OjzBk+dWYTj3ykdw8ofZ9AoSuCRlZzm1mqBskErIe7qiHA67w4m4tnykMRrM5RxNZOOkLmNpk8xEjDU1jbUc1Wavto4ZrwuQTZEOWVsyeROvZQYQdWKFaRm90gls4/9j1lDmVKmJhDE0j1fuCs+/x16ykt74JaoGxd7k9AAAgAElEQVQ5YYAGTc/J38+ZhfvtvFWlymf2JL472k5/3UqwAtK4g4vUpwsn5J4+FxqapjlJwgv2h2IIAd1Vy9lfvZodrVvoS84lZUV79FtgCGHC/pfkl+d+V1hp/z0RSXq/21EcM4NgRF1JnCnoGJnKIf/7jhbUtaqi5TJOjkKHGm9uWX3PU/D8LZ5EynYXjGYyaBqOX96oz3/b0aBFUgzFpDbHXctkxGBP3Xr6knMh2eASVggkGTE8ZstBXn3OQsCo1y/NPY/sqk/RkIrmTXlzGemDFqBBc5PVo2jkBTQZqbHIwoH/+nVfN8N9cqw1zRm2TM7kFSs5dq8VjKcYZseJ9CfayYp8GpGsEWdH6xYZvCjdLK2KANpPgEi6UKDCKzR6zdXzWtOwlbqjMhHjeDsJNxTkHyQ7SpUREL02ACWgzTZeugt+8lYOijSbuj/JggWLueovjlPCmWLiNK+C998oH9BXboG9L4x/jA8hBG9d3cKtf3cK6zpr+OKNT/HO7z/Iq/sHxj9YoZhF2JNI3cy/XAVmPhKhHiaka2QtgWZoNEu8ZQlVizcwmsyv7udMk4ztg2ZPAp/+lTSbDMIIiDaWy3gmwHaEQ9M0XQmmrZ+wBQ7/pKqiRZrz1C7wntpqZ8QYe3owPGpPMm2tXl5Ao7bLdUJL4Iqk5GSsXvqJuCdqx3fWsa6rmbaqOPPrkhxnTXCWt+a1DCFdY0/denbVn8xgrJEXtQ56U3KSXJCg2a91AumzNG+jZwJWEZHH9SfyJoIFGjTLBy1n+d2kDz4HL98td1p9GrG0LulYyDMVdmvHhBDEB14rrJfF6rYqTuisQUR9q+t+AU3ozsQtZ5pFNWilyGdnLWti/bwau4IMpzscASYazp83YujUJb3Xj1coNgu0U86YVM2BxpUyUp8WdnztFrdUO8Jhhc+ETJimd6I81ANZ13UdMIn2nkCX91ZGpqAwMTERgSaBuhCB+agyLpNaLTeKqRneNjcsly13Ik8WntsRKrJW+BRLSPFHZnTu0SLtSlXWkItWY7SshrknO4KfGVDvIMyABYtiGtZoSCca0lnZWsmSZstPz6qXyAzJ6KHWvbB+XnB0z6we9V2AwtFMFvpYyu35yLiuAu0nOH9GDI1oSGdf34iT7mQ8jFSNXFCyvsfCOmcsacjXu2EJNFmmiBWtsODMcVc23NdAPopqBs0R0CLeqJph3+JBZght210l1V8JaLOJZ2/G/Onb2WM0cnr3p1m9YgXffe9aqU5VKCZD43K4+Ca5unzlmyccOMSmIR3lRxcfy79fsIKnXu9h09fv5poHt5f8AlEoyh37ZezOtWPn8TKRE8qQLhiK1PKH7dJsaDiTg4pWjyFSNgdDcRl6PhsqEsTBjcu/yMHKKQbSZMd0aeSGfcmlnflEkNCSbiqYkNjJoINCd7uxJ5lZ3fKZ00OOEDu3azksPQ8alkp/Dms/y853wu5r1Z2saK3kmPYqalIxwuEwDekolbWNNFVIobSpIuaEFDd0IU2uwmnS8Qi96S4nGXiBWWAxDVq82uNb0rlsHQPxFkdrIY/1tjusuzVoAs0loNsT7nQ0xIauOpY0pTk46ErCLDRnkhXSxn4WttfEZaCSjhNlPYu1RdNdASkKBTR73PzXQRBhQ/NoJt2Tdr8/XUdtwvPdI6zoEUKaIKPnFxOcMYnXFGhMAWoq8+NQEQs7wUIEQgrTbi3Z87fAy3dC/145jrbWpzo4/5wj6Az3gB72he8v1KDZLGlKO38ft2oVByukCZ3IWWa77jb77qegab0TNCibtXLUFdOg2QKahq4J5tcnnUik/Yl2Qos2sWVlMyd11XrOW7KxilXXlsoYSxrTrG6rHLP4pqWNzHGPt52mYPeTQH5BqFiEUJnr0Ns+vxmrgzWWS5rTLGvOL8jQtCqv3UIKbrXJMAcHRxkaLVFAs+rn1tDGwwZ1qbE1gEGkoyFCuuYV0Ozz5jKOD6LmD+4S9j3jD75asvmxEtBmC49fh/nz97Aj0skZ3Z/k9ONW8o13rC5+0SsUpdKwFC7+NSDgR5vgxTsmdRohBBeubeO3H9vAMR1VfO6GJ3nfjx6WEd4UillOeOdWQtkBhMtUrjktJ1EhXYNEDUbnBnoqFjkTgnn1cvXUvU5hr/7KhK7CCctdlEiAv0J93vfCNjGMGFrBhGlNexXVMTuKYmkLedkADdru+hPZX7XKU86enPSm5rO/ahVmutUxPQrpQk4K6xeDUSSsdNMKwrrl+2VPWLrOhDkne4rplsmeW5BorvBqFQsmPEFmoXbUP9e+aEUD3dWrvEKOzxRUc5s0+fvQpfGoSoQLBMWcMJx3dFIbZWlzetzADoRijpYRKBw3zSBkhXDvqk8V1NdesB0YKS0QgcfU1NePdckIHTVyou6fi3uutYoWhhrX0pue72zyRxjMZL3SRNgXzc+TT/vAK9D7hvcHBw/AS3d6fLxoXl3YIPCE4EczCoNDFGmHO19hZTzMcET6KwlyaNlhb//4ow/66xCKOddlJiNNkjVLUPKbdDoaLqFx9opmljZXQCjGroYNHKhcWuijN0ENmn0NxcM68Ygh/VybVkHLmnyZjvVWw9sLj7dNHEfl/T0akoJssWAuWT3m0YSZQuT7zn9MRSukmwk1LvYuCvkTSSMXDQZGMry4pzBxeBBRQydiaCxtqRi/8DictqjeE00SyPsi7nuJ2h23Azhj7BCKQdQlEGeGQdPzVhRjoGb35Y5pwt3/CddfygvRZWze/3Eu2rCSK85bppJQK6aO+kVw6e3y4XzNBfCHKyd9qpbKGD+55Dj+5a3L+MP2bs782t3c8McdSpummL1kRzB6XqWy9zlP3qC6VJS1HdXOREyv8GqkkhF3/h3rVDnT8scQ1CXDkBvHH0EPEHB0A1JN1rmtQABCMKfGq+Voq47nI4YFadACsM293JPC0XAlg4kWT7ACZ5IpBIOJFkZyphMoqCDhcxBB5lzRdIFwZZs+xV2WIsmo99gCE8ciZqGysKs/hcampY1sXFQ/fn3xJtGWxxfv046aBKbQqLN8+TBzxEJG0cAO3vMW9+Gz+62tOk48wLWhOi7bV6qWwC3c+qcU6+fXOpH4/IKlXwAbTrU6dR2INefPa0VQ9GuNIj4NrROjZbxpzdCBwoP8uHN76WHnmh5Pg1ao2XCZbOaGx5xzeQ6tnmsJaPJrbN8TVgLsYrnDLBNS3zsyE5JmwWHf/VSfipKOhVjQWFqwiXAo/4yQFaqE2vlSA1nRCnULpT/W8gtkgAw/9j0ppMY+p0doTEepiIVYO6eaqrjvGSW80RVNK2IsBN2rESkc1nR5E8r7/XQpotVfdDYsPodTF9Rz6gLvfaxpgs3Lmrw5Dg8R97NtWZuVjqD75bz8H/ScTfqeL6mmcf34QAlo5U02Azd9FO74Zx5MbGTL/o/xoU2r+fRZiwLtnRWKQ6KiFS75rfTTuPEjcMvnvHb/E0AIwXtO6ODmj5zMgoYUH/35Y1z+34+yv7/0HCAKRdlg5uQkT486uZWK4Y7iZ09iq10TmGzOJBaWE/d0RPcErwik2IvctepuU5MMEOasBLwFvhBFsDUdibB38m9o3pVov9/Ont58v4T0Et5PJb7D7Ih16ViIhY0p1rRXkfIJaP74J4FaO1szUG0Fd0jUgSb9WlLRAI2bDxPyAVecHy7uC7WqrZJNy1vz/uGmrZksQVB2R9Ys8K8b2/8qFtbZsryJuT6TxGK4J8xjzSsKgkv60yS4BLDu6lXkFr9FTrwtk1a/wOj/rZx9nZY6tYlVFd/n1kA6AnmwD5pHQPVdR+5gN73JzkDTNMMJyuITqtuOd7ZlcybCpUFz01mbZOPiJkdTGYS/r8KGxmkL6z3PmrGYV5+mozqeL5/I5zmj/QTp6jAWRsQRnkbC8v/jO2sQQiZ4P76zeqyjAY2VbRUsaUpTkyzyPAtFYb4rEXjAdR4UYZNQFIwIFfEQFQHpIKYa27IgYuievHiOgBYN0NbZzx47fUe62RVcqThKQCtXhvvg2nfCH67khuRFvHPfJfzjuau4/LT5SjhTHD4iKXjntXDsB+H+b8rgId3bJ326jpoE1/3VOj65eSG3PrWLTV+/mzue2TWFFVYoJocQ4rNCCFMI8a1xC1sCmhmKec2nbFxmebafCORX5+vTUTYvayQVNWSUxRz53EQ53yJIzTxHOwYESB92A6zcPi4BLRkULKpmHnSeBqmGwn0B2IoxQxec0FlTtNwTOw56vvcMBkdcHBP3RLEItrBnmiaLGtO0VcdJBgiPYzLnZGg9Vv6drJOags5TvGXmvylYe+CiQIMWMOE+pqNKmqiBb5IpO3ZhY4oFDeNoPtzXRM5nqliCJrQkDaZ9Oo8GqXg5/7zDr02qdwtgQmCEIp7Q96loyNHAzqsrXCwQdij38SqcHZE+enM3FC8TikLYEnh0mUtrX80x5JKNAWkUxtCgue6tnsrFgffXwoYU7dVx77FCg3ACs1GmjsiZVhAV13XaWhWjoybBwsaU0/ZS82NNlJAw85pcCNROjUu0CkFwwBH3/ddVLxdRqMgndjaFIGLodI133bsJEND8i0Yzgb3AksnlPNr4mFW3QRGgrYtVymdO+3qZH62ijYF4gG+xDyWglSO9b8CVb8Z84Ta+Gb+cj+9/K1+9aDXvXTdnpmumOBrQDdjyFTj/B7D7Kfivk2Q+msmeThNcdup8brj8RKrjYS65ciufuf5xeodKCzWrUEw1QogTgEuBx8crC1gRCk05KQ3SoLkmE8WSPkcMmRsskzPJmqY0dzLNwgl40yqYc2IJjSgU0Nxmie3VcbsSkCguaPmxNWOGLmhIR52IiuOZ5vUMjaIJwZlLGksX0OZugGVvG7OIPflz+zD5fb3GlM/mnCSF0/EWNmNVwf43WHmohChJQGutissgD/79tuloopr0ss0ADEeKaB7c2iH/bxwGU3FbuAhK0pwvk/97dVtVgYA2ry7JyV15gTvIHNC+H5yw/YvPcSL12fs8RwVpC0f6pU9PkJ+hvxxArJIFDUk6OjppWP6mQh801/fCOnu/t9cU+g1FQ7qTkiJ/mBWx0HVhanj7NxYyWNVWKe9ZO1hQxrv4c9qiek+fThq/4BdkNj0eVqj7gnsAX6CVZrmIQutxJWvti/2enypfQvOmiqkzXSyVsK6R0yLSxNt1DTZVRKmMhWmpHsO3LFEj86MJwUCiVYb7HwMloJUbbzwJPzid3J7n+UTos3yndwM/eP9azlvdOv6xCsVUsuLt8Nf3SEf/X34AfvGBQsftCbC0uYJf/e2J/NWGTq595FXO/JrSpimmHyFEBfBT4ANA9zjFARkmu2dolP7h0WANWon+XYamkc2ZmKZpTdZMeMMlI7oDZizYLDVfxdA0oobuFdBcgtGqcSK1FcMW0Gzfr8pYaZO5vqEMiYg+sajC7vYWwdZwZMYIWTe2Bu0QLU4iabqWrmHJCZvJCWNivt/uyaxttti0AmJVvN58Jntrjw8+LpyQUTDbjvNGdATvZNvwCQaTxG7TWE1za4iCBBUhBNWuCXSQOaDw7zMijr+gGSQ0+HNIORUu4ZqssDQUkTSGrrG4KR0ogLq3+atcihlaINY9abjMWYUmPAs2nuvIjvSXGfKcJh0Nefp00vgDEQWZ4Y2LsDRoeuDYFqBpkJI52ypDk9AMjpdKAZzFo+nE0AR76tbJoEmu+8/QNebXJ51opFOBEtDKiad+BT88k5HREd6T/Sd+n13FtX95AqctLM2BWaGYcqrmwMW/gVM/K3M1fXMt3PNVaYI7CSKGzmfevJhffmg9yYjBJVdu5aPX/lH5pimmk+8BvzBNs+RwpSMjUmsmzFx+Zd5NCZMJyCfFzZlWwAAz501O7Z7QR5Jja76ExpLmNGvaZRnT9GrsJmsKf9L8WubXJ52Jq62Vq09FiY4RNbh/JEO0FP+qCdJVn6KjJkGHTyhw+6EFCk1O1LhDnOYIQbh1Nel0JTXpGPMDzPOK4r4ubE2prfnUQiA0lrVUBGtJNL2oRg+QkTwXn116XcaqphOgYywN2sSuJ38QkaLnStRA43K6K6UflDvZdan3VSCtx0oN3Tj1HnGlIyhoo6aTihisaK1k87J8yHcnb1tRQUeeJxoOOec8OGezRwD2/JQd1KaUtBuTwS3Uz9sozXwnirAFTDPQvzFi6DKqqJuKNla3VXLsoiLpEMb8vbGfJaeVGNhnqjF0jYwRYzDRIq/Pw+hypAS0ciCXgzu/DNe9lwOp+Wzs+QKvxxdx/WXrWTnJVVCFYsrQDTj1U3DZg9Lp+/YvwteXwz1fkQlEJ8Ga9ipu+vBJfPhNXdz0+E5O/+pd/OpPr6tIj4rDihDiUmA+8I+TOX5BbRQG90PNfFh6fn5HiRNJQxccGBhhNJuTWXNM06sFmYgwYeVCsxMpTxWV8XDehwop/Jy2qJ5jOqqoSUYK/HDckdMiQU78h0jY0FjVVllgNjm+6dcUCWj22YRgcXtzSQFFHNzXxbD9rPRO6ObUJCamJameJ68/X4LxQ8ExcRxrsjnBeWip5pIA1C1kmBAHKpfSn+jIby9mxljKmGpaPijDGLj73h9h8OzlzSxsTBPWNcf3aOOiek6cb/mZdp0Bi99StH5C6E5QCWFEi6dY0DSZ/24sv7pDwR10xq+RLRkZmfGUrlqWNqcL9m5e1sgS//ZELWLZ+YgS/V89FHmm2lYC/siW04X8WfvZIvJBmA4DSkCbaUb64RcXw51X8FLz2Zyw82PUNLbziw+tHzOqj0Ix7dTMg3dfBx+8HVrXwu3/nxTUfvcPsPvpCZ8uYuj83RkLuOnDJ9FWFePDP/sjl/5kq8qbpjgsCCEWAlcA7zZNc1yVrRDiL4UQW4UQWwcGpQatWuuTQpUVAdCZAJbo0+HWKhwczMqJk5mV95YQJZtKygpafi6+fEErWyulk/4Uko6GHC1VpS9SWsjw+tlNF+P6uTlRLqdwhbu2K//3vI3jlw+aZFr1sUOTTzhbjm5A86rAPFGTxe7LsTVo3rKHQpB2dzSToz85J59bCgr7r2GZ9CmsW3jIdXBOmc4vkJSS8DxlJSx2CIoYah+n5QU0/6kLQtOnmwqTGk8VTaukNczicyZ/DmvMDJ+p5rhM5JlWwnHr59UytzbhydM4bTQslTnk3O2PT+2z1o0S0GaSA6/ADzdhPn0jv2//MBtfeifrFrTws788wZMwUaEoK1rXwrv/By79vYyE9tB34TsnwPc3wiM/gJ6dEzrdosY01192Ip/bsph7X9jLmV+9m2se3O4kwlUopoh1QC3wpBAiI4TIAKcAl1nfPQ9d0zS/Z5rmWtM01xrpBkyhExraL3cmrBX0zlNh0ZYCh/ZiE5h0LC/cLGi0VpuzozKUvhEd16zHixVUwTFFlP/PqU1IJ/3DhC+dlWeyGpinaKYISENwyLi1naVoIjQN5vqiRVrXxgmdNZyyoK4sojLnQ8UXL2M/jgNDnU+QIKF0xH9hQaEGbHQAUoWRGA8bTnK2Cf6e69pz8n9pglhYZ+2cajYuqi85T92UEI7L93YJGsXiWH0xXVYuRfq8Ih5iRWvlzNw39YtlNEZcixlNq+SiAUwo+EpNIlIopPuY+ZiVRyvP3wbXfxAzl+FbDf/KV57r4J3HtfHP5y6bUIhchWLGaFkDF/4E+vfC49fBH6+GX39cfppWyiAH80+XD7CgVUYXuib44MmdnLGkgc9c/wSfu+FJ/mfrq/zzW5exolWZ+SqmhBuArb5tPwaeR2rWxtSqmWhSA1bRmp/oRILDRp/l8ldxU+ES0CpiYei1vuhhOfk3AyapxbDqEDZHWN5SQWPF1ASMGI+miqhHy+0W0Px+YjPK4dCgTeZcBf4+ed++8DjPxenCDrIyVgCUZMSgKh5mWcvkTboct8CAfpxTk+D53b3ejfVLpAtI/SJ4/bEpNet0s6ylgte6Ayw37AWT+iWlnUgPyyBCTkN1x2zSvk+mMmnytFIzTwYJq5470zWZcU5dUJ+PmKvpctGg9ViIlx4t152OpRhKQJtuclm468tw178zWruIyzMf5dbtKT63ZTEfOGluWaymKRQTIlEL6y6DEz4kTR2f/x08+1u4+z/ktW7E5Opdx4nQsU4+yMLB5rsdNQl++sHj+d/HXudffv005377Pt51XDuf2LSQynFWmxSKsTBN8wBwwL1NCNEP7DdN88mxjo2GNHQy8vlcwku4mBmYOwy+R7Ojh+R9lJlAsBxbOBzupbPlEMJZT5C26jgtlTFufPz1gn1TYf42ZUxVkJAgxkqSPB6Hoz6HiK1BG8vEUdcEGxZMQch3CPQ1XtKcLhTQQjFos/LXzT254JipYl5dMjA3G5rmaEzGpGWNDJzVt8snoGlkrQiKZaVdngyhGHSdPn65o4DAhNhVHYXbDhEloE0n/XtluPKX7mR/1wWc+/L57B/R+eH7V7Nx0SScKBWKckIIaFgiPyd9DPr3wfZ7YfsDsP0+uPvfpYZAM6BxucyT0mZ9Ktqcl5oQgreubmHj4nq+dutzXHX/Nn7zxE4+fdYi3n5M25jO5wrF4SAa0jnG9usqFvq7BIpGttMMqXWeCIk6aWbXsHzS9ZksmibYtLTRsXbqqk95oipOF0KI4oGFDpcgtPT8Q9PKleEibD5IyOH9naaKGM/v7vUuVLg4c0mjTADcuxwG9h/eykwl1VaUwudvk/87PmgaqViIvf0j1CTVAmNJLNgMo8oPHZSANn288hD8z8UwsI/H1/wLFz4yj5pElF9+cC2LGgsj4igUs55EDSw5V34Ahg7Cq49IYe21R6RJ5MPflfuSjXlhrfU4aFpJOhrl8+cs5cK1bfzT/z7Jp375BD97+FU+f84SVk9xEATF0YlpmqdO+KBDCIXtMSHLZfJ/TyacuKaXFqjiMOHWCBREb5smipmSAvlJ8kTMRkthzJxrRajuhP0vyb/LUINmX5eH24JncVOKefWJosFkZA49HaJTFwRkWvFrbYVOTSJCRTSEMZHon0czkaT8KJSAdtjJZeHer8GdX8KsaOVny3/IZ+/XWNNewffet1YFA1EcPUQrpImEbSaRzcCuJ6Ww9urD8OpDMtcaSFv+lmNg7iks7jyV6z64luv/tJt/++0znPed+3nrqmY+uXkRzbPVnl8xezkER/uQe3LvEdBmufnTDDGmSWXjcvlsCZfBZK+iLS+gHWri7MNAwkouPhWR8QJNBS2EENMa6XPacYLI2CEv5WKOiiugmAxKQDucdG+D6/8KXn2Q4YXn8uG+9/O7B4c4f3ULV5y/fPbbJCsUh4IdLrp5FRx3qdzWu8sS2B6Cbfda/pr/hggneVvHes7esIGf7FvIf2x9g9/++Q3+csM8/vqUTuJh9ShTTBMTiNTlx2Oem3Mlj1UC2tSTbs5r72caUcT3sEyYW5ugJhmZWD62IhxKEJFZjy2gZS1f0ilMhaA4+lBXz+HANOGx/4abPwVC8OJJX+VdD7XTPTjCFect553HtalgIApFEKkGWHy2/ID0Q9h2L7x8F7x0J5Hnb+FS4OL6BdzBWr5zxyKue3gxf795CeevblH+aYrDT7HkuRPFrUGbUHh9xazD/b4vw3e/oWtTIpwd9dja9cxwftucwxfcRHFkowS0qebAK3DjR+HF2zHb13N182f54h29tFeH+PFfHD9jtvoKxawkXg1L3iI/AN3b4dmbCT1zE5u2/5xNkSz7s9X85v+u5tN3nMrGM8/jzOWtSlBTlD+H6oOmmD14tGbq2XTEUjMP+nd7Q9GnVAA4xeRQb4WpIpeDR74Pt30RgO5T/5XLn13D/Xd2c/aKJr50/nJSyklUoTg0qjrghL+Wn4H98PytVD1zE+947lbe0387+6//Mrf8ej31x1/I6lPegjikxJwKxWEk2QC9VlJ3ZeJ4ZFPmJo6KKSIUm9HAPYojCyWgTQXbH4DffRZefxRz/un8pv2TfOr2A5hmD/92/nIuOlaZNCoUU068GlZehFh5EcbIANkXbqPvvms5eccdJO65hb57E+xvPZ2mdRcR6noThKYnka9CMRanLqyX4fYjBux8TG5UGrQjG4+ApuYCCoVifNRb4VDY/zLc9gV46gZINbP3zG/x8acXcNdv9nLc3Gq+8vaVtFVPPiSzQqEokXAcfclbaF/yFjLDA9xz6y/pf+x61r1yK6FX/5cRPUFm/pnEV54P80+HsLovFROk6wwolnNrAlTEAiwplAbtyMY9vkqDplAoSkAJaJPhwKtw3zfg0atAM8hs+DTfy2zhG7/ZQUg/wOfPWcL7181RfjAKxQxgROKcfPZ7yb35Pdz9zA4evuMG2t64jU3P3Er82f9LRo+Rm3c64RXnQdeZEEnNdJUVs4HoURydTnGIlHeQEIVCUX4oAW0i7HsR7v0q/OlaQJBb9S5urrmYK+45wI4Dr7JlRRP/dPYSGtLKlEqhmGk0TXDqklZOXfI3vLLvEq5+dBsvb72FY/rvYfOzd1P33I1kRJi+tlNIrn4bxsJN0mxSoZguFmyCwQMzXQvF4UZpzRQKxQRRAtp4ZEbg2d9IbdmLvwcjQm7tJdxV806+dH8vz92/k2Utab78thWc1FU707VVKBQBtNfE+cgZSzBPX8zjr13MNU/vZOcTd7K4+/ds3v4Ila/cSg7BG/GFHGxaT7jrTTQvP5VYogyS3CqOXCIppcE9GlAmrAqFYoIoAS2IXA5eexie+hU8cR3074F0K6MnfYIbQpv59sM9bNu3k7m1Cb79rjWctaxRmTMqFLMAIQQr2ypZ2VYJZy5mT+8l/GHbXm7+831EX7mLBf1bWfnCVYRe/BHDN4d4VF/Ea6kV9Netweg4ltbmVubVJahLRVTgH4VCURpKg6ZQKCaIEtBsBvbDKw/AC7fBM7+Gvl2ghzG7zuCF1rfx412d3HjPbnqH32BlWyXf2byITUsb0ZVgplDMWupSETYvb4HlFytfUhIAAA29SURBVAIXksuZvPrGHvb8+fdoL99J3b6HWXnwZ+gHfwovwIu5Ju42u3hWm09P5SJE/RKaGxvprEvQWZtkbm2CWFitlisUChdqMUehUEyQo1NAy4zAnmfgjSdkmOPt98OuPwMmhBJk5p3OM1WncEPfMn77Qj+vPTZILLSLs5Y38o5j2zl2TpVaPVcojkA0TdDRXE9H80XARXLjcB+5HY/S+8IDVGx/iC27H+WC0bvhIHAQXn22jqfNdm4123k610F3sotofSdz6tJ01iWl8FaXpCkdVZp2heJoJd0CI30zXQuFQjFLOPIEtFwOhg7AwD7o3wsDe6FnJ3Rvg+6XYf/LmPteQORGZXEjRnfNarZ1XsaDuUX8truFpx4fJpsziYe7WT+vlo+8qYuzljeRjBx53aVQKMYhkkTr3EBF5wb53TSh53W5qLPrCZpef4L6nU9wxoH/RZCDEci8prPj1TpeyjXwgtnIbWYjr2uNZCs70ataqauqoCkdpbEiSnNljIZ0hIpYmFTUIBpSGjiF4oijY92UpGlQKBRHB9MvcTx9E2y/z3pQma7/cwHb/P8DuVEYHYDRQRgdyv+dGYSRfmmqaGYLfzcUZzDRxgPdKZ7OnsXTuQ6eMjvYNtRIrk/ah9cmwyxrSXHZ0laOn1vDsXOriBhqsqRQKFwIARUt8rPgTAysB+noIOx+GnY/hbH/Jdr3vUjz3hfZsP9e9Ey/PLZXfva/kmZnroo3zGq2m1U8ZFaznxTdZoo+LcVwuIpstIpcrBphRAnpGoYuCOkaIV1g6Bq6pcX/zJsX0VQRm6HOUCgUJaMsbxQKRYlMv4D22iPw6NXWg0rI9CBCs/4W4/+vGRCKQygm/4/XQCia3xavgXgtJGrl34laSDZCsp693YPcf/82YmGdxSGdNSGdxoooLZUxWqpi1CTCynRRoVBMjlAMWtbID/LRFgK5wNS/B/a/JFN19OygumcHFQdfp+vADkTfnwgN7fOeKwcMyM+QiDIgEgyKKAPEGCAqPyLGIFESd3RAIg56CPSwfEbqYetjPeJNa6HLzOU/+L7nspDL5P83s9bf9nbftupOOO0z09a9CoVCoVAcLUy/gHbGF+VnBmirjvO5s5fMyG8rFIqjFCEgWS8/7Sc4m3XrA0B2FAa7pWn2wH7r/30wuJ/owH6iwz0w3CetBEb65Ge4W35/5iHIjshzBFkPTLi+ugwLrhnyIzTrb2ubvT+XOfTfUigUCoVCUYAwJ2ATLYTYA2w/fNWZUmqBvTNdiSlEtae8OZLacyS1BVR7yp1aIGGaZt1MV6QYQohe4NmZrscRwJF27c4Eqg+nBtWPh47qw0Ono9i7b0IatHJ+gfoRQmw1TXPtTNdjqlDtKW+OpPYcSW0B1Z5yx2rPnJmuxzg8eyT1+UxxpF27M4Hqw6lB9eOho/rw8KKyJyoUCoVCoVAoFApFmaAENIVCoVAoFAqFQqEoE45kAe17M12BKUa1p7w5ktpzJLUFVHvKndnQntlQx9mA6sdDR/Xh1KD68dBRfXgYmVCQEIVCoVAoFAqFQqFQHD6OZA2aQqFQKBQKhUKhUMwqlICmUCgUCoVCoVAoFGXCrBPQhBCXCyEeF0L0WJ8HhBBbxig/RwhhBnw2T2e9S0UI8Vmrft8ap9xyIcRdQohBIcQOIcQ/CSHEdNWzFEppS7mPjxDiCwF1e2OcY8pybCbalnIfGwAhRJMQ4iohxB4hxJAQ4ikhxCnjHFOW4wMTb085j5EQYluRuv16jGPahRA3CiH6hRB7hRD/RwgRns56B9TpMiHEy9Z4/EEIcfJM1qecEEJ8RgjxiPUu3mON3TJfGWE9e1637rk7hRBLfWWqhBBXCyEOWp+rhRCV09ua8iDovan6sDTGe36qfhwbIYQuhPhn1/PuZSHEvwghDFcZ1YfTxKwT0IDXgE8Ba4C1wB3ADUKIFeMctxlocn3uOJyVnAxCiBOAS4HHxymXBm4FdgHHAh8GPgH83eGuY6mU2hYX5Tw+z+Kt2/JiBWfB2JTcFhdlOTbWA/8+QABbgMXA3wK7xzimbMdnMu1xUY5jdCzeOq0BTOC6oMJCCB34NZACTgbeCVwAfGU6KlukThcB3wCuAFYD9wM3CyHaZ6pOZcapwHeA9cBGIAPcJoSodpX5JPBx5LV8LPJ6vlUIkXKV+W/k9XEW8lpeA1x9uCtfbozx3lR9OA4lPj9VP47Np4DLke/FRcBHrO+fcZVRfThdmKY56z/AfuCviuybg5wUrJ3peo7ThgrgReRL7k7gW2OU/RDQA8Rc2z4H7MAK/DKL2lLW4wN8AXhyAuXLdmwm0ZZyH5srgPsmeEw5j89k2lPWY+Sr6z8AB4B4kf1nATmgzbXtPcAQkJ6hOj8EfN+37XngSzPdn+X4AZJAFjjH+i6AncA/uMrEgF77nY2cSJvAia4yJ1nbFs50m6ax7wLfm6oPS+6/MZ+fqh9L6sObgKt8264CblJ9OP2f2ahBc7DUse9AvhTuH6f49UKI3UKI+4QQF0xD9SbK94BfmKZZysr3OuAe0zQHXdt+BzQjJ2wzzUTaYlPO49MppCncy0KIa4UQnWOULfexmUhbbMp1bN4KPCSE+LlVv8eEEH8jxJjmiuU8PpNpj025jhEgzWKADwDXmKY5UKTYOuBp0zRfdW37HRABjjnMVSxASNPKY4BbfLtuQWqMFIWkkJY53db3uUAjrj607r27yffhOqAP7zv8PqCfo6ufi703VR+WxnjPT9WP43MvcJoQYhGAEGIJcsHgN9Z+1YfTyKwU0IT0IekDhoH/As4zTfOJIsX7gL8HLgTeDNwO/FwI8Z5pqWwJCCEuBeYD/1jiIY1IEy03u1z7ZoxJtKXcx+ch4GLk6v6lyP69XwhRU6R82Y4NE29LuY9NJ3AZ8BKwCWmK9m9Ik4xilPP4TKY95T5GNmcgX+4/GKNM0NjsRWpkZmJsagGd4Otlpq+VcuUbwGPAA9Z3u5/G6sNGYI9pLbUDWH/v5ijp53Hem6oPS2O856fqx/H5MtIU8SkhxCjwZ6RG7TvWftWH04gxfpGy5FlgFVAJvA24SghxqmmaT/oLmqa5F68Pw1YhRC3Sjvaa6ajsWAghFiJV8yebpjkygUP9CexEke3TxmTaUu7jY5rmze7vQogHkS+A9wNfLXaY7/uMjw1MvC3lPjbIBaatpmna9vF/FEJ0IV/IYwXZKcvxYRLtmQVjZHMp8Ihpmo+NU67YGMzk2ARdLzN9rZQdQoivIk2ZTjJNM+vbPV4fBvXnUdHPE3hvqj4cm1Kfn6ofi3MR8D7gXUjhbBXwDSHEy6Zp/tBVTvXhNDArNWimaY6YpvmCaZr2zfgY8LEJnOIhoOvw1G7CrEOu1D4phMgIITLAKcBl1vdIwDFvULgSUW/971/ZmE4m05Ygyml8PJim2Yd8cBWrX7mOTQEltCWIchqbncBTvm1PA2MFcCjn8ZlMe4IopzFCCFEPnAt8f5yiQWNTTIs1HRTT3tUz89dKWSGE+BoyqMtG0zRfcu2yo8SO1YdvAPVuU17r7zqOjn4e870J7LPKqT4cm/Gen+paHJ//AP7TNM1rTdN8wjTNq5GLt7bQq/pwGpmVAloAGtJPoVRWIW/mcuAGZCS9Va7PVuBa6++gFbUHgJOFEFHXtjOA14Fth7Oy4zCZtgRRTuPjwerzRRSvX7mOTQEltCWIchqb+4CFvm0LgO1jHFPO4zOZ9gRRTmME0qx2GPkcGIsHgMVCiFbXtjOsY/9weKpWHEub8QerDm7OYHyf56MGIcQ3kCvuG03TfMa3+2XkhO0MV/koMkqn3YcPIP3I17mOWwckODr6ebz35nOoPiyF8Z6f6locnzhyUcpNlrysoPpwOpnpKCUT/SBtik9GOvQvB76EjPx1lrX/S8DtrvLvR748FiNv3r9HCgofm+m2jNHGO3FFPgxoUwXyJrkWWAacj4xM9/GZrvsk2lLW4wP8J3I1cy5wPDLKUQ/QMdvGZhJtKfexORYYRUYHnA+8HTgIXD7G9VbO4zOZ9pT7GAnkBPP7Afv+BnjG9V0HnkCmCFgNnI6MrvnNGaz/RVZ/ftDq428g/f46Zrpvy+EDfNu6fzYiV9XtT9JV5lNWmfOte+5a5IJIylXmZmvsT0BO5p4Abpzp9s1gv96J972p+nD8Pivl+an6cew+vBKZymoLco59HrAH+IrqwxkYj5muwCQvoO3IVdXdwG3AJt/+ba7v70eqvfuti2or8J6Zbsc4bfQ/nD1tsrYtR0bOGUKuln+eMgixP9G2lPv4uB4+I8jJ4i+BJbNxbCbalnIfG6uOW4A/WX39HDJ/i3DtnzXjM5n2lPsYAach/Q6OC9j3BSz/cde2duTCwQDStOubQGSG23AZUrtqa/I2zHS/lsvHGtugzxdcZYQ11jut6/ouYJnvPNVIn8ke63MNUDnT7ZvBfr0T73tT9WFp/Tbe81P149j9lwK+jpxjDyJ91K8AoqoPp/8jrM5UKBQKhUKhUCgUCsUMc6T4oCkUCoVCoVAoFArFrEcJaAqFQqFQKBQKhUJRJigBTaFQKBQKhUKhUCjKBCWgKRQKhUKhUCgUCkWZoAQ0hUKhUCgUCoVCoSgTlICmUCgUCoVCoVAoFGWCEtAUCoVCoVAoFAqFokxQAppCoVAoFAqFQqFQlAlKQFMoFAqFQqFQKBSKMuH/AUiIe/W5paJ/AAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 864x144 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"az.plot_trace(trace1, var_names=['ν']);"
]
}
],
"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.7.3"
}
},
"nbformat": 4,
"nbformat_minor": 4
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment