{ "cells": [ { "cell_type": "code", "execution_count": 1, "id": "4099c66d", "metadata": {}, "outputs": [], "source": [ "import GPy\n", "import matplotlib.pyplot as plt\n", "\n", "%matplotlib inline\n", "import seaborn as sns" ] }, { "cell_type": "code", "execution_count": 2, "id": "e355bc94", "metadata": {}, "outputs": [], "source": [ "import numpy as np" ] }, { "cell_type": "code", "execution_count": 3, "id": "c70ef438", "metadata": {}, "outputs": [], "source": [ "true_kernel_lengthscale = 2.0\n", "true_kernel_variance = 1.0\n", "true_noise = 0.5" ] }, { "cell_type": "code", "execution_count": 4, "id": "6efdbc47", "metadata": {}, "outputs": [], "source": [ "kernel = GPy.kern.RBF(\n", " input_dim=1, variance=true_kernel_variance, lengthscale=true_kernel_lengthscale\n", ")" ] }, { "cell_type": "code", "execution_count": 5, "id": "44f9be85", "metadata": {}, "outputs": [ { "data": { "text/html": [ "\n", "\n", "\n", "\n", "
rbf. valueconstraintspriors
variance 1.0 +ve
lengthscale 2.0 +ve
" ], "text/plain": [ "" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "kernel" ] }, { "cell_type": "code", "execution_count": 6, "id": "c5388fab", "metadata": {}, "outputs": [], "source": [ "X = np.linspace(-4.0, 4.0, 250)[:, None]" ] }, { "cell_type": "code", "execution_count": 7, "id": "7652eec0", "metadata": {}, "outputs": [], "source": [ "C = kernel.K(X, X)" ] }, { "cell_type": "code", "execution_count": 8, "id": "b059db35", "metadata": {}, "outputs": [], "source": [ "np.random.seed(0)\n", "f = np.random.multivariate_normal(np.zeros_like(X.flatten()), C)\n", "noise = np.random.normal(loc=0, scale=true_noise, size=f.flatten().size) " ] }, { "cell_type": "code", "execution_count": 9, "id": "27829513", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAD4CAYAAADxeG0DAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABBfklEQVR4nO2deXhc1Xm43zPSjGYkjSTL2mxLtow3jIwX8EJiQyBAgMQGGqWUmMQEWghJY6eUH1tp06YpSeM0tDEli0nSYMBOCEoKmMYYEjYbjBe8YHmV8SJ5k2xZ1kia0Yxm7u8Pe5TR6M6+Xvl7n4fnwZqZe8+955zvfOfbjtI0DUEQBMG4mDLdAEEQBCExRJALgiAYHBHkgiAIBkcEuSAIgsERQS4IgmBwcjNx07KyMq22tjYTtxYEQTAsW7ZsOaVpWnnw3zMiyGtra9m8eXMmbi0IgmBYlFKH9f4uphVBEASDI4JcEATB4IggFwRBMDgiyAVBEAyOCHJBEASDI4JcEAQhibhcLtra2nC5XGm7Z0bCDwVBEIYiTU1NrFy5Eo/Hg9lsZuHChYwfPz7l9xWNXBAEIQm4XC5WrlyJ3W6npqYGu93OypUr06KZiyAXBEFIAg6HA4/HQ0FBAQAFBQV4PB4cDkfK7y2CXBAEIQnY7XbMZjPd3d0AdHd3YzabsdvtKb+3CHJBEIQkYLVaWbhwIQ6Hg+bmZhwOBwsXLsRqtab83ioTR73NnDlTk1orgiAMRVwuFw6HA7vdnnQhrpTaomnazOC/S9SKIAhCErFarWnRwgMR04ogCILBEUEuCIJgcESQC4IgGBwR5IIgCAZHBLkgCILBEUEuCILhyUShqmxCwg8FQTA0mSpUlU2IRi4IgmHJZKGqbEIEuSAIhiWThaqyCRHkgiAYlkwWqsomRJALgmBYMlmoKpuQolmCIBieVBaqyiakaJYQNxfKJBGMSyYKVWUTIsiFsEholyBkP2IjF0IioV2CYAxEkAshkdAuQTAGIsiFkEholzAUGYrp/GIjF0LiD+1auXIl7e3t/TbyC9mpJBiboerzkfBDISIStSIMBVwuF0uXLsVut1NQUEB3dzcOh4OHHnrIMOM6VPihmFaEiFitVsrLyw0z2AVBj2T4fLLVLCOmFUEQLggCfT5+jTwWn082m2VEIxcEIetIheabSDp/tofiikYuCEJWkUrNd/z48Tz00EMx+3z0zDLt7e04HI6sMDkmRSNXSt2olNqrlGpSSj2SjGsKgnDhkQ7NNx6fT7aH4iYsyJVSOcBTwE3AJcAXlVKXJHpdQRAuPLI1CS1as0ymnKHJMK3MBpo0TfsYQCn1a+AWYFcSri0YHAldFGIhUYdkLMQ6NiOZZTLpDE2GIB8FNAf8uwWYk4TrCgYnm738QnaSriS0eMdmqCqLgSYh/wK0cuXKtMWoJ0OQK52/DcoyUkrdC9wLMHr06CTcVshmMj2wBeMSr0MyWlIxNjPtDE2Gs7MFqAn4dzVwLPhLmqYt1zRtpqZpM8vLy5NwWyGbyVZbp2AMUpmEloqxmWlnaDIE+SZgglJqrFLKAtwOvJyE6woGJtMDWxBCoTc24ZymHq+T0mq1Ul9fz4kTJzh48GDaj5xLWJBrmtYHfAN4DdgNvKBpWmOi1xWMjZylKGQrwWPz0KFDdHZ2snz5cpYuXUpTU1PM0SdNTU00NDRgMpnw+XzU19en1R8kRbOElCJRK0K24hfWy5cvp7S0tN9efujQIYqKigCicoSmsxiXFM0SMoIU3BKylcAIFL+93GKxsG3bNiwWS9QJSdngDxJBLgjCBYdfGzebzQPs5e3t7QCUlpYC0QnlbPAHSa0VQRAuKIJjyOfOncv69ev7hfj06dNxu90DhLPfEaq3s8yGA1jERi4IwgVDKHv2kiVL8Hg82O12Wlpa+gV9Z2cnAEVFRRHt5enwB4mNXBCEC55Q9myPx4M/v6W4uJglS5Zw7733UlRURG1tbVT28kz6g8S0ImQ1EvUiJJNwtVyCTS433HADQNaWrg1EBLmQcuIVxlKrRUg2oezZwKC0/dWrVwOkpYBXooggF1JKvMJYarUIqUKvlktbW5turZQFCxawZs2ajDkxo0UEuZAyEhHGmS5CJAxtgqsYhjK5TJ06lalTp2a9eU+cnUJMxJK6nEiiRDbE5goXDuFKSgQ7MTN1eEQ4RCMXoiZWM0kihwRkQ2yucGERTfncbPXbSBy5EBXx1pNIdOBL1IqQLiKNtXTWVAlFqDhy0cgziJGEVLw260QPCQh1IosgJJNoFI5s9tuIIM8Q2bpFC4WemQTCpy77EWEsZDPROuXTeZ5orIizMwMEDpxoK6xlmkBn0MGDB9m1axdnzpwZUMNZEIxItE75bK6xL4I8A8QTzZENnvLx48dTX1+P2+1mz549NDU1kZeXZ4iFSBBCEUuElN9UuGTJEh566KGs2UWLaSUDxLpFyxYzjMvloqGhgZKSEux2O4WFhWzcuJEbbrihfyHKBu1EEGIh1gipbDQViiDPALEMnHiTalLhSPXvJKqqqsjNzUXTNLxeb/8zmM1m2traDOG8FYRA4nHKZ1OwggjyDBHtwInHU54qDd6/k3C73cyePZt169bhdDpxu91cffXVLFu2LOO7BkGIl1g07WzZJfsRG3kGiabsZawZjql0pAY6e3p7e5k1axZLly7lgQceYP369YZy3gpCvGRjsIJo5FlOrPa7VMe6RltwqLW1lYMHDzJ27NiMbzuF7CZVJopUXTcb48lFkBuAWOx36Yh1jVRw6NChQ3zwwQcA2Gy2mCoeZovNUUgPqTJRpNL0kY3x5JKiPwTJhP3Of0+n08kHH3zAFVdcQW1tbdpS+f3IYmAcUpXyHu91XR4vH7d109TWxeFT3Zzq6uVUl5tTXb04XH30+Xx4vBoerw98fTg7O8jBi0VpTBk/mok1FYwotjKi2Mr4ikJGldhQSsX9HHpIiv4FRKJp8Ync8+DBgwDU1tYC0W07k1V7PNscUEJ4UmWiiOa6Pp/GvlYHmw+dYfOhdrY1d3C4vYfY9Npz10eDln1nWbPv7IBPC/NyGV9RyOQRdmbUDOOyMSVcVFaIyZRc4Q4iyIcsmYh1tVqtjB07FpvNFtO2MxkTWg6iMB6pMlGEui5mG6t3HOOPu1t5c28rHT2eJD2JPl29fWxr7mBbcwerNjYDUGTNZcboYcwYXcL8qSMYX5Ecc4xErQhJJZ405mTUHk+k9nkyyYYMXD0y1a5w901Gyrve9QeUkzjSwvYzuewsvYpP/uAdvrFyK7/fejTlQjwUna4+3t7Xxn+9sZ9dx5M3NkUjF5JOrKadwMic1tZWfD4fixYtimlC++PbOzo6KCkpyYgDKltNO5lqVzT3TcQMGO76Tms5josX8IfGk3T1eqGlO6nPlgyG+84CI5NyLXF2CllDY2MjK1aswGQyRR3t4nK52LFjB6+88grt7e3s3r2byZMnU1lZmVZBmg21qv3tCBSKmWpXrPeN1Umtd/2znQ5mLPgKz21sYdOhM3G1e2SxlXEVhYyvKGRksY0yu4XhBXkMy7dgzlX4PB5++YunsRUUkmstoLOnl47uXuZe9zlO9fRx/KyLfcc72Huikz5yQt6nwNTHbXnbefjh2PpBnJ1CRDIZ8eGv41JVVRW1jbupqYlnnnmGd955B5vNxrx58xgzZgynT59myZIllJSUpK392RBbrKehFhcXZ6RdsbyP4HbX19dTUVExYDEKHpeB1/dqsNdbzrvO0fzqhY+ibmO+JYfLxwzj8jHDmDmmlGk1xdit5oi/u+9LXzjX3s6TFJvNfG3RQIXB5XLx/e8vhYJhdOUUc6gLDneZaFd2+s6L3Gqrh74k1icSQT5ESFQIZ9osEKsg9Ds3LRYLBQUFAwp4+a+VTjIdWxzK2btkyZKMtCva9xHc7kOHDrF48WLmzJmDzWZj7ty5rF+/ftC4tNvtmHItbGjLZbNzOJ190YmyEcVWrp1cwXWTK7niouFYzaG15lBEMgdZrVbuuOOcqdDqaeVSs5l/WVTPb19swG0bTrupmLy+bswkrx9EkA8BknGcWjIiPvxmunhiZ2O1cUcq4JXu5IxMnzEaaiH0eDzU19cPMlmlul3Rvo/Advf09LBlyxaUUpSVleH1enn88ce58cYbqaqq6h+XDz74IGv3trMmZzYnHX0R25JvyeGzl46g/rJq5owtTUr4X3BUWLAipSfs8/LyWLlyJeWeE0kfHyLIDU4yhHA02rCmaRw63cP25g4OtHVxoK2Lox0uTjl6ae924/b68Po0lIJ8cw75eblU2PMYUWxjdGk+k6oKubiqiItH2MnLHagF+Reirq4uNm3aNMDGHeoZwhXwuvPOOzMScpiJ+H0/oTTg1tZWGhoaMJlM+Hw+6uvr07bTiuZ9+Nt96NAhtmzZwr59+8jLy+v/jcfjISfn3HgpKChgzyk3X/jZBzSeiOy8vLjKzl1za1kwbST5llxcLhenT59Ket+EUqSChX0qx4cIcoOTDNtsKCHQqyy8sKmZP+1pZeOhdtq73RGvpWnQ7fbS7fbS5uil8VjngM+tZhMzx5TyiXHD+cS44Uwqs/YvRFVVVVHbuAM1Po/Hw6xZs5g/fz7Tpk2LOXwtmRMrU7Wq9TTg+vp6Ghoa+t9td3c3DQ0NjBs3Lm01TSK9D6vVSn19PYsXLwbAYrEwatQoPvzwQ+bMmYPZbMbr9eLoM7HmZCH7vaMgghC/ZlI5f3PlRXxy3PD+3WG8u9ZoDmSORZFK1fgQQW5wkmGbDRQCJ06f5aipHMfwefzPD9fFmOkWGZfHx7qmU6xrOgVAWYGZMtdIZpfkkq+5KSkpidrGnaiGk2m/QLIJfh+pdsAGRgwBEd9hKKFYUVHBnDlzKCsro7Ozk61bt3LmzBk6Ozt59B/+gV+8vZ8t7eV4Ioir6y+p5JvXTmDKqOJB941n12qkA5lFkBuQ4AmRFNtsURXOS27mpW3Hcbl9cNyZmsYHcarbwylGsOcYFOf2UWc9w5hcW9QLUbwaTrZkgqZ6R5AqR6dexFBBQUHIdxhOKNrtdmw2GyaTiVGjRlFQUMDp06e5+Svf4N/fOMxHntFh23LNpHIe+MykQQLcTzzC1mgHMosgzxKindChJkS8mumeE508sXYfa3edTMZjJMTZvlze6ypnk6mC3pf28JW5tUytLknJvbJBk0r1jiCZDtjA8QmEjBjSO/IvklAMbie5eXim3MKXVuwIuyO8uLKQf1pQx9zxZWHbHo+wjXZ8RPOO0xHWK4I8AZLVQdFO6GgmRLTtPd7l5T9f38crO45FbT7JyzUxrbqES6uLmVBRSG1ZARX2PIYX5pFvySHXpPBp0O3uo9Pp4cRZF0c7nOw/2cWeE518eKQjKju7x6fxu61H+d3Wo1w1sZxvXjuey8eURtdIouuXTGtS6doRJMPBFjw+/QI72oihaISiv53r9hzjO2sPcbgxdEJPvsnLNHWI/1x4O1WV4YU4xLegxTI+wr3jdJnvEhLkSqm/BP4FmAzM1jTtgknXTGbZ1WgndCJapL+93W4fOxnNbm8lvigEeJFyUa3aGWN2sGThAiZPmhD2+zkKiqxmiqxmqoflE5iC5q849/6B07y7/xTv7m/D4w3fiHf2tfHOvjY+OW44f3/9RGbWhhfo0fZLLJM7leefpmNHkIiDTW98rl69GiDqiKFohKLL4+WJPx3k6Xc/DqtYzLB3M9t6HE/3WUqKi6J+huLiYpYsWYLH44m5bES8BzKn03yXqEa+E/g88LMktMUwJLODYpnQ8WqRLpeL559fyQnLCN7traTLGz4JojAvl5unVuLY8QZji3MoLCygu1vx29+sSmgQmkzqXAhiVRF3zR3L2R4Pr+06weodx1m3vy3swvLegdO8d+B9bqyr4uGbLmZsWYHuc8bSL9Foq6k+/zTTtlU/oRarUONzwYIFrFmzJqqIoUhCcUdLB/f/ZhsH2kJHo4wdZqGuZwfDXB14vNGbiPT6r7y8PKp3kuhuxuFw4HQ6+wu4pXKxTkiQa5q2G+JLADEy8WpTepMllgkdr82z6dhpXnON45izJOz3ygtyue+aidw2sxqXo4Nl285SWFgT0zPGQp7JyzVjrNw8ZSptPT6e/+AIv9l0hDNhKtOtaTzBG7tP8qUrxvB3102gJN/S/1k8/RJOW02lRpXpBKJAIjki9cbn1KlTmTp1atRCTk8o+nwaP3vnY364di99IVbxvFwTf3/9RO6eNxav58qY67Ek2n+J7GZaW1v54IMPyMnJwWazMWXKFKxWa0oWa7GRx0E82lS4pIFYJnSsWsLqHcf4h9/tpFMrCfmd3L4erig6y0W+49wx81NYrWbMnHuWkydPUlpaitvtTqrGqPc+HrnpYv7uugm8tO0oP37rAIdP9+j+ts+n8av3DvHy9mM8ctPFfOGyakwmlTQt17/gulyutJ9/mm5idUQGj89Y2hwoFE92urj/N9t478DpkN+fXlPCD2+bxrjyQgDMOX/+fbDzNZbdRDoc2v7aQVdccQU7d+7E6XSyYcMGnnzyyczEkSul3gCqdD56TNO0l6K9kVLqXuBegNGjw4cTZTuxCt9IkyWesq+RvuNwefjnlxr53dajob/U56bG0Ujh0U3MuOYqOjrc/YO8paWFzs5Otm3bBsD06dNZvHhx0uKPA99HR0cHTz/9NI8++iglJSX81azR1F9Wzcvbj/Hff2ri41P6W+72bjcPvbiDFzY1851bpzB5RFHCWm7gAgPQ2dmZMvNHNhxLF4sjMta2hnq+13ed5KEXt4fceVlyTNx//UTuuXIsuTmDj0wI7KPOznMJZ0VFRVHtJvztcrlcKX3n/vdaW1vLqFGjcLlcnDp1ioqKipTcL6Ig1zTtumTcSNO05cByOFfGNhnXzCShBnekSm2gP1mSmfG1+3gnX312C0fa9TVagOKuw/S8+wyFI0v7s+f8QsovaGtra5k4cSLt7e243W6qq6vjblPgewm0HR49erQ/AeR73/se99xzD+PHjyc3x8TnL6vmlumjaNjSwg9f38vJzl7da28+fIb5T67jK5+s5e+vnzigXwDa2tqiEkB6C25nZyft7e1JN39kIhkpEdNerONT7/mqx4zl8Vd38+yGwyF/N2VUET/8y+lMqtJfLAP7yGKxsGnTJgAWLFiA2+0Ou5vwC/3ly5en/J0Hv1e3243NFn1+RKyIaSUBggd3qMmZTsfWS9uO8nDDDlwen+7nBaY+Pmk5QvO2Vzlt8dHT08OUKVNwuVz9QqqtrW3AwlNZWdl/gksywtdmzZrFBx98AMDRo0eprq5m2LBhDB8+fJANM8ekuG1WDfOnjeDn7x7kZ28foNvtHXQPr0/jF+sO8ocdR/nurXVcfcnImIWlnnOqqKiIe++9t9+2mYodSTqSkZJl2osGved76rnfsbNoNvtbQzs0v3rVRTzwmUlYckMfXNbW1sbZs2cZPnw4LpcLk8mEUgqXy4Xdbg+5m2hra2P58uWUlpYOeuegb5pJZMeUbh9IouGHfwE8CZQDryqltmmadkNSWmYwwpUR9Vega2hoSFmnerw+vvt/u/mf9YdCfufK2kJqTq5Dc3VRVFfHrbfeyqWXXjooJCuZC4+eGeUHP/gBM2fOZPv27bjdbo4ePcrtt99OSUlJyAUj35LLkmsncPusGv7t1d28vP2Y7v2Odbr5yoqt3DRhP8Nb3mV4kb6w1JukoZxT5eXlSe2rdNtuk23ai0Tg82ka7POWs9ZVjNelL8RLbTksWziTeRPCx4T7s0k3b95MY2Mjc+bMwec7p7BYrdaodhPB73zHjh390TeBC1wydkzp9IEkGrXye+D3SWqLodGbnE1NTXzve9/DYrHoFswPRayawFmnh68/v4X1TfqOowJLDv9ycx1fuLyaXbvKWLFiBRaLhTfffJMRI0YMGqDJ1CaC30tOTg4ej4fKykpuuukmAPr6+rDb7VEtGBVFVpZ9cQZ/NauGf3ppJx+HCFn7w/4u8rXJXG9pZaJlYOhXS0vL4G1/dXXanFPpDj1Mt2nP/3ztDidv94xkd3d+yO/WWrr4BAeYWXP1oM/0sklLS0u5/vrrWbduHW+99RZTpkwhNzeXEydO9M8xh8PR/0zBbQq2l7/yyiuDtPQlS5YkbceUzPcaDjGtJInggdLR0cHu3bu58cYb++trNzQ0JKVQTyAtZ3q4+1eb2HeyS/fzSZV2fvrlyxlbVhDTKTzJ0Cb8TiWg/734bfFerxe73c7ll1/Ohg0bOHXqVEy1sueOL+MP37ySn797kGV/3Edv32C3S4+y8tLZ0ZQe2s2nS09QZDNjNpt1J+ldd92VNudUurfdyV44IikaVquVy6//PI+8vA+Hlqd7jRyl8enSs8ws6qalxTloNxIqm7SgoICCggIWLFjAgQMHuP/++ykvL8fhcPSX7NWbO3rvfP78+axevXrQAnf8+PGYDznJtNNaBHmSCB4obrebyZMn95dijadQT3A0RzAftZzl7mc20ebQdwKOzWnnP26c0J88E+uWPhFtIjiyoLOzk6Kic5l4X//619m+fXv/fZ988smodirB5OXm8LfXjOe6iaUsemotJ336gql92GR+7+rk29dOwOPx6L4DGFhgKtXOqerqau666y6ApJtugknmwhFJ0dC0c76K7685gieEEC9RTj477CS1xWb9LM8w2aSB/VNSUjLg3f3kJz8Jq0UHKycAr7322qAFbsSIEVEvfLGU10ilsBdBHoJ4XnzgQDGbzSxbtizuQj2tra1s3LhxUDSHn7f3tXHfs1twegY7/hQany7tpC73FL//7a+pOW+nN5vNCWlm0b4TvYnY3t7OZz7zGdauXcv27dsB4qofrsekUaU8e9flfGvFG2z2jNI99NabV8Q/vnGSxo5ccnPzBr2D8vLytGnJkc6oTAXJ2GH5FYvhw4cPOLHH7zBsbm3nu2808+a+UyGvMTGnjb/7VDWbN+yiudmj+54jZZPq9U8sRa4C/63X5yUlJVGNhWid1umIUBJBrkMiLz5woMRbqKejo4ONGzdiMpl0oznW7DzO4lVbdeuUKG8vtS1/ZGzhCApLKjhwYKCd3n8GYqzCKpZ3ojepWltbWb169QCzzmuvvca0adOiea0RmThhAr96rIY9Lad4fO0hNh05q/u9VZuPMsI+mxkdOxne3jzgHaTDORXpjMpYNbpYFI5Ed1hPP/007733HsOGDWP27NlUVFT0Owx/9eq7vNU7lh70tXCL8nJT+Vlq1Rk2bzgStu5JPNmk8ZqPQvV5NGMhmsUjXRFKIsiDSOaLjyfRZ+HChTz99NOcOXOmf8IERnP8365TPPjiDrw6Kc0FysPIQy9T4Olg48Zm5syZM8hOv379+piKB8XzTvQmlc/nw2QyJRypEek0munjq/nNRaP40atbeGr9MV3t/LjDwwkm8cXLK3ngxsmU2AsGXCOc6StRIR84+T0eDzt37iQnJ4eysjJMJlNMGl2ov4drZzzP4O//4cOHM2zYMEwmExs3bmTevHloGnz3f7ew3TcJDf1SHUXesyysdVJq9gJ/Pks0VM2TSKagUGbAeM8mDdXnkRa+aBaPdEUoiSAPItkvPlYtaPz48Tz66KN873vfY/jw4QMOIl695yzffnWv7u+G0c0dNZ24Csb3m2ROnDiha6cPN4n0iMe2HjwRFy1aRENDQ0IOt+Csy+uvv56xY8cOsjGbTIr7F8xkwcwzPPK7j9jc7Bh0LQ1YueUkG4508cO/nMaM0cOivnc0u7RQAjNw8vt8PpxOJzabDavVei7So72dtra2/rh1IGRYq97f/WGueu2Md6fp7/+qqipmz57dP76Onu7kQPk8dp4abN7zMzXvNI51z2Iqvx6iOFTbT6xKUFNTU9rPJo3G95CuCCUR5EGkOzRMj5KSEu655x5WrlxJc/O57X9O3Q0hhfjlo4uZ3LaVnN4CKioqmDdvHqdPn+b+++9n+fLlEZ8lkpYWzzsJd4p4PDbowF2B3yzz7LPPUlNTw2WXXcbixYsHTdzxI4bx3c+M4FsrGtnkGYVXRzv/uK2b+p+8x32fGsd980bT6+zWNV/EsiMJJzADJ7/T6cTr9TJlypT+99vZ2cny5csBBkVrQPjIitbWVlasWKEblQT6C0KoZwgcE4H97x9fH5328bZtGh0hhLjV5OVz5R1MKnCx/eKJnD59ut93lKimrNdW/7Ol8mxSPSItOOmKUBJBHkS6Q8NCEThAnt96mv/60wHd7105oYyffflyjh2pGtDme+65h6qqqojP4hc6TqcTn8/HokWLqKurG3CPeN9J8ESMVssKXlhcLhd79uzh1KlTFBUV8f7779Pe3o7NZqO0tJS9e/fyzDPP8Nhjjw0SwL9etYo5pXamW9p46UQxxz22QffzafDjtw6w6u2P+GROE5UWzwDhG8uOJBqhH/ge/CFzzc3N/dcIjGsOjtYIF1kRznwFRP0MeguRv/9bT3ewlVr2qApw6WcPl3rPcEvZGUYUnHMqV1ZWxmzOC0fw+MhkcSyIvOCkw/ciglyHbKhKB+cGyIqNR0MK8RvqKln2xRnk5ebE5bTxCx2Xy0VjY2N/Esy//du/MWvWrLiEcDTPFO63wUJk7ty5vPzyy3z44Yc0NzdTVVWFyWRC0zQsFgv5+fn09fXR09MzaOIGTvACvNxZ3c7aw33sVLW4dRzFZ7R8/tB3KbPzTvPc8yt55OGH+k0c0e5IYo2eKC8vZ9y4cf3VFpcvXx5VtIZeZEUk81U0zxBuIfrsl+7jgd9+xOEzrpD996XLyvnS9Im88OtVNJ/1DGhvMtCt4VJdnfFddCQScTRHgwjyECT7xcfjZPqf9Qf57v/t0f3s8zNGsfQLUwdUh4vVaeOvLdLY2Ehe3rlogwMHDvDwww9z7bXXcueddw4wV6R6MOrF0f/rv/4reXl5VFVVUVhYyN69e+nr6yM3N5cxY8bg9Xrxer3k5+cPmrjBArjzbAcTvKd57J5b+faaj9nRMjiyxYdiQ3cZ+5SVGw+e5IrJY2LakcRjhvK/V5fLFVO0Rqzmq2ieQW8hOnW6nWVv7GX5e0dD1g23mfq4yvwx/3jztVit1pQoQuEWmWzYRWcSEeRpIB4n03MbDvPtV3bpfjYpt437LpuoW+IzFux2e7/DraCggF27dmGz2SguLsZisaS8kFMweun8vb29mM1m8vLyyMvLY+LEiQwbNgybzcaBAwfo6elh+vTpukeMBQrgpqYmdu/ezeTJk3l15c9ZevsXeb25kh/9cb+ucGrXCvnyc4088Jk+7rnyoqh3JImY5uKN1oh25xTNMwQvRC1n3bzhq+PkupaQ7R6pOrgq7zB/c8dtA9qa7HETbreTLbvoTKG0aE/eTSIzZ87UNm8eOsd7Rgr3Wrp06QAtwuFw6ApI/3XW7nfw2Mu7de813d7NVbajdHXpXyOWtgE0NjayePFivF4vx48fZ9y4cf1OthMnTrBkyZKYIlwSIfhddXR08Oqrr5KXl0dxcTFKKTo7O5k1axYPPPBAv+03UmZkR0cH3/nOdygoKKCyshK3293fBwfae3nghe3sOTE4ssXPZaNL+H79VCZURr9VTyRUMRlhjolco6mpieeeX8mO3jK2+2rwoq8w5OWaePD68dw8uZiioqKUC89Y5tJQRSm1RdO0mcF/F408QSJp29HaTP3X2d9r513vBNCJyb20sIebyjpQqoAzZyI7c6LZCdTV1fHkk0/yy1/+ko0bN6JpGrNnz454IlA0J7TEip5G+q1vfYuXX355wAEXd955JyUlJVHbXffu3cvmzZspLCwkJyeH2bNn4/F4cDgc1I0s5+VvzOOJ13bz9LpD6J0F/eGRDm760bv8zZUXseTa8eRbIk+bRDTSRLXZRDMJe/PL2VhyNbuOh17cLhlRxI9unz5gcUt1Gnq2BCJkIyLIEyCaCIVobKb+67RbKlnfMwo9IX5RTjtX559CqehOOoklZK6uro7HH3+cHTt28Morr9Db24vP5ws5SaI9oSUe9LbIc+fOpa2tDYi9LonL5WL16tXYbDby8/NRSrFu3TpmzZrV3wdHDn2Me0sDnzFZeM83nrPa4MiWPp/GT98+wMvbjvKtBXXcUFeZlWfVxtLvfsFrNpvxeDzk5OXz8/ea+dk7H+smnAEoBV+9ahx/f/3EAXXDY1k8EhH4F7oJJRQiyBMg2vKgkbQIh8PB8V4LbzhH4tMR4jdNqeLvZk/kN79eRceZ6E46iSeJZ/bs2f1ONf/kDl4oYjmhJRTRVM8LjvCoqakJeb1w+M0v8+bNY+PGjfT19eF0OlmwYEG/g/HPMcgFXNTVwtvtdhr7qtATZcfOurjvuS1cM6mcB68bR7nVl1DafCji3fHEugM8efIku3bvpuTST7Mvv45uzRLq0tSU2vjBF6ZxxUXDB7U12sUjGXVHUu10NyIiyBMgnLYdOBEjaRFtLhN/8k7Go2OLvPbiCn50+wwsuaaIJ50EXjfexCar1apbr1svptrhcEQ8oSWYdB9x5n8PBQUF3HDDDf2VKadOnTroeQCKC/OZceYQ3/yra3j8jcMcCnEA9Jt723hrbysTTCe5PO8kf33HbUk7kCCRHU8sO0Cr1cq+th5ap3yRg7aR6K5cgEnBXXPH8sBnJuqalaJdPDJxMtKFQmJhDxc4fm3b4XD010JZuHAhLS0tLF26lGXLlrF06VKampr6Y4aDB+zxs07+5rlt9OqsqbNqh/HUHZf1b2EDNZHASeO390bTtmico/7JVlNTg91u7481h4GCwmq14vP58Hq9WK2hT2iJ9tqpIPA9nDhxAk3TBkS4BD4P0P8MV04eyY8XjGKG+Rgm9BNfNBT7fFW86LqUB595kyMnTiX8fMFZinv37mXv3r1UVVVFdb3A5z148CAnTpygvr5+0A7wtDuHtd017L/oC3Tnjwx5vYmVhTR87ZP80/xLQvoGQr3D4HGgJ/D1xq4QO6KRJ4hejeNgz3ooreNMt5sv/2Ijx84OnpgXV9n5+Z2zsJoHppXHomnHY0+MpF0Fmoo8Hg+TJk0C6D+hJdxikakMvHDvIZTpC+DFF37NVaV25lha+cNJO4fdBbrX92gmtnpGctNTmxjnreTachOgxfV8ie54/M8bWECqoaGBvLw8xo8fz6FT3Tz1Zgsv901D6wtt47fkmPja1eP422vGhz1DE6J3QmZD+YuhigjyJBCoKQcfXBxqMve4+7j7mU00tQ4+2aem1MaKu2dTbDPr3isWz32s9sRoJpve4hXNYpHJiRzuPegJ+sB+LMDLwlEdrD9ylD22KbR2uXWv0+3xsYNq9hzxMs3uZLK5Lebn09vx+Nsf7fvSOwnqv579X9wXXcXre9o4F3EcWohfe3EF/zT/EmrL9BcuPaJRGqzWwRUKQx3NJsSGCPIkE42w8vo0lqzaytYjHYN+X1Zo4dm751BRlJnaDdEuFHoOyWRdOxMEP09wP/b0dDM+z8Gyb36S5zcf56k/7afHo29ycWs5bOosZBOFzKqewrqPz3LVREtEzdbfjnh3PH78Wn2urZBtnTa2O8o52muB3W1hfzfSnsvj9TO45uL4jreLpDQEVyicNWtWyEqNQmxIQlAKCOfw0jSNf365kRXvHx70OzNelt16ETddMSUpBwkkQirvkw1nHEZDuH481dXLE6/t4TdbjoYM1QukJN/MTVNGcPO0kcyqHRYxKzfeqJUedx9v7jrOf/72TQ55h9EXhRusIMfLrPxTjOk7yqMPP5iSPtFL9lqzZs2AWvkXWnJPPEhCUBoJpzH//N2DukI8B42bhx1j09ptjLYrXU0lWTWxoyGVIV5GCR8L149lhXl8t34aX7tmAj95+wC/3dyse2KTn44eD6s2HmHVxiPYrbl84qLhXDmhjCsuGs5F5YXkmAaaOvR2PC6Xi7a2tgFt6e3zsutYJ9uaO3hnXxvvHThNb58PGBgiqIeZPuYO62ZmcTcWk0Zzsztl/gq98gsej4ecnHM+oHRXLBxqiCBPEXrC6tUdx3n8//RS7zVurmhnUqGJgwedujWlQx0kEE9NbCF6Ii06NaX5fPcvLuUb14znp28f4IXNzbhCmFz8OFx9rN11krW7TgLnUt0nVtq5uMrOqGE2ygrzKCvMY1i+mdwchVKKjw81s+aNN3F4TTiVjbIxEznda2L3cQdub/j7BVNZlMeiOTUcX/dbhlsKsJhChykma+cUbKryer2YzWa83nP1zMXxmRgiyNPE5kPt3P/CNt3PrixsY3KhJ2xNab2DBBKpiW0EjGKCARhZYuNfb5nCA9dP4oXNzazYcIjmdmdUv+3t8/HR0bN8dFT/nNE/E7AQH4g9ZO/iKjt/PW8st0wfhSXXRNOY20P6K5KtCOj5Rx577DHWr1/ff3hKtvhLjIgI8hQQLIA+buvinhWbcfcN1pxuvaSEYQe39J8oHqqmtN5BAonWxM5mjLqjKM43c89VF3H3vLG8va+VFzcd4Y97T503d6SfImsut84YxW0za6gbWTSgrEAo01GqFIFQ5ReMslhnMyLIk0ywAPrsX9zG3//hOGd6PIO+e/0llfzwS5fjcV8esaa03kECyayJnSySlZ5u9B1FjkkxOtdB1cE/8Bd4OW4pwzNiKjtOuOjq7UvpvcsKLVwzqYJrJ1dy9aTyQbkIgeiZjlKpCOjZ/o3Sp9mMCPIkEiyAznb18LVV22nzFQ767rSaEpbdPoMckyLHGl1N6XTUxE6EZGnRQ2FHMSBDs6CAyu5uHG1v87OH/h/7TrlYt/80Gz4+ze4TnXToLPKxUG7PY+pIO5PKrVw1qYrZ4yowmeIv6GU2m3G73XR0dAw4/NuI9msjmecSQQR5EgkUQD4NXu+qps03uJJeTamNX9w5E5slNk0p3N+DSXeVuGRq0UMhA1BvMWptbaXlyGHqxo7l8jGlfJMJaJrGyc5edp/o5OO2bk519XLK0cuprl4crj68mobPp5FnzkFzO2ltOUSh6qXE4uPmT8/jyukT6W5rYdWqVTg+9rD2QzNlCZih/ItxV1cXmzZtYvLkyVRWVhrSfm1U81w8iCBPIoEC6D3nSPb1DBbiJflmfnXXbMoK81LennRuW5OpRadrR5FKbS1wLFgsFvbv38/WrVsBsNls/UJFKUVVsZWqYivXTArf1qVLl2Kv8i+ULva+8xI3TF/CT1aswGKxUFVVFXUFylD3CKzzMmbMGE6fPs2SJUtiOnMzG7TgoWCeiwUR5EnEL4D+8ZnX2eQZbE6x5JhY/uWZjCsf/JnRSbYWncodhcvl6q+9DqREW/OPhSeffLL/4OjRo0efS/kvKIhZqIRaKD/44APeeecd3UMzYn1nwfcoKSnp/1u0ZIsWPBTMc7Eg1Q9jwJ+QEa763EG3nY19o3U/++Ft05g9tjRVzcso8VZbjHTNWA+SiERTUxOPP/44Dz74IJs2bSIvLy9lVRirq6spKipi7ty5XHzxxVRVVbFx40YsFkvMVf/0KgwCvPXWW/2HZlitVtatW9f//ViJpophuDmQieqWoYi2IuNQQTTyKIlG09jR0sGSVVvRy9h+5KaLWTDtXLnQZG89s2ErC9l/eotf0FgsFgoKCigsLGTjxo3ccMMNcWux4fAL6pqaGnbt2oWmaXi93n5zUSxCRc/cNH/+fFavXh3y0IxYiWTSStaxhukgm+v6pAIR5FEQjb2t5UwPd/9qM06Pd9Dv75gzmq9edRGQ/K1nY2PjgGpymXboZHM4mV/QVFVVkZubm5BgjQa/Vuh2u5k9ezbr1q3D6XTidrsH1ESPFr2qk6+99lrIQzPiIZHY8mxzUme7YpFMxLQSBZEK4p91erjrfzZxqqt30G+vnlTOt2+u668pncytZ2NjI4sXL+bDDz+ksbGx//qZ2MoagWDB6nA46Orqwu12p0RbCzQ39fb2MmvWLJYuXcpjjz0W92IbaG4KvL7eoRmJtDvYpBXNoRCpMK/Fi98EBLGf82pERCOPgnCahrvPx9ee28J+nbril4wo4r8XXtZf6S6ZW0+Xy8WKFSvIycmhsrISt9vNzp07qaurG7IOnUQJLhE7a9Ys5s+fz7Rp01L2vlKtFaZL64xW284GLThbHK7pRAR5FISyt+Xl5fHgizt478DpQb+pKrLyy6/MojDvz684mVtP/+kxNpuN3t5e8vLyOHPmDD6fb8g6dPwY7RT2VJub0mHOisXmnEnz2oUWduhHBHkYIh2gvOyP+3lxS8ug3xXm5fI/d82iqnjwYQzBJ6TEu/W02+3YbDamTJnCzp076ejowOv1smjRoiE9YIfyKezZ4rQORTZo25HIJodrOhFBHoJQAsM/GH6/tYUnXt836Hc5JsVTd1zG5BFFutcMPCGlvr4+IVupX0Oqq6vD5/OxaNEi6urq4rqeERjK2pZRzAHZugj6SZXDNdsX2YQEuVLqB8ACwA0cAO7SNK0jCe3KKJEExvsHTvPQizt0f/udW6bwqYnlYa9ZVVVFd3c3DQ0NjBs3LuzACDeAjKAhJZNs0rbindh6vxvKC1SixPqeUxF2aIRFNlGN/HXgUU3T+pRS3wceBR5OvFmZJZzAaOns46vPbtY9Deb2aaUsnKOfDBSPEIpmAGW7hpRMsiW8Lbhf6uvrqaioiChsQvVnNi1Q2US8AjSZCo5RFtmEwg81TVuraZq/JucGoDrxJmUWl8vVH74XnBXWo1m485cb6XQNLkM6Ma8T277XQ4b+xZpplk1ZctlCNoS3BfeLy+Vi8eLFPPHEEyxdupSmpqaofhfYnxdaFmI0JDr+k5UVHE3YZTaQTBv53cBvQn2olLoXuBdg9Gh9rTXTBGoAnZ2ddHZ2UlRUhNls5uYv/BVfXbmdox2DT32pzuvl1hEOjh8NnR0Y65ZPtDR9Mm1OCuwXj8fDzp07ycnJoaysDJPJFFJbC9ef5eXlMZsDst1mmyjZMv6zZRcYiYiCXCn1BlCl89Fjmqa9dP47jwF9wPOhrqNp2nJgOcDMmTMjHzueBoJPKg/eQrW3t3PvvfdSVFLKV1ftYM+JwatwcY6bL1S10+uM3MF6QijUhDTKAMoEmTQnBfaLz+fD6XRis9mwWq2YzeaQwiZSf8ayQBnBZpso2TL+jZLqH1GQa5p2XbjPlVJ3AvOBazVNywoBHQ3Bk8FfbyNYAzBb8nj05T1s+Lh90DVKrDlcyz5OH+uKuoMDhVC4CWmUAXShEdgvTqcTr9fLlClT+oVOKGETTX9Gs0Cl22abKc0/m8Z/pneB0aASkb1KqRuBJ4BPaZrWFu3vZs6cqW3evDnu+yZKf23nIO0boLS0tP9vnZ0O3JfewvMbB8eK5+Llv24Zy3UzJsQdvRDcBofDMWhCDvUttFHx90traysNDQ1Ra8eJ9mdbWxvLli2jpqam/2/Nzc0sWbKE8vLB0VKJkA2av4z/gSiltmiaNjP474nayP8byANeP3+o6wZN0+5L8JopJ5T9bcGCBaxZs6ZfA1B1N+gKcRManys5zoevb+O6GQ/FNYGitQFeSFEpRsLfL+Xl5YwbNy5qYZNof6bL5JAt0Roy/qMjIUGuaZohDXOhJsPUqVOZOnUqDoeDtU0OHntpt+7vP1d+hjq7ork5/tKn2WIDTBZDTXOK5XnSKWzSZXLIFmejEB0XZGZnpMnwp/1n+KeX9YX4vMI2LrW7Exa82WQDTJRs2IInk2x/nnTYbIeaojHUSchGHi/psJFHo1EFR604HA62nujlb3+9Qzfh5/NThlF84A36+pI3wY2uyUZr6zcKmX6ebBoP2b6gXYikykaeFQQP/mgHoH9L7P9+S6+VP3on49XJk1owbST/8VfTcbsv051o8U5Ao9sAh9oWPJPPk22C0wjRGsI5DC/Ig0/Iqa+vp6GhIWonjd+p020t5y1nta4Q/9TEcv7jL6diMildwZttEzCdDLUteKaeJ1uci8EYXdG4UDD0CUF6J+SsWLECp9MZMaXWf4JIS0sLR8728VLHKNza4Ncxe2wpP/3S5eTl5ui2IVwqcTSHNRudbEibTyahngeI2JeJ9LdRUsGF7MSwGnmoE3ImTZoEEFaj8mvQJ0+e5MOPT3Ji8m1ovsGCelp1Mb+4cyY2S05I00morfiOHTtYs2bNBaGlD7UtePDztLS0sHTp0rB9meiubKjtbIT0YliNPPiEHIvFgtPpxGQysWjRopAaol+Dtlqt7GvroXXKF9EsBYOuf3GVnWfuno3daqapqYmlS5eybNmyQYWR9AoeAbzyyisXVMGrZBUpyhb8zwOE3HH5SUaBs6G2sxHSi2E18nAn5NTV1YVM0vBr0F2W4RwYPR9fzuCJMrasgGf/eg4l+ZaItku9MML58+ezevXqIeMAvJCJ5Px0uVwcPHgQp9NJVVWV7neiZajtbIT0YVhBHihA9U7ICeWksdvtdJiKee1sDT4du3fNMBvP/c0cyu15QHRRDMETEOC1116TbfIQIJzJw29OcTqdfPDBByilqK2tTai/Y3EuZlOoopBZDCvIIT4NZueJHv7oq8Ot+QZ9NqrIzG+++glGltj6/xat7TJ4Ag6VZJ8LnVCJW/Bnk0tVVRVKKTZs2ICmaQmdxRotF3KklDAYQyUEJaqB/HH3Sb7+/If09g0W4rXD8/n1vZ8YdGAyxD9pRGMaOgT3pV7xqoMHD7Jo0SLGjh07qL+TORYynbQkZA7DJwQlqoH8dnMzj/zuI7y+wQvXuPICVt1zBRVF+pMgXtulxOAOHYL7Um+nZrPZdIV4srXnoZaEJSSOIaJWEokK0DSNn759gAdf3KErxCdUFPLMnTNQvY6Q1xPNWggm2iiTVBzZJ0fDCcEYQiOPVwPxeH1866WdrNrYrPv5paOK+eery/jlj38UUlsSW6QQimh2aqnQnodSwTUhORhCkMeTLNHR4+Zrz33I+x+f1v18pKmTn942h1/87MchQwuzNW1ayB4imc9SlegjoYpCIIYwrcSaLHGgrYtbn1ofUohPLujhatNuOtvbwqZFS9q0kCipTPQZaklYQvwYQiOH6DWQV7Yf45GGHXS7vbqfzyzq4pPWY3RpuYwYMSKstiRp00IyiFd7Ft+MEC2GEeQwcBsbPMh7+7w8/upuVrx/WPe3JgVzco8w3nmMrr5zNsWSkpKwtkaxRQrJItYIJvHNCLFgqDhyP8GD/KrPfYEn3m9nR8tZ3e/brbn8+I7LmFVjj6uWuGhGQjpJdZy4jGfjYvg4cj+BDkhbfgHr28w889sDunXE4VzdlKcXzWR8RSGA7sCNpC1JPLgQSKoFYSrjxEXTH5oYTpD7B3mvpYiG48NoceWF/O7nLh3Bv9dfit1qTmMLhaFMOgRhqnwzEoU1dDFE1EogymJjuzaGp5srQgpxc47i2zfX8d8LZ4gQF5JGKpJ79EhVpItEYQ1dDKOR+3wav9t6lO+v2UNb34iQ3xszPJ9lt89gWk1J+honXBCkMzU+FXHiEoU1dDGEIN9y+Az/+koj20M4M/185ZO1PHTjJPIthngswWCkWxAm2zcjUVhDl6yPWjnT7eYT//5HXJ7BFQv91JTaWFo/jU+MG677uXjphWQxFJyFMh+MS6iolawX5ABPrN3Lsj81Dfq7JdfEPVeO5etXj6cgT18LDyz+H3z4hCDEgwhCIVMYOvzwvqvH8cLmFk50/tmpdNOUKv7hs5OpKc0P+Tu/c8rlctHY2IjT6WTTpk08+eSTIsyFuJFwVCHbMETUSr4ll0duuhg4dyjyqnuu4CdfupzyfBNtbW0howYcDgdOp5OdO3eSl5dHZWUlOTk5rFixYkgfhCwIwoWFITRygFumj8ScY+LGKVXkmFRUtkq73Y7P58PpdFJUVITb7cZms2EymaQIvyAIQwZDaOQASik+N3UEOSYVdTyv1Wpl0aJFeL1eTp48icvlYsqUKdhsNgm5EgRhyGAYQR5ILIkNdXV1PPnkk1x22WXU1dX1h2CJNi4Ig3G5XGHNlUJ2YhjTSiCxxvPW1dXx7W9/WyINBCEMQyG08kLFkBp5PCnMUoRfEEKTrvIDQmowpEYOctSVICSTdJYfEJKPYQU5DI7nlUQNQYgPqcNibAwnyEMJa7HvCUL8SB0WY2OIFH0/oYR1qk9UEYQLBdnVZjehUvQTcnYqpb6jlNqhlNqmlFqrlBqZyPXCEc4ZI3WWBSE5SFCAMUk0auUHmqZN1TRtOrAa+FbiTdInnLAOtO8BKbfvSaytIAjZREI2ck3TOgP+WQCkzE4TzhmTTvue2OIFQcg2EraRK6UeBxYBZ4FrNE1rC/G9e4F7AUaPHn354cOHY75XJCGaavue2OIFQcgkcdcjV0q9AVTpfPSYpmkvBXzvUcCqado/R2pMvM5OiE1YJ1uwt7W1sWzZMmpqavr/1tzczJIlSygvL0/4+oIgCOGIux65pmnXRXmPlcCrQERBngjR1oJOhQlEYm0FQchGEo1amRDwz5uBPYk1JzmkKt04VaebC4IgJEKiCUH/rpSaBPiAw8B9iTcpcVKZbiylAQRByDYSjVqpT1ZDkkmqTSBy1JcgCNmEIasfRkJMIIIgXEgYKkU/ViTdWBCEoUTcUStGRkwggiBcCAw504qkzwuCcKExpDRySZ8XBOFCZMho5HJUlSAIFypDRpAHxo57PB58Ph9Op1NK2QqCMOQZMqYVf+z4oUOH2LlzJ06nE6/XS2trq9RBEQRhSDNkNHKr1Up9fT0bNmzA6XRis9m44ooraGhoGGBeEWeoIAhDjSGjkQNUVFQwZ84cysrKsFqtmM3m/oQgq9UqzlBBEIYkQ0Yjh3PmFZvNhslk6k/R96fmizNUEIShypAS5OFS8+VcT0EQhipDyrQCoasTSi1xQRCGKkNKI/ejdxK4FNISBGGoMqSLZukhhbQEQTAqF2TRLD2kkJYgCEONIWlaEQRBuJAQQS4IgmBwRJALgiAYHBHkOkgavyAIRuKCc3ZGQtL4BUEwGqKRByBp/IIgGBER5AFIGr8gCEZEBHkAgWn8gKTxC4JgCESQByBp/IIgGJELLkU/GiSNXxCEbERS9GNA0vgFQTASYloRBEEwOCLIBUEQDI4IckEQBIMjglwQBMHgiCAXBEEwOBkJP1RKtQGH4/x5GXAqic1JFtKu2JB2xYa0KzaytV2QWNvGaJpWHvzHjAjyRFBKbdaLo8w00q7YkHbFhrQrNrK1XZCatolpRRAEweCIIBcEQTA4RhTkyzPdgBBIu2JD2hUb0q7YyNZ2QQraZjgbuSAIgjAQI2rkgiAIQgAiyAVBEAyOoQW5Uur/KaU0pVRZptsCoJT6jlJqh1Jqm1JqrVJqZKbbBKCU+oFSas/5tv1eKVWS6TYBKKX+UinVqJTyKaUyHiqmlLpRKbVXKdWklHok0+0BUEr9UinVqpTamem2BKKUqlFKvamU2n2+D7+Z6TYBKKWsSqmNSqnt59v17Uy3KRClVI5SaqtSanUyr2tYQa6UqgGuB45kui0B/EDTtKmapk0HVgPfynB7/LwOTNE0bSqwD3g0w+3xsxP4PPBOphuilMoBngJuAi4BvqiUuiSzrQLgV8CNmW6EDn3AA5qmTQauAP42S95XL/BpTdOmAdOBG5VSV2S2SQP4JrA72Rc1rCAH/hN4CMgab62maZ0B/ywgS9qmadpaTdP6zv9zA1Cdyfb40TRtt6ZpezPdjvPMBpo0TftY0zQ38Gvglgy3CU3T3gHaM92OYDRNO65p2ofn/9/BOeE0KrOtAu0cXef/aT7/X1bMQ6VUNfA54OfJvrYhBblS6mbgqKZp2zPdlmCUUo8rpZqBO8gejTyQu4E/ZLoRWcgooDng3y1kgWAyAkqpWmAG8EGGmwL0my+2Aa3A65qmZUW7gP/inPLpS/aFs/aEIKXUG0CVzkePAf8AfCa9LTpHuHZpmvaSpmmPAY8ppR4FvgH8cza06/x3HuPclvj5dLQp2nZlCUrnb1mhyWUzSqlCoAH4u6AdacbQNM0LTD/vC/q9UmqKpmkZ9TEopeYDrZqmbVFKXZ3s62etINc07Tq9vyulLgXGAtuVUnDOTPChUmq2pmknMtUuHVYCr5ImQR6pXUqpO4H5wLVaGpMHYnhfmaYFqAn4dzVwLENtMQRKKTPnhPjzmqb9LtPtCUbTtA6l1Fuc8zFk2lk8F7hZKfVZwAoUKaWe0zTtS8m4uOFMK5qmfaRpWoWmabWaptVybgJelg4hHgml1ISAf94M7MlUWwJRSt0IPAzcrGlaT6bbk6VsAiYopcYqpSzA7cDLGW5T1qLOaVG/AHZrmvZEptvjRylV7o/KUkrZgOvIgnmoadqjmqZVn5dZtwN/SpYQBwMK8izn35VSO5VSOzhn+smKkCzgvwE78Pr50MifZrpBAEqpv1BKtQCfAF5VSr2WqbacdwZ/A3iNc467FzRNa8xUe/wopVYB7wOTlFItSqm/znSbzjMX+DLw6fNjatt5bTPTjADePD8HN3HORp7UUL9sRFL0BUEQDI5o5IIgCAZHBLkgCILBEUEuCIJgcESQC4IgGBwR5IIgCAZHBLkgCILBEUEuCIJgcP4/uVIo33QjqLIAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.plot(X, f, lw=5)\n", "plt.scatter(X, f + noise, s=20, alpha=0.4, c=\"k\")" ] }, { "cell_type": "code", "execution_count": 10, "id": "e454a83f", "metadata": {}, "outputs": [], "source": [ "k = GPy.kern.RBF(input_dim=1)\n", "m = GPy.models.GPRegression(X, (f + noise).reshape(-1, 1), k)" ] }, { "cell_type": "code", "execution_count": 11, "id": "9d491d23", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "{'dataplot': [],\n", " 'gpmean': [[]],\n", " 'gpconfidence': []}" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAagAAAEYCAYAAAAJeGK1AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABp7UlEQVR4nO2dd3xc1Zm/nzN9NDMqlrvlhrEMtgzG2JQAhhBCICEk2RTLwYSStUkhP+KSbMpmQ5JNNrtxgU21SIAEZy0IIYFNZ0PoxRWwjLENuHerzkjT7/n9MbqjO00aSaMp0nk+HydIunPvuVdX5z3nLd9XSClRKBQKhaLYMBV6AAqFQqFQpEMZKIVCoVAUJcpAKRQKhaIoUQZKoVAoFEWJMlAKhUKhKEqUgVIoFApFUTJoAyWEcAghNgshXhNC7BJCfCsXA1MoFArFyEYMtg5KCCEAl5TSJ4SwAs8Dd0kpX87FABUKhUIxMrEM9gQyZuF83V9au//1avVGjx4tp02bNthLKxQKhaIE2bZt2xkp5Zi+jhu0gQIQQpiBbcDZwI+llK/0dvy0adPYunVrLi6tUCgUihJDCHEwm+NykiQhpYxKKecBNcBFQoi6NANaLoTYKoTYevr06VxcVqFQKBTDmJxm8Ukp24CngevS/KxBSrlASrlgzJg+d3YKhUKhGOHkIotvjBCisvu/ncA1wJuDPa9CoVAoRja5iEFNAH7ZHYcyAY9IKf+Qg/MqFApFzgmHwxw5coRAIFDooQx7HA4HNTU1WK3WAX0+F1l8rwMXDPY8CoVCkQ+OHDmCx+Nh2rRpxKpkFEOBlJLm5maOHDnC9OnTB3QOpSShUChGFIFAgOrqamWchhghBNXV1YPaqSoDpVAoRhzKOOWHwT5nZaAUCoVCUZQoA6VQKBR5xmw2M2/ePObMmcP555/PunXr0DSt188cOHCA//mf/8nTCIsDZaAUCoWiF5L1SgerXwrgdDp59dVX2bVrF08++SR/+tOf+Na3etfZHmYGKivfnzJQCoVCkYG7776bFStWxI2SlJIVK1Zw99135+waY8eOpaGhgR/96EdIKTlw4ABXXHEF8+fPZ/78+bz44osAfOUrX+G5555j3rx5rF+/PuNxxYg/EOb4GS/vHGvjzYPNWOxl7mw+lxMtPoVCoRhuSClpa2vj3nvvBWD9+vWsWLGCe++9l7vuugspZc6SLc466yw0TePUqVOMHTuWJ598EofDwb59+1iyZAlbt27l+9//PmvWrOEPf4iVmXZ1daU9rhjQpCQa1YhqknAkysFTPux2K1abDasNRB+C4jrKQCkUCkUahBCsX78egHvvvTduqO666y7Wr1+f80xAfZcWDoe58847efXVVzGbzezduzft8dkely80TRKJamhSImXs+ZmEQAiBq8w+oHMqF59CoVBkwGikdIbCOL3zzjuYzWbGjh3L+vXrGTduHK+99hpbt24lFAql/Uy2xw0lUsZ2SIFQhGA4igRMQmA2CUw5eEQFMVCaJnMSaFQoFIqhRI85GTHGpHLB6dOn+cxnPsOdd96JEIL29nYmTJiAyWTioYceIhqNAuDxePB6vfHPZTouH2iaJBiOEAhFiWo9RinX1WUFMVChSJQ3D7XQ7lNaWAqFojjRjZMec9I0jbvuuot777130EbK7/fH08yvueYarr32Wr75zW8C8LnPfY5f/vKXXHLJJezduxeXywXAeeedh8Vi4fzzz2f9+vUZjxtKopoW3y1B7nZKmRh0y/eBUHfePPnnp16k0x/EYYbJ4yowm5W3UaFQDD27d+/m3HPPzerYu+++m7a2trhbTzdalZWVOc3kK3aiUY1wVENKMPVzp7Rnz5tUjk/U4ptRM2ZPKNB5Tl+fLWiShMtpJxLV2HO4hcljPHhcAwukKRQKxVBw9913J2Tr6TGpkSKVpGkaoYiGBMxCZFm9lDsKvm2xmE143GUcbe7k2Blv3x9QKBSKPJJsjEaCcdKkJBiKEAxrCCFixqkAFNxA6bjKHHSFNN460kI02rvkh0KhUChyj5SScDhKMBSFIUp86A9FY6AA7DYrZquNfUdaCEfyl5GiUCgUI52ophEMRYlKCm6YdIrKQEHM5ed0Otl3uJVAMFzo4SgUCsWwRspYyngorGEa4qy8/lJ0BgpiWSJut5O3j7UrI6VQKBRDRDSqEQj1pIwXG0VpoCAWiPS4nbxzvAN/IP8V0gqFQjFUCCG4+eab419HIhHGjBnDDTfckJfrSykJhSOEIlrRuPPSUbQGChKNVCCkdlIKhWJ44HK5aGpqwu/3A/Dkk08yadKkvFxb6y621WRx7pqMFLWB0vG4nbxztF0lTigUimHD9ddfzx//+EcANm3axJIlS+I/6+zs5Pbbb2fhwoVccMEFPP744wAZW2w8/fTTXHXVVXzsYx/jnHPO4aabbkqrdBGORAmGNUwmU1HFmjJREmrmQghcLidvH22jdvIoTKXwZBUKRdGzcPmvh+S8Wxpu6vOY+vp6vv3tb3PDDTfw+uuvc/vtt/Pcc88B8N3vfperr76a+++/n7a2Ni666CKuueaajK04AHbs2MGuXbuYOHEil112GS+88AKXX345oLv0YmKuxb5rMjLoHZQQYrIQ4h9CiN1CiF1CiLtyMbBkTCaB3WHnnWOtSmhWoVCUPOeddx4HDhxg06ZNvP/970/42d/+9je+//3vM2/ePK666ioCgQCHDh0iHA6zbNky5s6dy8c//nHeeOON+GcuuugiampqMJlMzJs3jwMHDgC6fl6srslUYkXGudhBRYBVUsrtQggPsE0I8aSU8o2+PthfLGYTmmbhyKkOJo+ryPXpFQrFCCObnc5QcuONN7J69Wqefvppmpub49+XUvLb3/6WWbNmJRx/9913x1tsaJqGw+GI/8xu75GKM5vNRCIRwpEokagsqV2TkUHvoKSUx6WU27v/2wvsBoYs2mezWugKSVo6/EN1CYVCocgLt99+O//2b//G3LlzE77/vve9jx/+8Idxb9GOHTuA7Fts6H2aolrpGifIcZKEEGIacAHwSpqfLRdCbBVCbG1taU75bH8oc9o42dqlaqQUCkVJU1NTw113pUZFvvGNbxAOhznvvPOoq6vjG9/4BpC5FYcRTYu1Wi9Fl14yOWu3IYRwA88A35VSPtbbsXq7jcEgpaSzy8+sydUqaUKhUGRNf9ptlBrF6NIbTLuNnOyghBBW4LfAr/syTrlCCIHdbufIqY58XE6hUCiKFtmtPl7qLr1kcpHFJ4BfALullOsGP6TssVrM+MOSVhWPUigUIxS98HY4uPSSycUO6jLgZuBqIcSr3f/e39eHckWZ08bxli4iqkWHQqEYYRgLb4eXaYox6DRzKeXz5L3PYiJlZXYOnmhnxqSqQg5DoVAo8kKpFt72l5KQOuoLs8mEJoVKPVcoFMOemAL58HTpJTMsDBSAw2HjhHL1KRSKYUq8b1NEwzxMXXrJlIQWX7aUldk5dKKds5SrT6FQZMnhU+0EQrlb2DpsJiaP7Vvp5sSJE3zxi19ky5Yt2O12pk2bxj333ENtbW3KsdGoRigSayhoTtpWPP/8c9z5+c9jtVr5/eNPsGrlChoffiTlHO+95mq+/5//xYUXLhjwveWbYWWgzCYTAQ06fAHK3Y6+P6BQKEY8gZCGw2Hv+8BszxcI9nmMlJKPfOQj3HLLLTQ2NgLw6quvcvLkyQQDJaUkFImiaZljTZs2bWLFypXccsutAGmNU6kybFx8Oi6nnaNnOpWgrEKhKFr+8Y9/YLVa+cxnPhP/3rx587j88sv50pe+RF1dHXV1c9n4602A4PnnnuG911xN/eJPMLduDrd86maklNx//y/47aO/4bvf/Xdu+dTNHDhwgAvmnQ+A3+9n6U2f5ML5F3DTJ5fEe08BPPnk31h0xWVcfNFCltQvxufzAVA7cwbf/tbdXHzRQuZfMI8333wTAJ/Px7J//jTzL5jHhfMv4HePPdbreXLFsDNQADa7lePNuX1QCoVCkSuampq48MILU77/2GOPsWPHDl7ZvI0//vkvfP1rX+HE8eNAbIe1Zu06Xnt9J/v3v8OLL77A7bd/mhtu+CDf/4//5Je/eijhXBs2/IyysjK2bd/BV77yVbZv3w7AmTNn+P5/fI8//+VvvLJ5CxdeeCH33rM+/rnq0aN5ZfMWlt9xB/esj5W2fu+7/055RQXbd7zKtu07uOrd7+7zPLlgWLn4dGxWC20+P6Mrotis5kIPR6FQKPokqmk8/cyzfOwTi7FaLUwYP54rrljE1q1bKS/3sGDhQmpqagA47/x5HDxwkMsuuzzj+Z5/7jk+f+edAMw97zzmzj0PgFdeeZndu3dz1ZWLAAiFQlxyySXxz334wx8BYP78+fz+978D4KmnnuKhjT29s6qqqvjjH//Q63lywbA0UACuMgeHT3Wo2iiFQlF0zJkzh0cffRSIJUCEoxp6VCJT6nhKO41opM/riDTnklLynvdck2Bw0l3HbDYTjUTin0k+V1/nyQXD0sUHsQaH4Sh4u/oOWCoUCkU+efe7300gEOAnP/0ZoaiGEIId27dSVVnJb37zCNFolNOnT/P888+xcOHCAV3j8iuuoHHTJgB2NTWxc+frAFx88SW89NKLvPXWWwB0dXWxd+/eXs91zTXX8NOf/iT+dWtr64DO01+GrYECcJXZOXZaxaIUCkVmHDYTgUAwZ/8ctszTqqZpBMMRgmGNxkce5am//505587igvPP4zvf+TaL65cwd+5cFlw4n/dd+16+973vM378+AHd1x13fAafz8eF8y9g7do1cUM3ZswY7vv5L/jUzUu5cP4FXHH5ZezZs6fXc331a1+nrbWVC+adz4IL5/PM008P6Dz9JWftNvpDLtptZEsgFKbSaWFMVWrfFIVCMfLId7sNTZNEolE0CVLGvDsjochWp+DtNooZh83KmXY/mqbSzhUKRX6Idu+UAqEIwXAUSUyWyDzCjNNgGfYGCsBmt3GiRbn6FArF0CClJBKJEgxF8AcjhCMaKKM0aEaGgbJaaO8MKZ0+hUIBkJNCfk3TCHXvkgKhKBEtljVnNsUMkzJK+nMe+JMYEQYKwOGwc+yMt9DDUCgUBcbhcNDc3NxvIyWlJBLV4rukYERLcN0N464XA0JKSWtrC2arbcDnGLZ1UMlYzCa8nUFCYVW8q1CMZGpqajhy5AinT5/O6nhNk2iaRBLbDQzzDhc5RGC22nBVjh3wGUaMgQJwOR0cPd3B9ImqeFehGKlYrVamT5/e6zHeriBn2roIhDSsVjNO+8B3AYqBM6IMlMkk6AzFsmvs1hF16wqFog/CkSinWjvx+sMIYcLpsOG2qe1SIRlxs7SrzMGx0161i1IoFAC0dwY40+YnGJE4HTZcZc5CD0nRzYgzUCaToCsi8QdCOB1q265QjEQ0TXK6rZNWXxCTyYTTYceWu5ZQihwx4gwUdPeMau7k7EnKQCkUI4lwJMrxZh+d/gg2u1XtloqcEWmghBBEomoXpVCMFIKhCMfOePGHNZwOO263tdBDUmRBTuqghBD3CyFOCSGacnG+fFDWvYtSKBTDly5/iLeOtvLO8XZMFhselxOLecSUf5Y8ufpNPQhcl6Nz5QV9F9XlDxV6KAqFIsf4/EH2HWnh0GkfNpsNt8uJSVXSlhw5cfFJKZ8VQkzLxbnySZnTzrEWFYsaThjVAdI1a1MMb7xdQU60dBLVBGVOu3oHSpy8xaCEEMuB5QATJ9Xk67K9IoQgrGJRJYGUklA4ij8QpisUIRyOEumu8Nck3ZX+KZ9CdFf+m7r10Sxmgc1qxuWwUuawYrUoVZHhQEdngJMtXUQRuJyOQg9HkSPyZqCklA1AA8T6QeXrun2hMvqKDykl/mCYDl8QfyhKOKoRjUqEyYTZbMJqNWO22hjobyykafjag0SbOxHQbbAsVHmcSgarxGj3BTjV2oWGiTJlmIYdIzKLz4jaRRUeKSVdgRCt3gCBkEY4qmEymbDZrFhstpy/pGaTCafdBPaeTK7OUJSW4+2YhMRltzCm0oXdNuL/PIqWdm+AU21daAhlmIYx6i+Q2C7qeHMnZ6ldVN7QNEmHL0CLN0AoooEw4XBYsTusFKJe0moxx919EU1j/4kOTEJS6bZTXV6GWWV+FQWtHX5Ot3WByYRTGaZhT04MlBBiE3AVMFoIcQT4ppTyF7k4dz4QQhCMxGT01ap5aOnoDNDc7icQ0rDYLDjsdqxFVsFvNplwlcUmP28gQnNHK2V2C+OrXUrDsQBIKWnp8HO63Y/JbMapimtHDLnK4luSi/MUkjKng2NnlEbfUBCJapxp7aStK4TJZMJhLx0RTpvVgs1qQdMk7xxrx2oxMXGUizKn2m0PNbocUYs3iNVqUaoPIxC1HOxGKZ3nnlC4W1YmEMFR4iKcJpPA7XIipeTQaR8Ws1CGaoiIRDVONHvx+iNYrVbcrlJ9b5K7yQ6uu+xIRM3EBpTSeW4IR6IcO+2lM6RR5rThGUayMkIIXGWOBENVM9qNwz587rFQBEMRTjT76ApFY7vskjVM0NCwAa/Xx6pVK4kZJcnatevweNwsX35HoYdXMqjIrwGTSeAPa4TC0UIPpSTRNMmRU+3sO9KGNFvwuByYTcPzFdMNlc1mY/+JDg4ebyMcUe/NQOjoDLDvSAtvH28Hc2zHZCnp+jSJ1+tj06ZNrF27Dt04bdq0Ca/XB2kq9hTpUTuoJMqcdo43+5g6vqLQQykpWjr8nGjpxOGw43GX7sq3vwgRc/1FNY23jrbhcVqYOLpcyer0gbHdhRCxGqbhk5MnundOsGnTJjZt2gTAkiVLDDsqRTYIozRMvqg7b57881Mv5v262eLt9DNzUqVSGciCcCTKwZMdaBKcjiJLxysA4UiUYDBEdbmdMZUuJbWTRDAU4URLJ52BCDabdZhnzUoWLFgY/2rr1i0o4xRjRs2YPaFA5zl9HTc8/S+DpMxp5/gZX6GHUfS0dPjZd6QNq9WmjFM3VosZt8tJhz/KnkMttPsChR5SwdE0yZm2LvYdbmH/iQ4wWfC4ncPeOMXcez3o7j5F9igDlQazyURnMKJiChnQNMmB422c6QjicSuV6HTYbVZcLicn2wLsO9yMPzDyVPO7/CEOHm+LGWp/BIfTgavMMQLel56Y05IlS9i6dQtLlixJiEkpsmM4L2EGhcMRi0VNGadiUUYCoTD7j7fjsNtxOpQLtC+cDhtSSg6c8uGwCCaNKR/Wen/hSJTmti46/CE0TDjtNtzukZaKL/B43AkxJz0m5fG4UW6+7FExqF7w+fzMnFylGpx10+YNcLy5E5fLoWIrA0DTJJ3+AG67hYljPMPmvdI0SavXT6s3QCgqcdptRZaFV6h6JFUHlYlsY1BqB9ULDqedY2e8ahcFnGz20doZwj2CMvRyjckk8LicRKIa+4604nJYmDi6NA2VrqXY7A0QCmtYbVbsDkdBdBR7o7D1SMnGSBmn/lJ6fxl5xGI20emPEIlqhR5KwZBScuhkOx3+SFyfTjE4LGZTLDUdM/uOtHLwRDvBcKTQw+oTTZO0dvh552grew61cNoXwma34y7ahAdVj1TqFONbVVSM5F2UlJL9x9qQJjMOh1JKyDW6odJ1/mwWE2MqnZS7imchEI5EafMGaO8MEo5KLNbiFPhNj6pHKnVUDCoLRmIsKjZptmKyqK6z+aTLHwSp4XJYGFOVf/V0KSWd/hAt3gCBUISIJrBZLUW6Q8oWVY9UbKg6qByi76JGCtGoxltHWjAXnXFKXkwNPxdNmdNOWZmTsDTzzvEO9h5u5sipdjr9QYZiMRmNang7Axw51c5bR1p581ALR1u6kMKC0+nE43KUvHFS9UilSym/eXnDYjbh6wwSjkSLbMLOPTHj1IrNYS+qHeNIE9+0mE24u2N+EU3j8OkupObDajFhs5hwl9lw2izYrOasmilGohrhcBR/KEJXIEQoIolENCKaxGIxY7NZCtYscuhIrEdatWpl/GtAuflKAGWgssTZrS4xZRhr9BWrcTIGu4GEiWbJkiUM9/TdWAPFHtMhpaTFFyYSCaBJDSHBJECYRPwpyO7/0aREk4AQmITAYjFjtViw2gTWYV+epOqRSh0Vg+oHXp+fsydVDstCS02T7DvcUoTGSUcmrH5BBbsV2ZJNPZKqWconKgY1BJSV2Tl2uqPQw8g5miZ560gxGycwrn51lHFSZEfv9UgNDRuS4lKxxVBDw4Z8DE7RC8U6GxUlZpMJf1gjEAwXeig5Q8pYtp7VZiti4wQq2K0YGgZTKzX8k3YKjYpB9RNXmYMjZ3ycPWl4dN3df6wNk8VSZNI0yahgtyIduXDLDaxWaqQl7RSKYl4yFyVCCKIa+PzBQg9l0Bw62Y40xYLmxU36YPeSJUtUsHuEklu3XH/dx0qhIl8U+8xUlJQ57Rw77aN2Sukm5R474yUYkTjspaEQEVuVGlfIQu2cRiy5zupM7z7O/H4phYp8kZMdlBDiOiHEHiHEW0KIr+TinEWPyUxLh7/QoxgQp1p9+AJRHPZSyzPur/imihEMT3p20Js2bWLBgoUJ7t+BGCf981u2bE7o3SRlJh3OvnZd6t3LBYM2UEIIM/Bj4HpgNrBECDF7sOctdpwOGydbu9C00nrxWjv8tHhDOB2lZpz6h8rMGu7kKquzx33sdrtZt249K1euSPg6/TuTOWlHvXu5Ixc7qIuAt6SU70gpQ0Aj8KEcnLfosdmsnGgpndbwHb4Ax1v9I0CVXMUIhj+5y+pcvvwOVq5cgc8Xe2d0I6V/nfrOZO6Yu2bNWjo6vOrdyxG5iEFNAg4bvj4CXJx8kBBiObAcYOKkmhxctvDYrBbafH5GV0SLvnjXHwhx5IwPj7us0EPJAypGMLzJfVanEKY+3pmEo3tVqFi2bBlCCPXu5YBBK0kIIT4OvE9K+c/dX98MXCSl/EKmz5SqkkQ6NE0SDoeKOu08FI7y1tHWEWKcjCgV6+HK0KV5p3tnYM2atZSXewzn1ufNxDT3hoYGvF4fK1euYOHCi+I/WbZsGXfcodLPpZS0+YIsOHdy3jrqHgEmG76uAY7l4Lwlgckk0KSgtcNPVXnxdZuNRDXePtqG21V8Yxta+puZpSglhiKrU0ot5Z1Zs2YtO3fuZNeuXdTX16MbpkzGUHctb9++LeH7zz77DMuXL0OIkVPZ0+4Lsv94G4dPejly2svhUx0cOeXF589e6CAXBmoLMFMIMR04CtQDn8zBeUsGp8PG8ZZOyl32rJSl84Uu/lpW5kCIkTQpq8LekUHuWqpv2LCBZ599hj179nYbImhsbKSxsTHl2DVr1tLY2GhIadevLVi5cgXbt29jz5698eNnzaplz569rFu3fti+e+2+IG8fa+PtI628fayNd462caY9fZZzWT+anw7aQEkpI0KIO4G/AmbgfinlrsGet9RwOh0cPtXBtAmVhR4KEHM9vn20FbvDjsk0/P4gekepWCv6g8Tn87Fnz15mzaqNL2iM1NfXJxis+vr6+Dtl3E0JYWLRoisTDNTGjRtZt279sHn3oprk4Il23th/hjcONPPWkda0xshhMzN9QiVTxpczeayHmjEeasaWU+Wxc3ZDdtfKSaGulPJPwJ9yca5SxWI20RmStHkDVHoKmyVXOvp6Q8fAXUBK1XrkkZhUY4wd9UVqgTD4fImZvXpWYKm69/QwwRsHzvDG/jPsPthCVyDRTeewmTlrYiVnTapkxsRKZkyqZMJoD+ZBLo4L0m5jeu0c+cwLm4el28nr8zOzprJgjQ1144TJgi3P7cKLk+wNjtJXG+kkJkjU19fj8Xh45pmn2bt3X8ZPGTP9MrmWSymLLxiO8tbhVnZ1G6Q9h1oIhqMJx4yrKmP29NHMnlbNrCmj+m2Msm23UZAZ7FSrn683PMtt75/LzMmjCjGEIaOszM7BE+2cXZP/+5JScuBYmzJO3fTP4IzspoiK1KQaKSUdHR1x41RXV0ddXV1KXMpoeErRtewPhtlzqIU39jeza/8Z9h1pJRJNVNCYNMbN7GmjmT19NHOmVTO6Mj8ZwQXZQZWPmyFnfejbAFxxXg1Lr5vDmDzdcD4IhsI4LFAzNn/dd6WU7D/WhjSZR4Bxyq4BXf9Xs6op4sgk9V3REyHq6+vjWXw6etKDTuI7kj79vJjeH58/xO4Dzbxx4Ay79jfzzrG2BEUcIWDq+ApmT6uOGaVp1TkPW2S7gyqIgZo993x581d+zv++8BbhiIbTbuG2D8zlPRdOHTZuv05/kHEVjryknmuaZP/xVoTZUgLK5IOjv7ui/hscVTs1Esn0Xrndbu64Y3nCOwE97j8982/JkiWsXLkingxRLC5hKSWn27rYfaCZNw+2sPtgM4dPdWCc9k0mwYyJlcyeXs2c6aM5Z2o1bufQSqEVtYvPJARL3zeHay+axgN/3MkrbxznJ4/t4OVdx/jcRy5gVBHWE/UXl9PO8dYurFYTbufQqZ5HoxpvH+1OiMh53KvYEgZS3XCZU35j/79q1coEA9WXcVK1UyOTdEk1brcbr9fLmjVrE46trq5m5cqVrF+/Pp75p+v2FdolHI5oHDrZwZsHm9l9sJk3DzbT0hFIOMZiNlE7uSrusps1ZRROe3EubAs6qrFVLr5808U899oR7nviNbbvOckX7/07n/un+VwyZ2Ihh5YTPC4nh056mT5eDIk4ayAUZv/xdhwOR86z9YozYSC9hFF9fT0rV64Akgso+2NwVO2UInEx5vV6E9LKIVYb1dzczLp16xLeRd3ll0+XcCSqcfiUl7ePtvL20TbePtrGgePtKfEjt9PKOVOrOWdqNedOq2bGxMqil2bTKbjZFEKwaN5k5kwfzU9/t4Pte0/yX79+hY+/exaL33NuydfweNxl7D/hZeo4D64cbps7fAEOn+nE43IOgVu0PwkD+d5lxVa2yaxdu46mpiaamppYsmQJUmoJK9q+DY6qnVIYEaxevSr+TiUbqvJyT1y/L/sd+sBp9wU5fKqDQyc7OHzSyzvH2zh4vJ1QJLUdyMTRbmonV3HutGrOnVrNxNGekp1HCxKDyqTFJ6Xkf194i1/9uQlNwsJzx3PXxxf0q/K4WPH6/Eysdg062Cil5OjpDnyB6BCrkvcdv+lvPCgXhkxKjaVLlyYEqY3U19ezevUqQLBhwwZ8vv7uAovNrakoLOm1+YzvU/q/Ef2Y/uHtCnH4ZAeHT3m7jVEHh0510NEZSnv8hGoXZ02sZEZNFWdPqmT6xEpcJTBfFnUMKhNCCG68fCZTxlWwrnEzW3af4Cs/fYavfeoSxlenrppLCY/byck2P96uIDVjywe06/H5gxw95cNis+ahZUZf8Zvsd1m5cxdK1q3r8funM1LGa/h8vu7dVn+KdXMnn6ModXpzEfcYJz39HGLuPn3Rnygum0i7L8iRUzF9usOnvLH/Pt1BmzeY9nin3cLksR4mjytnythypo4v56xJlUOezFBoispA6cybOZb/+txV/MdDL3P4lJevNTzLN2+7nKnjyws9tEFR5rQTjkTYc6iF8aOy3035AyFOtHQSjEhceRN97St+k21Li1zWF/W44ZLVonWWLl3KQw9tZP36HveelJqhil8ZnNIlV7vbwZUpQOw99njcCbGp+vp66uvrE1zN/mCYAyc62H8sJpqqC6Z2dKXfEdmtZmrGepgyLiYPNGVcOZPHlTO6Yihc+cVPUbn4kvEHw/znxld4/e3TuJ1WvnHbZcysKd62Fv3BHwiCplFeZqOqwonNYk54AYPhCG3eAB2dISIylhWYvxe0PzVE2aRl57a+SFedNhZMLl68mFdf3ZFSn9J36q9y6ZUCudqF9+c82R3bo27e+PBvsFVMwFY+gZlzL8FeMYFjZ3ykm2Kddgs1Yz0xjbqxMWM0eayH0RVlJRsv6g8l6eJLxmm38rVPXcqaTZvZ+uYJ7v7F83ztU5cyZ/roQg9t0DgdsdRzfyRK2/EOpJSYhAAp0SQIk8BqteBwFkLXL9uEgWyz5Pqb7t0bMTdfY2NjgmulsbGRxYs/kWCgdOOUabc2+ElPGbf8kKtdeP/O05eeY48Cwxla3AupuXoKwhSbUs+EgNM+LGbB5HHlnGUUTR1bTnX5SOswMDCK2kAB2KxmvnzTxfzwN9t47vUjfOeBF/japy7lvLPHFnpoOcFiMQ9B/dLg6VtstT9p2bmsL0o2nj00NTUlfK27ANPv1gY36RVnGn6pk8ng56pD8kDO0/M9KSX7j7ez7c0TbN97kn1HWhMUGBBmQh0nCbYfZYzHxHe+/kWmjPMkFc+rRUx/KHoDBbHCsv/3iQU47Bae3HKA7298hW/98+XDxt1XvPSWMNC/XVYu64tSjWdsktD9/snxqViNVOq9DHTSk1JLMG597dRyQ6F3a0N7/b4N/kB24enH3J/zRKMabxxo5oWdR9iy+wSt3p6iV5NJ4DYHOPb2DgIth7jh6ouwVGs0Nv6ZVuDbXzvA3Llz41mlIDN05lUGKxMlYaAAzCbBHR+aRzAc5dlXD/PvD77Id5cvomasp9BDG7Fk19JiIPVF2UyGidcoL/ckxJyMZG4Ul2myyoyeuq4XBg9uRZ8d+d2tpT57vY350F0/m90saXfhmdpY9CZdlNwOI3k3L6XkzUMtPPvqYV5uOkZ7Z09m3ahyBxfOGs+Fs8Yxd8YYnHYLt976V5rOvMWjj7wFxEodpJTs2rWLxsZGhIi9Z7feehtNTU1ZdeZVxCgZAwWxFcudH51Ppz/Etj0n+dYDL/C9OxYNK6HZ0qM3wxGb7HoMWc8xmSbygU7Gy5ffkVKca9zZQLqVsuTWW29LOI8ua5MuRdjYdRVgxYoVCcYt/U5tsORPZT3ds1+zZm08K23ort/XbranjcWsWbVs3Lgx/vX27dtYtOjKpHbqko6OHhUI45j18oR0u/nbl3+Op3cc5qmtBzl6pseITah28a65k7i0bhLTJ1SkxI4efPCBhESh2I6JhHEb35Pknym1/MyUlIGCmLtv9ZKL+PYDL7L7YDPfvv8FvveZRXjKhk7vTtF/Uie75NVi+p3TYCZjIUz92K1JfvCDNfGJV1/1JioGJLoRjV1X0006Q9PSO1fxl75I/+xjySeLqaurG+Lr9+56c7vdceNy6623MXfuXGprZ7Jnz14sFivPPvsMixZdyR139Cwq0o3Z7XYzf/6F8XN/ccUKWkJl7OuqYtl//jUeU6ryOLjygslccV4N09IYJeNz6y2+unJl4iImuTOvUsvvnaJOM++NTn+If73vOQ6e6KBu+mj+7fbLRmz32OJjIK0uUj+r0/8/4r5dhA0NG3jxxZeYM2cOQoiElPW6ujoefPCBNNdLHRvA5s2bE+quhmbCyYfKeur96ZmSq1atTIjrLVu2LMEY9Pc6qb8f+vy9S6lx6623JbS+WLz4E/z1r3+jra0tvrvSd826MdDpeWaSzkCEv289wJ9eeodTrV1AzEOzYNY43rNgGvNrx2FOO58kLlqMYsVutztB3XzFihXcfHOi6knmMY0ssk0zL9kZ3eW08a+3vItKj52m/We4/w+vF3pIijix1fCSJUvYtGkTCxYs7Mfk3bNj0On/hN+XGkRst9DU1BSPDxh54IH7M1xPxGNPRtavj7X0XrJkSYa4WvIisL+LwvSr9PTnGcy1Up+F3qBv6dKlCd9/9tlnkDJVB653JA0NGwxjj/1bs2Ytt956W/wd2bJlc/zdWbt2Xfw6Qgjmzp2bcMaHH36EtrY2Kisr2bNnLwsXXhQ3TsmsXbuO481efv6/r7Ps+3/hwT81caq1iwnVLj513Rzu+5f3YT31IjueeyKtcUoce4ympibq6upYuXJFwg7b5XLFjdOsWbVs2bI5xTjpY+r/+zByKFkDBVBd4eQrSy/BajHxl1f289dX9hd6SMOYbCa+xDhT8mSenaHpz2Q8UHoMaGNjY4oiRSzJIvV6ug6gEd3dt25dzEilK/ZMHH/s/hoaNmQ51sTd6NatWxImb+M4c3WtxHuWKbJS+teZnlM6Gho2sGbNWjo6vGzatIk1a9bGDZPRraq3rdANvv51X/fwsY99LOV7eq3c1q1b+MDHbuXJ3RE+v/ZJ/vTSOwRCESosfk7veITpcicfXjST+xt+wqZNm/B6fWnuq8cFqj9jXaC4rq4uvtBZsmQJe/bs5ec//3ncOG3cuDHhTHV1dSlGWBmp9JS0gQKonTyKz37kAgB+/r+v0fTO6QKPqBgZ3Ao+m4kv+Zh0k3nyOVLHlHkyjiUvyKTjB0OqAe190kjUAdSPNfYDigXpEz+TblLLPAmmH2e6uFrqbm2w10r/7B9++GEslkTx0Y0bN/ayW0w8p3FsRkOkx2H0rLYHH3yAVatW4vP5Egy+/nVHhzfuTlu8eHHKlR599NGEr5988m8A1MxayL/9/Hl2dkyibGwtArj6wims+8LVPPDtJXz4mgU0Nmazy8/sFehJI0/dgW7cuJH77ruvW83EQ319PQ88cD/r1q3H7XZn+RxHLiUbg0rmV39u4vfP7cNTZuMHn383Y6tUZh/kRi2h73hSYkbSypUr4orjyXEBfVVsVBnXpYvKy2MlAx0dXlatWhmf8PVMu564UC7SnNPHunqTRtKfpZ7arGf1XXHFIj7zmc9kGFdfMbVs64v6px+X/lq9k00Wn/E5pUvvznSu3lToE+Mwme+hoaGB9vaOFEkrh8NBIBCgsrKStrY2AJxjZjJ29nuQ9moAyuwWrr/0LD5w6QyDBmZs7utfbK+3WGDq2I1affpz05NPsnmOw5W8xKCEEB8XQuwSQmhCiAWDOddguel9c5hfOw5vV4h1jZtTmnaNTPq7qk73dTbxpMRjFi68KME46X1zjMZJ3xVt2LCB6667jsbGRtrbO1i2bBlSSm677fb4Dq2uro6mpqb4PaxZs7afu5Cku+o2iPp9bN7cs3PK5KqDWCq70XDqMYfOzs5enm3mmFr/XHJ9xdV6v1Y26PeXfLxeAG3c0Sa69xJ/B8ZCZv3+9N1nOhKfQeZ7WL58OUKIlLiObpw++k8fpWzcuYy/9J8Zc8EnkPZqystsfPLa2Wz4l+u46do5ceOkuxyTu+XGFkOZ3qne3M/pd6CNjY1s376N+vr6+N+G7nrseZcUmRjUDkoIcS6gARuA1VLKrdl8bih2UADeriArf/gPmtv9fOiKs7nl+rl9f2jYk92quu+dVnaisMZjtmzZnPQHqK/6Zbxo0Uht7UysVls8S6u+vh6Px4PX6wXg4d88hqWsCktZFTXTZnH9DR+mvTNEmy9IVyBMJKrF/kU0TCaBw27BabfgsFmocjsYO6qM17e9iAz5GFtuxd/li++YdMPZv11ZNs82/ap61aqVCdlmq1ev6kemY4bRGOrAMo+nf/T1XvRVEGscS23tTPbu3ZfwHKBnd5G8G093D8m7WJD8YM06OmQ1b7c7CWix9hORgJeZo4L8+5duwWm3YXz3gLi7EGJCw8ZMzsWLF/OlL61OemZ9exKSC5qNu8bkBAlj37KRSF7EYqWUu4GiET30lNlZVb+Qf73vOR5/7i1mTxvNwnMnFHpYBSa1viS5oDRZuif5j0+f+Iykaumlri5T64J6rqnviowYJ68bP34zF13zER554u+8c9SEZ/RkJr/nS/GfB4HfP/dWv55EjIrYv3bJ9Akz+dy3f8Ub23Zx7eVz+dqX+uNu6Us2JzF1Wp/EjEkBq1atjNdeDbYuxlhEbHSzDkZSCpLVQnruWZ+Ae3tvkmuA9u7dlyLwq7eo0Lsk659PNtr6+Y3jCYWj/H3bQfZpczjT5gfAbgozo9LPOzue5Jln3+Sfd/8j3n7FuAgxdst9+OGHgUQFiIaGhqSFSt+KKMnKKkKY4kXFydl7iuzIW6GuEGI5sBxg4qSaIbvOOVOrWXrtbH71l1388NFtrP3C1SNcaSLVcCxdujTueutZEXvi7htjYWOyzlx6Lb3Uibi/entWz3gco6Zgr5qCvXIyr7aX8eoj24BKHNWVhCVo0TARfyuRrlZmTB7He656FxUuOxVuOy6HNSa8azZhMQk0KfEHIzz2+yfwdga5/Kpr+L9nX6HDr9HcEUTYq9h/vB2ooLrug2xrg9u+9SjXXzGXhedOSKsYoGNcxacz3MZdmXFSk1LGDbOe4p7swRiYIUksItavb0zgGNxKXaTdKRkD/ZneGyOzZtXywAP3x5+rELFJX0ri8kPGHksNDQ3x52ZMJPAHI/x18wGeeH5fvMGfyxphor2Ff//Sp/nv/76XvXvejKeeG9O9dQMOsXICYwan/t42NTV1p7MnGubly5fHn4fxM4nPNvE5p3uH9N1UT4lDcSzwi5E+DZQQ4v+A8Wl+9HUp5ePZXkhK2QA0QMzFl/UIB8CNl89k1/4zbNtzkrWbNvPvyxeN0CJembKi1V0OupG67bbb45lUmXYDyerhiStH0hzTu4LD2rXrePg3v6X2wms447fjHHM2Zntix2SbiDDvnBqmT6zktVf+wd//9CjRQAeg/4H/mLPcLX0W/dojLTz22CaqbT4mWeC5P8ZWssJsxeYZj61yEtU15xCxj6EjZOPhv7/Jw39/k8ljPbxnwVSunDeFCrc94Zz6rmH79m3s2bM3HpyvrKxM+H5PDC42vttuux0pZXyCMrpDdQam8p6oNmFUcc9NEL53hY/knVJvixrjrjqdS2/VqpUJxa8bNjQAsHz5cjr9If740ts8+vddRGSsA8D0iRV87KpZXDx7Ivfd18C9996T8CyAeOzLKHGkX9eI/tx0I290We7du4fa2lmG+Fim3lCpBcjJngLQ3dcqe68vcpLFJ4R4miKIQRkxxqM+cfU51F9z7pBer1jJNpsq2Q8PyX7yRNFQKSVCCMMfaerqMnkFGo1qvPb2aX7+yFOc6rKhGXJ0Iv52wu2H6Tqzn0DrQaL+9rjR1I2ocVzQe0tt47XTqT8k8/HFS7jmxk+yZfcJXtl1PC4QajELFpwzgQ9efjbnTq2Onzv5nMm1QtXV1TQ3N8cnxOSYh+5W0j+bnOk40J3U0KlN9J71mPz9dJmaidmRiYoMmc69du06HvntE1zy/ltp0cbQFYwA4LEEmDkqgCPSwurVMaOhP+O6urqU3ZFOLC4q+jw22bBm0vDLFH9KzoBMF3ccyUkS2caghq2BAtj59mnuvv95hBB8745F1E4eNeTXLBy9pSEn/kxKLeEPMnnC7Ct4nS6OAJkFVs+0B6iZezXP7DhMm69HGTrYdhTReYTmg68T9p0CYoH0Cy6Yz65du7j00kvjgXbdWOoTS4/xhORJOJ1RTjdZGTEa40hUY9ueE/x960G27zmB3vJn1pRRfHjRTBaeMwGTKTE9OfkZZjL4UsqE4/RjY0H5gatb954gAbkxVKlJMJl2Sj27NxGfvPWvgZS2E8m/o61bt9Dc7uff7nmY4343dDcCDDTv58Lpdr715c+wfv09Cdc2LgIyoS8G9EVPut9T8mKjN0Pcl4tb/1vxeDxxFfyBJ+UMH/KSJCGE+AjwQ2AM8EchxKtSyvcN5py5ZO6MMXzwsrN54vm3uPeRray582qc9pLTx+2TvjPwEg1Xcmwg2U+uT/x6jED/vDGobxRZhR5jZjSG//nfv2Dr20Eizhq2dSc0aIE2Og6/yiRPiHfNnsHKld/joot6JqYrr7yKO+7o2Y01NGxIWKXq19RrptK5WZJ7NV133XV9PkPjDs1iNnHx7IlcPHsize1+/rp5P395+R32HGrhPze+wqQxbtz+txM+n2x09PiCceLTJy4jumETQrB69ape6qMyY0yQ0CfE7du3sWnTpniMq8cYpHdDZVtjZcQYg0rn2r3vvvvo6IhlYDY2NsbH8tRTf+fUqdNxQV5NkyxdelP8vFbPOO64+yHOhNxIKsEE/tP7aH/nBULtR/nOhi0J1zLGv5IXCtBjcIxySPqxunHKZISAhHvLJGabbizG2KPxvPnpHzY8GNT+Ukr5OylljZTSLqUcV0zGSeema2czdXw5x5s7efBPOws9nCGgP7VOqbUamfXBjFlTOiKhlQAQz0AzTlJ7DjXzjfue5ZUTo4i6JiMEdB7byYlXHuDIsz+mY/8LnD97BitXruTmmxPVJnp69cQm6Z7aHBLUCJYvX95r3ZGxJqu5uQWITVRGKisr49831lkZqa5w8Mn3zmbDl9/H7R+Yy+hKJ0dP+9jjG8e516/mZxv/l9qk8+rPNJ2Sxs6die+gEII5c+YYDGTPAiN9PVRqzZGeIFFbOzOunbd37z5qa2eyc+dOGhsb8Xp9bNiwIUmRI5bun1j7E3NNGa9trBszqmhs2rSJZ599BpfLhTFTM/beLM+oHnHqVI/ay89+toFrr72WvXv3MeP8K5j7/hVMuPSfOR3yoElJtdXLeeWHOb3jEULtR+PPMTGjMPF5Gpk1q5YrrljErFm1fOyjH034WXl5ebyPmO5uW7lyRcp7ousB9i7Blbl+y1gHqL+Tg3PljhyGvQPUajHzxU8swGI28eSWA2zZfbzQQ8ox/RFmzdwuXdcsS5T6SSZ1Fa2zatVKjp728V+/foWv/uxZdu1vpsxhYZKjjaPP/ojmpicItR8Deiariy7qKehdtmxZXNfOOCFs2LCh+x4S9fP6us/k+6uvr0+JEf3tb39jyZIlLFp0ZVLQOlZ4uWHDhvgYPnvHMt586fecLZpwe3dCxE9n1M73fvUy4fHvxl45uecpSRlvBaE/V/2ed+3axeLFixO+F4mEEwxkpgVGusJe4y5m7959CYuNvXv3sWvXLurr61m5cgVeb6xHkm6QjDES3XDpafCxnY+MX2PPnjcTNPj0idxYqKyPyfj7Sm4vYfx9NDY28vP7HyDsmsqkyz9LeNwi2iNlaJEQHQc3c+z5nzDLc4o/PPqrtDqE6YyG7v41SlE999yzsf9//rmEY5999hmWLVuWkpmox5qM19NT9tONw1ioayRbA6bIzPDzd6Vh6vgKbrp2Nr/8cxM/eWwH935xFOWu4dQ/Kvs21slunvJyTzxga/wjSk1N7pk4jW49AGG28blv/ZJT4Uo0TWKzmvngZTP48KJayuxmFj7x415Hr9ep6G4Yl8sVd4Hoq1m9tUN295nZkOr85S9/jq9sddfh8uXLaWjYQEeHl507X+fEiZM0NzdTUVFBe3t7fJdTWzuTw2/9CffkC6k6+wqi9tGMu+hT+E+/RWX4HR5++GGqq2Pxzrq6uu40aQ/V1aOYMGFiPN6k//zSSy9l/vwL07qHjHHE/mTSJWPcWTQ1NcVjScbfZbr4jfEaK1asiLcVSTfOVFdzZgJRMxUzrsRdMy+ewRkJdOA9tAXfkR3ISDA+1kzvZ3L8Sy/+1hN49GNdZbEyk0yCt8Y4XaZs1L179/SSpZpNmUX6jsDKSPXOsNHi6wtNk3zzF8+za/8Zrji/hhWLU1N8S5fBabDpbcwzZ1zF0Cdv6Fmp7th7is7yOkw2DyC5ZsE02t5+ltEVTpYtW5aSMZisJqCPNV3vHD1925jRlU4XLl3BcF/B83QB7mTDaxxDOurr63G6ytl7xsLp6GgCoShSSjqPvU77W8/w8Y98IKG2Rp9sjQWoxmSPbJQ6ss2ky0T6fkSJ104+Jlk7MP04s3vuZ5+/iDYxDueYs+MZbKGOE3QcfIWuE2+A1Hp5RolJP+lir7fccitCiLhuo/4u79nzJtu374iPQ0/wSJ+kkCnhKHMiUu9x4OUZDdhIdfPlJUmilDCZBJ//pwtY8d9P8dxrR3jX3ElcPHtioYeVAzJLsEDyDiN9gNwoSWNMhU0O4Oqp5A0NDXxs8U1YJl+Ov+UIJkAEW7h0msS/fz+PPRybVNauXZdQD2SsqNebAuq6ehDbSRkTJvQmdEaBTSBlAku8z97dmA8++EDK55LrZoxkMk5AQgp+uy/Eb/6xmz++8BbuSedTNn42Y+vO5T9/cA+PPhIz5np8KFU1ItvVdXpVEGPqMpDwnNJlDRpJ1qJLR/LCJdM40yUKAJw79wKqz7qYk8FygpqVMkBqUTpP7GLtvy7ne//2RU4cT9Xpy1xLp7+PqUkfc+fOpbGxMT6u5J24ju6mTJ/mnUn3MLMeYupYEotws68TVBgZMTsonT+++Da/+MPrVLrt3PvFa/CU2QoyjlzSdxZf9is8neSVnfHzuw+2sHbTZlo6ApiEZIqzhWcf3wDd75Ix686oTL5hwwYee+y3jB8/gV/+8kGAeBzEmFKeCT1DTTcM+oTpdru7XYCJemvGsXd0eBMMSrrsv3SFs72RXCd266238ebbR6mYeRWu8bMBiIa6KPO/wy/v+Qo2qyXN7iOzeygbXb/kmiP9mQNxHcN08R9jyn6m7xmvYdzh9jbOnnRxgaN6GrMv/SBt0Yq4eLOIdNL6ziv4jr2GFuqKu9r0ZAafz5dS3tC/CTxzjVrhdy/ZKtcPf9QOKgPXX3IWLzUd5Y0Dzdz/x9e56+MFFWHPCX2t3vqKX+jfyxzb6fn8UX8FhwOjiWqSYNsRLp2m8fUvfYGFv/9Z/LN6Ff7y5ctZs2Yt9913H8uXL8fn89Hc3MJ733stQHxXtHjx4vhEmlyDYuTBBx/oTlXuce/ok7PRWCXXl/Q8H53kgHXfMSsjupvSqKmnu/EArjnXxu/+/ACVtVfjGDWVoK2OW771O6Y4WxLOc+utt/Hggw+krK71PlWJSRskKCzoLjRjDCr2zGP3qqc268/3S19azS233MquXbvi6d6rVq2Mjzk5O1MvNNYTV4zqGPr43G63QUoplu33yds+S+XMd1M2YS4Wh4czITAJjQXnjKfr6A7+8tsHWLKkHre7Lp4aP2tWLb/61UPcc889ccmt5Hs3/q56n9RTd5mLFl3J/PkXFsHuJfMOTJGeEWegTCbB5/5pPit/+BTP7DjMu+omDRNB2d5e/vTul57gc18uJsFnPv8FdnvHcaDLBUg6DrzM++aPZfWqVM01XUZJnyCNRlAfQ7Kba8OGhoSV7ooVK7jiiisIBnsKe3V3lF70aKzwN06gPXprxhhBcqU/KarbvcWg7HY7wWCQ2tqZ8QmvqamJ8nIPQpjiySZGl9qprRtxjD6bcXXXEbRVsK8zlpo+tayZv/z2gXj2nNE4GHX+jFqJe/fuIRgMJfzOgO4iUHeSq0rEa9j0bq8bNjQwd+5c6urmsGvXLjZsaMDn83HppZd2152lJs0Y41pGEVpjfdOePXupO38BK79zH2+fkZgnXkd59ygcpjAn9r7IexZM52uf+jANDa+zZEl9ilbgFVcs4p577kkwtkKYBtjLLHWxoRfIZl7AKYqVEWegACaOdnPTe2fzwJ928rPfv8q506pxO0vf1dc7qStLIKGIcNasWhYtujIlJnW6zc+/P/giLSEXWjhA864/4D+1h9W/TlQSMGr9GfXgkgsa0+3U9Cy9+fMvjMecdOOku/aM7qjkuFF6vbXsdpDGYtOGhgbq6+vZuXMnQgjuv/9+1q9fj8vl6t6ZeQzuROLXWL78jrRpz7GOrT9mXO0lmMcvoNNRzhveiUxf9M9YWnfi8XjQJ9/k5oDGceqGV1fw0Hef9fX1BpmpRJYvX05HR+LOVP//cDiUIJza85me3Xi635WUsWf58G8e4+obP8X5H/gim1scCJMZswvQQrz34plcNX8q50ypYt26w3g8Doy7u+QFk/67S07G6H3X33vDxr7jsco4lQIjLgalE9Uk/9rwLHsOtfDu+VP4wscuLOh4hp70mnR6z6XkVg16htM1Nyzmu798iVZvAKc5xNtP30fE3wZk1lxLlqzpK4aSLpEjOWNQR090SJdRlv6ame8/eUJMHIOO8Wd9uZdS+1zpuzKPx8PPf/EAnikLKZ/+LkzWWOO8CmsXX7z5Gp78/UYaG9Pv4vqW2ul9TJm0CHv/fOq9fGzxElpDLtq1Sk77rcjuMkqpaQSa36Hz2Ov87Xf347BZE86TeXy9awcORL5pYLsuRb7JqxZffykGAwVw9LSXVT98ilBE4+u3XMqFs9KJtg8H+pt6HfvD3r7nJGs2bSYQilJu8bPrbz+m/uP/lFAPo7vjTKYet1l/9coyT5S9teQeyMSbHzHV5CzDZL03k9XJRdffxslABdHuST7Yfoy6SRa+/aVPY7eaM97zwMbfH0MeO15/P8x2D1dcv4TWcBktATvCbIkfE2g5RNfJ3XSdfBMt1AkkLyCye1466WqqVq5ckbDgWbZsGXfcsbwPo6OSEYqdvLR8L3UmjfGw5L2xbKuf/m4Hnf5QgUc0VKQWIK5evSq+WtcxTur/2H6Y7z30MoFQlDE2LxdP9FL/8X/C7XZz881Lcblc8R3U7bffzi233BoP2utN6Orq6hJUKZJ121auXJEUEE8lXWwsWXoHSPj/xAp/nb4q/QdLzzPWM/t0JYVdu3bFY3Fbt25h8cc+zMtP/IQLKw/Stu8fREOd2Csmss83lk//x5/5/Ld/idU9Nu0993/8mRNAjCoIOtGoxpsHWzgRHc85169m0pX/j3e6xtAadoHJjDl4hk/fMBf74f/l1NaN+A5vQwt1Uls7E4gV1v7gB2tIfvbpxpRJmcHYCPGmm25K+OQzzzwdf89Spbx0VDLCcGFExqCM3HDZ2bzUdJS9h1t58E9NfP6j8ws9pCEhvVBoInpixJNbDvKz3+9ASqhxtvLi4z+hNk0x7UMPbWTdunVpXVp6fEQP4jc0NCQFq3tkejKtgnuLJ+jGzu12x+NW+vnmz78wyehlG5sY3Mo7XTbl6tWraGiIJSgYs8iklDQ1vU7H/ia8BzdTNqGOqee9B1/AQRcVTHjXMiZUuxCdR3nsz4+zbfs29qZJlU4cf+ZnqMeeqqtH0dzcQm3tTDZt2kRUk3RG7YTMlbjGzmD3gWb8wQjgAsBhs9ByuAn/6bfwn3mLV174B+vWrWff7p3xcwFx7T+r1cauXbtSxpC42+m9f5iuiacnvhjZu3cfe/fuS1CXUAxfRryBMpsEd350Pqt+9A/+vu0gl86dxPzacYUe1hDRM2FnSlk+HvDwTucYAG6+bg4fWTSTtWWtCdl/+mSnF9UajZI+OSUXyiZPqH0HvHufxDK3Ik+9177OZXQpDT52kXofmRYHTU1NSQZnLefOexeeyfPxmUZzvLkTqGTcRbcQ1EKcf3aEmnnz2b7nJLct+xxSGptGph+Lft8ul5vZ5y0koFl57pXXiI6ZzbkzbuTlFhvoGYCtJwGYUO1ifu14LjxnHE8+vpGHX/tt/IxG7b9k95ueYGMslM30e+6rNEIIE4sWXZmx5CBRWUIxXBnRMSgjv3t2Lw/9ZRfVFU7u+X9X4xrGWX0NDRt48cWXDEWvsRTubfu7CFbMBeC298/lg5ef3f2J1D5AqYkQfUn1DFSOKftdTd9GJrOETb6laPoqnI5GJW8ebGHz7mNs3n2cky1dKeewmE1UVzjwdzTjtAnOnjEDi9mE2SwIhzV8gTCd/hA+f5jTbV2EI1rKOaSUlJnDXDZ/JnVnjWbO9DFUVzhJt4jJpMunky49Hfonu2UYWcYY48DPmXh+FacqHKpQt5/ceNnZvLzrGPsOt/LAn3Zy50eHa1ZfLH23qakpIWX5Dy/sY9S5sb5Jyz54Ptdfelb8+OQYRnIbiXRyOalSPdkL2iaSbTxhIGnJyav39P18BjdxpZ8I+9pBmM2COWeNZs5Zo7n1/XM53dbFvsOtbPr9k7xzvAP3qElEonQbLicdETi582ivI6l02xlf7aLjzFF2bX+eUMcJQh0n2Pzy89x22+1shu4Ehxj6O5K863S73QlJMqtWrYyroEO2gr69P7NM+oj61327NzOjMv1KB2WgujGbTXzhoxey6kdP8dS2Q1wyZyILzhkOBbzJJE7GjY2NlI2fzejzPgzAHR+ax/sunt59bOLOwljnlE5fr3eNvPRJCrnbnQzWyAzUgGZCJhUHx0icCLMzvkIIxla5GFvl4rLzbo2lf//1IYTZitnu4T3X3UgUM3aHk4A/iM3h4Jqr343LacPttOJ2WqmucOK0W5FSY+nSpXgP9rjObrrppriAr144rNdkJSbS9BjRhoYNCW7curo6mpqa2LlzZ8rvuacYOX2cLJ0BN7pk9do00BsvLk9pppk9A1nIKAqFMlAGasZ6+OR7Y205fvq7V7nnruoS1urrzYXRMxk7qs+iuu5GQLD0fXMMxil2XHLsZtGiK4FYvEEIE6tXr0qSy0knJdOfAsrBMBgjkzsDmlzDpbvR9SLcwU6EujGQ0TCRrhYO7HoxRUXjygsmpzm/TFBw0OM7e/fuo6Kiguuvvz5BzNaoe9hD+hjS6tWr4ioaehzSuPvpaX/REzPM1iWbKFXVU/A9sOc3lLtlRa4Z0Wnm6bjhsrM5Z+ooWr0B7v/D64UezoBoaEhtbLd2rbFDa+xrW8UkRs/7KMJkZqKjjY8sOttwlthnjR1tIdaXaePGjXHlB4i1NkiufVm5ckWCuyRda/C+Usz7z0BTyXtPe+5fKrrM2Ek2ubdRruhdRcNIz4Jj48aNCT9pb29POVp3AWcm0XAZ+1xBzGh5PB6DXl+PEdqwYUMWnaCTDWOqoRwYqnlgqaAMVBKxrL4LsVnNPPPqYV5541ihh9RPem8Br9fTPPrEk9RccjMms42xdi8vP/Hjbk09mcbAkWDgUlsU9EweDQ0buPXW2xJqbNasWctLL72U1ARRbw2eK5//YIxM+iy/gRnQxM6/ycXQg8s+kwkSR+noa6JdvvyOeEq+kdramWkbFvZ3bEbWrFmL1+tlz569+Hw+jL8jveRA/x313iE5l2TT/VZRLCgXXxomjnZz8/vm8Is/vM7Pfvcq506tLqEOvH27MMyOcqZe/mlC0szCc8fz5U/eyD3u0/GU5YH76CUdHd6UpoJGl1GyqzGX9z2Ynju9Jy30N+Mrve4hDDbuFrvHTMYpu/NLkjvRplMV0Xd+ia65TGR24epF25nexdzG/Xon2fWqP8empqYhcDcrcoFKM8+AsQPvu+ZOYvWSi/r+UMFJ1JNLl/bd5g3w9YZnOd7cyZzpo/nXW9+F3Wo2fDb2PgwsJTx23UwTXn5qV3KbPjxQRe3kZ2CMxwxul9Bz7mQBWGMMKrUZX+bur1JqXHfddfGiW+N4Y0kJfe9y+0qbz1auauh2UD3X0sWHIdFQZXuv+vlUmvrAUVp8OeBESycr//vvBEJRVi5eyOXn1xR6SBlJnCBIO0F+9s7/x7/9/Hn2H2tn+oQKvr3sClwOa8ZzZNOCPP0faar2W2517/LFQOqjBjsR9j3xJffC0muS9K/37HmTWbPO6cOoptaA9d1qve/nlTr29Isdo/jt4GvPsjUWvRlEsr6eSlMfPHmpgxJC/AD4IBAC3gZuk1K2DeacxcT4US5uuX4uGx5/lZ89/iq1U6oYW+Uq9LDS0BN30hccxiD9zp07aXzkNzR1TKQj4mRCtYvRgdf59a/2JUxYyecwkpwqnOmP1O124/V6+/x8aTCQjK/0Lef1tOhYa4zsC43XrFmbZNBkgjtyw4YNrF/f075c0zRuvnlpFi7aHldm7tqRpx6bjVzVYK7bP2ORC5eiSlPPJ4ONQT0JfFVKGRFC/CfwVeBfBj+s4uHai6axY+9JNu8+zrrGLfz78kVYzMWUW5LaX0dHzxhbs2YdJ63n0hFxMqrcwcTobn77yMaUCUvXiMvkmuqJR2SOU+ndZpM/q5+zVI1Ufya1dKoV2cRxkp+p3u4i9hx7diPGNOzk3l3r1/ekkfdtVNOlcffc8+B/T/2RqxrIdftrLHJRSqDS1PNJzlx8QoiPAB+TUt7U17Gl4uLT8XYFWfHfT9HSEeBjV83ik9fOLuBoEmMJHR3e+KSf3Itpy5bNgODHj23nqW2HiIb9nNryEGHf6axbUmzZshkhREI8obf4gdvt5qWXUmWUmpqaeNe7Li1BF0hh4iRGkl1vyW1RMrnQMvflyrebaijjNdn+fnItZzWUrVuGP4Vot3E78OdMPxRCLBdCbBVCbG1tac7hZYceT5mdFZ9YiEnAb5/Zw863TxdkHInp3zLeLfXWW2+LKwQYWbp0Kb/6SxNPbTuE3Wrm9PaHCftiY+8tfpJ8Dilld0r4ckO6efpakjvuuIMHH3zAsFOKKXo/+OADJW2cBl8flQ2pz1TPpsuchp36mXRp5Ml1cX3XIA2WdE0fjV/nimxrmnJZSqDS1PNFnwZKCPF/QoimNP8+ZDjm60AE+HWm80gpG6SUC6SUC6pGVedm9Hlkzlmj+di7z0FKuOeRrbT7gnkeQfKk0kNTUxMLF14Ud+1s3ryZWbNqORYcxePPvYXZJJBH/kGovUerLVO/JH0S3LJlczwrTDdS6eqpjPScM11RZT6y93r7eiDkclLLhsy9m3Qy7QyMLF26tA+j2nMfQ1GD1HeheC7J3lj0FJ0PphYv34uWkU2fBkpKeY2Usi7Nv8cBhBC3ADcAN8lCpATmkY+/exbnTq2m1Rtg7abNRKKp6tBDR+qkkq5gc+PGjZhMJm7+wneoqn0PUkpOvvoYb732HLNm1bJly+YMf1CJk7EQJjZu3Bg3UgsXXpSgyWfMwCr0H+lQToi5mdTSkfiMjE0Y9Weq756MpLvP5N+DvlDRe2+lN6pDpaaQj91Zz7X6byxSd1b9I9+LlpHNYLP4riOWFHGllDK1H8Aww2w2sbJ+IV/+ydM07T/D/X94neUfmpfjq2SnoZeJdevWs+Dqj3Pf/8Zkmlrf/CtdJ94AYsZLCFPGbKnkoLVupIyxDP2PMneZX4MlH1lVg53UEkkX/1m3bj179+5Jyf6Dnjbq6fQLM/0e3G43QmQ2RrlJGEhHPpMICvMeDj65Q5Etg0qSEEK8BdgBPaj0spTyM319rtSSJJLZc6iFb9z3HJGoxmc+PI9rL5o+wDMlCmDG9Mm8hvhNqnhmpiD6qlUrWbp0KYdaBWPmfRQQTC1r5rnf/yx+XGIBZ2Kzv0zj6z0AXSzFivlMZBgsvQfrjQW2yUkwmRMZUn8PiUrqqYWzuU0YSH+f+UsiKJb3UJEteUmSkFKeLaWcLKWc1/2vT+M0HJg1ZRSf/cg8AO574jXe2H+m3+fQNetivZQkUmo8++wz8aSHZNeIlFq86+3ixYsTXHu6e+jgmQhjzv8nQFDjbOW53/8sHk/SXR+xeJKGPmmtWbM2yRXWu/soOZaRSKEmhVIS/+w9/tOzeIit1PVMSONnY5mUiedMpnc3G0Pspsp3EkGxvIeKXKO0+AbIu+dPZf/xdv7wwtv81/+8wvc/exXjR2VbxJtes05XpW5qaoqvPvVJpKGhgaamJubMmYMQgpUrV7Jjx3b27NnLww8/gnPMTMZduBiJiQmOdmaPCTPVEC9yuVzxeNK6deuzqrMpHjdeXwx1r6lck+qq1eNFybuk1J1Qci1U5vNDZjfb0Lmp8tVaRTESKKaK05LjluvqmDdzLB2dIb51/wu0dASy/GQs9TpdIevixYsTjjQWxjY1NRGJhGlsbOTmm5fGDVrZ+NmMPv+jSEyMt7fzyhM/orPTl5DM0NnZyUMPbUxYuRuNI5ASzB66BIFcUopZVZnS+bWE34GmaQk7oeSf935v2ewqh2LnoZIIFLlDafENkq5AmG/+4nnePtrG1PHlfOvTl/dD+TxVsy45c8sYNE8Xf3JNmseo2e9HCMEkRyv3/ustrF9/T69xI+M1M1+vdCaS0tJGy9ylWEcveNZbUuiLDOPPs1GmyLaAdWjiNyoupMhMIQp1RyRlDiv/esu7mDTGzcETHXzzF8/T5s20k0qsO/rBD9akHKGnjqfuBEhZEXumXUL1nA8ghGBqWQsvPfET1q+/p9td1EPyxN0bpWacYChTwYeC9On8RlauXBGXMFq3bn3K79PoDkykP/HDoa5XUnEhxeBRBioHVLjtfOvTV8SN1Dfue46TLZ0JxySrQPzgB2t4+OGHAeI9c3T0Xa3uGok1+qPHuAgTo2Z/gKra9wAwvewMa792G0uWLMHlcqVVERhYnU0pUToTYqJBlSm/L90o6UbFmOav/3zDht6MSzZutnzWKykUA0MlSeSIUeUOvrNsEd+6/3kOnujgX376NF+9+RJmTakmXa3Orl27AJg1q7Z7JR2bUHfu3ElFRTn6BKurg+vtM2bOnkeH+3zslTVo0TCOlq08s+Mpbt79D8rKXHR1dSb0BNJVBfRz9bfORjFUJOropUsoWLFiRYKLbsuWzXF3n57wAulrv/pOglCip4riR8WgckxnIMya/9nMa2+dwmI2cfN1c7jhXTMQIjWGpNcv9TSWS6yLMk5gc+rqmDDzYvZ6R4HJRiTQgat1M1dePJfnnns2pXHdxo0bEyazRYuu5I47Eiet7OtsFENFb21LjErlkNhHKdPPleipohRQDQsLSCSq8cAfd/Lnl98B4Pyzx7LsxvOZONo1gMlA8t3/+m+e2xugbGwtAFXWTs52neK3j/xP/KhkCaLk7yd2V008vwpmF5rEZy6lliAllSwtZSy2zv59yr6ZoNpBKYYalSRRQCxmE8tuPJ9/WXoxbqeV1946xRfv/T/u/PaDmO2e+HHGmFQ6TjT7aHj8NV7zTqVsbC1aOEDLG3/i53d/kq9+OTURwihjlPz9zJRO7Gb4kvjMhTClJFIYY0i6ccq2GDZdMsSaNWu59dbbBpiePxTCvApFKioGNYRcPHsisyaPYuPfdvHUtkMci1YyadGdVNn8jLJ18pvHf4umSYSItQG/9bZ/5nizj6Z3zrBl93Fej7f1kPiO7aR93z+IBn1pM/HWrl2Xsc2CWhGXCj27nOXL7zAofkBiDKk/xbDptQobGxupq6uLu5mzLcQurZR+RamjDNQQU+lxcOdHLyRw7FX2Nls5E3LTGnbRGnYx8fLP8lKLRjQSxN5q4y/ffCLhszaLiUpzO68+9T98/Mb3surBf8STJSC1kd327dviCRKqgr+0yCQgmzjx9xir7FU++kqGSDxv7++JaneuyC8qBpVXJC0dfv7jhxtpOtCG1T0Wi7Mi/lOTSTCmwsmsqdWcN2MMF8+ewK8feiBl4opp9cGDDz6AcRW7d+8eamtnqdVtyTHQbq/9iR/mKhmilIR5808kEiUciaJpWsw7QuzvWkD88UgZKyXRpEQIgTCZsFrMWC1mgwL98EYlSRQ1PZOFMFl46ulniGrgcVoxm9PFi9IHuNNPTirpoTQZyok/1+dWmX8QM0bBUASpaVgsJmxmgcNmocxhwWazYMvC4IS7DVqnP4w/GCEU0QhHNITZhN1qwWIx5+lu8otKkihaEoPbUovQ8NMfUem2pTFOegKFSPN1puQGlfRQmgxdA8HcahWO3HbnUkq6AkF8Pj/BQACnBaaOdXPO1FHMrKli6oRKxlW78bgc2K2WrHZDVouZMoeNMVUupoyv4OyaKs6dVs20sW6cFggGAvg6/fiDIYZ5P9i0KAOVV7KfLJLbcRgzr4amdbaisAzVxJ9L8dZSFOYdHJom6fQH6ezyI6MRJlQ5OWfqKM6uGcX40R6cDuuQuOUcdivjR3s4u2YUsyaPYlylAy0Sxtfpp8sfHDHGSiVJ5JWeycLYjhvoljPq2Smla8ehJ0fU1dWhXHfDiaFtUZG71hrF1El56JBSEgiG0KIaZXYLU0a7KHPaCjYek0lQ4XJQ4XIgpcTbFaS1I0BnKIrJbMJptw3b2JWKQRWADRs24PP1JD7ohZmJyQwyIWNPR8/cGy6TgSJGaaVv9ycmWjqEwhFCwTAOm4nqCiflLkehh9QrUko6uo2VPxjFbDXjtBfOkPaHbGNQageVd2SCRM2qVSsTVAKMMabVq1elGChlnIYnQ9dAcChIHFPv7eWLzbgmIqXEHwiDjFLptjN6bFWGRKXiQ4jEnVWHL0hzh59gWMNsNeOwDY37MZ8oA5V3shXplN3xp0TWrFmrjNSwpRQTXEqzNiqqaXT5Q9gtgonVZXjKsu3hVpwIIajwOKjw9BirVl+AQCgKwoTDYcVsKg3Da0S5+PJOTzp4plRdvTWGsSAXSCnQLcY/fMVIpHRqo3Q3nstpYfwoNzbr8EzjNtLpD9LqDdAVjBLVJGazGbvNislUuN+NcvEVIXqcIZ0kUc/OKNbvZ9++vdTV1VFXVxf/PkBTUxPl5R6K7Q9fMZKJeQWMBqrYjFMgFCYSjjDKbWfMuFEFnZzzjctpx+WM7RD1rMT2ziDBoEYkohEFTDkuFpZSEo1qRDVJVIsiNdCkhpCxpA9Ni2rZnEcZqLzR4wrRJYlqa2eyd+8+gIRYU2Njo0G1uqfmqcdQjZw/LkUpkD5FvhiMlD8QBE2jusLJqPLyko/JDBaTSeBxOfAYEkCiUY1gOEKnP0wgFOk2KhJNk3HVC0lMAcOIEN3R8m6lDJNJIITALGL/XWY1Y7VasFsdWCwmLCZTPL4XDQe7shnvoAyUEOI7wIcADTgF3CqlPDaYcw5fEmNPAHv37ktw3+lGKrN7ZGBV/6WeXaUoZoY2RX6gdPlDmNAYV1lGhae4s/EKjdlsosxso8xRfBmAg4pBCSHKpZQd3f/9/4DZUsrP9PW5kR6DSo09MSTSMaWVujxwpIyt+KSMrfp0hIgFj80m04hy6eSbYnnPYkoPIcxCMq6qrOjTxEcyQohtUsoFfR03qB2Ubpy6cTEcS8lzSqorJF2mnt46I7GPU393PqWZXZWJUDhCOBxFkxomwGwSmEwi/v9OswmTWWA29fjQY4ZLIxSOoEUg0u26iEY1NAQWswmbzVKS2U3FRKFT5KWUdPmDmE1QM9qF21naGXmKHgadxSeE+C7wKaAdeLeU8nSG45YDywEm1Uy+8G/P7cDhsGEdpmKIqaS6QnprnaG3bTc2p+v/irR0squMhCNRQt0inFaLCavFRJndgttpw2Y156ROJRLV6AqE6OgMEgxrhMIaJsvwrsofbkgZC/hbzYKJowqr9qDoH9nuoPo0UEKI/wPGp/nR16WUjxuO+yrgkFJ+s6+LLliwQL7yymZOtnbS5gvhdNiGrWqvkXSukOTWGVJqLF26NG1fp4EZl+JXno5ENYLBMEgNu9WMy2mhwuXAbstfDo+UEl9XkBZvEH8wgjCZcDqUsSpGYploAewWExOrXTiLMHai6J2cGah+XHAq8EcpZV1fxy5YsEBu3boViGWQHD3jxReI4HI6RkCsoG+ZGF36aPA7n+LdQYUjUQLBMFYTuJxWRnkcOOzWgo7JiK8ryJl2P/5QFKvVmldjqUiPniLttAomVLuL6n1R9I9sDdSgfCVCiJmGL28E3uzvOcxmE1PGVTBjYgXhcAh/IDSYIZUA6TLzEr8nhCkHrReKT3k6qmn4ugL4/QFcVsHMSRXMnDyKiaM9RTfZuMvsTJtQSW1NFeVOE51d/ljKsiLvRKIa3k4/RMPMmFjO9IlVRfe+KIaGwS4Lvy+EmEUszfwg0GcGXybsVgtnT6qizRvgeLMPh8M+Itx+PST2fUpOnuh/XUnxKE/Hq/cdFs4aX15SuxGz2cTYKjdjq9y0+wKcau0iiqBMuf+GnHAkSjAYwuWwMHNS5QiKVyt0CiJ1ZHTxpSMa1Thy2os/rMUroIczibEp4skTdXV1PPjgA4OOQRWqDsofDKFFo4zy2Bld4Ro27ltvV5ATLZ1EpRgR72e+CYUjhEJhKlw2xlW5Ska8VZE9eUkzHyrMZhNTx1fQ5gtw7IyPsjLHME4F7kkHl1KyevWqeA+oOXPmAAxy55N/AdJAKEw0EmFsZRmjyp1Dfr184ymz4ymz4+0KcvyMDylMlClDNWj8gSCaplHtcTB6vEftUBXFuYMyEolqHDjeBqaYwOHwJJbNZ2xOuHjxYoQQlJd7WL78DqTUBlkXNfToLpnR5Q5GV5aNmAmm3RvgeEsnZqsFx7B9R4cGTYu1UbeYYGxlGRVuVVw7EshLkkQ+sJhNnF0zCrfNRGdXoNDDGTJiXXJ7EELQ2NiI1+tjw4YN3eKy+mIilgBRLK3fpZR4OwPYTBqzJo9iTJVrxBgngAqPg1lTRlHptOD1dRGORAs9pKInFI7g6058OGt8OTNrRinjpEih6A2UzvjRHiaMKqPD56cQu75809jYSH19PStXrog3ONSz7/SYlNfro9DiHaFwhIA/wFnjPdSMrRg2cab+IoRgTJWLc6ZU47CAz+cnqmUl2Dxi0BUfOrv8uG0mamuqmDqhsqSSZhT5paTejAq3A6fdytvHWnE4HFiGRfBUxns/1dfXp3TQFSLbBof5x9cZoMplZdz4USNqx9QbJpNg4mgPYyqjHDvtxeeP4ipzjOjno7t+Y4W1pd8cUJE/Sm6Gt1nN1NaMIhoOEQpHCj2cHBBLB9dVzY30xKREDuqickckEsXX6WfqWDfjR6tgdjqsFjNTJ1QybZyHYDBIl3+41/clohfVdnX5KbPAzJoqzppUpYyTol+UnIGCWJbfjElVWIUkECz9P/zly5cDPX2gtm7dQn19PU1NTaxduy7eYddIoYpt/cEQMhph1uRRSvssC5wOGzNrRjGhykFXpx//MHhfeyMQCuPzxWJLU0a7mDWlmvGjPcPE26HINyXl4jMihGDK+ApOnPHS5g+WeD1KLFvP6LZbvXoVQgjcbndc9qjQ/XZ8XQFGuWyMq3bn5XrDiXK3g3K3g5YOP6dau7BYLcMmK1WvW7JZTIz2OKj0qMaAitxQ9Gnm2XC6tZMz3iDuslLPAkpfVFvofjtSSnydASaOdlGpMq0GjZSSM21dnOkIxDqOlqChisWVwlhMUOG2UV1epgpqFVmTd7HY/pBrAwXQ5g1wrKUroZXx8KI3RYihU4vQNElnl5/p48uVanSO0Q1VizeAMJmL/vmGwhGCoTBWk6DcZWNUuVPJDykGREkrSQyESo8Dk4AjZzrxuIefekEmRYih3F1FNQ2/P8DZk6qwWdVElGv01PQxVS7avQFOtXcRjkrKnPaiUE6RUhIMR4iEIlgtMaM0dWyViicp8sawMVAQ8/NPNgkOn/INUyOVzNB1zY1ENUKBIDNrRqkJKQ9UeBxUeBwEQxFOtnbi6wxiMue/gWKsWWQYkDisZsZ6HHhc5SO2vk1RWIaVgYKYTtqUcXDopBePu6zQwxlihqZGKhKJEg6FOLumSsUV8ozdZmHKuIqYOkdXkNaOAF3hKBBrT59Ll5qUklA4SjgSQUiJzWrG7bRQUV2O3TrspgZFCTJsYlDJdPpDHDzpHTE7qVx1ze0xTqPUqrmICIYitHkDdAYjRCIaEU1iMpswm0xYzCbMZlPanZamSaKaRiSqoWkaWlTDLMBiMWE1m3CX2XA7bcqFq8grIy4GlYzLaWPqOA8HTnopH9ZGSqatkRrIDioS1ZRxKlLsNktCer+myZjMVDBCMBwlFImgSUDKuGNXCIHZBE67GbvVjt1qxma1qN+tomQYtgYKYkZq2jjPMN5JyZReUQOtkYpENcLBoDJOJYLJJHDYraqzrGJYM6wNFPTspA6e7BiGMancdM2NahrBQJCZNVXKOCkUiqJh2MagkvH5gxw6OVyz+wZeBxXTTPNTq7L1FApFnhg2/aByhdtpZ8o4N16fv9BDGQIG1jVXSklnp5+Zk1Rti0KhKD5G1KzkdtqZPGa4Gqn+4+sMcNakCqUGoFAoipIRZaAAPC47NaNdeDtHtpHy+vxMHedRLcoVCkXRkhMDJYRYLYSQQojRuTjfUFPudjCx2oW3c/i2kO8Nb2eAmtEuXKpdhkKhKGIGbaCEEJOB9wKHBj+c/FHpdjBxVBm+EbaT6vIHGV/lpFypkisUiiInFzuo9cCXKUT3vEFS6XEwvmrkGKlAMERlmZVR5cMxk1GhUAw3BmWghBA3AkellK/laDx5p6rcyfiqMjq7hre7LxgK47QK1WxQoVCUDH0W6goh/g8Yn+ZHXwe+BlybzYWEEMuB5QBTpkzpxxCHnqruHcXJNj+ukm96mEo4EsEiJDVjKws9FIVCociaARfqCiHmAn8Hurq/VQMcAy6SUp7o7bOFKNTNhtYOPydau3C7ho8LLBKJEo2EmTGpSrXhVigURcGQi8VKKXcCYw0XPAAskFKeGeg5C42+kzrR5h8G7eNjEka6+KsyTgqFotQYcXVQfVFV7mRCVVnJ10lpmiTgDzBjktLXUygUpUnOxGKllNNyda5CU+rt4zVN0tUVUA0HFQpFSaNmrwyUux1MGeemw+enEIK6A0XX15sxqVJJGCkUipJGGahecDvtzJhQjtfnR9OK30hJKfH6/JxdU6U6pCoUipJHGag+cNit1E6uwu/3E4lqhR5ORqSU+Hx+ZkyqUMZJoVAMC5SBygKrxczMmlFEwyFC4Uihh5OCpkl8nbGYkxJ/VSgUwwVloLLEbDYxY1IVVpPEHwgVejhxopqG3+9nZk2l2jkpFIphhTJQ/UAIwZRxFYxyW/EWQfJEKBwhFAgys2aUSohQKBTDDmWgBsDoShfTx3vo7AoQiUQLMoZAMIQZjZmTR6lUcoVCMSxRM9sAcTpszJo8CjManf5gXq/t7QxQ7rAwbUKlUohQKBTDlpwV6o5ETCbBlPEVtPsCHGvuxG63DamrLRKJEgiEmDLOjdtpH7LrKBQKRTGgdlA5oMLtYNbkUdhM2pDUTEkp8XUFMAuNWVNGKeOkUChGBGoHlSNMJkHN2ApC4SjHznjx+qOUOe2YTYNbA3T5Q5iIMm2sB6dDtWhXKBQjB2WgcozNambahEpC4SgnWnz4OoOYzWYcdmvW8aJIJIo/GMZqhonVLjxlasekUChGHspADRE2q5kp4yqQUtLhC9LmCxAIa0Q0iclkwmQ2Yeo2WFFNQ2oamiaxmgUep5VJ1UoRQqFQjGyUgRpihBBUeBxUeGL9pTRNEgxFCEejRKISE2AyC+w2CzaLWWXlKRQKRTfKQOUZk0ngdFhxoiSJFAqFojdUFp9CoVAoihJloBQKhUJRlCgDpVAoFIqiRBkohUKhUBQlykApFAqFoihRBkqhUCgURYkyUAqFQqEoSpSBUigUCkVRogyUQqFQKIoSUYi25UIIL7An7xcuHkYDZwo9iAIy0u8f1DNQ9z+y73+WlNLT10GFkjraI6VcUKBrFxwhxFZ1/yP3/kE9A3X/6v6zOU65+BQKhUJRlCgDpVAoFIqipFAGqqFA1y0W1P0rRvozUPc/ssnq/guSJKFQKBQKRV8oF59CoVAoipKCGighxBeEEHuEELuEEP9VyLEUCiHEaiGEFEKMLvRY8okQ4gdCiDeFEK8LIX4nhKgs9JjygRDiuu53/i0hxFcKPZ58IoSYLIT4hxBid/ff/F2FHlMhEEKYhRA7hBB/KPRYCoEQolII8Wj33/9uIcSlmY4tmIESQrwb+BBwnpRyDrCmUGMpFEKIycB7gUOFHksBeBKok1KeB+wFvlrg8Qw5Qggz8GPgemA2sEQIMbuwo8orEWCVlPJc4BLg8yPs/nXuAnYXehAF5F7gL1LKc4Dz6eVZFHIH9Vng+1LKIICU8lQBx1Io1gNfBkZcIFBK+TcpZaT7y5eBmkKOJ09cBLwlpXxHShkCGokt0kYEUsrjUsrt3f/tJTYxTSrsqPKLEKIG+ADw80KPpRAIIcqBRcAvAKSUISllW6bjC2mgaoErhBCvCCGeEUIsLOBY8o4Q4kbgqJTytUKPpQi4HfhzoQeRByYBhw1fH2GETdA6QohpwAXAKwUeSr65h9iiVCvwOArFWcBp4IFuN+fPhRCuTAcPqZKEEOL/gPFpfvT17mtXEdvqLwQeEUKcJYdRWmEf9/814Nr8jii/9Hb/UsrHu4/5OjHXz6/zObYCIdJ8b9i879kihHADvwW+KKXsKPR48oUQ4gbglJRymxDiqgIPp1BYgPnAF6SUrwgh7gW+Anwj08FDhpTymkw/E0J8Fnis2yBtFkJoxPSpTg/lmPJJpvsXQswFpgOvCSEg5t7aLoS4SEp5Io9DHFJ6+/0DCCFuAW4A3jOcFia9cASYbPi6BjhWoLEUBCGElZhx+rWU8rFCjyfPXAbcKIR4P+AAyoUQG6WUSws8rnxyBDgipdR3zo8SM1BpKaSL7/fA1QBCiFrAxggRT5RS7pRSjpVSTpNSTiP2S5s/nIxTXwghrgP+BbhRStlV6PHkiS3ATCHEdCGEDagHnijwmPKGiK3GfgHsllKuK/R48o2U8qtSypruv/l64KkRZpzonuMOCyFmdX/rPcAbmY4vlFgswP3A/UKIJiAE3DJCVtGKGD8C7MCT3bvIl6WUnynskIYWKWVECHEn8FfADNwvpdxV4GHlk8uAm4GdQohXu7/3NSnlnwo3JEUB+ALw6+5F2jvAbZkOVEoSCoVCoShKlJKEQqFQKIoSZaAUCoVCUZQoA6VQKBSKokQZKIVCoVAUJcpAKRQKhaIoUQZKoVAoFEWJMlAKhUKhKEqUgVIoFApFUfL/ATpVqCD/zHeDAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "m.plot()" ] }, { "cell_type": "code", "execution_count": 12, "id": "636332d8", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Optimization restart 1/10, f = 195.84754738377032\n", "Optimization restart 2/10, f = 195.84754738420352\n", "Optimization restart 3/10, f = 195.84754738488476\n", "Optimization restart 4/10, f = 195.84754738373357\n", "Optimization restart 5/10, f = 195.84754738844853\n", "Optimization restart 6/10, f = 195.84754738371487\n", "Optimization restart 7/10, f = 195.847547383823\n", "Optimization restart 8/10, f = 195.84754738461945\n", "Optimization restart 9/10, f = 195.84754738372928\n", "Optimization restart 10/10, f = 195.84754738422166\n" ] } ], "source": [ "_ = m.optimize_restarts()" ] }, { "cell_type": "code", "execution_count": 13, "id": "2627c258", "metadata": {}, "outputs": [ { "data": { "text/html": [ "\n", "\n", "

\n", "Model: GP regression
\n", "Objective: 195.84754738371487
\n", "Number of Parameters: 3
\n", "Number of Optimization Parameters: 3
\n", "Updates: True
\n", "

\n", "\n", "\n", "\n", "\n", "\n", "
GP_regression. valueconstraintspriors
rbf.variance 3.173279047762375 +ve
rbf.lengthscale 2.3422120232188184 +ve
Gaussian_noise.variance0.25032019995651095 +ve
" ], "text/plain": [ "" ] }, "execution_count": 13, "metadata": {}, "output_type": "execute_result" } ], "source": [ "m" ] }, { "cell_type": "code", "execution_count": 14, "id": "4a7e5a09", "metadata": {}, "outputs": [], "source": [ "k = GPy.kern.RBF(\n", " input_dim=1, variance=true_kernel_variance, lengthscale=true_kernel_lengthscale\n", ")\n", "m = GPy.models.GPRegression(X, (f + noise).reshape(-1, 1), k)" ] }, { "cell_type": "code", "execution_count": 15, "id": "100be479", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Optimization restart 1/10, f = 195.8475473936342\n", "Optimization restart 2/10, f = 195.84754738474803\n", "Optimization restart 3/10, f = 195.84754754618757\n", "Optimization restart 4/10, f = 195.84754738374266\n", "Optimization restart 5/10, f = 195.84754738772745\n", "Optimization restart 6/10, f = 195.84754738398067\n", "Optimization restart 7/10, f = 195.84754738390444\n", "Optimization restart 8/10, f = 195.84754738371876\n", "Optimization restart 9/10, f = 195.8475473837782\n", "Optimization restart 10/10, f = 195.8475473893405\n" ] }, { "data": { "text/plain": [ "[,\n", " ,\n", " ,\n", " ,\n", " ,\n", " ,\n", " ,\n", " ,\n", " ,\n", " ]" ] }, "execution_count": 15, "metadata": {}, "output_type": "execute_result" } ], "source": [ "m.optimize_restarts()" ] }, { "cell_type": "code", "execution_count": 16, "id": "ac762d83", "metadata": {}, "outputs": [ { "data": { "text/html": [ "\n", "\n", "

\n", "Model: GP regression
\n", "Objective: 195.84754738371876
\n", "Number of Parameters: 3
\n", "Number of Optimization Parameters: 3
\n", "Updates: True
\n", "

\n", "\n", "\n", "\n", "\n", "\n", "
GP_regression. valueconstraintspriors
rbf.variance 3.173275150587815 +ve
rbf.lengthscale 2.342210766061685 +ve
Gaussian_noise.variance0.2503201961365106 +ve
" ], "text/plain": [ "" ] }, "execution_count": 16, "metadata": {}, "output_type": "execute_result" } ], "source": [ "m" ] }, { "cell_type": "code", "execution_count": 17, "id": "65880e12", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 17, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAagAAAEYCAYAAAAJeGK1AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAB2cElEQVR4nO2deXiU1b34P2f2NQth30FZlIDIpraKrQVrN7vZGlqsy61oWxUDtLd7ue3tbWtZStdL7E9tixJ6W21ta6toVdwJKEJcQEEEEiB7MjOZ/T2/PybvZGYyk0ySSSaB83kenjAz75z3vMuc7/vdhZQShUKhUCiGGoZ8T0ChUCgUinQoAaVQKBSKIYkSUAqFQqEYkigBpVAoFIohiRJQCoVCoRiSmPKx05EjR8qpU6fmY9cKhUKhyDN79+5tkFKO6mm7vAioqVOnsmfPnnzsWqFQKBR5RgjxbjbbKROfQqFQKIYkSkApFAqFYkiiBJRCoVAohiR58UGlIxwOc+LECQKBQL6ncsZjs9mYOHEiZrM531NRKBSKjAwZAXXixAncbjdTp05FCJHv6ZyxSClpbGzkxIkTTJs2Ld/TUSgUiowMGRNfIBCgpKRECacBRghBSUmJ0lQVCsWQZ8gIKEAJp0FCnWeFQjEcGFICSqFQKBQKHSWgEjAajcyfP585c+ZwwQUXsGnTJjRN6/Y7R48e5YEHHhikGSoUCsXZw7AVUKmNFnPReNFut7Nv3z5ee+01du7cySOPPMJ//dd/dfsdJaAUCoViYBiWAmr9+vWUl5fHhZKUkvLyctavX5+zfYwePZqKigp++ctfIqXk6NGjXHbZZSxYsIAFCxbw/PPPA/D1r3+dZ555hvnz57N58+aM2ykUCoUColENg8lqz2bbIRNmni1SSlpaWtiyZQsAmzdvpry8nC1btrB69WqklDkLApg+fTqaplFXV8fo0aPZuXMnNpuNt956ixUrVrBnzx5+/OMfs2HDBv7+978D0N7ennY7hUKhONsJR6IcrmnBaDRmJXuGnYASQrB582YAtmzZEhdUq1evZvPmzTmPUNO1tHA4zG233ca+ffswGo0cOnQo7fbZbqdQKBRnE/5AiHdOtuF02gGy8skMSxNfopDSGQjhdOTIEYxGI6NHj2bz5s2MGTOGV199lT179hAKhdJ+J9vtFAqF4myhzRvgnVNtuFx2DIbs1+lhKaB0n1MiiT6pXFBfX8+tt97KbbfdhhCC1tZWxo0bh8Fg4A9/+APRaBQAt9uNx+OJfy/TdgqFQnE2Ut/so6bRh9vl6LUSMewElC6cdJ+TpmmsXr2aLVu29FtI+f3+eJj5smXLuPLKK/ne974HwJe//GV+97vfcfHFF3Po0CGcTicA8+bNw2QyccEFF7B58+aM2ykUCsXZxom6Vpq8YVzOrGIiujAsfVBFRUVJPifd3FdUVNQvM1932s6MGTPYv39//PWPfvQjAMxmM0888UTStum2UygUirMFTZMcPdmCJgw47JY+jzPsBBTEwswTo/V0IaVK+CgUCkV+0SP1rFYLVpOxX2MNSwEFXevJKeGkUCgU+cXrD3LstAeno3fBEJkYtgJKoVAoFEOHhhYfda0B3C5HzsZUAkqhUCgUfUZKyfG6NgJhibuPwRCZUAJKoVAoFH0iHIly9GQLRrMZuy33HbqVgFIoFApFr/G0Bzlelzt/UzpykgclhLhHCFEnhKjOxXj5QgjBddddF38diUQYNWoUH/3oR/M4K4VCoRhanG70cqLei9vlGDDhBLlL1L0PuCpHY+UNp9NJdXU1fr8fgJ07dzJhwoQ8z0qhUCiGBpGoxtsnmmgLRvucfNsbciKgpJS7gKZcjJVvPvShD/GPf/wDgO3bt7NixYr4Zz6fj5tuuonFixdz4YUX8te//hUgY4uNp556ive9731cc801zJ49m89//vM5LcekUCgUg4XHF+TQ8SZMFgs2S+79TekYNB+UEGIVsApg8uTJ3W67eNX9AzKHqorP97hNWVkZ3//+9/noRz/K/v37uemmm3jmmWcA+OEPf8gVV1zBPffcQ0tLC0uWLGHZsmUZW3EAvPLKK7z22muMHz+e9773vTz33HNceumlA3J8CoVCkWuklNTUt+EJRHMaQp4NgyagpJQVQAXAokWLhqwaMW/ePI4ePcr27dv58Ic/nPTZY489xsMPP8yGDRsACAQCHDt2jPHjx2dssbFkyRImTpwIwPz58zl69KgSUAqFYlgQCIY5eqoNs8WMy2Eb9P0PySi+bDSdgeTqq69m3bp1PPXUUzQ2Nsbfl1Ly5z//mVmzZiVtv379+niLDU3TsNk6L6TVao3/32g0EolEBv4AFAqFoh9IKalr8tHoDeJy2PJWqWfYVTMfDG666Sa++93vMnfu3KT3P/jBD/KLX/wi7kd65ZVXANViQ6FQnDkEgmHeOt6EN6ThdtrzWkYuV2Hm24EXgFlCiBNCiP/Ixbj5YuLEiaxevbrL+9/5zncIh8PMmzeP0tJSvvOd7wCZW3EoFArFcEH3NR051YbNbsNizr+BTeQjqmzRokVSDyLQeeONNzjvvPMGfS5nK+p8KxQKnVZvgNpGHxaLeVAE0zkTRx0MBXyze9ou/yJSoVAoFHkhFI5yvK6NiCYGJa+ptygBpVAoFGcZ0ajGyUYPbf4ITrsNywBWg+gPSkApFArFAJDYVDXd63ygaZK6Ji/N3hA2myXn1cdzjYriUygUihyzfv16ysvL4xG/UkrKy8tZv359XuajaZLTjV4OHm/CF5a4XHZM/ex2OxgoAaVQKBQ5REpJS0sLW7ZsiQup8vJytmzZQktLy6CWO4tGNU41eDh4vAlPUMPltA+J6LxsGT4zVSgUimGAEILNmzcDsGXLFrZs2QLA6tWr2bx586CY+YLhCKcbffiCEaxWy5AMgMiGIRtmfryulUBIy9k+bRYDk0YXdrvNqVOnuPPOO6mqqsJqtTJ16lR+9rOfMXPmzF7t65lnnuHWW2/FbDbzj3/8g9WrV/OnP/2py3bve9/72LBhA4sWLerV+LlAhZkrFAOLlBKDodNIpWnagAonKSVt3iB1re1EohK73YrRMDSNZMM+zDwQ0rDZrD1vmO14gWC3n0sp+eQnP8n1119PZWUlAPv27eP06dO9FlD3338/69at48YbbwRIK5wUCsWZi27WS6S8vHxANKhAMExdczu+YASjyYjdNvg18waKoSle88CTTz6J2Wzm1ltvjb83f/58Lr30Ur761a9SWlrK3Llz2bFjB5C5lcZvf/tb/vjHP/L973+fz3/+8xw9epTS0lIA/H4/ZWVlzJs3j2uvvTbedwpihWgvueQSFixYwGc+8xm8Xi8AU6dO5Xvf+x4LFixg7ty5vPnmmwB4vV5uvPFG5s6dy7x58/jzn//c7TgKhWJwSPQ5rV69Gk3TWL16dZJPqr8EQxFO1LVy6HgjR097kQYTLqcdu9WSgyMYOgxZDWqwqa6uZuHChV3ef/DBB9m3bx+vvvoqDQ0NLF68mKVLlwLpW2l88Ytf5Nlnn+WjH/0o11xzDUePHo2P9Zvf/AaHw8H+/fvZv38/CxYsAKChoYH//u//5vHHH8fpdPKTn/yETZs28d3vfheAkSNH8vLLL/PrX/+aDRs28Nvf/pYf/OAHFBYWcuDAAQCam5t7HEehUAw8QgiKioqSfE66T6qoqKhPGpSUEp8/RFNbAH8ogkRgs1mw2wenL1O+UAKqB5599llWrFiB0WhkzJgxXH755VRVVVFQUNDrVhq7du3ijjvuAGJtPebNmwfAiy++yOuvv8573/teAEKhEJdcckn8e5/61KcAWLhwIQ8++CAAjz/+eNwUCVBcXMzf//73bsdRKBSDw/r165PynnQh1RvhFAxHaPEE8PrDhMIaBpMBm8WCw3FmC6VElIDqYM6cOWl9Rd2p431ppZHuBpVSsnz5crZv397tfhL3kS7pr6dxFIqzkXwlzKbuo6d9BsMRWr0BvO1hwhENKUSsNp7ViiV37vhhhfJBdXDFFVcQDAa5++674+9VVVVRXFzMjh07iEaj1NfXs2vXLpYsWdKnfSxdupT77491C66urmb//v0AXHzxxTz33HO8/fbbALS3tyc1PUzHlVdeyS9/+cv46+bm5j6No1CcyQy1hFkdKSX+QJhTDR6O1Lbw5ruNvHOyDU9Aw2y14nDacTpsmIdBMu1AMmQFlM1iIBAI5uyfzdL9oQoheOihh9i5cyfnnHMOc+bMYf369Xzuc59j3rx5XHDBBVxxxRXcddddjB07tk/H9KUvfQmv18u8efO466674oJu1KhR3HfffaxYsYJ58+Zx8cUXx4MhMvHtb3+b5uZmSktLueCCC3jyySf7NI5CcaYylBJmNU3i8QU4UdfK2yeaefPdJt6t8+KPgNliwdkhkIZTEu1gMGTzoBQDizrfirOBRKGkMxgJs5om8fmDNHuDBMMakUinD8kwRAuzDibZ5kENWQ1KoVAo+ktiBJ3OQAmnQDBMbYOHt443cfBYE7XNfqQwYbNZcbnsOGxWJZx6iRJQCoXijCVTwmwuLEeaJmn1BDh6MuZDOnraSzAqsNltSiDlCCWgFMOG1EUlH+ZpxfBhIBJmNU3S3ObncE0zB481cbotgMFk7vAhWTEZ1ZKaS5RHTjEsWL9+PS0tLXHzjL74FBUV5T0iSzE0yVXCrJQST3uQxlY/wbCG0WzCdhaHfg8mSkAphjyJ0VgQ8yEkPhkPhUZwiqFJfxJmQ+Eodc1evP4IwmjAbrXgtKj7bDBRAkox5BkK7QsUw5dsEmYThVirJ0Bdi4+IBna7Fafz7KncMNRQBtMETp8+zec+9zmmT5/OwoULueSSS3jooYf4yle+wvz58zn//POx2+3Mnz+f+fPnd6k8UV9fz0UXXcSFF17IM888k6ejODMZzGgsxdnF+vXrufPOck41enjz3UZOt/j51a9/zQP3/37Itqs4W1BnvwMpJZ/4xCdYunQpR44cYe/evVRWVnLixAl+9atfsW/fPh555BHOOecc9u3bx759+7jmmmuSxnjiiSeYPXs2r7zyCpdddlmejuTMZCCjsRRnL6FwhFNNPv7w53/xw7t+htNp49e//iXbt2/H4/EC3d1fqZ+pezHXDDkTX8XeigHfx6qFq7q89+9//xuLxZLUbmPKlCncfvvtWY25b98+vva1r+H3+5k/fz4vvPACdvvw7GI51EiNxkr0QYHSpBS9JxSOcrLRiy8Q4T+/9V1szgK2b9/OH3fE6liuWLGCtWvXAOnvq4qKrXg83oRtJBs3bsLtdrFq1S2DdhxnOkqD6uC1116Lt7/oC/Pnz+f73/8+1157Lfv27VPCKYdkisZavXp1n9sXKIY/fUk7CIWjvHuyhbdrW5AGI26XHZPR2CFoOulOOIHE4/Gyfft2Nm7chC6cstO6FL1hyGlQQ4WvfOUrPPvss1gsFqqqqvI9nbOeXLQvUJw59DbtIByJUtvgxReM4LBbcTsTG/vJDkHTycaNm7oRUiIu0LZv3x7vHtCT1qXoPUqD6mDOnDm8/PLL8de/+tWveOKJJ6ivr8/jrBSJ9LZ9geLMpDdFYKNRLVagtaZDY3LaUwIfOrWfFStWsGdPFStWrEjSjtIjeql1KfpCTgSUEOIqIcRBIcTbQoiv52LMweaKK64gEAjwm9/8Jv5ee3t7Hmek6C+q8sSZSaKJd8uWLRgMhiT/pK5RnW70cuh4M2HNgKuLYIqPhtvtStB+YoJmxYoVuN2ubmaRXuvqFGgqgCIX9NvEJ4QwAr8ClgMngCohxMNSytf7Ml66AIbBQAjBX/7yF8rLy7nrrrsYNWpUvG26YvihKk+c2ehCKrFKuX6tm9v8nG5ux2wx43L17AuOBTVIKioq4oEPurBKH/iQrHWtXbsm/hrA5XLh9aoAilyQCx/UEuBtKeURACFEJfBxoE8CKp+MGzcuqY16KlOnTqW6ujrj5zfccAM33HDDAMxMoROORIlENMLRKFLr1IoMRoHJaMBoNGA2GlTliTOcdGkHt68u5/Z138FgMOJ09j5ISQ98AJKEzooVK4hpQPo9k6p1dZr7dOGU3TiKnsiFgJoAHE94fQK4KHUjIcQqYBXA5MmTc7BbxZmMpknaAyFavAFCEUkkohGVEoRAIBAGgUGIuHakSRkTVlKiaZJby7+DXzrYet+9/GrrPURDflV54gwhNe3gpxs2ctuab7Ljr3/H6Chh3bo1PQ/ShZ4CH5LRta5EoaULKym1bsZR915vyIWASnfGuxhcpZQVQAXEGhbmYL+KM4xIVKOprZ02X4hwVGIwGrBazJgtBsyWnr+fyre/8y3+8te/gChBk3DLnd/mSG0LBQ4LxQV2VXl6mKKnHdxxx2q+/p0f8HZNC9/49rexO10dfqO+CoGYkNGFChAXThs2bKSgwJ1gokunCYl4ftSaNeVJ47hc/ZnX2UsuBNQJYFLC64lAbV8GUqaXwWEoBQtIGeup09DmJxSVWC1mrDYb/S8UneDElhoGYOvW37B27Rpa/VHq21owG8BtN1NS5MBsMvZ7j4rBY826r1PT4MUXkrg6zHn911AkGzZsTHpnw4aNHDhwgNdee42ysjL0Z+9MvindTPjyy3uTxtm162lWrboZIdRDUW/IxdmqAmYIIaYJISxAGfBwbwex2Ww0NjYOqcXzTERKSWNjIzabLa/z0DTJqQYPB483UecJYrXZcDvtWMy5eGbqPnTYajHidtqw2W34I/BWTStvn2iisbUdTVP331AmHIlyuKaZ2mY/Lpcj5X7pv3CqrKykrKysQxhBZWUlr732WtKWmZNyBWvWlDNr1kwOHjwUf1d/vWnTZlQ0X+/o92ogpYwIIW4DHgWMwD1Sytd6+FoXJk6cyIkTJ1Te0SBgs9mYOHFiXvataZJTTV5avSEsVjNOx0BU3MjsxE41AZlMRtwd2lOzL0xdsx+nzcSYEU6sFpXHPlTQw8abfSEcdhuWHHeqraiooLq6mrKyMtatW4uUMilgqqysjMrKyvh7yT6lTnOfEAaWLr08SUBt27aNTZs299P8eHaSk1+glPIR4JH+jGE2m5k2bVoupqMYgkgpaWhpp6EtgNVqySr8tz9058TOhNVixmoxo2mSI6faMBlgdJGDQld+tc2znTZvgJpGHxaLOW7Oyy0x01x1dTWlpaVIKVm5cmW331izppz0IeQSr9ebtO2mTZtZs6Zcmff6gDpjigHH6w9y6HgTrf4ILqd9EP09XZ3YPSMxGAQuhw2bzcapFj8HjzVS3+xT5udBJhROMOflzPybDhFPzq2srGTx4iUcPHiIWbNmUlW1m2uvvbZL+snKlSuRUksy9yW+TjUtK/Ne31A2DMWAoWmSE/Vt+IJRXANiysuG1GirzHko6SpU/+ZXv8TtdnH9DTfR2NZEkdPC6BEuDDk2MSk6kVJyqtFLywCZ89LTNYJv27ZtVFTczb59rwBQWlrKPffcw5VXXsnBg4dYvHgJkGzuy9a0rMgOJaAUA4KnPcjxOi92mwWXIz8dSXvXEkF2m6hptZiwWsz4whHePNZEkcvCmGInRhWqnlNaPQFqG31YbQNlzstEptJFxLWpe++9h02bNtPS0pK0XaLpuC+mZUVm1K9LkVOklJyoa6WmwRdrZZC38O3etkToNPNs376dRYsWJ5Wy0RcYi9mE22UnEIFDx5s5UddKNKoN8rGdeQTDEd4+0cTp1gAulx2zaTCfndNHfepmvbKysrjGtH37dmbNmpn07a41+PpiWlakQ+TDrr5o0SK5Z8+eQd+vYmAJhaO8U9uCyWIeQH8BZG+2k0k10iCbjH7JokWL46/27KnqZluIRKL4AyEKHCbGlbiVRtVLolGN2gYPnkAEp92WN9Np99r2qqR7AmL3kcvlYteupzl48BArVqxgzZryeLSeqrnXPedMHHUwFPDN7mk79WtS5IRWb4C3a1qw2W0DKpwqKrZ2eWLduHETFRVb02zd25YIPVWo7orJFGt6F9YMHDqhNKpskVJS1+zl0IlmItKA22nPq19v1apbUu6N2L1z8803d7knZs2aSXl5OV6vN27+c7lcbNq0WTUtzDFKQCkykm27ilMNHk42teN2DfQi09Vst2HDxpRFQSZtn73A6WtfoBgmkxGXMyaoDnaY/iJKUKWl1Rvg4LEm2vxRXM58moFTSb53t26tYOXKlfF7oqpqdzzpdvPmWOj4ihUrOHjwEHfffXdak7Cif6ggCUVasmlXIaXk3ZMthDHgdAxGrlD6gp5lZWUdeSkkmWW6a4nQdRHJTQRWXKOKRHnrRDMum4lxI92q7h/Q7g9R0+BBir5VGx9cZJKGpJvvEjUmIQwZavcp4ZQrlIBSdCGxYymkb1cR1SRHapoxWczYkhzaA91OQHQU3kxm48ZNVFdXU11d3dHWgF4LnFxGYJlMRlwmO5GoxlvHm3HaTIwb6Tora/4FgmFONHiJRCUOu22Y1NtMfhhKDCnvTLrtbat4RW9RQRKKtCS2NNDR21UEwxHeqW1j2/2/w+f1DWpjNik1Vq5cmVRKJhG9VE1nCRrINg9qIIlqGu3+IHazgfGj3FgHNIhkaBAIhqlt8BKIaLgcw0UwpZIpaCbZJKxrWImvVeWIzKggCUW/0DuWJrJ582baA2EO17TidFrxeX29COPORQtsmWRmSUeqsKyoqEg9sj7st/8YDR2BACYLR2pbebumGU97MC9zGWja/SEO1zTzzqk2TBYLbqd92AqnzD7MTpOwHiCh+6T01+kDdxS94cx/jFP0iXQdS29bvZavlH+LArcDIK0/KJ2TuHcJs93RuSisWVMeN7sksnLlynhxzqHYxdRgELicdqSU1Db6EQ1eRhTYGVGQ3yi2XNDqDVDX0k5UA4fdisU6nI+n+7bua9euYdWqW5BSi99rQBdNKtO95w9GaGrz0+wJ4AuEaQ+E8fnDBEJRoppGVIs13tQ0SVRKBGAxG7GYjVhNxvj/HTYTbrsFtyP2z+WwnFH+TiWgFF1I7Vi6efNmbrtjLTv+8ihmZ3GCoMnU4C3ZpJZ9K+2e0ReF1Cfba6+9ln37XslQgqbbo+0y38EQZkIIHPZYF8ZWf4SG1iYcVhNjSpzDyvwXjWo0tLTT7AtiMBiwD2gbl1xdq2zGyS5oRg+UgK4Pal/6yu28/k4jNQ0eahu81DZ4Odngo7HNjz8Y6cO8s8NhNVHktjGiwNbx8GOjpMBOsdtGSWHsvWK3bVgIMuWDUqQlMYqvxRPgZFM7FVt/naL1ZJsI25eE2Ux0jlVaWkppaSkQ69tz7bXXsmPHjviWMX9BpuZyudDscivcopqGPxDCJKDQZaGkwDFkE3/bfAEaW/0EQhoW60AnZudOC+/9ONld42g0ymXLP4G1eBKWgvHMmHsRtQ3eLtvpWEwGigvsFLutuOwWnHYzTpsZu9WE0WDAIAQGo4j9NcSiaEPhaMc/jVAk9n9fIIynPYS3PRT76w+RTUszg4Ait42RhXZGFjoYWWTv+L+dkiIHowrtFDitA6bVZ+uDGj6PaopBZf369UgpafYEqGvx43bZUwRKzyaQ1Gi43ITjpj7ZdnLgwIGk13p31MrKyjTaWv80u9yZLTsxGgy4OsL12/wRGtuaMRsFRS4bRUPgidfnD9LYGqA9GEEYDLE6i5bBMOPlSgvvyzix15rU8IfbaY/4CEQChKMhahtb2X+kjrdONHHo3SYKL70KgzQhpJljjaewGG1MGzuSSaPdjB/pZvxIF+NKnIwscuCymwfEL6dpEl8gTLMnQFNbgKY2f9Lfxo6/LV798wCHjjenHctkNFBSaGdUoZ2SuABLFmYO28AcR3wOAzayYtjT1OanvjWQkOPU17yh3Ibjdg0Hjz0ydrblJqm5XFlZWdq8p2x9aF3JrdkyHbHitLGfZ0t7hLrWZkwGgcNqpNhtw2GzDHjgQSSq0eYL0uYL4g9GEAYjdpsFp9PMwJhCM2kr/blWiWQ3jj/cTl37aZr8DbQGm2kNtuANe9Gk5GSDl6MnW3n3VBttvoQgFzOYRZTpk0fx9uuv0HT6OBfOnc2SKz5MgbUIl8WL21qIxVaC2WIl1nw89xgMIu6PmjymION24YhGs8dPQ4uf+lY/ja3tNLT4aWjt+NfSjtcf5nSTj9NNvozj2CwmRhXZKSm0U+Sy4nZacTssFHT8czusuJ2x+ditJqxmY6/uW2XiU6Slqc1PXYs/iwTcnkwgmTWtXGbdJ2s09KKeXu9q7+nfAXJotsyG2HmNahrBYBgtqmE2GbCYDHHzkMVs6rNJJhLVCATDeNtDBDrMSBEZK477u/v+X861xVSy00h7e60y3ZvJ4zz/0nOc9NZS6z1One8U3nCnaU5KyekmH2/XNHP4REuS78hmMWHSvJx6902KHQb+48br2Pn4E1Tt3s2YsWOwWqyMGTuW5cuXxfe787Gd2Ow2Pvj+qyixj2SEfSTnlczFaBh6+XGBUCQutBpb2zsEly7E2jtMvNFejSlE7Ly9ePdNysSn6BvNbX5Ot/jj5qbu6aly8+D0yOnUquiFtpZes+suh2Xr1q14vbGFdLCqCCQu3kaDAYfdEl+8b755Fc2+CPWtQaTUMAAGowGjQWAQYBACRLKuKTWJJkHrSLiORiVSCIxGgcVkxGS2YO/okCKlNuDaYnYaafrrmulaZRJ4LpcLr9eLNEiirhDRgjC3/OJGli3/QNJx+PxhDh5r5M13m5I0pUKXlenji5g6rpBRxQ4MAnbu9FK1ezf/8z8/AmDxkiUsW/YBHu8QVgDLly/jvvt+R21NDYsXL6Y94qPd4+N3D97HRY7LhmRxWZvFxMTRbiaOdqf9XEqJ1x+OC7BWXwiPL0hbhz/M44v9bWsP4mkPEQhGCEW0XgWI5EVABcMRNE0O+7DaM5EWT4BTze057cXT+woNfQ8+yN4v1qnZzZo1Myk0/eWX97J06eXcckvyorF169Z49ep0locNGzYmJAnniu4XbyGSzYEZR+mYrxACYUyfAJluUd+0aTMulytel7Dv5rXu6Mn0Rpdrpb/Wr9WqVTcnCCpJW5snbuJNPGfTSidzsPENFn7uQpYt/wA7dz7O7t0vIZEsX/YBjtd5ee2det495YmfM6fdzLkTijl3YhGjihwxNSCB5cuXxQWR/hpEx1+o2r076XOdnTsfZ+9TL3P+ZRcylFIhskWITnPitHGFWX0nGtUIhCLM+312+8iLgJIa1DV5GTsyvWRW5Ic2b4DaJh/unAmnZKHU9XVX+hd80BttTeB0OuPFPzdtipVzevnlvRw8eIgFCxamzLezNtvMmTOSWoDPnDmDBQsWUllZiRC5blCXG/9Lz3b/zIJQ9+ENrLbYfSCNy+WKX6sbb7yJ0tJSZs6cwcGDhzCZzOza9XSXh4rS0lK2b9/OA5UPEC0IcfHnFyNtGguC8+Ma0/Lly9AktEsXlU+8SYsnpi0ZhGDq+CLOmzKCiWMKMAiBzWTDYXJiN9kxGy2YDCaMwsj27ZUYW62ARBolzz/2Ih/86JX4wj6WL/tAknBavHgxVVVVVFXFIkw/vPRjrF1z9pRGMhoNOO3Z+9/yY+IT0OwNUVIUPStrkw1FvP4gNY25E059EzT9Dz7IVlurqNiKz+fjD3/YxubNm7NY/JMFRSL3339/TDMRYoBae+cyCrL7fUCyICwtLUXKrqbQG264kfvuuzeHc+g+kOaWW25h1apY64vKykqqq6uB5Pw3gFWrbmbTps1UVlbymbJreOVkFZHiIBgkl135XhIflAKhCAferqM2PIZgKAoEcdrNzJk2iotnnsu0kkmMsJdQaC2i0FqMyZC6XMbm/PT257h+xX/E+0O9dP8rTNdmc+edd/K5G8qwnHAhLRqaNYpNOpJuz9v+444cnsMzj7zFrdodVmrq2vK1e0UC/kCIY6e9OTTr9babrU52XW0BfCEvp3wnOdx8iFfrXmbvqZfYXfs8L9Ts4qXa53j51G4O1L/CoaY3OOU9iS/kTTDLdc5Pb5uQSOx1ukVDdNkWYNOmWEkovbpAuvPR/eue6H2fqr7RtX/WnDlz2LFjB5WVlZSVlcWjJKurqzvC+PtWsirpVUfidWJbi8RWJ1LG2pYIYcDtTra67Nixg4MHD1FSUhJP0n5gxwNcsWIpB40HiJQEwBDb386djwMSfzDCi6/V8PtH9rP34GmCoSjjR4zg+qVXsAQn59YX8pnSz7Bo3MVML5pBiX1UgnBKnLuImz/XrEnuD+V0OvnCF67j8BvvcP6kObyy81W+8P6b2PPAfmxvF2E95sZ82sH23+7o4zk8O8hbkITRYMAXkfj8oV6pfIrcEgxHeOdUG26XI4ej9scslV5baPQ3cNJbQ0N7HQ3+eoLR3texMwkjJY7RjLKPZsWqz6IRTZqfzqZNm9POUy9Um8isWTNTfFzJ5CIZuHdtQ/qTPNxVEAohmDNnDq+99lpS2D5AQYG7F2PHxq+oqEiKtoTY+T506GCXunZA/HUsIORmdu16Ou3Iy5Yto/KPlUSKg0RGBKgz1FJVVcXiJUtYvnwZO3c+zp6XX6Uh6MCruYhGNYS0YGiMsGTcBH765a+xceMmHtn+z4zaerpr6fV6k1pvQOye1zU63WeWaGKdO2cu9913b/w6mrGoCugZyGsUn9Nho6bew8zJJfmcxllLOBLlcE1rlppTNgtfsmltzZryPpilOhfJqD1C1B3i5s3X874rL8/iu90TkVFO+05y2ncSgEkfGk3waS9Grxmj18yel/Yk1VVLDapILFSbGFSh9wdKdz76HwGXvV+tf8IwsyAsKyvjtddei28ZCwQhi7l3UlGxlbY2DxDLUdO12cQWKeXl5XFzKyTXtSsrK2Pjxk1xH+ChQ28lzFzyyAsPE5zWhjTFtK2amppYqLfVSiAUwTX2fIpnltAaFpg1NxdPvpAbrljKP/7v9zFt/f/+D+juISq7a5n6cBUTTrHWHAUF7oRq+2S8jopO8h5mLkwmGlp8jCxy5nsqZyxSyqQnOCljRSgP17TgzKINQjYLX+o26bSN5HDvdAIP7tp4Fw/8cxsXf34Rl115KTt3Ps7zu58jTIjlVy4nUWD070ct+ddjj6K5wmiuMGEJd/7sNv5z1TeQyLRBFYmFavUnZiklbrc7wTmfzv/VvwCH7Pxq/RWG6QWhlDLu79HpvI490Zl35PF4k8yEiUEmiYnU3Z2vioqKuPamCyjNFsF2jpG6cC0QC0IAYkEIwoh77GweeOx1QmGJNTKOy85ZzHUfWML08UUAzM7at5fNteyqgepRkHpqQuI5dLtdSnPqgbwLKJvFTH2rn2K3fcjWHRvOpOuMe+ed5Zhco7ntjjuyCPXPLkclcZs1a8rjPZtStQ0g5QcbE2Z3bboLv9NDnfMU8z8yh8uuvBQ9yqqmpoaa2tqkOe3c+Tg2m5XLLlvah7MS+37V7t0xE9CyD7Dz8Sd4fPejtIpmPvGxTzC7ZA7haBiz0Rz/Vqqg0Ft53HLLqvi4XTWWngIcsjXJ9Zxv1l9h2FUQxtA1nJ7Ni52kPrCsWROLkEwUTDrJofmZz5deKPjGG29CGjXCIwNEC4OUjBqPLWrBaDJ3CCaBpWgSrrGzaQpZsEXGcvmkxdx41ULOmVAcH1evRJ5I9xVOuruWmTVQPfpQ3z7596Pojn4JKCHEZ4D1wHnAEilln8pD2KwWaho8TB6TXSy9Ijsydca9+w//xxdvvhWTsTsTXefrbBa+dNt02t877fO6cNq+fTtSShwuO/965e8cbn2LxRctivsLdu16hqVLY8JnwvjxVFVVsXPn4/HPq16qYumSy5lcMBW3pQCr0YqxI+xXkxphLUQoGsIX9uIJteEJthLSwvEjs9mscf8EEP9rs1nxhX3sPbWbA3WvcG7xbM4bORebKbXcU6dWoIeWp9dYMkenJftjclGhIRfRfslaY0GBOwvzYvJ9ky65VzePpiNVs84czRfTSPYd28vCsgu44sr3xx80Fi1axJVXLuenP/8tjjHnYbS6sUXGcf6Iedz4wcXMnzE6aczEnDZdK165cmUPwrf7SMNMplin08mCBQsHMI/szKW/GlQ18CmgX525TCYjXl8IfyCE3aYCJnJFYtPBLVu2sGXLFqwFo7n51q/w1a8m/zh6MuP1vPB1Z3/v/Fwfe/+B/Wx75HeERwXAqGF3OZBSct/vfk9th/8AYOnSpXxg2TKi0siru9/hjef/goURLJv7n1w2/zJOvy04LcBqNlLsslHstlFcYGVkoQOLuTOFQUpJW6iVX933c+rb61j2oeUEowGAuDbWWZImRkgL83rjAQ41vc6skjnMLilNElTZmnzSPVXrPhjdH7Nu3docVGjIfQvynsyLvUnuTfUd6ea6dD6ndBrbF2+7iVPOd5l/dSlXLH8/uoYNIEw2/vC33bgnL8YSLcHaNJJ5I0388EtXYTQaSXxgAJIi7vT96q8z+ROz6Q+Veq70B7JUf2z6fShS6ZeAklK+AeSkaKXTYeN4nQqYyDW6kNqyZQsWxwgKikfyta+mVjvo3oyXnSkkvf09eZvY3+ZAI455JsKh9vi2/vZ29PqMY8aNp8mnUdtmpPLhQ0j/REyRjzNtaqe57bgfHtj5RsbjNhgEE0a6mDqukKnjCpk1aQQzJhVj9Tp5ZvuLjI9O4fovXcddW39E1UtVLL54UcaxIjLKaw37OwRVKeeXzO0w/fWksXR9qtYTTmN+q5hpMLGwbX/akGzYsDFeuX3t2jVs2LAxK3Ncz2QyL3Z/36QuyocOvdWlRYruk9Ij4RL9fHE/GBo+VyuPHP4r8y+bT6IQCEc0CifM4eWDpzBqbopD0zm/wMaxQ//i8V0HOf7arniumy4sOn0/MeGY2D8s5l9Mb97LLlgl+XelWwtefnlv0mi7dj2dUv1CkY5B80EJIVYBqwDGT5iY7nMwGqlv9jGqWAVM5Aq9+aDR6qJg1DgEWtqSPJm0gZ6eavXvZRMKrUmN6vpXea1hH36RXCHZYLZjcY/B7B5D2FHCKPcERNsETNIV2yAaIOitJxr0Eg16KS6wc/XVV8ePIRCK0OwJxP61Bahvaed4nYfjdR6eefUEADaLkfOmXsyMi3z8+eGd7KiM9Y66oWwVHq2VV585wAWXzct4LsNahOr6fRxuPsgFoxcwvWhGN4IbdL9JcrCAJ6FSBV1KJvVVkFRUVFBdXZ0QcBCjtLQ0h1FiXbW6nu6bRGbNmsm9994TFwB6YrOUMY0GZFKXWrfbxWe/8FnmfGIGTYEmtI58qNhUJG/XtPBCdQ3tfokrNItJhgJ+/NXPcO9vf8PjBw9SVFTEwYOHuO66Tn+obtLT59714SJz/7BVq3RfY+rvJtO57fS9JZo3E6uXKDNf9/QooIQQjwNj03z0LSnlX7PdkZSyAqgAKJ03P21mmt1qiQVMFNjz3vvmTEAXTr/41VZW3fE1vvmNr3HDDTfGn9b1cNcNGzZSUODOqA1k8+TY0zYtgSaer9lFc6Ax5jeoqmLRosW8evAY1hFTsDhHAgJbZAL29olceO5kzpsyknMmFLL5R9/m0Ov74/OaNWsmB186xOnp1ow/8GAowrE6D7/f8TCNXom5cCLH6zy8cqgOCucy/tK5hNpO0376Dfa/8S5vvPoSK1as4KPnfpI3G1/nnZa3icr0lZofe/JRHvb/FZt08NJfXmbenAuIRMKYTOYk811ndF9sflu3xoIqysrK0uZfQV9NcjFNprq6Oq6d6BUXysrKEhbWvtOdCTj1vunuoSZxUU59uIFOf9UD2x9g+efez9fvXUMgGEiqCP6Pfz1JY9hNe8SENTqaBSWLuPVjS5g9ZUR8DOg0HerCQRdO+ud6vzCdzP3D0uVvZecvFMLA0qWXJwkoPWhIhZf3TI8CSkq5bDAmomO3WzlxupWp44t73ljRLUIInAVFfPH2mHCC2BO1nnsCsZI11dXVXHvttV1+sPpimaoFpPqTumoKJG3zVtOb7D31UseCL7BYrJw7/3KOeM24J8U0CXNoBMZaCJ56nZaGh3Efm8L6Dh/W5e+9KElAZfMDt1pMzJhYRBEN/PNvsTyaMUVG/vXMq9hLpmMfNQNLwRgsBWPwAed/6FLmLb0YC24uGv9eLhi9gNca9vNW05spgkoSCASpqqpi/IQJzPvkLOrfquHUK02cd85srr322vj5nTVrZtyMI6UWd8rrvheddJGOvRNSyf6w/psLU+m+Tl/qfZPog8rGHJaqhWmWKBd9fgFzl8+Ot68AuOzy9/Hgo1W0RoowSQcTDKX8x4cv44qFUzEaRDzXat26tV2Epk6iUNQFOCSbWXW/lH5/b9iwMek3011ATNcgI107TD4/3VXMV3SSk35QQoingHXZRvGVzpsv//nv5zN+3u4PMn6EA7fT2u+5nc2EwlHeOtGM22VLqvSs+yoS0c0OqX8TF5neJoKGoiFeqn2WY23vAIJIVOO1I/W8cqiOQKij5L5fYquJ0PJmNTISC1ooKiqipaWli4lRJzEXKTPJUXRdFithxFYyDefY87CPnoXBFLvXbBYTl86bwIcuOYdp4wppD/uorn+Vwy2HEkxMnWHqOmPHjKWxuhVjsxWBoKSkhMbGxi4aRGqgAMTqyX31q+uAzOal7OhLb6vsx049j4mLe7peX8nXSF9nMofXL1y0iMiIIJESP9/89jdjn0uNnTufYP/bp7CPnonBaMUenshnFn6AUM1eAu2eJI2osrKSOXPmMHfu3LQh7Xqi7NatFQkFZlclnTegy/EkHmvqNpl+G4mCLZ25/GwWUtm2fO/X2RFCfFIIcQK4BPiHEOLR/oyn47BbOdHg7WKfV2Qm9VxFIlEO1zbjctpTfgQioRJAjLKysi5mkG3btrFixYqEJ97k+nqJ9dM8Hm+8XppOa7CFfx15mG1//z07H3ucN442sH3n67xQXUsgFMEQbed8s5OPFk2muXpPXDiVlpby2GOPJZW90X/Qe/ZUxc1jMf+GJF0ds4qKrQm16tLXz0NGCTS8TWP136h56mfYml/mvCklBEIRHt/zLmt/8W/W3/Msbx7xsGjsJXzs3E8zuWBq/Bzq0WM6EydNJDzKT3CKB80WYdmyZfG56jUFy8rK4r4nnWuvvZYdO3bEA0w6NdbeoV+PRDrPQS5+R13r9K1btzZtGLp+3+j3Xef1iM+WDRs2UlGxNf76hxv+m9BkL5GRfhAddfOkxt8efZqaQDHOcaWYDUXYT47i7i+t4uarFxBo9yTUe+wksSwTxKIHdSorK9mwYWPcF+jxeLpogHr5qsRakOvWre3yu0kURulqTyaWhkos4ZR4b3eeA0U6+hvF9xDwUI7mkoTFYqamvo2Jo1VuVE+kJuNqmsaX136bUaNHcestt6Zs3TXaLpVks0zPVRH0J0P9yb/Wc4LnTjzJ4089waEjxwlYxnPYdxwALeTBf/wwnzh/Gd9Z9x0A/m/HH+P71itkJz6ZJlYSgNgPPiY402kcskvS8FVXXZXxWBOfjBeeW8TP77yZR186yhN7j7L/7Xr2v13PxFFuPnbpOSydfzkzR5zH3lMv8se//TFpHD0CUVqjBCd7OGk8xo/u+GnSQrl27RpuvPGmpO8JIRKOp2/aTmJOjx4l9/LLe5P8YQUF7gyJuKmaDWm2ib3XXX5SpjD01JyoNWvKufHGm+IBHZoW5Rub/5O/7/0LzpFOfF4vY8aOYe/L+3mnUaBZSwCBxVtMsPotvKcO8af7Haxduybt/ZipSsXGjZvi2kxqTcF0GmAiiWbB9MfffcqBlDLJCpCqSfW/KsqZy5DVLy1mE23tEfyBUL6nMqRJTMYtLy9HSsmX7vwGf334b/i8PpKfnpNzOXSNJNUUkqh9JNNVG0n8sXk8Xt5oOMBTx3bSHgpQ22YkXHAeJnsR0XAAb80rePbuRlS34xCxSM105ZAS963nYAHxxNjY+6syVEhProi+ePESGhubgNiTcSKJfZ10ITFxdAH/8bF53P2fV/GFq+ZQUmjnRL2H3zy0j1vu+hfPVrWy/6Ej7PtbNZcsfg/f/OY34jlbOosXL+bRPY/wkXXLibg779+VK1fGF+Y9e6pYsWJFfP8333xzyrnOVutJ7lNVWlpKZWUlhw69xcyZMzhw4ACVlZV4PF62bt2aUoFccsMNN3LDDTcmvZes3ZCkLSfOffv27V20D/0aQExzStQa9Ouh+3KC+Fm3dTVPHnwMhOS82bNZuGgJrWEnhecuRbOWYNRsGN4M8cGRhTz79/uTqpxD+gK9XWbTkUgde/jpJJ0GuGZNeZf7ZMOGjfGQff13kzyP5IR2HX1MPVE98RwkBpAo4ZSZvJc66g6X08axOg8zJ43ISa7VmUhqMu7/3nM/BSNG87kVZWlu/tSIvE5KS0uTKixD+krZqaHDeg5J2YpreU/ZIl4+XUVdk4/H97xLW9iBEJJA01H8p97CdNKCxeOIC8VEp3Ril9TOfdOx/67O9O4DAbrmJ6UK4jlz5nDfffcmaH6rksZw2i3UvfkkM41efCNLONxso80H9z/2BkKbyMxZn6XUMoIZI87DaEjuaSalZMzYMZw8Vcvijy/huk98gZf+/DIPbn+IWbNmdnnqPnToYErIcXe+ve5DvRP9W/r/y8rKWLOmPKmX0n333Rv3kUBnN2A9aCamXcSEll5tXI8Q1Pep1+mrqKhI0mBTzV5AUk6URLJsxft4YM/vQMT2MXrMaGbMXcITuw/jGDMydhVPa3xs1jlwoY+1azvzooAkDToRPfAhMfG5urqaOXPmdFlDOntaEZ+znrCbqlGVlpamrRmYaALvPkl6MHp6nXkMaQElhMBkNnOy0ct41X03I7qQ+nXF73CPGINAy3jzp5p5kissd9eBNln7WrOmPC6cpJBc+OlSDje/zb63TlP1+ik0KRlRaOPy+ZO45+ePYalxYgilb065bds27r77biDmk0lcfFwuV0Ih1mx/5D2bMX/3u/sSjrezwoD+fSljbcN3VFbGA0ZKJpdiHHMhtuLJRIpK2XU8xIs/eoamA2Emz5/KmCmjgFih0vHjx8eraZ/ynWTah8ZzBUuZ5pyR5BNMDQLJHCHWu1DvVBIX5+rq6nhQQLootkTSBQl0dg2OjTV37ly6hmJ35gAlmr00S5Tw2HZOG2riwsloddFum8qjLx0FjEh/hOj+GsIna2Hq9JT+XF3D03VhogtX3aypC9EDBw6wY8eOLser97RKrM7eXbpE8kOM6HIdsskVTD23Skh1z5AWUBAz9bX4AhS7VBmkTEgpue2OtRSOHo+BWLBCT0UvdToFVudn6c0mydW8dU1KGjWCE73c++AfiDinUFMfC6mdd+4oLjp/HNu3VmJ914XQYotyart0fa4QWzRStSnd3HLLLbf0qqJF4mKRLmox8XuJC39FRQVtbR4OHDgAJGtejceq4Vg11uLJjCu9Es0+hpB7FtPe81XG2Jp5fudWLlw+k0WLFlF78iSnT51m8uQpIDWefuZZAoYAl37gYk56axjnGp9S0bqnIq+9C/VOJFH4pGqS+uKc6rdJX90iJujSz5Nu6+9JZDxCb8zY0VRVVSGMZuyjZmAtnhLznUbDRI/U88nSJfzn337djTafPvfuvvvu5frrb0hoDdJZI7G+vo66uvoee1pl25E5cR49CTbILpFd0ZUhL6AAXA4rx057mDlZmfpSkVJy++q17Pjrv1j5+b7d/F2LlqYPd07M8t++fTufWvFJTjmOs/+dw7zT5sLQ7sVmNVEsGok0e9j/eCN1L7Uwe9bsuNDRF4hUk+K1114b11R0zUyvBLBgwUI0TYtXBOj+R947M+aaNeXxhTWxPp5OONzVBxpsPsbRZ37LR665AdPYC3nlUB0nAyWcM//rHN37PE3+fxO2tnVWSt/5ODW1tdTW1IAQeJZ72P3YXp7944t87rOfj8+/e80wvSM+USNIDYUuKytDShnXHtLRnWBLN5fu5plufrNmzeSNo68TGtuOtEYZM2Y0p+vqsY6Ygn3UDAxGCyAJ1Z8ksv8EwhPBOPeSpPEyBZCkC/rQw8v1BxC9/XsmwdzTA1zmbXqaS7Jgy7anlyKZnORB9Zae8qDSEY5EsBikiupLwR8I8Z3/+RmRcDALH0b6CK5MT3fpnuJ1jaPeU8esD0/ngb/vpC7gRAgDTnOEUeZWXt7zEp9Y9BkmO6fFC2XqTduuv/4GhBDxaL3EsjY333xzXDjppGsOmNgELhtfTWICZ/rzkyFXqgeqqnYjhIG3jjfxf08eZM+bp2J7jkZort/Fjbct4smn/xnr7JrYp6iDi5dcwteu+wYzR5yXNtcrk/kyNdcp8QFD1wAhVsnC4/F00R51k26iZpnuvdS5pEajpZ9n5/ykQRIe6WfBlRd0COon2PfmMRxjZmG0xIJkbMYwTVUHmGWfxOVL3x+fb/IDRu+qanSXr5V53oNBf7odn1lkmwc1LDQoALPJhKc9gKc9iNuhEnghloj7zqk27rj9K8QUy8xmiZ78F9C9iSnx+9dc92keO/II9/71SdoiboQAf9M7NJ1+k+MaXLPoc/zXmu/H99M5RgVz585NmoNedUCPykslUaNKl9yYrCWlz+jv2WzT1beVDj2BWGflypVs27aNGZNG8M0vXMx1X1xNvTYW+5jZjBh7BQ/9MURjawHjJk6L5U0JkSSgrlj+fvacepGf/X4TL/5lDytXfKEHDbjnUG/9WHVhkhg0oPtodNatWxt/nZrjM3r0KK644gNA7L7QNI19+17pQYPtNNdGXSHCo/2MnjiK5cs+QE2DD49lMu5JsQCIaMhHoOYwE8U4PvieT7J27ZouNfr6tninv5aZQskHV0j1XjM72xk2AgrA5bBx/LSHWZPNZ31zw0hU4+14Im5PZoqemg72ZGLq/H47XqYtn8Lv//EC/qgdkLxvwWQe3PYIaALrCRf/tf37pHvyz9Q/SfehVCYEJKQjURjpAjOxnUFi99KuWlWq4Ep83XNQBUBLSwuzZs3kssuW8swzu+IFP/UIuTf2xawCH7tmJc8f9OAYez6jij9EBC//9/A+HNZTSePdd9/vuOGG65H2KBd8+nwWf+oC2oJtGcw/2Tni9WPTF/lUH41ed1EfNzmKLeaTufbaaxFCJFUbr66uTmpPoe8z1lHYFReI9/95G+/5/BIuuXIJO3c+zsv73+R3f6vCH435j7VIkNHOEHNKxnHkcDuHDr7NwtmxSvLJD0yZK050v7B3vZbJRXSVeW04MawEFIDdYePdU61MT+iMebYRjWocPtGMw95zu/YYmf0X2UUYxb7vw8OOV/6Iu3kRJlsBRqFx9WWz2P/yC4iIAcsJF4aQMSkyKrWnVOoc9AV069aKpPJK5eXlXHbZZQSDwficdH+J2+2OtzHQv6NXjE6uWN1zBFxiA8V0piC9NFFJyQgaG5uYP/9Cbr31Fm65ZVXcNCmEIR4NGfP5bAOg9fAuCqa9F9f4ueC9lGZ/LVPmFlFgaeXA3heoralh52M7O1rZw3HPMWq8J5hRPIvb7vwKVqPeeyp2HKn17XRB0bOPppNkTanrYp/4HT0gQkfXYO++++74wwHEcrh+sukn7DnxIkUX2QnbgtS3tBN1T6dw+kj8UTCbDBRb/JQQ4b9v+TF/+N/7OXRwZ7z1SHLiKt1es8yloDIL8cQQ+e4DHxRDiWEnoExGA4GIoK7Zy+hiV76nM+homuRwTTMWmxWjoTdaZHrTR6JfYdasmSxdenl8wYZOTeqU7xSjLx1HQdvFGM0OoiEfn/vIIl564Rn2PLeHz110A5Yl1i7hyqmVoTNpanoo+YIFC+MaiS6c9OoIif6SRGEHmSpW9xwBp5edSaxWUVZWxoEDBxBCcO+998Q1M4/HE9c+OkOtOxOKU8sNRdqbuXKui//7y68Zfd77sZXMQnpG0mg5wvi5LgpN7VhthqR5Pr3rKR4NPMpHPvhhZpbMYXbJ+fzyZ7/i0KGDzJw5K6mIaad5NHNpJN0/lb2mSXw+6a6VlJ3nctasmbz51pvsOvQEhz1vMWbcaBpag5hbrLz+ZOx6GI2COdNGsnj2BC6auIRZJedjEMak66dfu9Rmj91r/ekTybMPRlDCaTgwbIIkUvH4/Ewb4z6rQs+ljAkng8mM2ZQ+p6ibb6cNBIhVUHB3aX+dWLrolLeWHfse4sGnX0eTBiL+FjzH9rBo4XwKrUWMa5/M19d+HaCLkz110UmdQ3eBGLqzP3E8PRJPHy+1yCdkKpLa074TF7zksPtMvq2uyC5+Hl0rc7vd3PP7StxTLsI14QKi1iBeyyFsrgBzp4/ivCkjePqpJ+PRfno7+id2/pu9/9rH+SNLeevNw1kEsqQ/5myKAffmWFq8zTyw8w9EioJgBEvBWGwjpmGyFwExjWnOtJFccO5oZoyazqKxl+C0pD5Q9lTctrtrBpmvhQpGGOpkGyQxbAUUgMfrZ9ak4rPCHyWl5GhtC9JoxGzqreKbXW5QugW7zneKbXv/zF92vYFE4DAGWfHhxTz15L/Z88xernvPTXxtzdeSqlanX3Qym1+6WySBHsdLpTcRcLlbuGSSTy2xkkGiTwdAGC1Mv3AZPttktIIQPvMRpAgQaD7BtHEubIZIx9xiLF6yhKuu/CC7//UyT/3pmXjCc3aRaH07T6kRfhB7UNDMUd7/maUsvupCfvjjn2IpmoC1eApGc8wcabUYmTNtJPPOGcUo9wguHLOESQVTsppX+oeVdNe/v1XfFflmUKqZ5xu73co7J1vyPY1B4d1TrWiGvggnSGf6WLdubZe+RKmRbQ3tdfy+6s88/NwhJAKXKcDKjy7BbBRc+7Frue49N1HsLu7wg8WqZmeuqJ3e/JJcMb0r6cZLrA2nJ/Im/k2ukabT3dxyQefxJVblKCsri1fXjtc//OynObznEQrq/03b7l04TozBHp6OrXgqtYFimgwTsZVMx2C2A7B8+TKiUmPhBy8gOLWN4EQPkYIQd9x5Oz0L2EyJ191FsIm4T62srIztO7Zz/99/T3CCl5JL3RSWTuY3O56kcOb7cYyehdFsIxr04K2tRtS/wqWlU3jP5PfwkXM+1aNwylTbrqIiVjfwpz/dkPTN66+/IV4XL7kGo+JMZFhrUADBUBibiTM6P+rYqVZCGlgt5n6OlGzG6s4c1xxo5LcvVPK352Ka0+wpI1h64SSefeZZjr7+Lh+fcw1fvvXLQKymmZQyniSpC77EXjipUWad80l8nTzX7jQu3YeSnW+lL9pbX+lqXuqaCN0ZoKGXeDI5Srj56z/kH9VPUx96N35uIv4WRheauXLpQqpefDapB9VFSy7mP675IhPcUxjvmojFmM7cnY0GBanHH9EinPbVcrztKF9cdzOmgpFYCsbgHjmJUDhWrURKjbDnNHa81J94CySYWq2suOI6vr7u613GTKT7AIhVSffmtddeGw9x10ksz6UYfpwVJj4dnz/I6EIbIwrsORtzqHDsdCuhaC6EUyLdL9hfvO0mtj73AP948SCaJgk0v8uciU6WL/8Av9+6jYaXPMyeMatLdQi9URwkVzaIBRbQy4isbKK4Un1DmX1F2UWEDaTvInmsxIocOrEoQHjoqRcZf8kcmkythCOd3XyjoXaKHQYuWjCbt6r38vKel+K+KoGBkY7RjLKPpsQxihLbSBxmR1IAzMGDh+LRiHp0on6N7G4rn/jcJ2gKNHDSc5I3Tr7Difo29r9xFF/YgDB0au42YxSn0Y+n/iinao6BBGOrlbljLuTCORcmtLrv6Tx2d767+sASya15NvdIKQmFo0S1KFKTaFKClIjEe6DjrxAghAGj0YDZbOxl8NPw5KwSUBALmpg82oXLfuYk8eZOc+rK1q1b8Xo7F2x9wbS4TIRmWnn4hdfRNMm8c0fhO/lGrGpC0Ij1hIvZ587u8jQLmTL1Y/Rdg8ml0Mg8Vu9DmvtD9+HQpaWlrFu3luZ2D3979Wl2PPUP2iICRLJ51ySiOMwas8+dzIgCG26HFbfDgsUc81EZhIGXdr1ExK/htDkJ+kMIBE/8+3GmTpvKpKmTOFxzjJPNDZxzfiljJ06jsS1AQ4ufUDiatK/RIxxEvPUceb2KhfNLWb7sAxgNJn607i5MTTZExMCePVVs3VqRdF/17zymD4KBfFWCSE9U0wiFIkQiUUwGgclkwGwUOGxmbBYTJpMBk8GQ1lcupSQS1YhENIKhCO2hCMFQlHBUIxKVmMxGbBbzGVfi7YyrJNETbqedY6e9TB9vwDYAC/pgIqXk3VOtRKQYEOFUUbE1bhZLDFm2OM1o59n423PVaJpk7rmjeE/pBJg7gT3P78V6woWIGti2bVtSSaJ0BUdTF49sqlWkp7fZ990JtExj9SWkuT9kVzW72FHAFy75GGVLruTtprd47sge/vXc83iCBsJYiUShLWRk9+snk0a3WozYLSbMZiNmUxFmt4EWIGq0EdU0Rr/nchqCEZrazYgRMygYMYP6INQfboiPUeC0Yox4MRPgg++/GKfdAsxkp8VPoa2QBWMv4uHf/QNznSP+HT1PTb8P+nceu/oMU8te6fsYbCElpSQQChONRDEbDdgsBkqKbTjtVgyG3s1FCIHZZMRsMmK3mSlK+Cwa1fAGQrR4gvjDETRNYLdbzgoNS+eM0aAgduP4fAFmTCrGNEwj+6SUvFPbgjQYsZgH4vkh/dP7A3+8n/M/s5DjYTvRqGTOtJFcdsFEEPD0Y7vYu+MAIhI7p6nVHrKvdTaQUXT91YKyiSrLNb2vkvD9jf/Fn/79RxYtW8CCiy/hyWerePdkIyWjxmN1FNDWHiIa1bKegRYNMXJEASPcNkoK7ZQU2CkptHcIpE4fodlgYqJ7ClMLpzPONSEpsTa1Kgj0t+6d7OJv7C4Nom/07txrmiQQDCO1KHaLkaICG+4+CKT+EAiGaWhtxxuIIAwGHLbhay066zQoiD2NOBw23j7RxIyJI4Zd+LmmSY7UNmMwmbD0KVovG7pWlZBGjdLPLOJ4xE40qnHe1BIu7RBOTz32NHv/eIDPf2Yla9aUs3LlyngeTaoPKjW8GugiLBLJbT+c/mpBA9VQri8aXeY5fnft93Di5oEHHuCVvxwg6gxz0fLFvG/5YhAGkBr+YJRAKEI4ohGOxMxFIDAaBEYD/Otf/6LuVA1aOABSY8bixZyoruGEENxww/Xo1+upx55mlH0Mt6z4MhPckzAl+KFSNUCXyxWvCnHLLau61aZ7pquGuWrVzUmVO/pzbXrzIOMPhIhGozisJiaOtOPMowvBZjXHg8FavQHqWtqJaOCwDa6gHEzOKAEFYDAIbHYbb59o5tyJwydHKtJRvshisw6C9te5GEuDJDJdUKM5CEeizJxczNL5kxBCYDPZmOtYwOxrLoj/mJcuvRyApUsvRwhDmoKjmbuNDmw/nPTlnLJ/es+1AE2N3uscs39P/jFhIKTA6DNj9JnZuvpefrz5RwinpKG9DumIcsnlFxPSwl3mtPOxnZw6Fuu0m1RhXYIIGdnz6Ct86YYv88d7/sQrla+zYsUFTCmcSuo5SA1S0dvOX3jhgi4tPDpLX2UfKJFa4LerUOq7H7KnB5lIRMMfDGE1GRhTaKfAZR1yPqBCl41Cl41gKEJtgwdfWPbJxDjUOaNMfIlEohqhQHBYCKlAMMyR2lacTnsOb7DuI6Q2btzEA5UPEJkewTl7EUazgyljC/jgRdMwGAxYjBaWT/0wRbYRXcaSUkuqKJ4+XDz5O4MXhNAXM2Juw9BTq2BkDrvv/bWWUotrsTp6lfXU6hDBaAhvqI1gJEhYCxPWQuz4YyXtXj8CwT/+9g9E2ICIGCgpKuHK5R/spmpHjzPrkraQaO5LlwQ+uIEpyXNNZ8798m23o0WiuB1mRhc7+1CtJX8EwxFq6z0EIxKHfegJ1FTOShNfIiajAWG38laHJjVUfVItngAnG324XJmqkveennJMdOG0pGwRtZYCmloDRPzNGDwNGAzTsRjMXDHlyg7hBImLCpAinKDrApb60CM7HP+J2w5Ewc6+akG9qeHW8xz0qu16omvqot0frUzvUpvoB9RbgKSWLrIarVjto5JG+NYt30O/Pv/avjP+flNjc5e9JRdY7YlY8ndqU0C9CHCssGxyDlhqzceBDUxJnmuiOVfDyJe+9GVGuq0Uu7MtwDy0sJpNTBtfTLs/xIl6D8JoxGYd/mXghuaqnSOMBgN2u41Dx5sJhFLNHfnnZIOH0y3+nAqnRBOGnpWv//BjmffgdDm4eMUimpwjaGoNUOiycm6JhsNuxSRMHN91ivt+/Qc6BU1sjIqKrT3uvaJiKzfccGOHmSdWXWLDho3ccMONVFRUpGyde+HUXXWC7uistp4sQHv/JN9ZISO1cC7Qz+TSTkG6bdu2tFtkK/xSBfnMmTO6zLV3dH042LBhIx6Ph4MHD+H1ekm8RnoUqX6NFi1aPECJ013nGbsnN6FhwOf10ljzFr/Z9INhK5wScdgtzJxcwgiXBa/PT6QXATNDkTNWg9IxGAQup40jtW1MHOmkwGXr+UsDTCSqcfRkCxiMOHLudO3eF6NJyeyrZvDiv2uprWnBbjPx0feeg9txHkZh4vLJy7jfW9nHJ1tJW5snbsrSSWzzPnBPxrnQgnobtJB5nHSV46H/gSF61fRNmzan/bzn8bOry6hrfqlV27MdMzGqr6ysLKNfcGACU9JTUbGV1lYP4YjG7++7hxs+90lMBnjppZfYsmULAJs3bx72QgqgpNBBkcvG8dOt+ELkNbijP5yxPqh0eH1+RrisjCnJX5uOVm+AmgYvDodtgPMZuvpipITnTjzFA88+w4HD9VhMBq5eOpORhXYMwsDSSR9ggnsS/Qu57uqLgMEsTTMQ1SB6Hw6evT+m93PJpkJ5aufh7pKSpdS46qqraGxs6jLfggJ3VlpkT2blbAv+DpQGJaXGj3+ygd/fdw9zzp3ExRfFmiT+/Oc/54477gCguLiY9evX53S/Q4EWb4DaQVlzsuesqySRLf5gCIPUmDquaFD9UtGoxol6D/6wNoBPM8kdYpPK6Kwo49KyJfy56jlerK7FYBB85D3nMGGUG4HgvRPfx5TCaUlj9D1nqWsFgKFemiYT/Wmap1eFgK6ln/oTBNBTN2G9d1R2ZZ0yV2Hv/UNFOkFOWiGkz3kw6iO2+0MYhMa4EQ6+/c2vx7UlgNWrV7N5c0wb7Y3mJKVM2j719VAjEtU4dqqVqBTYhkCLorM+SCITdquFqKZx6HgzY4rtlBQ6ev5SP2lsbaeuxY/NZsVpH5gqF/qilfjDT2xA+Lsn/h8vaNXUBWN18a5YOJkJo2L/v2TCZUwpnNZljETSm47SL0ipYcb6e8OvuGcumuZ1vNvRgl2vEtHdPnvS1rrWISRJG0oN4U8/505fW+4CRLpu2116QWqH4P4EpqQ7Z/5gCBmNMm5Ep2l/8+bNSQKqLya99evX09LSEv+ulJLy8nKKioqGrAZmMhqYPqGY+mYf9W1+XI7h4W/rl4ASQvwU+BgQAg4DN0opW3IwrwHFaDDgdtlp8oRoavMzYaQbhz33TxWt3gCnGn0IkxGXcyAL2SYvpHrS5MGDh1iwYCEtzgZGlE6mLhgzbb537gTeeWMvJ9+xsvazX2da0blJY6S2T9dbWEDnQphOs9iwYSMHDhzgtddeA5LNWrq5a3gJqb7lVqXm8EB2vpVM5zRZ40oWMKm5VvqiU1pamvWcu843VxGWPZV06v9+052zu366CYvZzDe/ekfSA6guSBIpLy/vlZCSUtLS0pLksyovL2fLli2sXr16yGtSo4qduOxmjp7yYLMPRs5l/+ivBrUT+IaUMiKE+AnwDeA/+z+twcFmsyCl5Fi9F5NRMKbYgdvRP/OblJKmNj+NrX6kMOAYUMEE+g88cSHVWbFiBVd8/jJ+cv/PibimIRDMnzGaU++8StXu3Vy96NPMGDGrY+uuY6QrLaMvAuk0i8rKSkpLS+NVzfUafRDLAdLbpffraKUkGtWIahIpZedyLToqJRgNOV4g+urI723ARddzqlfzjgn6TnNZ4uKe6TqUlZUlBaqkn3NqVYvuqlz0jZ6FUH8CU5KP/87yO/nxj3/KvRW/4tb/uC6pu4EunHRBkihYIHtNSggRNwlu2bIl/n19zKEsnHTsNgszJxXzzskWolHjgNT7zBU580EJIT4JXCOl/HxP2+bTB5UJKSX+QBhBFJfdQkmBHaslO/mtaRKPP0hLW4D2YBSzxTSAFz3Z2d3W5olrJVJqSUVc7/3n3Txz5AUe2vUWgWCEYGsNvppXAfjYok/y4zU/Jd2ileg/qqra3eFs79rSovt23CSM3V3fp8yEwhFCoQhI2VERWmA0CqwmY6wtgRDxIaOaJBKNEgpHiUQlkWiHIANMRiMWi6mPDuLBrNGXvndTqm8oNem1i78xi3p4+UuSzTWxkPEHtlfia23G33aa1XfcnlZY5NI0J6XEkHA/aZo2LIRTKrUNHjyByKDX9Rv0IAkhxN+AHVLKtAkaQohVwCqA8RMmLnxx38Gc7HcgiEQ1gsEQSInZZMBsNGAxGTCZjAjRUV6/YyEMR2Jl8Y2DUBY/eVEhHilWWlrKvffek1RhIFwcYO6H5tEsxtPmCzJpjJtXn/pjbLFvtPHqo9WkE06pi51ec08XUqlP8Lks/hp7SAghNQ2LyYDLYabQacv6QSEd0ahGeyBEqy9IMKwRDGsYTQbsVksW12owGx127jPxnKYm+mZThFf/TuY596f9SfbHkftoyq7ofqbzpo1DdpR26k5Y5CK4IVEb0xlOGlQqzW1+Tja143YNXj+9nLV8F0I8LoSoTvPv4wnbfAuIAPdnGkdKWSGlXCSlXFQ8oiTb48gLJqMBp8OG02nHYrWC0YQ/Kmj1R2j2RfAGJVFMGM0WbHYbLpc9ywWvP6Qm4HZSXV3N4sVLOou4/us+Rs8v4XCjiTZfkFFFdkyed2LCqcmGudGesSW6vkhVVe2O+6BWrlyZ5Hz3eLzx14n0tYV6MBTG6/UjI2HGlziYNXkE0ycUM7rY1S/hBGA0GnA7bUwcXcg5E4o5b8oIxhXbkdEIvnY/Pn8QTcs05761qe87XZNdU0knnNLdD4kVK7rOufO9gUiSrajYmnIvZJ/onS3hSBSPt50Sl4XfbPpBXDhBzK+U6cE79TfaH+G0evVqNE1j9erVbNmypdv9DmWKC+xMGePG423v5reQH3r89Uspl3X3uRDieuCjwAfkcLw6WSCEwGQUkFeHYnqHfeoT9re3fIPdtS8iRszGFLWAFsQRamTvnt2YWqxct+zG+BiQuOB1XYy3bdsW18p002FPIcLJY2am06SqUVJgo3hsQb/qEGb7ZCyEoMBpo8AZi+pq94doaPXj9UcwGo3YU0JwBy6AIHOH3e4SaJOjKbtvfJhw1GnmPHDV2weyt5aUEm97gAK7iSmTRrB27Zp++5V6gxCCoqKiJI1J90kVFRUNSw0KwGm3MGNiMYdrmrHabEMmeKK/UXxXEQuKuFxK2Z6bKZ3tdN+eIVOFAoiZ9X7yh//BOGIG/qgFg9BoOvICTaF2jK1Wrrvipnjggh72nLivdBWkU5sT6otYf0KTff4gRiEZX+Lsd1AK9M+34LBbmGyPBcu0eYM0tPkJRiR2myXhR9ofR35X0vl/Nm3azKFDB9P48TqLtqZ7CMj+OqQ35yaSm/Yn/a0qnxl/IIgByTnjC7F29ErLh7BYv3590gOQvt/hKpx0zCYjMyaO4J2TLUhpwjxgLX+yp18+KCHE24AVaOx460Up5a09fW8oBkkMBXp2XGd2ol91/TI2bPsx1e+2YR95DmaTgasvPZe7f/0zjG0WzKcc7N2zp+MbycEOmReNngIEeudnCARCSC3KuJGunAgm6D46q69+gXAkyukmLx5/BIvFnOPGkd37tRIrQKQGwWQOZMh0HTK/Pxg+qFz5J2PBMmHGjXBS5O5aqmy4Jc0OdaSUHK1tISoMAxbsNSiJulLKc/vzfUVidFusjl2mdtmpJiCXy8WDD/6ZxsYm6kQtr57eQ0vYjn3kKECyfMlU9u15HqMnJpwEghtuuJHS0tIsFjzobhHT55etZhGORAkGgowZ4UwK/c0FAxH2azYZmTi6ECklDS3tNHn8CIMBe04inbLXMLI3L3Z93dPDTu6Sc9ORG+1M0yS+9gCFTjPTxozIaALur19JkYwQgmkTijl+upVAIJTXyhNDw9B4FtK16nfsyaWkpCSt41oIQzzzvry8HK/XS2NjE+MWjKTV3shD/3yGhpATgPcvnMLB/bt5+fF9fOF9/8HePXvieTGVlZXxfSZXOe90aMfITYCA1+fHLDRmTS7JuXDSSRRSOrkwuQghGFXsZNbkEsYU2Qn4A/j8wRw4wkUXM17y4i2Ttu3+dTp6qmgvc1i9veu++1NVXqfdHyIcDnHuhCImjOqff1LRNyaNKcRtM9LuD+VtDmddLb6hQXIx0XR5Kzq6aSS1/tqdd97Jx7/yYd5pexuTowT3lMUIYeCiOeO4cOYYXn/mIMW+0axb25ksu2HDxi6VxrsPI4/NtS/hwqFwhHAozKTRbpwDUKUjkcEM+/X6g5xqaicS7U9juPS5S7p/MPEa9D1faTDzt5LpT46Vbs4bX+KkcAh0HlBAfbOPJm8op50XVC2+IY2IL0apPYMSG9FBbLFas6Y8/kQ8a9ZM3jx0kCcO/5N3297BaHXjnrQAIQyUnjMSz6lD7H3nFBvLfxZfPPXFQd9nom/g4MFDbNq0OYuabZ1z74l2fxCHxcD0ySMG3NySqwoB2eKyWzl3gpV2f4jaJh/haKzVdvb7SC7OCl3vAT2HSdO0pIi41OjJ7h8WBreVRSJ9iXyMahr+9iDFLgvTxw78faPInlHFTgwGQV2Lf4BLtnVFCai80bUDKSR3RU18Al6zphyAByofIDTBy7stLRjMdlyTFyOMZoJtJ2k/Wc9IRvP0A8+zSXYVOlLKLkVg9Tp7uYi20n0GE0YO3tNvvsJ+HXYL506w4A+EqG30EYxoWRbgTF9MNlFI6T7GzZs3x697769Rtn6ggUqoze7BRkqJzx/EbhKcO7FoWLVZP5soKXRgNBiobWrH7Rw8zVaZ+PJG+r5J3ZncApF2Fn92PpotijBaKJh2CUaLk/EjXdj8R9n36P54nlNq9YGecpd0+hptFY5EiITCTBufn0Um35FcgWCY2gYvgawFVaIg6FpeKvFarVlTnhTq31l+KvPY2UTp5bvckT8QREhtwIo1K3JPmzfAiQZfv6tO5KyShCKXdAYipPqgdHPPwYOHEpzJnY7rlkAjt//iS2i2KBiMuCcvwmhxEgm0YW0/ytdWfp0vLL+JggJ3UpFWIB5kkRr0sGZNObNmzUzati/VIPyBECYhmTFpRN6egPMdyWWzmpk+oZjpYwuIhEJ4ff4esvI7hVOqprNp0+akduiJwkn/PHaNUsfvTYBLz4EUA0UoHKviMarAxoxJJUo4DSMKXDYmjXLh8foHZX/KxDdIpD6tut1uSkpGMG7c+B6qfguOtx3le/d8i5eqXmTM2LH4LJMx2YuIhtqx+4+y/+9N/FPbmVRnLRHdtJPaR2jTps0ZTYrZmvk8vgCjC62MLHL29xSdEeiCKhAMc7LRS3tYw2G3ZYhC6z6Uv7y8PEnDTdSsDhw40G3KQDZVxAcqoTYT4UiUYDDECJeV0crPNGxxO61MNsCx0x7croHtp6cE1KDQtfyLHia+fPmVHduIBEHV+cP9wdbvcdh/CKvVyuLFS6gPFxIOmUELM60wyAXjP4RhnInq6mpuueUWfD5fUpuMlStXZhA6/asGIaXE6wswaZQLt3NwKyEPB2xWM9PGFxMMR6it9+ALazjs1pSK6pmvge6DSkTXrKSU8ZQBIUQ/AlwGJ5AiEtUIBIIUOi1MHZ05n0kxfHDZrUwZI3j3tGdAi8wqH9Sg0buw33A0zEu1z/DbP99N1e7dLF68BPuY2VQfaUBqETyHd/PpuR/ju+u+FzcXJjYY3LZtW5fOurfckr5FeW+d5FFNo709yPTxBdiGcC+ZoUQoHI1pVMEINltqo7ju6/Kl+g9j1SZEDsLIBzYUPRLVCPiDuB0mxpa4h0x9N0Xu8PlDHD3toaCXQmrQ2230hrNTQEG25V9aAk08c/xJ2kKtgGTnzsc5cKQJx+iZSKkxzRbmPDGFBysfin8ndSFLfb97p3r2RKIaoUCQcyYWqwWnD0SiGqcaPXjaI5gt5ozV2rMLYOhNOaGei9PmqtxRJBIlEAzhspsYpwTTGY8/EOKdU73TpJSAGnL01OAvthi81fQmL596iYiMxrd77Ug9z7x6Aokk8NabPH/PP7AYLRkWp9z2aEokHImiRcJMH1+szDT9RNMkDa0+mjxBDBnLKHWn3Wav/WQSdocOHWTmzFk5i+LTk2wLnRbGFDsxKsF01qALKZczmwhWlag7xOhcTEpLS+OtELZv3x4vm2N3Wzn/qpnUemuSvvnmuw088+oJAMKHTiBfb+bnm3/RZQ96Qm9qnlNuKlTHFh+0COdMKFbO7RxgMAhGF7sYXeyi1RugvtVPsMNPlVnj6Cqceg5w6b79RSzPKvuEWn1MfRspJe3+EAY0it02Ro51q/vjLMRuszBtrLtXQioblIAaFGLO8MTKAXpoeXV1NfuO7WXhxy6gyOsgcXF463gTT718DIGJ2YYxVPz+513C0xNbgb/88t5+R+WlIxSOYNCiTB2vhNNAUOiysXnDj2lqbuGr3/wv2tuDSGHgN7/6JQUF6bSZ3rXZ6F20XvfXV9fGbrv9DsKhCFaz4JcbfsDIEYW9bpuuOLOw2yxMH1fA4dpW3C57TtYKJaAGic6w31iOzvbt29EsUcJj2ll49QKuWP5+EheHwzXN/HvvMYyaiznGCfzqm7ehR/rptfT0EGN9AUruJ5SbCtXBUBizkEyZUNy3A1f0iJSSlpYWfvHzLRhErDzT7avXcM/vK7nuxi/S7g9ityV3bO5dOaHcROuFIxEaGlu4/w+/J+hr4Reb7mLdurX86hexMlOqzYXCZjVzzvjCnAkp5YPKA+1hLxd/fAnRwiAI+Oa3vkniYnH0ZCuPvfQOlvA4bnjvR/n88lKSF5PENh2J7wn6EpWXiWAojMUAk8cW9un7iuzJVPB206ZNeP0hmj1BAqEIGgKrxdzLhOi+RetJKQmEwkTDUcwmgdtuobjAxte+um5QCvMqhi+BYJgjJ1txOdMLKRUkMQQJRAK83vAqv6jcwu7dL8XfX7xkCcuXLwMEx0638dgLJ3AEZnHtey/hug/OycsPXwmnwUdKiSEhT0rTtC7XPhiO0OoN4PNHCEU0ohJMRgNmsymD7yr75oShcIRwOIrUNExGgdVspNBlxe2wJgXFZDNPhaI7IaWCJIYQUS3K6437eaPhAI88+s9YXlOHUNq583Gqdu8GYPYFF/H08z4K2xfzsfecx3UfPJ/k65pOc8o9oXAEs5BMHls0oPtRdKJrUImUl5d30UysZhOji13QYXGNRDX8gRBef5hgOEQkKtE0SVSTsbtFCMxWO5+65rOsuvXLeH0BVt3yJfzBCCazBb8/iNEgMBsNFNhMuIrt2KymjAIn23kqFDarmenjCzlSm1mT6gkloAYBIQRHW44Q1mJJmokaU+wvIEp44akCHOEpfHDJNCI1L3DjjZsSytl09nN6z3suGbBinqFwBIPUmDK+aEDGV3SlPy1DTEYDbqctbYVpKWOC6n++XY6mxfxDQoDRYOB/N67v9YIx2K1NFMMfmyUmpA7X9M0npQTUIGAQBuaNXsCzJ57kssuWkliI02KwsGzxDTz0Tx9GTfCRS87hxo+UsmnT412aC+rRe7Ew9Vy1RegkHIkgtChTlXAaVAaqZYgQApMxd/dIvlqbKIY3NouZcyYUcqSmFVcvhZTyQQ0i/zz8V5oCjTzzzC5CgTB3rLiTU+8WcvdfX0NKmF7o46df+3zHBUzfjqOz82qO66VFokQjYZXnlEfy3TIkW4bLPBVDi1A4ytsnmnE67cyYPFq12xhqzB+zCIMQONoLeHX7G/zu7lep+EtMOLW89RTjTKcTfE6iS9sMYGCEU1QjHAoxXeU55ZV8twzJlnTzTH3QzceDr2JoYzEbOXdiMe3tgay/kxcBdaqpHU97MB+7zivjXBP4+IzP8tM7f8ZFH76Dd9tLAGh+8zE+cvHkLhUANmzY2GWM2Hu5+/FrmiQYCHLOBFW+SNE31q9fT3l5eVwo6b4qlbirSCUmpIp66JXWSV4ElD8Y4au/eop3TrbmY/d5xSit/PT+3ZwKFiK1CA2vPojnWFVa4ZSuoWFlZWXOhJSUkvb2AOdOLFZ10xR9Qk8y3rJlS1xI6YETLS0tSpNSdMFsMhIN+73ZbJuXIAmr2Uhdczvf+N+nue1TC7j0gon5mMag0+wJ8JNtL3LoeDNGEaV27wMEm48DXWvmFRS443X7um9o2He8Pj/nTlBVyRV9JzFQYsuWLfGIPpW8q8gFeQmSmDN3vvzIV37Bky8fA+Bjl57LyivnYDaduQvlwWNN/PSBl2hqC2A1hHnnmXv4zNXLcblc7Nr1dFKTwU2bNndUlF7V8e1OzSr5dd/xeNuZNrYAu02121b0H5W8q+gNQoi9UspFPW2XF4kgBNz26QX8x0fnYTAI/vbs23yrYhenGrPS+oYdO6uO8p27n6GpLUCByc8l41v4zNXLKS8vjwunWbNm4nK54j16PB79XCQX/syN5hRg0ii3Ek6KnJApeVeZ9xT9pV8CSgjxAyHEfiHEPiHEY0KI8b34Lh95zzn8982XMbLIztsnmln7yyfZte94f6Y0pAiGIvzmoVf4zUOvEIlqjLO1Uv3PnxH2e1izJpZ9rwungwcPcffdd+ekYVwyyYuErz3A2GK7atOuyAmpybuaprF69eokn5RC0Vf6ZeITQhRIKds6/n8HcL6U8taevpeaB+X1h/j1g6/w4mu1ACydP4mbPjKXgmG8iB6uaeFnf6yipt6L2WTglo/P54qFk9MW7VyzppzFi5fE38tVk8HURnX+QJDNd/2QMSNcKsJKkTPWr19PS0tL3OekC62ioiJ1nynSkq2JL2c+KCHEN4DJUsov9bRtukRdKSU7q45yzz8OEApHKXBa+OLHLuC9cycMK1t2VJP89Zm3qHz8dSJRycRRbsqvXcS0eHWG5I63VVW707Zp778GlVwk9LbbbmfDT37I1l9sUA5sRc5RybuK3jBoAkoI8UPgC0Ar8H4pZX2G7VYBqwDGT5i48MV9B9OOd7LRy28efIXqdxoAWHzeWL74sQsYVeTo1zwHg6MnW6l4eB9vvtsEwIcuns4XPlSK1ay3Ruja9kA376VWmu5aMaIvpY069le5A09LE/6Wk0o4KRSKvJMzASWEeBwYm+ajb0kp/5qw3TcAm5Tyez3ttKdSR5omeXzPUX7/z2ragxEsJgMfu/RcPnX5TOxWc0/DDzrtgTCVT7zBIy8cQdMkRS4rt316AQtmJZ62zG0PZs2aybZt2xDCAEhuuOFGAO677170Hk8bN27qiOzrXZHYSDTK3NJSfA2xiEkVXaVQKPJNzqL4pJTLpJSlaf79NWXTB4BP93XCSZMyCK5cMo0tdy7jvXMnEIpo/PmpQ3x5w04efekdIlEtF7vpN5GoxhN7jnLb5p38/bnDICUfung6vyhfliKcIFOb7hUrVrB06eUdwilGaWkp1dXVbNy4iUTBFovsy17j1TSNH/3ox/gaOgNPlONaoVAMF/obJDFDSvlWx/9vBy6XUl7T0/d6Wyz24LEm7nvkAAePxUxnI4vsXP3ec1m2eCo2y+DnGkeiGk+/cow/PXmQ083tAMycVMzNV8/nnAlFPXw7m463feuAmjrGf//wx1Rs+R9W33F7UmsEZeZTKBT5ZFB8UEKIPwOzAA14F7hVSlnT0/f6Us1cSsnzB2qofOINaupjOUIuu5kPXTydKxZOYcwIZ+8PoJe0eoPs2necR144HBdM40e6+Mz7Z3HZBZNyXMsuOZiit5F9Xq+fbff8Al+biq5SKBRDi0GP4usN/Wm3oWmSPW+e5KFdb8U1KoDzppSwdP4k3jN3Am5H7hJQw5Eo+w/X88Sed9nz5kki0dj5mjDSxTVXzObSeRMx5kQwJWpSXVtt9EaD8vr8TChxUuCyqegqhUIx5MhWQA27hoUGg2DJ+eNZfN443jjayKO73+Gl10/yxruNvPFuI3c/vI/pE4opnT6S0umjmDV5BE5b9oEVkajG8ToPBw7X8erb9bz+TgPBcDS2bwELZ43hAwunsPj88TkSTKn5SsSFU2lpKffdd2+Sua8nIdXuDzGq0EaBK9Zhdbi0cFAoFIpUhp2A0hFCcP60kZw/bST+YJiXXj/J068cp/pIPW+faObtE838ZddbABQ4LYwrcTG2xMkItw27zYzRIJASAqEIPn+YxlY/dS3tnKjzdAnCmDymgMsumMj7F0xmRIE9x0ci8Xi8SQJI76Ib65xLXHC53S66E07+YAi33cjIooE3dyoUCsVAM+xMfD3hD0Z4891Gqo/UU32kgXdPtRKK9C7qb0yxg/OmlnDBuaOZe85oRhTYBmSunXQNiigrK2Pt2jUJ0X3d50GFwhFMaEwZVzSgM1UoFIr+csaa+HrCbjVx4cwxXDhzDBDzWTV7Apxq9FHb6KXNF8QfjBDVNAQCq8WI02ZmRIGdUUV2Joxy4+iFSTA3CFwuV9I7a9euSahqfgvdCadIJIqMRpg8oXiA56lQKBSDxxknoFIxGAQlhXZKCu3MmT4y39NJi5Qau3Y9nfTeypUr4xUm0rfZiGlUmiYJBkPMnDRC+ZcUCsUZxZnbgGnYINm0qbOquY7+es2acioqKuJJu/p3Nm7cxNat/4vP5+ecCUWqXbtCoTjjUAIq73RWmNi2bVvSJ7EKEyIeRJFaWaK+oYVp4wowm4zph1YoFIphzBlv4hsOrFp1C1JqbNq0Oel9rzeWkKxH8W3fvj0eSPHJT3+Gn/7gG6rpoEKhOGNRGtSQQMZbbqxYsYI9e6pYsWJFgtbUKaQAolLwo+99HbdjoKMLFQqFIn8oATUkyFxINpb7RFxQaRhoqz/Bf333G6roq0KhOKNRJr4hQiyUPDHXScS1Jt3n9JnPlvGNr97JT3/4XbZs2QKgir4qFIozFiWghhSpgib22u128dlry1h35x1MGFXA5s0xX1VRUZESTgqF4ozljKskcSYSjkSR0TDTxxfHBVJi0VdVEFahUAwnctawUJFfIlGNaDhZOEFn0df169cnNSHUW2qodhoKhWK4owTUECaqaYQCQaaPT2/Kk1LS0tLCli1b4kJKb0rY0tKigigUCsWwRvmghiiaJvH7g8yYWIzRmP45QggR90dt2bIlHjihOuYqFIozAeWDGoJIKfF6/cyYVJxVlQgpJQZDpxDTNE0JJ4VCMWRRPqhhipQSr8/P9AmFWQun8vLypPcSfVIKhUIxXFECaojh9QWYNrYAm6Xnlh+JPqfVq1ejaRqrV69O8kkpFArFcEX5oIYQHm87U8YUZF1fTwhBUVFRks9J5UgpFIozBeWDGiJ4vH4mj3Hhslt7/V2VB6VQKIYTygc1jPB4/Uwa1TfhBHQRRko4KRSKMwEloPKMx+tn4kgnbmffhJNCoVCcqSgBlUc8vphwKnCpthkKhUKRihJQecLj9TOxRAknhUKhyISK4ssDus9JmfUUCoUiMznRoIQQ64QQUggxMhfjncl4vX4mjVbCSaFQKHqi3xqUEGISsBw41v/pnNl4vH6mjHHjtGeX56RQKBRnM7nQoDYDXyPWDlaRBiklbZ52po1VwkmhUCiypV8alBDiaqBGSvmqyr1Jj6ZJfD4/50wozKp8kUKhUChi9CighBCPA2PTfPQt4JvAldnsSAixClgFMH7CxF5McfgS1TTa/YGsq5IrFAqFopM+lzoSQswFngDaO96aCNQCS6SUp7r77tlQ6igSiRIKhjhnYjGmDP2cFAqF4mwk21JHfTbxSSkPAKMTdngUWCSlbOjrmGcK4UgELRJhxqQRGAzK9KlQKBR9QeVB5ZhAKIxZSKZNKFY18RQKhaIf5ExASSmn5mqs4Uq7P4TbbmT8SHe+p6JQKBTDHqVB5Qhve4BRBVZGFjnzPRWFQqE4I1ACqp/EWrQHVNFXhUKhyDFKQPUDPYx8+niV46RQKBS5RgmoPhIKR5DRCLMmjsCowsgVCoUi5ygB1Qf8gRB2s2DiWBWpp1AoFAOFElC9xOPzM7rQpoIhFAqFYoBRAipLoppGe3tQVSNXKBSKQUIJqCwIhsKgRZk1qVj5mxQKhWKQUAKqB3ztAYocZsaOLMz3VBQKheKsQgmoDEQ1DX97gImj3bgdqvutQqFQDDZKQKUhEAxhQDJzkgohVygUinyhBFQCelWIUYU2RhWrKD2FQqHIJ0pAdRAKR4iEwkwfX6CqQigUCsUQ4KwXUFJKvO0Bip0Wxo4doRJvFQqFYohwVguoYChMNBxh+tgCbFalNSkUCsVQ4qwUUJom8bUHKHFbGT2uQGlNCoVCMQQ56wRUuz+EySA5d0IRFrMx39NRKBQKRQbOGgEVCkcIhcKML3FSqPo2KRQKxZDnjBdQkahGwB+k2GVhugqCUCgUimHDGSug9OKuLruJGZOKMamEW4VCoRhWnHECKtblNojTYlR+JoVCoRjGnDECKhLVCASCOKwmzh2vBJNCoVAMd4a9gNKDH9x2E5MmKlOeQqFQnCkMSwGlaRJ/IIhBQLHLSskYNwaDCn5QKBSKM4lhI6CklPgDIaSm4bCamDrGrao/KBQKxRnMkBZQkUiUQDCMEBK7xcSEEgcu1ZtJoVAozgqGjICSUhIKR4lEIiAlZpMBp9XE+BEFWC1DZpoKhUKhGCT6tfILIdYDNwP1HW99U0r5SDbf9fr8IMFgEJiMArPRQLHTjMvhwGpWAkmhUCjOdnIhCTZLKTf05gtWs4lzJxRhMhpUZQeFQqFQpCUvqooQYDapPCWFQqFQZCYXSUO3CSH2CyHuEUIUZ9pICLFKCLFHCLGnvr4+02YKhUKhUAAgpJTdbyDE48DYNB99C3gRaAAk8ANgnJTypp52umjRIrlnz57ez1ahUCgUwx4hxF4p5aKetuvRxCelXJblDu8G/p7NtgqFQqFQ9ES/THxCiHEJLz8JVPdvOgqFQqFQxOhvkMRdQoj5xEx8R4Fb+jshhUKhUCignwJKSnldriaiUCgUCkUiqvS3QqFQKIYkSkApFAqFYkiiBJRCoVAohiQ95kENyE6F8AAHB33HQ4eRxPLHzlbO9uMHdQ7U8Z/dxz9LSunuaaN8VWU9mE2S1pmKEGKPOv6z9/hBnQN1/Or4s9lOmfgUCoVCMSRRAkqhUCgUQ5J8CaiKPO13qKCOX3G2nwN1/Gc3WR1/XoIkFAqFQqHoCWXiUygUCsWQRAkohUKhUAxJ8iqghBC3CyEOCiFeE0Lclc+55AshxDohhBRCjMz3XAYTIcRPhRBvdjS7fEgIUZTvOQ0GQoirOu75t4UQX8/3fAYTIcQkIcSTQog3On7zq/M9p3wghDAKIV4RQpyV7YmEEEVCiD91/P7fEEJckmnbvAkoIcT7gY8D86SUc4AN+ZpLvhBCTAKWA8fyPZc8sBMolVLOAw4B38jzfAYcIYQR+BXwIeB8YIUQ4vz8zmpQiQBrpZTnARcDXznLjl9nNfBGvieRR7YA/5JSzgYuoJtzkU8N6kvAj6WUQQApZV0e55IvNgNfI9au5KxCSvmYlDLS8fJFYGI+5zNILAHellIekVKGgEpiD2lnBVLKk1LKlzv+7yG2ME3I76wGFyHEROAjwG/zPZd8IIQoAJYC/w9AShmSUrZk2j6fAmomcJkQ4iUhxNNCiMV5nMugI4S4GqiRUr6a77kMAW4C/pnvSQwCE4DjCa9PcJYt0DpCiKnAhcBLeZ7KYPMzYg+lWp7nkS+mA/XAvR1mzt8KIZyZNh7QUkdCiMeBsWk++lbHvouJqfqLgT8KIabLMyjuvYfj/yZw5eDOaHDp7villH/t2OZbxEw/9w/m3PKESPPeGXO/Z4sQwgX8GbhTStmW7/kMFkKIjwJ1Usq9Qoj35Xk6+cIELABul1K+JITYAnwd+E6mjQcMKeWyTJ8JIb4EPNghkHYLITRiBRTrB3JOg0mm4xdCzAWmAa8KISBm3npZCLFESnlqEKc4oHR3/QGEENcDHwU+cCY9mHTDCWBSwuuJQG2e5pIXhBBmYsLpfinlg/mezyDzXuBqIcSHARtQIITYJqVcmed5DSYngBNSSl1z/hMxAZWWfJr4/gJcASCEmAlYOEuq+0opD0gpR0spp0oppxK7aAvOJOHUE0KIq4D/BK6WUrbnez6DRBUwQwgxTQhhAcqAh/M8p0FDxJ7G/h/whpRyU77nM9hIKb8hpZzY8ZsvA/59lgknOta440KIWR1vfQB4PdP2+apmDnAPcI8QohoIAdefJU/Rihi/BKzAzg4t8kUp5a35ndLAIqWMCCFuAx4FjMA9UsrX8jytweS9wHXAASHEvo73vimlfCR/U1LkgduB+zse0o4AN2baUJU6UigUCsWQRFWSUCgUCsWQRAkohUKhUAxJlIBSKBQKxZBECSiFQqFQDEmUgFIoFArFkEQJKIVCoVAMSZSAUigUCsWQ5P8DrpObHp8bw4UAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "m.plot()\n", "plt.plot(X, f, lw=5, label=\"GT f\", alpha=0.4, color='green')\n", "plt.legend()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.9.7" } }, "nbformat": 4, "nbformat_minor": 5 }