From 77b46a6515b3bdaa77079e8f964322b54dc6d841 Mon Sep 17 00:00:00 2001 From: Kees van Kempen Date: Tue, 8 Mar 2022 17:22:51 +0100 Subject: [PATCH] 09: Task 6 is donish, but AB4 still looks crappy. --- .../9 Ordinary Differential Equations.ipynb | 49 +++++++++++-------- 1 file changed, 29 insertions(+), 20 deletions(-) diff --git a/Week 5/9 Ordinary Differential Equations.ipynb b/Week 5/9 Ordinary Differential Equations.ipynb index 951e6bd..7a38d73 100644 --- a/Week 5/9 Ordinary Differential Equations.ipynb +++ b/Week 5/9 Ordinary Differential Equations.ipynb @@ -635,7 +635,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 9, "metadata": { "deletable": false, "nbgrader": { @@ -650,32 +650,41 @@ "task": false } }, - "outputs": [], + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtgAAAHkCAYAAAAEg1hOAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAACU7ElEQVR4nOzdd3xUVd7H8c+dPpNJ740kJJCEXkKvgiKuYHftoCiu3XXVZ8uzLvi467q7ruva1q5gAV1RVOxK79JLAoQSAkkgvWcmU87zx4RApAVMI/zer9d93Zlzy/wmy+I3h3PP0ZRSCCGEEEIIIVqGrr0LEEIIIYQQojORgC2EEEIIIUQLkoAthBBCCCFEC5KALYQQQgghRAuSgC2EEEIIIUQLkoAthBBCCCFECzK0dwEtLSwsTCUmJrZ3GUIIIYQQopNbv359sVIq/KftnS5gJyYmsm7duvYuQwghhBBCdHKapu0/UbsMERFCCCGEEKIFScAWQgghhBCiBXXogK1pWldN097QNO2j9q5FCCGEEEKI5mjzMdiapr0JTAIKlVK9jmmfCPwb0AOvK6WeUkrtBW6XgC2EEEKIk3G5XBw8eBCHw9HepYhOymKxEBcXh9FobNb57fGQ49vAC8DsIw2apumBF4GLgIPAj5qmfaaUymyH+oQQQghxDjl48CD+/v4kJiaiaVp7lyM6GaUUJSUlHDx4kKSkpGZd0+ZDRJRSS4HSnzQPBnYrpfYqpeqBucDlbV2bEEIIIc49DoeD0NBQCdeiVWiaRmho6Bn9C0lHGYMdCxw45v1BIFbTtFBN014G+mua9vuTXaxp2p2apq3TNG1dUVFRa9cqhBBCiA5GwrVoTWf656ujBOwTUkqVKKXuUkolK6X+eorzXlVKZSilMsLDj5vrWwghhBCiVen1evr169e4PfXUU6c8/+233+a+++5ro+pEW+soC83kAfHHvI9raBNCCCGE6PCsViubNm1qtfu73W4Mho4S28TpdJQe7B+BbpqmJWmaZgKuBz5r55qEEEIIIX6WxMREiouLAVi3bh1jx4497pyioiKuvvpqBg0axKBBg1ixYgUAM2fO5JZbbmHEiBHccsstbVm2+JnaY5q+OcBYIEzTtIPADKXUG5qm3Qd8g2+avjeVUtvbujYhhBBCiLNRV1dHv379Gt///ve/57rrrmvWtQ8++CAPPfQQI0eOJDc3l4svvpisrCwAMjMzWb58OVartTXKFq2kzQO2UuqGk7R/CXzZxuUIIYQQohN5/PPtZOZXtug9e8QEMGNyz1Oe83OGiHz//fdkZh6dmbiyspLq6moALrvsMgnX5yAZzCOEEEII0UoMBgNerxfgpNO8eb1eVq9ejcViOe6Yn59fq9YnWocEbCGEEEJ0GqfraW5riYmJrF+/nksuuYR58+ad8JwJEybw/PPP8+ijjwKwadOmJsNNxLmnozzkKIQQQghxzjoyBvvI9rvf/Q6AGTNm8OCDD5KRkYFerz/htc899xzr1q2jT58+9OjRg5dffrktSxetQFNKtXcNLSojI0OtW7euvcsQQgghRBvJysoiPT29vcsQndyJ/pxpmrZeKZXx03OlB1sIIYQQQogWJAG7hXg9Hj59+s9kLlvU3qUIIYQQQoh2JA85tpBti79n94+rKcjeSbchwzGazO1dkhBCCCGEaAfSg90CXE4Hq/77Hv6h4dSUl7Hpmy/wej3tXZYQQgghhGgHnaYHW9O0ycDklJSUNv/s3T+uprqslOtmPMWKD99l6btvsuy9t7EGBOAXGIQtKBh7cAi2oGD8AoPxCw7GLygYv6AQ7MHBmKy2Nq9ZCCGEEEK0jk4TsJVSnwOfZ2RkTG/rz04fOZaQmDgiu6Zw2cN/YOfKZdRUlFFT3rCVlVGSd4Da8nK8Hvdx1xvNFuwhIfgFhzSE7oYtJBR7cCh+Ib7XMuxECCGEEKLj6zQBu71FdvX1nFv9A+h38aUnPEcphaO66mjwLi+jurSkIYSXUl1WyuG92ewpLcVd7zzuerOfH/bgUPxDw3zhOyQM/5BQ7KGh+IeEYQ8Nw+JnR9O0Vv2uQgghhBDi5CRgtyFN07D6B2D1DyAsPuGk5ymlcNbWUFNWSlVpiS98l5ZQVVpCdcNWtH8fNRXl8JN5zA1mM/4hYfiHhuIfGoF/WBj+oWEEhIbjH+bbTBZrK39TIYQQQpytxMRE1q1bR1hYGMOHD2flypUAPProo3z55Zf84he/IDk5GZvNxpQpU5p9X7vdTnV19XHter2e3r17N76//vrrGxfKOZG3336bdevW8cILL5zBt2qqoKCA6dOns2DBAgDWrl3LI488wuHDh7HZbAwcOJDnnnsOm+3Ew2iPrWHmzJnY7XYeeeSRk37e/Pnz6d69Oz169DjpOS+88AI2m41p06ad9fc6QgJ2B6RpGhY/OxY/O6FxXU56nsftpqa8lKoSX+iuKilq2BdTVVLE/q0bqSkrQylvk+ssfvbGsB0QFkFAWDgB4RG+1+ER2AKDpBdcCCGE6ACOhGuAV199ldLS0pOuCHm2rFYrmzZtatF7HsvtdmMwNI2czzzzDNOn+0b1Hj58mGuvvZa5c+cybNgwAD766COqqqpOGrDP1Pz585k0adIpA/a0adMYMWKEBOzznd5gaAjIESc950gIrywuoqqkmMqiQl8ALy6kqqiQvKztOGtrmlxjMJrwD48gMNwXuAPCIwkMjyAwIorAiEisAYEdIoB7vR50upb9S0YIIYQ4Uzk5OUycOJGhQ4eycuVKBg0axG233caMGTMoLCzkvffeY/DgwZSWljJt2jT27t2LzWbj1VdfpU+fPpSUlHDDDTeQl5fHsGHDOHaV7SO9zpdddhnV1dUMHDiQ3//+92RlZTX22u7Zs4d7772XoqIibDYbr732Gmlpaezbt48bb7yR6upqLr/88jP+Xsf2pK9bt45HHnmExYsXNzmnqKiIu+66i9zcXACeffZZRowYwcyZM9mzZw979+6lS5cuzJkzp8l18+bN489//jMAL774IlOnTm0M1wDXXHMNwEl/Zifz2muv8eqrr1JfX09KSgrvvPMOmzZt4rPPPmPJkiX8+c9/5sUXX+Tee+9tvGbr1q3s3buXhIQEEhMTWbt2LYMHDz7jn9exJGB3cs0J4c7aGiqLCqksLmzYF1FReIjKoiIO7d2No6qyyfkGs5nA8EgCIyIJjIwiKCKKwMioxgBuNFta+2tRnJvDh//3B8bffg+pw0a2+ucJIYQ4R3z1Ozi0tWXvGdUbLnnqlKfs3r2b//73v7z55psMGjSI999/n+XLl/PZZ5/x5JNPMn/+fGbMmEH//v2ZP38+CxcuZMqUKWzatInHH3+ckSNH8qc//YkvvviCN95447j7f/bZZ9jt9sae5pkzZzYeu/POO3n55Zfp1q0ba9as4Z577mHhwoU8+OCD3H333UyZMoUXX3zxpLXX1dXRr1+/xve///3vue6665r1o3nwwQd56KGHGDlyJLm5uVx88cVkZWUBkJmZyfLly7Famw5N3bdvH8HBwZjNvskbtm3bxtSpU094/5P9zE7mqquuauwZ/+Mf/8gbb7zB/fffz2WXXcakSZMag/uRe7z44ossWbKEhATf0N2MjAyWLVsmAVv8fGabH+EJSYQnJJ3weH1dLRVFhVQWHaai8OhWWXiIA5nbcDnqmpzvFxRMYGQ0QZFRBB3ZR8UQFBWNxe7/s3u/lVIsmvUqdVWVLJ8zi26Dh6Fr4X8uE0IIIc5EUlJS4zjmnj17Mn78eDRNo3fv3uTk5ACwfPly5s2bB8C4ceMoKSmhsrKSpUuX8vHHHwNw6aWXEhwc3OzPra6uZuXKlVx77bWNbU6nb6KEFStWNH7eLbfcwm9/+9sT3uPnDBH5/vvvyczMbHxfWVnZOM77sssuOy5cg2/8dXh4eLPuf7Kf2cls27aNP/7xj5SXl1NdXc3FF1980nNXrFjBa6+9xvLlyxvbIiIi2LFjR7NqOxUJ2OK0TFYb4V0SCe+SeNwxpRR1VZVUFB6i/PAhKg4VUF54iIrDh8jdtpnMpQubnG/28yMo0he2g6NjCI6KISgqhuDoGKz+Ac2qZ/faVeRu20JyxhD2rFtD5rJF9Bp7YUt8VSGEEOe60/Q0t5YjvbEAOp2u8b1Op8PtPn6K3pbi9XoJCgo6aUD+OZ1aBoMBr9f3HJfD4Tjp569evRqL5fh/vfbz8zvhNVartcn9evbsyfr1689qGMtP3XrrrcyfP5++ffvy9ttvHzek5YiCggJuv/32xn8ZOMLhcJzwl4IzJQFb/CyapmELCMQWEEh0Supxx131TioO+8J3+aH8xv2h3TvZtWp5kwcwLXZ/X+iOjj1m84XvI8NOXPVOFr/zOmHxCUx+6Pe8/8eHWf3xXNJHjkVvkD/OQgghOq5Ro0bx3nvv8dhjj7F48WLCwsIICAhg9OjRvP/++/zxj3/kq6++oqysrNn3DAgIICkpif/+979ce+21KKXYsmULffv2ZcSIEcydO5ebb76Z995774zrTUxMZP369VxyySWNvcg/NWHCBJ5//nkeffRRwDf04tjhJifSvXv3xl59gPvuu4/Bgwdz6aWXMmTIEAA+/vhjRowYcdKf2clUVVURHR2Ny+XivffeIzY2FgB/f3+qqqoAcLlcXHvttfztb3+je/fuTa7ftWsXI0aMOGX9zSGJRLQqo8lMWHzCCacldLtcvp7vQ/mUFeQ37PPI3b7luJ5v/7BwQmLiUEpRWVTItY89id5gYPi1NzL/70+QuWwhvS+Y0FZfSwghhDhjM2fOZNq0afTp0webzcasWbMA3zjjG264gZ49ezJ8+HC6dDn5DGIn8t5773H33Xfz5z//GZfLxfXXX0/fvn3597//zY033sjf/va3U/YO/3QM9sSJE3nqqaeYMWMGt99+O4899hhjx4494bXPPfcc9957L3369MHtdjN69GhefvnlU9br5+dHcnIyu3fvJiUlhcjISObOncsjjzxCYWEhOp2O0aNHM3HixJP+zE7miSeeYMiQIYSHhzNkyJDGUH399dczffp0nnvuOZ599lnWrVvHjBkzmDFjBgBffvklMTExrFixosn49rOlqZ/Mo3yuy8jIUOvWrWvvMsTP5HI4KGsI3KX5BynLz6M0P4+ygjy6DR7GxHseAnxDVN77w2+oKS+l34RLCYmLJzS2C0GRUTIuWwghzhNZWVmkp6e3dxniDHzyySesX7++cSaRjmDjxo0888wzvPPOOyc8fqI/Z5qmrVdKZfz0XOnBFh2S0WIhIrErEYldT3mepmlcMHU6X77wNMvnzm5s1xsMBEfHEhLXhdDYeMLiuxAal0BQVLQMJRFCCCHa2ZVXXklJSUl7l9FEcXExTzzxRIvcS3qwRafhrK2lNP8ApXkHKck7QMmB/ZTmHaS88FDjipc6vYGQmFhC47oQFp9AaLxvHxgZJXNqCyHEOUp6sEVbkB5scV4y22xEp6Qe97Cly+loErqLD+ZyaM8udq5a1niOwWQmNC6esPhEwrokENYwa4pfUPOnShJCCCGEAAnY4jxgNFuI7JpCZNeUJu31jjpKDuRSfHA/JQf2U5S7n32b1rF9yfeN59gCgxrD9pG5wkPj4tEbjG39NYQQQghxjug0AVvTtMnA5JSUlNOeKwSAyWIlulsq0d2a9njXVlZQnJtD0f4cig/49pu//RK3qx4AnV5PSGw8EQ2BO7xhrHhz5/EWQgghROfWaQK2Uupz4POMjIzp7V2LOLfZAgLp0qsvXXr1bWzzej2UFeRTtH9f45a7bTOZyxY1nmMPDWt8MDMiKZnIxGT8w8J/9sqVQgghOj673d64giHA22+/zbp163jhhRfasSrRXjpNwBaiNel0ekJj4wmNjSdt+OjG9trKCgpz9lKUs5fChm3fhnWNC+hY7P6+sJ2U7Nt3TSEoMlpCtxBCCNGJScAW4mewBQSS2Kc/iX36N7a5nA6K9uf4Ave+3RTm7GXDl5/iaVgq12zz8/Vyd00hqmsKkcndJHQLIUQnVlRUxF133UVubi4Azz77LCNGjGDmzJnY7XYeeeQRAHr16sWCBQsAuOSSSxg5ciQrV64kNjaWTz/9tEWW8BZtQwK2EC3MaLYQ0z2NmO5pjW0et4uSgwc4vHe3b9u3m03fLMDjcgFg9vMjMskXtqMaNv9QGV4ihBDnip+uhlhaWspll10GwIMPPshDDz3EyJEjyc3N5eKLLyYrK+uU98vOzmbOnDm89tpr/PKXv2TevHncfPPNrfkVRAuSgC1EG9AbjI3js3uP8y3p7nG7KTmYy6E92Rzem82hPdmsXzAfr8fX020LDPKF7ZTuRCd3JzKlO1a7f3t+DSGE6PD+tvZv7Cjd0aL3TAtJ47eDf3vKc6xWK5s2bWp8f2QMNsD3339PZmZm47HKysom47VPJCkpqTGwDxw4kJycnLOqXbQPCdhCtBO9wXB0tcrxFwPgrq+nKHefL3TvyaZg9y72blzXuFBOUFS0b67vbr75vsMTk2TKQCGE6OC8Xi+rV6/GYrE0aTcYDHi93sb3Doej8bXZbG58rdfrqaura/1CRYuRgC1EB2IwmY5bLMdZW8Phvbsp2L2LQ7t3krt9C1nLFwOgNxqJSEomplsa0d3SiO6WSkBYeDtVL4QQ7e90Pc3tYcKECTz//PM8+uijAGzatIl+/fqRmJjYOOZ6w4YN7Nu3rz3LFC1IArYQHZzZ5tdk2kClFFUlxRzavZP87J0c2r2Tzd9+yfov5gNgDwn1Be7uacR0TyciKRmDUXq5hRCivTz33HPce++99OnTB7fbzejRo3n55Ze5+uqrmT17Nj179mTIkCF07969vUsVLURTDf/03FlkZGSoI2OehDhfeNwuivbnkL9rBwXZvq2i8DDg6+WOTEohJjWdmNR0YrunYwsMat+ChRCiBWVlZZGent7eZYhO7kR/zjRNW6+UyvjpudKDLUQnoDcYG2cf4ZLJANSUl5G/K4u8nVkU7NrBxq8+Y93nHwO+sdyxqT2JTetBbFoPgqNjZcYSIYQQooVIwBaik/ILCqbb4OF0Gzwc8D1AeXjfHvJ3ZpK3M4u9G9ayfcn3AFj9AxrCdk/i0noSkZSMTq9vz/KFEEKIc5YEbCHOEwaTidjUdGJT0xmEbyx3WUEeeTsyG7bt7P5xNeCbyzu6expxaT2JS+9JVLdUjCbzqT9ACCGEEIAEbCHOW5qmERITR0hMXOPc3NVlpeTt2M7BrO3k7djOyo/eB6XQGwxEpXQnLr03cek9iUlNx2SRFcWEEEKIE5GALYRoZA8OIXXYKFKHjQLAUVNN3o5MDmZtIy9rO2s//S9rPvkATacjKrkb8T16E9ejN7Gp6ZistnauXgghhOgYJGALIU7K4mcneeBgkgcOBqDeUUf+ziwOZm3nQOZW1i2Yz9pPP/IF7q7diOvZmy49ehOb1hPjTxZUEEIIIc4XErCFEM1mslhJ7DuAxL4DAHA5HOTv2sGBzK0cyNzK+gWf8OOnH6HT+4aUdOnVh/gefYjpnobBZGrn6oUQonV5PB4yMjKIjY1tXEDm2Wef5c4778Rm8/0rn91uP+0y6WcrMTGRdevWERYWdlbX33rrrSxZsoTAwECUUjzzzDOMHz++hatsniNLzb/wwgt4vV5uu+029Ho9b7zxxglnvSovL+f999/nnnvuASAnJ4eVK1dy4403tnXpgARsIcTPYLRYSOjTj4Q+/QBf4M7bmUnu9i0c2L6FNR9/yOp5czEYTcSkpvsWzOndl8ikFJmlRAjR6fz73/8mPT2dysrKxrZnn32Wm2++uTFgtxS3243B8PNinMfjQf+Tv4v/8Y9/cM0117Bo0SLuvPNOsrOzf9Zn/FxKKe666y5cLhdvvfXWSaeULS8v56WXXmoSsN9///12C9i6dvnUVqBp2mRN016tqKho71KEOG8ZLRYS+w5g9I23ctNfnuHeN+dwxf88Rp+LLqG2soLlc2fz/v8+zEt33Mj8fzzBhq8+p+TgATrbgldCiPPPwYMH+eKLL7jjjjsa25577jny8/O54IILuOCCCxrb//d//5e+ffsydOhQDh8+fNy9SktLueKKK+jTpw9Dhw5ly5YtAMycOZNbbrmFESNGcMstt1BSUsKECRPo2bMnd9xxR5O/S999910GDx5Mv379+NWvfoXH4wF8PegPP/wwffv2ZdWqVSf9PsOGDSMvLw/w9Sbfd999jccmTZrE4sWLG+93ou+zZ88ehg4dSu/evfnjH/+I3W5vvP4f//gHgwYNok+fPsyYMeOUP9cHHniAkpISZs+ejU6nY+bMmTz99NONx3v16kVOTg6/+93v2LNnD/369ePRRx/ld7/7HcuWLaNfv37861//Iicnh1GjRjFgwAAGDBjAypUrT/m5P1enCdhKqc+VUncGBga2dylCiAZmmx/JA4dwwdTpTP3HC9z96rtc+uD/0H3YSIoP7GfR26/w9sN38+o9t/L1S/8ic9kiasrL2rtsIYQ4Y7/+9a/5+9//jk53NFo98MADxMTEsGjRIhYtWgRATU0NQ4cOZfPmzYwePZrXXnvtuHvNmDGD/v37s2XLFp588kmmTJnSeCwzM5Pvv/+eOXPm8PjjjzNy5Ei2b9/OlVdeSW5uLuBbcfCDDz5gxYoVbNq0Cb1ez3vvvdf4+UOGDGHz5s2MHDnypN/n66+/5oorrjjt9z7Z93nwwQd58MEH2bp1K3FxcY3nf/vtt2RnZ7N27Vo2bdrE+vXrWbp06Qnv/f7777Nhwwbmzp172t76p556iuTkZDZt2sQ//vEPnnrqKUaNGsWmTZt46KGHiIiI4LvvvmPDhg188MEHPPDAA6f9bj+HDBERQrQZW2AQacNHkzZ8NAAVhYfYv3UT+7duZs+GH9m+5AcAwuITSOjTn8Q+/YlN74nRLA9MCiGa59CTT+LM2tGi9zSnpxH1hz+c9PiCBQuIiIhg4MCBjT27J2MymZg0aRIAAwcO5LvvvjvunOXLlzNv3jwAxo0bR0lJSeOwk8suuwyr1TdN6tKlS/n4Y98KvZdeeinBwcEA/PDDD6xfv55BgwYBUFdXR0REBAB6vZ6rr776pPU9+uij/OEPf+DgwYOn7OE+3fdZtWoV8+fPB+DGG2/kkUceAXwB+9tvv6V///4AVFdXk52dzejRo4+794ABA9ixYwdr165lxIgRp63lVFwuF/fdd1/jLxy7du36Wfc7HQnYQoh2ExgRRZ/xE+kzfiLK66UwZ68vcG/ZyKZvv2D9F/PRG43EpvUksU9/Evr0JzwhSZZ1F0J0KCtWrOCzzz7jyy+/xOFwUFlZyc0338y777573LlGo7Hx7zC9Xo/b7T6jz/Lz8zvtOUoppk6dyl//+tfjjlksluPGXR/ryBjs559/nmnTprF+/XoMBgNer7fxHIfD0fj6TL+PUorf//73/OpXv2rS/uKLLzb2fn/55ZcApKWl8X//93/88pe/5JtvvqFnz56nrOVU/vWvfxEZGcnmzZvxer1YWnmmKwnYQogOQdPpiOyaQmTXFAZffg0up4O8rO3kbNnI/i0bWfreW/DeW9gCgxpnMkno0x9bgAwLE0Icdaqe5tby17/+tTHMLl68mKeffroxXPv7+1NVVXVGM3uMGjWK9957j8cee4zFixcTFhZGQEDAceeNHj2a999/nz/+8Y989dVXlJX5htiNHz+eyy+/vHFoRGlpKVVVVSQkJDS7hvvuu48333yTb775hsTERF566SW8Xi95eXmsXbv2tNcPHTqUefPmcd111zF37tzG9osvvpjHHnuMm266CbvdTl5eHkajkXvvvZd77733uPsMHz6c//znP0yaNIklS5aQmJjYOEPLhg0b2LdvH3D053zET99XVFQQFxeHTqdj1qxZjWPSW4sEbCFEh2Q0W0jsN5DEfgMBqCotZv8WX+/23o3ryFy6EDSNyKQUkvoNILHvQKK7pcrsJEKIDuXOO+9k4sSJjWOxm2PmzJlMmzaNPn36YLPZmDVr1gnPmzFjBjfccAM9e/Zk+PDhdOnSBYAePXrw5z//mQkTJuD1ejEajbz44otnFLA1TeOPf/wjf//73/n+++9JSkqiR48epKenM2DAgNNef2T2lL/85S9MnDiRI8/ITZgwgaysLIYNGwb4HpJ89913G4ewnMjkyZMpLi5m4sSJfPfdd8yePZuePXsyZMgQunfvDkBoaCgjRoygV69eXHLJJTz55JPo9Xr69u3Lrbfeyj333MPVV1/N7NmzmThxYrP+JeDn0Drb0/sZGRlq3bp17V2GEKIVeb0eCvfuYd/m9eRs3kjBrh0o5cVs8yOhT3+SGoK5PTikvUsVQrSBrKws0tPT27sMcYza2lqsViuapjF37lzmzJnDp59+2t5l/Swn+nOmadp6pVTGT8+VHmwhxDlHp9MTldKdqJTuDLv6Bhw11eRu28y+jevJ2bSOXauXAxCRmExS/4Ek9cuQ3m0hhGhD69ev57777kMpRVBQEG+++WZ7l9SmpAdbCNGpKKUo2r+PfZvWk7NpPXk7M1FeL2Y/PxL7DKDrgEEk9hsoY7eF6ESkB1u0BenBFkKctzRNIyKxKxGJXRlyxbW+3u2tm9i7cR05m9azc9Uy0DSiU7rTtf8gkgYMIiKxq8xMIoQQosVIwBZCdGoWPzvdh46k+9CRKK+Xw/v2sHfDj+zb+CMrPnyXFR++iz04hKT+GXQdMJiE3v0wtvL0TUIIITo3CdhCiPOGptMRldyNqORuDL/2RmrKy9i3aT37NvzIzlXL2brwWwxGE/G9+pA8cDBdBwzGP7T5U2sJIYQQIAFbCHEe8wsKptfYC+k19kI8bhcHs7azd/1a9mxYy76N64CXiEhMpuvAwSQPHExkUjLaMcsgCyGEECciAVsIIQC9wUhC734k9O7H2KnTKc07wJ71a9mzfi1rPv6A1fPmYA8OoevAwaRkDCW+V18MRmN7ly2E6EDmz5/PlVdeSVZWFmlpaeTk5JCenk5qaipKKfz8/HjrrbdITU1l7dq13HnnnYDv4eyZM2dy5ZVXtvM3EC1FZhERQojTqK2sYN/Gdexdv5Z9mzfgctRhtFhJ7NuflIyhJA0YhNXu395lCnHe6iiziFx33XXk5+czbtw4Hn/8cXJycpg0aRLbtm0D4JVXXmHlypXMmjWL2tpaTCYTBoOBgoIC+vbtS35+PgaD9H12VDKLiBBCtCBbQCA9x4yn55jxuOvrObB9C7vXrWbP+rVkr1mJptMRm9aDlIyhJGcMJSgyqr1LFkK0serqapYvX86iRYuYPHkyjz/++HHnVFZWEhwcDIDNZmtsdzgcMpNRJyMBWwghzoDBZCKpfwZJ/TO48HYvh/fuZve6NexZt5rFs19n8ezXCeuSSMqgYXQbPIzwhCT5D6cQ54FPP/2UiRMn0r17d0JDQ1m/fj2hoaHs2bOHfv36UVVVRW1tLWvWrGm8Zs2aNUybNo39+/fzzjvvSO91JyJDRIQQooWUHz7EnnWryV67irydmaAUAeGRpAwaSrdBw4hJS0enk9UkhWhpx/7T/bIPd1F8oLpF7x8Wb2fUL7uf8pxJkybx4IMPctFFF/Hcc8+Rm5vLfffd12SIyAcffMBbb73F119/fVz9U6dOZenSpVhkmtAOS4aICCFEOwiKjGLgpVcw8NIrqK0oZ8/6tez+cRWbv/2CDV9+itU/gOSMIaQMGkZC734YTKb2LlkI0QJKS0tZuHAhW7duRdM0PB4PmqZx7733Njnvsssu47bbbjvu+vT0dOx2O9u2bSMj47isJs5BErCFEKIV2AKD6D1uAr3HTaC+rpZ9mzaw+8dV7Fq9gm2LvsNotpA0YBDdhwwnqd9ATFbb6W8qhDit0/U0t4aPPvqIW265hVdeeaWxbcyYMRw4cKDJecuXLyc5ORmAffv2ER8fj8FgYP/+/ezYsYPExMS2LFu0IgnYQgjRykxWG6nDRpI6bCQet4sD27aQvXYVu9etZteqZeiNRhL7DqDb4OEkDxyCxW5v75KFEGdgzpw5/Pa3v23SdvXVV/PXv/61cQy2UgqTycTrr78O+ML2U089hdFoRKfT8dJLLxEWJgtbdRadZgy2pmmTgckpKSnTs7Oz27scIYQ4La/XQ/6OLHatXUH22lVUlxSj0+uJ79mH7kNGkJwxBL+g4PYuU4gOr6NM0yc6tzMZg91pAvYR8pCjEOJcpJTi0J5dZK9ZSfaalZQfLkDTfNP/dRsynG5DhuMfIr1bQpyIBGzRFuQhRyGEOMdomkZ0SirRKamMuvFWinNz2LVmJdlrVrDo7VdZ9ParxHRPp/vQkXQbMpyAsPD2LlkIIcRJSMAWQogORtM0whOSCE9IYsQvb6Ik7wDZq1ewa80KFs9+jcWzXyO6Wyrdh4yg+9CRBIRHtHfJQgghjiEBWwghOrjQ2HhCr76eoVdfT2l+HtlrVrBr9QqWvPsmS959k6jkbnRrCNuyiqQQQrQ/CdhCCHEOCYmJZciVv2TIlb+k/FABu9asYNfq5Sx7/22Wvf82EUnJpA4bReqwkQRGSNgWQoj2IAFbCCHOUUFR0Qy+/BoGX34NFYWH2LVmJbtWLWsM21HJ3UgdNoruw0YSECbDSIQQoq3o2rsAIYQQP19gRBSDJl/FTU/+izuef51RN96K1+tlybtv8tq903j/sUfY8OWnVJUWt3epQnRa8+fPR9M0duzYAUBOTg5Wq5V+/frRt29fhg8fzs6dO5tck5ubi91u5+mnn26PkkUrkYAthBCdTGBEFIMvv4Zbnvo30559hZHXT8HtcLBo1mu8es9tzJ3xWzZ+/Tk15WXtXaoQncqcOXMYOXIkc+bMaWxLTk5m06ZNbN68malTp/Lkk082ueY3v/kNl1xySVuXKlqZDBERQohOLDj66JjtkrwD7Fq1nJ2rlrHwrVdY+ParxKf3InX4KLoNGYEtILC9yxXinFVdXc3y5ctZtGgRkydP5vHHHz/unMrKSoKDjy4eNX/+fJKSkvDz82vLUkUbkIAthBDnidDYeIZdcwPDrrmB4gP72dkQtr9//SV+ePNlEvr0J234aFIGDcVsk//gC3EmPv30UyZOnEj37t0JDQ1l/fr1hIaGNi6VXlVVRW1tLWvWrAF8gfxvf/sb3333nQwP6YQkYAshxHkoLD6BsPgEhl97I0X797Fj5VJ2rlzK1y/9C73RSFK/DNJGjKbrgEEYzZb2LleIZlv09qsU7t/boveMSOjKBbfeecpz5syZw4MPPgjA9ddfz5w5c7jvvvsah4gAfPDBB9x55518/fXXzJw5k4ceegi73d6itYqOQQK2EEKcxzRNIyKxKxGJXRl1w1QKsneyY+USdq1azu4fV2E0W0jOGELaiDEk9u2P3mBs75KF6HBKS0tZuHAhW7duRdM0PB4PmqZx7733Njnvsssu47bbbgNgzZo1fPTRR/zP//wP5eXl6HQ6LBYL9913X3t8BdHCJGALIYQAfGE7pnsaMd3TGDvlDg5mbmPHyqVkr17BjhVLsPjZSRk8nLQRo4nv2RudTt/eJQtxnNP1NLeGjz76iFtuuYVXXnmlsW3MmDEcOHCgyXnLly8nOTkZgGXLljW2z5w5E7vdLuG6E5GALYQQ4jg6nZ4uvfrSpVdfxk+7i/1bNvmGkaxaxrZF32ILDCJ12CjSRowmulsamqa1d8lCtJs5c+bw29/+tknb1VdfzV//+tfGMdhKKUwmE6+//no7VSnakqaUau8aWlRGRoZat25de5chhBCdkqveyb4NP7JjxVL2bvwRj8tFYEQkqcNHkz5iDGFdEtu7RHEeysrKIj09vb3LEJ3cif6caZq2XimV8dNzpQe7hZTk7GLB458QFO3g8idntHc5nd6+5QtZ+vpmug0xMfzue09/gRCiRRhNZroPHUn3oSNx1taw+8fV7FixhB8/m8fa+f8lLD6BtBFjSBsxWpZqF0Kct6QHuwXkrFvB9y/swWmJA8DuXIvB7EQD0BSapkADrWFDpxpfa5qGplMNe9B0mu+1XkOn87Xp9Dpfu05reK1Dr9eh0+vRGfTo9Dr0Bj2aQY/eYEBnMGAwGNEZDeiNRgxGI3qDEb3RhMFswWg2+16bzBgtVgwWKzqjGZ3ZhM7Y8X/n2vLxXFYvMOEyBaF5XXTrupmLfvc/7V2WEOe12opydq5ezo7lS8jflQVAdPc00keMofvQkfgFBZ/mDkKcPenBFm3hTHqwJWC3gK+ffoqcHb04bFxIbG0Y9eaBgA6laaCde4tlal4P4EVTRzYPoH6y9x0DLxpeUF7A43t9pO2Y15rWcK3W0Nb4i0fDa53vvU4H6Hx7nU75fsnQ0/ALhYbb5aWgdCA6r4O4/jvJXxeJ05JAjN8KEvp2wWi2YDQbG36RsGL288Ns88NgsaK32tDbrOitVnQWEzrdufe/jRDngorCw+xYuZQdK5ZQnJuDpuno0rsvaSPG0G3wMJljW7Q4CdiiLUjAbuOA/fGGg/z7/S/o0S+DdfvLKKpyAgrwgteDSXNjUB4MuDDjxqx5sWgKs86DRacwah4seDDpFEalMGlejHgw4DtmUF4MGuhQ6JVCjxedUhiUQkOhUwqdAk150QGaV6EpfJvX97+vpgCvBkqhKd8epYHSUEpDU6Aa3h9pO3aP0qHw7VEaCh2g871veN3Y1vhe73utHfNe06E0X7vS9A2vtcbXqhmzEpgcBQROUPx5Vzi399RhX7ALh/XM/2LVvB405UZTnmM2Nxqehl8iGo7h29AaXmu+XyR0mhdN50XTvOj0Cp1OodeDTg96PeiNOgwGHXqTAaNJj9FkwmQxYbRaMFutWGw2jDY/jHY7Bj8rBj87xoBADHYberNBHhoTnUJxbk5j2K4oPIzeaKRr/0GkjRhN0oBBGE3m9i5RdAJZWVmkpcnDtqL1KKXYsWOHBOy25PZ4+XRTPlcNiKWoysny3cW4PF5cHoW7Ye/yenF71E/avbi8vtf17p+8bjy36Wu3VzUc951X7/Hd1+1t2f8dTQYdZr0Os1GHSa/DbNQ37HWYDTrfcYMes0HXsOkbjzW2G3VYjPombY3vjXosjfcCg96LQedF8zrxOmtxOeqor6vB7XDgctTiqqvD7azH7ahjsyWNpxaXEGIzUVJTz60DAxh0aDnuuno8bi9elwevW/k2j0J5FcoDyqPh9fh+aVBeDbw6lFd39BcIpQfV8EuB0qMwAPqGzffa9wuBwbfXDKjGTY9XZ2jRf7HQvC50XjeacqFT7oZfBlxouNHhPrrXPOh1HnQ6LzqdF70eDAYwGMFg1GMy6TFZTJgtJix+Vsw2Gwa7FaPdjjHA37cFBmEK8MNgM6PTyX+gRMtTSjXOsb1z5TJqK8oxWW10Gzyc9JFjie8l0/6Js7dv3z78/f0JDQ2VkC1anFKKkpISqqqqSEpKanJMAnYn5/H6Qnh9Q/B2e1RDAPdQ71bUe5qG8nr3MZvHi9PlweVRON0e6t1enA3nON3H7j2N7x0uT8N1vusdLk9Duxen28PPyftGvYbFoMdi8oVwi0GP1aTHYtDjUYr1+8uY2DOKf13Xj398s5M3V+xjYs8oBiYEYzHpsRl951uNeixGPTaTbzv62oDFqPtZfwkrpXB5XdR76qn31lPvdlJfX4OjpoL6qkqc1ZU4q2uor63F5XDgqnXidjhxO524HR68LjeeeoWn3otyK7xuUG5QHh3Ko6G8ejiyHQn+yujb0ANG0AyAybfXDHg1E0pnwKszonRnvxiIL9i70Hnr0al6dMqFDt+m19zodB4MOg8GvRejUcNk1DCa9FisRsxWM0Y/CyY/K8YAP0yBAZiDgzAGBWMO8cfoZ8ZglBB1vvN6PORu38KO5UvIXruC+ro6/IKCG2ciiUzuJiFJnBGXy8XBgwdxOBztXYropCwWC3FxcRiNTf/7KgFbtCmXxxfKnS4PjiP7hvB9/N73unF/TJvveg+19b62OpeX0d3C+PWF3dHrNJRSvLBwN8/+kI3nDFK9poH1mMB9JIQ3eW024NfQ5mf+yd5kwGbWYzf7zrebDfiZDRj1bTeu26u81HvqcXqcONwOHB4HDlcdTkcFjppy6srLcFRW4KyuwVVdR31tHe5aFx6HC7fTg9ep8NYDLl/PPm49mscAXt+mKSOaMvn2GNEwgebbfGHehFdnwqs3nVHdmteNzluP3utEr+rR49sMOjdGnQejQWEygsWsx2g1YLKZMPlZMAX4YQ4KwBwagjkkGHNoAOZAG0azXsLYOezItH9Zy5ewb+OPeNxugqNjSBsxhvSRYwmOjm3vEoUQ4qQkYItOze3x4nB7qav3+DaXh9p69zGvPT957aa23kPNMa/rXB5qnL7Xvu3o6+YyGXQNYVuPn8mA3WzAbmnYm0/w3mLA32LEbjYQYDn63s/UcUKjy+uizl1HrauWOncdde46auprqKuvpLaulLrKEpxlFTgqKnFVOnDVOPHUuvDWefHWg3Jq4NKhuY1obgM6rxm914jea0GHCZ0yA2bQzKAz49WZ8eoseAyW5hWovOg9DgzKiV45MVKPSefGpPdgMiqMZg2T1YDZZsIcYMMcHIAlNARLeDCW0EAsoXbMNiP6NvzlSJyYo7qa7LUryVq+mAOZW0EpIrt2I33kWFKHj8IeHNLeJQohRBMSsIU4S16vorYhsNc6PVQ3hPAap5uaejc1TjfVTg+1TjfVDe9rnB6qHO7Gc6odbqqd7sZrT0fTaAjdRvwtR/f+FgMBVuPR9obXAVYDgY2vjQRYDBg6cGCs99RT7aqmxlVDjauG6vrqhn0VNY5SqmpKcFQU4yitwFVZi7uqHm+NF+oUOPRo9Qb0bhMmtxmjx4LRa0GvfJuGBbCCztIQ1K149Kd/kE7ncWL0OjDixKTVY9S7MRu9mMw6LDYDJn8LlmB/rGEhWCOCsUUEYg0LxOJnxGCSYS8traq0mJ0rlpK1YgmF+/agaTrie/UhfcQYug0ZLjORCCE6hE4fsDVNmwxMTklJmZ6dnd3e5QhxUh6votrppsrhosbpodrposrhpqohhFc5XFQ73FQ2tFU5XFQ6fOcc2Vc53KcdEnOkVzzAaiSwYQuyHX0daDMR1NAWZDX59jZfb3pH6T0/Fa/yNobzyvpKquqrqKyv9G3OSirrSqisK6G6qghneRmecieq0gM1Gvo6PZZ6M1aXDYvLgsljweSxofda0SsrOs0GmhWvzobbYMVtsJ7yAVad14XRW4cJJyZdPWajF7NVw+JnxBrkhzU0CGt4ELaIQGyRwVgDLFj8DOg68C9BHUnJwQPsWLmErOWLqTh8yDcTyYBBpI8cS1L/QRiMZ//MgRBC/BydPmAfIT3Y4nyglKK23kOlw0Vlnbth76KiYauscze+9r13NXlf5zp5L7pBpzWE7SMB3ESIn5Fgm4lgPxPBNt/rED/f+xCbiUCr8ZyafUQpRY2rhor6CiqcR7dyZ7lvc5RRVltERV0xVVUl1JdXoipdWGtM+Dus2B027E4bVpcNs9u3mZQNnfIDzYZX54fb4IfL6HfKB06NnlpMOLDoXZjNXqw2A5ZAK7bQQGzhAfhFBGKLCsYWZMXmb0JvPL8DuVKKQ7t3kbV8MTtX+WYiMfv5kTp0FD1GjyMmNf2c+OVQCNF5SMAWQjRyuj2+sF3rorzORXmti7Laeioa9mW1Lirq6imrOfLe11bv9p7wfjqNxgAe4mci9Jh9qN1MqN1EqJ+ZMLvvfdA5FsiPcHqclDnKKHWUUuoobfK61FFKWV0ppbWFlDhKqKwuw1TtIqjWRHCNjcA6OwEOP2z1fgTU+2F1+2Hy2NCUHaXzx62z4zL5+0K5duIhJ0ZvHWbNgcXowWrTYQu0YguxYwv3xx4djD0qGFugGb9Ac6cP416Ph9xtm8lctojstStxO50ERUbTY/Q40kddQFCkLNMuhGh9ErCFED+LUoo6l4fSGl/wLq2tp6ym3ve+tp6SmnpKq33vS2qclNbUU17n4kR/xeg0CGkI3GH2Y/b+ZsLsZsL9zYTbzUQEmAmxmc7JMK6UotpVTUldCSWOEorrihtfl9T53hfVFVFUc5gSRxlKefCrg8AajZAaCyE1fkTWBhDi8MPutGN2+aHzBqA0f1y6AFzGAOqN/riNthN+vlHVYdXXY7XgC90hftgj/PGPDcEeFYg9yIIt0NQppk2sd9SRvWYlmUsXkrt9CyhFbFpPeoweR+qwkTJeWwjRaiRgCyHanNvjpbzORUl1PSXVToprfPuS6nqKq50UN+59m8N1fA+5XqcR6mci3N9MhL8vfEcGWIjwNxPRsI8MsBDub27TaRJbksfrodRRSmFdIcW1xRTWFVJUW0RRXRGFtYUU1hZyuOYwZc4y3wVKYa2HoGqIrzGRUB1IRLUfgTU2LE4belcAHhWESxeE0xiI0xRIvSkApTMc99lGnFgNLvzsevyCrQRE2gmIC8U/Jgj/EAv2YAtG87kTwiuLi8hatojMpQspzT+I3mgkJWMoPcaMI7HPAHT6c+e7CCE6PgnYQogOTSlFTb2Hoipf2C6q+snW0Ha40kFxtfOEixmF+pmIDLAQGWAmKtBCZICFqAALkYEWIv0tRAdaCLIZz9lxug63g8O1hzlcc5hDtYc4VHOo8fWRfYWzosk1fm4dqbUWulbpiCszE1ppxlptRe/0x+MJwkkQdfoQnOYgnOZg6k0Bx32ukXpsFi9+dj32ECv+0YEExoXgH+5HQJgVe7C5wz2wqZTi8J5sti9dyI6VS3FUVWILDCJ95Bh6jB5PeELSOfvnQAjRcUjAFkJ0Gh6voqTGSWGlL3AXNgTvw5VH9o6GIF5/3LUWo47oQCvRgRaiAi3EBFp9+yAL0YFWYoOtBFjO3Vkpal21FNQUkFedR351PvnV+Udf1+RT6ihtcr5ZZyRW7093h4mUUjfxRU4CSwwYKy146kOoUaHUaqHUGXwh3GEOxmXyb3IPDS9Woxt/fx0BYTYC44IJ7BJCQKiVgDALfoFmtHYc5uNxu9i3cT2ZSxeyZ/1avB43YV0SfeO1R4zBHhLabrUJIc5tErCFEOedereXwipf2D5U4aSgoo5DFQ4KKhwUVNRRUOE79tPecH+zgdhgKzFBVmKCLMQG2YgJshAXbCUu2Ea43XxOjguHEwfwA1UHyK3K5UDVAercdY3nGjQDcZZQEoz+JLh0JJc4iCusJORQFbpSPVX14VS6w6kmghpDOA5LKHWWUOrNQU0+U4cXm8WLf5CRwEg7gfEhvn2ElaAIGybr8UNXWktdVSU7Vy0nc+kPFGTvRNN0JPTpR4/R40gZNBSjuZkLHAkhBBKwhRDihNweL0XVTvLLfaE7v7yOvLI68sodvtfldVTUuZpcY9LriA22NgbuuIbX8SG+1+F28zk5/EApRVFdEbmVueRW5R63bxK+dQbibdEkm0NI1llIqa0nMb+EyAMFuA9VUFEdQoU7gipPJNVaBHXmMOosITgsocf1gFuMHgJDDATHBhKcGEpwpI2gSBsBYVb0htYbelKan0fWsoVsX7qQquIijBYr3YeOoOeY8cSl9zon/zcUQrQtCdhCCHGWqp1uCsrrOFhex8GyOg6W1XKwtGFfVkdJTdOhKFajni4hNuJDbCSE+rYuIb4tLtiGqRVDY2tRSlFcV8z+yv0cqDrA/sr95FTmsKd8DweqDuBRvrnVdZqOeHscyfZYkg0BdPXqSHbUEl9wCF3uflx5h6ipMlHhiKLcFU0lMdSYI6m1RVBrjWgSvjUUdj8IjLASkhBCcLSdwEgbQRE27EEtN+xEeb0c3LGdzKUL2bV6OfV1dQRHx9Bz7EX0HDNelmgXQpyUBGwhhGglNU43eeW+wH2gtI7c0lr2l9SSW1pDbmltk9lRdBpEB1pJCLWRGOZHUqgfSWF+JIb50SXk3Azf9Z56cipz2Fu+lz0Ve9hT7ttyK3NxKzcAGhpx/nEkB3alqzWSZJ2F1HoXSRWFaDlZ1O/dTX1hLVXVAZQ5Yyn3xlKti2oM3rW2CLzHLHmv1ykCg/UExQQQnhRMaKyd0Fg7/qGWn9Xz7HI42LVmBVsXfkveju1oOh1J/TPofcEEkvpnoDe03XAWIUTHJwFbCCHagVKKoion+0tryS2pZX9pLQdKa8kpqSGnuIay2qPDT3QaxAX7gnfXMD8SG0J4crid2CDrOTfu2+VxkVuV2xi4j4TvnMoc3F5f8DbpTHQL7kZ6aDrpfnGkKxPd6qqxFGfjzcukfs8unEW1OCsNVNSEUe7pQoU3klqLL3TX2KJwWMMbP9NoUIREWQlLCiEs3t8XvGP8zmqcd2l+HtsXf8f2JT9QU16GLTCIHqPH0XvcBEJi4lrs5ySEOHdJwBZCiA6ovLaefcU15JTUsK+ohn0lteQU17CvuIZqp7vxPLNBR9dwOykRdpLD/UiJ8L1ODPXDco4tFuPyujhQeYAdpTvYUbqDzNJMskqyqKyvBECv6UkKTKJHaA/SQ9JJt0WT6nJjL90PRTtQh7Nw7dtJfWE1zgoD1RV2it3JlHljqfbzbTX2WNwGa+Nn2gN0hCUEERYfQGisnbA4OwHhzfulxevxsG/TOrYu/I69G9aivF5iUnvQ+4KL6D5sJCaL9bT3EEJ0ThKwhRDiHKKUorjaF773FlWzu7CaPUXV7C6q5mBZXeMKmToN4kNsJB8TvrtH+tM90h8/87kznEEpRUFNAVklWY2Be0fpDorqihrPSQhI8AXu0HTSgtNIt0QQXH4ACjZB/ia8uRtw5hbgLDNSV26gvK4LJa5YqswxVNtjqLbHUWuNAM03DEevh5BoG2EJgb7Q3TDMxGI/+TSNNeVlbF/yA9sWfUdZQR5Gi5W04aPodcEEorulyoORQpxnJGALIUQn4XB52FtUw+6iavYUVjfu9xbXUO8+Ot47PsRKamQAqVF2ukf6kxYVQFKY3zk1zruotoisUl/YzirJIqs0i7zqvMbjMX4x9I3oy8CIgfSP7E+KKQTdoS2QvwnyN6LyN+I6mI+j3IijzEhtbRgldbFUaBFU+8VQbY+l2j8el+Hocur+QUaiu4cQ1TWQ6JRAQmLsx/V0K6XI25nJtoXfsXP1MtxOJ6FxXeh1wUX0GD0OW0BgW/2IhBDtSAK2EEJ0ch6v4kBpLTsPV7HrUBU7D1ex81AVe4tr8DRM9m3QaSSH2+ke5U9qpJ3UqADSovyJC7aeM72vFc6KxsC9rWQbGw9vpLCuEAB/kz/9wvsxIHIA/SP60yusF2ZHdUMv98bG3m53YR7OMgOOciN1tSFUVYdS4gyl2hZLZUAiFcEp1Bt8M5oYjRqRyYFEJQcR3TWQyKQAzLajvdzO2lp2rlrKtoXfUbB7Jzq9geSMwfS+YAIJffuj051bQ3iEEM0nAVsIIc5TTrevx3tXQ+De2RC+D5Ydndc6wGKgR0wAvWIC6RkbQM+YQLqG+WHoYEugn4hSirzqPDYWbmRD4QY2Ht7Inoo9ABh1RnqF9aJ/RH8GRAygX0Q/As2BUFPs6+Uu2NjQ270Jb9lBnBVG6kos1NZEUFoSTIkWQ0VgVyqDUqiyRTcOLwmJshLVLZjoroFEdQ0kMML3C0rxgf1sW/QtmUsXUVdViT00jF5jxtNz7EUERUa13w9JCNEqJGALIYRoosrhIruwmqyCSrbn+7YdBZU4G4aZmA060qID6Bnj23rFBJIa5X9OPFRZ7ihnY+HGxtC9vWR748wlKUEpDIgYQP9IX+iO9ov29d5XF/l6uXNXwr5lqLyNuKqhrsRKbW0MVcU2imqDqQhIojIwmYqgZNw639SBFpueqJRgopN9gTs01sL+LevYuug7cjZvAKXo0qsPvS6YQMrgYRhN5lOVL4Q4R0jAFkIIcVpuj5e9xTVsy6toCN2+fZXDF071Oo2UcLsvcMcG0q9LED2iAzp86Ha4HWwr3saGwg1sKNzA5sLNVLuqAYi0RTIgYkDjsJJuwd3QaTpwVkHuGshZCjnLIX8THqeXulIbtY4u1BabKS4yUWHrQkVAVypDu1NjCgV8Hd3hXfyJ7hpEYLiH0rx17Fq9iMqiw5j9/EgfeQG9x00gIrFre/5YhBA/kwRsIYQQZ0UpxcGyusawvT2/km15FRRWOQEw6jV6RAfQNz6Ifg1bUphfhx7T7fF62F2+2xe4D/u2I+O4QywhjIodxei40QyPGY7dZPdd5KiE3NWQs8wXuAs2odxeHJU2ap1dqS3zoyLXSZkuoiFwd6PS3gUPvtlc/IJMBAQX46jazOF9G/G63UR27UbvcRNIGzEGs83WXj8OIcRZkoAthBCiRR2udLAxt5xNB8rZdKCMrQcrqKn3LZkeaDU2Bu7+8UH0jQ8ixM/UzhWfnFKK/Jp81h9ez/K85azIW0FlfSUGzcDAyIGMihvFmLgxJAYmHr3IUXE0cO9bBoe2oLxe6mut1Lm7U1sRSM3+akorTVQEdKUiOJnKkO44dHaUtw6bfR/1tVuoKcvHYDaTOnQUvcdNICY1vUP/ciKEOEoCthBCiFbl8SqyC6vYfMAXujfmlrPrcBUNE5iQEGpr7OHuFx9Ez5jADjtloNvrZnPRZpYeXMrSg0vZXb4bgC7+XRgdN5rRcaPJiMzAqD9mzuy68mN6uJdBwRZA4XZZqfWmUVcVRs2BWkry6ikO6UVx9EAqLDEoz2F0Wib1tVl43U6Co2PpPf5ieo4ehy0wqD2+vhCimSRgCyGEaHM1Tjdb8yp8vdwNvd2HKh0AWIw6BnQJZnBSCEOSQunfJajDjuXOq85j2cFlLDm4hLUFa6n31mMz2BgeM5zRcaMZFTeKMGtY04vqymH/St9wkpxlcGgroHCrQKrd/anOM1OSeZgiv+4UR/anNDABV/0elHsbnvp8NJ2ergMG0/eiiST06SfT/QnRAUnAFkII0SEcqnCwMbeMtTmlrN1XSmZBJUqBSa+jb3xgY+AemBDcIVejrHXVsvbQWpYeXMqSg0sorPWN3e4Z2pPRcaMZEzeG9NB034OSTS4s9QXunV/CjgXgqMCrD6BGP4TqQ/6UbcqlSEVSFNaXotA4nPW78NRngqrDYg+m59gLGTDxEgLCI9rhWwshTkQCthBCiA6pos7F+v2lrNlbypp9pWzNq8DjVeh1Gr1iAxmSFMKQpBAyEkMItJ58GfP2oJRiV9muxrC9pWgLCkWoJbRx3PawmGH4Gf2aXuiuh31LIHM+ZC0ARznK6I/DbyRVxWFUbjpAYbGiKLQXBWFh1Ln34HXvByA4Jo2+F02kz4WjMZo67rh2Ic4HErCFEEKcE2qcbjbklrFmr6+He9OBcuo9XjQN0qMCGJwUwtCuIQxKDCHU3rHmky51lLIibwVLDy5lRd4KqlxVGHQGMiIzGnu3uwR0aXqRx+UL29vn+3q268rAZKc+bCzV5XFUbsunKOsQ+aGp5If4UevNBVWFprcSlTyE/hMvIXVYj+OWcxdCtD4J2EIIIc5JDpeHTQfKfYE7p4T1+8twuHyL4XSPtDOqWzhjU8MZnBSC2dBxxim7vC42FW5qHLu9t2IvAN2Cu3FVylVM6jqJIEtQ04s8Lti39GjPdl0pmOx44sdTXZNMVVYJxWu2sTcwjoIgEw4KAS96UwwxqSPod9E4kvrFYDR3nJ+DEJ2ZBGwhhBCdQr3by9a8CtbsK2HVnhLW7Cul3u3FatQzLDmUsanhjOkeTkKo3+lv1oYOVh1kycElLNizgG0l2zDqjIzvMp6rul3FkOghx4/Z9rh8D0ce6dmuLQGjH6rrRdR4e1KdXcOh5WvYafDjUIAOt1YHGDFY0ohNG0H6yP4k9gnDL7Bj9fIL0ZlIwBZCCNEp1da7WbO3lMU7C1m8q4j9JbUAJIX5MaZ7OGNSwxnWNbRDzVCys3Qnn+z+hM/3fE5lfSWx9lguT7mcK1OuJMov6vgLPG5f2M78FLI+h9piMNpQKRfhtA2hco+TnGXL2eFwUeynoTQvmj4MvakXUd0G0y0jgbRh0dgCZMy2EC1JArYQQojzQk5xDYt3FrJkVxGr9pbgcHkxG3QM6RrK2O6+4SQdZaVJp8fJwtyFzMuex5qCNWhoDI8dzlUpV3FB/AVN59k+wuOG/SsahpF8DjVFYLRBt4twRYylZI+HrCWLyK4sp8qsA3TojN0w2frTY0wG/S/qQlCErBopREvo9AFb07TJwOSUlJTp2dnZ7V2OEEKIDsDh8rBmXylLdhaxeFche4tqAIgPsTK2ewRjU8MZlhyKzdT+0wEerDrI/N3z+WT3JxTWFhJiCWFS10lc1e0qkoOST3yR1+ML29vnQ9ZnvrBtsEK3i/AkTSR3h4OtS5ewr6IEt05Dp49EbxlAt0HDGfiLZCISAtr0OwrR2XT6gH2E9GALIYQ4mQOltSzeVcSSnYWs3FNCbb0Hk17H4KQQxjT0bqdE2Nu1d9vj9bAyfyUfZ3/M4gOLcSs3fcP7cnW3q7k48WJsxpP0Pns9vnm2M+dD5mdQUwgGC6RciCPuItYv2sOWLZuo1Sk0zQ+9pT+xqcMZfFlP4nuEdIgefSHONRKwhRBCiGM43R7W5ZSxZFcRi3cWsutwNQCJoTYm9YlhUt9oUiP92zV4ltSV8Pmez/l498fsq9iHzWDjkqRLuLLblfQJ63Py2rweyF11tGe7+jD4x+Adeh+ZO3Ss//47ir31gAG9uSeh0YMZfNUQumVEoNN3zOXrheiIJGALIYQQp5BfXseinYV8tfUQK/cU41WQEmFnUp9oJvWJISXC3m61KaXYXLSZednz+CbnG+rcdaQEpXBlypVMTp5MsCX45Bd7PbBnESx/xjecxBYKQ+/moDOdtXM/JKeqFKWBzpCEX9BABl85jp6jYmWqPyGaQQK2EEII0UzF1U6+2naIBZvzWZtTilKQFuXP5L4xTOoT3a5TAFbXV/N1ztd8kv0JW4q3YNAZGBc/jqu6XcXQ6KHodacIxvtX+YJ29rdgDoBBd1ARdSlrZ71LVs5eXDrf7CNmW3/6T7yI/henYPWXmUeEOBkJ2EIIIcRZOFzp4MutBSzYUsD6/WUA9I4NZFKfaC7tE01ccPvNyJFdls3H2R+zYO8Cyp3lRPtFc0XKFVyRcgUx9piTX1iwGZY945v2z2CGAVOp7zONTXM/ZcPa1dRoLtCsGM19SBt2AUOuHkhguLXtvpgQ5wgJ2EIIIcTPlFdex5dbCliwJZ/NBysA6Bcf1Bi2owPbJ4TWe+pZeGAhn2R/wqr8VQAMixnGld2uZFz8OEz6k/RCF2fD8mdhy1zf+77X4x32ILsXr2fNZ59SWF8Nmg6dKY2E7sMYNfVCwrv4t82XEuIcIAFbCCGEaEG5JbUs2JrPgs0FZBZUAjAoMZhJfWK4pHcUEf6Wdqkrvzqf+bvnM3/3fApqCggyBzE5eTJTe0wl0i/yxBeVH4CVz8GG2eB2Qs8rYNTDHD5Qy8rXX2df2WGU5kVniCcscgCjp11Ol57hMvOIOO9JwBZCCCFayd6iahY09GzvOlyNToMhSaFM6hvNJb2iCfFr+3HMHq+H1QWr+Tj7YxYeWIhBMzC151Ru63UbfsaTjCGvLoTVL8Ha16G+CrpdDKMepsYQz+oXnmP77p24NBeaLgj/gD4Mv+4K0sekyMwj4rwlAVsIIYRoA7sOV7Fgcz4LthSwt7gGvU5jeHIok/vEcHHPKAJtJ1idsZUdrDrIcxue46ucrwi1hHJv/3u5MuVKDLqTLLBTV+YL2atfgrpSSBgJox/GHT2UzW++zo/LV1CjakAzY7GkM3DCRAZeMwyjSWYeEecXCdhCCCFEG1JKkVlQ2dizfaC0DqNeY3xaJFOGJzCsa2ibD7HYUrSFf677JxsKN5AcmMxvMn7DqNhRJ6+jvgbWz/INH6kqgJj+MOphSL2UvV8uYPmHH1PkLAHAYEwmte9wxtx9BVa7zDwizg8SsIUQQoh2opRia14Fn2/O56P1BymrdZEa6c+U4Qlc2T+2TZdqV0qxMHchz6x/htyqXIZED+HhgQ+THpp+8ovcTtg8x/dAZNk+CE+Dkb+BXldTuiOLxS++Rk5JLgo3On0kcfH9ufCBmwiOPcX83EJ0AhKwhRBCiA7A4fLw2eZ8Zq3MYXt+JQEWA7/MiGfKsES6hLbdlH8uj4sPd33Iy5tfpsJZweTkydzf/36i/KJOfpHH7VuKfdk/oTATghJgxIPQ7yYclTUsf/ZFMndtw0UtmmYnNLgHF0y/kS4DUtrsewnRliRgCyGEEB2IUop1+8t4e2UOX287hFcpxqdFMHV4IiNTwtps+EhlfSWvb32d9zLfQ9M0pvSYwrRe07CbTrFypdcL2d/A0qchbx3YI2HYfZAxDa/ewobX3uLHZcuo9ZQCBuy2FIZdPpnel4+WmUdEpyIBWwghhOigDlU4eG/NfuaszaW4up7kcD+mDk/kqgFx2M1tM3wkrzqP5zY8x5f7viTEEsI9fe/hqu5XYdSd4qFMpWDfUl+P9r4lYAmCoXfD4DvBFsLeb75nyZx5lNblAV5Mxjj6DB7BiHtuxGCQByLFuU8CthBCCNHBOd0evthSwKyVOWw+WIG/2cDVA+OYOjyRpLC2WZ59e/F2/rHuH6w/vJ6kwCQeGvAQY+PHnr7n+cCPvmXYd34JJjtk3Obr1faPonRXNt8/9zoHi/agcKDXBdM1qQ8XPjIdW0hQm3wvIVqDBGwhhBDiHLIxt4xZK3P4YmsBLo9iTPdwbh2RyJhu4eh0rTvMQinF4gOLeWb9M+RU5pARmcEjGY/QM6zn6S8+vB2W/wu2zQOdEfrfDCMegOBEHBWVLP7nK+zctQm3qkDTrPRIH8XEGQ+06vcRorVIwBZCCCHOQYVVDuasOcC7a/ZTVOUkKcyPW4YmcE1GHAGW1p1T2+V1MW/XPF7a9BJlzjIu7XopD/R/gBh7zOkvLtkDK/4Nm94H5YXe18Ko30B4Kh6Phx9fm8PaJd/j8hYT5JfCDf98HFtwYKt+HyFamgRsIYQQ4hxW7/by1Tbf8JENueXYTHquHhDH1OEJpET4t+pnV9VX8ea2N3kn8x2UUtzc42bu6H0H/qZmfG5FHqx6Eda/Ba46SJ/km0s7pj/O6ho++M1MiiqyMOhCuOyu+0gaM7hVv4sQLUkCthBCCNFJbD1Ywdsrc/h8cz71Hi8jU8KYOjyRcWkR6Ftx+MihmkM8v/F5PtvzGUHmIO7qexe/TP3lqR+EPKKmGNa8DGteBWcFJI+DUY9A4ggWPf0aG378GlD0638R4393d6t9ByFakgRsIYQQopMpqXYy98cDvLNqP4cqHcSHWLllaALXZXRp1SXZM0sy+ee6f7L20FoSAhJ4aMBDjOsyrnlT8Dkq4cfXfb3atcXQ/xa49J/krN7Gpy8+h9tbQph/Gtc/+zhme9s82CnE2ZKALYQQQnRSLo+Xb7cfZtbKHNbmlGIx6riyv2/4SFpUQKt8plKKZXnL+Oe6f7K3Yi8DIgbwcMbD9Anv07wb1NfCsqd9U/zFDYbr3qXOY2PuwzMord6FURfGlb/+NfFD+rVK/UK0BAnYQgghxHkgM7+SWStzmL8pD6fby9CuIdw6PJEJPaJaZfYRt9fNx9kf8+KmFyl1lHJJ4iU8MOAB4vzjmneD7fNh/t2+ObSvfxdiB/LdX15gy5Yf0DQ9GUN/wehfT2vxuoVoCRKwhRBCiPNIWU09H6zzDR/JK6+jT1wgMyb3YGBCSKt8Xo2rhre2vcWs7bPwKA83pt3I9D7TCTQ3Y2aQQ1th7o1QdRguex76Xsfu75az4M1X8HjLiAzqyXX/ehyjzdIqtQtxtiRgCyGEEOchj1fx6aY8/vb1Dg5XOrm8Xwy/uySN6EBrq3ze4ZrDvLDpBT7d/Sn+Jn/u6nsX16dej1F/mjHhNSXw36mQswyG3w8XPk5taSXvPzKDirq9mPQRXPPoI0T379EqdQtxNiRgCyGEEOexGqeb/yzew6vL9qLXNO4Zm8z00V2xGFtnyfKdpTt5et3TrC5YTbx/PL8e8GsuSrjo1A9CelzwzR9g7auQPB6ueQOswXz5p2fI2rkUTTMybOyVDLvrxlapWYgzJQFbCCGEEBworeXJL7P4atshYoOs/O+l6VzSK6p5M4CcIaUUK/JX8M91/2R3+W76hvflkYxH6BfR79QXrp8FXzwMQfFww1wIT2XHgoV8/e4beFQlMeF9+OUzM9GbWnehHSFORwK2EEIIIRqt3FPM/32eyY5DVQxJCuFPk3vQM6Z1VlL0eD18uudTnt/4PMV1xVyUcBEPDXiI+ID4k1+Uuxo+uMW3OM3Vr0HqJVTlF/L+7x6n2rkfiyGKXz72B8LTurZKzUI0hwRsIYQQQjTh9niZ++MB/vntTirqXFw3qAuPTOhOqN3cKp9X66pl1vZZvLX9LVxeF9enXs99/e/Dz3iS+a4rDsLcm6BgM4z7I4x6GI/Xyxf/+w+y961C08yMnng9Gbde1Sr1CnE6ErCFEEIIcUIVtS7+/UM2s1flYDXpeXB8N6YMS8Rk0LXK5xXVFvHiphf5ZPcnpAan8p8L/0OoNfTEJ7vq4LMHYOuH0OMKuOIlMPmx9b9f8f282XhVDV2iB3DVPx9Dr2+d8eRCnIwEbCGEEEKc0u7CKh7/PJNl2cV0DffjT5N6MDY1otU+b+nBpTy8+GEibBG8OuFVYu2xJz5RKVj5PHw/AyJ6wvXvQXACZfvzmPu/T1DrOojVGMsNj/8vwcldWq1eIX5KArYQQgghTkspxcIdhTyxIJOcklrGpUXwx0vT6Rpub5XP21i4kXt/uBeL3sIrF71Ct+BuJz85+3v4aBroDfDL2ZA4Eo/Hw/xHnyQn70d0mo1xV9xM3+sntUqtQvyUBGwhhBBCNFu928vbK/fx3A+7cbg83DYikfvHdyPA0vIzd+wq28Vd392Fw+PgpfEvnXqWkeLdMPcGKN0LE5+CQXeAprHhnU9Y/MVclKojOWEIk//6OxkyIlqdBGwhhBBCnLGiKidPf7OTD9cfIMRm4pGLU/llRjz6Fl52Pa86jzu/vZPC2kL+OfafjI4bffKTHRUwbzpkfwMDpsIvngaDieKd+/jg/57E4S7Aborn+r/OJDAuskXrFOJYErCFEEIIcda25VXw+Ofb+TGnjB7RAcyY3IMhXU/yYOJZKqkr4e7v7ya7LJsnRj7BpK6nGOrh9cCiv8Cyf0L8ULjuHbBH4Kl38dHDT3CwcCM6zZ8J199GzysuatE6hThCArYQQgghfhalFAu2FPDXL7PIr3BwaZ9ofn9JGnHBthb7jOr6ah5Y9AA/HvqR3w76LTf3uPnUF2ybB/PvBVuI7+HHmP4ArHltLit++BilnKSmjGLSXx5psRqFOEICthBCCCFaRF29h1eW7uHlJXtQCn41Jpm7xyRjNbXMmGenx8lvl/6WH3J/YHrv6dzf//5TrzRZsAXm3gg1RXDZC9DnWgAObd3JR3/9O07PYQIsSdzwj5nYI1q2112c3yRgCyGEEKJF5ZXX8dRXO/h8cz7RgRZ+/4t0JveJbpFl191eN0+sfoKPsz/m2u7X8r9D/he97hQBvqYYPpwC+1fAiAdh/AzQ6XE5nXz468c5VLoFvS6IX0y9k+4TTzG+W4gzIAFbCCGEEK1i7b5SHv98O9vzK8lICGbG5J70jvv5y64rpfj3hn/zxrY3uCjhIp4a9RQmvenkF3hc8NVvYd0bkHIhXP0GWIMAWP78bNau+Ayl3PTsMY6JMx742fUJIQFbCCGEEK3G41V8tP4A//hmJyU19Vw7MI5HLk4lwt/ys+89a/ssnl73NEOihvDvcf8++dLqR6x7E758FIIS4Ia5EN4dgLz1W5n3z2dweYoIsqVwwzOPYwv++b8IiPNXpw/YmqZNBianpKRMz87Obu9yhBBCiPNSpcPFCwt389aKfZgNeu4fl8KtIxIxG37e+OzP9nzGn1b8idQQ39LqIZaQU1+wfyV8cAt46uGq1yB1IgDO6ho++M1MiiqyMOhCuOyu+0gaM/hn1SbOX2cdsDVNa+6ao+VKqcqzKa4lSQ+2EEII0f72FlXzly+y+GFHIYmhNv54aQ/Gp0f8rPHZiw8s5pEljxDtF80rF71CjD3m1BeUH4APbvI9BDn+MRj5G2j4/EVPv8aGH78GFP36X8T439191nWJ89fPCdiLAAWc6v8RCnhbKTX7Z1XZAiRgCyGEEB3H4p2+Zdf3FNUwuns4f5qUTkqE/1nfb8PhDdz3w31YjVZeufAVUoJTTn1BfS18dj9s+wh6XgWXvwAm3xCTnOXr+fTF53B7Swj1T+WGZ2ZiDjj72sT5p9MPETlCArYQQgjRsbg8Xt5ZtZ9/fb+L2noPU4Yl8Ovx3Qm0nd2y6ztLd3LX93dR76nnpQtfom9431NfoBSs+Dd8PxOiesH1cyAoHoC68irmPjyD0updGHVhXPnAg8QP639WdYnzT4sFbE3T/ACHUsrTUsW1JAnYQgghRMdUUu3kme92MWdtLoFWI7+ZkMqNg7uc1bLrB6oO8KvvfkVxXTHPjH2GkbEjT3/Rrm9h3h2gN8IvZ0PiiMZD3/3lBbZs+QHQMWjYpYz+9bQzrkmcf37OEBEdcD1wEzAIqAfMQBHwBfCKUmp3i1d8liRgCyGEEB1bZn4l/7dgO6v3ljKhRyTP3dAfi/HMH4Isrivmru/uYk/5Hv4y8i/8ousvmnFRNsy5Hspy4JK/w6DbGw/t+X4ln7/xHzzeMiKDenLdvx7HaPv5s6CIzutkAVvXjGsXAcnA74EopVScUiocGAmsBv6madpp1jEVQgghhPDpERPAnOlD+dOkHnybeZipb66l0uE64/uEWcN4a+Jb9I3oy++W/Y73s95vxkXd4I4fIHkcfPEb+PzX4K4HIPnC4dz5/HMEWrtyuHw7L99xDwUbM8+4LiGa04NtVEqd8k99c85pK9KDLYQQQpw7Pt2Ux8MfbqZ7pD+zpg0m3N98xvdwepw8uuRRFh1YxF197+KevvecfrYSrwcWPgHL/wVdhsEv3wF7eOPhL//0DFk7l6JpRoaNvYJhd910xnWJzu+se7CVUi5N09I0TRuvaZr9JzedeOSclitVCCGEEOeLy/vF8vrUDPYV13DNyyvJLak943uY9WaeGfsMV6ZcycubX+Yva/6Cx3uaR8V0erhwpm+1x/xN8OpY377BL/7vN1x6ywPoMLFy0Vzm3Pu/eOol7ojmOW3A1jTtAeBT4H5gm6Zplx9z+MnWKkwIIYQQ54exqRG8N30IFXUurn55JZn5Z76shkFn4PHhj3Nbr9v4YOcH/HbZb6n31J/+wt7XwLSvfa/fnAhbP2o8lDZpHLc/8y/s5i7kF2/m5dvupjCzwzx2Jjqw5ozBng4MVEpdAYwFHtM07cGGY2c/W7wQQgghRIMBXYL576+GYdBpXPfKKtbsLTnje2iaxm8G/oaHBz7MNznfcO8P91LrakaPeEw/uHOxbz/vdt90fg094P4xEdzx1nN0SxqJw13Eu//3B9a9/dEpbiZE8wK2TilVDaCUysEXsi/RNO0ZJGALIYQQooV0i/Tno7uHEx5g5pY31/Lt9kNndZ9be93KEyOe4MdDP3L7N7dT5ig7/UX2cJjyGQy8zTcue8714KgAQK/Xc9lTv2PCNXehoWfJV7P4769n4vF0yBmLRQfQnIB9WNO0fkfeNITtSUAY0LuV6hJCCCHEeSg2yMpHdw0nPcqfu95dz4frDpzVfa5IuYJ/jf0X2eXZTP16KgXVBae/yGCCyc/Cpc/AnoXw2njftH4Nel97Cbf+7WlsxlhyC9bxytR7KNuTe1b1ic6tOQF7CtDkV0illFspNQUY3SpVCSGEEOK8FeJn4v3pQxmREsb/fLSFl5fsOav7XNDlAl6+8GWKaou45atb2Fu+t3kXDrrd15tdVwqvjfMtUNMgOCGWO2e9SGLsEOpcBbz9v//D5rkLzqo+0Xk1J2DnKaVO+G80SqkVANpp58IRQgghhGg+P7OBN6YOYlKfaJ76agdPfpmF13tmq08DZERl8PbEt3F73Uz9eipbi7Y278LEEb5x2cEJ8P4vYfmzviXX8Q0ZufqZx7hg0m0oFN9/8iqfPPoXGTIiGjVroRlN0+7XNK3LsY2appk0TRunadosYGrrlCeEEEKI85XJoOPf1/dnyrAEXl26l0c/2oLL4z3j+6SGpPLOJe9gN9q5/dvbWZm/snkXBnWBad9Azyvg+xm+Zdbrjz40OeCWK5ny+FNYDJHszV3F67feR8XBw2dcn+h8mhOwJwIeYI6mafmapmVqmrYXyAZuAJ5VSr3dijUKIYQQ4jyl12k8fllPHrqwO/M2HOSud9ZTV3/mPcXxAfHMvmQ28f7x3PvDvXyd83XzLjT5wTVvwfg/wbZ58NZEKD86LjwsNYm73nqJuMgBVNcf5M1HHmL7/O/OuD7RuZx2JccmJ2uaEd/DjXVKqfLWKurnkJUchRBCiM7pndX7+dOn2xjYJZg3pg4i0GY843tU1ldy/w/3s7FwI38Y8geuT7u++Rfv/NrXi220+FZ+TBjW5PDa1z5g+Q/zUKqe6//nCWIHylwQnd1Zr+R4zA2+A3oopQo6argWQgghROd1y9AEXrhhAJsPlnPdq6s4XOk443sEmAJ45aJXGBM3hr+s+Qv/2fQfmt3ZmDoRpv8A5gCYNRnWvdXk8ODp13HZ7Q8Abpa8MfeMaxOdR7MDNvBb4FlN097SNC26tQoSQgghhDiZS/tE89atg8ktreXq/6xkX3HNGd/DYrDwrwv+xWXJl/HS5pd4cs2TeFUzx3aHp8L0hdB1DCz4NSz4DbiPrhiZctFILIYoDpXuxlV75r8AiM6h2QFbKbVBKXUBsAD4WtO0GZqmWVuvNCGEEEKI443sFsac6UOprfdwzX9Wsi2v4ozvYdAZeGLEE0ztMZW5O+fy26W/xeVxNe9iaxDc+CEMfwDWvQHvXAHVRY2H0/tkoFQNy16YfcZ1ic7hTHqwj0zHtxP4D3A/kK1p2i2tUZgQQgghxMn0jQ/iv3cNw2LUc/2rq1m5u/iM76HTdDwy6BEeGvgQX+d8zX0L72ve0uoAOj1MeAKueg3y1sNrF0DBZgBG3X8rmmYja/OPZ1yT6BzOZAz2CiAP+BcQC9yKb9n0wZqmvdoaxQkhhBBCnExyuJ15dw8nJsjCrW/9yFdbm7Fa4wlM6zWN/xv+f6wuWM30b6dT7ihv/sV9fgm3fQXKC29cDNvmYbRZiAxOweEuYN+StWdVkzi3nUkP9p1ArFLqIqXUY0qpBUqp3Uqp+4FRrVSfEEIIIcRJRQVa+PBXw+gVG8C972/g/TVnt3T5ld2u5Jmxz7CjdAdTv57KoZoTrrF3YrEDYPoiiO4DH02DH55g7O03ABrL3//orOoR57YzGYO9XZ38MdtLW6geIYQQQogzEmQz8d4dQxnTPZw/fLKVFxZmN39mkGOM7zKely96mcO1h31Lq1c0c2l1AP9ImPo59LsJlj1NbFglNmMMRRV7qSuvOuNaxLntjMZgn4xS6gz+BAohhBBCtCyrSc+rUzK4sn8sT3+7i8c/zzyrpdUHRQ3irYvfot5Tz9SvprKteFvzLzaY4RdPgyUQfnyNXoOHo5SDpc+9dfprRadyxgFb07TJrVGIEEIIIcTPYdTr+Oe1fZk2Iom3V+bw6w82Ue8+86XV00PTmX3JbPyMfkz7Zhqr8lc1/2KTDfrfAlmfM/ymceg0O7t2bDrjGsS57Wx6sP/S4lUIIYQQQrQAnU7jsUnp/M/EVD7bnM8ds9dRW+8+4/skBCQw+5LZxPnHcc8P9/BNzjfNvzhjGnjd6Le+R0xEd+o9hez8askZ1yDOXWcTsLUWr0IIIYQQooVomsY9Y1N46qreLM8u4sbX1lBWU3/6C38iwhbBWxe/Re+w3jy65FE+3Plh8y4MTYaUC2HdW1xw542AjlXzPj3jzxfnrrMJ2Gc+oEkIIYQQoo1dP7gL/7l5IJkFlVz7yiryy+vO+B6B5kBeuegVRsWN4onVT/Dy5peb9wDl4Duh+hARZGI3x1FSlUN1YclZfAtxLmqRhxyFEEIIITqii3tGMXvaYA5XOLjmPyvZXVh9xvewGqw8e8GzXJZ8GS9uepGn1j51+qXVUy6EoARY+xr9R40B6ln8b3nY8XwhAVsIIYQQndrQrqHMuXMo9R4v1768kk0Hys/4HkadkSdGPMGUHlN4f8f7/H7Z70+9tLpOD4Nuh9yVDLw0Hb0ukD17t579lxDnlLMJ2IdbvAohhBBCiFbUKzaQj+4ajt1i4MbXVrN0V9EZ30On6Xgk4xEeHPAgX+77kvsX3X/qpdX73wIGC/r1b9IlJh23t4QtH375M76FOFecccBWSl3UGoUIIYQQQrSmxDA/5t01nC4hNm6f9SOfb84/43tomsYdve9g5rCZrMpfxZ3f3UmFs+LEJ9tCoNc1sOUDLph+FWBg7Zdf/7wvIc4JMkRECCGEEOeNiAALH/xqGP3jg3lg7kZmr8o5q/tc3f1q/jnmn2SWZHLr17dyuOYk/8A/eDq4agkuX0KgNZ6KulzK9ued/RcQ5wQJ2EIIIYQ4rwRajcy+fTDj0yL506fbeea7XWe1tPqFCRfy8oUvU1BTwC1f3UJORc7xJ8X0g7hBsPY1MiZcBLhZ8uLbP/MbiI7ubFZy9NM0Td8axQghhBBCtAWLUc/LNw/g2oFxPPdDNo99ug3PWSytPjh6MG9e/CZOj5MpX01he8n2408aNB1K99B7oBWDLoScA1l4PJ4W+BaiozptwNY0Tadp2o2apn2haVohsAMo0DQtU9O0f2ialtL6ZQohhBBCtCyDXsffr+nDr8Z05d3VuTwwZyNO95kH3x6hPZg1cRZWg5VpX09jTcGapif0vAJsYeg3vElSYi883nI2viMLz3RmzenBXgQkA78HopRS8UqpCGAksBr4m6ZpN7dijUIIIYQQrULTNH5/STp/+EUaX2wtYNrbP1LtPPOl1RMDE5l9yWxi7DHc8/095Fcf8wClwQwDp8LOr7hg6kWAiQ0Lf2i5LyE6nOYE7AuVUk8opbYodXRWdaVUqVJqnlLqauCD1itRCCGEEKJ13Tk6maev7cvqvaXc+NpqSqqdZ3yPSL9I/nPhf/AqL+9mvdv0YMY00DT8D8wn2C+BKudBinbsbaHqRUdz2oCtlHJpmpamadp4TdPsxx7TNG3ikXNaq0AhhBBCiLZwzcA4Xr1lIDsPVXHty6s4WHaKOa5PIsovigmJE/g4+2Oq6quOHgiMg9RfwIbZDL1sIuBh8cuzW6540aE0Zwz2A8CnwP3ANk3TLj/m8JOtVZgQQgghRFsbnx7Ju3cMobjaydX/Wcmuw1Wnv+gnpvacSo2rhnm75jU9MHg61JXSI7EEoz6cg4d24amXPsrOqDlDRKYDA5VSVwBjgcc0TXuw4ZjWSnUJIYQQQrSLQYkhfHjXMJSCa19exfr9pWd0fY/QHgyKGsS7We/i8h4ToJPGQFh3WPsa3bv3xasqWfPa3BauXnQEzQnYOqVUNYBSKgdfyL5E07RnkIAthBBCiE4oLSqAeXcPJ9hm5KbX17BoR+EZXT+1x1QO1x7m25xvjzZqmm/KvvwNjL52AJpmYfOqFS1cuegImhOwD2ua1u/Im4awPQkIA3q3Ul1CCCGEEO0qPsTGR3cPJyXCzvTZ6/hk48FmXzsqbhSJAYnM2j6r6SI2fa8Hkx3brvcJ9U+i1pVH3voTzJ0tzmnNCdhTgEPHNiil3EqpKcDoVqlKCCGEEKIDCLObmTN9KIMSQ3jog828sXxfs67TaTqm9JxCVmkW6w6vO3rAEuAL2ds+ZuQVFwKKpW+93zrFi3bTnICdp5Q6dKIDSqkVAJqmyVARIYQQQnRK/hYjb902iEt6RfHEgkz+/vWOZi2tPrnrZEIsIcze/pPZQgbdAR4nyQHbMRuiKCjejct55tMCio6rWQvNaJp2v6ZpXY5t1DTNpGnaOE3TZgFTW6c8IYQQQoj2ZzHqeeHGAdwwuAsvLd7D7+Ztxe3xnvoag4XrUq9j8cHF7Ks4puc7Ih0SR8G6N0nv1R+lalj+vEzZ15k0J2D7JmuEOZqm5Tcskb4XyAZuAJ5VSr3dijUKIYQQQrQ7vU7jySt7cf+4FD5Yd4B73tuAw3XqpdWvS70Ok87EO5nvND0weDpUHGDkxbFomo3MjWtbsXLR1pqz0IxDKfWSUmoEkACMBwYopRKUUtOVUhtbvUohhBBCiA5A0zQenpDKjMk9+DbzMFPfXEul4+RzWYdaQ5mcPJnP9nxGqeOY6f5SLwX/GMzb3iYyOBmHu4Cc5evb4BuItnDKgK1pmv7Y90opl1KqQClV3qpVCSGEEEJ0YLeNSOLf1/dj/f4yrn9lNUVVJx9DPaXHFJweJx/s/OBoo97gWz597yJGXzUa0Fj+zoetX7hoE6frwf67pmkvAGia9r2maX3boCYhhBBCiA7v8n6xvD41g33FNVzz8kpyS068tHrXoK6Mih3F3B1zcXqOCeIDp4LOSLx7MVZjDIUVe3FWnvnKkaLjOV3ALgP2N7z+H+BZTdPe0jQtunXLEkIIIYTo+MamRvD+9CFU1Lm4+uWVZOZXnvC8qT2nUuooZcGeBUcb7RHQ8wrY9D69+g9EqTqWPDerbQoXrep0AXsw4ABQSm1QSl0ALAC+1jRthqZp1tYuUAghhBCiI+vfJZj//moYBp3Gda+sYs3ekuPOGRw1mLSQNGZnzsarjpl9ZNB0cFYybJiGTrOzM3NDG1YuWsspA7ZS6jLg3SPvG+a73gn8B7gfyNY07ZZWrVAIIYQQooPrFunPR3cPJzzAzNS31nKwrOlwEU3TmNJjCnsr9rIi75jl0eMHQ1RvjJveJDq8G/WeQnZ9vbSNqxctrTmziJQBaJq2AsgD/gXEArcCY4HBmqa92nolNo+maZM1TXu1oqKivUsRQgghxHkoNsjKO7cPweNVvLxkz3HHJyZNJMIWwazMY4aBaBoMvhMKMxl7WV9Ax8qPPm27okWraM482EfcCcQqpS5SSj2mlFqglNqtlLofGNVK9TWbUupzpdSdgYGB7V2KEEIIIc5TsUFWrhkYz4c/HuRQhaPJMaPOyE3pN7GmYA07SnccPdDrGrAEEVX2OXZTLCVVOdQWl7Vx5aIlnTZga5rWpWEVxyog/sj7YzfgF5qmBbR6tUIIIYQQHdw9Y5PxKMWrS/ced+ya7tdgM9iaLp9uskH/m2HHAvoOGQg4WfTvN9uuYNHimtODPQt4u2F/su0N4IpWqVAIIYQQ4hwSH2Ljin6xvL92P8XVTefHDjAFcFW3q/hq31ccrjl89EDGNPB6GJR2GJ0WwJ7dW9u4atGSmjMG+wKl1LiG/cm2cUqp2ae7lxBCCCHE+eDeC5Kpd3t5bdnxvdg3pd+EFy/v73j/aGNoMqRciH7TLLpEp+LyFrNt3tdtWLFoSWcyBlsIIYQQQjRD13A7k/rE8O6q/ZTV1Dc5Fucfx/gu4/nvrv9S6zpmtpHBd0L1YcZN7ALoWbPgq7YtWrQYCdhCCCGEEK3g3gtSqKn38NaKfccdm9pzKlX1VXyy+5OjjSkXQnAiwXkfEGDpQnltLhUHDx93rej4zihga5o27ti9EEIIIYQ4sdQofyb2jOKtlTlUOlxNjvUN70u/8H68k/kOHq/H16jTQcbtkLuKgcP6Ai4WPy8PO56LzrQH++mf7IUQQgghxEncNy6FKoeb2Stzjjs2tedU8qrz+CH3h6ON/W8Gg4W+UZsx6ILZl5uJx+Npu4JFizjbISJai1YhhBBCCNEJ9YoNZFxaBG8s30eN093k2AXxFxDvH9904RlbCPS+Bv32j0iKT8PjLWPz+5+1cdXi55Ix2EIIIYQQrej+cSmU1bp4d/X+Ju16nZ6b029mS9EWNhVuOnpg0HRw1TJmtBUwsu77HxDnFgnYQgghhBCtqH+XYEZ1C+O1ZXtxuJoO97gi5QoCTAHMzjxmtuOYfhA3mMC97xJkS6DKcYCS7Jw2rVn8PBKwhRBCCCFa2f3julFcXc+ctblN2m1GG79M/SU/5P7AgaoDRw8Mng6lexgyLAXwsOglWW7kXHKmAbu6YV/V0oUIIYQQQnRWg5NCGJwUwitL9uJ0N+3FviHtBnSajncz3z3a2ONy8Aunl305Rl0YBwp2yMOO55AzCthKqdHH7oUQQgghRPM8MK4bhyodfLT+YJP2CFsEv0j6BZ/s/oQKZ4Wv0WCGAVNh19ekJKXiVZWsfW1uO1QtzoYMERFCCCGEaAMjUkLp3yWI/yzeg8vjbXJsSo8p1Lnr+O+u/x5tzLgNNB1jM6oBM5tXLG/bgsVZk4AthBBCCNEGNE3j/nEpHCyr45ONeU2OpYakMjR6KHOy5uDyNCxKExgHab/Alv0+Yf6J1NTncWhzVjtULs7UGQdsTdP8NE3Tt0YxQgghhBCd2QWpEfSKDeClRbvxeFWTY1N7TqWwrpCvcr462jhoOtSVMmJYNOBl8Wvvt23B4qycNmBrmqbTNO1GTdO+0DStENgBFGialqlp2j80TUtp/TKFEEIIIc59mqZx3wXdyCmpZcGW/CbHRsSMICUohVnbZ6FUQ/hOGg1hqaS4P8esjyS/eBcup7MdKhdnojk92IuAZOD3QJRSKl4pFQGMBFYDf9M07eZWrFEIIYQQotOY0COS1Eh/Xli4G+8xvdiapjGlxxR2le1idcHqI42+KfsKNpGanIxSNax88d2T3Fl0FM0J2BcqpZ5QSm1RSjWOyFdKlSql5imlrgY+aL0ShRBCCCE6D51O495xKWQXVvPN9kNNjl3a9VJCLaFNF57pcx2Y7IzusQ9Ns7J9w5o2rlicqdMGbKWUS9O0NE3TxmuaZj/2mKZpE4+c01oFCiGEEEJ0Npf2jqZrmB/PL9x9dDgIYNKbuCHtBpbnLWdP+R5foyUA+t6Aefd8IgKSqHPlk7tiQztVLpqjOWOwHwA+Be4Htmmadvkxh59srcKEEEIIITorvU7jngtSyCyo5IeswibHrku9Dove0rQXe9Ad4Kln9BArAEvf+bAtyxVnqDlDRKYDA5VSVwBjgcc0TXuw4ZjWSnUJIYQQQnRql/eLIT7EyvOLmvZiB1mCuDzlcj7f8znFdcW+xog0SBpNl8p5WI0xFJbvwVld006Vi9NpTsDWKaWqAZRSOfhC9iWapj2DBGwhhBBCiLNi1Ou4Z2wKmw+Usyy7uMmxm9Nvxu11M3fHMas3DpoOFQfomRKHUnUs/ffbbVuwaLbmBOzDmqb1O/KmIWxPAsKA3q1UlxBCCCFEp3fVgFiiAy08vzC7SS92YmAiY+LH8MHOD6hz1/kaU38BAbEM77IOTfNj53YZh91RNSdgTwGaPOKqlHIrpaYAo1ulKiGEEEKI84DZoOeuMcn8mFPGmn2lTY5N7TGVcmc5n+/53NegN0DGbRhzlxATnIjTc5jd38ny6R1RcwJ2nlLq0IkOKKVWAGiaJkNFhBBCCCHOwnWD4gn3N/P8wuwm7QMjB9IztCfvZL6D98hMyQOmgs7I2P71gI4VH8xv83rF6TVroRlN0+7XNK3LsY2appk0TRunadosYGrrlCeEEEII0blZjHruHNWVFbtLWL+/rLFd0zSm9pxKTmUOSw4s8TXaI6DnlUSVfISfKYbiqn3UFped5M6ivTQnYE8EPMAcTdPyG5ZI3wdkAzcAzyql3m7FGoUQQgghOrWbhnYhxM/ECz/pxb4o4SKi/aKZlTnraOPg6eCspG+3cMDJ4uffattixWk1Z6EZh1LqJaXUCCABGA/0V0olKKWmK6U2tnqVQgghhBCdmM1k4PaRSSzaWcTWgxWN7QadgZvSb2L94fVsL97ua4wbBFF9GBz0LTotgN27trRT1eJkmrPQzDRN08wNb38BXAb0aNWqhBBCCCHOM1OGJRBgMRw3FvvqbldjN9qP9mJrGgy+E33pDuJD43F5i9k+/7t2qFicTHOGiDyolHJqmjYT+A2QBMzQNG2lpmlRrVqdEEIIIcR5wt9i5LYRSXybeZgdhyob2+0mO1d3u5pvc76loLrA19jrarAEcUGvQ4Ce1Z990T5FixNqTsCub9j/AhinlPqdUupi4C/AS61WmRBCCCHEeea2EYnYzQZeWLi7SftN6TcB8F7We74Gkw3630xo4ef4m+Mor9lPxcHDbV2uOInmBOwDmqa9DUQA1iONSqkv8PVmCyGEEEKIFhBkM3HLsAS+2FrA7sLqxvZoezQTEicwL3se1fUN7YNuB+VlYHcb4GLJC/KwY0fRnIB9K7AEuByYp2naQ5qmTdA07bcc7d0WQgghhBAt4I6RSVgMel5a3LQXe2qPqVS7qpmXPc/XENIVul1EP8N89Lpg9u3fjsfjaYeKxU81ZxaRSqXUW0qpzcC1gAFf6O4CXNe65QkhhBBCnF9C7WZuGtKFTzflk1tS29jeM6wnAyMH8l7We7i9bl/j4DvROwpJCovG7S1j85zP26lqcazm9GA3agjb/1BK3aiUulcpldNKdQkhhBBCnLfuHN0VvU47YS92QU0B3+1vmDUkeTwEJzE2ORMwsv67H9q+WHGcMwrYQgghhBCi9UUEWLh+UDzzNhwkr7yusX1M/BgSAxKZtX0WSinQ6WDQ7QSWLyPIEkulI5eyPbntWLkACdhCCCGEEB3Sr8YkA/DKkj2NbTpNxy09bmF7yXbWH17va+x3ExisDE72Ah4WvjjrBHcTbUkCthBCCCFEBxQbZOXqAXHM/fEAhZWOxvbJyZMJMgcdXXjGFgK9r6G37mOMulAO5O+Qhx3bmQRsIYQQQogO6p6xKXi8ileW7m1ssxqsXJd6HUsOLCGnIsfXOHg6uOtIDg/Foyr48fUP26dgAUjAFkIIIYTosLqE2ri8XwzvrdlPSbWzsf36tOsx6oy8m/WuryG6L8QP4YK4ZYCJzcuXtk/BApCALYQQQgjRod0zNgWn28vry/c1toVZw5iUPIlPd39KuaPc1zhoOrbaXYTaYqmuz+PQ1p3tU7CQgC2EEEII0ZGlRNi5tHc0s1fmUF57dI2/KT2m4PA4+GDnB76GHpeDXwTDk8oAL0tefa99ChYSsIUQQgghOrr7xqVQU+/hrRU5jW3JQcmMjB3JnB1zcHqcYDDBwKl09y7ApI8gv2gXLqfz5DcVrUYCthBCCCFEB5cWFcCEHpG8tWIfVQ5XY/uUHlMocZTw5d4vfQ0DbwNNR2qEH15Vzar/SC92e5CALYQQQghxDrh/XDcqHW5mr9rf2DY0eijdg7szO3O2b+GZwFhIu5QxoV+iaVa2rVvTjhWfvyRgCyGEEEKcA3rHBXJBajivL9tLbb0bAE3TmNpzKrvLd7Mif4XvxMHTMXuKCbdFUefK58CaTe1X9HlKArYQQgghxDnivnHdKKt18d7qo8uhX5J4CRHWCGZtb1h4JnEUhKcxKn4foFj69gftU+x5TAK2EEIIIcQ5YmBCMCNSQnll6V4cLt9qjUa9kRvSb2B1wWp2lu4ETYNBd5CoLcNiiOJw2R6c1TXtXPn5RQK2EEIIIcQ55P5x3SiudvLBjwca267tfi1Wg5XZmbN9DX2vB5M/PcL1KFXL8udnt1O15ycJ2EIIIYQQ55AhSSEMSgzm5SV7cLp9vdiB5kCuTLmSL/d9SWFtIZj9od8NjAz8BE3zI2vb+nau+vwiAVsIIYQQ4hyiaRr3j+tGQYWDeevzGttv7nEzXuVlzo45voZBd2DUnET7heN0H2LPwlXtVPH5RwK2+P/27iw4yvNQ8/jzdrf2Da0sEos2FoHZLDDYYATEcZwYLzhxbAJibAcSOyRT59xMztWpuZiai6mZmsRgZ/CxwxLwRkhs4jjYMQiwsTG72ZGQWCSxCNC+q/s9FzBO4nhhaentbv1/Va6Cr0T3o8L6+NdbLTUAAAgzMwszNGHoAL1QVqFuf0CSNDRpqOYOm6s3Tryhtu42KXOUlDtLJYP3SjL68NWNbkf3IwQ2AABAmDHG6OezC1Rd3663DtR+fr20qFRNXU36Y8Ufr12YukSDfUcV7xukK81VaqtvdDO4nyGwAQAAwtDcMVkqGpysF7ZWyB+wkqSJWRM1IXOC1h5dK3/AL418QErO0YSMDlnboe2/+q3j1f0DgQ0AABCGrr0Wu0CVl1v1zqHzn18vLSpVdUu1tp7bKnl9UvFTuit5kzwmSSdPHnS4uP8gsAEAAMLU/WMHqTArUcu3lCtw/RR77rC5yk7M/tsbz0xeLG90lHISU9Xtr9PRtz9wuLh/ILABAADClMdjtGxOgU5ebNF7Ry9IkrwerxYVLdKBugM6WHdQSsyUxj6qksxtkrz65I9/cju6HyCwAQAAwtiD44coNyNBz2+pkLXXTrEfLXhUSdFJfzvFnrJEmTHVSooaqPrW02quveRwceQjsAEAAMKY12P0bEm+jtQ2aeuJa+EcHxWvH4z8gT44+4Gqm6ulnGJp8ARNTquT1K2tz69yujnSEdgAAABh7tFJ2cpJjdOvP/jbKfaC0QvkkUfrjq2TjJGmLtWklC3ymgGqOn3Y8eLIRmADAACEuSivR8+W5OvAuQZ9VHFFkjQwYaAeyH1AG8s3qqmrSRr3mLyJqRqRFKeewFUdeHWT49WRi8AGAACIAN+/M0eDkmP16y3ln19bPHax2nratOHkBikqTpq0ULPTN0vyac/m99yNjXAENgAAQASI8Xn1k1l5+rTqqnZVXjvFHpU2SncNvkvrjq1Tt79bKn5GKbFNSokaqMb2c6qvOud4dWQisAEAACLEk1OHKSMxWsu3Vnx+bXHRYl1qu6S/nP6LlJYrFX5bU9JOSerR1uWr3Y2NYAQ2AABAhIiN8mrJzDztKL+s/WfrJUn3ZN+jvJQ8rTm65to3QE5dqglpu+XzpOls7TH5/X7HqyMPgQ0AABBBFk4brtT4KD2/5doptsd4VFpUquNXj+vTC59K+XOktDzlJxj5A43a+8oGx4sjD4ENAAAQQRJifHpmRq62HL+kwzWNkqQH8x9UWmzatTee8Xik4mdUkrlZUrT279jmdnAEIrABAAAiTOndI5QU69Py66fYMd4YPTH6Ce2o2aHKhkpp0o+UmCClRWeopbNal46cdLw4shDYAAAAESY5NkpP3T1CfzlyQScuNEuSfjjqh4rxxmjN0TVSXKo0/ge6e8BBSQGV/b91bgdHGAIbAAAgAj11T64Sor1acf0niqTFpumh/Ie06dQmXWm/Ik1ZolHpxxXtyVDNpZPyd3U7Xhw5CGwAAIAIlJoQrYXTh+tPn9Wqsq5FkrSoaJG6Al16/cTr0uDx0tBpGhnfroBt1s7frHe8OHIQ2AAAABFqycw8Rfs8WrH1lCQpNyVXJTkleu34a+ro6ZCmLtG9Az+QMbE6/OlOx2sjB4ENAAAQoTISY7Rg6nD98UCNzl1tkySVji1VfWe9NlVuksY8pLjUVGVEpaqtu1bVn37meHFkILABAAAi2NJ78+Q1Ri+UXTvFLh5YrKL0Iq05skYBr0+6879oZuonkqy2r3rN7dgIQWADAABEsEEpsXp8So427D2n2oZ2GWO0uGixTjed1o7qHVLxU8pNr1GsJ0sXrlaou63D9eSwR2ADAABEuJ/Oype10srtlZKk+0bcp4HxA7X66GopeYg05kGNib8qa9u04/lVbsdGAAIbAAAgwuWkxmv+5Gy9+ulZXWruUJQnSgvHLNTuC7t19MpRacoSzRy8TcbE69hne1zPDXsENgAAQD/wXEmBuv0BvXT9FPuxkY8pISrh2tunj5ihqMGjNCgqUR09F1S17VPHa8MbgQ0AANAPjMhI0MMTs/W7T87qamuXkqKTNL9wvjaf3qwLbRelqT/WrPSPJBntWPem67lhjcAGAADoJ342O18dPX69/OG1U+yFYxZKktYdWyeN/6Gys1oV78nS5aYqtTc0u5wa1ghsAACAfqIgK0nfHTdYq3eeUWNbt4YkDtF9w+/ThpMb1OrxSBMX6I6EalnboW2/+q3ruWGLwAYAAOhHls0pUEtnj1btPC1JWjx2sVq6W7SxfKM05ceaPniXPCZR5ScOON0ZzghsAACAfmTM4GR9a8xAvfJRlZo7ujUuY5wmZ03W747+Tj3pefIWzlK2L0pd/ks68U6Z67lhicAGAADoZ34xt0CN7d1a+8kZSddOsWtba/XXs3+VpixRycCPJHm0c+PbboeGKQIbAACgnxmfM0CzRmbqP3ZUqa2rRyVDSzQ8ebjWHFkjW3i/srLjlOjN1NWW02q5UOd6btghsAEAAPqhX8wt0NXWLq3fdVYe49GiMYt06PIh7b9ySJrytCbFn5LUpTLe2fGmEdgAAAD90J3D0zQ9L10rt1eqo9uvhwoeUkpMyrU3npm8WHfmHJHXpOhU5WHXU8NOSAe2MSbBGLPaGPOSMeZHrvcAAABEkp/PLdCl5k69ueec4nxxenzk49p6bqvO+FvlHT9fw6MC6glc0cHX33E9Naz0eWAbY14xxlwyxhz+wvXvGGNOGGMqjDG/vH55vqQN1tolkh7q660AAACRbHpeuoqHp+rFslPq6glowZgF8nl8Wnt0rTR1iUoG75Dk0+53N7ueGlZcnGCvkvSdv79gjPFKWiHpAUlFkp40xhRJypF07vqH+ftwIwAAQMQzxmjZnALVNnZo475qZcRl6Ht539NbFW+pIT1fqfn5SvGmqbH9rOrP1LieGzb6PLCttdslXf3C5amSKqy1ldbaLkmvSXpYUrWuRbYU4i9nAQAACEezRmZqfE6KXig7pR5/QKVFperwd+jN8g3S1KWaknRcUo/Klq9yPTVshEq0ZutvJ9XStbDOlrRR0mPGmBclbfqqP2yMWWqM2WOM2VNXx4+SAQAAuFHGGP18TqHOXm3T2wdrVZhaqHuG3KP1x9era8yDmpBXK59J1ZnqY/L7eUHBjQiVwP5S1tpWa+1T1tpnrbXrvubjVlpri621xZmZmX05EQAAIOx9a0yWRg9K0vKtFfIHrErHlupy+2X9uXqrNGmR8qLa5A80aN/qja6nhoVQCewaSUP/7vc5168BAACgl/3/U+zKula9e/i8pg+erpGpI7Xm6BrZ4qc1O+djSdHaV1bmempYCJXA3i2p0BiTa4yJlvSEJN6bEwAAoI88MG6QCrIStXxLhayVSotKVV5fro/bzyvxjplK8ySrpfOc6o5Xup4a8lz8mL5XJX0saZQxptoY84y1tkfSMkmbJR2T9Ia19khfbwMAAOivPB6jn83O1/ELzXr/2EV9N/e7yozL1Oqjq6WpSzQ99ZCkgMpeXON6ashz8VNEnrTWDrbWRllrc6y1L1+//mdr7Uhrbb619n/09S4AAID+bt74IRqeHq/lWyrk8/j05OgntbN2p06mDtXo0QFFm3RVXzwhf1e366khLVReIgIAAADHfF6PflZSoEM1jSo7WafHRz2uOF+c1hxbK035sQqjGxSwzfpk5auup4Y0AhsAAACfe2RStrIHxOn5D8qVHJ2sh/Mf1jtV76hu1Lc1K3e/jGL12Sc7Xc8MaQQ2AAAAPhft8+inJfnad7ZBH5+6okVFi+QP+PXq6XcVd+cjyvDEq627RjV7+Xa5r0JgAwAA4B/84M4cDUyO0a+3lGtY8jDNGTZHr594XW2TF+mejAOSrLa9vN71zJBFYAMAAOAfxEZ5tfTefH1SeVW7T1/V4rGL1dTVpLdbTil/0hDFmgxduFqh7rYO11NDEoENAACAf7Jg6jBlJEbr+S0Vmpg5UeMzxmvt0bXyFz+t0bF1srZVH65Y63pmSCKwAQAA8E/ior368cw8bT9Zp4PVjSodW6qzzWdVlpSiGaNPyZh4HT3wqeuZIYnABgAAwJdaOG24BsRHafmWcs0dNlfZidlac3y9Yu4q1UATrY6e86ravtv1zJATMYFtjJlnjFnZ2NjoegoAAEBESIzx6el7cvXXY5d04kKrFo5ZqH2X9ulQ7jTNGrRPktGH6950PTPkRExgW2s3WWuXpqSkuJ4CAAAQMRbfPUJJMT6t2FqhRwsfVVJUklaf+bNypk9RnElXXWOlOpuaXc8MKRET2AAAAAi+lLgoLb57hN49fEG1VwP6/sjv6/0z76tm/HzdEV8razu07VerXM8MKQQ2AAAAvtbTM3IVF+XViq0VWjBmgTzy6HdNx3X3pGZ5TKJOHNvvemJIIbABAADwtdISorVo2nC9fbBW7e1Juj/3fm0s36jWKQs1xONRl/+STry7zfXMkEFgAwAA4Bs9MzNXUV6PXthaocVFi9XW06bfx3o1e9ghSR59/Pu3XE8MGQQ2AAAAvlFWUqyenDpMf9hfo0QzXFMHTdW68g1Kvfc7SjRputJ8Wi2XrrieGRIIbAAAANyQn8zKk8cY/WbbKS0eu1gX2y7qvezRmph4RlKXyn79W9cTQwKBDQAAgBsyOCVO3y/O0Zt7qlWQWKzclFytPrtZd96dKK9JVuWpQ64nhgQCGwAAADfs2Vn5Clirl3ZUqbSoVMeuHtP+cfdpmM+v7sAVHXrzXdcTnSOwAQAAcMOGpsXr0UnZWr/rrKZnfVtpsWla3XRCswsrJfn06Z//4nqicwQ2AAAAbspzswvU7Q9ozc4aPTHqCW2r2a76GQ8o2QxQQ9sZNZ4773qiUwQ2AAAAbkpuRoLmTRiitZ+c0f1DH1W0J1pro7o1ZUClpB5tfb5/f7MjgQ0AAICbtmx2gdq6/PrD3kbNy5+nTWffV/acfPnMAJ0+d1R+v9/1RGcIbAAAANy0woFJemDcIK366LTm5y9Qp79TGwblKDeqQ/5Ag/av7b9vPENgAwAA4JYsm1Og5s4elR02ujfnXr12frvunlQvKUr7tn7gep4zERPYxph5xpiVjY2NrqcAAAD0C2OHpOhbY7L0ykdVerxwoa52XFXZhDuV6klWc0e1Lp+ocj3RiYgJbGvtJmvt0pSUFNdTAAAA+o1lcwrV0Nato5WZGpM2RmubjuuuzNOS/Cr7zRrX85yImMAGAABA35s4dIBmFmbo5Q+r9OSohapsqtKV+0YpyqTq3PkT8nd1u57Y5whsAAAA3JZfzC3U5ZYuXb44RgPjB2qNp1WF0S0K2Cbtevl11/P6HIENAACA2zJlRJruyk3Tf+w4ox+OfFK7Ln+mIfcGJMXos48+dD2vzxHYAAAAuG2/mFuoi02dUss0xfvi9dqQNGV4EtTaXavz+4+6ntenCGwAAADctrvz0zV52ACt2nFRjxTM118uH9C4EbWSAip7eZ3reX2KwAYAAMBtM8bo53MKVdPQrrSeOQoooG3TByvGpOv85Qp1d3a6nthnCGwAAAAERcmoTN2RnaL1O1s0d+i3tKGlQoVxjbK2VR8tX+t6Xp8hsAEAABAUxhgtm1OgM1falBv9XTV3t+jiXI+MidPRfbtcz+szBDYAAACC5r4xAzVqYJLe2uXTxMxJWh/doCwTq/ae8zr70T7X8/oEgQ0AAICg8XiunWJXXGrRuMR5qmmvU/LYq5KMtq/pHz8Tm8AGAABAUH33jsHKy0zQln2ZGpo0VK+PDijOpOtSY6U6m5pdz+t1BDYAAACCyusx+llJgU5caNWU1Ef0WWu1slMaZG27tj+/2vW8XkdgAwAAIOgenjhEw9LitfdIvpKjk7VtVpeMSdCJI5H/OmwCGwAAAEHn83r0XEm+Dld3aFrmPP21u1pZXp86/Zd0cvMO1/N6FYENAACAXjF/co6GpMTqVPkEeT1eXbyzXpJHO9/4o+tpvYrABgAAQK+I9nn005J8HTgT0JSMuXo986ISTJqutJxW2+V61/N6TcQEtjFmnjFmZWNjo+spAAAAuO7x4qHKTIpRXc00tQe6FJ3VKKlTZb/6retpvSZiAttau8lauzQlJcX1FAAAAFwXG+XVT+7N0/6KOI1NLdbvimvlMcmqqDjoelqviZjABgAAQGhacNcwpSVEq/PyTNXZNqVEG3UHrujwH95zPa1XENgAAADoVfHRPv14Zq72n8xSTkKePiyuleTVrrf+5HparyCwAQAA0OsWTRuulLhoRbeWaHfKZSV4UtXQfk6N1RddTws6AhsAAAC9Lik2Sk/dM0IHj+dpQHS6rgytl9StsudfcT0t6AhsAAAA9Imn7s5VYnSsBvSU6LXRp+Q1Kao6e0R+v9/1tKAisAEAANAnUuKjVDp9uA4fL1J0VKyi4iV/oEEH17/telpQEdgAAADoM8/MyFWsJ1lZZoZ+P6lcUpT2vv++61lBRWADAACgz6QnxmjhtGE6eXKSLid3KM6boqbOGl0pP+16WtAQ2AAAAOhTS2bmyRvIUpZ3so7kXpTkV9mKVa5nBQ2BDQAAgD6VlRyrJ6cM1dmqKdo+slo+k6qzF05GzDc7EtgAAADocz+ZlS/bMUIDvPnqSOpRwDbp05dedz0rKAhsAAAA9LkhA+L0/TuHqq76Lv1h8nFJMTr44TbXs4KCwAYAAIATz84qUHfTONmEAYryJam1+7wuHDrhetZtI7ABAADgxLD0eD0ycZiaL03TrsJzkgIqe/G3rmfdNgIbAAAAzjw3O1/t9cU6ntskn0lX7dUqdXd2up51WwhsAAAAOJOfmagHx+Wps+EuXU5tl7Wt+njFGtezbguBDQAAAKeWzS5Qa910bZ5ULmPidHjPLteTbguBDQAAAKdGDUrSt0eNUnv3ONmoBLX7L+rcx/tdz7plBDYAAACc+/mcQrXWzdDO0ZWSrLa9stb1pFsWMYFtjJlnjFnZ2NjoegoAAABu0rjsFJWMmKxTGSnyeTJ0qfmcOltaXc+6JRET2NbaTdbapSkpKa6nAAAA4BYsm1Og5rp7VZveLGvbteP/vuR60i2JmMAGAABAeJs8LFXTBs3Q1nF1MiZBx48cdD3plhDYAAAACBm/mDNS9Q0l6omOVWegTqfe/9D1pJtGYAMAACBk3JWXrgmpc/XxmDOSPPpw/RuuJ900AhsAAAAh5b/OKdKxqPHyetJ1ua1GbfXh9UMsCGwAAACElBkFGSqIf0DnBtZL6tS2/7XC9aSbQmADAAAgpBhj9K9zJmvL0Dh5TJLKq064nnRTCGwAAACEnNmjspQZ86g6Y33qDlzRkd+/43rSDSOwAQAAEHKMMfrXkpnaWdggyatdb/3J9aQbRmADAAAgJH27aJAakx6U15um+s4Laq656HrSDSGwAQAAEJI8HqN/mfk9nRnYKKlbW//3864n3RACGwAAACFr3oRs7Rs+SR6Tosra0/L7/a4nfSMCGwAAACHL6zF67p6F6og38tsGHVrzuutJ34jABgAAQEibP3m4dhcmS/Jp11+3up7zjQhsAAAAhLQor0ffm/UzeX2paumpU315petJX4vABgAAQMhbOHW0qgd1S+rRe/9nues5X4vABgAAQMiL8Xk19MFn5TWpqq2/ENLf7EhgAwAAICz8ZMY0dcYbBWyTPvn1C67nfCUCGwAAAGEhNsorz/1zJEVr/+79rud8JQIbAAAAYeO5RxbK60tRp/+yzu/d63rOlyKwAQAAEDYSYnwKFA2SFNA7K15yPedLEdgAAAAIK0//y7/L60lXU1uD/J1druf8EwIbAAAAYWVAfKx8A2JlbYs2/Pd/dz3nnxDYAAAACDsP/vK/yZhYXay64HrKP4mYwDbGzDPGrGxsbHQ9BQAAAL1sxPA8RfsGqMe2qaaq3PWcf2Csta43BFVxcbHds2eP6xkAAADoZVVHDmpA1mClZmY5eX5jzF5rbfEXr/tcjAEAAABuV+7YCa4nfKmIeYkIAAAAEAoIbAAAACCICGwAAAAgiAhsAAAAIIgIbAAAACCICGwAAAAgiAhsAAAAIIgIbAAAACCICGwAAAAgiAhsAAAAIIgIbAAAACCICGwAAAAgiAhsAAAAIIgIbAAAACCICGwAAAAgiAhsAAAAIIgIbAAAACCIjLXW9YagMsbUSTrTx0+bIqmxj5+zL4Tq5+VqV28/b7AfP1iPdzuPc6t/NkPS5Vt8Tty8UP1av12h+nlxD+vbx+MeFvlcfq0Pt9ZmfvFixAW2C8aYldbapa53BFuofl6udvX28wb78YP1eLfzOLf6Z40xe6y1xbfynLh5ofq1frtC9fPiHta3j8c9LPKF4tc6LxEJjk2uB/SSUP28XO3q7ecN9uMH6/Fu53FC9f8h/KNI/XsK1c+Le1jfPh73sMgXcn9PnGADCDmc/gAIZ9zDwAk2gFC00vUAALgN3MP6OU6wAQAAgCDiBBsAAAAIIgIbAAAACCICGwAAAAgiAhtAWDHG5BljXjbGbHC9BQBuhDEmwRiz2hjzkjHmR673oPcR2AD6jDHmFWPMJWPM4S9c/44x5oQxpsIY88uvewxrbaW19pneXQoAX+8m72fzJW2w1i6R9FCfj0WfI7AB9KVVkr7z9xeMMV5JKyQ9IKlI0pPGmCJjzB3GmD994b+svp8MAF9qlW7wfiYpR9K56x/m78ONcMTnegCA/sNau90YM+ILl6dKqrDWVkqSMeY1SQ9ba/+npAf7eCIA3JCbuZ9Jqta1yD4gDjf7Bf6SAbiWrb+d7EjX/iHK/qoPNsakG2N+I2mSMebfenscANyEr7qfbZT0mDHmRYXg23oj+DjBBhBWrLVXJP3U9Q4AuFHW2lZJT7negb7DCTYA12okDf273+dcvwYA4Yb7GSQR2ADc2y2p0BiTa4yJlvSEpLcdbwKAW8H9DJIIbAB9yBjzqqSPJY0yxlQbY56x1vZIWiZps6Rjkt6w1h5xuRMAvgn3M3wdY611vQEAAACIGJxgAwAAAEFEYAMAAABBRGADAAAAQURgAwAAAEFEYAMAAABBRGADAAAAQURgA0A/ZowZYYw57HoHAEQSAhsAAAAIIgIbAOA1xrxkjDlijHnPGBPnehAAhDMCGwBQKGmFtXaspAZJj7mdAwDhjcAGAFRZaw9c//VeSSPcTQGA8EdgAwA6/+7Xfkk+V0MAIBIQ2AAAAEAQEdgAAABAEBlrresNAAAAQMTgBBsAAAAIIgIbAAAACCICGwAAAAgiAhsAAAAIIgIbAAAACCICGwAAAAgiAhsAAAAIIgIbAAAACKL/BJPY4ZYb/Dp2AAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "# We assume the definitions of the previous block as not to copy everything.\n", "\n", - "fig = plt.figure(figsize=(12, 8))\n", - "gs = gridspec.GridSpec(2, 1, height_ratios=[2, 1])\n", - "\n", - "ax0 = plt.subplot(gs[0])\n", - "ax0.plot(x, y, label=\"exact solution\")\n", - "ax0.set_yscale(\"log\")\n", - "ax0.set_title(\"Solutions to $y' = y - x^2 + 1.0$\")\n", - "\n", - "\n", - "ax1 = plt.subplot(gs[1])\n", - "ax1.set_ylabel(\"$\\delta(x) = |\\~y(x) - y(x)|$\")\n", + "h_list = 2/np.arange(1, 100)\n", + "y2 = np.zeros(len(h_list))\n", "\n", + "plt.figure(figsize=(12, 8))\n", + "plt.ylabel(\"$\\delta(2) = |\\~y(2) - y(2)|$\")\n", "\n", "for alg, alg_name in algos:\n", - " eta = integrator(x, y0, ODEF, alg)\n", - " ax0.plot(x, eta, label=alg_name)\n", - " ax1.plot(x, np.abs(eta - y), label=alg_name)\n", + " for i in range(len(h_list)):\n", + " x = np.arange(0, 2, h_list[i])\n", + " y = np.zeros(len(x))\n", + " y[0] = y0\n", + " \n", + " eta = integrator(x, y0, ODEF, alg)\n", + " y2[i] = eta[-1]\n", + " plt.loglog(h_list, y2 - y[-1], label=alg_name)\n", "\n", - "ax0.legend()\n", - "plt.xlabel(\"x\")\n", - "ax0.set_ylabel(\"y(x)\")\n", - "ax1.legend()\n", + "plt.xlabel(\"h\")\n", + "plt.legend()\n", "plt.show()" ] },