From debd58eb5602182b465a963425a72b9b71cca9fe Mon Sep 17 00:00:00 2001 From: Kees van Kempen Date: Tue, 8 Mar 2022 16:37:05 +0100 Subject: [PATCH] 09: Do task 5. Looks good except for AB4. --- .../9 Ordinary Differential Equations.ipynb | 77 ++++++++++++++++--- 1 file changed, 68 insertions(+), 9 deletions(-) diff --git a/Week 5/9 Ordinary Differential Equations.ipynb b/Week 5/9 Ordinary Differential Equations.ipynb index 242a58b..aaa82c5 100644 --- a/Week 5/9 Ordinary Differential Equations.ipynb +++ b/Week 5/9 Ordinary Differential Equations.ipynb @@ -472,7 +472,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 7, "metadata": { "deletable": false, "nbgrader": { @@ -490,12 +490,20 @@ "outputs": [], "source": [ "def phi_ab3(x, y, f, i):\n", - " # YOUR CODE HERE\n", - " raise NotImplementedError()\n", + " # The first three y values are to be calculated using the Runga-Kutta method.\n", + " if i < 3:\n", + " return phi_rk4(x, y, f, i)\n", + " \n", + " phi = (1/12)*( 23*f(x[i - 1], y[i - 1]) - 16*f(x[i - 2], y[i - 2]) + 5*f(x[i - 3], y[i - 3]) )\n", + " return phi\n", "\n", "def phi_ab4(x, y, f, i):\n", - " # YOUR CODE HERE\n", - " raise NotImplementedError()" + " # The first three y values are to be calculated using the Runga-Kutta method.\n", + " if i < 4:\n", + " return phi_rk4(x, y, f, i)\n", + " \n", + " phi = (1/24)*( 55*f(x[i - 1], y[i - 1]) - 59*f(x[i - 2], y[i - 2]) + 37*f(x[i - 3], y[i - 3]) - 9*f(x[i - 4], x[i - 4]) )\n", + " return phi" ] }, { @@ -522,7 +530,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 8, "metadata": { "deletable": false, "nbgrader": { @@ -537,10 +545,61 @@ "task": false } }, - "outputs": [], + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtMAAAH0CAYAAAD/taEiAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAACicklEQVR4nOzdd3yV9fn/8dcn55zkZE/2CpuwRxCU4aAyFHBvUWtdVVvqt1ptawu21trqz7rrqBYXuIuIGwUVUCDBoECQvVf2HifnfH5/nBDDCISQcDLez8cjD3Lu+z73fd137oQrn1z39THWWkRERERE5PgFBToAEREREZGmSsm0iIiIiEgdKZkWEREREakjJdMiIiIiInWkZFpEREREpI6UTIuIiIiI1JGSaRERERGROlIyLSIijZ4x5hRjzDfGmK+MMXOMMa5AxyQiAkqmRUSkadgBnGWtHQtsBc4LbDgiIn5KpkWkWTPGbDXG/OwE3r/GGHNG/UUUWMaYycaY8EDHcbystXustSWVL8sBXyDjERE5QMm0iDQJxpjRxpilxpg8Y0y2MWaJMWZ4PR/jsMTbWtvPWruoPo9zvDHUs0cBZwPuv0EZY7oA44H362FftxtjUowxZcaYWbXYPs4Y8z9jTJExZpsx5soTjUFEmr4m+wNVRFoOY0wUMB/4JfAmEAyMAcoCGVdTY4zpA3xtrc0LdCx1UXkfvAJcZ631HGPbmQDW2plH2Ww3cD8wAQitRQhP4R8VbwMMBj4wxqyy1q6pxXtFpJnSyLSINAW9AKy1c6y1XmttibX2U2vt9wDGmCRjzCJjTG5lWcbUmnZkjLHGmB7VXs8yxtxvjHkF6Ay8b4wpNMb8rnJ91UjxsY5Tue2dxpjvK0fQ3zDGuCvX3W2M2WWMKTDG/GiMGXeE2GqKoVbnZ4yJMMZ4jTHtqi3rb4zZY4yJBM4Cnq3NBT/K9TvWMU5k3/80xsyt9vohY8znxphgY4wTeB24z1r744kc5wBr7bvW2rlAVi1iCwcuAv5krS201i4G5gHT6iMWEWm6lEyLSFOwHvAaY14yxkwyxsQeWFHZ1eF94FOgNfAr4DVjTO/jOYC1dhqwHZhirY2w1v6z+vrjOM6lwESgKzAQuK5ym9uB4dbaSPwjoVtrE8PxnJ+1thBYBwyttvhB4AFrbQEQaa399tD3GWPmVybqR/qYf5zHOBH/AM40xgwxxtyC/zpeaK0tB64ARgB/qvzF4rITPNbx6gVUWGvXV1u2Cuh3kuMQkUZGybSINHrW2nxgNGCB54EMY8w8Y0wbYCQQATxorS231n6BvyTkinoOo7bHedxau9tam40/CR4MeIEQoK8xxmWt3Wqt3VTPxz1gBZWJrjFmLNCXytFoa+0/jvQGa+1ka21MDR+Tj+cYJ8JamwX8C3gJ+D1wzoGSFGvtK9baeGvtGZUfb5zo8Y5TBJB/yLI84IRG40Wk6VMyLSJNgrU23Vp7nbW2I9AfaI//Ybr2wA5rbfXuDtuADvUcQm2Ps7fa58VAhLV2I/AbYCaw3xjzujGmfT0f94CqRBf4J/6yhPJaHqu2anWMyhFkW8PH4hr2/R0wAPi9tXbH8QRVfYQduAe4p6YR9jooBKIOWRYFnOhovIg0cUqmRaTJsdauA2bhT6p3A52MMdV/nnUGdtXw9mIgrNrrttV3fZTDHu9xDo15trV2NNCl8jhHHCU+QgzHe9wVwFBjzEWAG5h9rNiMMR9V1mgf6eOjuh6jcgTZ1PAx+ghxDAD+jX9k+vpjxX2E41WNsOMvPXnwGCPsx2M94DTG9Ky2bBCghw9FWjgl0yLS6Blj+hhjfmuM6Vj5uhP+ModvgWX4E+TfGWNcxt8Tegr+h9WOJA240hjjMMZMBE6vtm4f0K2G9x3vcarH39sYc5YxJgQoBUqouU/yoTEc73FX4f8F4f/hH9092i8IAFhrJ1XWaB/pY1J9HONYjDEd8JfF3ALcCgwwDdzf2xjjrHxA1AE4jDHuygcdD2OtLQLeBf5ijAk3xozCP3HMKw0Zo4g0fkqmRaQpKMD/8NkyY0wR/iR6NfDbyvKCKcAkIBN4GrimcvT6SKZXbp8LXAXMrbbu78C9lWUBd1Z/Ux2OU10I/pHSTPxlIK3x1wQfyUExHO9xrbVlwA/AVmvtkUaVT1h9H8P4W959CDxirZ1nrS0GHgL+dqL7PoZ78f9icw9wdeXn91aL6yNjzB+qbX8r/hZ6+4E5wC/VFk9ETD0MKIiISCNhjAkGNgKXHqlzR1M5hohIU6GRaRGR5mUGsKSBk9yTcQwRkSZBybSISDNgjBlqjMkDxuLvRd0kjyEi0tSozENEREREpI40Mi0iIiIiUkdKpkVERERE6qjRJNPGmG7GmBeMMW8HOhYRERERkdpo0JppY8yLwGRgv7W2f7XlE4HH8DfK/4+19sFq69621l5cm/0nJCTYxMTE+g1aREREROQQqampmdbaVocuP+JMT/VoFvAk8PKBBcYYB/AUcDawE1hhjJlnrV17vDtPTEwkJSWlnkIVERERETkyY8y2Iy1v0DIPa+1XQPYhi08BNlprN1fO7PU6/ilZRURERESalEDUTHcAdlR7vRPoYIyJN8Y8AwwxxtQ0zS7GmJuMMSnGmJSMjIyGjlVEREREpEYNXeZRa9baLOCWWmz3HPAcQHJysppki4iIiEjABCKZ3gV0qva6Y+WyeuHxeNi5cyelpaX1tUuRKm63m44dO+JyuQIdioiIiDQCgUimVwA9jTFd8SfRlwNX1tfOd+7cSWRkJImJiRhj6mu3IlhrycrKYufOnXTt2jXQ4YiIiEgj0KA108aYOcA3QG9jzE5jzC+stRXA7cAnQDrwprV2TX0ds7S0lPj4eCXSUu+MMcTHx+uvHiIiIlKlQUemrbVX1LD8Q+DDhjquEmlpKLq3REREpLpGMwOi1M7cuXNZu/a4W3IfZtGiRUyePPmo2+Tm5vL0009Xvd69ezcXX1yr+XREREREWgQl001MfSXTtXFoMt2+fXvefluzvYuIiEhgNOTM3XWlZLoBvPrqq5xyyikMHjyYm2++Ga/Xy4oVKxg4cCClpaUUFRXRr18/Vq9eTWFhIePGjWPo0KEMGDCA9957r2o/L7/8MgMHDmTQoEFMmzaNpUuXMm/ePO666y4GDx7Mpk2bDjruW2+9Rf/+/Rk0aBBjx44F/DXkP//5zxkwYABDhgxh4cKFh8U7c+ZMHn744arX/fv3Z+vWrdxzzz1s2rSJwYMHc9ddd7F161b69+9/1P3OmjWLCy+8kIkTJ9KzZ09+97vf1fv1FRERkZbB+izluwsp+HonmbPWsP+J7wId0mEaTZ/p42GMmQJM6dGjx1G3u+/9NazdnV+vx+7bPooZU/rVuD49PZ033niDJUuW4HK5uPXWW3nttde45pprmDp1Kvfeey8lJSVcffXV9O/fn4qKCv73v/8RFRVFZmYmI0eOZOrUqaxdu5b777+fpUuXkpCQQHZ2NnFxcUydOpXJkycfsdziL3/5C5988gkdOnQgNzcXgKeeegpjDD/88APr1q1j/PjxrF+/vlbn+uCDD7J69WrS0tIA2Lp1a9W6o+03LS2N7777jpCQEHr37s2vfvUrOnXqdIQjiIiIiPzEWkvF/mLKNuVRtimXsi15+IorAHAmhBLSPRpb4cM4G894cJNMpq217wPvJycn3xjoWA71+eefk5qayvDhwwEoKSmhdevWAPz5z39m+PDhuN1uHn/8ccB/0/zhD3/gq6++IigoiF27drFv3z6++OILLrnkEhISEgCIi4s75rFHjRrFddddx6WXXsqFF14IwOLFi/nVr34FQJ8+fejSpUutk+mjOdp+x40bR3R0NAB9+/Zl27ZtSqZFRETkMNZavFmllG7KpWyzP4H2FXoAcMSE4E6KJ6R7NCHdY3BGhwQ42iNrksl0bR1tBLmhWGu59tpr+fvf/37YuqysLAoLC/F4PJSWlhIeHs5rr71GRkYGqampuFwuEhMT69x67ZlnnmHZsmV88MEHDBs2jNTU1Fq9z+l04vP5ql6faOu3kJCfbnaHw0FFRcUJ7U9ERESaj4rc0p9Gnjfl4c0rAyAoMhh3jxhCuscQ0i0aZ3xogCOtncYzRt5MjBs3jrfffpv9+/cDkJ2dzbZt2wC4+eab+etf/8pVV13F3XffDUBeXh6tW7fG5XKxcOHCqm3POuss3nrrLbKysqr2AxAZGUlBQcERj71p0yZGjBjBX/7yF1q1asWOHTsYM2YMr732GgDr169n+/bt9O7d+6D3JSYmsnLlSgBWrlzJli1bjnms2uxXRERExFtQTnHafnLe3cCeh1aw98EV5Ly1ntIfswnuHEnM+d1p89thtPvDKcRd3ofw4W2bTCINzXxkOhD69u3L/fffz/jx4/H5fLhcLp566im+/PJLXC4XV155JV6vl9NOO40vvviCq666iilTpjBgwACSk5Pp06cPAP369eOPf/wjp59+Og6HgyFDhjBr1iwuv/xybrzxRh5//HHefvttunfvXnXsu+66iw0bNmCtZdy4cQwaNIg+ffrwy1/+kgEDBuB0Opk1a9ZBI8cAF110ES+//DL9+vVjxIgR9OrVC4D4+HhGjRpF//79mTRpErfddlvVe2699dZj7ldERERaHl+xh7LNef7SjU15VOwvBsC4HYR0jSbi1PaEdI/B1SYME9T0528wjbHFSG0lJyfblJSUg5alp6eTlJQUoIikJdA9JiIi8hNfWQVlW/OryjY8uwvBgnEFEdw1mpBu0bh7xOBqH9Gkk2djTKq1NvnQ5RqZFhEREZFasx4fZdt/Sp7LdxSAz4LDENIliqifdSGkezTBHSMbVdeNhqJkWkRERERqZL0+yncWVibPuZRty4cKC0EQ3DGSyNM7+jtudInCuByBDvekUzItIiIiIlWsz+LZW+RPnDfmUrYlH1vuBcDVLpyIke0J6RFDSGIUQW6lkroCIiIiIi2YtZaKzJKqso2yTbkHTZQSNrS1f+S5WwyOcFeAo218mmQyXdsZEEVERETkcBW5ZT+NPG/KxZtfDoAjOhh3nzj/yHMjniilMWmSyXRjngFRREREpLHxFnl+qnnelEdFZgkAQeFOQrrFENIjBnf3GBzxboxpuh03AqFJJtONncPhYMCAAVWvL7/8cu65554at581axYpKSk8+eSTJyM8ERERaeZ8ZV7KtuRVjTx79hQBYEL8vZ7DR7QjpEfz6fUcSEqmG0BoaChpaWkNtv+KigqcTn3pRERExM9W+Cjfnk/pxkPa1TkNIZ2jiBrfhZDuMf52dQ4lz/VJGdlJlJiYSEpKCgkJCaSkpHDnnXeyaNGig7bJyMjglltuYfv27QA8+uijjBo1ipkzZ7Jp0yY2b95M586dmTNnTgDOQERERBoD67N4dvvb1ZVuzKV8az7W4wMDro6RRI6tbFeX2DLb1Z1MzTuZ/uge2PtD/e6z7QCY9OBRNykpKWHw4MFVr3//+99z2WWX1Wr306dP54477mD06NFs376dCRMmkJ6eDsDatWtZvHgxoaFNZ756EREROXEHddzYkEvp5jxsSWXHjdZhhA9vS0j3GEK6RRMU2rzTu8ZGV7sBnEiZx4IFC1i7dm3V6/z8fAoLCwGYOnWqEmkREZEWwptfTumBjhsbc/DmHei4EUJo33jclR03HFHBAY60ZWveyfQxRpBPNqfTic/nA6C0tPSI2/h8Pr799lvcbvdh68LDwxs0PhEREQkcX2lFVZ/n0o25VOwvBiAozOkfde7u77rhVMeNRqV5J9ONTGJiIqmpqUyaNIl33nnniNuMHz+eJ554grvuuguAtLS0g0pGREREpHmwFT7KtuVXjjznUr6zACwYVxDBXaMJH9bG33GjXbg6bjRiSqYbwKE10xMnTuTBBx9kxowZ/OIXv+BPf/oTZ5xxxhHf+/jjj3PbbbcxcOBAKioqGDt2LM8888zJCVxEREQajPVZPHuKKNuYc/BDg0EQ3DGSyDM74e4RQ3DnKIwzKNDhSi0Za22gY6iz5ORkm5KSctCy9PR0kpKSAhSRtAS6x0REpLYqskr87eoq+z1XTdPdOsxf89yj8qFBt8Y3GztjTKq1NvnQ5U3yK6fpxEVERKQxqpppcKO/7tmb7X9GyhH10zTd7h4xOKI0TXdz0SSTaU0nLiIiIo2B9Xgp21I5WcrGHDy7K2cadDsI6RZD5OgO/ocGW4XqocFmqkkm0yIiIiKBYH0Wz67CquS5bFs+VFhwGII7RxF1dhdCesYQ3EEzDbYUSqZFREREjqKq7nlDDqWbfposxdU2nIiR7XH3jCG4azRBwZppsCVSMi0iIiJSja/YU/XQ4EF1z9HB/slSelZOlhKpyVJEybSIiIi0cFX9njfkUroxB8+uQn+/5xAHId2i/XXPPWNwJqjuWQ6nZFpERERaFGstFfuKKd2QQ+mGXMq35P3U77lTFFHjOhPSM5bgjqp7lmNTMt3IJSYmkpKSQkJCAqeddhpLly4F4K677uLDDz/knHPOoXv37oSFhXHNNdfUer8REREUFhYettzhcDBgwICq15dffjn33HNPjfuZNWsWKSkpPPnkk8dxVgfbs2cPN954I/Pnzwdg+fLl3Hnnnezbt4+wsDCGDRvG448/TlhY2DFjmDlzJhEREdx55501Hm/u3Ln06tWLvn371rjNk08+SVhYGNdff32dz0tERBoPb345pRtyKks3cvAVeABwtgolLLkN7p6x6vcsdaI7pgk5kEgDPPfcc2RnZ+Nw1O/DDqGhoaSlpdXrPqurqKjA6Tz4tnvkkUe48UZ/l8N9+/ZxySWX8Prrr3PqqacC8Pbbb1NQUFBjMn285s6dy+TJk4+aTF9//fWMGjVKybSISBPlK/dStiXPX7qxIYeKfcUABIU7CekR658wpWcszhj1e5YT06yT6X8s/wfrstfV6z77xPXh7lPurnH91q1bmThxIiNHjmTp0qUMHz6cn//858yYMYP9+/fz2muvccopp5Cdnc3111/P5s2bCQsL47nnnmPgwIFkZWVxxRVXsGvXLk499VSqz1B5YDR56tSpFBYWMmzYMH7/+9+Tnp5eNRq7adMmbrvtNjIyMggLC+P555+nT58+bNmyhSuvvJLCwkLOO++84z7v6iPkKSkp3HnnnSxatOigbTIyMrjlllvYvn07AI8++iijRo1i5syZbNq0ic2bN9O5c2fmzJlz0Pveeecd7r//fgCeeuoprr322qpEGuDiiy8GqPGa1eT555/nueeeo7y8nB49evDKK6+QlpbGvHnz+PLLL7n//vt56qmnuO2226re88MPP7B582a6dOlCYmIiy5cv55RTTjnu6yUiIifXgam6SzfkULYhh7Kt+eC14DSEJEYTPrQ1IT1icbULxwSpdEPqjyZ+bwAbN27kt7/9LevWrWPdunXMnj2bxYsX8/DDD/PAAw8AMGPGDIYMGcL333/PAw88UFWicd999zF69GjWrFnDBRdcUJWYVjdv3ryqEeTLLrvsoHU33XQTTzzxBKmpqTz88MPceuutAEyfPp1f/vKX/PDDD7Rr167G2EtKShg8eHDVxxtvvFHr854+fTp33HEHK1as4J133uGGG26oWrd27VoWLFhwWCK9ZcsWYmNjCQnxjwysXr2aYcOGHXH/NV2zmlx44YWsWLGCVatWkZSUxAsvvMBpp53G1KlTeeihh0hLS2PUqFGkpaWRlpbGjTfeyEUXXUSXLl0ASE5O5uuvv671+YuIyMnlzSujKGUfWXPWsedv37L/ie/I/3grvqIKIk5rT8L1/ekw41Ra3TCAyNM7EdwhQom01LtmPTJ9tBHkhtS1a9equuN+/foxbtw4jDEMGDCArVu3ArB48WLeeecdAM466yyysrLIz8/nq6++4t133wXg3HPPJTY2ttbHLSwsZOnSpVxyySVVy8rKygBYsmRJ1fGmTZvG3Xcf+dqcSJnHggULWLt2bdXr/Pz8qrrsqVOnEhoaeth79uzZQ6tWrWq1/5quWU1Wr17NvffeS25uLoWFhUyYMKHGbZcsWcLzzz/P4sWLq5a1bt2adevq9y8bIiJSd1WlG+v9Dw5W7K8s3Yhw4e4VR0gvf/mGWtbJydQkk2ljzBRgSo8ePQIdyhEdGGUFCAoKqnodFBRERUVFgx3X5/MRExNTYzJ8Iu18nE4nPp8PgNLS0hqP/+233+J2uw9bFx4efsT3hIaGHrS/fv36kZqaWqdSlENdd911zJ07l0GDBjFr1qzDylIO2LNnD7/4xS+YN28eERERVctLS0uP+AuAiIicHNZn8ewt8k+Wsr566UYQIV2jCE9ug7tXLM42YWpZJwHTJMs8rLXvW2tvio6ODnQodTZmzBhee+01ABYtWkRCQgJRUVGMHTuW2bNnA/DRRx+Rk5NT631GRUXRtWtX3nrrLcDf+mfVqlUAjBo1itdffx2g6rjHIzExkdTUVICq0eFDjR8/nieeeKLqdW1GuHv16lU1Wg9w++2389JLL7Fs2bKqZe+++y779u2r8ZrVpKCggHbt2uHxeA4658jISAoKCgDweDxccskl/OMf/6BXr14HvX/9+vX079//mOcgIiL1x1tQTtHKfWS/vo49Dyxj/+PfkfdRtdKNX/Snw4yRtPrFACLHdsTVNlyJtARUk0ymm4OZM2eSmprKwIEDueeee3jppZcAf13wV199Rb9+/Xj33Xfp3Lnzce33tdde44UXXmDQoEH069eP9957D4DHHnuMp556igEDBrBr164a339ozfSBtngzZsxg+vTpJCcn19hB5PHHHyclJYWBAwfSt29fnnnmmWPGGx4eTvfu3dm4cSMAbdq04fXXX+fOO++kd+/eJCUl8cknnxAZGVnjNavJX//6V0aMGMGoUaPo06dP1fLLL7+chx56iCFDhrB06VJSUlKYMWNG1Tnv3r0b8Jd+nH322cc8BxERqTvr8VG6MYfcD7ew79GV7PnbMnLeXE/phlzcPWKIvaQX7f4wgja/GUrMud1w94zFuDRttzQepnq3iKYmOTnZpqSkHLQsPT2dpKSkAEUkdfG///2P1NTUqo4ejcF3333HI488wiuvvHLYOt1jIiJ1Z62lIqPEX7axIYeyzZUTpjgMIV2i/HXPPdV1QxofY0yqtTb50OVNsmZampcLLriArKysQIdxkMzMTP76178GOgwRkWbBV+yhdFMuZetzKV2fgzfP/3C8MyGU8OFtCekVS0jXaIJCNOIsTY+SaWkUqrfRawxU3iEiUnfWZynfWeDvurE+h/IdBWDBuB24u8cQclYn3D1jccYd/sC6SFOjZFpEREROmDevjNLKrhulG3KxJRVgwNUxksgzO+HuFUtwpyiMQ6Ub0rwomRYREZHjZit8lG3N89c+r8/Bs7ey53NkMKF943H3iiGkRyyOcFeAIxVpWEqmRUREpFYqMv0PDpauz6FsU+5PDw4mRhE9KZGQXnG42qrns7QsSqYbQERERNXMfwCzZs0iJSWFJ598MoBRiYiIHB9fuZeyzXmU/pjtf3Awyz/JliPeTVjlhCkh3WL04KC0aEqmRUREBKhsW7e/2D/6/GMOZVvywGsxriBCuscQOaqDf8bBBM0OKydHRUUF6at+YO2yFLK3bsdbXMrtTz4U6LAOomT6JMvIyOCWW25h+/btADz66KOMGjWKmTNnEhERwZ133glA//79mT9/PgCTJk1i9OjRLF26lA4dOvDee+9pmmsREakXvtIKyjblUvqjv3zDm1vZtq51GBGntfePPidGY1ya503qn8fjYeW3y9mUmkbu9l1U5BZgy8qwPg/WluGzxUBFtXc48Xg8uFyNpxa/WSfTex94gLL0dfW6z5CkPrT9wx+Ous2BWQQPyM7OZurUqQBMnz6dO+64g9GjR7N9+3YmTJhAenr6Ufe3YcMG5syZw/PPP8+ll17KO++8w9VXX33C5yIiIi2PtRbP3mJKf8ymbH0OZVvzwWcxIQ5CesQQeVYn3L3icMaEBDpUaQY8Hg+rv1vFumUp5GzZgSc3Hw5KlosA7yHvCiHIhOIwwYQEBRPsKCfCVUhCaBaJCR5oZBMONutkOlBCQ0NJS0uren2gZhpgwYIFrF27tmpdfn7+QfXVR9K1a9eq5HzYsGFs3bq1vkMWEZFmzFdaQemG3KoE2ptfDoCrXTiRYzrg7h1LcJcojEOjz3L8dm7bzoovvmT/j5soy8rBlpRivdWT5YpD3hGCozJZdjtcBAeVExVcQOuwLHq28dC+Y1uIbg0xnSC600//RncEZ+P7Ja9ZJ9PHGkEOBJ/Px7fffovbfXCjeqfTic/nq3pdWlpa9XlIyE83jsPhoKSkpOEDFRGRJqv66HPpj9mUbyvwjz67Hbh7xuLuFYu7dyyOqMaXmEjjU1ZaRuripWxIWUnBzr14C4qwnnLwleOjBGsPzUtcBJkwHMZFSFAMwQ5/stwqNJNebSro0LktxLSFmC4Q0/mnj8h2ENT0HmZt1sl0YzR+/HieeOIJ7rrrLgDS0tIYPHgwiYmJVTXSK1euZMuWLYEMU0REmhhfWQVlGyprn3/MPnj0eWxH/+hzZ02aIkeWmZHBtwsWsmf1Okr2Z1WOLpfXUIphMCYcByEEm3CCnSGEBxcRH5JD94RiEjvH4IpvDbFd/AnzgX8j24Gj+aWeTfKMjDFTgCk9evQIdCjH7fHHH+e2225j4MCBVFRUMHbsWJ555hkuuugiXn75Zfr168eIESPo1atXoEMVEZFGrKrzRmXyXLY13995I8SBu2cM7t5xGn2Wg+zctp3lCxay78eNeDJzKx/0K8dnS7C26JCtnQSZcH8pRpALt7OU6JB8OkRm06+zi4h20RDbCWITISax8t9OjbIMo6EZ28iKuI9HcnKyPVCLfEB6ejpJSUkBikhaAt1jIhIovnLvT5031mVXdd5wtQ0jpHdcZecN1T63ZDu3bWfZZ1+wf91GPNmVCbO3HB/FRyjHcFfWLjtxBVUQ5iomwZ1Nt/hiuneNwZXQA+ISIbar/yOuK4TGQgudlMcYk2qtTT50eZMcmRYREWkpKrJKKF2XTcmPOZRtzoUKiwkOIqRHLJFndsLdW503WprMjAyWfPwZe1evozQzB0pLKxPmEqwtPmhbY0IJIpRgE0Gw0014cAFtQrPp2bacLl07QlzbykS5mz9ZjukCLncNR5YjUTItIiLSiNgKH2Vb8ihdl03pjzlUZPpHE52tQokY2R53n8q+z06NPjdnxSUlfLNgIVtWfEfx3gxscYk/Ya4qyaheWeDGYcIINuG4nCFEVCbMfTpW0LFLJ4jr6E+U47r5P6I7Ncva5UDRlRQREQkwb34ZpetyKPkxm7INudhyLzgNId1iiDi1He4+cTjjNVlXc5T+w1pWLviC3E3bqSgowJaXY20pPlvIwQ/9uXCYcFwmFJfDRbirmDZhmfRuV06Xbp0gvoM/UY7vDnHd/W3kmmBnjKZIybSIiMhJZn2W8p0F/tHnddl4dvsf/nJEBxM2pBXu3nGE9IghKFjJUHNQVFzMko8/Y1vqKkr2ZWJLKkeZKcLa0mpbBlU+9BdCiCOGUEcJ8aE59GhdSM+uCbjatIH4Hv6PuG7+kgyNMAecvgIiIiInga+0gtL1OVXlG74iDxgI7hxF1IRE3H3icLUNw7TQh7uag13btrP4g0/J/HEjFbn52PKyGkaZQ3CYcIJNBCGuEKJD8ukYlc2A7qFEtu/gH11O6PlT0uzSXyUaMyXTIiIiDcSTUexPntOzq6btDgpzEtIrltA+cYT0jMUR7gp0mHKcVq9MY+Xni8jbvANfQWFlLXPxIe3lDEEmAocJwe2IIcxZTEJoFn3al9Gte0dMqy6VCXNP/79h8S22S0ZTp2S6gXi9XpKTk+nQoUPVZCyPPvooN910E2FhYQBEREQccyrxukpMTCQlJYWEhIQ6vf+6667jyy+/JDo6GmstjzzyCOPGjavnKGvnwHTsTz75JD6fj5///Oc4HA5eeOGFI47g5ObmMnv2bG699VYAtm7dytKlS7nyyitPdugi0sLYCh9lW/MoTfeXb1Rk+f+E72wTRuTYDrj7xBHcSROnNAUVFRWkfrOMtYuWULBjN7aoGOstO0JphpMgE4HLhBHsDCYyOJ8OkVkM6h5MTKeO/kQ5oZf/IzYRnMGBOiVpIEqmG8hjjz1GUlIS+fn5VcseffRRrr766qpkur5UVFTgdJ7Yl9Lr9eJwHFyb99BDD3HxxRezcOFCbrrpJjZs2HBCxzhR1lpuueUWPB4P//3vf2v8U2hubi5PP/30Qcn07NmzlUyLSIPwFpZX9X0uXZ+DLfM/POjuHkPEaH8C7YxVq7HGyuPxkPL1UtZ+tZSiXXuxxf6k2WuLgLJqWwbjMBEEm0hCXCHEhuTQJT6PQb1iCW7X7qeEuVVviGijUeYWRMl0A9i5cycffPABf/zjH3nkkUcA/8yHu3fv5swzzyQhIYGFCxcC8Mc//pH58+cTGhrKe++9R5s2bQ7aV3Z2Ntdffz2bN28mLCyM5557joEDBzJz5kw2bdrE5s2b6dy5M08++SRXXHEFu3bt4tRTT6X6ZDyvvvoqjz/+OOXl5YwYMYKnn34ah8NBREQEN998MwsWLOCpp55i9OjRRzyfU089lV27dgEHjxIDTJ48mTvvvJMzzjiDiIgIpk+fftj5bNq0iauuuoqioiLOO+88Hn300aoR+Yceeog333yTsrIyLrjgAu67774ar+uvf/1rsrKyeOONNwgKCmLmzJlERERw5513AtC/f3/mz5/PPffcw6ZNmxg8eDBnn302X3/9Nenp6QwePJhrr72WCy64gGnTplFU5P9z3JNPPslpp5123F9nEWmZrLVU7CumJD2b0vQsyncUgIWgyGDCBrXC3UcPDzZG/pHmb1mzcAmFO3ZXS5oLgfJqW4bgNOG4gyJxO4OJc2fRo00RfXq1xdUuERJ6+xPmhJ7gjg7Q2Uhj0qyT6a/fXE/mjvoto0joFMGYS48+1fdvfvMb/vnPf1JQUFC17Ne//jWPPPIICxcurCq9KCoqYuTIkfztb3/jd7/7Hc8//zz33nvvQfuaMWMGQ4YMYe7cuXzxxRdcc801pKWlAbB27VoWL15MaGgov/71rxk9ejR//vOf+eCDD3jhhRcA/2x9b7zxBkuWLMHlcnHrrbfy2muvcc0111BUVMSIESP4f//v/x31fD7++GPOP//8Y16bms5n+vTpTJ8+nSuuuIJnnnmmavtPP/2UDRs2sHz5cqy1TJ06la+++oqxY8cetu/Zs2eTlJTEokWLjjkK/+CDD7J69eqq67Ro0SIefvjhqnKb4uJiPvvsM9xuNxs2bOCKK67g0Jk0RUSqsxU+yjbnUZKe5Z95MKdy5sEOEUSN6+x/eLB9BCZIo5GNwdq0Vaz4+HPytuzAFhZVS5qrjzQfSJqjcDvLiA/LolfbUnr2bI+rXU9o1eenpFkPAMpRNOtkOhDmz59P69atGTZsGIsWLTrqtsHBwUyePBmAYcOG8dlnnx22zeLFi3nnnXcAOOuss8jKyqoqHZk6dSqhof5v8K+++op3330XgHPPPZfY2FgAPv/8c1JTUxk+fDgAJSUltG7dGgCHw8FFF11UY3x33XUXf/jDH9i5cyfffPPNMc+9pvP55ptvmDt3LgBXXnll1Ujyp59+yqeffsqQIUMAKCwsZMOGDUdMpocOHcq6detYvnw5o0aNOmYsR+PxeLj99ttJS0vD4XCwfv36E9qfiDRP3iJP5cODWZSu9/d+Nq4gQnrEEHlmJ0L7xOGI0syDgbR7xy6+fm8+Ges24s0rwFaU4bNFh8wC6C/PCAmKItRZSnxoFr3aldCzZwdc7XtVJs19/EmzU19POX7NOpk+1ghyQ1iyZAnz5s3jww8/pLS0lPz8fK6++mpeffXVw7Z1uVxVdb8Oh4OKiorjOlZ4ePgxt7HWcu211/L3v//9sHVut/uwOunqDtRMP/HEE1x//fWkpqbidDrx+XxV25SW/vQQxvGej7WW3//+99x8880HLX/qqad4/vnnAfjwww8B6NOnD3/5y1+49NJL+eSTT+jXr99RYzmaf/3rX7Rp04ZVq1bh8/lwu1XLKCKV5RsZJZSmZ1GyNpvy7fn+8o2oyt7PfeJw94jBuFS+cbIVFRfz5fsfsW1FGuUZWdjy0sqZAAuqbeXAYSIJNuGEOF3EunPo3rqA/n3a4GrfDVr3rUyae2m6bKlXzTqZDoS///3vVYnrgfKCA4l0ZGQkBQUFx9VhY8yYMbz22mv86U9/YtGiRSQkJBAVFXXYdmPHjmX27Nnce++9fPTRR+Tk5AAwbtw4zjvvPO644w5at25NdnY2BQUFdOnSpdYx3H777bz44ot88sknJCYm8vTTT+Pz+di1axfLly8/5vtHjhzJO++8w2WXXcbrr79etXzChAn86U9/4qqrriIiIoJdu3bhcrm47bbbuO222w7bz2mnnca///1vJk+ezJdffkliYmJV6cbKlSvZsmUL8NN1PuDQ13l5eXTs2JGgoCBeeuklvF4vItIyWa/9qftGelZV9w1Xhwgiz+pMaFIcrg4R6v18Eq1OTWPFR59RsG0ntrgYn7cUny0AfhqgCTKROI2bYIeT6JB8OsXkMrRvNGEdO0ObvtA6yV/bHBIRuBORFkPJ9El00003MXHiRNq3b1/1AOKxzJw5k+uvv56BAwcSFhbGSy+9dMTtZsyYwRVXXEG/fv047bTT6Ny5MwB9+/bl/vvvZ/z48fh8PlwuF0899dRxJdPGGO69917++c9/smDBArp27Urfvn1JSkpi6NChx3z/gS4mf/vb35g4cSLR0f4HNsaPH096ejqnnnoq4G8V+Oqrr1aVoRzJlClTyMzMZOLEiXz22We8/PLL9OvXjxEjRtCrl/8vEfHx8YwaNYr+/fszadIkHnjgARwOB4MGDeK6667j1ltv5aKLLuLll19m4sSJtRrhF5Hmo2rylLVZlPyYgy2pAIchpHsMEWM64k6KwxmtP/c3tIL8Ar743zx2pa3Gk5WL9ZQeVqJhjBsH4YQ4ogl3FdM2IoshPR207t4dWverTJz7QlhcAM9EWjpTvetDU5OcnGwPfXAsPT2dpKSkAEUkR1JcXExoaCjGGF5//XXmzJnDe++9F+iw6kz3mEjTU5FbSml6NiVrsyjbnAdeS1C4E3efeEKT/JOnBIWofKOhpH+/hm/nf0TBlh3YomJ8vpLK0eYDfxkMqhxtDsbtLCfOnU1Sh2L6JCUS1K4/tOnvT5pjOqvlnASMMSbVWpt86HKNTEuDS01N5fbbb8daS0xMDC+++GKgQxKRZs5ai2d3ESVrsyhNz8Kz298K09kqlIhRHQjtG0dw5yh136hnHo+HLz/8hA1LllG+LwvKS/HaooNmBjww2ux2xhDhKqRjdDbDkiKISuwCbfpB2/7+Eg3VNUsT0SSTaWPMFGBKjx49Ah2K1MKYMWNYtWpVoMMQkWbOeivb163NonRtNt68MjAQ3DmK6EldcfeNw9WqfifNasmyM7P5/J257PthLd7cAmxFKV5bwE89mw1BJhKXCSPE6SIuNJuk9iX0SkrA1XEQtBngT5w1wYk0cU0ymbbWvg+8n5ycfGOgYxERkcDxlVZQ+mOOP4H+MRtbWtm+rmcsUWf7+z87IjR984navHETi9+dR+6GbdiiosoyjXzgQEclJw4ThdsRRZirmHYRWQztHUzrHh39CXPbAdAqSaPN0iw1yWRaRERaLm9+OSXpWZSsyaJsU25V/XNovwRC+8YT0lOzD56IH1K/Y/n8TyjatgtbWozXV3xQC7qfyjRiiQrJp0tcDsMGtCa0cx9oO9CfPMd00WiztBhKpkVEpNHz7C/2jz6vqZy+G3DEu4k4rT2hfeMJ7qL657pYsfgb0j5eQMmuvdjSkiPUN4fjNKGEOJ3EunP8DwX2bY2r05DKxHmAOmlIi6dkWkREGh3rs5TvLPC3r1uTRUVGCQCujhFEje9CaL94nK3D1P/5OCz78mtWfbKQ0t378JWV4LOFB7Wh89c3hxLidJEQnsXARA89+nWDdoP8H637QrBqzkUOpWS6gcydO5cLLriA9PR0+vTpw9atW0lKSqJ3795YawkPD+e///0vvXv3Zvny5dx0002A/wn0mTNncsEFFwT4DERETq6qBwjXZFGyNgtffjkEQUi3GCJOa487KR5njPo/18aKxd+Q9tFnlOyqKXGOwmXCCXU5aRWexdCehk5JnX5KnFv1BocrgGcg0nQomW4gc+bMYfTo0cyZM4f77rsPgO7du5OWlgbAs88+ywMPPMBLL71E//79SUlJwel0smfPHgYNGsSUKVNwOvXlEZHmzVfupWx9jj+BTs/GllZgXEG4e8Xi7p9AaO9YgsKU1B3N6tQ0vp33AcXb9/hrnI+SOLcJz2RYbyftk7pAu8HQfjDE94Ag1ZiL1JWytQZQWFjI4sWLWbhwIVOmTKlKpqvLz88nNjYWgLCwn/5sVlpaqj9bikiz5iv2UJKe7X+AcEMO1uMjKMxJaN84Qvsn4O4Zg3EpuTuSzRs38eWb/6Ng0zZsSTFeXxHWFlat95dq+BPn1hFZJPdy0D4p0Z80txtcmTgHBSp8kWapWSfTC2c9x/5tm+t1n627dOPM62466jbvvfceEydOpFevXsTHx5Oamkp8fDybNm1i8ODBFBQUUFxczLJly6res2zZMq6//nq2bdvGK6+8olFpEWlWvPnllKzNomRNJmWb8sBncUQHE5bchtB+CYR0jcY4NJBQXWZGBp/NeYuMNRuwBYX4fEWV7ej8jInAady4nU5ah2cxpBd06dsF2g/xf8R314izyEmgjK0BzJkzh+nTpwNw+eWXM2fOHG6//faDyjzeeOMNbrrpJj7++GMARowYwZo1a0hPT+faa69l0qRJuN3qxykiTVdFdiklazIpWZ1F+fZ8sOBMCCVybAdC+yXg6hihv8RVKistY8G777H1m1QqcnLxVRTjs3kc6ONsjBuniSDMGUecO5tBiWX0GtwLOgz1J84JvcGh/9JFAqFZf+cdawS5IWRnZ/PFF1/www8/YIzB6/VijOG22247aLupU6fy85///LD3JyUlERERwerVq0lOPmz6dxGRRs2zv5iSHzIpWZNZNYW3q304UT/rQmh/deA4IHXpMlLe/5jS3Xux5SV4ffn8NHOgC6eJwu2MJTokj97tCxk4uAOuLsOh/VD/lNtOPYgp0lg062Q6EN5++22mTZvGs88+W7Xs9NNPZ8eOHQdtt3jxYrp37w7Ali1b6NSpE06nk23btrFu3ToSExNPZtgiInVircWzu4iS1ZmUrM6samEX3DmS6HO6+lvYxYcGOMrA2rd3L5+99ibZ6ZuwxYWH1DkbHCYatyOKCFcRXRKyOWVILGFdh0GHYf7OGiERAY1fRI5OyXQ9mzNnDnffffdByy666CL+/ve/V9VMW2sJDg7mP//5D+BPrB988EFcLhdBQUE8/fTTJCQkBCJ8EZFjstZSvqOgMoHOwptdCgZCukb7J1HpF48jqmWOnHo8Hha+/zEbv1qCJzPnCOUa/klQQl0O2kZkMWJAMK379PEnzh2GQUTrwJ6AiBw3Y60NdAx1lpycbFNSUg5alp6eTlJSUoAikpZA95i0RNZnKd+e7y/hWJ2JN68cHIaQ7jGE9U/A3TcOR0RwoMM86bZt3sIXs9+mYPM2fCXF+Gw+1pZWrnXiNFEEOwyxoTn071xK/2EDoEMydEyGuO7qrCHShBhjUq21h9XgamRaRESOyPosZVvyqkagfQXl4DS4e8YSNT6B0KS4FtUD2uPx8MV7H7Dxq2+oyMrB5z0w6uwflAoyUQSbSMJDHHSKzeHUoTGE9zjFnzi3GwSull3uItJcKZkWEZEq1msp25zrT6DXZOEr9PgnUekdS+iABNx94ggKaRn/dezZtZvPXn2d3B+3YEuK8B406ux/SDDUGUtcWDbDevnoPqg7dBzu/4hqF9DYReTkaRk/EUVEpEbW66NsU+UI9JpMfEUVmOAg3H3i/Al07ziCgpt/v+JvF37Jd/M/o3xfxmG1zkEmsmrUuUt8NiOHtyKs20h/4tx2gKbeFmnBlEyLiLRABxLo4u8zKF2bha+4AhPswJ0UR9iABNy9Y5v1LITFJSV89Mob7EpdhS0owOsrwNqiyrUOnCYatzOGuNBchvSooNfQHtDpFOh4CkS2CWjsItK4KJkWEWkhqkagK/tA+4orMCEOQpMqR6B7Nd8Eevf2HXzyyhvkb9yCLS3G68vjQF9nY0JxmnDCXA46xGQyalgronqPgE4joO1AcLa8BytFpPaUTIuINGNVNdCVXTgOjECH9o0jdECrygS6+XWU+G7ZCpa9+z5lu/djPcV4bS4/lWz4+zpHBhfQu0MBQ4d3x9VtlD95jukMmlRGRI6DkukGMnfuXC644ALS09Pp06cPW7duJSkpid69e2OtJTw8nP/+97/07t276j3bt2+nb9++zJw5kzvvvDOA0YtIU1bVheP7DH8CXVRZwtG3soSjmY1AV1RU8OUHH7NuwVdUZOfg8xbhs/mVa/0lG2HOWOLDshmeZOk6uC90HuGvd3ZHBzR2EWn6lEw3kDlz5jB69GjmzJnDfffdB0D37t1JS0sD4Nlnn+WBBx7gpZdeqnrP//3f/zFp0qRAhCsiTZz1Wcq35VP8fQYlP2T+1IUjKY6wga2aVQ10eWk5H73xJtu+ScXm51fWOxdXrg3GFRRFuDOG9lGZjEqOJDZpBHQe6S/Z0IOCIlLPlEw3gMLCQhYvXszChQuZMmVKVTJdXX5+PrGxsVWv586dS9euXQkPDz+ZoYpIE1Y1E+EqfwLtzS/3J9B94ggd2Hy6cBTkF/DBS7PZv2oNtqiQCl8eUAaAMWG4TBjhriASW+UwemR7gnuMgS6nQVw3lWyISINr1sl07vubKN9ddOwNj0Nw+3BipnQ/6jbvvfceEydOpFevXsTHx5Oamkp8fHzVdOIFBQUUFxezbNkywJ98/+Mf/+Czzz7j4Ycfrtd4RaR5sdbi2V3kH4H+PgNvThk4DO7ecUQPTMCdFE9QSNNOoLOzsvlw1qtkr1mPr6QIry8X8ABgTCQhQVFEhhTSu0Mew0b0wtV9DHQ+FSLbBjRuEWmZmmQybYyZAkzp0aNHoEM5ojlz5jB9+nQALr/8cubMmcPtt99+UJnHG2+8wU033cTHH3/MzJkzueOOO4iIiAhg1CLSmHn2F1O8KoOSVRlUZJZAkMHdM4aon3UhtF88Qe4m+eMcgMyMDD588RVy1m3Clh5InisA/6yCbkc00e58BnYtZeDIXtBllL9NXWhMIMMWEQHAWGsDHUOdJScn25SUlIOWpaenk5SUFKCIIDs7m44dO9KqVSuMMXi9XowxfPnll0yZMoXVq1cDUFJSQnx8PMXFxYwZM4YdO3YAkJubS1BQEH/5y1+4/fbbA3YeUrNA32PSclRkl1Yl0J69RWAgpFs0oQNbEdo/AUd406z/zczI4IP/vEzej5vxlRXh9eUAXsDfaSMkyEVMaC7Dennpfcpp/pKNjsMhOCywgYtIi2aMSbXWJh+6vOkOZTRSb7/9NtOmTePZZ5+tWnb66adXJcsHLF68mO7d/eUiX3/9ddXymTNnEhERoURapIXyFpT7SzjSMijfUQBAcOdIoqd0I2xAKxxRTa/ncXZmNvNffInc9I34Sg9NnmMIdcQSG5bD8CRDj+GDocto6DAUnCEBjVtEpDaUTNezOXPmcPfddx+07KKLLuLvf/97Vc20tZbg4GD+85//BChKEWlMfMUeStZkUbwqg7JNuWDB1S6cqImJhA1shTPOHegQj0tBXj7zXniJrB9+xFdaeFDZhqMyeY4Ly+GUfkF0Sx4CiWOg3WBNjiIiTZLKPESOk+4xqQ/W46UkPZvitAxKf8wGr8UR7yZsUCvCBrfG1brplDQUFRczf9ar7Fu5Gl9RwUEPDB4o24gNy2V4EvQ4ZQwkjoX2g9WmTkSaFJV5iIgEmPVayjblUpy2n5I1WdgyL0GRwUSMbOdPoDtGYJpAKzePx8PHb7zDlq+XYfPzDmpVd+CBwRh3PsP6VNBn5ODK5HmIRp5FpFlSMi0i0oCqekGnZVD8fYZ/MhW3g9ABCYQNbk1It2hMUONPoL/6+DO+n/8p3uwcvL78qklSjAknJCiKKHc+g7uXMXDUAOg6xv/AoGqeRaQFUDItItIAPBnFFKdlUJK2n4qsUnAaQvvEETa4Ne4+cRhnUKBDPKq1363iq9feomzPfnzeAnzW/zCkMW5cJpKIEENSpwKGjeqJq+eZ/j7P6rYhIi1Qs0ymrbVN4k+l0vQ05WcMpOEd6MRR/N1+PDsL/a3suscQeWYnQvsnNOpe0Pv37uOD518if8MWfOWF+GxO5RoXrqBoIl2Grq1zGDO2M8G9x/l7PavPs4hI80um3W43WVlZxMfHK6GWemWtJSsrC7e7aXVWkIblK/dSujaL4u/2U7ohB3zgah9O9DldCRvUCkd04yx1KCst4/2XZ7N72Up8xQXVOm4E4TQxhLniaB+TxRmnJRDZ/2zoOlYzDIqIHEGzS6Y7duzIzp07ycjICHQo0gy53W46duwY6DAkwKzPUrYxl+Lv9lOyJhNb7sMRE0Lk2E6EDWmFq014oEM8om+++JKV735ARXY2Xl8u1pYC/ocGQx0xxIfnMGpoKB2HjYFuZ0J8d9CghIjIUTW7ZNrlctG1a9dAhyEizVD5niKKv9tH8XcZ+ArKMW4HYYNbEza4NcGJUY3uQcI9u3bzwbP/pWjLdnyeQnw2D/ip7jnKDUN6ljFwzBDodha0GwSOZvffgohIg9JPTRGRo/Dml1Octp/ilfv9U3oHGdy9Ywkb2prQPvEYV+N5kNDj8fD+y7PZ+U1KZb/nAzMNOnCaGCKCo0lsncPpp3cluM/ZkDgKghvnKLqISFOhZFpE5BAH6qCLVu6nbEOOf0bCTpHEnNed0IGtcIQ3nslG1nyXxpevvIln7/6DWtYdKN1oHZHDmJExtBk6AbqdAVHtAxuwiEgzo2RaRAR/HXT51jyKVu6n5IdMbJnXXwd9RifChrbG1apxtH0rKi7mvWdfJHNVOr7SArw2B7BACMFBUUS6YUD3EoadMQy6j4O2AyGo8Yyei4g0N0qmRaRFq8gqoWjlfopX7sObU4YJrpxQZWhrQro2jglVli9eyoo35lKRlUWFN5cDsw06TCzhzhg6xudw5hkdCe83ERJHQ0hkQOMVEWlJlEyLSIvjK62g5IdMilL3Ub41398PukcM0eMTcfeLJyjYEdD4ioqK+N8zL5D1/Tp8ZQVVPZ+NCSUkKIro0HxGDAii16ix0GMcxHULaLwiIi2ZkmkRaRGsz1K2OZfi1P2UrM7Eenw4W4USNTGRsCGtcQa4H/SqFSksfvUtKjKyqPDm4B99NjhMLBGuWBJbZXPmuB4E953on23QGRzQeEVExE/JtIg0axVZJRSl7qN45X68uWX+dnZDWxM2rA3BnSIDNrmTx+Ph3ednsW95WmXtczbgb1sXEhRFTGg+pw1y0m30WdDjZxCt/uYiIo2RkmkRaXZ8Zd7KMo69lG+pLOPoGUv0pERC+yYErJ3dlg2b+OT5lyjbtQevN6+q84Z/9DmOTgnZ/OxnPQnuNwk6jdTos4hIE6BkWkSaBWst5dvyKUrZR8n3mdhyL854N1ETuhA2tE3Ayjg+f28+a+Z/hq8wD68vG/ABwQQHRRPlhlP6+Ug6fSz0HA+xXQISo4iI1J2SaRFp0rz55RSt3Edxyj4qMkswwUGEDmxFeHIbgrtEnfQyjqKiIt596nmyV/+IrzwPn80HKvs+O2NpF53NuDM7EjVoMnQdA67QkxqfiIjUr6Mm08YYNzAZGAO0B0qA1cAH1to1DR+eiMjhrNdH6bpsilbso3R9NvggODGK2DM6EjqgFUEhJ7cbx6Yf1/Ppcy9TvmdvtYcHg3CaWCKDo0hKLGTUhOHQayK06gMBqtMWEZH6V2MybYy5D38ivQhYBuwH3EAv4MHKRPu31trvT0KcIiJ4MoopWrGP4pX78BV6CIoMJnJsR8KGtTnpk6osnP8xP8z7GF9BbrXyDf/EKbGh+YweHkriqHP9Dw+GxZ3U2ERE5OQ52sj0cmvtjBrWPWKMaQ10boCYRESq+Mq9lHyfSdGKvZRvy4cgcPeJJ3x4G9y94jCOkzPKe6D7xt5l3+Ery8dncwF/+UaYM5b2sdn8bFw3wgdP8beuczSeKcdFRKTh1JhMW2s/AH+ph7W2tPo6Y0yCtXY//tFqEZF6V76rkKLleyhOy8CWeXEmhBI9KZGwoW1wRJ6cLhfZWdn874lnKdy4BW9FHtYWAQaniSUiOJo+XQoZPWE4ps+5kNBT5RsiIi1QbR5AXGGMudFa+y2AMeYi4O/4yz1EROqNr7SC4rT9FC3fi2d3ETiDCBuQQPgpbQlOPDkPE27esJFPnvkv5Xv2UeHNBsoBJ66gGKLclpEDDX3OmujvvhGe0ODxiIhI41abZPpK4EVjzCL8DyHGA2c1ZFAi0nJYaynfXkDR8r2UfJ+B9fhwtQsn5rzuhA1uTVBowzcdSl36LUtffQtfTjYVvmzAC7gJCYomPiKPs8a0ps3IC6Dr6eByN3g8IiLSdBzzfylr7Q/GmL8BrwAFwFhr7c4Gj0xEmjVfsYfi7/ZTuHwvFfuKMcEOwoa0JvyUtrg6RDT4KPTn8+azZt6n/v7PNgsAYyIIdcTQPi6Hsyf0InzwBdBhGAQFZpIXERFp/I6ZTBtjXgC6AwPxl3bMN8Y8Ya19qqGDE5HmpWpileV7Kf4+Eyp8uDpGEHNhD8IGtSIopOFGoSsqKvjwtTfYunAp3tK8ag8QRhPujKV7+zzOOGcYrv7nQStVsYmISO3U5n+uH4AbrLUW2GKMGQE80rBhiUhz4iup8I9CL9vjH4UOcRA+rDXhp7QjuENEgx3X4/Hw7nP/9XfgqDaBisPEERUcTd9uxYw69wzocy5EtW+wOEREpPmqTZnHo4e8zgN+0VABiUjzYK3Fs7OQwmV7KFlVWQvdMYLYC3sSOqjhJlYpKynl7aeeIzNtTWUHjkL8HTjiiHZHkNzPx8Dxk6DXBPV/FhGRE3a0SVveB54DPrbWeg5Z1w24DthqrX2xQSMUkSbFV+71d+T4dg+e3UWY4KCqWujgjpENcszikhLeeuzf5K5eh7ciF2uLgSBcQXFEuS2nDgmm98/Ohx5nQ0jDjYSLiEjLc7SR6RuB/wP+ZYzJATLwz4CYCGwCnrTWvtfgEYpIk+DZW0Thsj0Ur9yPLfPiahtOzPmVHTnc9V8LXVRUxJv/eor8dRvxVuRgbQngwBUUS0yoZczISLqecQl0O1MdOEREpMEcbdKWvcDvjDE7ga/xJ9IlwHrrH/YJGGPMFGBKjx49AhmGSItnK3yUrM6k8Ns9lG/NB6chbEArwke2I7hzZL135KhKoNM34vUeSKD9PaDjwixnjImn45hLIHEsOE/OxC4iItKy1Wa4qDXwFrASeBFY1aAR1YK19n3g/eTk5BsDHYtIS1SRW0rRsr0UrdiLr9CDI95N9DldCRvWBkd4/U6jXVRczNuP/ZvcNT9WG4H2J9Dx4ZazTm9Fu9GXQ5dR4Gj4ntQiIiLV1eYBxHuNMX8CxgM/B540xrwJvGCt3dTQAYpI42B9lrJNuRR+s4fSdH9fZnefOCJObU9IjxhMUP2NQpeXlvPmk/8mK21NtRron0agx53Z+qcEOqhhHmQUERGpjVoN41hrrTFmL7AXqABigbeNMZ9Za3/XkAGKSGD5SisoStlH0bd7qMgsISjcReTpnQg/pS3OuPqrRfZ4PLz97xfYvyINrycHa4s4UAMdG2Y56/QEOoy9Qgm0iIg0KrWZtGU6cA2QCfwHuMta6zHGBAEbACXTIs2QZ28Rhd/spvi7/dhyH8GdI4m9rDdhAxIwzvqbEXDurFfZtnAJ3rJcrC3gQBeOmFDLGaNj6HzmldBltEo4RESkUarN/05xwIXW2m3VF1prfcaYyQ0TlogEgvVaStZmUfTNbso25/kfKBzUmohT29VrW7tP353Huvc/wVuSi8/mUdUHOjScsadG0m3c5dD1dHDUb/21iIhIfatNzfSMo6xLr99wRCQQvEUeilbspeibPXjzynDEhBA9KZGw5Lb19kDh8q+X8O3Lb+IryMZrcwBwmjii3JGcMsTFgEmXQ/ezwBlSL8cTERE5GfR3U5EWrHx3IYVLd1OclgEVPkK6RxMztRvupPh6eaBw07p1fPzkC3iyMvD6sgCLw8QQ6YphSF8Pw6eeB70mQnD4iZ+MiIhIACiZFmlhDpRyFC7dRfmWfIwriPBhrYk4tT2utiee1GZnZvPWQ49Run0nFb5MwIsxEYQ7Y+ndpZAzL5kIfSZDaMwJH0tERCTQlEyLtBC+kgqKVuylcOluvLllOGJDiD6nK+HJbQgKO7FSjrLSMt547Glyvk+noiILKMMYN25HLF1a5zHholG4Bl8CkW3r52REREQaCSXTIs2cJ6OYwiW7KV65z9+Vo2s0MZO74e574qUcH8x+k40fL8RbnoO1hYCT4KBYWkfnMWnqQKJGXAXx3evnRERERBohJdMizZC1lrKNuRQu3kXpjzngMIQNbk3EqPYEt484oX2vWr6CL//zGt78LHw2hwOdOGLCffzszHZ0OOsaaD8U6nkqcRERkcZIybRIM2I9PorT9lO4ZBeevcUERbiI+llnwke2wxERXOf9Zu7fz1v/fJzyXbsr66B9BJloooKjGTHEwcAp06DbGeoFLSIiLY7+5xNpBryF5RR9u4fCb/fgK/TgahtO7MW9CBvcqs4TrHg8Ht564lkyUr+vVgcdRpgzlh6dijj78smQNBlC6q//tIiISFOjZFqkCfPsL6bw610UfbcPKizuPnFEjO5ASPdoTB3LLBZ/uoCVr8/FW5yNz+YDDoKD4mgbk8+5Fw4n7JQrIbpD/Z6IiIhIE6VkWqSJsdZStimPwq93+uuhnUGED2tDxOgOuFqF1Wmfe3fv4d1/Pkb5vr14fZkAOEwccaERnHV6K7pMuB7aDVYdtIiIyCGUTIs0Edbro+T7TAq+3olndxFB4SdWD+3xeHjrqWfJWPE9FRWZQDnGhBPujGNA73JGXXoZ9DgbnHWvtRYREWnulEyLNHK+sgqKlu+jcMkuvLllOFuFEnthT8KGtMa4jr8eOnXptyx5YTbeoszKMg5/O7t2cflMuWw0IcOuhPD4+j8RERGRZkjJtEgj5c0vo3DJbgqX7cGWegnuGkXM1O64+8Qdd3/ogvwCXn/wEYq3bKvqxuEwccSGRnLmGQl0nXgTtB3QMCciIiLSjCmZFmlkPPuLKfhqJ8Xf7QefJbR/AhFjOhDSOeq49/XJO3NZ994neMuzsLYYY9yEOmPp3bWEcVddAj0nqIxDRETkBCiZFmkkyrbnU7BoJ6XpWeAIInx4WyLHdMAZH3pc+9m3dy/vPPhotYcJDa6geOIjLedeOISYUT+HqHYNcxIiIiItjJJpkQCy1lL6Yw4FX+6gfEs+JtRJ5JmdiDit/XE/VDh31qtsW/A1FZ5M/D2hI4hwxXLKEAdDLr4ZOo9UNw4REZF6pmRaJACs11LyQwYFi3bi2VuEIzqY6HO7EX5KW4JCHLXez+7tO/jfPx/Hk7kPr80GgggOiqdtXAFTrxpHyNBLwR3dcCciIiLSwimZFjmJrMdH0cp9FHy5E292Kc7WocRe0ouwQcc3U+H/XnyZ7V8srhyFLifIRBEVEs2Y0yLpc94vod2ghjsJERERqaJkWuQk8JVVUPTtXgoW78RX4MHVMYKYc5NwJ8XXujPHnl27effBR/Fk7D1oFLpDq0KmTJuIa9DFEFy3SVtERESkbpRMizQgX7GHgiW7KVy6G1tSQUiPGCIv60hI95haT/f9wZy32PjBgqpaaP8odAxjRkfT57zboE3fhj0JERERqZGSaZEG4C0op2DxLoq+2YMt9+LuG0/UmZ0I7hRZq/fn5WQz+/7/R9nuXZUdOfyj0O1bFTD1mkn+UWjX8XX5EBERkfqnZFqkHlXkllH41U4Kl+8Fr4/Qga2IOrMTrrbhtXr/0s++YMWrb+Ety8TaEowJJ9IVy4gRYQy65FfQtn8Dn4GIiIgcDyXTIvWgIruUgkU7KErdBxbChrYm8oxOuBKOPXrs8Xh49cFHyF/7IxW+DMDiDEqgTQxMnfYzwpKvgJCIBj8HEREROX5KpkVOQEVmCfkLd1D83T4wxj/Ryhkdcca4j/neTevW8eEj/6YiPwOfzQdCCHPGMag/nDbtdugwTH2hRUREGjkl0yJ14MkopuCLHRSn7QdHEBGntidybEcc0SHHfO8Hs99g04ef4/FkAB4cJpb48EjOuaA/rc+6BcLjG/4EREREpF4omRY5Dp79xRR8sZ3iVRkYZxARozv4k+jIo89WWFRUxGv3P0Tx1q0HPVCY2LGYKTdOgx4/g6Da95kWERGRxkHJtEgtHJZEj+1I5JgOx5zy+8fVa/j0X8/gKcrA2kKMCSPCFcvoUdH0u/QOiO9+ks5AREREGoKSaZGjqGsS/cHsN9n04QI8nv1ABQ4TR3xUOFOvPp3oU6+F4Np19xAREZHGTcm0yBF4Mkso+Hw7xWn7a51EezweXvrLPyncuBGvLwMIIiQonu6dS5h0y42QOFoPFIqIiDQzSqZFqqnILiX/8+0Uf7cP4wgiYkxlTfRRkug9O3fxzt/+H56cvfhsPsaEEuGK47RR0Qy4/LcQm3jyTkBEREROKiXTIvgnWylYuJ2iFfsgCH93jjM6HfXBwmWLvmbZi7PxlO3HP813DPHhkZx/1anEjL5evaFFRERaACXT0qJ5C8opWLiDwmV7AAg/pS1RZ3Y6aou7t5/7L7sWLaHCux/w4QpKoGNrF1NuvgpXnwnqyiEiItKCKJmWFslX7KHgq10ULtmF9foIG9qGqHGdccYeebIVj8fDS3/9J4UbNlS2tnPidsQxsH8QY35xJ7Tpe3JPQERERBoFJdPSovjKvBQu2UXBVzuxpV5CB7Ui6uwuNU77nZWRyZz7HsSTuQefzfPXQwfHctb4DvS84LcQ0eokn4GIiIg0JkqmpUWwFT4Kv91DwcId+Io8uJPiiBqfSHC7I7eo27B2LR89/DQVxfuxtpggE01ceCTnXz2G2DHXg+vY04WLiIhI86dkWpo167MUf7ef/M+24c0tI6RbNFETEwnpHHXE7ZcuWEjKS6/jKd8PeHCaeNq2Mlz4yytxJU1SPbSIiIgcRMm0NEvWWkrTs8n7ZCsV+4pxdYgg9qKehPSIwRyh1/O8l15jyycLKx8qtAQHJdA90cM5v7oD2g8+6fGLiIhI06BkWpqdsq155H20lfJt+TgTQom7sg+h/RMwQYcn0a898gSZK76jwrefAw8VDhkawmm/+L36Q4uIiMgxKZmWZsOzv5i8j7ZQmp5NUFQwMRf0IDy5DcZxcGmGf6bCf1C4cQNeXxYQQrgrnrHj2tH30rsgPD4wJyAiIiJNjpJpafK8+WXkL9hO0Yq9mGAHURMSiRjVnqBgx0HbFRUV8fKf/kbp7u34bC7GhBEZEsOkiwfSaeLtEBwWoDMQERGRpkrJtDRZvtIKCr7cSeHiXVifJeK09kSe1RlHuOug7fJysnnl3r/hydqNzxYQZCKJDYvm/OvOJG70teBw1XAEERERkaNTMi1NjvVailbsIf+z7fiKPIQOakX0+C444w/uFb1vzx7emvkPyvP2YG0RQSaa1lFRXHjrpYQPmqrOHCIiInLClExLk2GtpXRdNnkfbqEio4TgrlHEnNOP4E6RB223c+s2/vfXh/EU7cXaEhwmljbxDi6+4xe4epwJR+jmISIiIlIXSqalSSjfVUjeB5sp25yHMyGU+Gl9cfeNO6jN3eYNG5n/wL/wFO8DSnGaeDq0dXHeb27F1fXUwAUvIiIizZaSaWnUvHll5H2yleLv9hMU6iRmanfCR7Q9qEPHpnXr+PDBxykv2QeU4TQJdO7k5oLf3gVt+wcueBEREWn2lExLo+Qr91L49S4KFu3wP1w4piNRZ3YiKPSnW3bd6jV8+vCTeEr2AeW4ghJI7BLO1N/+AVr1ClzwIiIi0mIomZZGxVpLyaoM8j7aijevjND+8URP6nrQw4Ub1qbz0T8fPyiJ7t4Nzv2/P0F898AFLyIiIi1Oo0mmjTHhwNNAObDIWvtagEOSk6x8RwG572+ifHsBrvbhxF3Wi5BuMVXrN6/fwAcP/KuqnMMVlED37oZzfztDsxWKiIhIQDRoMm2MeRGYDOy31vavtnwi8BjgAP5jrX0QuBB421r7vjHmDUDJdAvhzS8n7+MtFK/cT1Cki9iLexI2tE3V9N/bNm7mvb89XPlgYWUS3S2Sc+/6C8R0DmzwIiIi0qI19Mj0LOBJ4OUDC4wxDuAp4GxgJ7DCGDMP6Aj8ULmZt4HjkkbAVvgoXLKL/M93YL0+Is/oSOSZnQgK8d+We3bu4u0Zf6e8cC9QijMogW6JkUz53UyNRIuIiEij0KDJtLX2K2NM4iGLTwE2Wms3AxhjXgfOw59YdwTSAM2m0cyVrMsmb/5mKjJLcCfFEXNuN5wJ/rro7MxsZv/hPsrzd2NtCU4TT2KXMM67ewbEdQ1w5CIiIiI/CUTNdAdgR7XXO4ERwOPAk8aYc4H3a3qzMeYm4CaAzp31J/6mxpNZQt77myj9MQdnq1ASft4Pd+84AAry8nn5D/dRlrULawtxmDg6dnBz8R//rAcLRUREpFFqNA8gWmuLgJ/XYrvngOcAkpOTbUPHJfXDV+6lYOEOCr7aiXEGEX1OVyJOa49xBlFWWsaLf5hJ6e5t+Gw+DhNDm9YxXHzPnbg6DAx06CIiIiI1CkQyvQvoVO11x8pl0gxZayldk0Xu/M14c8sIG9Ka6EldcUQFU1FRwUt//jsF69fhtTkEmShax8Zw+R9+g6tzcqBDFxERETmmQCTTK4Cexpiu+JPoy4ErAxCHNDBPZgm58zZRtj4HV9sw4m4aSEi3aABmP/Ik+1eswOvLwphwYiNiuOx3NxLe+/QARy0iIiJSew3dGm8OcAaQYIzZCcyw1r5gjLkd+AR/a7wXrbVrGjIOObmsx0v+wh0UfFlZ0jG5GxGntsc4DO+/MofNH35GhW8/xriJcsdy4e0XEJ98ARgT6NBFREREjktDd/O4ooblHwIfNuSxJTBKf8wm571NeLNL/SUd53TFERnMskVfs+z5l/BU7AOchLvimXT1aXQZfyMEqXmLiIiINE2N5gFEadq8eWXkzt9MyQ+Z/i4dNwzA3SOGzes3MP+2/4enbC/gw+1I4PTJPel/2V3gcAU6bBEREZETomRaToj1Wgq/2U3+p9uwPkvU+C5Eju1IfkEu/7np15W9oksJDmrFwFNiOP22+yE4PNBhi4iIiNQLJdNSZ+U7C8h5dwOe3UW4e8cSM7U7NsrJ83/4M0XbN+GzBf4JV7pHc94f/wHhCYEOWURERKReNclk2hgzBZjSo0ePQIfSIvnKvOR/to3CJbsIiggm7qo+hPZP4M2nnmXP4m/w2iyCTBRt4mO5bOa9uFr3DHTIIiIiIg3CWNt05z1JTk62KSkpgQ6jRSlZl03u3I14c8sIH9GW6IldWfzlF6S9/AYe7z6McRMZGsplv7ueqKQzAx2uiIiISL0wxqRaaw+bCKNJjkzLyectKCf3/U2UfJ+Js3UYrW4ZSAa5vHLzbZUPF0KoM4EJV4+m+8RfqM2diIiItAhKpuWorLUUp+4nd/5mbIWXqLO7EDwigRfu+ROlmTuwtpjgoFYMGd2O0TfPBGdwoEMWEREROWmUTEuNKrJLyfnfBso25BKcGEXsRT158+XnyHg+Fa/NwWFi6dw1igtnPAShsYEOV0REROSkUzIth7G+ynZ3n2wFDDHndSctbw0rpz9eWRcdSmxkHFfN/B0hHfsHOlwRERGRgFEyLQfx7C8m550NlG/LJ6RXLIyN5b9/+yvlxXsAS5gzgXNumkCX0484uaWIiIhIi6JkWgD/5CsFX+8kf8E2goIdxF7aizlvPknBpxvx2QJcQQkMGdWRMb+cCQ7dNiIiIiKgZFqoHI1+az3lOwpw94vnO5PO2oeeo8KXQZCJon27VlzxwP+DsLhAhyoiIiLSqCiZbsGsz1L49S7yPttKULADxiUw5+VHKS/fAziICInn0t/fRGzSqECHKiIiItIoNclkWjMgnjhPRuVo9PYCQpJi+WD122S/kI61RQQHteK0yf0YduVv1S9aRERE5Cg0A2ILc6BTR95HWzGuIHa0y2HZwrcqSzqi6dg5jEvufxSCwwMdqoiIiEijoRkQhYrcMnLeXk/ZxlxM13A++PZlCtZuAYKIcCdw+Z9vJ7r7YfeIiIiIiNRAyXQLYK2l+Lv95M7bBD7LWuePrFn4OT5bQHBQK06Z1JcR0+5USYeIiIjIcVIy3cx5izzk/m8DJauz8MQZPv/hLfLKNxFkImnbrg1X/eNRCIkMdJgiIiIiTZKS6Was5Mdsct5aj6+kgrUlP7A69Qss5YS6Erjodz+nzcDTAx2iiIiISJOmZLoZsh4feR9toXDpbkpCylm8fT7Z5RtwmDiShvVgwm8fgKCgQIcpIiIi0uQpmW5mPHuLyJqzjop9xWwsWsd3Wz/EZ4OIjkhg2oN/IaRV50CHKCIiItJsKJluJqy1FC7ZTd7HWyj3lrNs36fsLl5NcFArxlx6KoMvvCnQIYqIiIg0O0qmmwFvQTnZb62nbH0Oe0p2sGz/XMqtoXWbtkx7+HEIDgt0iCIiIiLNkpLpJq50Qw7Zb/yIp7CUtKyv2Zi/jBBHK87/1SV0O+2cQIcnIiIi0qw1yWRa04mD9frI/2wb+Yt2UlCew9L971JQUU7n7l245K+PgaNJfmlFREREmhRNJ94EVeSUkj1nHeXbC9hUsJrvsj4l2BnDxff+ktZ9Tgl0eCIiIiLNjqYTbyZKVmeS+UY6FWXlrMj8iF3Fe+jatxvn/+lhzWAoIiIicpIpmW4ibIWPvA/9vaNzyvaxdP97+IyTafffTUKPgYEOT0RERKRFUjLdBFTklLL7PysJyvLyY94Kfsj5jsR+bTn/3n9qNFpEREQkgJRMN3IlP2az979pWJ+PFRkfkeXJ5pq/301c1/6BDk1ERESkxVMy3UhZn2X3u6uxKbkUleXwTcZnxHYP49aZL2o0WkRERKSRUDLdCHkLykl/6DNiyqPYXLCKNbmrOO+eG2g/YESgQxMRERGRapRMNzI5a/eQ+d/viTChLM9eQElMJr984lkICgp0aCIiIiJyCCXTjciKp+fTZlsEXl85SzM+ZcQN4+gzVrMYioiIiDRWSqYbgfLiMlLufZ3Owd3YU7KV9PJvuPa/j4HDFejQREREROQolEwH2I9LvqPs7S10DulGeu53uEc7uebqpwMdloiIiIjUQpNMpo0xU4ApPXr0CHQoJ2T+zMfpU9SbSFcUy7M/Z8KDvyQ8rk2gwxIRERGRWmqST7VZa9+31t4UHR0d6FDqpDC/gM9ufYKBJYMo95axNvwbLnxuphJpERERkSamSY5MN2XL536Ia1EhSVGD2V64kXbX9ubc4X8KdFgiIiIiUgdKpk+i9373MP3sAJzBCXxfuJhJj92FcehLICIiItJUKZM7CbL2ZvDDX99hWOQIcsszyei6h3N+9ftAhyUiIiIiJ0jJdAP76tW3iUk19Ijsx+aCdHpNH0v/XhcEOiwRERERqQdKphvQvLseoh+DMU5IK/2KyU/+AYwJdFgiIiIiUk+UTDeA/OwcUv40h6GRI8kpyyCnbw6Tb/xjoMMSERERkXqmZLqerZj3Ec4vCugVOYAtBT/S8zdjGdCzV6DDEhEREZEGoGS6Hr3/x/9Hn/J+uIITSCtawuQn71ZZh4iIiEgz1iQnbWlsysvK+OTWRxlUkYzHV87GdulMfuIeJdIiIiIizZxGpk/QzvR0dj69gn5Rw9hZtIXWP+/Hz4ZODXRYIiIiInISKJk+AYv+8zKt1kbRPqwLawtSOOORXxIc4g50WCIiIiJykiiZrqP3f/cP+tmhWIdllfmWKU/dHeiQREREROQkUzJ9nEqLivj6rhcYHHkquZ4scgcVM+U6JdIiIiIiLZGS6eO04A9PMjjqNLYXbqLLr05jQG+1vRMRERFpqZpkNw9jzBRjzHN5eXkn/dg/e+B2vitfSvL/u4IOSqRFREREWjRjrQ10DHWWnJxsU1JSAh2GiIiIiDRzxphUa23yocub5Mi0iIiIiEhjoGRaRERERKSOlEyLiIiIiNSRkmkRERERkTpSMi0iIiIiUkdKpkVERERE6kjJtIiIiIhIHSmZFhERERGpIyXTIiIiIiJ11KRnQDTGZADbAnDoBCAzAMdtynTNjp+u2fHTNTt+umbHT9fs+OmaHR9dr+N3Mq5ZF2ttq0MXNulkOlCMMSlHmk5SaqZrdvx0zY6frtnx0zU7frpmx0/X7Pjoeh2/QF4zlXmIiIiIiNSRkmkRERERkTpSMl03zwU6gCZI1+z46ZodP12z46drdvx0zY6frtnx0fU6fgG7ZqqZFhERERGpI41Mi4iIiIjUkZJpEREREZE6UjJ9CGPMRGPMj8aYjcaYe46wPsQY80bl+mXGmMRq635fufxHY8yEkxp4ANXimv2fMWatMeZ7Y8znxpgu1dZ5jTFplR/zTm7kgVOLa3adMSaj2rW5odq6a40xGyo/rj25kQdOLa7Zv6pdr/XGmNxq61rcfWaMedEYs98Ys7qG9cYY83jl9fzeGDO02rqWeo8d65pdVXmtfjDGLDXGDKq2bmvl8jRjTMrJizqwanHNzjDG5FX7/vtztXVH/Z5ujmpxve6qdq1WV/7siqtc11LvsU7GmIWVecQaY8z0I2wT2J9n1lp9VH4ADmAT0A0IBlYBfQ/Z5lbgmcrPLwfeqPy8b+X2IUDXyv04An1OjeSanQmEVX7+ywPXrPJ1YaDPoZFes+uAJ4/w3jhgc+W/sZWfxwb6nBrDNTtk+18BL1Z73RLvs7HAUGB1DevPAT4CDDASWFa5vEXeY7W8ZqcduBbApAPXrPL1ViAh0OfQCK/ZGcD8Iyw/ru/p5vJxrOt1yLZTgC+qvW6p91g7YGjl55HA+iP8nxnQn2camT7YKcBGa+1ma2058Dpw3iHbnAe8VPn528A4Y4ypXP66tbbMWrsF2Fi5v+bumNfMWrvQWltc+fJboONJjrGxqc19VpMJwGfW2mxrbQ7wGTCxgeJsTI73ml0BzDkpkTVS1tqvgOyjbHIe8LL1+xaIMca0o+XeY8e8ZtbapZXXBPSzDKjVfVaTE/k52GQd5/Vq8T/HAKy1e6y1Kys/LwDSgQ6HbBbQn2dKpg/WAdhR7fVODv+CVW1jra0A8oD4Wr63OTre8/4F/t8eD3AbY1KMMd8aY85vgPgao9pes4sq/1z1tjGm03G+t7mp9XlXlhF1Bb6otrgl3mfHUtM1ban32PE69GeZBT41xqQaY24KUEyN1anGmFXGmI+MMf0ql+k+OwpjTBj+pO+daotb/D1m/KW1Q4Blh6wK6M8zZ33vUKQmxpirgWTg9GqLu1hrdxljugFfGGN+sNZuCkyEjcr7wBxrbZkx5mb8fw05K8AxNRWXA29ba73Vluk+k3pjjDkTfzI9utri0ZX3WGvgM2PMuspRyJZuJf7vv0JjzDnAXKBnYENqEqYAS6y11UexW/Q9ZoyJwP/LxW+stfmBjqc6jUwfbBfQqdrrjpXLjriNMcYJRANZtXxvc1Sr8zbG/Az4IzDVWlt2YLm1dlflv5uBRfh/42zujnnNrLVZ1a7Tf4BhtX1vM3U85305h/xptIXeZ8dS0zVtqfdYrRhjBuL/njzPWpt1YHm1e2w/8D9aRpnfMVlr8621hZWffwi4jDEJ6D47lqP9HGtx95gxxoU/kX7NWvvuETYJ6M8zJdMHWwH0NMZ0NcYE47+ZD33yfx5w4GnQi/E/HGArl19u/N0+uuL/zXv5SYo7kI55zYwxQ4Bn8SfS+6stjzXGhFR+ngCMAtaetMgDpzbXrF21l1Px14gBfAKMr7x2scD4ymXNXW2+NzHG9MH/kMk31Za11PvsWOYB11Q+BT8SyLPW7qHl3mPHZIzpDLwLTLPWrq+2PNwYE3ngc/zX7IjdGloaY0zbyueKMMacgj/vyKKW39MtkTEmGv9fcN+rtqzF3mOV988LQLq19pEaNgvozzOVeVRjra0wxtyO/0I78HcDWGOM+QuQYq2dh/8L+ooxZiP+hwgur3zvGmPMm/j/k64Abjvkz8zNUi2v2UNABPBW5c/U7dbaqUAS8Kwxxof/B+yD1tpmn+TU8pr92hgzFf+9lI2/uwfW2mxjzF/x/0cE8JdD/gzYLNXymoH/+/H1yl9wD2iR95kxZg7+TgoJxpidwAzABWCtfQb4EP8T8BuBYuDnleta5D0Gtbpmf8b/jMzTlT/LKqy1yUAb4H+Vy5zAbGvtxyf9BAKgFtfsYuCXxpgKoAS4vPL784jf0wE4hZOqFtcL4ALgU2ttUbW3tth7DP8AyDTgB2NMWuWyPwCdoXH8PNN04iIiIiIidaQyDxERERGROlIyLSIiIiJSR0qmRURERETqSMm0iIiIiEgdKZkWEREREakjJdMiIiIiInWkZFpEREREpI6UTIuIiIiI1JGSaRERERGROlIyLSIiIiJSR0qmRURERETqSMm0iIiIiEgdKZkWEREREakjJdMiIiIiInWkZFpEREREpI6UTIuIiIiI1JGSaRERERGROlIyLSIiIiJSR0qmRURERETqSMm0iIiIiEgdKZkWEREREakjJdMiIiIiInWkZFpEREREpI6cgQ7gRCQkJNjExMRAhyEiIiIizVxqamqmtbbVocubdDKdmJhISkpKoMMQERERkWbOGLPtSMtV5iEiIiIiUkdKpkVERERE6kjJtIiIiIg0CdbaQIdwmCZdM30kHo+HnTt3UlpaGuhQpBlyu9107NgRl8sV6FBERERaDJ/Xy8qP5rH9hzQuuHsGJqjxjAc3u2R6586dREZGkpiYiDEm0OFIM2KtJSsri507d9K1a9dAhyMiItIi7Nuyic+ee4J9mzfSbehwyktLCQkLC3RYVZpdMl1aWqpEWhqEMYb4+HgyMjICHYqIiEiz5ykrZelbs0n9YC6hkVFM/s3d9Bo5utHleM0umQYa3UWW5kP3loiISMPbumolC/7zFHn79zHgrPGMvep63BERgQ7riJplMh1oDoeDAQMGVL2+/PLLueeee2rcftasWaSkpPDkk0+ejPBEREREGqWi3By+fPVF0r9eSGy7Dlw64+906jvg2G8MICXTDSA0NJS0tLQG239FRQVOp750IiIi0jxYn4/vP/+Er+fMwlNaxsgLL2PEBZfhDA4OdGjH1HgehWwBEhMTyczMBCAlJYUzzjjjsG0yMjK46KKLGD58OMOHD2fJkiUAzJw5k2nTpjFq1CimTZt2MsMWERERaTD7t25mzp/vYsF/nqJ1l25c89ATjLpsWpNIpEEj0w2ipKSEwYMHV73+/e9/z2WXXVar906fPp077riD0aNHs337diZMmEB6ejoAa9euZfHixYSGhjZE2CIiIiInTXlpCUvfms3KD9/DHRHJpNv+j6QxZza555OadTJ93/trWLs7v1732bd9FDOm9DvqNidS5rFgwQLWrl1b9To/P5/CwkIApk6dqkRaREREmjRrLeu/XcKiV/5DYVYmA8dNZPSV1xIaERno0OqkWSfTjY3T6cTn8wHUOKmMz+fj22+/xe12H7YuPDy8QeMTERERaUhZu3bwxX+fZfsPabRK7Mbk6XfToXdSoMM6Ic06mT7WCPLJlpiYSGpqKpMmTeKdd9454jbjx4/niSee4K677gIgLS3toJIRERERkaamvLSEb999g9T5c3GFhHDWz29m0NnnEORwBDq0E6YHEBvAgZrpAx8H2uLNmDGD6dOnk5ycjKOGm+fxxx8nJSWFgQMH0rdvX5555pmTGbqIiIhIvbHW8uM3i/nv//2SFe+9TdLoM7j+0WcZMnFKs0ikAYy1NtAx1FlycrJNSUk5aFl6ejpJSU37zwXSuOkeExERObaMbVtYOOs5dqz9gVaJ3Rh3/S+bdEmHMSbVWpt86PJmXeYhIiIiIidXSUE+S958je8/+4iQ8HDG/eJWBo6b0GxGog+lZFpERERETpjP62XVgo9Y+uZrlBUXMWj8OZx26VVNtktHbSmZFhEREZETsu2HNBa9/B8yt2+lc/+BnHntTSR0Tgx0WCeFkmkRERERqZPs3Tv58pUX2LxyBVGt2jDl/35Pz1NOa3ITr5yIE06mjTGda7lprrW2fmdQEREREZGTrqSwgG/ens2qTz/EGRzMmCuvY+ikqU1mCvD6VB8j0y8BFjjaryAWmAW8XA/HExEREZEA8FZ4WPXph3zz9hzKiosZcNZ4Trv0KsJjYgMdWsCccDJtrT2zPgIRERERkcbJWsuGZUv4es5L5O7dQ+cBgznjmhto1ULqoo+m3idtMcaEG2OaZ++TAEhMTCQzMxOA0047rWr5XXfdRb9+/bjrrrt45plnePnl4xv0j4iIOOJyh8Nx0IQzDz744FH3M2vWLG6//fbjOvah9uzZw+TJk6teL1++nLFjx9K7d2+GDBnCDTfcQHFxca1imDlzJg8//PBRjzd37lzWrl171G2efPJJXnzxxeM4CxERkeZp14/pzPnzXbz/rwdxOF1ccPcMLv7jX5VIV6qPmukg4HLgKmA4UA6EGGMygA+AZ621G2u5LzfwFRBSGdvb1toZJxpjc7F06dKqz5977jmys7NrnEmxrkJDQ0lLS6vXfVZXUVGB03nwbffII49w4403ArBv3z4uueQSXn/9dU499VQA3n77bQoKCggLC6uXGObOncvkyZPp27dvjdtcf/31jBo1iuuvv75ejikiItLU5OzZxdezX2LD8qWEx8Zx9k2/ov8ZP2u2/aLrqj5GphcC3YHfA22ttR2tta2A0cC3wD+MMVfXcl9lwFnW2kHAYGCiMWZkPcR40mzdupU+ffpw3XXX0atXL6666ioWLFjAqFGj6NmzJ8uXLwcgOzub888/n4EDBzJy5Ei+//57ALKyshg/fjz9+vXjhhtuoPoMlQdGk6dOnUphYSHDhg3jjTfeOGg0dtOmTUycOJFhw4YxZswY1q1bB8CWLVs49dRTGTBgAPfee+9xn1f1EfKUlBTOOOOMw7bJyMjgoosuYvjw4QwfPpwlS5YA/tHiadOmMWrUKKZNm3bY+9555x0mTpwIwFNPPcW1115blUgDXHzxxbRp06bGa1aT559/nuHDhzNo0CAuuugiiouLWbp0KfPmzeOuu+5i8ODBLFmy5KCReIfDwbZt2wgLCyMxMbHq6yUiItJSFOXmsOCFfzPrt7ey9fvvOO3Sq/jFo88164lXTkR9PID4M2ut59CF1tps4B3gHWOMqzY7sv7MsbDypavyo+7znX90D+z9oc5vP6K2A2DS0UsfNm7cyFtvvcWLL77I8OHDmT17NosXL2bevHk88MADzJ07lxkzZjBkyBDmzp3LF198wTXXXENaWhr33Xcfo0eP5s9//jMffPABL7zwwmH7nzdvHhEREVUjyDNnzqxad9NNN/HMM8/Qs2dPli1bxq233soXX3zB9OnT+eUvf8k111zDU089VWPsJSUlDB48uOr173//ey677LJaXZrp06dzxx13MHr0aLZv386ECRNIT08HYO3atSxevJjQ0NCD3rNlyxZiY2MJCQkBYPXq1Vx77bVH3H9N16wmF154YdWI97333ssLL7zAr371K6ZOncrkyZO5+OKLAar28dRTT/Hll1/SpUsXAJKTk/n666855ZRTanX+IiIiTVlZcREp779Lygdz8VVUMOCsCZx68RUt+uHC2qiPBxA9AMaYx4Df2OpDqYdsUxuV9dapQA/gKWvtskPW3wTcBNC5c2278p1cXbt2ZcCAAQD069ePcePGYYxhwIABbN26FYDFixfzzjvvAHDWWWeRlZVFfn4+X331Fe+++y4A5557LrGxtb+BCwsLWbp0KZdccknVsrKyMgCWLFlSdbxp06Zx9913H3EfJ1LmsWDBgoNqkfPz8yks9P9uNHXq1MMSafDXS7dq1apW+6/pmtVk9erV3HvvveTm5lJYWMiECRNq3HbJkiU8//zzLF68uGpZ69atq0b2RUREmquK8nLSPv2AZXPforQgn96njmHU5dOIbds+0KE1CfU5aUsBMM8Yc7m1tsgYMwH4s7V21PHsxFrrBQYbY2KA/xlj+ltrV1db/xzwHEBycvLRR62PMYLcUA6MsgIEBQVVvQ4KCqKioqLBjuvz+YiJiakxGT6RBupOpxOfzwdAaWlpjcf/9ttvcbvdh60LDw8/4ntCQ0MP2l+/fv1ITU3lvPPOq3OsB1x33XXMnTuXQYMGMWvWLBYtWnTE7fbs2cMvfvGLqhH/A0pLS4/4C4CIiEhz4PN6WfvVFyx9ezYFmRl0GTiEMVdcS5tuPQIdWpNSb908rLX3AnOARcaYJcD/AfecwP5y8ddjT6yXABuZMWPG8NprrwGwaNEiEhISiIqKYuzYscyePRuAjz76iJycnFrvMyoqiq5du/LWW28B/jY2q1atAmDUqFG8/vrrAFXHPR6JiYmkpqYCVI0OH2r8+PE88cQTVa9rM8Ldq1evqtF6gNtvv52XXnqJZct++oPEu+++y759+2q8ZjUpKCigXbt2eDyeg845MjKSgoICADweD5dccgn/+Mc/6NWr10HvX79+Pf379z/mOYiIiDQl1udj3ZIvmfXbW/nkmccIj47hkj/9jYv/+Fcl0nVQb8m0MWYccCNQBCQAv7bWfn2c+2hVOSKNMSYUOBtoln9nnzlzJqmpqQwcOJB77rmHl156CfDXBX/11Vf069ePd99997hLWV577TVeeOEFBg0aRL9+/XjvvfcAeOyxx3jqqacYMGAAu3btqvH9B2qmD3zcc889VXFNnz6d5OTkGjuIPP7446SkpDBw4ED69u3LM888c8x4w8PD6d69Oxs3+hu+tGnThtdff50777yT3r17k5SUxCeffEJkZGSN16wmf/3rXxkxYgSjRo2iT58+Vcsvv/xyHnroIYYMGcLSpUtJSUlhxowZVee8e/duwF/6cfbZZx/zHERERJoCay0bU5bxyt2/5oPHH8LhdHLenfdy5d8eoXP/QYEOr8kyRyhxrtuOjPkCf1nHYmPMAOAV4P+stV8cxz4G4p9R0YE/0X/TWvuXmrZPTk62KSkpBy1LT08nKSmpLqcgAfK///2P1NRU7r///kCHUuW7777jkUce4ZVXXjlsne4xERFpSqy1bPv+O5a8+Sp7N64ntl17Tr3kKvqcOgYTVO9TjjRbxphUa23yocvrrWbaWntWtc9/MMZMwt/N47Sa33XYPr4HhtRXTNI0XHDBBWRlZQU6jINkZmby17/+NdBhiIiI1Jm1lu2rV7H0rdns/nEtkfGtGH/zr+l3+ji1uKtH9TFpi6mhg8eeytKPGrcROeCGG24IdAgHUXmHiIg0ZdtXf8/St15j17o1RMTFM+76X9L/rPE4XbXqVizHoT5GphcaY94B3rPWbj+w0BgTDJxqjLkW/4OEs+rhWCIiIiJSgx1rf2DpW6+xc+1qImLjOOvnNzPgrAk4g4MDHVqzVR/J9ETgemCOMaYrkAu48dc9fwo8aq39rh6OIyIiIiKHsNay/YdVfPvu6+xMX014TCxnXncTA8dNVBJ9EtTHpC2lwNPA05UzHSYAJZWt7URERESkAVhr2ZKWwrfvvM6eDT8SERfPmdfdxIBxE3AFhxx7B1Iv6u0BRGPMZ8Cd1tpV9bXPpioiIqJq5j+AWbNmkZKSwpNPPhnAqERERKQ5sD4fm1KX8+27r7Nv80YiE1rxsxtupd8ZZ6smOgDqcwbEu4FHjTFbgT9Ya/fU475FREREWjSf18uPS79i2dy3yNq5neg2bRl/86/pO/ZMHE4l0YFSn63xVgJnGmMuAj42xrwL/NNaW1Jfx2gOMjIyuOWWW9i+3f+s5qOPPsqoUaOYOXMmERER3HnnnQD079+f+fPnAzBp0iRGjx7N0qVL6dChA++9956muRYREWkhKsrLWfPlAlbMe4e8/fuI79iZc27/Lb1PG6sWd41AfY5MY4wxwI/Av4H7gRuNMb+31h4+80UzdmAWwQOys7OZOnUqANOnT+eOO+5g9OjRbN++nQkTJpCenn7U/W3YsIE5c+bw/PPPc+mll/LOO+9w9dVXN+QpiIiISICVFRfz/YKPSP1gLkW5ObTt0YszrrmR7sNO0WQrjUh91kwvAboCa4BvgevwTwU+3Rgzxlp7U30dq7b+sfwfrMuu39nI+8T14e5T7j7qNqGhoaSlpVW9PlAzDbBgwQLWrl1btS4/P/+g+uoj6dq1a1VyPmzYMLZu3Vqn2EVERKTxK8zJZuVH81j16YeUlxTTecBgzvnVnXTqNxD/uKU0JvU5Mn0TsPYIk7P8yhhz9KHXFsTn8/Htt9/idrsPWu50OvH5fFWvS0tLqz4PCfnpiVyHw0FJiSpnREREmpvs3TtZMe9d0r/+Ap/XR88RpzF8yoW07dEr0KHJUdRnzfSao6w+t76OczyONYIcCOPHj+eJJ57grrvuAiAtLY3BgweTmJhYVSO9cuVKtmzZEsgwRURE5CSw1rL7x3RS5r/LxpRlOJ0u+p85nmGTzye2bftAhye1UK810zWx1m4+GcdpCh5//HFuu+02Bg4cSEVFBWPHjuWZZ57hoosu4uWXX6Zfv36MGDGCXr30W6iIiEhz5fN62bB8KSnz/8fejetxh0cw8sLLGDJhMmHRMYEOT46DObwq4wR3aMwUa+379brTGiQnJ9sDtcgHpKenk5SUdDIOLy2U7jEREamr8pJifvjiM1Z+9B75GfuJaduOYeecT7/Tx+E6pARUGhdjTKq1NvnQ5Q0xMv034KQk0yIiIiJNQd7+fXz3yXxWf/EpZcVFdOjTjzOu9XfmCApSe7umrCGSaT1mKiIiIi2etZad6atZ+eE8NqUsAwO9Roxi2OTzadejd6DDk3rSEMl0/daNiIiIiDQhFeXlrFv6FSs/mkfG1s24IyIZft5FDB5/LpHxCYEOT+rZSXkAUURERKS5y8/cz6rPPuKHzz+hpCCf+I6dOfum20kafQauENVDN1dKpkVERETqyFrLjjXf893H8/2lHED35FMYPGEynfsP0iQrLUBDJNP7GmCfIiIiIo1GWXExa7/+grRPPiB71w7ckVEMn3ohg84+h6hWrQMdnpxE9T6xu7X27PreZ1Pk9XoZMmQIkydPrlr26KOPUlxcXPU6IiKiwY6fmJhIZmZmnd9/3XXXVU1jPmjQID7//PN6jO74zJo1i9tvvx3wzyB57bXXcv3111NTW8fc3Fyefvrpqtdbt25l9uzZJyVWERFp3vZv3cxnzz3Js7dcwxcvPoMrJISJt97BzU/PYsyV1ymRboFU5tFAHnvsMZKSksjPz69a9uijj3L11VcTFhZWr8eqqKjA6TyxL6XX68XhOLg1z0MPPcTFF1/MwoULuemmm9iwYcMJHeNEWWu55ZZb8Hg8/Pe//63xT2cHkulbb70V+CmZvvLKK09muCIi0kx4ystY/81iVn36IXs2/ojTFUzvUWMZfPY5mupb6n9kWmDnzp188MEH3HDDDVXLHn/8cXbv3s2ZZ57JmWeeWbX8j3/8I4MGDWLkyJHs23d4hUx2djbnn38+AwcOZOTIkXz//fcAzJw5k2nTpjFq1CimTZtGVlYW48ePp1+/ftxwww0Hjdq++uqrnHLKKQwePJibb74Zr9cL+EfGf/vb3zJo0CC++eabGs/n1FNPZdeuXcDBo8QAkydPZtGiRVX7O9L5bNq0iZEjRzJgwADuvffeg0bkH3roIYYPH87AgQOZMWPGUa/rr3/9a7Kysnj55ZcJCgpi5syZPPzww1Xr+/fvz9atW7nnnnvYtGkTgwcP5q677uKee+7h66+/ZvDgwfzrX/9i69atjBkzhqFDhzJ06FCWLl161OOKiEjLlLljGwtnPcdzt1zLx0//i9LiIs689kZufuZlJv7yN0qkBWiAZNoYE26MadHdx3/zm9/wz3/+k6Cgny7vr3/9a9q3b8/ChQtZuHAhAEVFRYwcOZJVq1YxduxYnn/++cP2NWPGDIYMGcL333/PAw88wDXXXFO1bu3atSxYsIA5c+Zw3333MXr0aNasWcMFF1zA9u3bAf9sfW+88QZLliwhLS0Nh8PBa6+9VnX8ESNGsGrVKkaPHl3j+Xz88cecf/75xzzvms5n+vTpTJ8+nR9++IGOHTtWbf/pp5+yYcMGli9fTlpaGqmpqXz11VdH3Pfs2bNZuXIlr7/++jFH4R988EG6d+9OWloaDz30EA8++CBjxowhLS2NO+64g9atW/PZZ5+xcuVK3njjDX79618f89xERKRl8JSVsubLz5nz59/x0p23kfbph3QZOIRL/vQAP3/k3ww95zzcDVimKU3PCZd5GGOCgMuBq4DhQBkQYozJBD4AnrXWbqzFfjoBLwNt8Peqfs5a+9iJxLb3gQcoS193Irs4TEhSH9r+4Q81rp8/fz6tW7dm2LBhVSO2NQkODq6qqR42bBifffbZYdssXryYd955B4CzzjqLrKysqtKRqVOn8v/bu/PwOKoz3+PfV2trl6zNtmxZssG2bINjYzZjMNsFO4lhiOcmzsIAvjPOM8FzuZObfWDi5CaQZPJMZu5MyJiHgUAWwVxIgMkQlsRAYgI2Fhi8yLstb/KmfVer+9w/qtW0ZNla3HZr+X309FNVp05VvXW61P2qdKoqJSUFgD/84Q/86le/AuBjH/sYOTk5APz+97+noqKCyy+/HIC2tjYKCrz+XPHx8SxfvvyM8X35y1/mG9/4BocPHz7rmev+9uett97iueeeA+Azn/kMX/rSlwAvmX7llVeYN28eAM3NzezevZvrrrvutHXPnz+fHTt2sHHjRq655pp+Yzkbv9/P6tWrw39c7Nq165zWJyIiI9+JA/vYsu4VKv/4Gh2tLeRMKGLx51Yya/FNpGZmxTo8Gcai0Wf6NeB3wNeBrc65IICZjQNuAL5vZr92zv28n/V0Af/bOfeumWUAFWb2qnNuexRivGDefPNNXnjhBV588UXa29tpbGzkc5/7HD//+em7n5iYGO73Gx8fT1dX16C2lZaW1m8d5xx33XUXDz300GnzfD7faf2kI3X3mf6Xf/kXVq5cSUVFBQkJCQSDwXCd9vb28Phg98c5x9e//nU+//nP9yj/8Y9/HD6r/eKLLwIwc+ZMvv3tb/PJT36Sl19+mdmzZ581lrP50Y9+RGFhIe+//z7BYBCfT/f+FBEZi9qbm6l883W2rnuVEwf2Ep+QwMVXXsOlNy9hUtkc3dZOBiQayfTNzjl/70LnXC3wLPCsmSX2txLnXDVQHRpvMrNKoAgYcjJ9tjPI58tDDz0UTlxff/11fvjDH4YT6YyMDJqamsjLG/jTj6699lp+8Ytf8MADD/D666+Tl5dHZmbmafWuu+46fvnLX3L//ffz29/+lrq6OgBuuukmbr/99nD3htraWpqampgyZcqAY1i9ejWPPfYYL7/8MiUlJTz88MMEg0GOHDnCxo0b+13+qquu4tlnn+VTn/oUTz31VLj81ltv5YEHHuCzn/0s6enpHDlyhMTERO69917uvffe09azcOFCfvKTn/Dxj3+cN954g5KSEn7zm98A8O6777J//37gw3bu1nu6oaGBSZMmERcXxxNPPBHuQy4iIqOfCwY5uO0Dtr72Krs3/omA309+yVRuvOfzzFx0PSnpGbEOUUaYc06muxNpM/tn4H+5Pu5X1leyfTZmVgLMAzaca3zDyapVq1iyZEm47/RArFmzhpUrV3LppZeSmprKE0880We9b37zm3z6059m9uzZLFy4kOLiYgBmzZrFd77zHW655RaCwSCJiYn8+Mc/HlQybWbcf//9/OAHP+B3v/sdpaWlzJo1i7KyMubPn9/v8t13Mfnud7/LkiVLyMry/l12yy23UFlZydVXXw14FzD+/Oc/D3dD6cuyZcs4deoUS5Ys4dVXX+XJJ59k9uzZXHnllUyf7l0IkpubyzXXXMOcOXNYunQpDz74IPHx8cydO5e7776bL3zhCyxfvpwnn3ySJUuWDOgMv4iIjGx1x46y/Q/r2PbG72k6dZLktDQuufEW5txwC4Wl02IdnoxgdqZ79Q56RWbfAeYCK5xzLWZ2K/D3zrlBdXA1s3TgDeC7zrlf9TF/FbAKoLi4+LKqqqoe8ysrKykrKxviXsj50NraSkpKCmbGU089RXl5Oc8//3yswxoyHWMiIiNDR2sru95ez7Y3fseRHdvBjJJL5zFr8U1cdPlVJCYlxzpEGUHMrMI5t6B3edTuM+2cu9/MPgO8bmadQDPwtcGsI9Qd5FngF30l0qHtPAI8ArBgwYLo/CUg51VFRQWrV6/GOUd2djaPPfZYrEMSEZFRKhgIULVlM9v/sI4977xNV2cHORMnsejTdzHr2hvIyB14V0uRgYhaMm1mNwF/BbQAE4CVzrmdg1jegH8HKp1z/xituCT2rr32Wt5///1YhyEiIqOUc44T+/ey/Y+vsePNN2htqMeXls7sxTcye/HNjL9oui4mlPMmmk9A/DvgAefcejO7BHjazL7onFs3wOWvAe4EtpjZ5lDZN5xzL0YxRhERERklGk4cZ8ebb1C5/nVqDh8kLj6BqfMvZ9Z1N1A673ISEvu9/4HIOYtmN48bI8a3mNlSvC4bCwe4/HpAfzaKiIjIGbU21LPz7fXsWP8GR3dVAjBxxixu/st7mX71It2NQy64aDy0xc5wB4/qUNePM9YRERER6U9Hawt73nmbHX/6A1UfvIcLBskrLmHRp+9i5sLryCoojHWIMoZF5aEtZvYs8Lxz7mB3oZklAVeb2V14D3b5aRS2JSIiImOAv72dvRUb2PnWH9m/uYKA309GXj6XL/sEMxddT35xSaxDFAGik0wvAVYC5WZWCtQDPiAeeAX4J+fce1HYzojy3HPPcccdd1BZWcnMmTM5cOAAZWVlzJgxA+ccaWlpPP7448yYMYONGzeyatUqwLuIYs2aNdxxxx0x3gMREZELy9/Rzv7NFex8az37KjbS1dlBWs445t68lBkLr2XCxTN1IaEMO+eUTJtZvHOuHXgYeDh0a7s8oM05Vx+F+Eas8vJyFi1aRHl5Od/61rcAmDZtGps3bwZg7dq1PPjggzzxxBPMmTOHTZs2kZCQQHV1NXPnzmXZsmUkJETz+lAREZHhx9/ezv7Nm9j59pvsf/cd/B3tpGRmMXvxTcxYeC1FM2cRFxcf6zBFzuhcs7UfmFmyc261mb0KfMk5N+bvgdbc3Mz69et57bXXWLZsWTiZjtTY2EhOTg4Aqamp4fL29nb91S0iIqNaZ3sb+9/bxK631rNv8ya6OjpIycyi7NrrmX7VIibPuoS4eCXQMjKcazJdB3SExr8K/JOZHcC7pV31Oa57xHr++edZsmQJ06dPJzc3l4qKCnJzc9m7dy8f+chHaGpqorW1lQ0bPnxa+oYNG1i5ciVVVVX87Gc/01lpEREZVdqam9hXsZHdG//EgfffJeD3k5qVzezrbmL6VYuYNGu2zkDLiHSuGdsVwKsAzrl3gRvMbDnwkpn9CviBc67tHLcxZH/8j12cOtQc1XXmTU7n2k9OP2ud8vJy7rvvPgBWrFhBeXk5q1ev7tHN4+mnn2bVqlW89NJLAFx55ZVs27aNyspK7rrrLpYuXYrP54tq7CIiIhdSU+0p9m7yEuhD2z7ABYNk5OYz9+alXHzFQibOLFMCLSPeOSXTzrnbzCynezr0FMOdwE+A7wB/ZWZfd8797NzCHDlqa2tZt24dW7ZswcwIBAKYGffee2+Perfddhv33HPPacuXlZWRnp7O1q1bWbDgtMe/i4iIDFvOOWoOH2Tvpg3seectju3dDUDOhCIuv205F1+xkMKpF6k7o4wq59yXwDlXB2BmbwKlwDbgbeBuYAdwn5ld65xbda7bGqz+ziCfD8888wx33nkna9euDZctXryYQ4cO9ai3fv16pk2bBsD+/fuZPHkyCQkJVFVVsWPHDkpKSi5k2CIiIkMSDAQ4urOSve9uZM87b1F/zOvlOf6i6Sxa8RdcdPlVjCuarARaRq1odsxdBWzv4+Esf2NmlVHczrBWXl7OV7/61R5ly5cv56GHHgr3mXbOkZSUxKOPPgp4ifX3vvc9EhMTiYuL4+GHHyYvLy8W4YuIiPSro7WFA++/y96Kjex/bxPtzU3ExSdQPOdSFnz8DqZddiXp43JjHabIBWHn+mBCMyseQLV4oMY513hOG+tlwYIFbtOmTT3KKisrKSsri+ZmRHrQMSYiY1Fd9RH2vbuJfe9u5HDlVoKBAL6MTKbOW8C0y65gyqXzSY64O5XIaGNmFc650/rgRuPM9BOAA872/xuH9wTEJ6OwPRERETnPuvx+DlduZf+777DvvXfC3TfGFU3mso/9GVMvu4KJ02fqAkIZ86LRZ/qGaAQiIiIisdVw4jgH3q9g/+YKDm55H39HOwmJSUyefQnzP3o7U+ctIKtgfKzDFBlWdDNjERGRMaqrs9M7+7y5ggObK6g9ehiAzPxCZl13I1PnX87k2ZeQmKxbtYqciZJpERGRMcI5R82hKg588B5VH7zH4e1b6fJ3Ep+YyORZl3DpzUspnXcZOROKdPcNkQFSMi0iIjKKtdTXcXDr+1R98B4HPniPlrpawOv7fOnNS5gydx6TZ+nss8hQRTWZNrMbnXPruofRXLeIiIj0r7O9jcOVWzm4ZTNVW97n1MEDAPjSM5hyyUeYMnceUy6ZR2ZefmwDFRklon1m+ofA/IjhmPXcc89xxx13UFlZycyZMzlw4ABlZWXMmDED5xxpaWk8/vjjzJgxI7zMwYMHmTVrFmvWrOFLX/pSDKMXEZGRosvvp3r3Dg5t+4CDWz+gevcOgoEA8YmJFM2czbWfuZviOXMpKJ2qO2+InAfnq5vHmO9oVV5ezqJFiygvL+db3/oWANOmTWPz5s0ArF27lgcffJAnnngivMwXv/hFli5dGotwRURkhAh0dXF83x4ved72AUd3VtLV2YFZHAWl01jw8TsovuQjTJxRRmJScqzDFRn11Gf6PGhubmb9+vW89tprLFu2LJxMR2psbCQnJyc8/dxzz1FaWkpaWtqFDFVERIa5cPK8fQuHt2/hyM5K/O1tAOQXl3DpTbcyec5cJpXNxpeWHuNoRcYeJdPnwfPPP8+SJUuYPn06ubm5VFRUkJubG36ceFNTE62trWzYsAHwku/vf//7vPrqq/zwhz+McfQiIhJLXZ2dHNu7iyM7tnNo+xaO7qzE39EOQO6kYmYvvpFJZZcwefYlpGZmxThaERnVyfRrP32EE1X7orrOgilTueHuVWetU15ezn333QfAihUrKC8vZ/Xq1T26eTz99NOsWrWKl156iTVr1vC3f/u3pKfrjIKIyFjT2dbK0V07OLJjG4crt1G9ZycBvx8IJc/X38TkWZcwqWwOqVnZsQ1WRE4T7WS6OTRsivJ6R4za2lrWrVvHli1bMDMCgQBmxr333tuj3m233cY999wDwIYNG3jmmWf4yle+Qn19PXFxcfh8PlavXh2LXRARkfOoqfYUR3dWcmTHdo7s2M7Jqv04F8Ti4igsncZHbvkYk8rmUDRzFikZmbEOV0T6EdVk2jl3XeRwsMzsMeDjwAnn3Jxzjae/M8jnwzPPPMOdd97J2rVrw2WLFy/m0KFDPeqtX7+eadOmAfDHP/4xXL5mzRrS09OVSIuIjALBYIBTB6s4umsHR3dVcnTndhpOHAcgITmZiRfP4MpPfJKiGbOYOH0mSSmpMY5YRAZruHXz+Cnwr8CTMY5jyMrLy/nqV7/ao2z58uU89NBD4T7TzjmSkpJ49NFHYxSliIicD23NTRzbvZOju3dwdGcl1Xt2hS8WTMvOYeKMMuYtuY2iGWXkl0wlPmG4fQ2LyGCZcy7WMfRgZiXAbwZyZnrBggVu06ZNPcoqKyspKys7T9GJ6BgTEU8wEODkwQNU795J9e4dVO/eSV31EQAsLo78KaVMnF7GxOkzmTi9jMz8Aj2iW2QEM7MK59yC3uVR/5PYzNKAdudcINrrFhERiQXnHI0nj1O9ZxfH9uykes9uTuzfS1dnBwCpWdlMuHgms6+/mYkXz6Bw2sUk+VJiHLWIXAjnnEybWRywAvgscDnQASSb2Sngv4C1zrk957qdiO2tAlYBFBcXR2u1IiIiYS31dRzbu5vj+3ZzbO9uju3ZRVtTIwAJiUkUlE7j0ptuZcLFM5hw8UyddRYZw6JxZvo14HfA14BtzrkggJmNA24Avm9mv3bO/TwK28I59wjwCHjdPKKxThERGbtaGxs4vm9P6LWbY/v20FxzCgCzOMYVTWLqZVcw4aLpjJ82nbziEvV1FpGwaHwa3Oyc85vZPzvn7usudM7VAs8Cz5pZYhS2M2DOOZ0hkPNiuF1jICKD01xXy4n9ezm+fw/H93nD7sQZIGdCEZNmzmb8tOkUTruIgpKp6q4hImd1zsm0c84fGm0ys/8EVjjnWszsVuDvnXPXRNQ5KzMrB64H8szsMPBN59y/DyYen89HTU0Nubm5Sqglqpxz1NTU4PP5Yh2KiPTDBYPUnzjGif37OHFgLycO7OPE/r20NtR7FczCiXNh6TQKp15EQek0klPTYhq3iIw8Ufs/lXPufjP7DPC6mXXiPcDla4Ncx6fPNY5JkyZx+PBhTp48ea6rEjmNz+dj0qRJsQ5DRCL4O9o5dbCKk1X7OXFgHyer9nPy4IHwLeni4uPJLZpM6Ucuo6BkKgUl0ygonap7OotIVEQtmTazm4C/AlqACcBK59zOaK1/oBITEyktLb3QmxURkfPMBYM0njrByaoDnDy4n1NVBzh5qIr66qOELtchKSWV/CklzF58I/lTplJYOo3cScUkJCXFOHoRGa2ieQXF3wEPOOfWm9klwNNm9kXn3LoobkNERMaA1sYGTh2s4tShKmoOVXHy0AFqDlXR2dYWrpNdOIG84hJmLryW/CmlFJRMJTO/UF38ROSCimY3jxsjxreY2VK8CxAXRmsbIiIyurQ1N1FzqIqaw4eoOXyQmsNVnDp08MO+zYAvPYO84inMuu5G8qeUkl9cSu7kYl0YKCLDQjTuM22uj1scOOeqQ10/zlhHRERGP+ccrQ311Bw+RO2RQ9QcOUTtkYPUHD5ES31duF6iL4XcSZOZOv9y8iZPIXfyFPKLS0jNytbZZhEZtqJyn2kzexZ43jl3sLvQzJKAq83sLrx7Uf80CtsSEZFhKhgI0HDiGLVHD1N75DC1R49Qe8RLoNtbmsP1klJSGFc0mZK5l5E7uZi8ScXkTi4mIzdfSbOIjDjRSKaXACuBcjMrBeqBFCAOeAX4J+fce1HYjoiIxJhzjramRuqOHqGu+gi11UfC43XVRwkGusJ1U7OyyZlQxPSrF5FbNJlxk4rJLZpM+jjdulRERo9o3Ge6HXgYeDj0cJY8oM05V3+u6xYRkdjoaG2l/thRL0k+dpS66qPUVx+l7thR2pubwvXi4hPILhxPzsQipl52BeMmFDGuaBI5EyeRkp4Rwz0QEbkwotFneiXwC+dcB/BRYDywBfjTua5bRETOD+cc7S3N1B87Sv3xY97wWLX3Ol7d4wJAgIzcfHImTGD6VdcwbuIkciYWkTOhiKz8QuLi42OzEyIiw0A0unnc55x7zMzWADcAbwGfMLMM4BPOuWNR2IaIiAxSMBCgqeYk9ceP0XDimDc8foz649U0HD9GR2tLj/rpuXnkFE5g2mVXkD1+IjnjJ5I9YSLZheNJTNaTP0VE+hKNZLozNPwocLVzLgBgZh/D6/7xiShsQ0REenHBIC31dTScPEHjyeM0nPBejSePUX/8OE01J3HBYLh+XHwCWQUFZBWMZ8LFM8kuHE9W4Xhyxk8kq3A8iUnJMdwbEZGRKRrJ9CEz+ylQgHfhYTOAc+6/zOw7UVi/iMiYFAwEaK6rofHkCRpPnQwNT3jD0HjA7++xTGpWNlkFhUycPpOsguvJKigkq6CQ7MIJpOfmEhenLhkiItEUjWT6bmA58CPgWTN7CdgGzOPDs9YiIhLBOUdHawvNNadoqjnlJcunTtB06iSNp07SVHOSpppTPc4sA6RkZJKZX0B+cQnTFlxJVn4hmQUF3jC/QN0xREQusGjczaMReBzAzP478Hm8BLsO+NS5rl9EZKTpkSjX1tBUc4rm2lM01dSEk+SmmlP429t6LBcXH0/6uDwy8/IpmjmbjNw8MvMKyMwv8IZ5+ST6lCyLiAwnUXucOIQT63+I5jpFRIaTQJeflvo6mmtraa6robm216uulqbaU3R1dPRc0Iy0rGwycvPILZpMyaXzyMjNIz03j4zcfDLz80nLzlE3DBGRESaqybSIyEjl7+ygpa6OlrpaWhpCw/o6mutqvfG6Wprramlrajxt2fiEBNJycknPGUf+lFKmzl9A+rg8L1nOySUjN4+0nBziExJjsGciInI+KZkWkVGry++ntaE+/GppqKO13hu21NfTWl9HS4M37H2bOACLiyMtK5u0nFwyCwqZOKOMtJxxpGWPIz1nHOnjckkfl0tKRqae6CciMkYpmRaREcMFg7Q1N9HW2EBrY4M3bGigtbG+17CB1oY6OlpOT5ABklJSScvOIS07h/wppaTNnUdaljftJcs5pOeM85LkuLgLvJciIjKSKJkWkZhwztHV0UFbcyNtjY20NTWGE+W2pkZv2NhIa1No2NhAe1MTzgX7XJ8vI5PUzCxSs7LILy4hNWsuqVnZpGXnkJqVQ1pWNqlZ2aRmZemOFyIiEjVKpkXknAW6/LQ3N4deTbS3NNHW1OSNNzfR1tRIe1MTbc1NtIeS5vamJrr8Z7h7phkp6RmkZGSSkplFzoQiimbMIjUri5RM75WakfXhdEYm8Qn6OBMRkQtP3z4iAnj9iztamulobaG9uZmOlmbaW1voaA6VtTR/WN7cFJ5ub2k+7RZvkSwuDl96BinpGfgyMsksKKSg9CJSMjO9ZDkjE19GKHFOzyQlMxNferruaiEiIiOCkmmRUSDQ1UVnWyudba10tLbS0drSczw07Ghp8Ybdr5YPh2c8SxySkJRMcloavrR0fOnpZOYXUFAyFV96Oslp6aSke0mwLz0j4pVOcmqaLs4TEZFRS8m0SIwEAwE629vobGsLJcJtdLa34Q8Nw2VtrXS0teEPDbuT5u5kubOtja7Ojn63F5+YSHJqmvdK84YZufnhcV9aemh+Ksnp6aHp9HBCnJCUdAFaRUREZGRRMi3Sj0BXF/6Odu/V3kFXZwf+9nb87W34Ozrwd7TT2d49v3vY5pV1DzvaQ0lyd1kbAb9/QNuPi08gKSWFpJRUb+hLITUzi+zCCSSnppGU6pUnp6SSlJoWGqaGE+OklFSS09JJSNQ9jkVERKJtWCXTZrYE+GcgHnjUOfe9GIckw5RzjmCgi67OztCrg67OTvwdHeHxrs5O/J2h6Y6Iso72cL3Iob+jPVQvlCx3duBv7yAY6BpUbAlJyST6fCT5fCT6Ukj0+UhM9pGWlU2SL4XElFSSfD5v3JcSSpRTQvO8YWTiHJ+YqG4SIiIiw9SwSabNLB74MfDfgMPAO2b2gnNue2wjk764YJCuLj8Bf+jV5aer0xt2l3X5O73pTr9Xt7OzR72uzk4C/k66/H5v2On3Et/wdKe3jtDQH5E4Bzr9Z7xFWn8SEpNISE4mISmJxORkEpKSSUhOJjHZR2pmlpcMJ/u8ecnJJCZ3T/tCiXGv6VCynOTzkZCcrAvnRGREcs7hHLjucQhNe+Uf1utZ1rs+EfMj50UuS3he9zp61Xcf1u9zfniei1hvz9hOi7OfuOgx7/Q26N7Z07YT2X7hOpExnB5z5Hb62r/T4u0d82n7cno8RMbfvU44+3b6en/dWbbTx3Fy5ve6n22csT0il3X4EuP59u1zGE6GTTINXAHscc7tAzCzp4DbgTGRTHsHSpBgV4BgoItAIECwq4tAVxfBQIBAl59gj7LQsMurG/D7vTO1/u4yb36gyx9epvsV7PKfYbq7bmja313mD493bycYCERlv+MSEohPTCI+MZG4hEQSQuPxSUnEJyQSn5xKQnoWqYlJXr2k0Pzu6cQk4pKSek6HhglJScQnJXuvBG/ZuMRELC6ux5dAd/v3nA6/M6fVDTpox9HeXc+BawPaunD4ey3fc33dH8g9vph6bT9y2z2X7bWuM8Tc/eHVe329t33aF0CosM8Ye6yv71jdILd5+rp7f+n1H19f+957W/TRTn3tQ+S+RS56pu313U6nt8NpMfW57bO0X19t38++R7bb2fYh8hjr6z3rsXw/x0D4yzIcQMQXYh8xnPnY6bnu3u9L7wSjr/ajV/mZ1tv7vTh9vT3b/Gzr7JEwnCXW3u/N6evs9b70tc7u6T7avs/3u8/2751IigyNGRhgZqEheGNAeJ5X1rsukdMR46FZWKhSnH24jrTk4XfCajgl00XAoYjpw8CVMYrljJ686yucbN8VxTU6IEjEV/95FId3KMaFxuMwixjHvHGL61mHZLAUIC78YwmR83uvK94rC68nvsf8nmXm7Xpn6AUEQq+zi1hA5DyyXkMRGUn0mzvadCQcAa6PdRg9DKdkekDMbBWwCqC4uPiCb98X5yM5fnLU1hdOYMNJrjdtWEQyGlHWXbdHYmofjltEsmy9l9GHykijd0xERORDqTb87iw1nJLpI0BkljopVNaDc+4R4BGABQsWXPB/Tn3y8W9f6E2KiIiIyDAVF+sAIrwDXGxmpWaWBKwAXohxTCIiIiIiZzRszkw757rMbDXwMl6H2secc9tiHJaIiIiIyBkNm2QawDn3IvBirOMQERERERmI4dTNQ0RERERkRLHe91UdSczsJFAVg03nAadisN2RTG02eGqzwVObDZ7abPDUZoOnNhsctdfgXYg2m+Kcy+9dOKKT6Vgxs03OuQWxjmMkUZsNntps8NRmg6c2Gzy12eCpzQZH7TV4sWwzdfMQERERERkiJdMiIiIiIkOkZHpoHol1ACOQ2mzw1GaDpzYbPLXZ4KnNBk9tNjhqr8GLWZupz7SIiIiIyBDpzLSIiIiIyBApme7FzJaY2U4z22NmX+tjfrKZPR2av8HMSiLmfT1UvtPMbr2ggcfQANrsi2a23cw+MLPfm9mUiHkBM9sceo2Zx8cPoM3uNrOTEW3zlxHz7jKz3aHXXRc28tgZQJv9KKK9dplZfcS8MXecmdljZnbCzLaeYb6Z2f8NtecHZjY/Yt5YPcb6a7PPhtpqi5n9yczmRsw7ECrfbGabLlzUsTWANrvezBoifv/+PmLeWX+nR6MBtNeXI9pqa+iza1xo3lg9xiab2WuhPGKbmd3XR53Yfp455/QKvfAeY74XmAokAe8Ds3rV+QLwb6HxFcDTofFZofrJQGloPfGx3qdh0mY3AKmh8b/ubrPQdHOs92GYttndwL/2sew4YF9omBMaz4n1Pg2HNutV/2+AxyKmx+Jxdh0wH9h6hvkfBX4LGHAVsCFUPiaPsQG22cLutgCWdrdZaPoAkBfrfRiGbXY98Js+ygf1Oz1aXv21V6+6y4B1EdNj9RibAMwPjWcAu/r4zozp55nOTPd0BbDHObfPOdcJPAXc3qvO7cATofFngJvMzELlTznnOpxz+4E9ofWNdv22mXPuNedca2jybWDSBY5xuBnIcXYmtwKvOudqnXN1wKvAkvMU53Ay2Db7NFB+QSIbppxzfwBqz1LlduBJ53kbyDazCYzdY6zfNnPO/SnUJqDPMmBAx9mZnMvn4Ig1yPYa859jAM65aufcu6HxJqASKOpVLaafZ0qmeyoCDkVMH+b0NyxcxznXBTQAuQNcdjQa7H7/D7y/Hrv5zGyTmb1tZn92HuIbjgbaZstD/656xswmD3LZ0WbA+x3qRlQKrIsoHovHWX/O1KZj9RgbrN6fZQ54xcwqzGxVjGIarq42s/fN7LdmNjtUpuPsLMwsFS/pezaieMwfY+Z1rZ0HbOg1K6afZwnRXqHImZjZ54AFwOKI4inOuSNmNhVYZ2ZbnHN7YxPhsPKfQLlzrsPMPo/335AbYxzTSLECeMY5F4go03EmUWNmN+Al04siiheFjrEC4FUz2xE6CznWvYv3+9dsZh8FngMujm1II8Iy4E3nXORZ7DF9jJlZOt4fF//LOdcY63gi6cx0T0eAyRHTk0JlfdYxswQgC6gZ4LKj0YD228xuBv4OuM0519Fd7pw7EhruA17H+4tztOu3zZxzNRHt9Chw2UCXHaUGs98r6PWv0TF6nPXnTG06Vo+xATGzS/F+J293ztV0l0ccYyeAXzM2uvn1yznX6JxrDo2/CCSaWR46zvpzts+xMXeMmVkiXiL9C+fcr/qoEtPPMyXTPb0DXGxmpWaWhHcw977y/wWg+2rQP8e7OMCFyleYd7ePUry/vDdeoLhjqd82M7N5wFq8RPpERHmOmSWHxvOAa4DtFyzy2BlIm02ImLwNr48YwMvALaG2ywFuCZWNdgP53cTMZuJdZPJWRNlYPc768wLwF6Gr4K8CGpxz1YzdY6xfZlYM/Aq40zm3K6I8zcwyusfx2qzPuzWMNWY2PnRdEWZ2BV7eUcMAf6fHIjPLwvsP7vMRZWP2GAsdP/8OVDrn/vEM1WL6eaZuHhGcc11mthqvoePx7gawzcy+DWxyzr2A94b+zMz24F1EsCK07DYz+w+8L+ku4N5e/2YelQbYZv8ApAP/L/SZetA5dxtQBqw1syDeB+z3nHOjPskZYJv9TzO7De9YqsW7uwfOuVoz+z94X0QA3+71b8BRaYBtBt7v41OhP3C7jcnjzMzK8e6kkGdmh4FvAokAzrl/A17EuwJ+D9AK3BOaNyaPMRhQm/093jUyD4c+y7qccwuAQuDXobIE4JfOuZcu+A7EwADa7M+BvzazLqANWBH6/ezzdzoGu3BBDaC9AO4AXnHOtUQsOmaPMbwTIHcCW8xsc6jsG0AxDI/PMz0BUURERERkiNTNQ0RERERkiJRMi4iIiIgMkZJpEREREZEhUjItIiIiIjJESqZFRERERIZIybSIiIiIyBApmRYRERERGSIl0yIiY4CZXW5mH5iZL/Q0tW1mNifWcYmIjHR6aIuIyBhhZt8BfEAKcNg591CMQxIRGfGUTIuIjBFmloT3WN12YKFzLhDjkERERjx18xARGTtygXQgA+8MtYiInCOdmRYRGSPM7AXgKaAUmOCcWx3jkERERryEWAcgIiLnn5n9BeB3zv3SzOKBP5nZjc65dbGOTURkJNOZaRERERGRIVKfaRERERGRIVIyLSIiIiIyREqmRURERESGSMm0iIiIiMgQKZkWERERERkiJdMiIiIiIkOkZFpEREREZIiUTIuIiIiIDNH/B4Cb54k2bhN4AAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ - "# YOUR CODE HERE\n", - "raise NotImplementedError()" + "algos = [\n", + " (phi_euler, \"Euler\"),\n", + " (phi_euler_modified, \"modified Euler (Collatz)\"),\n", + " (phi_heun, \"Heun\"),\n", + " (phi_rk4, \"4th order Runge-Kutta\"),\n", + " (phi_ab3, \"AB3\"),\n", + " (phi_ab4, \"AB4\"),\n", + "]\n", + "\n", + "def ODEF(x, y):\n", + " return y - x**2 + 1.0\n", + "\n", + "h = .02\n", + "x = np.arange(0, 2, h)\n", + "y = (x + 1)**2 - 0.5*np.exp(x)\n", + "eta = np.zeros(len(x))\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", + "\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", + "\n", + "ax0.legend()\n", + "plt.xlabel(\"x\")\n", + "ax0.set_ylabel(\"y(x)\")\n", + "ax1.legend()\n", + "plt.show()" ] }, {