{ "cells": [ { "cell_type": "code", "execution_count": 16, "id": "4236a4d0-f996-40b3-832a-263753d3707b", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAfAAAAGiCAYAAAAY6pQaAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/NK7nSAAAACXBIWXMAAA9hAAAPYQGoP6dpAABmXklEQVR4nO29fXQc1Znn/63qbrXkF8nYYMkKslFYT0x4GYwBI2CTXdDGywBrBp8QZ51dB9gwSWyC7d+G4FlsWPNiYM6AQ+LYgcM4MEAYmAkkkBOzRNk4S2KMcYYEQmKciQcrGMnDiyRbdre6u+r3h0L3vU9131J1dUsq9fdzTp/Tt++91bdLL7fr+db3eSzXdV0QQgghJFLYY70AQgghhASHGzghhBASQbiBE0IIIRGEGzghhBASQbiBE0IIIRGEGzghhBASQbiBE0IIIRGEGzghhBASQbiBE0IIIRGEGzghhBASQbiBE0IIIQH42c9+hssvvxytra2wLAvPPPOM1u+6LtavX49Zs2ahoaEBnZ2d2Ldvnzbm/fffx7Jly9DY2Ihp06bh2muvxZEjRwKtgxs4IYQQEoDBwUH8+Z//OTZv3ly0/5577sH999+PrVu3YteuXZg8eTIWLVqEVCqVH7Ns2TL85je/wQsvvIDnnnsOP/vZz3DdddcFW4hbJb75zW+6c+bMcZPJpHvuuee6u3btqtZbEUIIIWMCAPfpp5/Otx3HcVtaWty/+Zu/yb/W19fnJpNJ97vf/a7ruq77xhtvuADc3bt358f86Ec/ci3Lct9+++0Rv3e8/O8gpfmHf/gHrFmzBlu3bsXChQuxadMmLFq0CHv37sXMmTONcx3HwcGDBzF16lRYllWN5RFCCKkiruvi8OHDaG1thW1XL9CbSqUwNDQU+jiu63r2m2QyiWQyGfhY+/fvR09PDzo7O/OvNTU1YeHChdi5cyeWLl2KnTt3Ytq0aTj77LPzYzo7O2HbNnbt2oW//Mu/HNF7VWUDv/fee/GFL3wBV199NQBg69at+OEPf4i/+7u/w0033WSce/DgQbS1tVVjWYQQQkaR7u5unHjiiVU5diqVwvENDRiswLGmTJni0Z9vueUW3HrrrYGP1dPTAwBobm7WXm9ubs739fT0eC5m4/E4pk+fnh8zEiq+gQ8NDWHPnj1Yu3Zt/jXbttHZ2YmdO3d6xqfTaaTT6Xzb/VN58pU2kPzTN6L/03hupZdJCCGkSuTcLN7o24OpU6dW7T2GhoYwCGBV3ELw6+QCaQCbjhxBd3c3Ghsb86+Xc/U92lR8A3/33XeRy+WKfvv43e9+5xm/ceNG/O///b89ryctK7+Bx6yqBAoIIYRUkdGQQZMoXOyVxZ8uGhsbG7UNvFxaWloAAL29vZg1a1b+9d7eXpx55pn5MYcOHdLmZbNZvP/++/n5I2HMd8a1a9dizZo1+fbAwIAnhL76UEFD+dq+v9X6Tt51ltaevF/RQ959T+tz+/5NazuDfXp/5lihL9Wv9w0d1du5IeV5RhyntBVAnVe8nUYp5Fhvf+m5nrHOyHWjIMcdjeMUw6nisUeKm0v5D4oQ1fx5EVIpbAuIhdi/K63Qt7e3o6WlBV1dXfkNe2BgALt27cKXvvQlAEBHRwf6+vqwZ88eLFiwAADwk5/8BI7jYOHChSN+r4pv4McffzxisRh6e3u113t7e4t+syj3RgFCCCEkbg0/yp5fxpwjR47g97//fb69f/9+vPrqq5g+fTpmz56NVatW4fbbb8fcuXPR3t6OdevWobW1FVdccQUA4JRTTsF//s//GV/4whewdetWZDIZrFy5EkuXLkVra+uI11Hx2wPr6uqwYMECdHV15V9zHAddXV3o6Oio9NsRQgipYWJW+EdQXnnlFcyfPx/z588HAKxZswbz58/H+vXrAQA33ngjrr/+elx33XU455xzcOTIEWzfvh319fX5Yzz22GOYN28eLr74YvzFX/wFLrzwQjzwwAOB1mG5H941VkH+4R/+AcuXL8e3v/1tnHvuudi0aROefPJJ/O53v/No45KBgQE0NTXh/4sVNPDbHvnzfP+jV03Rxn998L9q7XcOnJd/Pr37JK3vuHd0faPhXf2jW0eUEOhhEQZP6eFRN3VUea7fB6mG4gEAmcJcT7hdhuYdvR/KeM9cT/hdGesXIjeE4/3C/HpfgLC9jwRQ7nGBYJJAmPcZrWOpjJU8MJEkAcoBo0vOzeK1D3ahv7+/IrpyMT7cK25NWqgPoYGnXBe3pt2qrrVaVEUD/8xnPoN/+7d/w/r169HT04MzzzwT27dv9928CSGEkCCUexWdn1+5pYw6VbuJbeXKlVi5cmW1Dk8IIYQgZrkhN/CKB6FHDeZCJ4QQQiLImNvIRsLg/+nOP7/01Ll656mPa81//LO3889fPf4cre/9ttlae+oh/W6/xvea8s8n9etZchKH9W9p9tGs8lzoa2nRVvXzjNCXh8TYjNDaM0qSG6GBe8Ya9HJI/VzV2j3aullrN441adG+NriR6/Le/jLtdwk90UQYrdSNle+mML1vzHDcqlrzqnbk0lRLd7dC/Gz8oL4+tsRC3oXOEDohhBAyBsTskBp4dCPoDKETQgghUYRX4IQQQiIL70If57z/rwVPdtuvdY91x3F1+uDWHfmnx5/Qp3Xtrp+vtd9t1FOrHjm+kClu8gfHa32T+nWttOFIwZCfPKx70+NHJ2ttO11QEy2hl1uZrNaWfnMro2jMfvq5kx1RHyC0a4OWXqytedMNvnVJ1bR1IJCv3dxXviYexudu2XUlRpo/u5+2G0afrVZ4zuRrt2L1JfuCMlo+9nL1dWrnlaGWN3CG0AkhhJAIEokrcEIIIaQYtXwFHokNfKAvl3+e++0HWt+Jcz6itc+oL4TYczN+rfXVN+ohq18n+rT22/UDhfecrPcdbZquH+vwtPzzSQN62DV5NCnGJvLP4+lJWl8sJexpaT38bqXUULceBrekXW2oMNYS4XY4Oa2phdhleD1I+D2XHfHYSlrbJKbqcFZCjFVC0p6/e08YvCCP+FvZRp5y1orVif6Rh8mDhF5laL7clLNhkOu1Q1i6gqSVDROOH43wezXlj1qilu9Cj8QGTgghhBTDRshyohHewKmBE0IIIRGEV+CEEEIiS2gNPMTcsSYSG/hgqmDDSh3US3dOfksvx/mR4xryz48m9LFo3Ks16xt0PfA3scKxuhN6OdH+er19ZNLh/PP05GlaX3JQ18TrpxR07Xqhj9cd1TXKRFq/pSJ+tNC2M6L8aUo/lpVRdG5pT8sIXVizp8k+WU5UaO9Kelepgbs5XWu3NGubKLMaQC/3atVh9PPJSp/4rAY7F2Jm/dgylGH107xNeqhJPw9iXRstvOsduebtp/sG0c/DlGENop+PRfpX6uMFbHv4Ufb8yi1l1Iny2gkhhJCaJRJX4IQQQkgxGEInhBBCIgg38HHOsXRB+z32nq79TDqoa+BNswoa+IkN+sdLx/WxuUlvae1YsqDfTrb1sW/FdP32/UShP5XUx2brm7T2UENj4fmgnna17miD1pYaeWKS4iFP6YpHPK0bnGOpwue1RZ+VEdquqpH76eVDQkN2FC+79JtLX7iil1sJ6Wc26OUez3j5fnPLoJdbOWESF+jHTZTsAwBI/dygmUu9XD+u1MdL6+dhPOKVRPPWB/R5q2ueaHr5WJVHpUZeG0RiAyeEEEKKEbOHH2XPr9xSRh1u4IQQQiKLbQ0/wsyPKrwLnRBCCIkgkbgCTyn+52OHdd0x06N7vesOFTy+x03R9eWWBkdrZ2zd2436P+afxhO6Phu39Hadom++K7TOo3FdfxpKFHSwXJ2uiaUb9NzoQ5P0tqqRe/RxoXPHk4XvY3HhJ7czelvVy5ERWrRJLwd0jTwjtPacI8Yqaxb52GHIuW75+MvdjH5PgmUXPo8V89Gq1baUwOXYWGm93D8/e0J5btbLdQ1ZnP8Ker1V7XQ86OVA+X7nSurAJr08avr48PuO/c95tLBtC7EQl9HD/zWjmU81Ehs4IYQQUoxaDqFzAyeEEBJZeBPbOEeN8KZTIgz+gR4equsthFbrj9NDkdOn6j+qdEIP0w7FCuH4XF2vvghxpuJWYW7S0tdwyNZDhEeUMGymToTX6/S0q05cD7VmlXYmqYfjksf0dl2yELL1hNelBS1RaNsZvc/O6B/WSunnzVL7RQgdjgyhK/3CjmYlZKheLQkq0rdKy1lMX5MWYhcWMxkW18LzplKpcq5dOrxeDMsQQvdaxQxj5XHVEHsAy1klCZUK1pCutlLh9bBzVfzsaOWG2GV4fTRKmJKJRSQ2cEIIIaQYDKETQgghEaSWQ+i0kRFCCCERJBJX4FqVzCFdERwa0DWz+vcLOlLsA91GNmWq3p6eFJp4vKCHZizdnpZL/Ju+KGVqvE7XZ2PCcvaeoon3Cx01LXTUjLCkqbazbGqy3pfIiHbh82WFDS6RSIh24Ucfy+gxpHhab9sJmbJV+fAZqY8Lq5ha508cR6ZshV04liUtZ7LEqdDALZMGLu1qsdIpW1U7GmDWpj2WM4N2La1tQSxpJlucn/llLKKDQfTx4fFKKlVTOVcES9k6Fha0SlrOKqWJj5d0u9XCtizYYWxkLkAbGSGEEDLK1LIGzhA6IYQQEkF4BU4IISSyhL6JLZrRcwAR2cBzygkeyggN/Jiuleb6CxpZfEDXFesO6/ralMn6T32akoo0Fdd1VMc6rL9PXJnrcxuj6hmX+nifrbdTlu6jzqk+cDHWFW11bF1cT8maFZ7rnNJOpPVfA6kvx4RP3LVVD7nQy4UmrmrkRn0c0DVygz4OmDVyNyfGCl0bGUVblH2x0r5wkz4uxwLQNXFfb3fhWJXSx4u9jzbX0FdNouYhHyvPODXxkWHZ4TRwK8IbOEPohBBCSASJxBU4IYQQUoxavoktEhu4GkLPZPV4Rzath5yzg4UQYqJfD6nFBvTwW8NkPUQ1taEQek0n9PeR4eshq1DJLCsEmGys/NQAAyLEfkwJqedskUZWpF114gWbmRMT4fWECMfFCsdyRFgyF9fX70nLqoS+pQVNhtRjqiVNnCdLCldqJTMZXpdtkwVNWM4gzoUqEVgZET4U4XctxG4IrwOVs6B5wuuybTquYawcX0kLmhoW97NojfQ4xY5ltIYZKrx5xgY57jixnKkhdaZdLWDb3n8PgeZHOIQeiQ2cEEIIKYYdUgOP8gZODZwQQggJwOHDh7Fq1SrMmTMHDQ0NOP/887F79+58v+u6WL9+PWbNmoWGhgZ0dnZi3759FV8HN3BCCCGR5cMQephHUP7H//gfeOGFF/D3f//3eO211/CpT30KnZ2dePvttwEA99xzD+6//35s3boVu3btwuTJk7Fo0SKkUpWVPiIRQlel0WxOaNNDQhdWNPBkv64xxQd0LSveqOuDkwYLP8nGOv2nmhapSTNqiVCZdlVYnoJky4/ZpdOyHhV9GaGX55S2tJi5Qj93lP5sXNfHE2ldS3SEUdKJKWlYh6ReLixnsUJoS1rMZBlTK1NYo0cf97SlXS1TemxGatXK3Jiw0Im0q5aqe5v0ccCrkQ8VStv6BvhUrVrq2iL9ruso1ja/44ZA+6n7lC01HieAru1XDrVcy5lJD/dbkx+VspyNlcWs3PWPF0ZbAz927Bj+6Z/+Cd///vfxiU98AgBw66234tlnn8WWLVtw2223YdOmTbj55puxePFiAMAjjzyC5uZmPPPMM1i6dGn5i5Vrr9iRCCGEkIgyMDCgPdLp4l9mstkscrkc6uv1L1ENDQ148cUXsX//fvT09KCzszPf19TUhIULF2Lnzp0VXTM3cEIIIZHlw5vYwjwAoK2tDU1NTfnHxo0bi77f1KlT0dHRgdtuuw0HDx5ELpfDo48+ip07d+Kdd95BT08PAKC5uVmb19zcnO+rFJEIoRNCCCHFsK2QIfQ/KXfd3d1obGzMv55MlpZr/v7v/x7XXHMNPvKRjyAWi+Gss87CZz/7WezZs6f8hZRBJDZwzQcuMmhmpS88VRiQO6ZrkrEjGdHW+5OTCz+whnr9N2JKUn+fdLyg12YsPdSSS+hpV7NW4VhSH8/Gy/eMHxVt9dM5ss+jiU9WnuujXY9vWurnShpWW5QIFb9SMcXeEfNYPYSH3BAQ8mi9Jp+49IhLvXxI6RcpWeX7aGlZpc8buu7o8W/XFcq7+nnGPSVQTccdYV8lCeIZ9+rWI9e1PTq2z9xSx5HHCuIR91vTSOcFnUvGlsbGRm0DN3HyySdjx44dGBwcxMDAAGbNmoXPfOYz+OhHP4qWlhYAQG9vL2bNmpWf09vbizPPPLOia2YInRBCSGSpVAi9HCZPnoxZs2bhgw8+wPPPP4/Fixejvb0dLS0t6Orqyo8bGBjArl270NHRUYmPnCcSV+CEEEJIMayQd6FbZcx9/vnn4bouPvaxj+H3v/89vvrVr2LevHm4+uqrYVkWVq1ahdtvvx1z585Fe3s71q1bh9bWVlxxxRXlL7QIkdjATalUc6LtKNXKciJk7sgQ+lERYk8VwsH1R/Wf6mQZQldSrWZkilNbD5sNxQs2s5wlQujSciZ+IllXecHHuaOmXZVrckRb/eTSchYGGY5X07CqlrJhZJUzxXIm/iI935IzhtC3x0ZWOjyNnPmvV7WcuZDhUP0uVEvYyDy2M4UwVcLUuV5pQaROdUTo3nRcT2hbyiPqcWW62sIvZ9DrGVMaVr+QeqnjyGMFqWoG6J+vUhYzv7njxVZG/Onv78fatWvxxz/+EdOnT8eSJUtwxx13IPGnioo33ngjBgcHcd1116Gvrw8XXnghtm/f7rlzPSyR2MAJIYSQYli2BStMOdEyrsCvuuoqXHXVVaWPaVnYsGEDNmzYUPa6RgI3cEIIIZEldCKXCN8Jxg2cEEJIZAldzIQbeHVRU6k6QjzMCtkxo5QXzaWEHUrYynBMpC09WmjXCRuZtJVNrissJB3Xj5OK6zpwo1XQo3Ky9KiPJp6Nl/4RybKfJmTRSXWFfgq4LE2q6ulSW3fSpbVFx6fMakw5N1Ivl21b6unKXEt6DU34/fWqZUpln7TbGZxJXq26dBrWMGVKPXPF2xr18wrh0fdDpGE1vk8Ay1kl7V2jZRULo4mT2iASGzghhBBSjLHQwMcL3MAJIYREllrWwCO8dEIIIaR2idwVeE7kCJXlRXNZRQMXWqhfalXraEFbjE3RT01C6JsNqULIZkpC/x6UEalHM1ZWea7rWGlbL0WasnQdT9XEs0JDzrmlNWXZ57pijY6iN8u5JY8ajjppTTVaIv2+W4ZIw6pN9Hsf9XdEaLcZoe3WSW209M8niDZdKc+4HO+XhlXV0/10a22sTFsqBxs08SBpWIN4xoOUNAXKL0Xqp4ePVulO1RcexBMexVSwDKETQgghEcS2Qt6FXq27OUeBCH/3IIQQQmoXXoETQgiJLJYdLgzOEHqVUWVu6QPPOaVzoztDokzmkPAzS01c8Y3bwiOeqNdPVbKu8FOvq9fXMGlIj8mkYmqOcl1PGxI+8LR1RGtrGrj4cWUhPOOK7m3SxyXSIz5amng8U1rddcQfVSKAJm7Sw/WRgJ4VfvidSiPPlNCFDaVJA5XjNPT5HStMHnWjXi795cac5NLHPnJN3L/0aGldOwyVKkUaREMOqjervnB6wgvYMcubFyLIfFl7OUJE+LsHIYQQUrtE4gqcEEIIKQbvQo8Q0kbmGNo5EaKVqVXjQzLVqhJC94zVQ9KqrUy1lAHAUEJvT8oUFiXTrE62REg9dkxfsxImTxssZgCQU9Ku5sSPVobUHaXtOOI40HFFPNsUUrcN1ik/7NzI/wjNIfUwFrMwIXU5tXSJzSBh8CBWML9jjfS4st9kMQOChdQ9qHMD2MgkxjB4gOP4ESSkTqqPFbNghQihWxEOoUduAyeEEELyhLwCj7KQHOGlE0IIIbVLoA1848aNOOecczB16lTMnDkTV1xxBfbu3auNSaVSWLFiBWbMmIEpU6ZgyZIl6O3treiiCSGEEKBwF3qYR1QJFELfsWMHVqxYgXPOOQfZbBZ//dd/jU996lN44403MHnyZADA6tWr8cMf/hBPPfUUmpqasHLlSlx55ZX4+c9/XpUPIFOpZhUbmZpWFShmKxPih2orExq4ndbHxtMFrTdRr/8C1GX0dn280J4sdPmMpb9Pg7CZpZVSpE0GixkApNyCBuiXdlVtOz5lST1pWZXnUj/PJkunbbSFli7b6if3pF31HMvzimF0pTRxkx4OGDVx8bOqlCYe5DjFxpvmVivtqud9VQ1ZzjVq3n6Ws9JpS8PMDUK1bGVBSo2qaVWHjzvy1KpRgBr4CNm+fbvW/s53voOZM2diz549+MQnPoH+/n489NBDePzxx3HRRRcBALZt24ZTTjkFL730Es477zzPMdPpNNLpwi/fwMBAOZ+DEEIIqSlCaeD9/f0AgOnTpwMA9uzZg0wmg87OzvyYefPmYfbs2di5c2fRY2zcuBFNTU35R1tbW5glEUIIqSEsq2AlK+sR3Qh6+Ru44zhYtWoVLrjgApx22mkAgJ6eHtTV1WHatGna2ObmZvT09BQ9ztq1a9Hf359/dHd3l7skQgghNcaHIfQwj6hSto1sxYoVeP311/Hiiy+GWkAymUQymTSOUWVuIXl7fOCuklrV4wPPmFOpqqlW7Yx+YI8GPlT47pMQqVOTQ6V94Umhy9dLX7ij66hDsYI2l4buEZelR1O2ooEL3Tpt6+c4p2jkjsEjDnhLkeYU7dpxq6OnZXPm75Z2Qvar53GU0q46fp73AKVIledVTbvqSJ+7gVhB8w+TdlVixUrfS+DxVAutWn+f6pQeDTy3SqVHCfGjrA185cqVeO655/Czn/0MJ554Yv71lpYWDA0Noa+vT7sK7+3tRUtLS+jFEkIIISqhb2KL8BV4oBC667pYuXIlnn76afzkJz9Be3u71r9gwQIkEgl0dXXlX9u7dy8OHDiAjo6OyqyYEEII+RNWzA79iCqBrsBXrFiBxx9/HN///vcxderUvK7d1NSEhoYGNDU14dprr8WaNWswffp0NDY24vrrr0dHR0fRO9DLQYbMTf2uqFTmivi7I1OpKrYyV1Qjs2RIXYlEJtLCRmawlamWMgBIZ31sZVYhHDlk6yE3GVJXU61Ki1lahAFNlcs8tjERQneV8LW0kbnCkqb258RYS9jILCUNq+UJkevIudAqs+k/K0uE49WPIz6aJ6OTq7StRJjEhSKsbAi/ey1apZPX+obMZTjbLnwG38plfilQVQzhdu9xDdawQJW/xqZymYlKplkNWq2M1B6B/iNt2bIFAPAf/sN/0F7ftm0bPv/5zwMA7rvvPti2jSVLliCdTmPRokX41re+VZHFEkIIISq1HEIPtIG7rt/tNUB9fT02b96MzZs3l70oQgghZCRwAyeEEEIiiGWH3MDDFEIZYyK/gZtsZdmstJGJVKpismojg9DHTalV7XpZalTYyhQbmSw1Km1lqbi+pnq3ILYPubqeNtnSLVyqrUxazCbbR7W2mmo1SOlRQNe13YS+Bll61IkX2lK3znlSq8aU5z42MkO/lROfJ2EyPZn/eG1FW3eFDdGyxRpkW8Xv5g0TMT+7moLQy6VlK0jKVtM8j/bujFwvH++lR+WxTGlWi80dC4KkViUTh8hv4IQQQmqXsHeSWwG+H483uIETQgiJLLWsgUfXAEcIIYTUMJG4Ag+SStVR7pT3+sDF2IwsN6qkUhWlRmVqVdUXHsvq34NiIltlPGspz/W+hFh/vdDtU0rdzDrhEU9aul44RfGFp2WaVUvXyFKKXzXp6HrZkPCyZkSJU7X8qPSIO3X6sVy1bKnHBy7857mRa+BZQ7/tyG/UBs+1GOt62upY8Z6OVIbFn1NG+WGb9HEAyJXWyC1HP65Rxxa/fEFKhHqOZeirVOlRz3sa0qwCo1N6dLi/MhqyyRc+Vlp6kPKi41Hv9xCzhh9h5keUSGzghBBCSDEYQieEEEJIpIj8FXhOhDHVkLoMr+eyPjYyLVYvJg+VDqHbQ/pxYhkRUleOm8gJG1lGb6eljUz5EBlbD49KW1naKoQfk5beN0XYyFJuITSWlilaRdhsSKSDVCuZ5RxhU5JpV2OFcKkpzerwsZRUqoY0qwCQS+jnwlbSpebiuZJ9gB4WdxJ+377LrVwGQP0MfpXLTHfRirS4ekUuoQvF4iXHAnp42xvaDlCpLAgiLG4MzfukPA1W+asyqVSDhN+rGWJmJbPi8C50QgghJIrYITXwCCdyYQidEEIIiSC8AieEEBJZavkmtsht4NJG5ulXNHHHp5xoLiPLieaU52a9HKoGntH7bIONrE6kWU0IDTYh3kdtJ2x9TSZb2RRZahS6jqdq4ilX75O2sly8dLnRrLDJOLYsL6po4IY0q4Cue5vSrA639f6sool7+0R6USXVqmXr59tbTrT48+JjLdGvTPCzkZlSrZr08yBpVgFPqlX9WKX/HfiWCDWMrVSaVcAn1WoFbWVjQSQsW+ORmGW+h8R3vn+RLpVcLodbb70Vjz76KHp6etDa2orPf/7zuPnmm2FZw7/9ruvilltuwYMPPoi+vj5ccMEF2LJlC+bOnVv+OovAEDohhJDI8uEVeJhHEO6++25s2bIF3/zmN/Hb3/4Wd999N+655x584xvfyI+55557cP/992Pr1q3YtWsXJk+ejEWLFiGVKu25L4fIXYETQgghlWZgYEBrJ5NJJJNJz7hf/OIXWLx4MS699FIAwEknnYTvfve7ePnllwEMX31v2rQJN998MxYvXgwAeOSRR9Dc3IxnnnkGS5curdiaeQVOCCEkunyYiS3MA0BbWxuampryj40bNxZ9u/PPPx9dXV148803AQC/+tWv8OKLL+KSSy4BAOzfvx89PT3o7OzMz2lqasLChQuxc+fOin70yF+By2yWJh+4bJt84K7HBy50VDW1qlhELCvbhedx4QOXbekTTyjHTrj6mhIiN2wdCm+UhK4VSl94PQr6mvSIp239W2dWlBOtU3TvTEykWfWUHi38irlC83YNZUpNaVYBr86ttrMej7h+TnN1Sspc8T6mVKqmvuE1m1KtGtKsAmaNPEyaVUd/H0vRzH1TnpZeEeBX9lM9rm0uaar72oMct3SaUs/YgJq3yXNtOpZvmdIAPvYwqOVFJ3pp0fA+8OHfyO7ubjQ2NuZfL3b1DQA33XQTBgYGMG/ePMRiMeRyOdxxxx1YtmwZAKCnpwcA0NzcrM1rbm7O91WKyG/ghBBCSFgaGxu1DbwUTz75JB577DE8/vjjOPXUU/Hqq69i1apVaG1txfLly0dhpQW4gRNCCIkuo5zI5atf/SpuuummvJZ9+umn46233sLGjRuxfPlytLS0AAB6e3sxa9as/Lze3l6ceeaZ5a+zCNTACSGERJcKaeAj5ejRo7CF5BWLxeD8SaNtb29HS0sLurq68v0DAwPYtWsXOjo6wn9ehQl9Be4pJ+rjC9c08CGZSN0t2bZEqVE7q+uqMSXfuSXWkJCecWHTVS3MdWINdbbQwF1FAxelRush8qYr7TiEn1zmRnd0LahO0b1lnvSsKGNqK6Z4V2qhwsfuKp/HFZ/NiYm28MRnlTqtVs7HB676zRPCB54r7QuXPvAgvnBL2q+lZqdq8fJmDc9Y5WDCdw+Zz9weubfbdUaeC91U9tOjnQfQtb3vU9oX7quXq3N9Ne/x5wsn45PLL78cd9xxB2bPno1TTz0V//zP/4x7770X11xzDQDAsiysWrUKt99+O+bOnYv29nasW7cOra2tuOKKKyq6lgm9gRNCCJngjHI98G984xtYt24dvvzlL+PQoUNobW3FX/3VX2H9+vX5MTfeeCMGBwdx3XXXoa+vDxdeeCG2b9+O+vp6w5GDww2cEEJIdLGtcAVJAs6dOnUqNm3ahE2bNpUcY1kWNmzYgA0bNpS/rhEQiQ3cVOVTtqWtTOvzlB4tHWKXNjJX2siUELsMu8pF2EoqVTWtarG20UYmPmvCKR1Sj1n64KQIoceUmG69pYfM6y09ZO4pN6qUIq2z9eNKW5labtSJ6SFPk61MlhoNYitzY/pnd0RI2lXC74749p0TqW1txTombWOBbGWmX0xAt5H5pV2toK1MX0PpUqS+ySYrZCvzljitjq0sCOMxxel4XBMZfSKxgRNCCCFFidmjmgt9PMENnBBCSHQZZQ18PMENnBBCSHQZZR/4eGLCbeCOVk7UPFZkIhWpVEduOYO0kcnyokp3TLp8hI5qSq1qu8KC5mkX3ki1lAFmW1laplkVmnjKLl1udEjYxhJCE1fLjXpKjRpsZdJGFsRWplrKALOtTKY/tUU4zVHKjZpKjRZr6xY0UWrUZCuT5UNNtjJHHGiC28pM7xlIL6+gbUw9ljyOSav20+xHS+e2YoU7o91cZatlkeoy4TZwQgghNQRD6IQQQkgEGWUb2XiCqVQJIYSQCBL5K3Cz79vcNqZWNaROHW47pfs8PnBFRxWad0zMFRZmTT8XGUERE3ptTNPA9cFx4QtXNfBByLGiLfpVD3md8H0POTKVaq7oc6CIL1zVxO3S3u1i7Vy8cCxZatTkC5fHkb5wN6b6wM2+b1O/udSoQPrAg5QazYl8BTGDL9xP41b0cq8/e+T6eBiM2noQ3VqmZA2QWpVpVcc3bszS/k7LmR9VIr+BE0IIqWEYQieEEEJIlJhwV+Am65gMmXvmKuFITypVT6hSDYuLClbCVqamAbWl60fayAypVWXqVFNq1ZgIDddBD3nGlLB4UtrIIFKrCrtL2i30p0XIPCbC76qtTFYqs0RY31LC4K4IOTsi7aoMSattWalMti0tlar42cl0u8rbWLYYG8BWZqpUBhSxlanEDLYyU1WzYqihb4OlbHisYVGxkdvRJPKv0JiyNUwlM+X31i+tahBbGdOYji8YQieEEEKiiI1wseQIx6EjvHRCCCGkduEVOCGEkMjixsKFwV0f1Wk8E7kNXDq2TDiuj+ZtLCfqZyMrbTkzlRe1hcZtG2xjsh0TerlMrWorCmJCeJxioq2mVpU2Mb92zCroqNJyZrKVxWWp0ZywCKnH8jkxJluZI+1qYqwdc0r3Sa1a0ciDlhN1DaVIjbYyT588F0q/1Kk9qVRL28pck8YNwFK092qlQw2D1K0lQVKrVopKpmitFLbQ7J2JptlbIe9Ct6iBE0IIIaOOa1uem0KDzo8q1MAJIYSQCMIrcEIIIZHFtb0WzqDzo0rkN3BpzzbhSaVqKBEapJyoSfOW/VLatTxSr9S5S5cTFRZmrS017xj0dlxpq2lVgWCpVZO2rqelXV1vU33htjiOJ7WqcnL8yoma2h59XJYeVT3jnjSrwn+u9Fs5XV/284UHQQ3jWX6pVFXtWnrE5R+Ex0NumCsxpEv1lPLU3qOCaVbF+2ie8TAad4DUqpUtPVr4+5D+cb/yoiM9brFjT2RqeQOP8NIJIYSQ2iXyV+CEEEJqF2Zim0Co0WxTWtWgeFOpGqqRSVuZVpVKWL9kCF2GxR3DWJlKVQ3Vi9hKTIbflUXFRMhcplaNu3pItF4Jm3tC5jJlqxI2l2lWPalUlbYlw+DWyEPojgyZ21JOKN0ncZQfiC01i4z5D18NzckwnSfcbkylakiXKn/JPWMNNjO/imJqulQfy5n+Hvq/lSCpVT1Vz5zKWMHChKfHC6bwey3DEDohhBBCIsWEuwInhBBSO9TyFTg3cEIIIZGFG3iNIFOnmmxlfqVHTZYz7xurpUf1LkvYxkwaeczVx3p0bbd0X8IV6VAVG1ncYzHTx9YLTTyl6N5qWlXAazkz2shMmrisrxkgtaq0kZlSq9rSRmZMrarbrlyhiZtSq/qlXdVSq8rfPdPNHCaLGWBMrWqJkqCm1KqWsJy5lbSKVQnV6hbUcqZax/xsY9SmyVhRUxs4IYSQiYVjecsHBJ0fVbiBE0IIiSwMoRNCCCERxLFDXoFzA58YOIZyokad2+MRL51K1ZM61SP1ynKjSipVqY+LJampVm1RADIm24oIK/3ZvuVFlfHeUqOirXzAmPAky1Sqaltq0QGqyBbxiBtSq2b1PwGZWhWZWMk+v9Sqtq8DeoSYUqtK3TpIalXpEfekXTXo3Aavt+vnL68S47GUJyHVhBs4IYSQyMIQOiGEEBJB3JAh9Chv4BFeOiGEEFK7TLgrcNUym/Pzfft4vTV88p2rSC+36X285URLt2VfzCntC5dpuyVBfOAmPKVGrZHnQjf6wCWByovqJSglqr4utXaZg92JqfcvyPsKzJhyoZtyo7siT7q0xOuLkLp1BQsAqDp3EF075vNvJVQZ0MLP1pM3fZTKiwaBunz1cW0XrrwhKOD8qDLhNnBCCCG1Qy3fhR7hpRNCCCG1SySuwNVodYBIdij80qO6AcqJaogwrDeVqj5cs5GJVKq2W/p9YjLcK3J3auF2U+gaRWxmbulwqieVKkqPlcdV8ZQaDVBeVIbXZXnRWLYQhjWVGpX9Md9UqmKyT7nREeMpEaqmXfUJmZvKi3rSrBrC5DIsHqS8qEBNcSqRv9HVKi8aBobFxxe8AieEEEIiiGO7oR9BOOmkk2BZluexYsUKAEAqlcKKFSswY8YMTJkyBUuWLEFvb281Pjo3cEIIIdHlw1zoZT8CBsp2796Nd955J/944YUXAACf/vSnAQCrV6/Gs88+i6eeego7duzAwYMHceWVV1b6YwMIuYHfddddsCwLq1atyr82mt8+CCGEkNHkhBNOQEtLS/7x3HPP4eSTT8YnP/lJ9Pf346GHHsK9996Liy66CAsWLMC2bdvwi1/8Ai+99FLF11K2Br579258+9vfxhlnnKG9vnr1avzwhz/EU089haamJqxcuRJXXnklfv7zn4debKWR5UVdg6znypSUKj56uZZKNYRoL7N8+lTY1PvE28ZQenAsSCpVy2wvigewkaltj6XM6KXS8aRONZQXjQnLmetJpVr6OH5GMlfplvq4qbyoJUU5mfLURJDyoqJvopcXVbFsoWM7o69jq2VIAZYiLZdywuByPgAMDAxoryeTSSSTyWJT8gwNDeHRRx/FmjVrYFkW9uzZg0wmg87OzvyYefPmYfbs2di5cyfOO++8stdZjLKuwI8cOYJly5bhwQcfxHHHHZd/vZxvH+l0GgMDA9qDEEIIGQmhwufKDXBtbW1oamrKPzZu3Oj73s888wz6+vrw+c9/HgDQ09ODuro6TJs2TRvX3NyMnp6eCn/yMjfwFStW4NJLL9W+ZQDw/fZRjI0bN2onra2trZwlEUIIIWXT3d2N/v7+/GPt2rW+cx566CFccsklaG1tHYUVegkcQn/iiSfwy1/+Ert37/b0lfPtY+3atVizZk2+PTAwwE2cEELIiHCs4DeiyfkA0NjYiMbGxhHPe+utt/DjH/8Y3/ve9/KvtbS0YGhoCH19fdo+2Nvbi5aWlvIXWYJAG3h3dzduuOEGvPDCC6ivr6/IAkaiM1SKIKlTHZPm7UcAndtPxw4i7ai2ZI/mbfKMe1KpmlOrqm1vOVGD71vWTpXyrOw3EeQGAAOu8Igby4tmYqX7AMQyfslVq4Cf5h0GrfSoj8atpV0de6+2RHrPQ6VdJeOK4TB4GA28vHnbtm3DzJkzcemll+ZfW7BgARKJBLq6urBkyRIAwN69e3HgwAF0dHSUvcZSBNrA9+zZg0OHDuGss87Kv5bL5fCzn/0M3/zmN/H888+P6rcPQgghZLRxHAfbtm3D8uXLEY8XttGmpiZce+21WLNmDaZPn47GxkZcf/316OjoqPgNbEDADfziiy/Ga6+9pr129dVXY968efja176Gtra2Uf32QQghpLYZi3KiP/7xj3HgwAFcc801nr777rsPtm1jyZIlSKfTWLRoEb71rW+Vv0ADgTbwqVOn4rTTTtNemzx5MmbMmJF/fTS/ffhRycJMHoLYwRy1opXeJStcSdTxnmpkrifRpNJnXpIpDWslMaVLNfXJcLpMeSrTpZr6vO3Svxim1KqOOKkxEVWW/0QsJaxne5KCjhxPdTI1bO4n9chqZep4TyWzyoTfveHqaFnMJEFSp1bSGjYerG5RwLFcOFaIEHoZcz/1qU/BLfE/tL6+Hps3b8bmzZvLXtNIqXgu9NH89kEIIYTUKqE38J/+9KdaezS/fRBCCKltcvbwI8z8qBKJamSEEEJIMWq5Ghk3cAWTzcxUXlT2hSkiKcuLegsslsbkpLANx/GUGpXlQ002Mp+xpuN61hggXapE1bktlC5XKfErJ2qeKzRxw1h5o4ypbQnN23BKvXh07QrdCGLLcqIBdG1ZinSU0q6q2jVLfk5cxkIDHy9E+LsHIYQQUrvwCpwQQkhkoQZOCCGERBAHIVOpVmwlo8+E28Bzmo5t/qlWzScu9fIQJURVbPFbahs0e6/PW58bC6Cth8FUbtSkiXvKiVYQR0mfGsua9XKpkUeemHK5IX3ftvSMV0mr9ujplUlrWkvpUlmKlAATcAMnhBBSOzghQ+i8C50QQggZAypVjSyKcAMvgV/lMpOtzIg4rowUywyheti8OmFvv2pkJuLQw6zSVlYpZEjdkqlVlff1S6UKg80sjK3MlflrM6PwnyEmLh+qmj9YfV8ZBi//5y5D3yrGMLicV8mQuZo+lRY0Mk7hBk4IISSy5KzhR5j5UYUbOCGEkMhSy5nYIrx0QgghpHbhFXgl8CvpOErEDHq5X3lREzI9qmb/CnDcGMy2JL9Uq2ONVx83K+SaJl5JPVyzgo2P372qIXRu7Tfcmbg2MTJyGEInhBBCIgjvQieEEEIiSM6ykLPK34XDzB1rqIETQgghEaSmrsAnmlwYJsunLCE63ghTWjQIjvCIxwKUIg2CK+XyIPKtx+sd4NzYYm4IvzYh4xGG0AkhhJAI4oQMoTsMoRNCCCFkNOEVOCGEkMhCGxnxJUju87LzpFeQML7vSlKt3OhRR2ririLEydsTLJ+8/IRI1HKjE73U6LAGHiaEXsHFjDIMoRNCCCERhFfghBBCIkst+8C5gUeIMLax8ch4T53qh+v5gQQpPjr+sZSSoS7tZxXBUsqUuixTWhFysJBDiA08xNyxhiF0QgghJILwCpwQQkhkcSwr5E1s0b0C5wZOCCEksuRgIxcimBxm7ljDDZwQQkhkYSY2QgghhEQKXoETQgiJLA5sOCGuRcPMHWu4gRNCCIks9IGTSOBE94tiUXKeGpvRwppoPxCBm8uO9RImHPR+k0rCDZwQQkhk4V3ohBBCSAThBk58sWIj10nk2LGoJTVeSuRlIx4mrxYyi6ylVLCzxkE1OxJtJnoFMjIMN3BCCCGRJWvFkLXKv1DIWtH9whzd2AEhhJCax4GVD6OX83DKKGby9ttv43Of+xxmzJiBhoYGnH766XjllVfy/a7rYv369Zg1axYaGhrQ2dmJffv2VfJjA+AGTgghhIyYDz74ABdccAESiQR+9KMf4Y033sDf/u3f4rjjjsuPueeee3D//fdj69at2LVrFyZPnoxFixYhlUpVdC01FUK3J9jXlTAuppw1vk+GM0raue2MzvuEqpyaC1FH1plgNWgJEWRhIxviWjTo3LvvvhttbW3Ytm1b/rX29vb8c9d1sWnTJtx8881YvHgxAOCRRx5Bc3MznnnmGSxdurTstUrG939xQgghxEAOsdAPABgYGNAe6XTxGwF/8IMf4Oyzz8anP/1pzJw5E/Pnz8eDDz6Y79+/fz96enrQ2dmZf62pqQkLFy7Ezp07K/rZuYETQgiJLFnEQj8AoK2tDU1NTfnHxo0bi77fH/7wB2zZsgVz587F888/jy996Uv4yle+gocffhgA0NPTAwBobm7W5jU3N+f7KkVNhdAJIYSQYnR3d6OxsTHfTiaTRcc5joOzzz4bd955JwBg/vz5eP3117F161YsX758VNb6IdzAK0FsfAQycnZpO0QYX/iH31DzxypTn875/LqN99SqdsCbDqxqmfHDaOJRI5fRmq6TKTGQ1CqV0sAbGxu1DbwUs2bNwsc//nHttVNOOQX/9E//BABoaWkBAPT29mLWrFn5Mb29vTjzzDPLXmcxxsfOQwghhJSBE1L/dhDswuGCCy7A3r17tdfefPNNzJkzB8DwDW0tLS3o6urK9w8MDGDXrl3o6OgI/4EVeAVOCCGEjJDVq1fj/PPPx5133omrrroKL7/8Mh544AE88MADAADLsrBq1SrcfvvtmDt3Ltrb27Fu3Tq0trbiiiuuqOhauIGXwLLN4c8gqVU1xHFdEQORribHEBavFDIXcJBwVFb8ClUrdaorTpQrw9nK+1riJMq2iaBhcu19xiJ/7ViF0ytYqcxVwuRuLkCIPMjYoLBqWGTIujHEQvzfybrB/obOOeccPP3001i7di02bNiA9vZ2bNq0CcuWLcuPufHGGzE4OIjrrrsOfX19uPDCC7F9+3bU19eXvc5icAMnhBASWbKwEQsYBtfnB0/ScNlll+Gyyy4r2W9ZFjZs2IANGzaUva6RQA2cEEIIiSC8AieEEBJZsoiFvAIf3+4XExNuA4/5aNcqVUutKvXxcvVygdTDTXKtY5nfM1dGAv9yyLqlf8VMtjGpeVcSOzc6mvi4xKSZO2HyvQbAqZx+rhJIP484LBdaIKckYyl3flSZYP+dCCGEkNpgwl2BE0IIqR2ybhy2IdLnPz+6iZG4gRNCCIksWcRgUwMnJu+3yfddtie8CG4I37djWIapaL0sLSo1IfkLrral79v0x+CXKjVMCdEgXm+VMBq3bTrhAssJ0HZCeP9zVdKxw+jWFfSMB8Gll7smyLox2KF84NHdwKmBE0IIIRGEV+CEEEIiSxZ2yBB6dK9juYETQgiJLDk3brSr+s8fJftkFZjQG3jVfN5AMG+3oq1Le7Pr0d1d0V94LuXanFV++VA/n3ilMOnepj5HaNqe3OcGzdsvF7pl0L1Nmrjtc1JtR7Yrc44tqYk7Ae6aNWniVdLLJ5ofO4iWXkl/tutQwydmJvQGTgghZGKTdWOwavQmNm7ghBBCIksOcdghtrJcGcVMxgs1tYH7lQhVsWMh4u8BwuueyLBsB4jCqhFeOS9nCJn7lRM12shEX86UOtXH6iXD5ubBPiduhFg+aVXtnCmkPg5ufgkSTg9KkBB7ldKjVoqJFtYnBKixDZwQQsjEgiF0QgghJIJkXTvkBj4OImllEnjlb7/9Nj73uc9hxowZaGhowOmnn45XXnkl3++6LtavX49Zs2ahoaEBnZ2d2LdvX0UXTQghhNQ6ga7AP/jgA1xwwQX4j//xP+JHP/oRTjjhBOzbtw/HHXdcfsw999yD+++/Hw8//DDa29uxbt06LFq0CG+88Qbq6+sr/gHCYAtN3DJ8ibNMmriP5u0q/W6ItKtScg0iA3s0ccN3N99Uqq6aStX8K6SOlbYxmTpVbXvKiQb4hi1tYrJtmyxoBl3bZD8rfqziz4fbbum2qeSnH36auKFkqBsg5albrZSto8R4sGixJGhlyCEOK9RNbNENRAda+d133422tjZs27Yt/1p7e3v+ueu62LRpE26++WYsXrwYAPDII4+gubkZzzzzDJYuXeo5ZjqdRjpd+EUeGBgI/CEIIYTUJrmQGrhfjYbxTKBLih/84Ac4++yz8elPfxozZ87E/Pnz8eCDD+b79+/fj56eHnR2duZfa2pqwsKFC7Fz586ix9y4cSOampryj7a2tjI/CiGEkFoj68ZCP6JKoA38D3/4A7Zs2YK5c+fi+eefx5e+9CV85StfwcMPPwwA6OnpAQA0Nzdr85qbm/N9krVr16K/vz//6O7uLudzEEIIITVFoBC64zg4++yzceeddwIA5s+fj9dffx1bt27F8uXLy1pAMplEMpk0jlFlYykhV7CSp4ZfiVBNEw+yKKG7e1Or6m1HKS/qiNSpJo94TryPp2So4gs36eGAN8SUNfzayG+zJn3JFLqSGrgnlapbOl2qTJ1qG7zefuVETf2WSK0q2xVDauJBNHLTWIMe7j1O5Xze0pOttWWfM/782yxTOr7IuXFYoXKhR1cDD3QFPmvWLHz84x/XXjvllFNw4MABAEBLSwsAoLe3VxvT29ub7yOEEEIqBUPoI+SCCy7A3r17tdfefPNNzJkzB8DwDW0tLS3o6urK9w8MDGDXrl3o6OiowHIJIYQQAgQMoa9evRrnn38+7rzzTlx11VV4+eWX8cADD+CBBx4AAFiWhVWrVuH222/H3Llz8zay1tZWXHHFFdVYvwe1AllMhJFldbIgqVWDhMk9FcYM7+MNmZdue6qR2XpIXY3g+kVz1bC5X+pUE/Lbq7SVqWFyPxuZxzqmDfarMDbyNathcRki97SVExm0upjlFH9erA2l4pin+piJatq5yk2PKsPtlUyzqoTYK5oetUphcYbbq4/j2qHuJHcinMgl0AZ+zjnn4Omnn8batWuxYcMGtLe3Y9OmTVi2bFl+zI033ojBwUFcd9116Ovrw4UXXojt27ePOw84IYSQ6JNz40CNauCBV37ZZZfhsssuK9lvWRY2bNiADRs2hFoYIYQQQkoT3a8ehBBCap6cGwuUqbHo/IjCDVxBTa0qbWRGW5lMsyrGaqlUpcYtfnccoWtrNjKfUqOOYg1zIGxksq3YynxTp5pSqXrKiYq2Y9DAndJtj20sAF59XKRWNdnKApQPtYWOHVQjHzEyPaopXaqfJm7qD1E+dDyU66TeXJvU8gYeXfWeEEIIqWF4BU4IISSyOCFzoUs3TJTgBk4IISSy5Fw7ZAg9uoHomtrAZflQjy9c0ar9POLaWKl5e964tAbuCs3b5AvPWaV937KtpkoFgIyolRrEB55y60r2SwuGJ5Wq0jaVDwWED1z+QfrUTlV1br/yoarXW5YPlT5wUwlRbypV2a8+N5QPBfSUp0FSpbJ8aEnC6PJB9HSWBR1bHDcWKA9EsflRJbpfPQghhJAapqauwAkhhEwsnJB3oUf5CnzCbeBqNFuGyMNgeaxiI69GpkWGRWjeYxuTYXFDKlXZzijv6zmOJW1lpW1kaU/IXP81STmF6nF+1cdMqVQ9FceUdpDqY7Itq4+ZKor5VyNTpJKA1cZMqVQRKF2qrEZmsoKVH1L3HqvMFKg+FjNPaFtNj1ql6mOuE32LGUP1xXFc26s9Bp0fUaK7ckIIIaSGmXBX4IQQQmoH143BDREGDzN3rOEVOCGEkMjiOLHQjyDceuutsCxLe8ybNy/fn0qlsGLFCsyYMQNTpkzBkiVL0NvbW+mPDWACXIFLadqEyTYm236pVLW21MCFzq2mUvVIux5dW2riSipVg21MttVUqYCueQO6dSyF0jYxoEjJUKU/rejhQJFUqiYbmfjD0XTvgOVD1bZf6lStnKi0kYmTqvb7pU716NwB0EqIjlXq1ABlP706dgVLhmrH9dHPyz5uEJtY5fRzk44dRqenPj66nHrqqfjxj3+cb8fjha109erV+OEPf4innnoKTU1NWLlyJa688kr8/Oc/r/g6Ir+BE0IIqV1c1/bcEBt0flDi8ThaWlo8r/f39+Ohhx7C448/josuuggAsG3bNpxyyil46aWXcN5555W9zmIwhE4IISSyOI4d+gEAAwMD2iOdLh3V2LdvH1pbW/HRj34Uy5Ytw4EDBwAAe/bsQSaTQWdnZ37svHnzMHv2bOzcubPin50bOCGEkJqnra0NTU1N+cfGjRuLjlu4cCG+853vYPv27diyZQv279+Pf//v/z0OHz6Mnp4e1NXVYdq0adqc5uZm9PT0VHzNkQuhm6p6SmzLPFimVrUM5USNXm9D+dA/LST/1BFn3M/braVS9ejjwlOulAzNSA1ctNNI5J/7lg/1lAwtfAipjw/ldD1d1b2zos8TulKPFSB1qmybUqcCevpUv9Spqvc7SOrU4bZb9DmAIt7uAKlUTZq4j8+77PSpAfTx0UJq0yxpWhxnguvjwyH0MHehD//Nd3d3o7GxMf96MpksOv6SSy7JPz/jjDOwcOFCzJkzB08++SQaGhrKXkc58AqcEEJIdPkwkUuYB4DGxkbtUWoDl0ybNg1/9md/ht///vdoaWnB0NAQ+vr6tDG9vb1FNfOwcAMnhBASWVzHDv0Iw5EjR/Av//IvmDVrFhYsWIBEIoGurq58/969e3HgwAF0dHSE/ageIhdCJ4QQQsaK//k//ycuv/xyzJkzBwcPHsQtt9yCWCyGz372s2hqasK1116LNWvWYPr06WhsbMT111+Pjo6Oit+BDkzADdyU/9yvRKitmMpl7nNvLnS1RKjQwBMix3dpe7PH952N6+1MrNDOyNznnnZhHdL3PSRzlKtebpH7POXqoaOUU7ptyn0OABmncGxT7nMAcLMFXd4Sermd9WkrXm+pa5vaHt+3I9ul+6Tv29gWuc+tQLnQA+Q+93i7DXOlrh3Eyy3HKscKkvscEPnPK6hjB/FVs3xohHFi3n+sQecH4I9//CM++9nP4r333sMJJ5yACy+8EC+99BJOOOEEAMB9990H27axZMkSpNNpLFq0CN/61rfKX5+BCbeBE0IIqSFCFjMJOveJJ54w9tfX12Pz5s3YvHlz+WsaIdTACSGEkAgS+StwU1RchtM9qVSljcyUHtXTNpQTlSVD42q4XQ+diiyfEJk9NfdUxjNWlAhVrGJDlrSC6QdW06cGSZ0K6GFyaRuTIXQ1XaondWouAX1w6dqpJtsYAMS08LtIs+pJl6pIJTK8LsbqNjKIPtl2S7dNtjFAt4b5pVLV7F0BbWN+qVe19y0dFh8tvOH4Mm1aAeeNRzsYKcEoh9DHE5HfwAkhhNQwbsgNnNXICCGEEDKa8AqcEEJIdHEsb+bGoPMjyoTbwNX0qCZLGQAImbjscqLSNuYkhAaupkP1pFIVNrJYaRuZTJ2a8bRVDVx/IzV1KqBr4CkfG5ksGZpW+k22MQDIOYX3dUQeWU8CBeVYfuVDbUOJ0Hg2XrIPAOIZ1XIm06OWtpUFso0BmnXM1zampVL1KxFafipVzTom9HE3iD4uULVqP9tYpQijy1e2RGjpY4UpHzpadjU3lxqV96kWlhPz/H8IOj+qMIROCCGERJAJdwVOCCGkdrDcGKwQN6KFmTvWcAMnhBASXRw7pAYe3UB0JDZwzZ4tM5pKr7fRF24Z23o5UeEHrhPf0uoK/abyoYBeQlSmSjWlTgX0dKmm1KmA7v2WqVTTKO3X9qROFZq41MA1H7jQvGXJUNX77fF9Z/W2pbRt0SfbMTlX0cRNvm9AlB4VfbFM6ZKhgXzfgLlEqMnrbfJ9A7rOLfoq5fsePpaa4rR06tSgaKlTAU0jD6VrB0idGui44zB16nhc01hBDZwQQgghkSISV+CEEEJIMaxczJN9Mej8qDKhN3BPqlRT6lTRtupkrL50alWPjUyc1VyiEFqVefMzYmxW/C6p6VOHYqVD5sPtwsFMtjHZlqlSZchcplZV06dK25jjSaVaWJPJNgbooawgtjFAt46ZbGPD/YXzKEPmpopjQWxjw/0G65hfFTHTWO09A9jGAC0UXinbmKddyYpiwqIVKMQeqMIYU6dGFduxPX/zgYiwBh7dlRNCCCE1zIS+AieEEDKxsRw75E1s0b2O5QZOCCEksliuHWoTtsLUEh9jIr+BS9uYmj7Vt5yoSQOX/jRhI9N0b1nWMy7bBS1Upko1pU4FdKuYmip1uC20aZTWwNOGdKlHnEn6WKGBDwmbmWodk7axnLCKadYxg20M0K1ifraxWEbayJRUquLGAk/p0SHVclbaNibbgWxjsp2RNiyDjcxkGxP9QW1jmu5tso0BunWsSraxou874uNWT+M2pkAtM3UqUD2rm8ShzawmiPwGTgghpHaxnRhs1gMnhBBCosWwBh4ihB5hDTy6KyeEEEJqmMhdgUs7tqdfKydq9n3HEsKHrOjc0gcu56oauKd8qPR2Kxr4UJ2P5u3xgSsauBDxpQ9c1b2PoEHrkz5wVfc2lQsFgGNZ/Viq91stFwoUKRmqaNeW0MvtbOm21Lhl2+v1VnzgQtf2+sDV56V937Id2PcdJD2q2vbziBtSqZp8374Y0qUafd+Ang5Vat4BqGXfN9OjloeVi3nyRATBZSIXQgghZPQJm8jFk2AqQnADJ4QQEllqWQOP/AYe81QYK/4cAGJxc1jcMpU9EyF1LYReJ2xjCT2UqkZoZMg8nShtGwOAlPIhUpb+4zpmiaphSghd2sakVUytOOaxjYlQt0yXqlrHjLYxAFDaJtsYoIfJ43KstIIJq5jaHxP5aL1WMUVmERFab7vw87EyemjbysgwuMFG5gl1G8ZK25jjU3FMO460hpUOsYep/BUIn/B7kHB2Na1j5R5ntELfDLETSeQ3cEIIIbWLnbNDauC8AieEEEJGnVrWwKO7ckIIIaSGicQVuCZNG1KnAoBtKelQPTYyMVaUAbUNNjJZMlRtm1KnArqNLCvOeEasPyWOpVrHTLYxQLeOHXGFjUxYw1TdW9rGhhyzBp7J1uefOzlhGxvSj2UZ0qOa0qV6bWT6+8h0qTHFKqamSh3us0q2valTDelSTRo3YE6X6mcjM1jDjOlSPelQzbYyvexnaduYHGuyjQHCOhYiVapHHzdp3hW0jVVLXzatf6w0bTeXCjB2/OvutmPDDhEGj/IVeCQ2cEIIIaQYtXwXenRXTgghhNQwvAInhBASWeyc5cm+GAQ3xNyxJvIbuEkTj8dl6lQfH7haMlSUD0W93naSig9c6toiXWpaaaeFPp6WmrclfOCK91v6vgdRr7cV3Vtq3oPCBz6YK7RlqtRUTqROzervo6ZLdTN6nyldqkydakqXGkTzHu4vtG0hN0sNXPV6qz7vYm0oXm9jqtRi7Zwh5WlGasqqrm1IswroWrVPqtQgJUKN6VJNmndAqpUuNZCf3K/sZ4iSoWNBLZcP5V3ohBBCCIkUkb8CJ4QQUrvYjuUpShQEN8TcsYYbOCGEkMhi5cLZyBxmYhs9pO/b1O/1gQsttE6WEy1MthqE5i30c1X3ziR1bXQoUbqdEhp4SuRnT8X09z1mF3ThQei6tiwZesQt6Noy97mqeQPAMaUtNe8hoYnLfOdORlmH6DPlO5eat8x3HldM8lLzliVBZVvVuWPCXF+xEqHS5y3bJp07QIlQj+/bUCLUz/ftIUiJ0CAYPONjlet8LEqGBlm/77FqWNcOgpWzPPUOgs6PKtH96kEIIYTUMJG7AieEEEI+xHbDaeC2G90r8Ehs4EFSqaph81hC2sj0sKvdoH981UYmU6eqtjEAyColRDN1eng0LW1kcdVGVrpcKAAM2npY+SgK1qtBETIfFOlS1bC5tI0ddWQIvTA37Rcyz+qhe6jpUYd0G1lMtOPpQjsh0qwm0npbtY55QuRZv/SohfNY0RKhapjcLwzusYYp4zMitGoqEZoRoVNTiVAfK5ifzcyIGhaXtrExSpdqtneVb/0KNTdQ2J9h8UrDEDohhBBCAnPXXXfBsiysWrUq/1oqlcKKFSswY8YMTJkyBUuWLEFvb2/F35sbOCGEkMjyoY0szKNcdu/ejW9/+9s444wztNdXr16NZ599Fk899RR27NiBgwcP4sorrwz7UT0E2sBzuRzWrVuH9vZ2NDQ04OSTT8Ztt90G1y2EJF3Xxfr16zFr1iw0NDSgs7MT+/btq/jCCSGEECsX/gEAAwMD2iOdNssdR44cwbJly/Dggw/iuOOOy7/e39+Phx56CPfeey8uuugiLFiwANu2bcMvfvELvPTSSxX97IE08LvvvhtbtmzBww8/jFNPPRWvvPIKrr76ajQ1NeErX/kKAOCee+7B/fffj4cffhjt7e1Yt24dFi1ahDfeeAP19fU+71DGBxCiuJo+NSYsWrYoESpLhkLVxA2pUwEgq1jHMgbbGKBbxwaFtp6yhW1MpEs9qqRL7XenaH0eq5jSHsg16scVOnda0cBlqtScSHnqKRGqpE+VJUFVzRsIViK0LqWMlZq3p0So+Nlmij8fbgtrmKJ7GzVvQOjY4sAmzRvQdW6T5g3ourdJ8waMVrDxUCLUL1XqeCgROlrlRYMcJ8jYIKlTg5QPrWXa2tq09i233IJbb7215PgVK1bg0ksvRWdnJ26//fb863v27EEmk0FnZ2f+tXnz5mH27NnYuXMnzjvvvIqtOdAG/otf/AKLFy/GpZdeCgA46aST8N3vfhcvv/wygOGr702bNuHmm2/G4sWLAQCPPPIImpub8cwzz2Dp0qWeY6bTae2bzsDAQNkfhhBCSG0RNgz+4dzu7m40NhYufJLJZKkpeOKJJ/DLX/4Su3fv9vT19PSgrq4O06ZN015vbm5GT09P2essRqAQ+vnnn4+uri68+eabAIBf/epXePHFF3HJJZcAAPbv34+enh7tm0dTUxMWLlyInTt3Fj3mxo0b0dTUlH/Ib0GEEEJIKSoVQm9sbNQepTbw7u5u3HDDDXjssceqElUOQqAN/KabbsLSpUsxb948JBIJzJ8/H6tWrcKyZcsAIP/torm5WZtn+uaxdu1a9Pf35x/d3d3lfA5CCCGk6uzZsweHDh3CWWedhXg8jng8jh07duD+++9HPB5Hc3MzhoaG0NfXp83r7e1FS0tLRdcSKIT+5JNP4rHHHsPjjz+OU089Fa+++ipWrVqF1tZWLF++vKwFJJNJY6hCEhNfOaQPXG17fOD1MnWqaCvpU3NirOr7BvSSocfqS/u+AeCoso60x/et680Dwuvdh4LuraZKBYA+Z6rW7s8W2kFKhHo0b1kiNFPa6y193zGhiate70Rafx9PiVBF95aadyJdWvOWbZPmPdxWNGWT5g3oOncQzRvQvN9GzRsIlB5VL/M5cs3bO9enRKhJLxeomnLglKyjUCK0kmlVK5kulYTHcoqkRA44PwgXX3wxXnvtNe21q6++GvPmzcPXvvY1tLW1IZFIoKurC0uWLAEA7N27FwcOHEBHR0f5Cy1CoA38q1/9av4qHABOP/10vPXWW9i4cSOWL1+e/3bR29uLWbNm5ef19vbizDPPrNyqCSGEEAzXK7Byrv9Aw/wgTJ06Faeddpr22uTJkzFjxoz869deey3WrFmD6dOno7GxEddffz06OjoqegMbEDCEfvToUdjiCjIWi8Fxhr/CtLe3o6WlBV1dXfn+gYEB7Nq1q+LfPAghhJAPr8DDPCrNfffdh8suuwxLlizBJz7xCbS0tOB73/texd8n0BX45ZdfjjvuuAOzZ8/Gqaeein/+53/Gvffei2uuuQYA8tlobr/9dsydOzdvI2ttbcUVV1xR9iJVp5goMIaYeCGm2Mi8tjFzKlVXCZs7ohpZpkGEyZOFn7q0jR0VqVSPxgvHkiHzw5Yegu6HbhXrdwrtvpweMh/ICquYUmFsMKsfJ52ZrLUzQ4WQuidkntZD9Z4wudJOCNuYTI+qhs1VmxjgtYolUkoI3WATK96uUHpUk1XMM3YcpEcNEjIHKpYe1c8qpvXJkHOIkHm1KowZw+9VSpUa1KoWxDpGRpef/vSnWru+vh6bN2/G5s2bq/q+gTbwb3zjG1i3bh2+/OUv49ChQ2htbcVf/dVfYf369fkxN954IwYHB3Hdddehr68PF154IbZv3z7md+sRQgiZeFi5kCH0EHPHmkAb+NSpU7Fp0yZs2rSp5BjLsrBhwwZs2LAh7NoIIYQQM447/AgzP6IwFzohhBASQSJRTlRF2shkKlU1faosHxqT5UOniLKZkwr9OYNtDNCtY0eSujZ6WGjvg7HCcY/YukY8AF2blulS+5yCzt0v0qMezuqauKp7DwkbWVaUBFV17yCaN6Dr3ibNGzCnR1U1b6D8kqBAwPSomgZewZKgQ0KjVPXpKqVHDaJ5A9VLj6q/x8g172LHHnlf9EqEjlY50XLTp0ax3ClD6IQQQkgUYQidEEIIIVGCV+CEEEKiS87xWjyDzo8okdjAVd1b+sDj4hMklLKfntSpQgOH8HrnJhXaQ7qEjGMN+g95sL7QPipSth4VQv1ArKAZfyB83u+5TVpbpkd9L1uoM3s4V1rzBoBjmUI7k9E/gCN0blX39mjcx3RdPjEkdW4lPWpKpkeVKVDVcqKi9GvKVBLUFX0izCV0bk33dqRfW2rgqg7sp4EX+l2Pl9uc8lRLjyrmVszb7ZfiNEB6VLNmXLmSoF49vXytuhrH8X+f0dGJ6fseIY7j/ZsPOj+iMIROCCGERJBIXIETQgghRXFChtAjfAUeiQ1cdYol4qVTpwKArYSzY8ImZou2O0n/+Ln6QkAiNUmEzEUI/UiyENI9LOxqAzE95PyBYhWTIfN3nWlaWw2ZA3rYfGBIn3t0SLeVaelRZchcpEuNpwprist0qIaQOaCHzdUQ+XCfqDCmhM1lCD0mIoSB0qF60pqqaUtF2Fj+gQ4p/YYKYoBIh+oJt6fE2NG3hvmlQ5W4WgjdL9Stvo+PrUqrRhbGNla+NSzIsaplEwszN0zIvFzb2ETAchxYITbhMHPHmkhs4IQQQkhRavgmNmrghBBCSAThFTghhJDoUsN3oUdiA9c0cF1uRlxo4HHFOuaXOjU7Re9PTy5osNI2dqS+dLrUvriuEb9nCSsYCtq1n+b9gWirurdJ8wZ03dtO61YwVfMGgLhqBfOUBBXpUEU7kS6ct9iQtI2ZrGEQfRVKhwqUbQ2rVjrU4WMdK91nTFtanXSow/2qVu1nQVN04iqlQx3uN2jVY6B5+x2rGvNICHKO928+6PyIwhA6IYQQEkEicQVOCCGEFIUhdEIIISSC5HIhQ+gh5o4xkdjA1cyksnxoQpTuTEwu6Nx2kyih2Sg08Mm6fnt0cuGb2IDwgfcn9fftixeO9Z6t68u90HXs3tyMwvPs8Vrfu0N6eyCje72PpQu6dy5TWvMGgJjq7U7pfVLnTh4rHCtIOtThYyvpaoW3W+rcsXSZ3m6Txl2sHaTsp6ple9KjGnRuPx17DLzdQdKhyvFBvN3e44w8Hap3bvmad6DUquPQ660yXnzf1O2jSyQ2cEIIIaQYrpPz3lAacH5U4QZOCCEkujg5b0bFoPMjCu9CJ4QQQiJIJK7ATbnQ40n9O0hc08B1bTcrNPBjU3Uf8uHJhW9iH0zS+/rr9FP1bqxQuvOQQfMGdN1bat59Q9O1diqte8hzitfbOqb3JTze7oLOXXdMauClvd2qr3v4uLoGHsvo5zimyG9Byn76ervVu0Glxj0k2oYc5q7Uog05zH293FrpztIat2esZ67Byy3x0cvHquznyPtGJ5+5Z2yVNO9KasTjRfeeUPAmNkIIISSCOE7IEDptZIQQQsjok8t6o2NB50eUSGzgagi9LqGH0OsaRLhXCZtL29jQVH3skan6t7Y+xUb2flI/Nb1xPXx9EIUw+du5mVrf25kWrf3eUGFsf1oPmaeHpmht55jeVlOiynSoMkxusoaZyn5KK1g87WMNS6mhbv0cGsPk8puuKUwepMwnINKjypSmpdOjVssKBpQfJh8rK1iQ/kqlP/U7lmfsBLaGeddQnZA5bWMTh0hs4IQQQkgx3FxIGxk1cEIIIWQMcLKeqFvg+RGFNjJCCCEkgkTiClwtIZqoExp4o671xqYXrFS54/S+I8fpGuz7jXro5N8mFb7PHEroenM3TtDab+Vm5Z8fHNI170NpXRM/ki7YzDJpXeOGtIYdFW1F5/amQ9XbdamEMrZ0+lNA17l9y3ym9G+oms5tsoIBus5dJSvY8PgyU55WywomxkfNCjbcX/0yn56x41zjLkaldG9q3mXihLyJLcJX4JHYwAkhhJBiuLkcXLv8TTjKGjhD6IQQQkgE4RU4IYSQ6FLDN7FFYgNXNfBkvSgfepxeMhTNBS906jjd9/3+NF1bPKTL3OhRdO9/RbPW9y/ZNq3dPfSR/PN3U7o+fjQ1TWvnjhZKgsaONWp9suxnUni7644WPk/9Uf2zmnTuqnm5AV3XNpX1BHT/tizrafBnG8t6AnAzx8T7lqlrS0w6d5DjQNeJfccG0LnD6Nr62ImV8rSSc1Uq6e1WYWrUCpHLAiFC6FFO5MIQOiGEEBJBInEFTgghhBTDzWXgyhBjwPlRJRIbeL2SPrVhqkgJ2qLHwbMzC9aqD47Xw649jbo96kC9bunaa52Yf/4vGT1k/lZqttZ+Xwmbp47p1cgw2KQ165QQev3gFNHXoLVNYXKTFQwAYqnC57PTAcLiphA5YLZ/+aU4zSjhR5+weLkpTuX4QClOfcPgpVOcVjIsrh+3OilOR3LsIMfSxgYIi4d5n0rMK8ZEC5NPeOuYgpvLhtzAoxtCj8QGTgghhBSlhm9iowZOCCGEjJAtW7bgjDPOQGNjIxobG9HR0YEf/ehH+f5UKoUVK1ZgxowZmDJlCpYsWYLe3t6qrIUbOCGEkMjiOplhHbzch19GRcGJJ56Iu+66C3v27MErr7yCiy66CIsXL8ZvfvMbAMDq1avx7LPP4qmnnsKOHTtw8OBBXHnlldX46NEIoTckFQ18hq4RW6267apfcX/9cYYeGvnXSbpe/gbmaO3fpP9d/nn3MV0D7z+qp0fNDRZ078SRaVpfclBPh1p/pPC+UuOuOyrLfgqd+2jhuTfFqf6LZ0xxKnVtkxWsWrr20LGSfcNzDTq2pEJ2ryC6dtDSnaNh9xovpTvLfY9KzpVUS9eWjIXOXUsaty+5DGCFuBEt4E1sl19+uda+4447sGXLFrz00ks48cQT8dBDD+Hxxx/HRRddBADYtm0bTjnlFLz00ks477zzyl9nEXgFTgghpOYZGBjQHum0/5ekXC6HJ554AoODg+jo6MCePXuQyWTQ2dmZHzNv3jzMnj0bO3furPiauYETQgiJLG4uGy6E/qcoYltbG5qamvKPjRs3lnzP1157DVOmTEEymcQXv/hFPP300/j4xz+Onp4e1NXVYdq0adr45uZm9PT0VPyzRyKETgghhBTDdYbgOjH/gYb5ANDd3Y3GxoLlN5lMlpqCj33sY3j11VfR39+Pf/zHf8Ty5cuxY8eOstdQLpHYwCcr6VPrW3UdOz1H18Dfbi7odnsb9bG78TGt/etj8/S5gwVN/NhhXfOOH56utaccnpZ/PmlA17yTQueuP6x4uQ3ebaCIf1vVuaWuLUM8ql9b6NjGNKZSizalOIUssTk6/uyxSlta6j2C9o+VP9szd5TSmFbjOMDE1rW9a6DOPZp8eFf5SKirq8O/+3fD90wtWLAAu3fvxte//nV85jOfwdDQEPr6+rSr8N7eXrS0tJQ4WvkwhE4IISS65DLhHyFxHAfpdBoLFixAIpFAV1dXvm/v3r04cOAAOjo6Qr+PJBJX4IQQQkgx3FwGrhUihB5wA1+7di0uueQSzJ49G4cPH8bjjz+On/70p3j++efR1NSEa6+9FmvWrMH06dPR2NiI66+/Hh0dHRW/Ax2IyAbeOK3ww4mdoqct/deT9JDtr6cXUpP+P/cMrW/34fla+90B3Spm9xVCHI0fHK/1TerXw+QNRwopW5OH9V+e+FEZFi9U6LKOirCYDIunROguUzosbgp1BwqDZ1Kl+4q01W+sHg+l0bIlj2sKOcvwejAL18jfxy8sHsB2NYHSlIada6KWwuAmGCKPJocOHcJ//+//He+88w6amppwxhln4Pnnn8d/+k//CQBw3333wbZtLFmyBOl0GosWLcK3vvWtqqwlEhs4IYQQUgzXGYKbK18NDvoF+aGHHjL219fXY/Pmzdi8eXPZaxop3MAJIYREFjeXgYvRC6GPJ7iBE0IIiSxubghuiPuxg2YzHE9EYgOfflKhBOcHZ+jlN3e26h/hH3OfzD9/9f1ztL7su3pJ0MZDrXr7vUIZ0En9IsXpYaFrHy1oytZRkSJU2rtUXdtPx5Z6tJK21PNN0aBdV9KyVTGtOpQNq3ytejxatCr1vtXUUUdLq1YZ77p1Mahlk7EiEhs4IYQQUgzXGYJrWf4DDfOjCjdwQgghkWU4hB5iA49wCJ2JXAghhJAIEokr8MmfKvi1nz5VX/LXB/+r1n7nQMEsP737JK3vuHf0NHkN74rynEcU/e1wn74I4c92U4U6n05qUO/LCE1c0ao9mvHQUb1t8FX7adOaBl5B33SYFKEjPU6Y4wJjU/qymsdSGQstGoimHl0K6tQTl+G70MNcgfMudEIIIWTUcZ00XLj+A0vOZwidEEIIIaMIr8AJIYREl1warlv+FTgifAUeiQ380asKPvCv7ftbre/kXWdp7dP3Kz+Md9/T+ty+/VrbGezT+xXt2kn1631SqzaUr3QzR1AKP305jJ+5Wvm1x2NZSclY6cQqE0kzBqgbk2jg5oZCbeCe+44iBEPohBBCSASJxBU4IYQQUgw3l4brOv4DS82P8BV4JDbw+2YWfjgtWK31DYqxsk0IIWTiwg2cEEIIiSDDGniYDTzrP2icQg2cEEIIiSDj7gr8w7sJ08pdhTk3ut+QCCGk1vjwf3Yoe9cIyeQG4Tjl1wPPubkKrmZ0GXcb+OHDhwEA33QAfJhd54NdY7YeQggh5XH48GE0NTX5DyyDuro6tLS04I2eV0Ifq6WlBXV1df4DxxmWOxpfkQLgOA4OHjwI13Uxe/ZsdHd3o7Gx0X9ijTIwMIC2tjaeJx94nkYGz9PI4Hky47ouDh8+jNbWVth29ZTaVCqFoaHwiVjq6upQX19fgRWNLuPuCty2bZx44okYGBgAADQ2NvIPZATwPI0MnqeRwfM0MnieSlOtK2+V+vr6SG68lYI3sRFCCCERhBs4IYQQEkHG7QaeTCZxyy23IJlMjvVSxjU8TyOD52lk8DyNDJ4nMh4YdzexEUIIIcSfcXsFTgghhJDScAMnhBBCIgg3cEIIISSCcAMnhBBCIgg3cEIIISSCjNsNfPPmzTjppJNQX1+PhQsX4uWXXx7rJY0ZGzduxDnnnIOpU6di5syZuOKKK7B3715tTCqVwooVKzBjxgxMmTIFS5YsQW9v7xiteHxw1113wbIsrFq1Kv8az9Mwb7/9Nj73uc9hxowZaGhowOmnn45XXinklHZdF+vXr8esWbPQ0NCAzs5O7Nu3bwxXPPrkcjmsW7cO7e3taGhowMknn4zbbrtNK9DB80TGFHcc8sQTT7h1dXXu3/3d37m/+c1v3C984QvutGnT3N7e3rFe2piwaNEid9u2be7rr7/uvvrqq+5f/MVfuLNnz3aPHDmSH/PFL37RbWtrc7u6utxXXnnFPe+889zzzz9/DFc9trz88svuSSed5J5xxhnuDTfckH+d58l133//fXfOnDnu5z//eXfXrl3uH/7wB/f55593f//73+fH3HXXXW5TU5P7zDPPuL/61a/c//Jf/ovb3t7uHjt2bAxXPrrccccd7owZM9znnnvO3b9/v/vUU0+5U6ZMcb/+9a/nx/A8kbFkXG7g5557rrtixYp8O5fLua2tre7GjRvHcFXjh0OHDrkA3B07driu67p9fX1uIpFwn3rqqfyY3/72ty4Ad+fOnWO1zDHj8OHD7ty5c90XXnjB/eQnP5nfwHmehvna177mXnjhhSX7HcdxW1pa3L/5m7/Jv9bX1+cmk0n3u9/97mgscVxw6aWXutdcc4322pVXXukuW7bMdV2eJzL2jLsQ+tDQEPbs2YPOzs78a7Zto7OzEzt37hzDlY0f+vv7AQDTp08HAOzZsweZTEY7Z/PmzcPs2bNr8pytWLECl156qXY+AJ6nD/nBD36As88+G5/+9Kcxc+ZMzJ8/Hw8++GC+f//+/ejp6dHOU1NTExYuXFhT5+n8889HV1cX3nzzTQDAr371K7z44ou45JJLAPA8kbFn3FUje/fdd5HL5dDc3Ky93tzcjN/97ndjtKrxg+M4WLVqFS644AKcdtppAICenh7U1dVh2rRp2tjm5mb09PSMwSrHjieeeAK//OUvsXv3bk8fz9Mwf/jDH7BlyxasWbMGf/3Xf43du3fjK1/5Curq6rB8+fL8uSj2N1hL5+mmm27CwMAA5s2bh1gshlwuhzvuuAPLli0DAJ4nMuaMuw2cmFmxYgVef/11vPjii2O9lHFHd3c3brjhBrzwwgs1XWLQD8dxcPbZZ+POO+8EAMyfPx+vv/46tm7diuXLl4/x6sYPTz75JB577DE8/vjjOPXUU/Hqq69i1apVaG1t5Xki44JxF0I//vjjEYvFPHcG9/b2oqWlZYxWNT5YuXIlnnvuOfzf//t/ceKJJ+Zfb2lpwdDQEPr6+rTxtXbO9uzZg0OHDuGss85CPB5HPB7Hjh07cP/99yMej6O5uZnnCcCsWbPw8Y9/XHvtlFNOwYEDBwAgfy5q/W/wq1/9Km666SYsXboUp59+Ov7bf/tvWL16NTZu3AiA54mMPeNuA6+rq8OCBQvQ1dWVf81xHHR1daGjo2MMVzZ2uK6LlStX4umnn8ZPfvITtLe3a/0LFixAIpHQztnevXtx4MCBmjpnF198MV577TW8+uqr+cfZZ5+NZcuW5Z/zPAEXXHCBx4b45ptvYs6cOQCA9vZ2tLS0aOdpYGAAu3btqqnzdPToUdi2/i8yFovBcRwAPE9kHDDWd9EV44knnnCTyaT7ne98x33jjTfc6667zp02bZrb09Mz1ksbE770pS+5TU1N7k9/+lP3nXfeyT+OHj2aH/PFL37RnT17tvuTn/zEfeWVV9yOjg63o6NjDFc9PlDvQnddnifXHbbYxeNx94477nD37dvnPvbYY+6kSZPcRx99ND/mrrvucqdNm+Z+//vfd3/961+7ixcvrjl71PLly92PfOQjeRvZ9773Pff44493b7zxxvwYnicylozLDdx1Xfcb3/iGO3v2bLeurs4999xz3ZdeemmslzRmACj62LZtW37MsWPH3C9/+cvucccd506aNMn9y7/8S/edd94Zu0WPE+QGzvM0zLPPPuuedtppbjKZdOfNm+c+8MADWr/jOO66devc5uZmN5lMuhdffLG7d+/eMVrt2DAwMODecMMN7uzZs936+nr3ox/9qPu//tf/ctPpdH4MzxMZS1gPnBBCCIkg404DJ4QQQog/3MAJIYSQCMINnBBCCIkg3MAJIYSQCMINnBBCCIkg3MAJIYSQCMINnBBCCIkg3MAJIYSQCMINnBBCCIkg3MAJIYSQCMINnBBCCIkg/z+gmrTr2zWdvgAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "# Initialize the grid with a 100 degree wall on the left side\n", "nx = 100 # Number of grid points in x direction\n", "ny = 100 # Number of grid points in y direction\n", "dx = 1.0 # Spacing between grid points in x direction\n", "dy = 1 # Spacing between grid points in y direction\n", "dt = 0.1 # Time step\n", "k = 1 # Thermal conductivity\n", "\n", "grid = np.zeros((nx, ny))+25\n", "grid[:, 0] = 100.0\n", "\n", "# Time-stepping loop\n", "for t in range(10000):\n", " grid[1:-1, 1:-1] = grid[1:-1, 1:-1] + dt * k * (\n", " (grid[2:, 1:-1] - 2 * grid[1:-1, 1:-1] + grid[:-2, 1:-1]) / dx**2 +\n", " (grid[1:-1, 2:] - 2 * grid[1:-1, 1:-1] + grid[1:-1, :-2]) / dy**2)\n", "\n", "# Plot the final temperature distribution\n", "plt.imshow(grid, cmap='turbo', interpolation='nearest')\n", "plt.colorbar()\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": null, "id": "4086c295-a895-40e1-add9-36600ce6559b", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python", "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.8.10" } }, "nbformat": 4, "nbformat_minor": 5 }