Files
ANum/QUESTION18.ipynb
2024-12-03 20:04:03 +01:00

148 lines
55 KiB
Plaintext

{
"cells": [
{
"cell_type": "code",
"execution_count": 2,
"id": "23bf2de4-1ad7-4bcc-b8a4-75737b31ec3a",
"metadata": {
"tags": []
},
"outputs": [],
"source": [
"import numpy as np\n",
"import matplotlib.pyplot as plt"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "87e8765d-2b61-478d-9ad1-2a8d4cb925a3",
"metadata": {
"tags": []
},
"outputs": [],
"source": [
"def f_edo(t, y):\n",
" return 1/(1+y)"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "b342537d-ad37-42db-aef3-b33ac342695e",
"metadata": {
"tags": []
},
"outputs": [],
"source": [
"def ptfixe_eulerimp(fct, x0, h, ti, yp, maxiter=10):\n",
" x = x0\n",
" for k in range(maxiter):\n",
" x = yp + h * fct(ti, x)\n",
" return x"
]
},
{
"cell_type": "code",
"execution_count": 24,
"id": "be078e7a-bfaf-4b71-82c7-3609dfed1ad3",
"metadata": {
"tags": []
},
"outputs": [],
"source": [
"def eulerimplicite(fct, a, b, y0, h):\n",
" t = np.arange(a, b, h)\n",
" y = np.zeros(len(t))\n",
" \n",
" y[0] = y0\n",
" \n",
" for i in range(1, len(t)):\n",
" y[i] = ptfixe_eulerimp(fct, y[i-1], h, t[i], y[i-1]) \n",
" y[i] = y[i-1] + h*fct(t[i], y[i])\n",
" \n",
" return t,y"
]
},
{
"cell_type": "code",
"execution_count": 29,
"id": "fcd619b9-0f98-4dfc-9413-58f430c5991f",
"metadata": {
"tags": []
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAhYAAAGdCAYAAABO2DpVAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABJNUlEQVR4nO3dd3QUZcPG4d+WJARSIIGQIAFC771LVVqkdxWVrigiyAsiooKKgCIKiiKgUhQVG4gFFJAqIogiiKiAKCgl0hIIkGRn5vuDl3xvpCVhk9kk93VOzjGb3Z0748DezMzzPA7LsixEREREvMBpdwARERHJPVQsRERExGtULERERMRrVCxERETEa1QsRERExGtULERERMRrVCxERETEa1QsRERExGvc2b1B0zQ5dOgQwcHBOByO7N68iIiIZIJlWZw+fZpixYrhdF75vES2F4tDhw4RHR2d3ZsVERERLzh48CDFixe/4s+zvVgEBwcDF4KFhIRk9+ZFREQkExISEoiOjk79HL+SbC8WFy9/hISEqFiIiIjkMNe6jUE3b4qIiIjXqFiIiIiI16hYiIiIiNdk+z0W6WEYBikpKXbHkBzO5XLhdrs1rFlEJBv5XLE4c+YMf/31F5Zl2R1FcoH8+fMTFRWFv7+/3VFERPIEnyoWhmHw119/kT9/fooUKaJ/aUqmWZZFcnIy//zzD/v376dcuXJXndBFRES8w6eKRUpKCpZlUaRIEQIDA+2OIzlcYGAgfn5+/PnnnyQnJ5MvXz67I4mI5Ho++U84nakQb9FZChGR7KW/dUVERMRrVCx8xNq1a3E4HJw6dcruKJk2f/58ChYsmPr9hAkTqFmzZrpf73A4WLp0qddziYhI9lGx8IJ+/frhcDgu+WrXrp3d0YAL+bp06ZLt2x01ahSrV69O9/MPHz5MbGwsAH/88QcOh4Pt27dnUToREckKPnXzZk7Wrl075s2bl+axgICAbM2QnJzsU8Mqg4KCCAoKSvfzIyMjszCNiIhkB52x8JKAgAAiIyPTfBUqVAi4/L++T506hcPhYO3atVd8z02bNtGsWTMCAwOJjo7mgQceIDExMfXnpUqVYuLEifTr14/Q0FAGDx6crqwtWrRg2LBhjBgxgkKFClG0aFHmzJlDYmIi/fv3Jzg4mDJlyrB8+fLU11y8VPPZZ59Ro0YN8uXLR4MGDdi5c+cVt3O5SyFvvPEGVapUISAggKioKO6///7Un/3vpZCYmBgAatWqhcPhoEWLFqnPmzdvHpUqVSJfvnxUrFiRV155JV2/t4hIrpZ0Gr6dg/HhPbbG8OliYVkWZ5M9tnzZPUHXzp07adu2Ld26dWPHjh0sXryYjRs3pvkgBpg6dSpVq1Zl27ZtPPbYY+l+/wULFlC4cGG2bNnCsGHDuPfee+nZsyeNGzfm+++/p23bttx5552cPXs2zetGjx7Nc889x9atW4mIiKBTp07pniV11qxZDB06lLvvvpudO3eybNkyypYte9nnbtmyBYBVq1Zx+PBhPvroIwDmzp3LuHHjePrpp9m9ezeTJk3iscceY8GCBen+3UVEcpXj+7A+fwjPcxVh+WhcO9/F+PsH2+L49KWQcykGlR//wpZt//xkW/L7p3/3fPrpp5ec9h8zZkyGPuz/19SpU7n99tsZMWIEAOXKlePFF1+kefPmzJo1K3VOhptuuolRo0Zl+P1r1KjBo48+CsDYsWOZMmUKhQsXTj3r8fjjjzNr1ix27NhBw4YNU183fvx4WrduDVwoJ8WLF2fJkiX06tXrmtucOHEi//nPfxg+fHjqY/Xq1bvsc4sUKQJAeHh4mkskTz31FNOmTaNbt27AhTMbP//8M7Nnz6Zv374Z2QUiIjmXacK+rzA2v4pr30ocXPhA32dGMc9oR9ezhaljUzSfLhY5ScuWLZk1a1aax8LCwjL9ftu2bWPv3r0sWrQo9THLsjBNk/3791OpUiUA6tatm6n3r169eup/u1wuwsPDqVatWupjRYsWBSAuLi7N6xo1apT632FhYVSoUIHdu3dfc3txcXEcOnSIm2++OVN5Af755x8OHjzIwIED01z28Xg8hIaGZvp9RURyjKTTsP0dUr55Fb9T+3ABpuVgjVmTt4klslY77mpcmgqRwbZFzHCx+PvvvxkzZgzLly/n3LlzlC9fntdff506dbzfjQL9XPz8ZFuvv296t50RBQoUuOJp/YuTNP3v5ZVrXT4wTZN77rmHBx544JKflShRIs12M8PPzy/N9w6HI81jFycpM03zmu+VngnNvDGT6sUsc+fOpUGDBml+5nJl7P+XiEiOcnwf1rezMX5YhDvlDH5AghXI+0YLVgZ15ObGjXi+bjSh+f2u+VZZLUPF4uTJk9x44420bNmS5cuXExERwb59+9LMXeBNDocjQ5cjfNXF0/qHDx+mVq1aANccRlm7dm127dp1xbJil82bN6cWm5MnT/Lbb79RsWLFa74uODiYUqVKsXr1alq2bHnN518c3WIYRupjRYsW5YYbbuD333+nT58+mfwNRERyCNOE37/C+ObSyx3zjbYcKtmZW5tUYVHFCFxO35mxOkOf2s888wzR0dFphlWWKlXK25lypKSkJI4cOZLmMbfbTeHChQkMDKRhw4ZMmTKFUqVKcezYsdT7G65kzJgxNGzYkKFDhzJ48GAKFCjA7t27WblyJS+99FJW/ipX9eSTTxIeHk7RokUZN24chQsXTvccGRMmTGDIkCFEREQQGxvL6dOn+frrrxk2bNglz42IiCAwMJAVK1ZQvHhx8uXLR2hoKBMmTOCBBx4gJCSE2NhYkpKS+O677zh58iQjR4708m8rImKDy1zuAFht1OJtYomq1Y67bixN+aL2Xe64mgyNClm2bBl169alZ8+eREREUKtWLebOnXvV1yQlJZGQkJDmKzdasWIFUVFRab6aNGmS+vM33niDlJQU6taty/Dhw5k4ceJV36969eqsW7eOPXv20LRpU2rVqsVjjz1GVFRUVv8qVzVlyhSGDx9OnTp1OHz4MMuWLUv33Bl9+/Zl+vTpvPLKK1SpUoUOHTqwZ8+eyz7X7Xbz4osvMnv2bIoVK0bnzp0BGDRoEK+99hrz58+nWrVqNG/enPnz56cOTxURybFO/I61fAzGf0d3+J3aR4IVyBuedtyW72X2t3mD58eOZGK3Gj5bKgAcVgbGVV4ciTBy5Eh69uzJli1bGDFiBLNnz+auu+667GsmTJjAE088ccnj8fHxhISEpHns/Pnz7N+/n5iYGK1E6WPWrl1Ly5YtOXnyZJZd+soKOqZExKdZFvyxEfObV3D8thwHFz6S/3254yYfuNyRkJBAaGjoZT+//1eGLoWYpkndunWZNGkScGHyol27djFr1qwrFouxY8emOUWdkJBAdHR0RjYrIiKSu6Sch58+xLPpZdz/7Eq9fLDGqMFbxBJRI5a+TUpTMfLKH+C+KkPFIioqisqVK6d5rFKlSnz44YdXfE1AQEC2T20tIiLik87EwdbX8Wx5Dfe5Y7iBc5Y/HxjNWJavI82b3ciz9UsQHpRzPzczVCxuvPFGfv311zSP/fbbb5QsWdKrocT3tGjRwvbZSEVEcqzDO7A2z8La+T5OMwU3cMgKY6GnDTuKdqF3s+osqhqFv9unJ8ROlwwViwcffJDGjRszadIkevXqxZYtW5gzZw5z5szJqnwiIiI5k2nAbyvwbHoF94GNOAAH8L1ZlvnmLVgVO9KvaVnGlCiUrvmAcooMFYt69eqxZMkSxo4dy5NPPklMTAzTp0/XnAIiIiIXJZ2GH94iZdMs/BL+xA14LCefmw14z9WeKo1aMaZRKW4oeP0TB/qiDM8+1aFDBzp06JAVWURERHKuk39cmB1z28LU2TFPWQV427iZdaGd6NC0PnNq35ArJn68mtz924mIiGQly4I/N2F88wrOXz/HgYkb2GsWY57RjqMxXbijaSWGlCuC04dmx8xKKhYiIiIZ5UmCnz66MFw0bmfq7JjrjOq8yS0UrXUL/W4sTTkfnsgqq6hYiIiIpNfZE/DdG3g2z8Z9Ni51uOhHRlM+DexEsxZNea5+NAXzp29G4txIxcJHTJgwgaVLl15zcbKs5q0ZNvv168epU6dYunSp17KJiNjm+D6szbMwv38Ll3EON3D4v8NFf4zoTO/mNVhYLQo/V84fLnq9VCy8IC4ujscee4zly5dz9OhRChUqRI0aNZgwYQKNGjWyO16W+uOPP4iJieGHH36gZs2aqY/PmDFD816ISM5mWXDwW8yvX8Lx62c4sHABu8ySvGa052y5TgxoVp6HYsJy1XDR66Vi4QXdu3cnJSWFBQsWULp0aY4ePcrq1as5ceKE3dFsExoaancEEZHMMTzwyyd4Nr6E+/C21Om2vzJqsoAORNduy7AmpSldJMjWmL5K52yu06lTp9i4cSPPPPMMLVu2pGTJktSvX5+xY8fSvn371OcdOHCAzp07ExQUREhICL169eLo0aPp3o5hGAwcOJCYmBgCAwOpUKECM2bMSPOcfv360aVLF5577jmioqIIDw9n6NChpKSkpD7nrbfeom7dugQHBxMZGcntt99OXFzcZbeZmJhISEgIH3zwQZrHP/nkEwoUKMDp06dTVxWtVasWDoeDFi1apMnyv+911113ERQURFRUFNOmTaNFixaMGDEi9TkOh+OSSycFCxZk/vz5qd///fff9O7dm0KFChEeHk7nzp35448/0rcTRUSuJuk0bJ6FZ0ZNeL8f7sPbSLL8eMfTkp6u6fzc8jVeeHg4E7tWV6m4Ct8+Y2FZkHLWnm375Yd0nNoKCgoiKCiIpUuX0rBhw8uui2JZFl26dKFAgQKsW7cOj8fDfffdR+/evVm7dm264pimSfHixXnvvfcoXLgwmzZt4u677yYqKopevXqlPm/NmjVERUWxZs0a9u7dS+/evalZsyaDBw8GIDk5maeeeooKFSoQFxfHgw8+SL9+/fj8888v2WaBAgW49dZbmTdvHj169Eh9/OL3wcHBbNmyhfr167Nq1SqqVKlyxSXUR48ezZo1a1iyZAmRkZE88sgjbNu2Lc3lk2s5e/YsLVu2pGnTpqxfvx63283EiRNp164dO3bsSPfy7SIiacT/DVtm49n6Bu7k07iB41Ywbxmt2VSoC92b1+atmsUIcLuu+Vbi68Ui5SxMKmbPth85BP4Frvk0t9vN/PnzGTx4MK+++iq1a9emefPm3HrrrVSvXh2AVatWsWPHDvbv35+6suubb75JlSpV2Lp1K/Xq1bvmdvz8/NIsPx8TE8OmTZt477330hSLQoUKMXPmTFwuFxUrVqR9+/asXr06tVgMGDAg9bmlS5fmxRdfpH79+pw5c4agoEsb+KBBg2jcuDGHDh2iWLFiHDt2jE8//ZSVK1cCUKRIEQDCw8OJjIy8bPYzZ87w+uuvs3DhQlq3bg3AggULKF68+DV/7//17rvv4nQ6ee2111KvZ86bN4+CBQuydu1a2rRpk6H3E5E87vAOzE0vwU8f4bQ8uLmwXPnrxi0cKdWFfs0r8UC5wrp/IoN0KcQLunfvzqFDh1i2bBlt27Zl7dq11K5dO/UU/u7du4mOjk6zXHzlypUpWLAgu3fvTvd2Xn31VerWrUuRIkUICgpi7ty5HDhwIM1zqlSpgsv1/606KioqzaWOH374gc6dO1OyZEmCg4NTL138+30uql+/PlWqVGHhwoXAhUJUokQJmjVrlu7c+/btIzk5Oc2NrGFhYVSoUCHd7wGwbds29u7dS3BwcOqZorCwMM6fP8++ffsy9F4ikkeZJvz2Jca8DjC7Kc6d7+G0PGw2K3GPZzSzq75D32FP8MbgZjQrX0SlIhN8+4yFX/4LZw7s2nYG5MuXj9atW9O6dWsef/xxBg0axPjx4+nXrx+WZV324LzS45fz3nvv8eCDDzJt2jQaNWpEcHAwU6dO5dtvv00b288vzfcOhwPTNIEL9zm0adOGNm3a8NZbb1GkSBEOHDhA27ZtSU5OvuK2Bw0axMyZM3n44YeZN28e/fv3z9AftvSODnE4HJc893/vDzFNkzp16rBo0aJLXnvxzImIyGWlnIcdi/F8PRP3id9w8f/rd7zr7kydJjfxVKOSRATnsztpjufbxcLhSNflCF9UuXLl1BsRK1euzIEDBzh48GDqWYuff/6Z+Ph4KlWqlK7327BhA40bN+a+++5LfSyj/0r/5ZdfOHbsGFOmTEnN8d13313zdXfccQcPPfQQL774Irt27aJv376pP7t4X4NhGFd8fdmyZfHz82Pz5s2UKFECgJMnT/Lbb7/RvHnz1OcVKVKEw4cPp36/Z88ezp79/3tsateuzeLFi4mIiCAkJCSdv7WI5GmJx2Hra3i+nYP73DHcwGkrkHeNlqwO6Ur7Zg14vXZxAv11/4S3+HaxyAGOHz9Oz549GTBgANWrVyc4OJjvvvuOZ599ls6dOwPQqlUrqlevTp8+fZg+fXrqzZvNmzenbt266dpO2bJlWbhwIV988QUxMTG8+eabbN26NXVURnqUKFECf39/XnrpJYYMGcJPP/3EU089dc3XFSpUiG7dujF69GjatGmT5t6IiIgIAgMDWbFiBcWLFydfvnyXDDUNCgpi4MCBjB49mvDwcIoWLcq4ceNwOtNeibvpppuYOXMmDRs2xDRNxowZk+YMTJ8+fZg6dSqdO3fmySefpHjx4hw4cICPPvqI0aNHZ/ieDRHJxU7+gbVpJub3b+IyzuMG/rbCmedpx57i3bijeTXerhiRZ9bvyE66x+I6BQUF0aBBA1544QWaNWtG1apVeeyxxxg8eDAzZ84E/n8YZaFChWjWrBmtWrWidOnSLF68ON3bGTJkCN26daN37940aNCA48ePpzl7kR5FihRh/vz5vP/++1SuXJkpU6bw3HPPpeu1AwcOJDk5Oc3Nn3Dh5tUXX3yR2bNnU6xYsdQy9W9Tp06lWbNmdOrUiVatWtGkSRPq1KmT5jnTpk0jOjqaZs2acfvttzNq1Cjy5///S1L58+dn/fr1lChRgm7dulGpUiUGDBjAuXPndAZDRC44/CPm+wOwZtTCsXUuLuM8P5mleCDlfiaXe4cOQyax4N5WtK5cVKUiizisbJ4eMSEhgdDQUOLj4y/5MDh//jz79+8nJiaGfPl0ncuXLFq0iOHDh3Po0CGvDets0aIFNWvWZPr06V55v8vRMSWSB1gW/L4WY+N0XPvXpj683qjG63QiunY7BjUtQ6nCOfPSuq+42uf3/9KlELmqs2fPsn//fiZPnsw999yjuSJExHcYHvh5KZ4N01NXGPVYTj41G/KOuwsNmrZgWuNSFA66dH4hyToqFnJVzz77LE8//TTNmjVj7NixdscREYHkRPhhEZ6vX8SdcBA3cNYKYLHRgk/zd6Nj8wbMqxdNfn99xNlBl0IkV9MxJZKLJB6HLXMuLFmedBK4MEPmfE9bthbpxu0ta3FL1UjcWmE0S+hSiIiI5A6XGeHxpxnBXKM9h0t1pX+LKowsG67JrHyEioWIiPimQ9sxv56BY9dSHJi4gB1mDHPNjrirdGJQ8/JUKaaVlH2NTxaLbL46I7mYjiWRHMay4Pc1GBum4/pjXeqcCOuM6sx3dKZU3XaMaVqa4oUyNjuyZB+fKhYX17hITk4mMDDQ5jSSG1ycufPfU52LiI8xDdj9Ccb6abiO7kgd4fGJ2Yj3/LtyY7MWvNCwJAXza2Sar/OpYuF2u8mfPz///PMPfn5+l8zMKJJelmVx9uxZ4uLiKFiwYJqF2UTEh3iSL6zhseEF3Cf34eLCCI93jZZ8VqArXVo2Yl6d4uTz05/hnMKnioXD4SAqKor9+/fz559/2h1HcoGCBQtecTl3EbFRciJ8vxDPxhdxnzmEGzhlFWCB0ZaNhbrR56baLK4epREeOZBPFQu4sKhVuXLlrrrapkh6+Pn56UyFiK85dxK2zMXzzSu4z5/EDRy1CjLX056fo7oy4KbqDNMaHjmazxULAKfTqTkHRERyk9NH4JuZGFvewOVJxA38YRZlttGBuNLduLtlJcbFhGnIaC7gk8VCRERyiRP7sb6egfXDWzjNFFzAbrMEs4xOWJW7cE+L8lS9QUNGcxMVCxER8b4jP2FueB7HriU4MHEAW8wKzDW7ULhWex5sXpYYLQqWK6lYiIiI9xz49sKQ0b1fpM5BscaowRuOrlRs0JaJTUtTNESXunMzFQsREbk+lgX7VmOsm4br4CZcgGE5+NxswFvubtzY4iZeaqQ5KPIKFQsREckc04RfPsFY91zqpFbJlosPjWZ8FNid2OZNmFdfq4zmNfq/LSIiGWN4YNcSPOum4j7+Ky4g0QrgbeNmvgzpQc+b6rOo5g34uzUHRV6kYiEiIuljpFyYJXPdc7hP7ccNJFj5mWe0Y0Oh7tzVqjbvVovCpTko8jQVCxERubqU87D9LTzrX8B9+i/cwAkriNc8t7AtogcDWtVgWKWimtRKABULERG5kuSzsG0+no0zcCcewQ38Y4Uy29OBn2/owd03V2V0+SKa1ErSULEQEZG0kk7D1tfwbHwJ9/njuIFDVhivejpyoFR37rmpKuNKa5ZMuTwVCxERueDcSfh2Dp5vXsadFI8bOGAW4RWjM8fLdmfITZWoU7KQ3SnFx6lYiIjkdYnHYfPLGJvn4Eo5jRvYZ0bxstGF5ErdGNKygqbdlnRTsRARyatOH4FNL2FsfR2X59x/1/GI5mWjK/7VunBvy/KUKxpsd0rJYVQsRETymoRDsHE65nfzcJrJuIAdZgyzzK4UrNmZUS3KUUrreEgmqViIiOQVCYewNr6A9d18nGYyTmCbWY5XzG5E1u7AozeV44aCgXanlBxOxUJEJLdLOIS14XmsbQtwmsmpK42+bPagRJ1YnmpZlmIqFOIlKhYiIrnVZQrFt2ZFXjZ7UKpuO6a0LEtUqAqFeJeKhYhIbhP/94VLHtvm4zRTUgvFTKsnZeq249kWZYkM1dLlkjUyVCwmTJjAE088keaxokWLcuTIEa+GEhGRTLhCoXjZ6knpuu14rmVZioaoUEjWyvAZiypVqrBq1arU710ul1cDiYhIBl2lUJStF8tzLcoQoUIh2STDxcLtdhMZGZkVWUREJCPi/8La8ALW9wsuKRTlGtzCc81LExGsQiHZK8PFYs+ePRQrVoyAgAAaNGjApEmTKF269BWfn5SURFJSUur3CQkJmUsqIiIXXKZQbDYr8YrVgwoNb2FaszIUCQ6wO6XkURkqFg0aNGDhwoWUL1+eo0ePMnHiRBo3bsyuXbsIDw+/7GsmT558yX0ZIiKSCQmHsNZPw/p+Yeooj2+MysyiBxUb3cK0pqVVKMR2DsuyrMy+ODExkTJlyvDQQw8xcuTIyz7ncmcsoqOjiY+PJyQkJLObFhHJO878AxtfwNwyF6eZDFwoFK/QnYoNb+Ge5mUoHKRCIVkrISGB0NDQa35+X9dw0wIFClCtWjX27NlzxecEBAQQEKADXkQkw86egK9nYHw7G5fnHE4uTGz1ktmLMvVjmaabMsUHXVexSEpKYvfu3TRt2tRbeURE5Hw8fPMyxqaXcaWcwQVsN0szw+xFsdrtefbmcprYSnxWhorFqFGj6NixIyVKlCAuLo6JEyeSkJBA3759syqfiEjekXQGvn0V4+sXcSXF4wJ+NkvygtGTsJqdePLmckSH5bc7pchVZahY/PXXX9x2220cO3aMIkWK0LBhQzZv3kzJkiWzKp+ISO6XfBa+ex3P+udxnz+BC9hj3sB0owf5qnVhXKsKWm1UcowMFYt33303q3KIiOQ9niTYNh/Puudwn43DDew3izLD6I5RuRsPtq5E2Yggu1OKZIjWChERyW5GCmxfhGfNM7jPHMIN/GUVZoanG4kVevJAm4pUjNSoOcmZVCxERLKL4YGd7+FZMwV3/J+4gSNWIWZ6uhBXthcPtKlC1RtC7U4pcl1ULEREspppws9LML6ahOvEXtzAP1YIszyd+TOmN8PaVqNmdEG7U4p4hYqFiEhWsSzYuwpj5RO44nbiAk5aQcz2dGBX8Vt5ILYG9UqF2Z1SxKtULEREssKBzZgrJ+A8+A0u4LQVyGueW9gU0Zv7Y2szplxhHA6H3SlFvE7FQkTEm478hLn6CZx7vsQJJFl+LDRa81nobdzdrh4jqkaqUEiupmIhIuINJ37H+moS/PQBTiw8lpP3jea8V+B2bmvdiA9q3YDb5bQ7pUiWU7EQEbkeCYex1j17YcVRywPAp0ZD3vC7jY5tmvNugxIEuF02hxTJPioWIiKZcfYEfD0dY/NsXMZ5HMA6ozovO2+nWctWvHljDAUC9Fes5D066kVEMiI5ETbPwtg4A1dyAi5gm1mO6dZtVGncnjnNS1Mwv7/dKUVso2IhIpIenuT/Tr/9LO6z/+ACdpvRvGD0pmjdLky7uZyWMBdBxUJE5OpME376EGP1k7jiD+AG/jQjeMHogbNaDx5tXYkS4VpxVOQiFQsRkSvZtwZj5XhcR37EBcRZBXnR05UjZXoy6pZqWs9D5DJULERE/u3wDsyV43H+/lXq5FavejqyNepWHoytSaMy4XYnFPFZKhYiIhedOoC1+inY+T5OLJItF4uMVnwc2od72tVnlCa3ErkmFQsRkbMnYMM0zG/n4DSTAVhmNOJ1/zvoGduU9+tF46fJrUTSRcVCRPKulHPw7asY65/HlZyAE/jaqMKLzjto0qI1bzfRXBQiGaU/MSKS95gG/PgOntVP4z5z6L9DR0sw1byNEvU68crN5QgPCrA7pUiOpGIhInmHZcGeL/F8OR73sd24gb+tcKal9MSs1pMJbSpr6KjIdVKxEJG84a9tmF8+hvPA17iBU1YBXvZ0Zl/M7YyMrU7VG0LtTiiSK6hYiEjudvJPrFUTcOz6KHUZ83lGW9YUvoP729dlXLkidicUyVVULEQkdzoff2GkxzezcJrJmJaDj8ymvJmvD3e1a8I7tW7A6dTQURFvU7EQkdzFSIFt8zG+moTr/InUkR7PO+6iZcubebdJaQL9tYy5SFZRsRCR3MGy4LcvML54FNeJPbiAvWYxphi3E1GnM6+2rkCRYI30EMlqKhYikvMd3oHxxThcf6zHBRy3gnnB04MjZXoxpn01yhUNtjuhSJ6hYiEiOVfCIayvnoLt7+DCIsly84YRy6rCd/Bg+7o0KVfY7oQieY6KhYjkPMmJ8PWLGF/PwOU5B1yYgvuNfHfRp21T3qtdHJduzBSxhYqFiOQcpgHb38az6incZ4/iAraZ5ZjKXTRuEcvbTWPI76+/1kTspD+BIpIz/L4WY8U4XHE/4QYOmEV4xriNkNo9ebFNeSKC89mdUERQsRARX3dsL9YXj+DY8wUuIMHKz4ueruwvfTsPta9BhUjdmCniS1QsRMQ3nY+Hdc9ifjsbp5lCiuXiLaMVy0Lv4IFODXm0QoTdCUXkMlQsRMS3mAb88BbGqidwnTuOE/jKqMkMVz+6xLbgvYYl8XM57U4pIlegYiEivuOPrzE+H4MrbicuYJ8ZxSTjTorX78z8VuUpVMDf7oQicg0qFiJiv1MHLqw8+vPS1Psopnu683vMbTzSsTrlNcGVSI6hYiEi9klOhI3TMb+egdNIwrAcvGPcxIehfRnWsSGPVYjA4dB8FCI5iYqFiGQ/y4Kd7+P54nHciYdxAt8YlXnO1Y/Ytq1Z3KgU/m7dRyGSE6lYiEj2+msbxvKHcP39HW7goFmESUYfwut2Z26bioTpPgqRHE3FQkSyx+mjWKvG4/jxHVxAohXAy54u7C51Jw93rKn5KERyCRULEclaRgpsmYPx1SRcKWcA+NBoysIC/bivYxNGVy6q+yhEchEVCxHJOr+vw/hsNK7jv+ICtpuledrsT5MWsSxuXpp8fi67E4qIl6lYiIj3xf+F9cWjOH5eggs4bgXzjOdWzlTszQsdqlC8UH67E4pIFlGxEBHv8STBNzMx1k3F5TmHYTl402jN0kL9GNWpAU3KFbY7oYhkMRULEfGOPSvxfDYa96n9uIAtZgWecQwktm1r3m9cStNwi+QRKhYicn1O7MdcMRbnb8txA3FWQZ5OuR13jd7MuqWiljMXyWOu658QkydPxuFwMGLECC/FEZEcI/ksrJmEObM+zt+Wk2K5mO1pz/1hc7nrnoeY1rumSoVIHpTpMxZbt25lzpw5VK9e3Zt5RMTXWRb88inG8rG4Eg7iBDYaVXjePYie7VvxTt1oXE4NHxXJqzJVLM6cOUOfPn2YO3cuEydO9HYmEfFVx/dhfj4a577VuIC/rXCe9txJ4Xo9eKNNBQrm16yZInldporF0KFDad++Pa1atbpmsUhKSiIpKSn1+4SEhMxsUkTslHIevp6BueE5nEYySZabOUYHNkTexeNd61L1hlC7E4qIj8hwsXj33Xf5/vvv2bp1a7qeP3nyZJ544okMBxMRH7HvK4xP/4Pr5O84gfVGNaa6BnNbx5a8Wy8apy57iMj/yFCxOHjwIMOHD+fLL78kX7703ZQ1duxYRo4cmfp9QkIC0dHRGUspItkv4TDWF+Nw7PoQF3DUKshTKXcSWLMH82+pRHhQgN0JRcQHOSzLstL75KVLl9K1a1dcrv+fhtcwDBwOB06nk6SkpDQ/u5yEhARCQ0OJj48nJCQk88lFJGuYBmyZi7H6KVwpZzAsBwuNNnwS3p9HujagbqkwuxOKiA3S+/mdoTMWN998Mzt37kzzWP/+/alYsSJjxoy5ZqkQER/31zaMT4bjOrrzv2t7lGEig2jXth2LNcmViKRDhopFcHAwVatWTfNYgQIFCA8Pv+RxEclBzp3EWvUkbJuHC4t4Kz/Pem4lvvLtzOxQjchQzUchIumjmTdF8jLLgh2L8awYh/vcMQA+NJqwKHgwI7rcSLPyRWwOKCI5zXUXi7Vr13ohhohku39+xfx0JM4/N+IG9prFeMIaSP0WnXmneWkC3Lq0KSIZpzMWInlNynnYMA1z4ws4zRTOWf685OnKL6X78nSXWpQI15LmIpJ5KhYieckfGzE+fgDXyX04gdVGLWb4D2Zwt5aMrh6Fw6E5KUTk+qhYiOQFZ09grXwcxw9v4uLCCqTjU/pSsE4P3oytRGh+P7sTikguoWIhkptZFvz0IcbnY3D99+bMRZ6beb/QQMZ1b0Q9zUkhIl6mYiGSW538E/OzkTj3rsIF7DFv4HFzMDfe1IH3mpXB3605KUTE+1QsRHIbwwPfvorx1URcnnMkWW5merrwY8m+TOpWh5jCBexOKCK5mIqFSG5yaDvGx8NwHd2BC/jWrMgU973c2bkVI2vdoJszRSTLqViI5AbJibBmEtbmV3BZJvFWfp729MGq0YfX21chrIC/3QlFJI9QsRDJ6fatuTCENOEADmCZ0Yh5QXczunszGpctbHc6EcljVCxEcqpzp7C+fDR1COlfVmEe9/Sn3I3debtVeQL9NXOmiGQ/FQuRnOjX5RjLRuBKPALAAk9rloYP4omeDalevKC92UQkT1OxEMlJEo9hLR+D46cPcAG/m5E8at5Do5YdWdxcQ0hFxH4qFiI5wcWJrj4bjev8CQzLwVyjA2ujBvBUj3qUKxpsd0IREUDFQsT3JRzG/PRBnL8txwX8YkbzmDWEW2Lbs6hRKVxODSEVEd+hYiHiqywLfngLY8UjuJITSLZcvOzpwo8xA3i+W22iw7QKqYj4HhULEV908k/MZcNx7l+DC/jRLM2Tzvu4rVssI2proisR8V0qFiK+xDThu9cxv3wcp+cs5y0/nvf04K8KA5jVtToRwfnsTigiclUqFiK+4tQBzKVDcf6xHicXpuN+2nUvg3q1YWz1KJ2lEJEcQcVCxG6WBd8vvHAvRcoZzln+POO5lSMV7uT1rjUoEhxgd0IRkXRTsRCxU8IhzGXDUpc2/84szxPO+xnUsxWdahTTWQoRyXFULETsYFmwYzHGZw/hSo4nyfLjOU9P/ijXj9e71SAiRPdSiEjOpGIhkt3OxF0Y8fHb57iA7WZpJjjvp2+PtjxSUyM+RCRnU7EQyU4/fYTx6Uhc50+SbLmY7unOnrIDmd29JkV1lkJEcgEVC5HskHgc87P/4Px5CS5gl1mSx533c1u3WxiteSlEJBdRsRDJar98hvHxA7jOHcNjOXnZ6MKO0oN4uXsdIkN1lkJEchcVC5Gscj4Ba8XDOLYvwgX8ahbnUYbSvVMHXqsXrbMUIpIrqViIZIU/v8H48G5cCQcwLQdzjA6sLTaIab3rUyJca3yISO6lYiHiTZ5kWDsZ6+vpuCyTg2YRRpv30bxVZxY1K62VSEUk11OxEPGWuF8wPhyM6+gOHMD7nma8Ez6Uib0bU7lYiN3pRESyhYqFyPUyTdgyB3Pl47iMJE5YQTziGUTJG2/lnTblCXC77E4oIpJtVCxErkfCIcwl9+LcvxYnsNaowfP5H2Bc75Y0KB1udzoRkWynYiGSWT99hPHJCFxJ8Zyz/JnkuZ2kmv1Z1LEKwfn87E4nImILFQuRjDqfgPX5KBw7FuMCdpgxTHCPYEivdrSpEml3OhERW6lYiGTEwa0YHwzEFf8nhuXgZaMz35e6m1d71dHCYSIiqFiIpI9pwMYXsNZMwmUZ/GUVZpQxjFZtO/HGjTE4NYxURARQsRC5toRDmB/ejfPPDTiAZUYjXi84jEm3NaFKsVC704mI+BQVC5Gr+eUzjCX34Uo6RaIVwHhPPwLr3sm77SsT6K9hpCIi/6ZiIXI5KeewvhiH47vXU2/QfNz9IEN7t6N15aJ2pxMR8VkqFiL/dnQXxvsDcB37BYBXPR3YXOpe5vSqqxs0RUSuQcVC5CLLgi1zMb8Yh8tMJs4qyEPGfTRp21M3aIqIpJOKhQjA2RNYS+/D8dtynMBqoxYvB4/gqTta6gZNEZEMULEQObgF471+uE7/TZLlxyTP7Zyu3p83u1SjQID+iIiIZIT+1pS8yzThm5cwVz2Jy/LwuxnJf6wH6dOtAz3qFLc7nYhIjqRiIXlT4nHMpUNw7vkSJ/Cx0Zj5YcOZ2qcJZSOC7E4nIpJjOTPy5FmzZlG9enVCQkIICQmhUaNGLF++PKuyiWSNA5vxzGqCc8+XJFl+PJwyiC21nuGd+1urVIiIXKcMnbEoXrw4U6ZMoWzZsgAsWLCAzp0788MPP1ClSpUsCSjiNaYJm2Zgrn4Kt2Wwz4ziIceDDOjdifbVo+xOJyKSKzgsy7Ku5w3CwsKYOnUqAwcOTNfzExISCA0NJT4+npCQkOvZtEj6JR7H/OhunPtWAbDEuJHFRUbwbJ8mlAjPb3M4ERHfl97P70zfY2EYBu+//z6JiYk0atTois9LSkoiKSkpTTCRbPXnpgsTXp05zHnLj8c9/Qhp2J+FsZXwd2foaqCIiFxDhovFzp07adSoEefPnycoKIglS5ZQuXLlKz5/8uTJPPHEE9cVUiRTTBO+no711URclsFesxhjnCO5t08nWmlabhGRLJHhSyHJyckcOHCAU6dO8eGHH/Laa6+xbt26K5aLy52xiI6O1qUQyVrnTmEtGYLjtws3F39kNOGdIiN4/o4biQ7TpQ8RkYxK76WQ677HolWrVpQpU4bZs2d7NZhIph35CePdPrhO/UGS5cd4T19cdfryWMcq5PPTiqQiIpmR5fdYXGRZVpozEiK22v4O5icjcBnn+csqzHDzP/Tp3olutTXhlYhIdshQsXjkkUeIjY0lOjqa06dP8+6777J27VpWrFiRVflE0seThLX8YRzb3sAJrDVqMD10FM/c0ZIKkcF2pxMRyTMyVCyOHj3KnXfeyeHDhwkNDaV69eqsWLGC1q1bZ1U+kWs7dRBj8V24Dn+PaTl40ejK3kpDeatHTYK01oeISLbK0N+6r7/+elblEMmcvavxfDAQ9/mTnLIKMNJzP01vuY2XGpfC4dAy5yIi2U3/nJOcyTRhwzSsNU/jxmKnWYrH8z3MowPaUadkIbvTiYjkWSoWkvOcO4n50T0493yBA3jH05IvS/6H126rT3hQgN3pRETyNBULyVnidmO8fWvqUNJHPf0JbzKQ19pWwOXUpQ8REbupWEjOsfsTjA/vxuU5e2EoqfUfBvTuqgXERER8iIqF+D7ThHVTYN0zuIBNRmWmBD/M1Ltu0lBSEREfo2Ihvu18AuaHg3HuuTBXyhuedmwqPZw3b61HaH4/m8OJiMi/qViI7zq2F8/bt+I+sYcky49HUgZyQ4sBzGlVHqfupxAR8UkqFuKbfvsS44OBuJMTOGyF8SCj6N+nO22rRNqdTERErkLFQnyLZcHG57FWP4ULi61meaYEP8IzfVtTNiLI7nQiInINKhbiO1LOYS29D8euj3AAizw3szrmP8zrU5+QfLqfQkQkJ1CxEN+QcBjjndtwHf6BFMvFeE8/CjQexNzYSpqfQkQkB1GxEPsd2n5h0qszhzlhBTHMGEnXbr3pUUdLnYuI5DQqFmKvnz++MOmVcZ495g2M8hvL4wPaU6dkmN3JREQkE1QsxB6WBRueg68m4gLWGdV5KXwcr/Rrzg0FA+1OJyIimaRiIdkv5TzmsmE4d74HXJj0aluF/7Cwd23y++uQFBHJyfS3uGSvM3EYb9+G69B3qTdpFmkxhJduLqdJr0REcgEVC8k+R37CWNQL1+m/OWUVYLj5ID1796FD9WJ2JxMRES9RsZDssXcVxuK7cKUkss+MYrTfWB7t24naJQrZnUxERLxIxUKy3rYFmJ8+iMsy2GRU5tmCj/JS/5ZEh+W3O5mIiHiZioVkHcuCr56CDdNwAh8aTfi4xMMsuKMhoYGaSVNEJDdSsZCs4UnCXDoU50/vAzDD041DNUbwerdq+LmcNocTEZGsomIh3nf2BMY7fXAd3ESK5eIRz0BiWt/DlOZlcDg08kNEJDdTsRDvOvkHnjd74D6xhwQrkAfMkfTsfSftq0fZnUxERLKBioV4z9/b8LzVC/e5YxyywnjQNY4xA7tr5IeISB6iYiHe8cvnmO/3x22cZ5dZkvFB43l+YCwlwjXyQ0QkL1GxkOu3bQHWpyNwWiZrjRq8GvEYs/s3IzwowO5kIiKSzVQsJPP+ZyExB/CupwUrSz/MG3fU05ofIiJ5lP72l8wxTawVY3BsmQPATE9nDtQYyexu1XFrOKmISJ6lYiEZ50nC/OgenD8vAWBCyl2ENL+fZ1qX13BSEZE8TsVCMibp9IU5Kv5YR7LlYrTnXup3ups+DUranUxERHyAioWk35l/8LzVA/eR7SRaAQwzR3Lr7f1oUyXS7mQiIuIjVCwkfU7+gWdBV9ynfue4FcwDzkcYOeBW6pQMszuZiIj4EBULubYjP+FZ2A332aP8ZRVmpP94nh7UlXJFg+1OJiIiPkbFQq7u4BaMN3vgTo5ntxnNowWeYMbdsRQvpImvRETkUioWcmW/r8N4+1ZcnrNsNcvzTKEneXXQTRQJ1sRXIiJyeSoWcnm/rsBcfCcuM5n1RjVmRz3J6/2bEhroZ3cyERHxYSoWcqmfPsT88G6clocvjLq8V+oJXruzEYH+LruTiYiIj1OxkLS+X4i17AGcWCwxbmRNxQnM6l0Xf7dm0xQRkWtTsZD/t3kWrHgYB7DIczO7ao3nha7VcTk1m6aIiKSPioVcWExs/XOwZiIAczztiWswjqc7VNYU3SIikiEqFnmdZcGap2H9VACeT+mB0XQU49pWVKkQEZEMU7HIyywLa/WTODY+D8DElD6E3vQgI28uZ3MwERHJqVQs8irLwlo5HsemGQA8kXInkW0e5J7mZWwOJiIiOZmKRV5kWVhfPIpj80wAHk/pS+lbHqTfjTE2BxMRkZwuQ2MIJ0+eTL169QgODiYiIoIuXbrw66+/ZlU2yQqWhbXi4dRS8WhKfyp2GqVSISIiXpGhYrFu3TqGDh3K5s2bWblyJR6PhzZt2pCYmJhV+cSbLAtr+UM4vn0VgEdSBlKz2yhub1DC5mAiIpJbOCzLsjL74n/++YeIiAjWrVtHs2bN0vWahIQEQkNDiY+PJyQkJLOblowyTazPR+P47jVMy8FYzyAa9hhB11rF7U4mIiI5QHo/v6/rHov4+HgAwsLCrvicpKQkkpKS0gSTbHbxTMV/S8VDnrtp1O0BlQoREfG6TM/TbFkWI0eOpEmTJlStWvWKz5s8eTKhoaGpX9HR0ZndpGSGZWF9MQ7H1rmYloPRnnto0HUY3euoVIiIiPdl+lLI0KFD+eyzz9i4cSPFi1/5Q+pyZyyio6N1KSSbWKufwrHhOQDGpAymTpfh9KqnciciIhmTpZdChg0bxrJly1i/fv1VSwVAQEAAAQEBmdmMXCdr3bOppeKxlH7U7PyASoWIiGSpDBULy7IYNmwYS5YsYe3atcTEaIiir7K+fhHHmqeBCzNqVuw0ktvqa/SHiIhkrQwVi6FDh/L222/z8ccfExwczJEjRwAIDQ0lMDAwSwJKJnw7B8fKxwCYmtKLkh3H0KdBSZtDiYhIXpCheyyutCjVvHnz6NevX7reQ8NNs9i2+fDJcABe9HQhqN0EBjTRmSUREbk+WXKPxXVMeSHZ4cfFWJ+MwAHM9rTH0WKcSoWIiGSrTA83FR/z63LMpffiwGK+pw0nGj3K/VqlVEREspmKRW7w5yaM9/ritAw+MJqxt/ajPHxLpSteuhIREckqWt00pzu8A89bPXEbSaw06rCp8uM816W6SoWIiNhCxSInO76PlIVd8Us5w7dmRZaUeYoZvergdKpUiIiIPVQscqqEwyTP74z/uWPsMkvyWvHJvNSnIX4uXd0SERH7qFjkROdOkrKgC/6nD7LfLMpzEZN4uV8z8vm57E4mIiJ5nIpFTpOciOetnvgd/4UjViEeC5nIzAFtyO+v/5UiImI/fRrlJEYKxuK7cP+9lVNWAUb6Pc5zgzpSML+/3clEREQADTfNOSwL85MRuPat4qwVwP2OsYwf1JNiBTWVuoiI+A4VixzCWvcMzu1vYVgORpgPMOyu26kQGWx3LBERkTRULHKCHxbhWDsZgMc8A+jWeyANSofbHEpERORSKha+bu9qzGUPAPCypxOVOg6nXdUom0OJiIhcnoqFLzu8A8+7d+K0PCwxbuR803Hc2VDLn4uIiO9SsfBV8X+R8mZ33J5ENhmV2VLtSUa2qWB3KhERkavScFNflHSGlLd64Xc2jl/N4iwoMZGZ3Wtr/Q8REfF5Kha+xjTwfDAQv3928Y8VwlOhE3jlzmaaqltERHIEfVr5GHPleNx7VpBk+THa/TCT+7cnJJ+f3bFERETSRcXCl3y/EOc3LwEw1hzC8L63Ex2W3+ZQIiIi6adi4Sv2b8D85EEApnu60arXUGqVKGRzKBERkYxRsfAFx/eR8s4dOC0PnxgN8b/5EW6pprkqREQk51GxsNv5BJLf6oVf8im2m2X4uupT3NuirN2pREREMkWjQuxkmqR8MBj/k3s5bIUxo8gTvNq9joaViohIjqUzFjYy1z2D394LI0Ae8RvDM31bE+B22R1LREQk01Qs7PLL5zjXTQHgcXMgI/reSkRIPptDiYiIXB8VCzv88xspHwwCYL6nDY26PUCN6IL2ZhIREfECFYvsdj6epLd64+dJ5FuzIkcbPk6XWjfYnUpERMQrdPNmdjJNUt4fRED87/xthfNm9JPMuKWq3alERES8RsUiG5kbpuG370uSLD/G53uYaX1uwuXUCBAREck9dCkku+xfD2smATDeHMiIu3oTml9rgIiISO6iYpEdTh8heXF/nJi852lO7U73U/WGULtTiYiIeJ2KRVYzPCS92w//88fYbUazs8aj9KoXbXcqERGRLKF7LLKY56unCfj7G85Y+ZgeNo4ZXerYHUlERCTL6IxFVvrtS9xfPw/Ak44hPHpXZ/L5aWZNERHJvVQsssqpgyT/dxKsBZ7WxN46lOiw/DaHEhERyVoqFlnB8JC0uD/+yfH8aJbmaKPHaFkxwu5UIiIiWU73WGQBY91UAg5vJcEK5OXCj/Jyu2p2RxIREckWOmPhbQc241j/LABPM5jH7ojFz6XdLCIieYM+8bzpfDznFw/AiclHRhOadLtX91WIiEieomLhLZZF0tIR5Ev8mwNmEX6s9igdaxSzO5WIiEi20j0WXmL++C4Bv3yEx3IyNWg0z3apb3ckERGRbKczFt5wYj+eT/4DwEyzO0Pvuo1Af81XISIieY+KxfUyPJxfPAB/I5EtZgXC2o2lYmSI3alERERsoWJxnYyN08l39HsSrPy8fcOj3NGotN2RREREbJPhYrF+/Xo6duxIsWLFcDgcLF26NAti5RBHf4a1UwB4xtGPMbe2xul02BxKRETEPhkuFomJidSoUYOZM2dmRZ6cw/Bw9v17cFkprDZqUa/TUKJCA+1OJSIiYqsMjwqJjY0lNjY2K7LkKCkbXiD/sR3EW/lZVXYsk2rdYHckERER2+kei8w4+jPOdRcugTzvGsjoHi1xOHQJREREJMvnsUhKSiIpKSn1+4SEhKzeZNYyUjjz3mCCLA+rjFo063U/YQX87U4lIiLiE7L8jMXkyZMJDQ1N/YqOjs7qTWaplPUvEHT8J+Kt/Gyu/Dg3V460O5KIiIjPyPJiMXbsWOLj41O/Dh48mNWbzDpHf8a5/hkApvsNYnjXpjYHEhER8S1ZfikkICCAgICArN5M1jMNEj+4jwKWh5VGbZr2HkpwPj+7U4mIiPiUDBeLM2fOsHfv3tTv9+/fz/bt2wkLC6NEiRJeDedLPFtep8A/P3DaCmRD+bE8WUmXQERERP4tw8Xiu+++o2XLlqnfjxw5EoC+ffsyf/58rwXzKQmHMVZOwA287LyN4V2b251IRETEJ2W4WLRo0QLLsrIii886/fEogo1EtptlqNT1QcKDcsGlHRERkSygeSyuwfp1BcH7PsVjOfnohtF0qpmzR7WIiIhkJRWLq0lO5NzSEQDMt9ozuGdnTYQlIiJyFSoWV3F+zXPkP3eYv6zCGM3GEB2W3+5IIiIiPk3F4kpO/I5780sAzAkcTP8WVWwOJCIi4vtULK4gfulo3FYK641qtOsxEH+3dpWIiMi16NPyMsxfvyT0wCpSLBdflxtN47JF7I4kIiKSI6hY/JsnicRlowB4y3ELg7q2szmQiIhIzqFi8S9JG2cSnPgncVZBnM0fokiw5qwQERFJLxWL/3XmH1j/HACv5+vHbU2r2hxIREQkZ1Gx+B8JXzxFgHmWH83SNOgyRDdsioiIZJA+OS/65zcK7HwTgM8i76NlRS0yJiIiklEqFv91/OOxuDBZZdamV49bNcOmiIhIJqhYAOb+jYT/tQqP5eSXKqMoGxFsdyQREZEcScXCNIn/+GEAPuRm+nRobXMgERGRnCvPF4uUnz6i0KmdnLHykdh4NIUK+NsdSUREJMfK28XCNEj8YiIAb7s6c2vLOjYHEhERydnydLE49/27FEzcz0kriII3jSC/v9vuSCIiIjla3i0WRgpJqyYB8EFAV7o1qmhzIBERkZwvzxaL09++ScHzf3HMCqFk7Ajcrjy7K0RERLwmb36aepIx1j4DwMdBvWhds4zNgURERHKHPFkszmx+g4LJR4izClI2drgmwxIREfGSvFcsPEmY6y4sNLYkqDfNqkTbHEhERCT3yHPF4szWRYSk/MMRqxDlY+/X2QoREREvylvFwjRIWvcCAJ8V6EYLna0QERHxqjxVLE7/+DHh5w8Qb+Unpu1Qna0QERHxsrxTLCyL06su3FuxPLAjLauXtjmQiIhI7pNnisXZ39ZQLHEX5y0/irbWSBAREZGskGeKxfEvpgLwhX8rmteqbHMaERGR3ClPFIuUv38k+sQmDMuB48ZhOJ06WyEiIpIV8kSxOLL8WQBWOm+kzY0NbU4jIiKSe+X6YmGdiSPyrxUAnKxxD/n8XDYnEhERyb1yfbE4uOpV/PCw3SpL29bt7I4jIiKSq+XuYmF4KLBzIQB7SvQmrIC/zYFERERyt1xdLA5vXUK48Q/HrWDq3jLA7jgiIiK5Xq4uFolfzwZgS8H2xEQVtjmNiIhI7pdri8W5w79Q9vRWTMtBeIt77I4jIiKSJ+TaYnHwi5cA+MZVh7o1atmcRkREJG/IncUiOZEb/lwCwMkqfTUhloiISDbJlcXi76/foYCVyAErgoZtetodR0REJM/IlcUi+bs3AdheuCOFgwNtTiMiIpJ35LpicfbIHmISt2NaDqKa9bc7joiISJ6S64rF/lVzAfjOXYM61aranEZERCRvyV3FwjQp8vtHACRU6KWbNkVERLJZrioWB7avJsL8h9NWILVa97E7joiISJ6TqWLxyiuvEBMTQ758+ahTpw4bNmzwdq5MObb5bQB2hjQjvFBBe8OIiIjkQRkuFosXL2bEiBGMGzeOH374gaZNmxIbG8uBAweyIl+6mSnJxMStAsCvhoaYioiI2CHDxeL5559n4MCBDBo0iEqVKjF9+nSio6OZNWtWVuRLt1+++YRCJHCCEKo17WRrFhERkbwqQ8UiOTmZbdu20aZNmzSPt2nThk2bNl32NUlJSSQkJKT5ygpnv38PgF/DbyZfQECWbENERESuLkPF4tixYxiGQdGiRdM8XrRoUY4cOXLZ10yePJnQ0NDUr+jo6MynvYLz5xKpcHIdACF1b/X6+4uIiEj6ZOrmTYcj7TBOy7IueeyisWPHEh8fn/p18ODBzGzymnl+qzueraFtqFS/tdffX0RERNLHnZEnFy5cGJfLdcnZibi4uEvOYlwUEBBAQBZfmgjIl586He8F7s3S7YiIiMjVZeiMhb+/P3Xq1GHlypVpHl+5ciWNGzf2ajARERHJeTJ0xgJg5MiR3HnnndStW5dGjRoxZ84cDhw4wJAhQ7Iin4iIiOQgGS4WvXv35vjx4zz55JMcPnyYqlWr8vnnn1OyZMmsyCciIiI5iMOyLCs7N5iQkEBoaCjx8fGEhIRk56ZFREQkk9L7+Z2r1goRERERe6lYiIiIiNeoWIiIiIjXqFiIiIiI16hYiIiIiNeoWIiIiIjXqFiIiIiI16hYiIiIiNeoWIiIiIjXZHhK7+t1caLPhISE7N60iIiIZNLFz+1rTdid7cXi9OnTAERHR2f3pkVEROQ6nT59mtDQ0Cv+PNvXCjFNk0OHDhEcHIzD4fDa+yYkJBAdHc3Bgwe1BkkW077OHtrP2Uf7OntoP2efrNjXlmVx+vRpihUrhtN55Tspsv2MhdPppHjx4ln2/iEhITpgs4n2dfbQfs4+2tfZQ/s5+3h7X1/tTMVFunlTREREvEbFQkRERLwm1xSLgIAAxo8fT0BAgN1Rcj3t6+yh/Zx9tK+zh/Zz9rFzX2f7zZsiIiKSe+WaMxYiIiJiPxULERER8RoVCxEREfEaFQsRERHxmlxTLF555RViYmLIly8fderUYcOGDXZHylUmTJiAw+FI8xUZGWl3rFxh/fr1dOzYkWLFiuFwOFi6dGman1uWxYQJEyhWrBiBgYG0aNGCXbt22RM2B7vWfu7Xr98lx3jDhg3tCZuDTZ48mXr16hEcHExERARdunTh119/TfMcHdPekZ59bcdxnSuKxeLFixkxYgTjxo3jhx9+oGnTpsTGxnLgwAG7o+UqVapU4fDhw6lfO3futDtSrpCYmEiNGjWYOXPmZX/+7LPP8vzzzzNz5ky2bt1KZGQkrVu3Tl13R9LnWvsZoF27dmmO8c8//zwbE+YO69atY+jQoWzevJmVK1fi8Xho06YNiYmJqc/RMe0d6dnXYMNxbeUC9evXt4YMGZLmsYoVK1oPP/ywTYlyn/Hjx1s1atSwO0auB1hLlixJ/d40TSsyMtKaMmVK6mPnz5+3QkNDrVdffdWGhLnDv/ezZVlW3759rc6dO9uSJzeLi4uzAGvdunWWZemYzkr/3teWZc9xnePPWCQnJ7Nt2zbatGmT5vE2bdqwadMmm1LlTnv27KFYsWLExMRw66238vvvv9sdKdfbv38/R44cSXN8BwQE0Lx5cx3fWWDt2rVERERQvnx5Bg8eTFxcnN2Rcrz4+HgAwsLCAB3TWenf+/qi7D6uc3yxOHbsGIZhULRo0TSPFy1alCNHjtiUKvdp0KABCxcu5IsvvmDu3LkcOXKExo0bc/z4cbuj5WoXj2Ed31kvNjaWRYsW8dVXXzFt2jS2bt3KTTfdRFJSkt3RcizLshg5ciRNmjShatWqgI7prHK5fQ32HNfZvrppVvn3EuyWZXl1Wfa8LjY2NvW/q1WrRqNGjShTpgwLFixg5MiRNibLG3R8Z73evXun/nfVqlWpW7cuJUuW5LPPPqNbt242Jsu57r//fnbs2MHGjRsv+ZmOae+60r6247jO8WcsChcujMvluqTpxsXFXdKIxXsKFChAtWrV2LNnj91RcrWLI290fGe/qKgoSpYsqWM8k4YNG8ayZctYs2YNxYsXT31cx7T3XWlfX052HNc5vlj4+/tTp04dVq5cmebxlStX0rhxY5tS5X5JSUns3r2bqKgou6PkajExMURGRqY5vpOTk1m3bp2O7yx2/PhxDh48qGM8gyzL4v777+ejjz7iq6++IiYmJs3PdUx7z7X29eVkx3GdKy6FjBw5kjvvvJO6devSqFEj5syZw4EDBxgyZIjd0XKNUaNG0bFjR0qUKEFcXBwTJ04kISGBvn372h0txztz5gx79+5N/X7//v1s376dsLAwSpQowYgRI5g0aRLlypWjXLlyTJo0ifz583P77bfbmDrnudp+DgsLY8KECXTv3p2oqCj++OMPHnnkEQoXLkzXrl1tTJ3zDB06lLfffpuPP/6Y4ODg1DMToaGhBAYG4nA4dEx7ybX29ZkzZ+w5rrN1DEoWevnll62SJUta/v7+Vu3atdMMt5Hr17t3bysqKsry8/OzihUrZnXr1s3atWuX3bFyhTVr1ljAJV99+/a1LOvC8Lzx48dbkZGRVkBAgNWsWTNr586d9obOga62n8+ePWu1adPGKlKkiOXn52eVKFHC6tu3r3XgwAG7Y+c4l9vHgDVv3rzU5+iY9o5r7Wu7jmstmy4iIiJek+PvsRARERHfoWIhIiIiXqNiISIiIl6jYiEiIiJeo2IhIiIiXqNiISIiIl6jYiEiIiJeo2IhIiIiXqNiISIiIl6jYiEiIiJeo2IhIiIiXqNiISIiIl7zf2+qatpbr71DAAAAAElFTkSuQmCC",
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAj4AAAGdCAYAAAASUnlxAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABQ10lEQVR4nO3df1xUdaI//tcwzA9+DvJzIGEABU2xUjDEQqzdUNxvjyw/xW7fj+ntXu9yqy3k483S7ZPb7Yq2bXX7+mu7a97cdlt313LdtE3clCyx1MUfKSoqCgIjDMgMP2dgON8/hjk6MSKDzDn8eD0fj3kIZ95zzvucxt4v3+/3eR+FIAgCiIiIiEYBH7krQERERCQVBh8iIiIaNRh8iIiIaNRg8CEiIqJRg8GHiIiIRg0GHyIiIho1GHyIiIho1GDwISIiolHDV+4KDCXd3d2oqalBUFAQFAqF3NUhIiKifhAEAc3NzYiJiYGPT999Ogw+N6ipqUFsbKzc1SAiIqIBqKqqwtixY/ssw+Bzg6CgIACOCxccHCxzbYiIiKg/LBYLYmNjxXa8Lww+N3AObwUHBzP4EBERDTP9mabCyc1EREQ0ajD4EBER0ajB4ENERESjBuf4EBER3QZBENDV1QW73S53VUY0pVIJX1/f215uhsGHiIhogGw2G2pra9HW1iZ3VUYFf39/REdHQ61WD3gfDD5EREQD0N3djYqKCiiVSsTExECtVnPxWy8RBAE2mw319fWoqKhAUlLSLRcqvBkGHyIiogGw2Wzo7u5GbGws/P395a7OiOfn5weVSoXLly/DZrNBq9UOaD+c3ExERHQbBtrzQJ4bjGvN/1pEREQ0ajD4EBER0ajB4ENERESjBoMPERHRKLN48WIoFIper7lz58pdNa/jXV10W7rs3ahuasdVixV1zR24arGi1doFW1c3Ou3d6LQLUPkqoPVVQqtSIlCjRGiABqEBaoQHqqHXaRGkVcl9GkREo87cuXOxZcsWl20ajcZt2c7OTqhUqltuGyw2m+221urpC4MPeaTW3I79Z+tx4ooZp2vMOGNshrWr+7b2GRqgRmyoP+JC/REX6gdDWAAmRAUhOSoIfmrlINWciMi7BEFAe6c8qzf7qZQeryGk0Wig1+vdvqdQKLBx40Z89tln2Lt3L5YtWwaFQoEdO3bg+eefx+uvv45Lly7BbrfDYrHg3//937Fjxw50dHQgLS0Nb7/9Nu6++24Ajt6lpqYm7NixQ9x/fn4+jh07hv379wMAZs+ejZSUFKjVamzduhWTJ09GcXHxgK7FrTD40C2VX23GrpO12Ft2Fd9VW3q9r1X5ICpYi6ggLSKCNQjWqqDx9YFKqYCv0gedXd3o6LKjo7MbLR1daGy1oaHVClOLDeb2TjS22tDYasPxqiaX/SoUQHxYACbqgzBBH4S7xupwT+wYhAZ4518BRES3o73Tjkn/93NZjn36tTnwVw9uk/7qq6+isLAQb7/9NpRKJbZs2YLz58/jj3/8I7Zv3w6l0vEP0x/96EcIDQ3F7t27odPp8Otf/xo/+MEPcO7cOYSGhvb7eB988AH+7d/+DV9//TUEQRjUc7kRgw+51d0tYM9pI97/6hK+vdQoblcogKmxIUhPDMPkmGBMjtHBEOoPH5+BrVbaYu1CZUMbKhvbUNXo+POiqQVnapvR0GpDhakVFaZWfPadUfxMXKg/psaF4J7YEEyNG4M7o4Og8WXPEBGRJz799FMEBga6bFu+fDleeeUVAMCTTz6Jp59+2uV9m82G3/72t4iIiAAAfPHFFzh58iTq6urEYbI333wTO3bswJ///Gf867/+a7/rM378eLzxxhu3c0r9wuBDLgRBwN6yOvzy8zM4d7UFAKD0UeCBCRHInqzHgxMjER7ofgx4IAI1vpgUE4xJMcG93qtvtuKssRlnjBacrrXgxBUzzte1oLInIP3lWA0AQK30wZSxOsxIDMWMxDCkGsYM+r98iIhuxU+lxOnX5sh2bE898MAD2Lhxo8u2G3to0tLSen3GYDCIoQcAjh49ipaWFoSFhbmUa29vx4ULFzyqj7vjeQNbBxJVNrRh5Y6TOFBuAgAEaXyxaGY8/vcMA/S6gS0NfjsigjSICNLg/qRwcZu5vRMnrjShtLIJx6ocr8ZWG45evoajl69h/b4L8PVR4K6xOsxIDBODUICGX3Ui8i6FQjGs/tEVEBCA8ePH9/n+rbZ1d3cjOjpanKtzo5CQEACO1Za/P3TV2dnZr+N5w4BuZ9+wYQMSEhKg1WqRmpqKAwcO9Fm+uLgYqamp0Gq1SExMxKZNm1zeP3XqFBYsWID4+HgoFAq88847vfbhfO/7r2effVYs4+72vBkzZgzkFEcVQRDw+28qMfe/vsSBchPUSh/82+xx+OqlB7FszgRZQs/N6PxUyEyKwPM/SML7i6fj6M9/iOJ/n41f/q+7sGDaWNwR4oeubgH/qGzChv0X8NT73+LuX+zBoxu+xi8/P4NDFxtgu83J2ERE5DBt2jQYjUb4+vpi/PjxLq/wcMc/WiMiIlBbW+vyuWPHjslQWwePo+m2bduQn5+PDRs24L777sOvf/1r5OTk4PTp04iLi+tVvqKiAvPmzcOSJUvw4Ycf4uuvv8YzzzyDiIgILFiwAADQ1taGxMREPP7441i6dKnb4x4+fBh2+/XZ8t999x0eeughPP744y7lvn97nrduhxsp2m12vPzxCezoGTZKTwjFG//rLhjCpEnet0uhUMAQFgBDWAAeT4sFAFQ1tuHQxQYcutiIbyoacOVaO0orHb1E6/ddgL9aiYzEMGQmhSMzOQKJ4QF8ojIRjTpWqxVGo9Flm6+vrxhY+uOHP/whMjIyMH/+fKxduxYTJkxATU0Ndu/ejfnz5yMtLQ0PPvggfvnLX2Lr1q3IyMjAhx9+iO+++w5Tp04d7FPqF4+Dz1tvvYV//ud/xr/8y78AAN555x18/vnn2LhxIwoLC3uV37RpE+Li4sRenDvvvBNHjhzBm2++KQaf6dOnY/r06QCAl156ye1xbxxTBIA1a9Zg3LhxyMrKctne1+155MrUYsU/bTmMk9VmKH0UWD53Av7l/sQBT1QeKmJD/REb6t8rCH193oSvzptgarHh72fq8PczdQCAO0L8HCEoKQL3jQ9DiD/DMhGNfH/7298QHR3tsm3ChAk4c+ZMv/ehUCiwe/durFy5Ek8//TTq6+uh1+sxa9YsREVFAQDmzJmDV155BS+++CI6Ojrw9NNP46mnnsLJkycH9Xz6XWfBg3vGbDYb/P398ac//QmPPvqouP2FF17AsWPH3N5zP2vWLEydOhX/9V//JW775JNP8MQTT6Ctra3X4kfx8fHIz89Hfn5+n/WIiYlBQUEBVqxYIW5fvHgxduzYAbVajZCQEGRlZeE///M/ERkZ6XY/VqsVVqtV/N1isSA2NhZmsxnBwb0n244kV6614anN3+KiqRWhAWqsf3IaMsaF3fqDw1x3t4AyowUHyk04UF6PwxXXYLNfH/ryUQB3jQ3BrOQIPDAhAneNDYFymAdBIvKOjo4OVFRUiFM/yPtuds0tFgt0Ol2/2m+PenxMJhPsdruY4pyioqJ6dZc5GY1Gt+W7urpgMpl6pc3+2LFjB5qamrB48WKX7Tk5OXj88cdhMBhQUVGBV155BQ8++CCOHj3qdjXKwsJC/OIXv/D4+MOd0dyBH793CFeuteOOED98+C/pSAgfHkNbt8vHR4HJMTpMjtEhL2sc2m12fFPRIAahc1dbxEnT7/69HKEBamQlR2D2hAhkJUewN4iIaJgb0PTz78+HEAShzzkS7sq7295fmzdvRk5ODmJiYly25+bmij+npKQgLS0NBoMBu3btwmOPPdZrPy+//DIKCgrE3509PiNZU5sNT73/Da5ca4chzB9/+NcZiNb5yV0t2fiplZg9IRKzJzh6BY3mDnxZXo/is/X4srweja02fFJajU9Kq+GjAKbGjcEDEyIwe0IkJscEc24QEdEw41HwCQ8Ph1Kp7NW7U1dX16tXx0mv17st7+vr2+u+//64fPky9u7di48//viWZaOjo2EwGFBeXu72fY1Gc9PnkoxEXfZu5H14FOeutiAqWIMP/zl9VIced/Q6LZ5Ii8UTabHotHfjH5evYd/Zeuw/W4czxmbxtvk395xDZJAGsydE4IEJkbgvKRzBfOYYEdGQ51HwUavVSE1NRVFRkcscn6KiIjzyyCNuP5ORkYG//vWvLtv27NmDtLS0AT3cbMuWLYiMjMSPfvSjW5ZtaGhAVVXVgIbTRqLVu8/g0MVGBKiV2Pp0OmJD/eWu0pCmUvogPTEM6YlheClnImqaHM8p23e2Dl+fN6Gu2Yo/HrmCPx65Al8fBdLix+CBCZF4cGIkxkcGsjeIiGgI8nioq6CgAAsXLkRaWhoyMjLw3nvvobKyEnl5eQAcw0fV1dXYunUrACAvLw/r1q1DQUEBlixZgpKSEmzevBkfffSRuE+bzYbTp0+LP1dXV+PYsWMIDAx0WVypu7sbW7ZswaJFi+Dr61r1lpYWrFq1CgsWLEB0dDQuXbqEFStWIDw83CWkjVafnqjB+19XAAB+9cQ9mKAPkrlGw09MiB+eTI/Dk+lxsHbZcbjiGvadrcO+s3W4WN+KQxcbcehiIwo/O4M7QvzE3qCZ48OG1aJmROQZbz5XilwNxrX26K4upw0bNuCNN95AbW0tUlJS8Pbbb2PWrFkAHHdWXbp0yWUVx+LiYixduhSnTp1CTEwMli9fLgYlALh06RISEhJ6HScrK8tlP3v27MGcOXNw9uxZJCcnu5Rtb2/H/PnzUVpaiqamJkRHR+OBBx7Af/zHf/R73o4ns8KHk6uWDmS//SXM7Z14ZvY4vDh3otxVGnEuN7SKvUElFxpcnliv9vXBjMQwPNAThOJHyURyopHObrfj3LlziIyMHNDUDfJcQ0MD6urqkJycLD4kFfCs/R5Q8BmpRmLwEQQBi7Ycxpfn6jHlDh0+fmYmVMoBLdhN/dRus+PQxQZ8ccbRG3TlWrvL+wnhAWJv0L0JodAO4Bk7RDQ01NbWoqmpCZGRkfD39+cQt5cIgoC2tjbU1dUhJCSk1xQWBp8BGonB5+N/XEHBH49D4+uDXc/fj/GRHOKSkiAIuFDfgn1nHL1Bhy81otN+/a+cn0qJ+8aH4YGJjjvL7gjhZHOi4UQQBBiNRjQ1NcldlVEhJCQEer2+V8Bk8BmgkRZ8LB2dePDNYpharHhx7gQ8M/vmD6MjaTR3dOLr8w3Y3zM36KrF6vL+hKggzJ7o6A1KNYxh7xzRMGG3290+eJMGj0qlchneuhGDzwCNtODzH5+exuavKpAYHoDP8jOh8eWQylAiCALKapsdE6TP1OEfldfQfcPfxiCNLzKTw3vWGYpAZBBXhiUicofBZ4BGUvC5ZGrFD94qhr1bwAdP34us5Ihbf4hk1dRmw5flJuw/U4f95xyLJ95oyh06x+KJEyNxNx+lQUQkYvAZoJEUfAq2HcPHpdXISo7AB0/fK3d1yEPd3QJOVJux70wd9p+tw/ErZpf3x/irkJUcgQcmRmJWUgTGBPBRGkQ0ejH4DNBICT7n61qQ/XYxugVg53P34a6xIXJXiW5TfbMVxeccE6S/PFeP5o4u8T0fBXBPbAgemBCJBybyURpENPow+AzQSAk+P/uoFH89XoPsSVF476k0uatDg6zL3o1/VDaJc4POGJtd3uejNIhotGHwGaCREHyqGtuQ9ct96BaA3c9nYlLM8DwP6r9ac8+jNM7U4avzJrTZ7OJ7vj4KTIsbg1nJ4ZiVHIGUGB18ODeIiEYYBp8BGgnBx3knV2ZSOH77z+lyV4ck5u5RGjcKDVDj/vHhyExyBKGoYN4pRkTDH4PPAA334NPc0YmMwi/QYu3C//zTdMyeECl3lUhmlQ1t+LK8Hl+eq8fBCw1osXa5vD9RH4RZyRHITArH9HiuIk1EwxODzwAN9+DzmwMX8fquMoyPDETR0lmc4EouOu3dKK1swoGeIHSi2owb//ZrVT5ITwjDrOQIzEoK5xPmiWjY8KT95iOjRwhBEPD7bysBAP90XzwbLOpFpfTBvQmhuDchFP8newIaW2346rwJX56rx4Hyely1OO4cKz5XDwCI0WmRmRSBWckRuH98OHT+nCRNRMMfe3xuMJx7fI5evoYFGw/CT6XE4Z//EIEaZlrqP0EQcO5qC748V48vy+vxTUUjbDc8Yd5HAdwdG4LMJMew2D2xIXycBhENGezxGYX+fLQKAJAzRc/QQx5TKBSYoA/CBH0QlsxKRLvNjm8vNTqC0Ll6lNe1oLSyCaWVTXj37+XwVyuRnhCK+8aH477x4ZgQFcS7xYhoWGALOQK02+z46/FaAMDjqbEy14ZGAj+1ElnJEeKjTmqa2h1zg8pNKLnQgMZWG/adrce+s45hsbAANTLGheH+niAUG+ovZ/WJiG6KwWcE+PyUES3WLsSG+iE9IVTu6tAIFBPih9zpccidHofubgFnjM34+rwJX18w4duKRjS02vDpiVp8esIRwONC/XHf+DDcNz4cM8eFI5SP1CCiIYLBZwT49EQNAODRqWM53EBe5+OjwKSYYEyKCcaSWYmwdXXjWFWTIwidN+FYVRMqG9tQ+W0bPvrWMQQ7KTpYDEL3JoTCX83/9RCRPDi5+QbDcXJzc0cnUl/fC1tXNz7Pn4UJ+iC5q0SjXIu1C99WNODr8w34+ryp1yM1VEoFpsaOQca4MMxIDMPUuBCuH0REt4WTm0eRL87UwdbVjcTwACRHBcpdHSIEanzx4MQoPDgxCoDjAasHL5hw8HwDvjpvQnVTO7691IhvLzXiv/5eDrWvD6bFhWBGoiMI3RPLIERE3sPgM8z97TsjAMfdXFy7h4aiiCANHrnnDjxyzx0QBAGVjW34+nwDDl10vOqarTh0sRGHLjYCcA1CGYlhuCcuBBpfBiEiGhwMPsNYm60L+87WAQByUqJlrg3RrSkUChjCAmAIC8CT6XEQBAEVptae4NOAkosNqL8hCL2Dcmh8fTAtbkxPj1AogxAR3RYGn2HsQLkJHZ3dGDvGD5P5FHYahhQKBRIjApEYESgGoYum1p7eIEcYqm+2oqQnFAFwCUIZ48Jwd6yOQYiI+o3BZxhzPlrgBxMjOcxFI4JCocC4iECMiwjE/5tuEINQyYUGMQyZWq4Hobf3AmpfH9wzNgTTE8ZgenwoUg1jEKTl4zWIyD0Gn2FKEAR82RN8siZEyFwbIu+4MQj97xmOIHShvlWcH+QMQs7J0sAF+CiAO6ODMT3e8Vyy6fGhiAjSyH0qRDREMPgMUxWmVly51g610gczEsPkrg6RJBQKBcZHBmJ85PUgVGFqxeFLjfi24hoOX2pEZWMbTtVYcKrGgv85eAkAkBAegOnxY8QwFBfqz15SolGKwWeYcvb2pMWP4WJwNGrdOEcod3ocAOCqpQPfVjT2hKFGnL3ajApTKypMrfjjkSsAgMggDaYnhOLeeEeP0AR9EJRc/JNoVGCLOUx9WW4CAMxK5jAX0Y2igrV4+O4YPHx3DADA3N6Jo5ev9widuNKEumYrdp2oxa6eR2wEaX2RahiD1LgxSDWMwd2xIQjgw36JRiT+zR6GrF12lFxw3OGSxeBD1Cedn8plQcWOTjuOVTXhcIVjXtA/Ll9Dc0cX9p+tx/6eh676KICJ+mBHGOp5jR3jx+ExohGAwWcYOnnFjPZOO8ID1ZjIR1QQeUSrUoqrRANAl70bZbXN+EflNRy97HhVN7XjdK0Fp2st+O2hywAcCzE6e4SmGcYg5Y5g3kZPNAwx+AxDjrtXgOnxofwXKNFt8lX6YMpYHaaM1WHRzHgAgNHc4RKETtWYUd9sxd9OGfG3U47V0tU9n0s1jMG0uBBMM4xBZJBWxjMhov5g8BmGDlc4gk9afKjMNSEamfQ6LeZNica8KY4V0Ts67ThZbRaD0D8uX0NDq0383Sk21A/T4sbg7rEhuCcuBJNj2CtENNQw+Awz3d0CjvT8j/ZeBh8iSWhVSkzvuQMMcKyjdbmhzaVX6OzVZlQ1tqOqsR1/OVYDwPEk+knRwbgn1hGE7okdg/gw3kpPJCcGn2Hm7NVmNHd0IUCtxJ3RnN9DJAeFQoH48ADEhwfgsWljAQCWjk4cq2zCsarrr8ZWG45fMeP4FTM+KHHMFdL5qXB3bAjuiQ3B1NgQ3B0bgtAAtZynQzSqMPgMM0d65vdMM4yBr9JH5toQkVOwVoVZyRHiEhOCIKCqsR2lVddwrKoJx6ua8F2NBeb2Tnx5rl5ciwsADGH+uCc2hENkRBJg8Blmvr3kGOZKM3CYi2goUygUiAvzR1yYPx655w4AgK2rG2eMFkePUE/v0EVTKy43tOFyQ5vbIbK7Y0Nw11gdEsIDucgi0SBg8Blmjop3dI2RuSZE5Cm1rw/uGhuCu8aG4KkMx7amNsdwmCMIOXqHrrV1ikNk6BkiC1ArMfkOHe66w3EH2l1jQ2AI9YcPwxCRRwY0VrJhwwYkJCRAq9UiNTUVBw4c6LN8cXExUlNTodVqkZiYiE2bNrm8f+rUKSxYsADx8fFQKBR45513eu1j1apVUCgULi+9Xu9SRhAErFq1CjExMfDz88Ps2bNx6tSpgZzikFTfbEWNuQMKBXB3bIjc1SGiQRDir0ZWcgRe+GEStvzTvfjHKw/hy39/AP/143vwT/fFI80wBlqVD1ptdnxb0YjffFWBF/5wDA+8uR93/2IPfvLeIRTuLsOnJ2pQ2dAGQRDkPiWiIc3jHp9t27YhPz8fGzZswH333Ydf//rXyMnJwenTpxEXF9erfEVFBebNm4clS5bgww8/xNdff41nnnkGERERWLBgAQCgra0NiYmJePzxx7F06dKbHnvy5MnYu3ev+LtS6ToG/sYbb+Ctt97C//zP/yA5ORmvv/46HnroIZw9exZBQcN/IvB31WYAwLiIQC6nTzRCuRsi67J340J9K05cacLJajNOXDHjdK0FzdYulFxsQMnFBvHzOj8V7hqrw5Q7dI4/x4YgRqflnWREPRSCh/88SE9Px7Rp07Bx40Zx25133on58+ejsLCwV/nly5dj586dKCsrE7fl5eXh+PHjKCkp6VU+Pj4e+fn5yM/Pd9m+atUq7NixA8eOHXNbL0EQEBMTg/z8fCxfvhwAYLVaERUVhbVr1+KnP/3pLc/NYrFAp9PBbDYjODj4luWl9u7fy/FW0Tk8OvUOvJ17j9zVISIZddq7ce5qM77rCUInq80oq7Wg0977f+lhAWrHIo136JByhw6TY4JxRwgfwUEjhyftt0fdBjabDUePHsVLL73ksj07OxsHDx50+5mSkhJkZ2e7bJszZw42b96Mzs5OqFSqfh+/vLwcMTEx0Gg0SE9Px+rVq5GYmAjA0bNkNBpdjqXRaJCVlYWDBw+6DT5WqxVWq1X83WKx9LsucjhxxdHjk3KHTuaaEJHcVEofTI7RYXKMDrnTHdusXXacM7bgRHUTTl5xBKJzV5vR0GpzeRYZ4OgZmhwT3PNyhKHECE6gppHPo+BjMplgt9sRFRXlsj0qKgpGo9HtZ4xGo9vyXV1dMJlMiI6O7tex09PTsXXrViQnJ+Pq1at4/fXXMXPmTJw6dQphYWHi8d0d6/Lly273WVhYiF/84hf9Ov5Q4BzqmsLgQ0RuaHyV4uM3kO7Y1tFpR1mtRRwiO1VjQfnVZpjbO3HwQgMOXrg+TOanUmJidJBLGEqOCoJWxVvraeQY0ESR73ePCoLQZ5epu/LutvclJydH/HnKlCnIyMjAuHHj8MEHH6CgoGBAdXv55ZddPmuxWBAbG9vvOkmprrkDRotjYvPkmKE3DEdEQ5NWpcTUuDGYGnf9TlBrlx3lV1twqsaM76otOFVjRlltM9o77SitbEJpZZNY1tdHgfGRgZgUE4yUnjA0KSYYQdr+99YTDSUeBZ/w8HAolcpevTt1dXW9elqc9Hq92/K+vr4ICwvzsLrXBQQEYMqUKSgvLxePAzh6mG7sReqrbhqNBhqNZsB1kBInNhPRYNH4KpHSM9/HOUxm7xZQYWrFqRozTtdYcKrGEYiutXXijLEZZ4zN+Pgf1eI+DGH+Ys/QpJhgTIoORmSQhvOGaMjzqAVVq9VITU1FUVERHn30UXF7UVERHnnkEbefycjIwF//+leXbXv27EFaWppH83u+z2q1oqysDJmZmQCAhIQE6PV6FBUVYerUqQAcc5KKi4uxdu3aAR9nqDh5xTH/iMNcROQNyp6enfGRgeLdZIIgoNbcIYag76otOF1jRo25Q1x0cffJ6/+wDQ1QY6I+CHdGB4t/jo8M5FAZDSkedx0UFBRg4cKFSEtLQ0ZGBt577z1UVlYiLy8PgGP4qLq6Glu3bgXguINr3bp1KCgowJIlS1BSUoLNmzfjo48+Evdps9lw+vRp8efq6mocO3YMgYGBGD9+PABg2bJlePjhhxEXF4e6ujq8/vrrsFgsWLRoEQDHEFd+fj5Wr16NpKQkJCUlYfXq1fD398eTTz55e1dpCPiuhhObiUhaCoUCMSF+iAnxw0OTrvecX2u1iWHI+WeFqRWNrbZe84aUPgokhgdgYnQw7owOwp36YEyMDoI+mLfYkzw8Dj65ubloaGjAa6+9htraWqSkpGD37t0wGAwAgNraWlRWVorlExISsHv3bixduhTr169HTEwM3n33XXENHwCoqakRe2kA4M0338Sbb76JrKws7N+/HwBw5coV/OQnP4HJZEJERARmzJiBQ4cOiccFgBdffBHt7e145plncO3aNaSnp2PPnj0jYg2fs8ZmAOCDSYlIdmMC1Lg/KRz3J4WL2zo6HfOGyowWlNVacKa2GWVGC5raOlFe14Lyuhb89fj1fYT4q8ReIWcY4kRqkoLH6/iMZEN1HZ92mx2TXv0bBAE48vMfIjxweMxLIqLRTRAEXLVYUWbsCUK1FpwxWnChvhX27t5Nj48CSOjpHZrUM1yWHBWEO0L8+GgO6pPX1vEheZyva4EgOMbPGXqIaLhQKBTQ67TQ67R4YEKkuL2j047zdS04Y7wehspqm9HYasOF+lZcqG/FrhO1YvkAtRLjo4IwISoQyVFBmNATiDiZmgaCwWcYOHfVMcyVFBkoc02IiG6fVnX9rjInQRBQ32xFmbEZZ2p7hsuMzbhY34pWmx3Hq5pwvKrJZT86PxUmRAUhWe8IRMlRQZgQFYQxAWqJz4iGEwafYeBcnSP4JEdxfg8RjUwKhQKRwVpEBmuRlRwhbu+yd+NSQxvOXW3GWWOz48+rzbhkaoW5vRPfXmrEt5caXfYVEaRBsrN3KCoIyfogJEUGcu0hAsDgMyyUX20BACRHsceHiEYXX6WPeJv9vCnX12jr6LTjYn2rGITKe/6samxHfbMV9c1WfH2+wWVfd4T4OQKRPgjJkY4eonGRAfBXsykcTfhfexgQh7rY40NEBMAxXDapZxXpG7Vau1Be14JzN/QOnbvajKsWK6qb2lHd1I59NzyzDHAEIme4Gh8ZiKSeP0P8OWQ2EjH4DHGt1i5cudYOgENdRES3EqDxxT2xIbgnNsRlu7mtE+fqbhguMzbjQn0LTC02MRAVn3MNROGB6uuBKCIQSVFBGB8ZyEnVwxyDzxB3od4xzBUeqEYoJ+wREQ2Izl+F6fGhmB4f6rL9WqsN5+tbUH61BefrWnC+vgXnrzajxtwBU4sNppZGHLroOocoSOt7QxhyBqMgjB3D2+6HAwafIe5cz/yepEj29hARDbYxAWpMD+gdiFqtXbjgDET1PaGorgWXG1rR3NHV62GuAKBV+SAx/PqQWWJEABLDA5EQHgA/NRdmHCoYfIa48jrn/B5ObCYikkqAxhd3jQ3BXWNDXLZbu+yoMLWKQai8rgUX6lpwsb4VHZ3dOF1rwelaS6/93RHih8SIAIyLuB6IEiMCoA/WspdIYgw+Q9xlUxsAIDE8QOaaEBGRxleJifpgTNS7Tqrusnej6lq7GIjO17XgoskRiMztneI8ogPlJpfP+amUSAgPcIShiECMc/YSRQQgUMMm2ht4VYe4Sw2tAAADgw8R0ZDlq/RBQngAEsIDXB7oKggCGlttuGhqxcV6RxC6UN+Ki6YWVDa0ob3TftNeoqhgjdgzlNjTUzQuPBB3jPGDkr1EA8bgM4QJgoDLDY4eH0Oov8y1ISIiTykUCoQFahAWqOk1j6jT3o2qxjZc7AlCF+paxV6ihlYbrlqsuGqxouSi63pEal8fxIf5iz1D8WH+iA9zhK4I3nF2Sww+Q1h9sxXtnXb4KICxYxh8iIhGEpXSp6cnJxBAlMt75rZOXOgJQc6eooumFlwytcHW1Y1zV1vEm19uFKBWwhAWgPhwRxiK7+mFig8LQHigmqEIDD5D2qWe3p47xvhB7esjc22IiEgqOn8VpsWNwbS4MS7b7d0Cqq+1i6HokqkVlxocr+pr7Wi13XzoLFDjC0OYvyMM9YSi+J7fwwJGTyhi8BnCnPN74sM4v4eIiACljwJxYf6IC/PHAxNc37N22VHV2I7LDa2ocAYiUxsqTK2oMbejxdqFUzUWnKrpHYqCNL6OIBR+fejM2Vs0xl81okIRg88QVumc3xPGYS4iIuqbxlcpriH0fY5Q1IYKU5tLL9ElUxtqzO1otnbhZLUZJ6vNvT4brHWEIkNYAAyhjtBlCPWHISwAkUGaYXc7PoPPEMYeHyIiGgyOUBSE8W4Ww+3odIYiRxhyhqPLDa2oMXfA0tGFE1fMOHGldyjS+PogLtQfhjB/xIUGOP7sCUZjx/gPyWkaDD5DmHhHF4MPERF5iValRFJUkNsHYXd02nG5wRGKKhtbcbmhDZWNbbjc0IbqpnZYu7pR3rOQ4/f5KIBonR8MYd8LRqH+uDM6WLZb8hl8hihBEK6v4cOhLiIikoFWpcQEfRAm6HuHok57N2qa2nG5oQ2XG9tQ2eAIRs5w1N5pFxduPHjh+i35CgVw5j/mQukjz2M8GHyGqGttnWju6AIAxHENHyIiGmJUSh/HvB83oxKCIKC+2YrLPb1DlQ2t4s9d3d3Q+Mr37DIGnyHK2dsTrdNCq+LD7YiIaPhQKBSIDNYiMljba+FGuQ29WUcE4PodXeztISIiGjwMPkNUdVM7AMfihURERDQ4GHyGqFpzT/AJYfAhIiIaLAw+Q1RNUwcAx62ARERENDgYfIaomp6hrpgQrcw1ISIiGjkYfIYoZ/DhUBcREdHgYfAZglqsXbD0rOETzeBDREQ0aBh8hqDant6eYK0vAjVcaomIiGiwMPgMQTVmx8TmGPb2EBERDSoGnyHo+sRmBh8iIqLBxOAzBPGOLiIiIu9g8BmCuIYPERGRdzD4DEG8lZ2IiMg7GHyGIOfjKqJ1HOoiIiIaTAMKPhs2bEBCQgK0Wi1SU1Nx4MCBPssXFxcjNTUVWq0WiYmJ2LRpk8v7p06dwoIFCxAfHw+FQoF33nmn1z4KCwsxffp0BAUFITIyEvPnz8fZs2ddyixevBgKhcLlNWPGjIGcomwEQeBdXURERF7icfDZtm0b8vPzsXLlSpSWliIzMxM5OTmorKx0W76iogLz5s1DZmYmSktLsWLFCjz//PPYvn27WKatrQ2JiYlYs2YN9Hq92/0UFxfj2WefxaFDh1BUVISuri5kZ2ejtbXVpdzcuXNRW1srvnbv3u3pKcqqodUGW1c3FApAzx4fIiKiQaUQBEHw5APp6emYNm0aNm7cKG678847MX/+fBQWFvYqv3z5cuzcuRNlZWXitry8PBw/fhwlJSW9ysfHxyM/Px/5+fl91qO+vh6RkZEoLi7GrFmzADh6fJqamrBjxw5PTklksVig0+lgNpsRHBw8oH3cru+qzfh//r+vEBGkweGVP5SlDkRERMOJJ+23Rz0+NpsNR48eRXZ2tsv27OxsHDx40O1nSkpKepWfM2cOjhw5gs7OTk8O78JsNgMAQkNDXbbv378fkZGRSE5OxpIlS1BXVzfgY8ihvtkKAIgM0shcEyIiopHHo+chmEwm2O12REVFuWyPioqC0Wh0+xmj0ei2fFdXF0wmE6Kjoz2ssmMeTEFBAe6//36kpKSI23NycvD444/DYDCgoqICr7zyCh588EEcPXoUGk3vIGG1WmG1WsXfLRaLx3UZbM7gE8HgQ0RENOgG9CAohULh8rsgCL223aq8u+399dxzz+HEiRP46quvXLbn5uaKP6ekpCAtLQ0GgwG7du3CY4891ms/hYWF+MUvfjGgOnhLfUtP8Alk8CEiIhpsHg11hYeHQ6lU9urdqaur69Wr46TX692W9/X1RVhYmIfVBX72s59h586d2LdvH8aOHdtn2ejoaBgMBpSXl7t9/+WXX4bZbBZfVVVVHtdnsIlDXcEMPkRERIPNo+CjVquRmpqKoqIil+1FRUWYOXOm289kZGT0Kr9nzx6kpaVBpVL1+9iCIOC5557Dxx9/jC+++AIJCQm3/ExDQwOqqqpuOpym0WgQHBzs8pKbONTFHh8iIqJB5/Ht7AUFBfjNb36D999/H2VlZVi6dCkqKyuRl5cHwNGL8tRTT4nl8/LycPnyZRQUFKCsrAzvv/8+Nm/ejGXLlollbDYbjh07hmPHjsFms6G6uhrHjh3D+fPnxTLPPvssPvzwQ/z+979HUFAQjEYjjEYj2tsdi/21tLRg2bJlKCkpwaVLl7B//348/PDDCA8Px6OPPjrgCyS163N8eCs7ERHRoBMGYP369YLBYBDUarUwbdo0obi4WHxv0aJFQlZWlkv5/fv3C1OnThXUarUQHx8vbNy40eX9iooKAUCv1437cfc+AGHLli2CIAhCW1ubkJ2dLURERAgqlUqIi4sTFi1aJFRWVvb7vMxmswBAMJvNHl+TwTL7l/sEw/JPhW8uNshWByIiouHEk/bb43V8RrKhsI5Pyqufo8XahX3LZiMhPECWOhAREQ0nXlvHh7yrzdaFFmsXAN7OTkRE5A0MPkOIc36Pn0qJALVS5toQERGNPAw+Q8iNixcOdI0jIiIiujkGnyGEqzYTERF5F4PPEOJctZnP6SIiIvIOBp8hhD0+RERE3sXgM4Rw1WYiIiLvYvAZQtjjQ0RE5F0MPkOI+GR2Bh8iIiKvYPAZQuosDD5ERETexOAzRAiCgIZWR/AJ5xwfIiIir2DwGSKarV3otDsemxYaoJa5NkRERCMTg88Q0dTaCcDxuAqtio+rICIi8gYGnyHiWpsNADDGXyVzTYiIiEYuBp8hotEZfDjMRURE5DUMPkNEk9jjw+BDRETkLQw+Q0RjzxyfEA51EREReQ2DzxDh7PHhHV1ERETew+AzRDgnN4dwqIuIiMhrGHyGiGs9Q12hHOoiIiLyGgafIeIa7+oiIiLyOgafIaKxlUNdRERE3sbgM0Q0tTmHuhh8iIiIvIXBZwgQBOGGyc2c40NEROQtDD5DQHunHdaubgCc40NERORNDD5DwLWeYS610gcBaj6glIiIyFsYfIaAa63Xh7kUCoXMtSEiIhq5GHyGgGtctZmIiEgSDD5DgHOoixObiYiIvIvBZwhwDnXxyexERETexeAzBHDVZiIiImkw+AwB13t8ONRFRETkTQw+Q4Bzjg+HuoiIiLyLwWcIEIe6GHyIiIi8isFnCGjiXV1ERESSYPAZApo7HMEn2I/Bh4iIyJsGFHw2bNiAhIQEaLVapKam4sCBA32WLy4uRmpqKrRaLRITE7Fp0yaX90+dOoUFCxYgPj4eCoUC77zzzoCOKwgCVq1ahZiYGPj5+WH27Nk4derUQE5RUs0dXQCAYC2DDxERkTd5HHy2bduG/Px8rFy5EqWlpcjMzEROTg4qKyvdlq+oqMC8efOQmZmJ0tJSrFixAs8//zy2b98ulmlra0NiYiLWrFkDvV4/4OO+8cYbeOutt7Bu3TocPnwYer0eDz30EJqbmz09TckIggBLT49PkNZX5toQERGNcIKH7r33XiEvL89l28SJE4WXXnrJbfkXX3xRmDhxosu2n/70p8KMGTPcljcYDMLbb7/t8XG7u7sFvV4vrFmzRny/o6ND0Ol0wqZNm255XoIgCGazWQAgmM3mfpUfDO22LsGw/FPBsPxTwdJuk+y4REREI4Un7bdHPT42mw1Hjx5Fdna2y/bs7GwcPHjQ7WdKSkp6lZ8zZw6OHDmCzs7OQTtuRUUFjEajSxmNRoOsrKyb1s1qtcJisbi8pObs7VEogAA1e3yIiIi8yaPgYzKZYLfbERUV5bI9KioKRqPR7WeMRqPb8l1dXTCZTIN2XOefntStsLAQOp1OfMXGxvarPoPJ0u6Y3xOk8YWPD5/MTkRE5E0DmtysULg20IIg9Np2q/Lutg/GcT2p28svvwyz2Sy+qqqqPKrPYGgW5/dwYjMREZG3eTS2Eh4eDqVS2asHpa6urldPi5Ner3db3tfXF2FhYYN2XOekaKPRiOjo6H7VTaPRQKPR9KsO3uK8o4sTm4mIiLzPox4ftVqN1NRUFBUVuWwvKirCzJkz3X4mIyOjV/k9e/YgLS0NKlX/ejn6c9yEhATo9XqXMjabDcXFxTet21Bg4Ro+REREkvG4m6GgoAALFy5EWloaMjIy8N5776GyshJ5eXkAHMNH1dXV2Lp1KwAgLy8P69atQ0FBAZYsWYKSkhJs3rwZH330kbhPm82G06dPiz9XV1fj2LFjCAwMxPjx4/t1XIVCgfz8fKxevRpJSUlISkrC6tWr4e/vjyeffPL2rpIXXV/Dhz0+RERE3uZxa5ubm4uGhga89tprqK2tRUpKCnbv3g2DwQAAqK2tdVlbJyEhAbt378bSpUuxfv16xMTE4N1338WCBQvEMjU1NZg6dar4+5tvvok333wTWVlZ2L9/f7+OCwAvvvgi2tvb8cwzz+DatWtIT0/Hnj17EBQU5PGFkQrn+BAREUlHIThnGhMsFgt0Oh3MZjOCg4MlOeabn5/Fun3nsSjDgF88kiLJMYmIiEYST9pvPqtLZuzxISIikg6Dj8x4VxcREZF0GHxkxru6iIiIpMPgIzMLe3yIiIgkw+Ajs+tDXezxISIi8jYGH5lZ2nuGutjjQ0RE5HUMPjLjXV1ERETSYfCRUXe3gBYrV24mIiKSCoOPjFptXejuWT6Sd3URERF5H4OPjJwTm1VKBTS+/E9BRETkbWxtZXTjHV0KhULm2hAREY18DD4yEhcv5PweIiIiSTD4yIh3dBEREUmLwUdGfE4XERGRtBh8ZHR98UL2+BAREUmBwUdGfE4XERGRtBh8ZMTndBEREUmLwUdG1yc3s8eHiIhICgw+Mmqz2QEAgRoGHyIiIikw+Miotec5Xf4apcw1ISIiGh0YfGTk7PEJULPHh4iISAoMPjJqtfX0+KjZ40NERCQFBh8ZtVkdPT7+7PEhIiKSBIOPjMQeH87xISIikgSDj4w4x4eIiEhaDD4yauMcHyIiIkkx+MjE3i2go7MbABDAdXyIiIgkweAjE2dvD8AeHyIiIqkw+MjEOb9H6aOAxpf/GYiIiKTAFlcm4qrNaiUUCoXMtSEiIhodGHxk4uzx4TAXERGRdBh8ZOLs8eGt7ERERNJh8JFJW2dPjw8XLyQiIpIMg49M+LgKIiIi6TH4yMT5uIoAzvEhIiKSDIOPTNqcd3Vx8UIiIiLJDCj4bNiwAQkJCdBqtUhNTcWBAwf6LF9cXIzU1FRotVokJiZi06ZNvcps374dkyZNgkajwaRJk/DJJ5+4vB8fHw+FQtHr9eyzz4plFi9e3Ov9GTNmDOQUva5VfE4Xe3yIiIik4nHw2bZtG/Lz87Fy5UqUlpYiMzMTOTk5qKysdFu+oqIC8+bNQ2ZmJkpLS7FixQo8//zz2L59u1impKQEubm5WLhwIY4fP46FCxfiiSeewDfffCOWOXz4MGpra8VXUVERAODxxx93Od7cuXNdyu3evdvTU5TE9ed0sceHiIhIKgpBEARPPpCeno5p06Zh48aN4rY777wT8+fPR2FhYa/yy5cvx86dO1FWViZuy8vLw/Hjx1FSUgIAyM3NhcViwWeffSaWmTt3LsaMGYOPPvrIbT3y8/Px6aefory8XFwAcPHixWhqasKOHTs8OSWRxWKBTqeD2WxGcHDwgPbRX6t2nsL/HLyEZ2aPw4tzJ3r1WERERCOZJ+23Rz0+NpsNR48eRXZ2tsv27OxsHDx40O1nSkpKepWfM2cOjhw5gs7Ozj7L3GyfNpsNH374IZ5++uleqx7v378fkZGRSE5OxpIlS1BXV3fT87FarbBYLC4vqTh7fPiAUiIiIul4FHxMJhPsdjuioqJctkdFRcFoNLr9jNFodFu+q6sLJpOpzzI32+eOHTvQ1NSExYsXu2zPycnB7373O3zxxRf41a9+hcOHD+PBBx+E1Wp1u5/CwkLodDrxFRsbe9NzH2xcuZmIiEh6A+pu+H4viyAIfT5vyl3572/3ZJ+bN29GTk4OYmJiXLbn5uaKP6ekpCAtLQ0GgwG7du3CY4891ms/L7/8MgoKCsTfLRaLZOGnTZzczB4fIiIiqXjU6oaHh0OpVPbqiamrq+vVY+Ok1+vdlvf19UVYWFifZdzt8/Lly9i7dy8+/vjjW9Y3OjoaBoMB5eXlbt/XaDTQaDS33I83iA8p5crNREREkvFoqEutViM1NVW8o8qpqKgIM2fOdPuZjIyMXuX37NmDtLQ0qFSqPsu42+eWLVsQGRmJH/3oR7esb0NDA6qqqhAdHX3LslJjjw8REZH0PL6dvaCgAL/5zW/w/vvvo6ysDEuXLkVlZSXy8vIAOIaPnnrqKbF8Xl4eLl++jIKCApSVleH999/H5s2bsWzZMrHMCy+8gD179mDt2rU4c+YM1q5di7179yI/P9/l2N3d3diyZQsWLVoEX1/XwNDS0oJly5ahpKQEly5dwv79+/Hwww8jPDwcjz76qKen6XWt4u3s7PEhIiKSisfdDbm5uWhoaMBrr72G2tpapKSkYPfu3TAYDACA2tpalzV9EhISsHv3bixduhTr169HTEwM3n33XSxYsEAsM3PmTPzhD3/Az3/+c7zyyisYN24ctm3bhvT0dJdj7927F5WVlXj66ad71UupVOLkyZPYunUrmpqaEB0djQceeADbtm1DUFCQp6fpdXxWFxERkfQ8XsdnJJNyHZ8pqz5Hc0cX/v5/sjAuItCrxyIiIhrJvLaODw0OQRDQzjk+REREkmPwkYHN3o2ubkdHG+/qIiIikg6Djwyc83sAwF/F4ENERCQVBh8ZOO/o0vj6wFfJ/wRERERSYasrA3ENHz6ni4iISFIMPjIQV23mGj5ERESSYvCRAR9QSkREJA8GHxlcDz4c6iIiIpISg48M2nomNwfwVnYiIiJJMfjIoLXndnY/FXt8iIiIpMTgI4P2zp7gwzk+REREkmLwkUFHT/DR+vLyExERSYktrwyszuDDVZuJiIgkxeAjA2tXNwBAq+LlJyIikhJbXhl0sMeHiIhIFgw+MujodPb4MPgQERFJicFHBh1djh4fDSc3ExERSYotrwycQ10a9vgQERFJisFHBuJQF3t8iIiIJMWWVwac3ExERCQPBh8ZdHRxcjMREZEcGHxkcH0BQ15+IiIiKbHllQGHuoiIiOTB4CMDceVmXwYfIiIiKTH4yKCDQ11ERESyYMsrA67cTEREJA8GH4kJgnB95Wb2+BAREUmKLa/EbPZuCILjZw3n+BAREUmKwUdizmEugHN8iIiIpMaWV2LONXwUCkCt5OUnIiKSElteiV1/TpcSCoVC5toQERGNLgw+EnNObOYwFxERkfTY+kqMqzYTERHJh8FHYlY+oJSIiEg2DD4Sc/b4aHx56YmIiKTG1ldiXLWZiIhIPgMKPhs2bEBCQgK0Wi1SU1Nx4MCBPssXFxcjNTUVWq0WiYmJ2LRpU68y27dvx6RJk6DRaDBp0iR88sknLu+vWrUKCoXC5aXX613KCIKAVatWISYmBn5+fpg9ezZOnTo1kFP0Gj6ni4iISD4et77btm1Dfn4+Vq5cidLSUmRmZiInJweVlZVuy1dUVGDevHnIzMxEaWkpVqxYgeeffx7bt28Xy5SUlCA3NxcLFy7E8ePHsXDhQjzxxBP45ptvXPY1efJk1NbWiq+TJ0+6vP/GG2/grbfewrp163D48GHo9Xo89NBDaG5u9vQ0veb6UBd7fIiIiKSmEATnAxT6Jz09HdOmTcPGjRvFbXfeeSfmz5+PwsLCXuWXL1+OnTt3oqysTNyWl5eH48ePo6SkBACQm5sLi8WCzz77TCwzd+5cjBkzBh999BEAR4/Pjh07cOzYMbf1EgQBMTExyM/Px/LlywEAVqsVUVFRWLt2LX7605/e8twsFgt0Oh3MZjOCg4NvfTEG4LeHLuOVHd9hzuQo/HphmleOQURENJp40n571ONjs9lw9OhRZGdnu2zPzs7GwYMH3X6mpKSkV/k5c+bgyJEj6Ozs7LPM9/dZXl6OmJgYJCQk4Mc//jEuXrwovldRUQGj0eiyH41Gg6ysrJvWzWq1wmKxuLy8zcrb2YmIiGTjUfAxmUyw2+2Iiopy2R4VFQWj0ej2M0aj0W35rq4umEymPsvcuM/09HRs3boVn3/+Of77v/8bRqMRM2fORENDg7gP5+f6W7fCwkLodDrxFRsbe6tLcNvEOT4c6iIiIpLcgGbYfv9RC4Ig9Pn4BXflv7/9VvvMycnBggULMGXKFPzwhz/Erl27AAAffPDBgOv28ssvw2w2i6+qqqqbnsNguX5XFyc3ExERSc3Xk8Lh4eFQKpW9elDq6up69bQ46fV6t+V9fX0RFhbWZ5mb7RMAAgICMGXKFJSXl4v7ABw9P9HR0f3aj0ajgUajuekxvIErNxMREcnHo24HtVqN1NRUFBUVuWwvKirCzJkz3X4mIyOjV/k9e/YgLS0NKpWqzzI32yfgmJ9TVlYmhpyEhATo9XqX/dhsNhQXF/e5H6k5V27WMPgQERFJzqMeHwAoKCjAwoULkZaWhoyMDLz33nuorKxEXl4eAMfwUXV1NbZu3QrAcQfXunXrUFBQgCVLlqCkpASbN28W79YCgBdeeAGzZs3C2rVr8cgjj+Avf/kL9u7di6+++koss2zZMjz88MOIi4tDXV0dXn/9dVgsFixatAiAY4grPz8fq1evRlJSEpKSkrB69Wr4+/vjySefvK2LNJi4jg8REZF8PA4+ubm5aGhowGuvvYba2lqkpKRg9+7dMBgMAIDa2lqXNX0SEhKwe/duLF26FOvXr0dMTAzeffddLFiwQCwzc+ZM/OEPf8DPf/5zvPLKKxg3bhy2bduG9PR0scyVK1fwk5/8BCaTCREREZgxYwYOHTokHhcAXnzxRbS3t+OZZ57BtWvXkJ6ejj179iAoKGhAF8cbOpzP6uLkZiIiIsl5vI7PSCbFOj5Lth5B0emrWP3oFDyZHueVYxAREY0mXlvHh24fh7qIiIjkw9ZXYtae29n5yAoiIiLpMfhIrKOLPT5ERERyYesrMa7jQ0REJB8GH4lx5WYiIiL5sPWVmLPHh3N8iIiIpMfgIzHnys0c6iIiIpIeg4/EeDs7ERGRfNj6SkgQBPb4EBERyYjBR0LO0AMw+BAREcmBwUdCzmEuAND68tITERFJja2vhJy3sit9FPBV8tITERFJja2vhMSJzeztISIikgVbYAldf1wF5/cQERHJgcFHQtZO3tFFREQkJwYfCdnsjuCj5lAXERGRLNgCS8jWczu7mhObiYiIZMEWWELOHh+Vr0LmmhAREY1ODD4SYo8PERGRvNgCS6jT2ePD4ENERCQLtsASEnt8OLmZiIhIFmyBJeTs8eFQFxERkTzYAkvIZhcAsMeHiIhILmyBJeQc6uIcHyIiInmwBZZQJxcwJCIikhVbYAmxx4eIiEhebIEldH1yMxcwJCIikgODj4R4OzsREZG82AJLyMYFDImIiGTFFlhC7PEhIiKSF1tgCfGRFURERPJiCywhZ4+Phj0+REREsmALLKHOnpWb2eNDREQkD7bAErJyjg8REZGs2AJLiHN8iIiI5MUWWEK8q4uIiEheA2qBN2zYgISEBGi1WqSmpuLAgQN9li8uLkZqaiq0Wi0SExOxadOmXmW2b9+OSZMmQaPRYNKkSfjkk09c3i8sLMT06dMRFBSEyMhIzJ8/H2fPnnUps3jxYigUCpfXjBkzBnKKXsGVm4mIiOTlcfDZtm0b8vPzsXLlSpSWliIzMxM5OTmorKx0W76iogLz5s1DZmYmSktLsWLFCjz//PPYvn27WKakpAS5ublYuHAhjh8/joULF+KJJ57AN998I5YpLi7Gs88+i0OHDqGoqAhdXV3Izs5Ga2ury/Hmzp2L2tpa8bV7925PT9Fr+JBSIiIieSkEQRA8+UB6ejqmTZuGjRs3itvuvPNOzJ8/H4WFhb3KL1++HDt37kRZWZm4LS8vD8ePH0dJSQkAIDc3FxaLBZ999plYZu7cuRgzZgw++ugjt/Wor69HZGQkiouLMWvWLACOHp+mpibs2LHDk1MSWSwW6HQ6mM1mBAcHD2gffZn7zpc4Y2zGb//5XmQmRQz6/omIiEYjT9pvj7oebDYbjh49iuzsbJft2dnZOHjwoNvPlJSU9Co/Z84cHDlyBJ2dnX2Wudk+AcBsNgMAQkNDXbbv378fkZGRSE5OxpIlS1BXV3fTfVitVlgsFpeXN10f6mKPDxERkRw8aoFNJhPsdjuioqJctkdFRcFoNLr9jNFodFu+q6sLJpOpzzI326cgCCgoKMD999+PlJQUcXtOTg5+97vf4YsvvsCvfvUrHD58GA8++CCsVqvb/RQWFkKn04mv2NjYvi/AbRKf1cWhLiIiIln4DuRDCoXr5FxBEHptu1X572/3ZJ/PPfccTpw4ga+++sple25urvhzSkoK0tLSYDAYsGvXLjz22GO99vPyyy+joKBA/N1isXg1/HR2Oc6bPT5ERETy8Cj4hIeHQ6lU9uqJqaur69Vj46TX692W9/X1RVhYWJ9l3O3zZz/7GXbu3Ikvv/wSY8eO7bO+0dHRMBgMKC8vd/u+RqOBRqPpcx+DycbJzURERLLyqAVWq9VITU1FUVGRy/aioiLMnDnT7WcyMjJ6ld+zZw/S0tKgUqn6LHPjPgVBwHPPPYePP/4YX3zxBRISEm5Z34aGBlRVVSE6Orpf5+dtnV1cwJCIiEhOHrfABQUF+M1vfoP3338fZWVlWLp0KSorK5GXlwfAMXz01FNPieXz8vJw+fJlFBQUoKysDO+//z42b96MZcuWiWVeeOEF7NmzB2vXrsWZM2ewdu1a7N27F/n5+WKZZ599Fh9++CF+//vfIygoCEajEUajEe3t7QCAlpYWLFu2DCUlJbh06RL279+Phx9+GOHh4Xj00UcHen0GlZU9PkRERPISBmD9+vWCwWAQ1Gq1MG3aNKG4uFh8b9GiRUJWVpZL+f379wtTp04V1Gq1EB8fL2zcuLHXPv/0pz8JEyZMEFQqlTBx4kRh+/btLu8DcPvasmWLIAiC0NbWJmRnZwsRERGCSqUS4uLihEWLFgmVlZX9Pi+z2SwAEMxmc/8vRj91d3cL8S99KhiWfypctbQP+v6JiIhGK0/ab4/X8RnJvLmOT6e9G0krHesUHf+/2dD5qwZ1/0RERKOV19bxoYFzruEDACpfPrKCiIhIDgw+EnE+oBTg7exERERyYQssEeet7AoFoPRhjw8REZEcGHwk4uzxUSt9+lzskYiIiLyHwUcinXau2kxERCQ3tsIS6eQaPkRERLJjKywRG1dtJiIikh1bYYnwOV1ERETyYysskes9PpzYTEREJBcGH4k45/hwqIuIiEg+bIUl4uzx0XCoi4iISDZshSXCHh8iIiL5sRWWiLWLk5uJiIjkxlZYIs4FDNnjQ0REJB+2whKxsceHiIhIdmyFJSKu3MweHyIiItmwFZYIe3yIiIjkx1ZYIjY7FzAkIiKSG4OPRNjjQ0REJD+2whLhOj5ERETyYysskU4+pJSIiEh2bIUlIg51sceHiIhINmyFJWLrWcCQwYeIiEg+bIUl4uzxUXGoi4iISDZshSXCyc1ERETyYyssEd7OTkREJD+2whK5/sgKLmBIREQkFwYfidh4OzsREZHs2ApLRJzczDk+REREsmErLBEbn85OREQkO7bCEhHv6uJQFxERkWzYCkvEOdSlYY8PERGRbNgKS6SzZ+Vm9vgQERHJh62wRPisLiIiIvmxFZaIjSs3ExERyW5ArfCGDRuQkJAArVaL1NRUHDhwoM/yxcXFSE1NhVarRWJiIjZt2tSrzPbt2zFp0iRoNBpMmjQJn3zyicfHFQQBq1atQkxMDPz8/DB79mycOnVqIKc46Dq5jg8REZHsPG6Ft23bhvz8fKxcuRKlpaXIzMxETk4OKisr3ZavqKjAvHnzkJmZidLSUqxYsQLPP/88tm/fLpYpKSlBbm4uFi5ciOPHj2PhwoV44okn8M0333h03DfeeANvvfUW1q1bh8OHD0Ov1+Ohhx5Cc3Ozp6c56DjURURENAQIHrr33nuFvLw8l20TJ04UXnrpJbflX3zxRWHixIku2376058KM2bMEH9/4oknhLlz57qUmTNnjvDjH/+438ft7u4W9Hq9sGbNGvH9jo4OQafTCZs2berXuZnNZgGAYDab+1XeE+NX7BIMyz8VapvaB33fREREo5kn7bdH3Q82mw1Hjx5Fdna2y/bs7GwcPHjQ7WdKSkp6lZ8zZw6OHDmCzs7OPss499mf41ZUVMBoNLqU0Wg0yMrKumndrFYrLBaLy8sburuF63d18VldREREsvEo+JhMJtjtdkRFRblsj4qKgtFodPsZo9HotnxXVxdMJlOfZZz77M9xnX96UrfCwkLodDrxFRsbe9Nzvx2d3d3iz7ydnYiISD4DaoUVCtdeC0EQem27Vfnvb+/PPgerjNPLL78Ms9ksvqqqqm56DrdDAQWef3A88rLGQeur9MoxiIiI6NZ8PSkcHh4OpVLZqwelrq6uV0+Lk16vd1ve19cXYWFhfZZx7rM/x9Xr9QAcPT/R0dH9qptGo4FGo+nznAeD2tcHBdkTvH4cIiIi6ptHPT5qtRqpqakoKipy2V5UVISZM2e6/UxGRkav8nv27EFaWhpUKlWfZZz77M9xExISoNfrXcrYbDYUFxfftG5EREQ0yng6c/oPf/iDoFKphM2bNwunT58W8vPzhYCAAOHSpUuCIAjCSy+9JCxcuFAsf/HiRcHf319YunSpcPr0aWHz5s2CSqUS/vznP4tlvv76a0GpVApr1qwRysrKhDVr1gi+vr7CoUOH+n1cQRCENWvWCDqdTvj444+FkydPCj/5yU+E6OhowWKx9OvcvHlXFxEREXmHJ+23x8FHEARh/fr1gsFgENRqtTBt2jShuLhYfG/RokVCVlaWS/n9+/cLU6dOFdRqtRAfHy9s3Lix1z7/9Kc/CRMmTBBUKpUwceJEYfv27R4dVxAct7S/+uqrgl6vFzQajTBr1izh5MmT/T4vBh8iIqLhx5P2WyEIPTONCRaLBTqdDmazGcHBwXJXh4iIiPrBk/ab91YTERHRqMHgQ0RERKMGgw8RERGNGgw+RERENGow+BAREdGoweBDREREowaDDxEREY0aDD5EREQ0ajD4EBER0ajh0dPZRzrnItYWi0XmmhAREVF/Odvt/jyMgsHnBs3NzQCA2NhYmWtCREREnmpuboZOp+uzDJ/VdYPu7m7U1NQgKCgICoViUPdtsVgQGxuLqqoqPgfMy3itpcHrLB1ea2nwOktnsK+1IAhobm5GTEwMfHz6nsXDHp8b+Pj4YOzYsV49RnBwMP9CSYTXWhq8ztLhtZYGr7N0BvNa36qnx4mTm4mIiGjUYPAhIiKiUYPBRyIajQavvvoqNBqN3FUZ8XitpcHrLB1ea2nwOktHzmvNyc1EREQ0arDHh4iIiEYNBh8iIiIaNRh8iIiIaNRg8CEiIqJRg8FHAhs2bEBCQgK0Wi1SU1Nx4MABuas04qxatQoKhcLlpdfr5a7WiPDll1/i4YcfRkxMDBQKBXbs2OHyviAIWLVqFWJiYuDn54fZs2fj1KlT8lR2GLvVdV68eHGv7/iMGTPkqewwVlhYiOnTpyMoKAiRkZGYP38+zp4961KG3+nB0Z9rLcf3msHHy7Zt24b8/HysXLkSpaWlyMzMRE5ODiorK+Wu2ogzefJk1NbWiq+TJ0/KXaURobW1FXfffTfWrVvn9v033ngDb731FtatW4fDhw9Dr9fjoYceEp99R/1zq+sMAHPnznX5ju/evVvCGo4MxcXFePbZZ3Ho0CEUFRWhq6sL2dnZaG1tFcvwOz04+nOtARm+1wJ51b333ivk5eW5bJs4caLw0ksvyVSjkenVV18V7r77brmrMeIBED755BPx9+7ubkGv1wtr1qwRt3V0dAg6nU7YtGmTDDUcGb5/nQVBEBYtWiQ88sgjstRnJKurqxMACMXFxYIg8DvtTd+/1oIgz/eaPT5eZLPZcPToUWRnZ7tsz87OxsGDB2Wq1chVXl6OmJgYJCQk4Mc//jEuXrwod5VGvIqKChiNRpfvuEajQVZWFr/jXrB//35ERkYiOTkZS5YsQV1dndxVGvbMZjMAIDQ0FAC/0970/WvtJPX3msHHi0wmE+x2O6Kioly2R0VFwWg0ylSrkSk9PR1bt27F559/jv/+7/+G0WjEzJkz0dDQIHfVRjTn95jfce/LycnB7373O3zxxRf41a9+hcOHD+PBBx+E1WqVu2rDliAIKCgowP3334+UlBQA/E57i7trDcjzvebT2SWgUChcfhcEodc2uj05OTniz1OmTEFGRgbGjRuHDz74AAUFBTLWbHTgd9z7cnNzxZ9TUlKQlpYGg8GAXbt24bHHHpOxZsPXc889hxMnTuCrr77q9R6/04PrZtdaju81e3y8KDw8HEqlste/Eurq6nr9a4IGV0BAAKZMmYLy8nK5qzKiOe+c43dcetHR0TAYDPyOD9DPfvYz7Ny5E/v27cPYsWPF7fxOD76bXWt3pPheM/h4kVqtRmpqKoqKily2FxUVYebMmTLVanSwWq0oKytDdHS03FUZ0RISEqDX612+4zabDcXFxfyOe1lDQwOqqqr4HfeQIAh47rnn8PHHH+OLL75AQkKCy/v8Tg+eW11rd6T4XnOoy8sKCgqwcOFCpKWlISMjA++99x4qKyuRl5cnd9VGlGXLluHhhx9GXFwc6urq8Prrr8NisWDRokVyV23Ya2lpwfnz58XfKyoqcOzYMYSGhiIuLg75+flYvXo1kpKSkJSUhNWrV8Pf3x9PPvmkjLUefvq6zqGhoVi1ahUWLFiA6OhoXLp0CStWrEB4eDgeffRRGWs9/Dz77LP4/e9/j7/85S8ICgoSe3Z0Oh38/PygUCj4nR4kt7rWLS0t8nyvJb2HbJRav369YDAYBLVaLUybNs3lVj4aHLm5uUJ0dLSgUqmEmJgY4bHHHhNOnTold7VGhH379gkAer0WLVokCILj9t9XX31V0Ov1gkajEWbNmiWcPHlS3koPQ31d57a2NiE7O1uIiIgQVCqVEBcXJyxatEiorKyUu9rDjrtrDEDYsmWLWIbf6cFxq2st1/da0VM5IiIiohGPc3yIiIho1GDwISIiolGDwYeIiIhGDQYfIiIiGjUYfIiIiGjUYPAhIiKiUYPBh4iIiEYNBh8iIiIaNRh8iIiIaNRg8CEiIqJRg8GHiIiIRg0GHyIiIho1/n8+GkWZFv+VbwAAAABJRU5ErkJggg==",
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"a, b = 0, 25\n",
"h = 0.1\n",
"y0 = 0\n",
"t, y = eulerimplicite(f_edo, a, b, y0, h)\n",
"\n",
"plt.plot(t, y, label='Euler Implicite')\n",
"\n",
"def f(t):\n",
" return np.sqrt(2*t+1)-1\n",
"\n",
"plt.plot(t, f(t), label='Sol analytique')\n",
"\n",
"plt.legend()\n",
"plt.show()\n",
"\n",
"plt.plot(t, np.abs(f(t)-y), label=\"Erreur\")\n",
"plt.legend(loc='upper right')\n",
"plt.show()"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "cd3c87b5-2624-4270-aa03-0a52c6ad2495",
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.11.5"
}
},
"nbformat": 4,
"nbformat_minor": 5
}