{ "cells": [ { "cell_type": "code", "execution_count": 2, "id": "6939d4fc", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "82.65368591439173 44.342458599117634 21.767499999999995\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAf0AAAFzCAYAAAA0dtAgAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABaiklEQVR4nO3dd3zV1eH/8de5Nzd7swkZoICylSXiQMFRv+6tuC1Uq3W0Wq102EFrrdvagdXW/hpHWxeodaE4gcg0gCyBBEIgELLvTXLH+f1xLyEhCQRIcjPeTx553JvP53M/n3MvkHfO+ZxhrLWIiIhI1+cIdwFERESkfSj0RUREugmFvoiISDeh0BcREekmFPoiIiLdhEJfRESkm4gIdwHaWs+ePW1WVla4iyHSraxbtw6AoUOHhrkkIt3P0qVLd1trezW1r8uHflZWFkuWLAl3MUS6lSlTpgCwYMGCsJZDpDsyxuQ1t0/N+yIiIt2EQl9ERKSbUOiLiIh0Ewp9ERGRbkKhLyIi0k0o9EVERLoJhb6IiEg3odAXERHpJhT6IiIi3YRCX0REpJtQ6IuIiIRBdm42WU9k4filg6wnssjOzW7za3b5ufdFREQ6muzcbGbMnYHH5wEgryyPmfNmAjB95PQ2u65CX0REpA0FbIAtpVtYs2sNy7YtY9G3i3j/2/fxR/rBCfgBJ7i9bmbNn6XQFxER6egCNkBeaR5rdq0hd2cuC79dyKqdq9havRUv3rrjerh64I/w77vB7tx3jvyy/DYto0JfRERkP9m52cyaP4v8snwykjKYPXV2XQ08YAPkl+Wzumg1q3et5qstX7Fy+0ry3HnUUlt3jhRnCkN7DOWcEecwqu8ohvcezrBew0iOTibriSzyyhqvgJuRlNGm70uhLyIiUk92bjYz583E7XUDwfvtN75xI88sfIaSyhI2VWxqEO5JjiSGJA9h6pCpHDfgOIb3CoZ7SkxKs9eYPXV2g2sAxLpimT11dtu9MRT6IiLSzVlr2VG5g1VFq1i9azU/nf9T3G53sNk91PTuDXhZlL+I0T1Hc+2IaxmXOY4RvUcwrNcwUmNSD/mae1sNmmtNaCsKfRER6TaK3cWsKlpV97V061K+Kf6GSn/lvoNqAS/BDnZR7Lv37oLldy1vtbJMHzm9zUN+fwp9ERHpdA50zx2gvKac1UWr62rvy7YtY1XRKkq8JXXHRNkoMqIzOCPtDE4afBLHpR3H8N7DmfDshLDcb28PCn0REelUmrrnftObN/Fi7ov4fD5WFq5kZ/XOuuNd1kW/iH5M7DWRSUdNYmLWRIb3Hk5aQhrGmEbnD9f99vag0BcRkU5hR+UOlhcu5/a3b8dd44YAwXvuBmr9tbyz7h36OftxbMqxXDfsOk4afBIj+4wkMzkTh2n5BLThut/eHhT6IiLSoQRsgI17NrK8cDnLdyxn0ZZFrNy5klJf6b6Dagjee3cCMaFHB2z/+fZWKUM47re3B4W+iIi0uoPdc9/L4/WwqmgVK3asYHHeYhblLWJD+Ya6IXEO6yDNlcaJvU/klMGnMGngJK557Rq2sjUY9vVkJmW2x1vr1BT6IiLSqpq65z5z3kwqaioYnDqYr7Z9xafrP+Xroq/ZXrsdaywA0UQzOGEw1464lsmDJnNcv+MY1msYkc7IBuf/3bTfddl77m1NoS8iIq1q1vxZwUAOANWAH9wuN7e+fWvdMT0iejCq9yhuGHQDx/c/njF9xzAweWCTHev215Xvubc1hb6IiByxrdu38p+F/+Hd9e+StysPEgiOb48BfNRNcvPBtR8wpu8Yesb2PKLrddV77m1NoS8iIi1WWVnJ6tWr+frrr1mSv4TNZjObzWbyHfnUUovDOIhMjaTWH5qm1gCu4NPMpEymDZoWtrKLQl9ERJrg9XrZsGEDubm5rFu3Dq/XSznlFMYUsithF6s9q9kVsQuAIT2GMGPQDKYNmsaUrCm8veFt3XPvoBT6IiLdmLWWbdu2kZuby6pVq6isDE5HGxERwYCjBuDu5aZobBEfbfmIb3Z/A17oVd2LaUdNY9qg4Nf+M9XpnnvHpdAXEekmSkpKWLVqFc+9+xxvfPUGZTVlJEUncekJl3Lz2Tdz84ybWVuxlg82fcD7m95n8abF+Db6iImI4ZTMU7j5uJuZNmgaI/uMPOhkN7rn3jEp9EVEupjq6mrWrl1Lbm4u3377LYFAAIDk5GSK4ov4t//feE70gIEyyvin458s+3oZGxdspKK2AoNhXP9x/PjEHzNt0DROTD+RqIioML8raQ0KfRGRTioQCLBp0yZyc3P55ptvqK6uBiAqKopjjz2WE044gauvvhqn01n3mqwnsvDEeBqcxxvwkluUy3eP+y7TBk3jtIGnHdZysdLxKfRFRDqBnTt3kpubS25uLiUlwZXiHA4HAwcOZOTIkZx99tnExMQc8BwrdqxocvU4AH/Az5/P/XOrl1s6FoW+iEgHsndIXG5uLlu3bgWCne169+7NyJEjue666+jRo0eLzxewAd5e/zaPLXqMBVsWYDBYbKPjusKysXJwCn0RkTDYf0icz+cDIC4ujuHDh3PGGWeQkZHRohnqmlJVW8U/VvyDJxc/yYY9G0hPTOfhaQ+TFJ3E3e/dreF03ZRCX0SkDR1oSNyQIUMYMWIEF110EZGRkQc5U8tsK9/G04ufZs6yOZRWlzIxbSIvX/IyFx97MS5ncJacuMg4DafrpsIa+saY54FzgSJr7YjQtlTgFSAL2AJcbq0tCe37CXAz4AfusNa+F4Zii4g0ae+QuNzcXHbs2FG3PT09nZEjR3LrrbeSkJDQJtfOKcjh8UWP85/V/8FiufjYi/nhCT9kUvqkRsdqOF33Fe6a/j+APwL/rLftfmC+tfYhY8z9oe/vM8YMA64EhgP9gQ+NMUOstf52LrOIdHPV1dV88803rFq1qtGQuJEjR3LxxRfTp0+fw26abyl/wM8ba9/g8UWP88XWL0iMSuTOiXfyg4k/ICs5q02vLZ1TWEPfWvupMSZrv80XAFNCz18AFgD3hba/bK2tATYbYzYCE4CF7VJYEel29g6JW7VqFWvWrGnRkLj2UF5TznPLnuOpnKfYUrqFgckDeeKsJ7jxuBtJjEps17JI5xLumn5T+lhrCwGstYXGmN6h7WnAonrHbQtta8QYMxOYCZCRoR6pInJwe4fErVq1ij179gANh8SdddZZBx0S19a2lG7hqcVP8bdlf6OitoKTMk7i0TMf5YKhF+B0tO8vHtI5dcTQb05T7WSNx50A1to5wByAcePGNXmMiHRP9YfE5efn123fOyTu2muvPaQhcW0hOze7rqNdemI6Nxx3A6uLVvP62tdxGAeXD7+cu0+4m3H9x4W1nNL5dMTQ32mM6Req5fcDikLbtwHp9Y4bAGxv99KJSKfg8/lYv359mw2JayvZudkNVqjLL8/nV5/8ijhXHD8+8cfcNuE2BiQOCHMppbPqiKE/F7geeCj0+Ga97S8aYx4j2JFvMJATlhKKSIdxoCFxgwcPZuTIka06JK6tzZo/q8EY+r1SY1L53bTfhaFE0pWEe8jeSwQ77fU0xmwDfkEw7P9tjLkZyAcuA7DWrjbG/BtYA/iA29RzX6R7KSkpqQv3wsLCulp6eno6I0aM4JZbbiExsfN2ZKusrWx2mtxt5dvauTTSFYW79/5Vzeya2szxswFNGyXSxe0dEpebm8umTZvCNiSuPb3/7fvMnDez2f2aJldaQ0ds3heRbqL+KnFr1qyhpqYG2DckbtKkSUyfPr3dh8S1pz2ePfzo/R/xjxX/YGiPofz81J/zyJePaJpcaRMKfRFpc9ZaioqK6laJqz8kbtCgQS1eJa6reXXNq9z2zm0Ue4qZdfIsfnrKT4mOiGZIjyGaJlfahEJfRNrcl19+yfLlyw9rlbiuqLCikNveuY3X177O8f2O591r3mVM3zF1+zVNrrQVhb6ItLnJkyczefLkcBcj7Ky1/H3F3/nhez+kxl/Dw9Me5u5JdxPh0I9iaR/6lyYi0g42lWxi5ryZzN88n1MzT+XZ855lcI/B4S6WdDMKfRGRVrb/jHonZ57M62tfx2mc/OX//sKMsTNwGEe4iyndkEJfRKQVNTWjXnZuNsf1PY65V83VbHoSVvpVU0SkFTU3o16xu1iBL2Gn0BcRaUXNzai3tXxrO5dEpDGFvohIK6moqSA6IrrJfZpRTzoChb6ISCvYUbmDU/9xKrX+WiKdDRf30Yx60lEo9EVEjtD64vVMem4S64vX89ZVb/H8Bc+TmZSJwZCZlMmc8+Zosh3pENR7X0TkCCzatohzXzwXh3Hw8fUfMz5tPIBCXjok1fRFRA7TvHXzOP2F00mOTmbhzQvrAl+ko1Loi4gchmeXPsuFr1zI8N7D+fLmLzkq9ahwF0nkoBT6IiKHwFrLgwseZOZbMznrqLP4+PqP6R3XO9zFEmkR3dMXEWkhX8DHLW/dwnPLn+OmMTfxl3P/gsvpCnexRFpMoS8i0gJVtVVc8d8reHvD2/zslJ/xyym/xBgT7mKJHBKFvojIQeyq2sW5L53Lku1L+Mv//YXvjfteuIskclgU+iIi+6m/Sl6/hH4EbIDS6lJeu/w1LjjmgnAXT+SwKfRFROrZf5W87RXbAfjFqb9Q4Eunp977IiL1NLdK3j9W/KP9CyPSyhT6IiL15JflH9J2kc5EoS8iUk//hP5NbtcqedIVKPRFRELKa8qbHIanVfKkq1Doi4gQnHjniv9ewY7KHdx/0v1aJU+6JPXeF5Fuz1rLHf+7g3c3vsucc+cwY+wMfjf1d+EulkirU01fRLq9pxY/xZ+X/Jl7T7yXGWNnhLs4Im1GoS8i3dq8dfO4+727ufCYC3lo2kPhLo5Im1Loi0i3tWLHCq569SqO73c8/7roXziMfiRK16Z/4SLSLW2v2M65L55LSkwKc6+aS1xkXLiLJNLm1JFPRLqdqtoqznvpPMpqyvj8xs+bHZsv0tWopi8i3Yo/4Gf6a9NZsWMFL1/yMqP7jg53kaS72pwNb2TBi47g4+bsNr+kavoi0q3c/+H9vLnuTZ46+yn+b8j/hbs40l1tzoacmeAPrfPgzgt+DzCw7eaEUOiLSLcxZ+kcHln4CLePv50fTPxBuIsjHdnmbFg5C9z5EJsBo2cfehhbCzW7wb01+FW1FW/5ZjzF6wgUvE+E8RITGTwswknwF4CVsxT6IiJH6sNNH/L9t7/Pd47+Do+f/Xi4iyMdWUtr4bVl9QI9H3/lFjy71+Mr34KzpoDoQDEuh6/BqY11YugBAS84g9sinPUOcLftwk4KfRHp8tbsWsOl/76UYb2G8fKlLxPh0I8+OYCVs/YFfkjA66b6s5nUrvwjDs82ovy7iHLUNDjGWAMkg6sfNvk46DEEkgZBXDrEpkNsBhHRvUgwoXv47rzG145t24Wd9C9fRLq0oqoizn3xXKIjonnr6rdIjEoMd5GkNRxJ83vAC57tULUV3PnYqnxqSjZSU/ItVOUTWZ1HTGTDlzgcgNcN1bux8YOxqedgU47CxGWGAj0dR0w/4lv6C+Xo2Q1bEwCcscHtbUihLyJdVrWvmgtfvpAdlTv45IZPtDxuV3Gg5vesq6C6qEGze23pt1Tv2UCgMh+Xt5AYynAYW3c6AxCIAWcvbEw61ruTQKAqGPT1xKZkEnvhhtZ5D3t/QTnSfgOHSKEvIl1Kdm42s+bPIq8sj1hXLG6vm/9e9l/Gp40Pd9HkSFkL3lJYfm9d4Pv84PGC3+/G+fF1xERdT4TxN3iZCbjA0ROi+mNTp2F7DoX4rGANPdT0Hu1KIHrvC/b/pQLaphY+cHqbh/z+FPoi0mVk52Yzc95M3N7gD2u3143L4aLaXx3mknVDh9P87nM36OkeqNyCp3g93rJNGE8BMYFdRDpqG7zEmFBNHbA2gB3yw+B99Nh0iMuA2HRckam4jGl8veaEqRbeHhT6ItJlzJo/qy7w9/IGvMyaP4vpIzv/D+xOo6nm98UzoHoX9BgL7q3YqjyqizdQU/otuLcS7dtJtKPh350DIJAIrr6QNAybejQkHw2rfwe1uwFwOiB+bxU9NhPGP9o67yEMtfD2oNAXkS4jv6zp4U7NbZdWFPAFa8WVm2DpHdRUu6n2BlvkXREQG+nBLL+77nADmEAsxtkbGzuAQMqp2NTBmFDtnLgMiEkjzhnZ+FrRfcLSCa4rUOiLSJeRnphOfnnjgFcHvv0cbs/32hKo3ETtnrWUb19BTfE3mKotxPgLSYooadA5LsIJDh/4ASxYQs3wp71X19s92hW/7z76oejCze9tTaEvIl3G6L6jG4V+rCuW2VNVA6xzoJ7vmZfjK99Mcf4SKrYvx7tnHQ7PFhLMblJcJcQ4g30jIoGeQK0jmUBCJq7UaTgSB0P8oODXF9fgrC4gIWa/a8dmQr8zW+d9dNHm97am0BeRLuHTvE95a/1bnJp5KltKt5Bflk9GUgazp87W/fyQQCDArk/uZ/t2N+7QvDJRLkiJc5Py4bUkxV5LhLH0AfoAAePC9szCmThmX6DXfQ0k0pXQ9IWO+72a3zsohb6IdHpl1WVc+/q1DEoZxFtXv0V8ZHy4i9TurLWUlZWRn5/P1q1b2bp1Kzt27MBFDf2iC0mL3k5a9HbSI7ZxXFb91wV7wIOF4T8NBnrCURA/CEdMfzCHsRirmt87LIW+iHR6t71zGwXlBXxx0xddNvDdbnddmG/dupVt27bh9Xox9YaipSbFMWJAgGHJxZw2JJ+YAasxFd+ADQQPiMuE6tgGNfC6l8dmwuhft16B1fzeISn0RaRTeyn3JbJzs/nllF8yccDEcBfnsHi9XgoKCvbV0lfOxZP3HtZbBhFJ0Pc0YvtPJD09nfT0dKZMmUJa/75Eer6F4q9gz1fBx9KvwVcLu4Ho3pA6HjIvCz72GBfc1l4Tz0iHpNAXkU4rvyyfW9++lUkDJvHAyQ+EuzhNCgQC7Nixo0EtvaSkpEEN3eVy0b9/f9LT0xnXdzsX+uYRN8YT2lsGjvdg9MkQ44Q9b8LGr+CrZeCrCp0gEVLHwtC7oMd46DEh2EO+qQlp1PTerSn0RaRT8gf8XPf6dfitn39d/K+wrJxnrWXPnj11YZ6fn09RUVGDYxwOB3369CE9PZ0hQ4YwdepUkpOTG4R+A2/MAJen4baAB5b/KHTCKEg5DgbdFAz41PGQOOTQ7r2r6b3bUuiLSKf06MJH+STvE/5+wd8ZlDKoTa5RWVnZINC3b9+O399wXvfU1FQyMjJIT09nwoQJ9OrVC8f+K7W0hA3Ari+aXm51r7OXQfIIcLgO/fwiKPRFpBNaXricn370Uy459hKuH339YZ2jpqaGbdu2NegYV1NTg7X7JpiJj4+vu49+5pln0r9/fyIiWvHHprVQvBjyXoH8/4CngOAUNrbxsbGZkHpc611buiWFvoh0Km6vm+mvTadXXC/+eu5fm2wm9/v9FBYWNriPXl5e3uCYqKgo0tLSSE9P54QTTmDAgAHExOw/m0wbsBZKloWC/t9QlQeOSOh3NmT+IbjozNI71NFO2oRCX0Q6lXs/uJdvdn/DHyf+kc/e/4ytW7eye/fuBsc4nU769u1LRkYGw4cP5+yzzyYxMbH9Crn/NLejfgMpoyD/Fcj7N1RuBBMBfc+Akb+EARdAZHK9NxCtjnbSJhT6ItKhlJWVNbiPXlhYSCAQwBjDBjbwkuMlpkRNYUDNANIz0pk8eTI9e/ZsvmNce2tqmttF1wE22Nmuz+kw7D5IvwiiejR9DnW0kzbSYUPfGLMFqCC4XoPPWjvOGJMKvAJkAVuAy621JeEqo4gcGo/HU3cfPT8/n4KCAmpraxvcR09KSqq7jz5y5Ej69euH0+mkqKqIUX8exci4kfxvxv+IjjispVra3ooHGjbNA2DBlQrnfRMcKy8SJh029ENOs9bWb7e7H5hvrX3IGHN/6Pv7wlM0ETkYay2//OUv60I9OjqaAQMGkJGRwSmnnEJaWhpRUVEtOs+MeTMorS7lg2s/6JiBby1sfwc8zSzj6y1R4EvYdfTQ398FwJTQ8xeABSj0RTosYwwPPvjgEZ/n2WXPMnfdXB4/63FG9hl55AVrTdZC4buQ+yAU5wTv1Vtf4+NitbyvhN9hDCZtNxZ43xiz1BgTWveRPtbaQoDQo35tFuni1hev5+737mbaoGncMfGOcBdnH2uh8H14/0RYcA5U74QJz8KEvwV729en3vfSQXTkmv5ka+12Y0xv4ANjzNqWvjD0S8JMgIwM/XYt0ll5/V6uee0aoiOieeHCF3Aczopvrc1a2Dkfvv4F7P4yON3thL/CwBvAGRk8xhGh3vfSIXXY0LfWbg89FhljXgcmADuNMf2stYXGmH5AUTOvnQPMARg3blwTs1yISEeWnZvNrPmzyCsLzk5358Q76Z/QP8ylAnZ+HAz7XZ9BTBqM/1NwOlznfv0S1PteOqgOGfrGmDjAYa2tCD0/E/gVMBe4Hngo9Phm+EopIm0hOzebmfNm4vbu6wH/7LJnGZ82nukj2ylI9x9nn3UN7P4cij6BmP4w9mk4+rvB8fQinUiHDH2gD/B6aNxtBPCitfZdY8xXwL+NMTcD+cBlYSyjiLSBWfNnNQh8CM7CN2v+rPYJ/abG2a+ZDa4kGPskHD1TYS+dVocMfWvtJmB0E9uLgantXyIRaS/5ZU0PeWtue6tbOauJcfZARCIM7UAdCUUOQwfoFSMisk96YnqT2zOS2qlTbnOr3Hm2tc/1RdqQQl9EOpTLR1zeaFusK5bZU9t4yJu3Ahbd1Px+jbOXLkChLyIdytLtS0mJTiEjMQODITMpkznnzWnb+/m7FsL/xsDmFyDtAnDst9qextlLF9Eh7+mLSPe0vHA5H2/5mIenPcy9k+9t+wsGfLDqN7D6NxA7AKZ+Ar1Patx7X+PspYtQ6ItIh/HowkeJj4xnxtgZbX+xio3w5TVQvBiyroVxT0NkUnCfxtlLF6XQF5EOYWvZVl5e9TJ3TLyD5OjktruQtbDp77D0DjAumPwyZF7RdtcT6UAU+iLSITy1+CkgOPtem6kpDo7B3/oa9DkNTngB4poeLSDSFSn0RSTsymvKmbNsDpcNv4zM5My2uUjhB7DoeqjZDWMehmN/BB1hLn+RdqTQF5Gw+9uyv1FeU86PJv2odU5YvyNeTDokDwsuf5t4LEx5B1LGtM51RDoZhb6IhJXX7+XJxU9ySuYpjOs/7shPuP80up784FffM+GUNyAi5oAvF+nKFPoiElb/XfNf8svy+eN3/tg6J2xuGt3ydQp86VBys3OZP2s+ZfllJGUkMXX2VEZOH9mm11Toi0jYWGt5dOGjDO0xlP8b8n+tc9LmptF1t9Pc/SItkJudy7yZ8/C6vQCU5ZUxb+Y8gDYNfoW+iITNp3mfsrRwKX899684WqNTXcAHzjjwVzXep2l05RAcaS3cV+3Ds8eDp8QTfNzjobqkGs8eD1XFVXz55Jd43B4AIonEgQOv28v8WfMV+iLSNT2y8BF6xfbi2lHXHvnJ/NXwxZXBwDcusN59+zSNbpfTlk3jTdXC37z5TXZ+vZM+Y/pQXVKNu9hNRVEFZTvLqNxVibvYTU1pDd4KL4HKAPgOcAEDXuvFYLDYBrvK8sta5T00R6EvImGxdvda3lr/Fr849RfEuI7wXru3Aj69AHZ+DGOfhsgUTaPbhR1q03hdrXu/mrdnj4eq3VWU7SijoqiCql1VeEo9lH5biiPgwIkTgwHAX+Pni4e/aHjiSHDGO3EluIhOiSZ1cCoJvRNI6ptEQq8EYnrEEJMSQ0xqDNEp0cSkBp9HJUTx5KAnKctrHPBJGUmt/Gk1pNAXkbB4bOFjREdE8/3x3z+yE1XvhgXfgZLlMOlf+8JdIR9WrV0TtwFLdVmwefy9e96rC3wAi8Xj9vDfmf/ly+e+DNa6S0K17qoAeA9wYgMmzhCREEFUUhRxPeMo2VDSqAa+1/fXfD8Y3ikxOCOdh/1+ps6e2uAXFwBXrIups6ce9jlbQqEvIu2uqKqIf678J9ePvp7ecb0P/0RVW+HjM6FqC5z8Ogw4r9XKKIevuZq4tZZjLjqmwf3tvTXv6pJqKndVBmvduyqo2l1FdUk1tWW1+Kv8WHfTIQxgQn98bh+lm0qJTommx5AexPeKD9a6eyfU1bL31rz3fkUmRGKMaXC+J7KeaLoWnplEr2N7tcpntPcXIPXeF5Eu709f/Ykafw0/nPTDwz9J+Xr46AzwlsJp70HvU1qtfF3dkdbCrbXUVtYGQ7vYg7vY3eD5pw9/isftaVBb9rq9vHrtqzgOtKK7A0xsqNadHEVc7zh6D+tNYp9EEnsnEtczjuiUaD649wPcuxoOy3TiJDUzlbu23HWoH0cj7VULHzl9ZJuH/P4U+iLSrjxeD8989QznDTmPoT2HHt5J9iyHj88KPp/6MaQe33oF7OKaqoXP/e5cyvLLSDshrS68K3ZWULK9hPKd5VTuqsSzx0NtaS2+Ch+2ykKg+Wv48dc9t9i6DmsWy7TfT9t3j3u/+92R8Y1r3U1xRjjbNJTDVQtvDwp9EWlX/1z5T3a7dx/+lLtFn8En54IrCU7/ABIP8xeHDqw17od7PV7KtpWxY/0Odm7cSfHmYkq3lbJh3gZsrSWCCCKJxIkTX7WP+Q/Mb3wSF0QkRBCZFElMagw9MnoEa919E4ntEUtsj9hgM3mPmH3PU2N4esjTzTaPT/7x5MP9WOq0RyiHoxbeHhT6ItJuAjbAY4seY1z/cZySeRjN8QVvw+eXQlwmnPZBl1whr6ma+GszXqOopIg+U/pQvKuYXVt2sXvzbsq3lePe6YYKGn6Vg6luXGN2RDmIro0OPg/9qe+6j64jJjUU4D1icMW4Dus9tEfzeFcN5bam0BeRdvPW+rdYX7yely55qUXNuA1seREWXg/Jo+C0dyG6dTpUtRe/309ZWRl79uxhz549lJSUNHj0eoMBufiJxTjcDuKII4YYoogi2hPN+3e9T2JSIr5SX4OmdYPBOAzxfeNJSEsgcWQi8f3jSUxLJKF/AglpCST0TyAxLZGopCieHNjMULHMJAaeNrBV3mtXbh7v7BT6ItJuHl34KBlJGVw67NJDe+H6P8GS24Od9U6dC67EtingQVhr8Xg8TYb2nj17qKioAGBn7k42f7SZmrIaopKiGHj6QPqP6U9SUhKpqamkpqaSnJhMKqn4a/14dnvYs34PRauKmFQ+qemL+2HE+SOCIV4vyBP6JxDXJw6Hs2UzGnblTmpycAp9EWkXXxV8xad5n/LYmY8R4Wjhjx5rYfVs+PpnkHYeTH6lVRbN8fv9lJaWNlvr9vman04tNjaW1NRUUlJSSE1NJTMzk+OOO47U1FTi4+NZ9eIq5j0yj37ufsEXlEHE/yKYcPQE4mLjKJpfRMGqAlasWYHPs+86yQOT6T2iN2X5ZdSU1zS6blJmEhf8/YIjfu+qhXdvCn0RaRePLnyUxKhEbj7+5gMfuDk7NJteHkQkgK8Csq6FE54Dx757zNZa3G53o9r23ueVlZXNXsLpdJKcnFxX605JSWHgwIGkpKSQkpKCy3V497IBPrj/gwa1aAjOCPflH74EIKF/Ar2G92LcLePoPaI3vUf0ptewXkTGRwKN7+mD7odL61Hoi0ibys7N5r4P7qOgooDEyETmrZ/H9JH7Zsvz+XyUlJQEw3rVS5Tk/I6SihqOy4Jj0ypYtNHw3uIK7PuN587fW+veG9xHHXVUXQ08Li7u0PsNHAb3bjd5n+ax+ePNbPl4CxXbKpo99sfFPyYm9cAtFaqJS1tS6ItIm8nOzebmP99MTW0NpEH5gnJu+OAGXh3yKqP6jgIgIiJiX6177Z9Ija3hpMGQEAP+AJxwtOWEUcvhwtfbrdwHGjLnKfGQ90keWxZsYcvHW9j59U4gWBvPOCmD8oJyakqbbp4/WODvpZq4tBVjbfNTG3YF48aNs0uWLAl3MUS6lSlTpgCw5cIt5JXkgQHqVbozkzLZcteWxi98sbmauYGrDzAbTCtqqnndGelk4OkDqdxRyY6VO8BCREwEGZMzyJySycDTBtJ/XH+ckc5mm+fPm3OeglzahTFmqbV2XFP7VNMXkTaTX5ZPU7Ou5pflN/2CiETwlTfeHpvRugU7gA8f+LDRPXl/rZ+N724k67QspvxyCllTskibkEZEVOMfoWqel45MoS8ibSYjKYO8srwmtzfyzSPBwDcRYOv1nnfGBpfGbWN7Nu4h98VcyvOb+KUj5PqPrm/RudQ8Lx3VIYe+MSYaOM9a+582KI+IdCGzp87mutevI2D3Nc3HumKZPXW/EP/2eVh+L2RcDv3PDQ7Rc+cHa/ijZ7fZMrmVOytZ/cpqcrNzKcgpAAPOKCf+Gn+jY5My23adc5H20KLQN8Y4gTOBq4CzgM8Ahb6IHNCVw69kxtwZOB1OqmqryEjKYPbU2Q1677P1NciZAX3PhEn/D5yRMOjaNitTTUUNa19fS252Lps+3IQNWPqO6cu0h6cx8qqR5H2SF5Z1zkXawwFD3xhzCnA18H9ADjAZGGitdR/odSIiACt3rsTj85B9cTZXj7y68QE7PoQvroIeE+GU14KBf4Sa6nk/7LJhbHxvI7nZuaybuw6fx0dyVjKT75/MqOmj6DVs35S+uicvXVmzoW+M2QbkA38G7rXWVhhjNivwRaSlPtnyCQCnZp7aeOfuxfDphZAwBE59CyLijvh6TS1W88YNb/DWLW9RW1lLbM9Yxtw4hlHTRzFg0oBmx/Hrnrx0VQeq6b8KXAhcAfiNMW8CXXt8n4i0qk/yPuGolKNIS0xruKN0NSw4B6L7wOnvQ1Rqq1xv/qz5jXreB3wBbMBy9dtXM+iMQThdzla5lkhn1OwKDdbaO4Es4DHgNGA90MsYc7kxJr59iicindln+Z81ruVXboGPzwRHJJz+AcT0a5VruXe7m1w9DoJryw8+Z7ACX7q9A97Tt8GZez4CPjLGuICzCXbm+xPQs+2LJyKdVVVtFXs8ezgl85R9Gz074aMzwOeGMz6F+EFHfJ3ygnIWPrqQpX9d2uwxSRnqeS8CLei9b4xxWmv91lovMA+YZ4zRzS4ROaDS6lIATs0K1fRrS+Hjs8CzHU7/EJKP7MdI8YZivnj4C1a+sBIbsMEOecN78ckvP1HPe5FmtGTI3lvGmEustW5jTCTwC4I1/rFtWzQR6cxKa0rJSMogKzkrWLP/5DwoXwOnzINezawZ3wI7Vu7g8999zpr/rMHhcnD8jOOZfO9kkrOSAUhMS1TPe5FmtCT0XwA+MMY8DPyK4Pj8iW1aKhHp9Mo8JXzH4QnOp++IgUA1TH4Z+p91WOfb+uVWPvvtZ2x4ewORCZFMumcSk+6eRHzfhl2M1PNepHkHDX1r7cvGmCLgNeBqa+07bV8sEenM3JVb8Qb8nOqqCm4IeMBEQsB74BfScJx9YnoiI68cybbF28j7JI+YHjGc9uvTGH/beGJSWrZinYjs05J7+k+Fnq4AnjfG/BvAWntHG5ZLRDqx0orggjoX1h96b2th5awDTqm7/zj78vxyvnj4C6JTojnr8bM4fsbxRMYd+QQ+It1VS5r3l+73KCJyQKU+H0kOSN1/hJy7mdX1QpoaZw8QmRDJCXed0IolFOmeWtK8/0J7FEREugZrLWUBGB0FjSa8O8ASucXri5sdZ1++tfmV70Sk5VrSvL9p/00Eh/Af+QBbEelyvi35lh5OiN1/6q9mlsj17PHwya8+4atnvgr9dGl8To2zF2kdLWne/wroA7xIcJx+bZuWSEQ6tZyNbzLIBX5HNMT0Bs/WJpfI9Xv9LPnzEhY8uICashqO++5x9B3dlw/u/UDj7EXaSEua968wxqQQXG3v78BCa+0v27xkItIpHbX5aRyAI2UUXLS40X5rLevfWs8H93xA8fpiBk0bxJmPnUmfkX0AiE6K1jh7kTbSkpo+QAAttiMiB7NtHhP9eew2sfR2xjRa5nbszLFs/mgzm+dvpsfQHlz11lUMPmdwg9XuNM5epO205J5+NtAfeAm4Aag1xqRaa/e0cdlEpDPxVuLLuYW1NeCL7kfVzqpGy9x+NOsjIuMj+c7T32Hs98ZqARyRdtaSmv5kgrX8nwD3s6+rjTryicg+X/8MZ3UhM4ogKTqFkjUleGsaD7+LTo5mwu0TwlBAEWnJPf2sdiiHiHRmxV/B+qf4OGIo600RIyPj2FWzq8lDyws0/E4kXFp6T19EpGkBL553r8VfnsCNG3bQv3QAezY0f/dPw+9Ewmf/kbQiIi2zORveyKLwoQxizDr++9Ip7Ih2c8q8UyjfXo7T5SQiumG9QsPvRMJLoS8ih25zNuTMBHcePfvvZu2SoXy8NZlb/nILvXb3AsA4DOf/7XySMpPAQFJmEufNOU8980XC6JCb940x/YA91tqaNiiPiHQGK2eB3w1A8fYevPrHSxnkdVGSXIINje711fg0/E6kgzmcmv7/A9YaYx5p7cKISCfhzqPGE8mHL09jzk+/h9cbwReTvuDds97FEfqx4o/wh7mQIrK/Q67pW2unmeBMGsPaoDwi0sHZgGXlFycy/8VJVJYmkDF0C4tqAnxw1gecO+/c4DHGUpJcEuaSisj+WjI5T5q1tmC/zd+z1v6ljcp0sPKcDTwJOIG/WWsfCkc5RLqjbYu28e6tL1Gw4kzSjtrGlT98mbSjCvjHqgTiPDGMWjmK0qRSil3F+JJ84S6uiOynJTX9t40xV1pr1xpjhgJzgG/auFxNMsY4gWeAM4BtwFfGmLnW2jXhKI9Id1GxvYIP7/+Qr//f18SnVHLhj9YzauZ5mDUfsKYUPo2qYEQk/PanvwXA8YKDISlDwltoEWmkJaF/FfCyMeZj4DTgDmvtp21brGZNADZaazcBGGNeBi4AFPoiraT+fPmJ6YkMmDiADe9sIOD1c9IlKznp4mVEXfgFxKbBkO/y+NwZuLa+QGlUH0xtARlJGcT0iKFPXJ9wvxUR2Y+x9uDr6Bhj+gP/A35nrX25zUvVfDkuBc621n439P21wERr7e3NvSYhIcGOHTu2vYoo0qlV7ayieH0xgUCgwfao+Eh6pe0iwlEOqcdBRDwAtf5aFm1bRN/4vgzpsa9mv2LFCgDGjBnTXkUXkZBPPvlkqbV2XFP7Dtp73xiTC7wLJAL/zxjztTHm61YuY0uZJrY1+q3FGDPTGLPEGLPE620897eINK1kc0mjwAdISikkwlEKScfWBT5AQUUB1lrSE9PbsZQicrha0rx/bpuXouW2AfV/ugwAtu9/kLV2DsG+B4wbN84uWLCgXQon0pnt+XYPTx/9dKPtp168gCmX5MGYh2DYfXXbq2qrSH88nYuyLuK1K15r8JopU6YAoP97Iu2v/lLV+ztoTd9amwckA+eFvpJD28LhK2CwMWagMSYSuBKYG6ayiHQJlTsreef2d3jmmGcataWNmJTLlEsWsHrpBDj2xw32/X3F3ympLuGeE+9pv8KKyBFpyZC9O4EZwN5f5f9ljJljrW1cJWhj1lqfMeZ24D2CQ/aet9aubu9yiHQFNRU1fPnIlyx8dCG+ah/Hzzie3sN78+GP38XrsaQdvZULZr5B/vpM7PF/hnq1B3/Az2MLH2PSgEmcmH5iGN+FiByKljTv30yws1wVgDHm98BCoN1DH8Ba+w7wTjiuLdIV+Gp8LP3rUj79zae4d7kZdtkwTv/N6fQY0gM2ZxOzbi45743mirtfobI8nsqyZEac+A1wfN05Xl/7OptLN/PImZqYU6QzaUnoG6D+fJp+mu5QJyIdmA1YVr28io9++hGlm0vJOi2LaQ9NI21C2r6DVs5i5KQ8Rk5aWrcpuefK4Fz7A6cHz2Mtf/jyDxyVchQXDL2gvd+GiByBZkPfGBNhrfUBfwcWG2NeD+26EHiuHcomIq3AWsu373/L/Pvns2PFDvqM7sP0d6dz1JlHNe7w426mu447v+7p5/mfk1OQwzPnPIPT4WzDkotIaztQTT8HON5a+5gxZgFwEsEa/o3W2uXtUTgROTT1J9ZJykjiuJuOI++TPDZ/tJnkgclcnH0xI64cgXE00VjnrQBHFASaWEAzNqPu6SMLH6FHTA9uGHND270REWkTBwr9up8K1tplwLK2L46IHK7c7FzmzZyH1x2cm6Isr4wFv1hAVGIUZz91NuO+Nw5nZDM1c58bPjkPAl5wREKgdt8+ZyyMng3Aut3rmLtuLj8/5efEumLb+i2JSCs7UOj3Msb8sLmd1trH2qA8InKYPnzgw7rAry8qKYqJP5jY/At9Hvj0Atj1GZz4L7CB4D18d36whj96dt39/McWPkaUM4rbJtzWVm9DRNrQgULfCcSjTnsiHZqv2seqV1ZRnl/e5P7ybU1vB8BfA59dDDvmwwl/h6yrgttDIV9fUVURL6x8getHX0/vuN6tUXQRaWcHCv1Ca+2v2q0kInJIygvKWfLnJSydsxT3LjcOl4OAt4kpdDOSmj6BvxY+vxwK34UJz8Kg6w94vWdynqHGX8MPJzXbACgiHVyL7umLSMdgrWXrl1vJeSqHb177hoA/wNDzhjLhjglUFlby1vfeatDE74p1MXX21MYnCvjgy6uhYC6M+yMc/d0DXtftdfPMV89w/tDzGdpzaGu/LRFpJwcK/SZ+UohIOOxtws95KofCZYXB+/R3TmT898eTMiil7jhjTIPe+1NnT2Xk9JENTxbww8JrYeurcPxjMOTg9+f/seIfFHuKuWeSptwV6cyaDX1r7Z72LIiINB5yd+KPTqRyZyVL/7oU9243vYb14v/+8n+MumYUkXGRjV4/cvrIxiFfnw3A4psh7+XgAjrH3H3QMu2dcndC2gROyjjpSN6eiIRZS2bkE5F20NSQu//d8T8wMPT8oUz4wQQGnj7wgCtoHZANQM73YPMLMPKXDVbMO5A3173JtyXf8tC0hw7/2iLSISj0RTqAgD/Aez96r8khdwn9ErjyjSuP7ALWwpIfwLd/g+EPwIiftfilj3z5CAOTB3LRMRcdWRlEJOwU+iJh4qv2sfmjzax9Yy3r5q6jamdVk8dVFFYc2YWshWU/gg1/gmPvgVG/abBi3oF8ufVLFm5byNPfeVpT7op0AQp9kXZUXVrNhnc2sPaNtWz830ZqK2uJjI9k8DmD2fTRJjy7PY1e0+yQu5awFlb+BNY9DkPugDEPtzjwAf7w5R9IiU7hxjE3Hn4ZRKTDUOiLtKL9O+JNnT2VzFMzWTd3HWvfWMuWj7cQ8AWI6xPHiKtHcMyFxzDw9IFEREU0uqcPBxhy1+ICPQhrfg9Hfw/GPnFIgb++eD1vrn2TB05+gLjIuMMvg4h0GAp9kVbSVEe81697HRuwAPQY0oMTfngCx1x4DAMmDmi06M3eXvcHHXJ3IJuz902h60oCbykMuhHG/+mQAh/g8YWP43K6uH3C7Yf0OhHpuBT6IkeotqqWwqWFvHP7O4064tmAJSo5iu8u/C49j+l50HMddMjdgWzOhpyZ4HcHv/eWgnFC79PAOFp8muzcbO7/4H62VWwj3hXP/M3zmT6y8bS8ItL5KPRF6mmqeb5+CAf8AXZ/s5tti7dRsLiAgpwCilYVYf222XPWlNW0KPCP2MpZ+wJ/L+uHr38Gg65t0Smyc7OZOW8mbm/wPJXeSmbOmwmg4BfpAhT6IiFNNc/P/e5cCr4qICI6goLFBWxfsp3ayuCys9HJ0aRNSGPo+UNJm5DGW7e+RcW2xj3tj6gj3qFw5zWzPf+gL/V6vVRVVXHfq/fh3uOGZKAGKAB3rZs7V91J6eRSqqqqcLvdWBv8JccY0+D50UcfzTXXXNM670dEWp1CX7o9r8dLRUFFk+PkfdU+Fj+5GIfLQd8xfRl9/WjSJqYxYOIAUo9ObXBf/oyyM1q/I1491lo8Hg9VVVVUVlZSVVVV99xXsZVxVRGkxvr4fB188g349q6940rErvsFQIPJdeqHdUREBHFxcRQsLYBMwAGUANFAIhRHFnPZZZcRHx9PTEyMJukR6aQU+tJpHKzpvSlej5fybeXBr63Bx7KtZVRsq6h77iluPExufz8p/wkR0Qf+77K3LO8/8D6783cTnRbNyLtG4hvm47PPPmsQ0nsf69ea96pfe66/DSAmJoa4uDji4+OJi4sjLi6OPpHbGVk7i4h4B8UVLkake5l4NERGAM5YmPBMk0vlNuVh/8PsqNoR/CZt3/bMpEx699ZyuiKdnUJfWs3hhPKhnLtR0/uMuVQUVtDv+H51AV4/3Mu3lTcZ6DE9YkhKTyKuXxw9xvQgsmckESkRfPzIx+wp2UM11dRQQy21ePHiSHbw69/9Gmi+plx/m/NmJ73iexEfH095XDl5eXnEx8eTlJREWlpaXVjvrTU7HC3vZNfI1tfhyzshqgecupiepav39d6PzYDRs1sc+L6Aj6iIKAwGy75fOmJdscyeOvvwyygiHYZCvxtp71CeN3MeQKNrWGvxur3UVtRSU1FzwMfqsmqqSqpY9e9V1HpqsVgMBidO8MAH937Q8NyxFhIJfg0ERoNNtJAAJslAIrhdbty4iY6OblBrHnjzQEqfLiWuJo4UUogkkriYOC564iKOv/74VvmcWo21sOYhWPkA9JgIp7wBMX0hZUyLQ35/Tyx6gryyPO6ceCdvrH2D/LJ8MpIymD11tjrxiXQRZv9mxK5m3LhxdsmSJeEuxkG1ZSDvPX9T95vPm3Nes9ex1uKv8VNbVYvX7cVb5Q0+3/u4d1tlLR/O+pDK0koCBLD1/jhdTpIzkvG5ffjdfmy1xdZajG3hPWEXmGiDp8LT6NwABsOV/72SXkf1ovfRvYmOjz7iz6kt/x5ahb8mODRv8z8h8yqY+BxExBzRKTeXbGb4n4ZzxlFn8MYVbxzxPfspU6YAsGDBgiM6j4gcOmPMUmvtuKb2qabfQh2llry/gD/QOIyrvMGadL1tH973ITXuGgL1/njcHl6++WXmPzIfr9uLz+PD7/ETqAlgaw8xnIFIIuvOXRfNXktyv2SiE6OJSY4hNiWWmKQYohKjiEqIIjIhstnHyPhIHM5g0/cTWU9QllfW6JpJmUmMvWRsi8t4MEc0Tr49VO+Czy6CXV/AyF/BiJ8e8qQ7+7PW8v13vo/T4eSP3/mjOumJdGEK/RY4klCGUDC7gzXivTXjukCurOV/d/6vQQ08QAC3282rM18l5//lBJu6Q8f73D781X4C1cFgxtfy9+HAgSH4A91gguFcY8EPcSlxRGVEEZ0QTUxSDDGJMUTGRxIZF4krztXwMdbVaNuzE56lfGs5jtCfvZIyk/jeZ99reSGbMXX21DbtGd8plK6CT86D6h0w+RXIvLxVTvvK6ld4d+O7PHn2k6QnpbfKOUWkY1Lot8D8WfMbDeWqddfy5i1vsuadNXjKPdSU11BTWVNXy/Z7/Pirg83ZhxLMQF3Tda27lqJVRUTEROCKcxHXO46ohFAwJ8YQmxxLdGJ008G83+Nzk56jfFs5BkNEvb/2pMwk7vr6riP9iJj2u2ltGsqtMkVtZ1bwDnxxJbjiYdqn0GN8q5x2j2cPd757J+P7j+e28be1yjlFpONS6LdAWX6wWdmLFxcuIBjMNZU1rH9nPc5oJ84YJ644F1HJUSQMSAg2ZyfGEJMUE6xFJ0QRGRdssnbFuRrUorPPyaZye2Xd9ZyhP0mZSdy15a5WeQ/THur8odzhm97bgrWw7glYfg8kj4ZT50LsgFY7/X0f3Eexu5j3rnlPS+eKdAMK/RZIykiiLK+sLvAh2FSekpnSKqF85sNntnnTtUK5k6i/YE5MOiQcBUUfQ/rFMOmfENF6q919mvcpf1v+N+498V7G9B3TaucVkY5Lod8CbX0/ub2arhXKHdz+C+Z48oNfaefDSf85pEVzDqbGV8PMeTPJSs7iF6f+otXOKyIdm0K/BVRLlnbR1II5ACUrWzXwAR76/CHWFa/jf9P/R1xk67UeiEjHptBvIYWytLkjWDDnUKzdvZbffv5brhpxFWcffXarnltEOrbWrT6IyOEp+gxMM7+Dx2a02mUCNsDMeTOJc8Xx+FmPt9p5RaRzUOiLhFNtKeR8Dz48BSJTwBHVcL8zNjh/fit5fvnzfJb/GX844w/0ie/TaucVkc5BoS8SDtZC/n/hrWPh27/BMT+CCzYHp9SNzQRM8HHCnMOeS39/Oyt3cu8H93Jq5qncdNxNrXJOEelcdE9fpL1VbYUlt0HBPEg5Dqa8DamhBX0GTm+1kN/f3e/djdvr5q/n/lVT7Yp0Uwp9kfYS8MOGPwVXxrMBOO4RGHonONr+v+G7G9/lpVUv8eCpDzK059A2v56IdEwKfZH2UJoLi2dA8WLodxaM/zPED2yXS1fVVnHr27dyTM9juP+k+9vlmiLSMSn0RdqSzwOrfg3f/AEik2HSvyDr6iNeGe9QPLjgQbaUbuHTGz4lKiLq4C8QkS5LoS/SmupPoxvVGwxQvRMG3RBszo/q0a7FWbFjBY8vepwZx8/g5MyT2/XaItLxKPRFWsv+0+jW7AQMDLsfxvyu3YvjD/iZMW8GPWN78vtpv2/364tIx6MheyKtZfm9TUyja2HLS2Epzh9z/siS7Ut44uwnSIlJCUsZRKRjUU1f5EgEfFAwF9b/EaoLmz6mlafRPZDs3GxmzZ9FXlkeBsPoPqO5YvgV7XZ9EenYVNMXORzVRbBqNswdCJ9dAhXfgiu56WNbcRrdA8nOzWbmvJnklQXn8LdY1hWv48VVL7bL9UWk41NNX6SlrIXinGCtPv/fEKiFvtNg3B+h/7mQ93LDe/rQ6tPoHsis+bNwexveXqj2VTNr/iymjzz0CX+stZSWlrJ9+3a2b99OQUEB27dvp7a2FmstAMYYbr75ZtLT01vlPYhI21Loi9RXv/d9bEYwsDMugbxXgmG/ZwlEJMDR34PB34ekY/a9du9Mevu/vo1m2NtfflnTtxGa2u52u+vCfG+gl5eXY4xpEOjJycn079+ftLQ0Tj31VPr160d0dHSbvg8RaTsKfZG99u99786DRTfAV7eArxISj4Vxz8DAa8GV0PQ52nAa3YNJT0zfF/A7gFKgHBIDiTz44IMNwjwmJqYuzMeMGcM555xDQkKCpucV6eIU+iJ7rZzVuPe99QWnzJ36EfSe0q6T6tQXCATYvXt3g5r5jh078Pv9dWEeUxUDCcAmwAekQMxRMTx9+dNcO/rasJRbRDoWhb50b9W7YdfnsOuzYM2+KX4P9DmtTS5vraW8vLzRfXOPx9Ogqd3hcNCzZ8+62vmoUaPo06cPLpcLgE/zPuU3L/yGk9NPJj8pn/yyfDKSMpg9dfZh3c8Xka5JoS/dS1U+FH0WDPldn0HZmuB2R1TwK1DT+DWH2fve4/FQWFjYIMxLSkoaNaEnJSXVhfnkyZPp168fsbGxLb5OiaeEa167hkEpg3j76rdJiGrm1oOIdHsKfem6rIXytcFw3xv0VaHavCsRek6GrGug18nQY3xwffsW9L73+XwUFRXVBXlBQQFFRUVYaxvcN4+KiqoL8xEjRnDWWWeRlJTUqvfNrbXMfGsmhZWFLLx5oQJfRA5IoS+dR1M96+t3mgv4oGRFvZD/HGp2BfdF9wmG+zE/DD4mjwKHs8HpbdbV7CmpYPtnv2J7YSEFlT0ojD6d2jXrgV/UHed0OunTpw9paWmkp6czceJEevXqhdPZ8Hzt4bnlz/HfNf/l99N+z7j+49r9+iLSuSj0pXNoqmf94hlQ9g1ERAdDfveXwV72APGDoP850Ptk6HUylaZfMMgLCti+ZA3bt39IVVVVXc0cgrXz1NRU+qc9RdqENI7p35++ffsSFdUxV6Zbu3std757J1MHTuWeE+8Jd3FEpBNQ6EvHF/DB8vsa96wPeGBNsOm9ynUUhYHJbCjvy+pdKVT4EkMH5QPZxMfH1zW1jx8/nn79+pGQ0Hmbwmt8NVz16lXERMTwz4v+icNock0ROTiFvrSOgzW9H0jAB54CqNwCVcGvQMVmvKUboGoLLu8OHATqDvcHYHsJ5O0OPtaOeZKeqUNIS0tjQv/+nJ2a2uXHm/9k/k9YsWMFc6+cS/+E/uEujoh0Egp9OXJNNb3nzAw+HzgdAv5gqFdtwVZupnr3WqqL12IrN+Oq3U4sxTjNvlAPWKj0JVBleuOLGowz5Sxii14nwVGK0wlOB6T3CH4RmwkX3tH+7zmM3t34Lo8vepzbxt/GeUPPC3dxRKQTUejLkbEWVtwPfjfVtVDmBo8X/AE3EW/dQELC7SS5yutC3QAxgCUFX9QATK9J2B7HQtLREJ8FcVk4YtNJdEaRWP86m6eGdV77jmJn5U6uf+N6RvQewR/O+EO4iyMinUyHC31jzIPADCDU7ZoHrLXvhPb9BLgZ8AN3WGvfC0shO6tDaYK3AajeBdWF+Cq2Ur5zHVXFG6ktywNPIZH+YuKdlcRHVOBy+AGIjgx+AXj9YPARMeicujDf95VBrPMQO8eFeV77jiBgA9zw5g2U15Qz/7r5xLhiwl0kEelkOlzohzxurX2k/gZjzDDgSmA40B/40BgzxFrrD0cBO50mmuADC2+i8tu3cAeSqC3fgnUX4vLtJj6igjhnJU4T7NkeAaSGvnzOREjphzN+KCamH8T0g41/A29Jg8u5nASb3idnt957COO89h3BU4uf4t2N7/LMOc8woveIcBdHRDqhjhr6TbkAeNlaWwNsNsZsBCYAC8NbrFZyJB3hrIXaEqjeAZ4dWE8hNWV5VO3+ltryfKxnBxGVa4iLDBBXr4LtoJbEopdJBLzOFGxSXyISRuKI7R8M8+hQqO/9iu5LREQTtcvk0Wp6b2PLC5dz34f3cf7Q87l13K3hLo6IdFIdNfRvN8ZcBywBfmStLQHSgEX1jtkW2tb5NdcRzl8LfaeAZ0cw0Kt34KvYhnvPJrwV28CzA5e/mFhTToRjX4OHAaIBFy68zh4Q34eIQICI5uaOubIWl8N1+OVX03ubqqqt4urXrqZnbE+eO/+5Lj8yQUTaTlhC3xjzIdC3iV2zgD8DvwZs6PFR4CaCWbY/28Q2jDEzgZkAGRmHN296I0dSE6/PWvBVgHs7VBcGH5fcTsDrptoHPn/wjUZGuInKuanRy53W4CIBE9ETE9+XyKRxOBPTIaY/RPeFmL51j05XEs69AfFGVtMLysRmwpEE/l7dvOm9Ld393t2s272OD679gJ6xPcNdHBHpxMIS+tbaaS05zhjzLPBW6NttQHq93QOA7c2cfw4wB2DcuHFN/mJwSA42JC14UfCWg2c7eArrHq2ngNqyPLzl+ZjqHUT6d+MytY0u4QuANxT6DhMclhYIgGPS8w2C3ET1IsZxGH9to2erCb4TenXNqzy77Fnum3wfUwdNDXdxRKST63DN+8aYftbawtC3FwGrQs/nAi8aYx4j2JFvMJDTLoVqap11vxvvwhl4cx8P9mb37SaiiTD3Bly4bTK+iF6Y2GH4kzIxPY4mIj4dYvsH75t/dCaRnq1E7v+3EZsJR93YOu9BTfCdztayrcyYN4Px/cfz69N+He7iiEgX0OFCH3jYGDOGYNP9FuB7ANba1caYfwNrAB9wW7v13HfnA+CphZjQkLQaL1TVePDihphjcPXMJK7H0UQlZwWb2mP6QUx/Il0JRB7s/GN+1z61cDXBdxr+gJ9rXr8Gb8DLi5e8iMvZCrdgRKTb63Chb6299gD7ZgPt3x4dmwHuvLrAB4hyQVRSJly45sjPr1q4hGTnZjNr/izyyoL9L7439nscnXp0mEslIl1Fhwv9Dqk97oerFt7tZedmM3PeTNzeff/O/rnyn5yceTLTR+rfhogcOS3N1RIDp8OEOcF77Jjg44Q5CmlpVbPmz2oQ+AAen4dZ82eFqUQi0tWopt9SqolLG9vbpL+//LL8di6JiHRVqumLdABVtVVEOpvu8pmR1EpzTYhIt6fQFwkzr9/Lpf+5FK/f2yj4Y12xzJ6quRREpHUo9EXCKGAD3Pjmjby78V3+eu5fef6C58lMysRgyEzKZM55c9SJT0Raje7pi4SJtZYfvvdDsnOz+e3pv2XG2BkACnkRaTOq6YuEye8+/x1PLn6Suybexf0n3R/u4ohIN6DQFwmDZ5c+y6yPZnHNqGt49KxHtXKeiLQLhb5IO3vtm9e45e1b+M7R3+H585/HYfTfUETah37aiLSjjzd/zFWvXsXEtIn857L/aE59EWlXCn2RdrKscBkXvHwBR6cezVtXv0VcZFy4iyQi3YxCX6QdbCjewNn/OpuUmBTeu+Y9UmNSw10kEemGFPoibWx7xXbO/NeZWCzvX/M+AxIHhLtIItJNaZy+SBsq8ZRw1r/OYrd7Nx9f/zFDew4Nd5FEpBtT6Iu0EbfXzXkvncf64vW8ffXbjOs/LtxFEpFuTqEv0oqyc7OZNX8WeWV5xETEUO2r5pVLX2HaoGnhLpqIiEJfpLVk52Yzc95M3F43AB6fB5fDRW2gNswlExEJUkc+kVYya/6susDfyxvwMmv+rDCVSESkIYW+SCvJK8trcnt+WX47l0REpGkKfZFWsG73OiIcTd8ty0jKaOfSiIg0TaEvcoTe2fAOE/42gVhXLFHOqAb7Yl2xzJ46O0wlExFpSKEvcpistTz8xcOc++K5DEoZxNe3fM1zFzxHZlImBkNmUiZzzpvD9JHTw11UERFAvfdFDovH62HGvBlk52Zz+fDLef7854mLjCMzOVMhLyIdlkJf5BBtK9/GRa9cxNLtS5l9+mx+ctJPMMaEu1jtzlrLli1byMnJYc2aNQQCAWbNmkV0dHS4iyYizVDoixyCL7d+ycWvXEyVt4o3r3yT84aeF+4itZuSkhJycnJYtmwZbndwaGJWVhYTJ07k0ksvxel0hrmEInIwCn2RFnp++fPc+vatpCemM/+6+QzvPTzcRWozNTU1rFy5ksWLF1NUVARASkoK48eP54477iAuTssCi3RGCn2Rg/AFfPzovR/xVM5TTBs0jVcufaVLLY1rrWXDhg3k5OSwfv16rLVERkYyZswYLr/8cvr06RPuIopIK1HoixxAsbuYK/57BfM3z+euiXfxhzP/0Ox4/M5i165dLF68mBUrVlBTU4MxhqOPPpqJEydy9dVX43BoUI9IV9W5f3qJtLK9C+bkl+XTN74vfuuntLqUv1/wd24Yc0O4i3fIPB4Py5YtIycnhz179gDQq1cvJkyYwD333KNOdyLdjEJfJGT/BXMKKwsB+MWpv+gUgR8IBFi7di2LFy9m06ZNAMTExHD88cdz3XXX0aNHjzCXUETCTaEvEtLUgjkA/1jxDx6c8mD7F+ggCgsLWbx4MV9//TVerxeHw8ExxxzDKaecwg033NAthxGKyIEp9EUIdtbryAvmVFZWsnTpUr766ivKysoA6NevHxMnTuScc84hMjIyzCUUkc5AoS/dWsAG+M/q//Czj3/W7DHtvWCOz+dj9erV5OTkkJ8f/IUjPj6esWPHMmPGDJKSktq1PCLSdSj0pVuy1vL2hrf56Uc/ZeXOlYzoPYK7T7ibvyz5Cx6fp+64tl4wx1rL1q1bWbx4MatXr8bv9+N0Ohk+fDhnnXUW6enpaqYXkVaj0JduZ8GWBTww/wEWblvI0alH8+LFL3LFiCtwGAdj+4+t672fkZTB7KmzW3Uu/dLSUpYsWcKSJUuoqqoCICMjgwkTJnDRRRcREaH/kiLSdvQTRrqNnIIcZn00iw83fUhaQhpzzp3DDWNuwOV01R0zfeT0Vgv52tpavv76a3JycigsDI4ESEpKYvz48dx+++3Ex8e3ynVERFpKoS9d3qqiVfzs45/xxto36Bnbk8fOfIxbx99KdETrjVG31rJp0yZycnJYu3Ytfr+fyMhIRo0axYUXXkj//v1b7VoiIodLoS9d1rd7vuXBTx4k++tsEqIS+PVpv+bOiXeSEJVwxOcuLi4mJyeH5cuX4/EE+wAcddRRTJgwgSuuuEKz2olIh6TQly6noLyAX3/6a55b/hwuh4sfT/4xP57848OeL7+6upoVK1awePFidu/eDUCPHj2YMGECd999NzExMa1ZfBGRNqPQly5jt3s3D33+EM989Qz+gJ9bxt7CAyc/QL+Efi0+RyAQYP369eTk5LBhwwYAoqOjGTNmDFdffTW9evVqq+KLiLQ5hb50euU15Ty28DEeXfgobq+b60Zfxy9O/QVZyVkHfe3OnTtZvHgxK1eupLa2FmMMQ4YMYdKkSVx77bUaLiciXYpCXzqN+ovhZCRl8PNTf06xu5iHvniIPZ49XDrsUn415Vcc2+vYJl9fVVVVt/hMaWkpAH369GHixImcddZZREVFteO7ERFpfwp96RT2XwwnryyP7879LhbLd47+Dr85/Tcc3+/4uuP9fj9r1qwhJyeHLVu2ABAXF8fxxx/PTTfdREpKSjjehohIWCn0pVN44MMHGi2GY7H0ievD21e/TUFBAa+++iqrVq3C5/PhdDoZNmwYU6dOJTMzU830IiIo9KUDq/HV8EneJ7y59k3yy+stehMAdgLfws7anfy87OcMGDCAiRMncv755+NyuZo7pYhIt6bQlw5lj2cP72x4h7nr5vLuxnepqK0g2hmNq9yFt8ALuwEH0B8YD5m9M/n1Xb8Oc6lFRDoHhb6E3cY9G5m7bi5z183l8/zP8Vs/8TaewXYwwxzDOPfYc9mVtIv7c+5v0MTf1ovhiIh0NQp9aXf+gJ+cghxeWfkKr695nXxPsOm+t+3Nd5K+w/Rx07lk0iW4Iho206f2SW3TxXBERLo6hb60iz0Ve3ju4+d4bfVr5NbkUmWqcODguJTj+MGkH3DJiEsYmDLwgOdozcVwRES6I4W+tIr6Y+jTE9L5weAfELU7ive2vMd61rPZsRkfPhIjEzl/5PmcP/R8zj76bJKjk8NddBGRbkOhL0fsTwv+xA+f/yE122ogGvJ75vPjnT/GxlpwQFZyFrcNvY3zh57PyRknN1jKVkRE2o9CXw6J2+1m2bJlLF68mE2lmyg0hbxT/A41aTUwCAgNh7dYkqOS+fTGTxnRe4TGyYuIdAAKfWmW3+9n7dq1LFq0iOVbllNoCtnl2kVFfAWbvZspiygLHti76deX1ZQxss/I9iuwiIgckEK/G9l/7vr9e78XFBSwePFiPsv9jG3+bRQ6CimPK2erbyulEaUAuHAxOn40Vw2+irH9xzKu/zgueOmChpPnhGQkZbTXWxMRkRZQ6HcTjeauL8rjpidv4t+J/yYyLpJCU0hJTAnb2U4ppeCECEcEIxNHcln/yxjXfxxj+41lRO8RREU0XJjmt9N+2+DcoDH0IiIdkUK/G/D5fNzzz3twb3KDBRKAZKjtV8vciLkAOI2T4SnDubjfxcGA7z+WUX1GER0RfdDz720t0Bh6EZGOLSyhb4y5DHgQOBaYYK1dUm/fT4CbAT9wh7X2vdD2scA/gBjgHeBOa61t35J3fNZaVm9czRtfvMGizYvYYXewy7GLHTE7YFTTr/nypi8Z3Xc0sa7Yw76uxtCLiHR84arprwIuBv5af6MxZhhwJTCc4OzqHxpjhlhr/cCfgZnAIoKhfzbwv/YsdEfjD/hZnrecNxe9yZfffsk27zaKTBGlpjR4gAMSoxIZ1WcUe3bsobK2stE5MpMymZQ+qX0LLiIiYRGW0LfWfgM0NYzrAuBla20NsNkYsxGYYIzZAiRaaxeGXvdP4EK6UOgfrJPdbvdulhUs473l77Fw00K2VG9hF7vwGR8ADuNgaO+hnNnnTEb1HsWoPsGvjKQMjDGN7umD7ruLiHQ3He2efhrBmvxe20LbvKHn+29vkjFmJsFWATIyOn4P8kad7MryuPGNG/nbF39jT8UeNldtpsJU1B3fI7oHYwaNqQv2UX1GcWzPY4lxxTR7Dd13FxGRNgt9Y8yHQN8mds2y1r7Z3Mua2GYPsL1J1to5wByAcePGddj7/ruqdrFy50p+8M4PcLvdwXXiXYADvAEvCwoXcEzqMVx49IWM7jO6LuD7xPc5rOvpvruISPfWZqFvrZ12GC/bBqTX+34AsD20fUAT2zsFf8DP+uL1rNy5kpU7VrKscBnLCpaxu2b3voMswdCv/yuKA76585t2Lq2IiHRVHa15fy7wojHmMYId+QYDOdZavzGmwhhzArAYuA54uj0LdrB77nuV15Tz9c6vWbFjBSt3rGTlzpXk7syl2l8NgMM66G16M6rHKE4deyqTj5rMDW/cwLaKbY3OlZmU2ebvS0REuo9wDdm7iGBo9wLeNsassNaeZa1dbYz5N7AG8AG3hXruA9zKviF7/6MdO/E1dc99xtwZFFUVMTB5ICt3rGTFzmDIby7dXPe6GBtDH/pwWsppnHbsaZwx8gyG9R5GpDOywfkfOuMhdbITEZE2Z7r6UPdx48bZJUuWHPzAA8h6Iou8srxm9xsMqTaVPrYPR8UfxenDTueCEy4gKzWrxQvNtLQlQaQzmDJlCgALFiwIazlEuiNjzFJr7bim9nW05v0OKb8sNK/83jaHCqAS8MBTlzzFxSddTFqvZgcTtIg62YmISFtT6LdARlJGsKbvBaqBJCA5eM/9Bxf9ILyFExERaSFHuAvQGcyeOjs4RW00kAwY3XMXEZHOR6HfAtNHTmfOeXPITMrEYMhMymTOeXPUHC8iIp2KmvdbSPfcRUSks1NNX0REpJtQ6IuIiHQTCn0REZFuQqEvIiLSTSj0RUREugmFvoiISDeh0BcREekmFPoiIiLdhEJfRESkm1Doi4iIdBMKfRERkW7CWGvDXYY2ZYzZBeSFuxztrCewO9yF6OT0GR45fYZHTp/hkeuOn2GmtbZXUzu6fOh3R8aYJdbaceEuR2emz/DI6TM8cvoMj5w+w4bUvC8iItJNKPRFRES6CYV+1zQn3AXoAvQZHjl9hkdOn+GR02dYj+7pi4iIdBOq6YuIiHQTCv0uwhjzB2PMWmPM18aY140xyfX2/cQYs9EYs84Yc1YYi9mhGWMuM8asNsYEjDHj9tunz7CFjDFnhz6njcaY+8Ndns7CGPO8MabIGLOq3rZUY8wHxpgNoceUcJaxozPGpBtjPjbGfBP6v3xnaLs+xxCFftfxATDCWjsKWA/8BMAYMwy4EhgOnA38yRjjDFspO7ZVwMXAp/U36jNsudDn8gzwHWAYcFXo85OD+wfBf1/13Q/Mt9YOBuaHvpfm+YAfWWuPBU4Abgv9+9PnGKLQ7yKste9ba32hbxcBA0LPLwBettbWWGs3AxuBCeEoY0dnrf3GWruuiV36DFtuArDRWrvJWlsLvEzw85ODsNZ+CuzZb/MFwAuh5y8AF7ZnmToba22htXZZ6HkF8A2Qhj7HOgr9rukm4H+h52nA1nr7toW2ScvpM2w5fVatq4+1thCCgQb0DnN5Og1jTBZwHLAYfY51IsJdAGk5Y8yHQN8mds2y1r4ZOmYWwSau7L0va+L4bjtkoyWfYVMva2Jbt/0MD0KflYSdMSYeeBW4y1pbbkxT/yy7J4V+J2KtnXag/caY64Fzgal231jMbUB6vcMGANvbpoQd38E+w2boM2w5fVata6cxpp+1ttAY0w8oCneBOjpjjItg4Gdba18LbdbnGKLm/S7CGHM2cB9wvrXWXW/XXOBKY0yUMWYgMBjICUcZOzF9hi33FTDYGDPQGBNJsAPk3DCXqTObC1wfen490FxrlAAmWKV/DvjGWvtYvV36HEM0OU8XYYzZCEQBxaFNi6y1t4T2zSJ4n99HsLnrf02fpXszxlwEPA30AkqBFdbas0L79Bm2kDHmHOAJwAk8b62dHd4SdQ7GmJeAKQRXhdsJ/AJ4A/g3kAHkA5dZa/fv7CchxpiTgM+AXCAQ2vwAwfv6+hxR6IuIiHQbat4XERHpJhT6IiIi3YRCX0REpJtQ6IuIiHQTCn0REZFuQqEv0okZYyr3+/4GY8wfw1WejsYYk2WM8RhjVhhjVhpjvjTGDA13uUTCRaEvIl3dt9baMdba0QQXW3kg3AUSCReFvkgXZYw5zxiz2Biz3BjzoTGmT2j7g8aYglDtd4UxpswYMyW0r9IY86gxZpkxZr4xpldo+4RQTXlF6LUPhrb/0xizdxKofxhjLjXGOIwxbxpjLgttP8oY864xZqkx5jNjzDH1j69X3lWhmnnW3jXljTEuY8ymva0XxphexphXjTFfhb4mH+LHkgiUHPaHKtLJae59kc4txhizot73qeyb9vZz4ARrrTXGfBf4MfCj0L7HrbWPABhj3qr3+jhgmbX2R8aYnxOcFe52glM8/9pa+19jzD1AfOj47wLvGmM21TvHo0COtfY/oe/nALdYazcYYyYCfwJOb+H7mwnUv4XxZKjsnxtjMoD3gGMPco6jQp9RAhALTGzhtUW6HIW+SOfmsdaO2fuNMeYGYFzo2wHAK6EFRiKBzS04XwB4JfT8X8DeBUv8BEOzAWttrTHm9wSnOP2G4AqGJwD9Q+WJB04E/lNvpbOoeqf4gzHmp6HnR9U/tzEmFrgR+DMwPLR5GjCs3rkSjTEJobXTm/Pt3s/IGHMFwV9Czj7A8SJdlpr3Rbqup4E/WmtHAt8Dog/jHHvn6X4QuCe0xsPde3caY5zArwiG8zEEWwoeAmaFDnEApaF76nu/6tfM7927Hfh2v2vfRTCgPfW2OYBJ9c6VdpDA399c4JRDOF6kS1Hoi3RdSUBB6Pn1BzqwHgew9z771QRvEQDsINjMfgrweL3jbwMWWGtfJ7hy2ZPAI8CZxpih1tpyYHO9+/vGGDO6hWW/EHh+v+3vE7zdQOh8Y0KPE4wx/2zBeU+i8S8XIt2GmvdFuq4HCTarFwCLgIEteE0VMNwYsxQoA64ILVf6D+ABa+32vU3rxpi+wAyCzfl1rLU+Y8xdwB+BM4DpwJ9Dzfgu4GVg5UHKMQC4J3Su+tvvAJ4xxnxN8OfXp8AtBFdP8zQ6S9Dee/oGqCXYD0GkW9IqeyJSxxhTaa2NP/iRHYsx5g/A/7PWfh3usoh0ZAp9EanTWUNfRFpGoS8iItJNqCOfiIhIN6HQFxER6SYU+iIiIt2EQl9ERKSbUOiLiIh0Ewp9ERGRbuL/A+qFUzngka9dAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "import matplotlib.pyplot as plt\n", "import numpy as np\n", "\n", "#5 mA\n", "#V\n", "voltage_5 = [25, 22, 19, 16, 13, 10, 8, 6.04, 3.95, 2.07, 0.48, -0.38, -2, -4, -6.04, -8, -10, -13, -16.07, -19, -22, -25]\n", "#muA\n", "current_5 = [120.02, 116.93, 113.78, 110.29, 105.35, 97.18, 89.20, 78.43, 63.32, 46.90, 25.17, 3.59, -18.41, -37.95, -55.38, -68.57, -78.62, -88.28, -94.05, -97.70, -101.01, -104.42]\n", "\n", "#3 mA\n", "#V\n", "voltage_3 = [25, 22, 19, 16, 13, 10, 8, 6.03, 4.05, 2, 0.29, -0.73, -2, -4, -6, -8, -10, -13, -16.07, -19, -22, -25]\n", "#muA\n", "current_3 = [69.05, 66.95, 64.84, 62.73, 60.27, 56.44, 52.25, 46.12, 37.51, 25.78, 14.42, 4.48, -4.80, -17.45, -27.78, -35.60, -40.79, -45.21, -47.58, -49.35, -51.09, -52.93]\n", "\n", "#1.5 mA\n", "#V\n", "voltage_15 = [25, 22, 19, 16, 13, 10, 8, 6, 4, 2, 0.36, -0.30, -2, -4, -6, -8, -10, -13, -16, -19, -22, -25]\n", "#muA\n", "current_15 = [36.20, 34.93, 33.65, 32.39, 31.05, 29.13, 27.02, 23.91, 19.54, 13.68, 7.97, 5.06, -1.01, -7.63, -12.91, -16.76, -19.39, -21.44, -22.38, -23.16, -23.95, -24.77]\n", "\n", "k_5_r, b_5_r = np.polyfit(voltage_5[:4], current_5[:4], deg = 1)\n", "k_5_l, b_5_l = np.polyfit(voltage_5[17:], current_5[17:], deg = 1)\n", "k_3_r, b_3_r = np.polyfit(voltage_3[:4], current_3[:4], deg = 1)\n", "k_3_l, b_3_l = np.polyfit(voltage_3[17:], current_3[17:], deg = 1)\n", "k_15_r, b_15_r = np.polyfit(voltage_15[:4], current_15[:4], deg = 1)\n", "k_15_l, b_15_l = np.polyfit(voltage_15[17:], current_15[17:], deg = 1)\n", "\n", "current_5_r = []\n", "current_5_l = []\n", "current_3_r = []\n", "current_3_l = []\n", "current_15_r = []\n", "current_15_l = []\n", "for i in range(11):\n", " current_5_r.append(k_5_r * voltage_5[i] + b_5_r) \n", " current_5_l.append(k_5_l * voltage_5[11 + i] + b_5_l)\n", " current_3_r.append(k_3_r * voltage_3[i] + b_3_r) \n", " current_3_l.append(k_3_l * voltage_3[11 + i] + b_3_l)\n", " current_15_r.append(k_15_r * voltage_15[i] + b_15_r) \n", " current_15_l.append(k_15_l * voltage_15[11 + i] + b_15_l)\n", "\n", "plt.figure(figsize=(8, 6))\n", "plt.scatter(voltage_5, current_5, color = 'green', marker = 'o')\n", "plt.scatter(voltage_3, current_3, color = 'orange', marker = 'o')\n", "plt.scatter(voltage_15, current_15, color = 'purple', marker = 'o')\n", "plt.plot(voltage_5[:11], current_5_r, color = 'black', linewidth = 0.5)\n", "plt.plot(voltage_5[11:], current_5_l, color = 'black', linewidth = 0.5)\n", "plt.plot(voltage_3[:11], current_3_r, color = 'black', linewidth = 0.5)\n", "plt.plot(voltage_3[11:], current_3_l, color = 'black', linewidth = 0.5)\n", "plt.plot(voltage_15[:11], current_15_r, color = 'black', linewidth = 0.5)\n", "plt.plot(voltage_15[11:], current_15_l, color = 'black', linewidth = 0.5)\n", "plt.plot(voltage_5, current_5, color = 'green')\n", "plt.plot(voltage_3, current_3, color = 'orange')\n", "plt.plot(voltage_15, current_15, color = 'purple')\n", "plt.axhline(y=0, color ='k')\n", "plt.axvline(x=0, color ='k')\n", "plt.xlabel('Напряжение, В')\n", "plt.ylabel('Ток, мкА')\n", "\n", "print((b_5_r - b_5_l) / 2, (b_3_r - b_3_l) / 2, (b_15_r - b_15_l) / 2)\n" ] }, { "cell_type": "code", "execution_count": null, "id": "e4a91142", "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.9.12" } }, "nbformat": 4, "nbformat_minor": 5 }