{ "cells": [ { "cell_type": "code", "execution_count": 26, "metadata": { "scrolled": true }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAEGCAYAAABsLkJ6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO3dfZRV1X3/8fd3BgbkSZBBDDMgoKgBNaAjaEyxiWmC+SXgSrSR2qqrZrFiYruSmLTaB/MrNWlMuhrbxJWG+JBoomhME/klGIxR81CNgkpQYtEBEQbQeAEnKA/DMN/fH+fccr3cmTn3+WF/Xmuddc89Z5999p4L93v33uecbe6OiIiEp6naBRARkepQABARCZQCgIhIoBQAREQCpQAgIhKoIdUuQD5aW1t96tSpBR3b29vLkCF1Vd2iqc5hUJ0bX7H1feqpp1LuPiF7e139BadOncqaNWsKOjaVStHa2lriEtU21TkMqnPjK7a+ZvZyru3qAhIRCZQCgIhIoBQAREQCVVdjACIildTT08PGjRvZu3dvtYvCyy/n7MZ/ixEjRnDCCSfQ0tKSKE8FABGRfmzcuJGxY8dy8skn09RU2x0mfX19vPrqq3R2djJz5sxExySqkZktMLMNZtZpZtfm2P8ZM/udma0zs5+b2fEZ+y43sxfj5fKM7Wea2bNxnv9hZpaoxCIiFbJ3714mTpxY81/+AE1NTUycOJG9e/fyyiuvJDtmsARm1gzcDFwAzAQWm1l2eHkG6HD304H7gC/Hxx4DfB6YB8wFPm9m4+JjvgEsAWbEy4JEJRYRqaB6+PJPa2pqwsy4//77OXTo0ODpE+Q5F+h0903u3gMsBxZlJnD3R9w93Un2G6A9Xn8/8DN33+Xuu4GfAQvM7G3AGHd/3KPnUd8BXJikgoUa8sQT8Mwz5TyFiEhN2L9/P/v27Rs0XZIxgDZga8b7LqJf9P25EnhggGPb4qUrx/YjmNkSopYC7e3tpFKpBEU+UusHPwhA6rXXCjq+HnV3d1e7CBWnOochpDq/853v5LHHHsvrmAMHDrBr1y72798/YLokASBX33zOWWTM7M+BDuC8QY5NnKe7LwOWAXR0dHixd/+1jhsHzc1F5VFPQrpbMk11DkMl6pzkyptyy/fLH2D48OEcc8wxjBo1asB0SbqAuoDJGe/bge3ZiczsvcDfAwvd/cAgx3ZxuJuo3zzLosAWhIhINaS/xN944w3OP/98zjjjDE477TTuv//+ovNO0gJYDcwws2nANuAS4M8yE5jZHOCbwAJ3/33GrlXAFzMGft8HXOfuu8xsj5mdDTwBXAZ8rbiqDKCv7/D6jh0wcWLZTiUiDepTn4K1a0ub5+zZcNNNiZIOHz6cH/7wh4wZM4ZUKsXZZ5/NwoULKeYCykEDgLv3mtnVRF/mzcBt7r7ezJYCa9x9BfAVYBTw/bgwW9x9YfxF/89EQQRgqbvvitevAr4NHEU0ZvAA5bJnz+H1HTuiP7qISB1xd/7u7/6OX/7ylzQ1NbFt2zZeffVVjjvuuILzTHQjmLuvBFZmbbs+Y/29Axx7G3Bbju1rgFMTl7QYb755eH3HjoqcUkQaTMJf6uXyve99j9dee42nnnqKoUOHMnXq1EEHeQdTPxe4FiPzelgFABGpQ93d3Rx77LEMHTqURx55pCQD1GE8CiJ7DEBEpM5ceumlfOhDH6Kjo4PZs2dzyimnFJ2nAoCISA174403gOiy18cff7ykeasLSEQkUGEEgHQLwEwBQEQkFlYAaGuLAoDnvOlYROQIfZldyDUu37KGEQDSXUDt7XDgALz+enXLIyJ1YcSIEezYsaMugkBfXx+vvPIKBw8eTHxMWIPAbfHz5nbsgHHj+k8vIgKccMIJPP3004mfr19tBw8eZMuWLbh7osdYhxsAEs6YIyLhamlpobe3l8cee4yxY8cW9diFYuzbt4+jjjoqUdo9e/YwceLEROnDCACZXUCggWARSeycc87BzNi0aVOiSVbKwcwSBR8zY8aMGZx33nmJ0ocRAHK1AEREEmhububcc8/l3HPPrVoZUqlUWR5/HcYgcDoAHH00jBwJ2yvz5GkRkVoWRgBIN9uam2HSJAUAERFCCQDpFkBzc9QNtG1bdcsjIlIDwgoATU0KACIisTACQLoLKB0Atm/X3cAiErxEAcDMFpjZBjPrNLNrc+yfb2ZPm1mvmV2Usf3dZrY2Y9lvZhfG+75tZi9l7CvfNF2ZXUCTJkFPD+zcWbbTiYjUg0EvAzWzZuBm4E+IJnNfbWYr3P13Gcm2AFcAn8081t0fAWbH+RwDdAIPZiT5nLvfV0wFEsnuAoKoG6gMl1WJiNSLJC2AuUCnu29y9x5gObAoM4G7b3b3dcBAD8y4CHjA3fcWXNpCZXcBgcYBRCR4SW4EawO2ZrzvAuYVcK5LgH/L2vYFM7se+DlwrbsfyD7IzJYASwDa29tJpVJ5n3jo7t0cDby+Zw99EydyDLBnwwYOzJ2bd171pLu7u9pFqDjVOQyh1blc9U0SAHLdT5zXCKqZvQ04DViVsfk64BWgBVgG/C2w9IgTuS+L99PR0eEF3Q03ejQAY485BmbNijZ1dzM6gC6gctw9WOtU5zCEVudq3QncBUzOeN8O5Hsn1Z8CP3T3/31Oqbvv8MgB4HairqbyyOwCGjoUjj1WN4OJSPCSBIDVwAwzm2ZmLURdOSvyPM9i4O7MDXGrAIueWHQh8FyeeSaXeRUQ6F4AERESBAB37wWuJuq+eR64193Xm9lSM1sIYGZnmVkXcDHwTTNbnz7ezKYStSB+kZX198zsWeBZoBW4ofjq9CPzKiBQABARIeHTQN19JbAya9v1GeuribqGch27mWggOXv7e/IpaFEyu4AgCgC/+U3FTi8iUovCuBM4uwto0iRIpaLpIUVEAhVWAMhsAYDmBRCRoIURAHJ1AYHGAUQkaGEEgFxXAYECgIgELawAoBaAiMj/CiMAZHcBjRsHw4YpAIhI0MIIANldQGa6F0BEghdWALCMxxpNmQJbt+ZOLyISgDACQHr2r6aM6k6ZAlu2VKc8IiI1IKwAkGnKlKgLqLe38uUREakBYQSAtOwuoL4+PRVURIIVRgDI1QI4/vjo9eWXK1sWEZEaEUYASMtuAYDGAUQkWGEEgFwtgMnxHDcKACISqDACQFpmC2DkSBg/XgFARIIVRgDI1QIAXQoqIkFLFADMbIGZbTCzTjO7Nsf++Wb2tJn1mtlFWfsOmdnaeFmRsX2amT1hZi+a2T3xdJPlZVnz2ysAiEjABg0AZtYM3AxcAMwEFpvZzKxkW4ArgLtyZLHP3WfHy8KM7TcCX3X3GcBu4MoCyp9Mfy2A44+PrgLqb7+ISANL0gKYC3S6+yZ37wGWA4syE7j7ZndfB/QlOWk8Efx7gPviTd8hmhi+vHK1APbsge7usp9aRKTWJJkTuA3IfGhOFzAvj3MMN7M1QC/wJXf/ETAeeD2ecD6d5xHzBgOY2RJgCUB7ezupVCqPU8cF2LOHUcDOnTvxjF/7LWPHMgbY/dvfcmjWrLzzrXXdAQY21TkModW5XPVNEgAsx7Z8+kymuPt2M5sOPGxmzwJ/SJqnuy8DlgF0dHR4a2trHqeOjRoFwPjWVjjmmMPbTz0VgHF79kAh+daBgv5edU51DkNodS5HfZN0AXUBkzPetwOJn5/g7tvj103Ao8AcIAWMNbN0AMorz7wNdBUQaCBYRIKUJACsBmbEV+20AJcAKwY5BgAzG2dmw+L1VuBc4Hce9cM8AqSvGLocuD/fwuctewxg4kRoaVEAEJEgDRoA4n76q4FVwPPAve6+3syWmtlCADM7y8y6gIuBb5rZ+vjwtwNrzOy3RF/4X3L338X7/hb4jJl1Eo0J3FrKiiXS1BTdEawAICIBSjIGgLuvBFZmbbs+Y301UTdO9nGPAaf1k+cmoiuMym+gyzynTNED4UQkSGHcCZyW3QUEMG0avPRS5csiIlJlYQSAgVoA06fDjh2wd2/lyiMiUgPCCABpuVoA06dHr5s3V7QoIiLVFkYAGKgFMG1a9LppU2XKIiJSI8IIAGkDtQAUAEQkMGEEgIFaABMmRHMDKACISGDCCABpuVoAZlErQAFARAITRgAY7HHPCgAiEqAwAkBarhYAHA4AmhdARAISRgBI0gLYtw9efbUy5RERqQFhBIC0gVoAoG4gEQlKGAEgSQsAFABEJChhBIC0/loAU6dGrwoAIhKQMALAYC2A4cOhrU0BQESCEkYASOuvBQC6FFREgpMoAJjZAjPbYGadZnZtjv3zzexpM+s1s4syts82s8fNbL2ZrTOzj2bs+7aZvWRma+NldmmqlEOSyzunT4eNG8tWBBGRWjNoADCzZuBm4AJgJrDYzGZmJdsCXAHclbV9L3CZu88CFgA3mdnYjP2fc/fZ8bK2wDokN1AL4MQTYft2ePPNshdDRKQWJGkBzAU63X2Tu/cAy4FFmQncfbO7rwP6sra/4O4vxuvbgd8DE0pS8nwkaQGcfHL0+uKL5S2LiEiNSDIlZBuwNeN9FzAv3xOZ2VygBcjsZ/mCmV0P/By41t0P5DhuCbAEoL29nVQqle+pOerNNxkJpHbuhGHDcqZpnjCBccAf1qyhp/2I2S3rUnd3d7WLUHGqcxhCq3O56pskAOTqN8nrmQlm9jbgTuByd0+3Eq4DXiEKCsuIJolfesSJ3JfF++no6PDW1tZ8Th0ZMQKA1tbWfgMAZ50FwJgdO6CQc9Sogv5edU51DkNodS5HfZN0AXUBkzPetwPbk57AzMYAPwH+wd1/k97u7js8cgC4nUpMED/QGMDIkTB5MrzwQtmLISJSC5IEgNXADDObZmYtwCXAiiSZx+l/CNzh7t/P2ve2+NWAC4Hn8il4XpI+5O2kkxQARCQYgwYAd+8FrgZWAc8D97r7ejNbamYLAczsLDPrAi4Gvmlm6+PD/xSYD1yR43LP75nZs8CzQCtwQ0lrlstALQCIAsCGDXoqqIgEIckYAO6+EliZte36jPXVRF1D2cd9F/huP3m+J6+SFiPpF/rJJ0N3N7z2Ghx7bHnLJCJSZboTONNJJ0Wv6gYSkQCEEQDyGQOAqBtIRKTBhREA0gZrAUydCkOHqgUgIkEIIwAkbQE0N0ePhFAAEJEAhBEA0gZrAcDhK4FERBpcGAEgn8s6Tz4ZOjuht7d85RERqQFhBIC0JC2AU06Bgwc1N4CINLwwAkA+LYBZs6LX9esHTiciUufCCABpSVoAM+OpDhQARKTBhREA8mkBjBoFxx+vACAiDS+MAJCvWbMUAESk4YURAPJ9uNupp0aXgupKIBFpYGEEAMCT9P+nzZoFPT3R5aAiIg0qmACQl/SVQM+Vb4oCEZFqCyMA5NsF9Pa3R1cMaRxARBpYGAEAkl0CmjZiBEybpgAgIg0tUQAwswVmtsHMOs3s2hz755vZ02bWa2YXZe273MxejJfLM7afaWbPxnn+Rzw1ZHkUMsPXqacqAIhIQxs0AJhZM3AzcAEwE1hsZjOzkm0BrgDuyjr2GODzwDyiSd8/b2bj4t3fAJYAM+JlQcG1SCLf+DJrVvRU0J6e8pRHRKTKkrQA5gKd7r7J3XuA5cCizATuvtnd1wF9Wce+H/iZu+9y993Az4AF8YTwY9z9cXd34A6iieHLo5AWwKxZ0WWgejS0iDSoJHMCtwFbM953Ef2iTyLXsW3x0pVj+xHMbAlRS4H29nZSqVTCUx82Yu9ejjLL69jmKVMYB+z51a84cNxxeZ+zFnR3d1e7CBWnOochtDqXq75JAkCuvpOkP6n7OzZxnu6+DFgG0NHR4a2trQlPneGoo3Agr2PPOQeGD2d0ZyejCzlnjSjo71XnVOcwhFbnctQ3SRdQFzA54307sD1h/v0d2xWvF5JnYfIdAxgyBE47DdauLU95RESqLEkAWA3MMLNpZtYCXAKsSJj/KuB9ZjYuHvx9H7DK3XcAe8zs7Pjqn8uA+wsofzKFjAEAzJ4NzzxT+PEiIjVs0ADg7r3A1URf5s8D97r7ejNbamYLAczsLDPrAi4Gvmlm6+NjdwH/TBREVgNL420AVwG3AJ3ARuCBktYsWyFXmc6ZA7t3w5YtpS+PiEiVJRkDwN1XAiuztl2fsb6at3bpZKa7Dbgtx/Y1wKn5FLZghf6CnzMnel27NnpEtIhIA9GdwAM5/XRoaoq6gUREGkwYAaDQFsCIEdEk8QoAItKAwggAUFgLAA4PBIuINJgwAkAxV/HMmQNbt8LOnaUrj4hIDQgjAEDhLYD0QLBaASLSYMIIAMW0AM48M3pdvbo0ZRERqRFhBAAovAUwbhycdBI88URpyyMiUmVhBIBi7+SdOzcKALojWEQaSBgBoFjz5sErr0BX1+BpRUTqRBgBoNhf7vPip1+rG0hEGkgYAQDwYmacfMc7YNgwBQARaShhBIBiWwAtLdHloAoAItJAwggAUPhVQGnz5sFTT0XTRIqINIAwAkAprt6ZOxf27oXnnis+LxGRGhBGAIDStABA3UAi0jDCCAClaAFMnw4TJ8Kvf118XiIiNSBRADCzBWa2wcw6zezaHPuHmdk98f4nzGxqvP1SM1ubsfSZ2ex436Nxnul9x5ayYjkqUfzx8+fDL39ZmvKIiFTZoAHAzJqBm4ELgJnAYjObmZXsSmC3u58IfBW4EcDdv+fus919NvAXwGZ3z5xl/dL0fnf/fQnqk1up7uCdPz+aHvLll0uTn4hIFSVpAcwFOt19k7v3AMuBRVlpFgHfidfvA86PJ3vPtBi4u5jCFqXYFgBEAQDUChCRhpBkTuA2YGvG+y5gXn9p3L3XzLqB8UAqI81HOTJw3G5mh4AfADe4H/lT3cyWAEsA2tvbSaVS2UkGNXLvXlrcCzr2LY47jmOOPpqeBx/kjQsuKC6vCuju7q52ESpOdQ5DaHUuV32TBIBcP52zv6gHTGNm84C97p55DeWl7r7NzEYTBYC/AO44IhP3ZcAygI6ODm9tbU1Q5CxHHUVfUxMFHZtt/nyGP/kkw0uRVwWUpM51RnUOQ2h1Lkd9k3QBdQGTM963A9v7S2NmQ4CjgV0Z+y8hq/vH3bfFr3uAu4i6msqjlE/xnD8fXnghejiciEgdSxIAVgMzzGyambUQfZmvyEqzArg8Xr8IeDjdnWNmTcDFRGMHxNuGmFlrvD4U+CBQ3jusSjEGAIfHAX71q9LkJyJSJYMGAHfvBa4GVgHPA/e6+3ozW2pmC+NktwLjzawT+AyQeanofKDL3TdlbBsGrDKzdcBaYBvwraJr038lSpfXnDkwciQ8+mjp8hQRqYIkYwC4+0pgZda26zPW9xP9ys917KPA2Vnb3gTOzLOsxSlVC2DoUDjvPHjoodLkJyJSJboTuBDve180DrB5c2nzFRGpoDACAJSuBQBRAAD42c9Kl6eISIWFEwBK6ZRToK0NHnyw2iURESlYGAGg1F1AZlEr4Oc/h0OHSpu3iEiFhBEAoLRdQBAFgN27o0liRETqUBgBoNQtAID3vjcKKuoGEpE6FUYAgNK3AFpb4Ywz4IEHSpuviEiFhBEAytECAPjQh+Dxx+H35XuStYhIuYQRAMpl4cIouPz4x9UuiYhI3sIIAOVqAcyeDZMnw4rsRyOJiNS+MAIAlH4MIJ3nwoXRQPDevaXPX0SkjMIIAOVqAQAsWgT79kX3BIiI1JEwAgDg5WgBQPRguDFj4P77y5O/iEiZhBEAytkCaGmBD3wgCgC9veU7j4hIiYURAKA8YwBpH/0opFLqBhKRupIoAJjZAjPbYGadZnZtjv3DzOyeeP8TZjY13j7VzPaZ2dp4+c+MY840s2fjY/7DrIzf0OVsAQAsWBB1Ay1fPnhaEZEaMWgAMLNm4GbgAmAmsNjMZmYluxLY7e4nAl8FbszYt9HdZ8fLxzO2fwNYAsyIlwWFVyOBcrYAhg+HD38Y/uu/YP/+8p1HRKSEkrQA5gKd7r7J3XuI5vZdlJVmEfCdeP0+4PyBftGb2duAMe7+eDx38B3AhXmXPqlytwAAFi+GP/wBfvrT8p9LRKQEkgSANmBrxvuueFvONPEcwt3A+HjfNDN7xsx+YWZ/lJG+a5A8S6ucLQCA97wHJkyAu+8u73lEREokyZzAub45s39S95dmBzDF3Xea2ZnAj8xsVsI8o4zNlhB1FdHe3k4qlUpQ5LcatW8fQ/r6Cjo2HyM/+EGG3303uzZtwseMKeu5kuju7q52ESpOdQ5DaHUuV32TBIAuYHLG+3Zgez9pusxsCHA0sCvu3jkA4O5PmdlG4KQ4ffsgeRIftwxYBtDR0eGtra0Jipxl+HAONTdT0LH5+MQn4PbbGf/gg/Dxjw+evgLKXucapDqHIbQ6l6O+SbqAVgMzzGyambUAlwDZD79ZAVwer18EPOzubmYT4kFkzGw60WDvJnffAewxs7PjsYLLgPLdSVWJMQCAM8+Ed7wDbrmlMucTESnCoAEg7tO/GlgFPA/c6+7rzWypmS2Mk90KjDezTuAzQPpS0fnAOjP7LdHg8MfdfVe87yrgFqAT2AiU98H65R4DSJ/jYx+LZgl75pnyn09EpAhJuoBw95XAyqxt12es7wcuznHcD4Af9JPnGuDUfApbsEq1AAAuvRQ++1m49Vb4+tcrd14RkTzpTuBSGzcOPvIR+O539YRQEalpYQSASrYAIBoA7u6OgoCISI0KIwBA5VoAAO96VzRf8E03VT74iIgkFEYAqPSXsBl86lPw/POwalVlzy0iklAYAQAq2wKA6Amhxx0XtQJERGpQGAGgGt0wLS3wyU9GLYDnnqv8+UVEBhFGAIDKtwAArroKRo2CG26o/LlFRAYRRgCo1kDs+PFw9dVw773ReICISA0JIwBU0zXXwIgRagWISM0JIwBU81LM1tboIXHLl8P//E/1yiEikiWMAADVGQNI++xno1bAdddVrwwiIlnCCADVvhnr2GPh2mvhRz+CX/yiumUREYmFEQCgui0AgE9/GtrbozGBvr7qlkVEhFACQLVbABB1AX3xi9GjovWMIBGpAWEEAMCr3QKA6FHR8+ZFrYCdO6tdGhEJXBgBoBZaAABNTbBsGbz+Onzuc9UujYgELlEAMLMFZrbBzDrN7Noc+4eZ2T3x/ifMbGq8/U/M7CkzezZ+fU/GMY/Gea6Nl2NLVal+KlHW7BM7/fSoBXD77fDII9UujYgEbNAAEM/pezNwATATWGxmM7OSXQnsdvcTga8CN8bbU8CH3P00ojmD78w67lJ3nx0vvy+iHgOrlRZA2vXXw/Tp8Jd/Gc0bICJSBUlaAHOBTnff5O49wHJgUVaaRcB34vX7gPPNzNz9GXffHm9fDww3s2GlKHjeaqUFANGA8J13wpYt8Fd/Ve3SiEigkswJ3AZszXjfBczrL42795pZNzCeqAWQ9hHgGXc/kLHtdjM7RDRv8A3uR/5UN7MlwBKA9vZ2UqlUdpJBjT5wAOvrK+jYsjnpJEZccw0jvvIV/vCud9Hz4Q+X/BTdAbYuVOcwhFbnctU3SQDI9dM5+4t6wDRmNouoW+h9GfsvdfdtZjaaKAD8BXDHEZm4LwOWAXR0dHhra2uCImdpaaG3qYmCji2nL34Rfv1rxlxzTTSL2MzsnrXi1VydK0B1DkNodS5HfZN0AXUBkzPetwPb+0tjZkOAo4Fd8ft24IfAZe6+MX2Au2+LX/cAdxF1NZVPLXUBpQ0ZAt//PowcCRdeGF0dJCJSIUkCwGpghplNM7MW4BJgRVaaFUSDvAAXAQ+7u5vZWOAnwHXu/t/pxGY2xMxa4/WhwAeB8s2aUmuDwJna2qIg8NJL8Gd/Br291S6RiARi0ADg7r3A1cAq4HngXndfb2ZLzWxhnOxWYLyZdQKfAdKXil4NnAj8Y9blnsOAVWa2DlgLbAO+VcqKHaEWWwBpf/RHcPPN8MADsGRJbQcsEWkYScYAcPeVwMqsbddnrO8HLs5x3A1Afw/CPzN5MYtUD1+oS5bA9u3wT/8UPTzuS1+qdolEpMElCgANoZZbAGmf/zy8+irceGN0qeg//mN9lFtE6lIYAaAeWgAQfdl//euwb18UDPbuhX/5FwUBESmLMAIA1M+XaHMz3HZb1AK48cboTuGvfS26YkhEpITC+FaplxZAWlNTNCg8ejR8+cuwaRPccw+MHVvtkolIAwnjaaBQPy2ANLOoBXDLLfDww3DOOfBc+a6UFZHwhBEA6q0FkOnKK+Ghh2DXLjjrLPjGN+q7PiJSM8IIAFB/LYBM550H69bBH/8xfOITsGgRdHVVu1QiUufCCACN8It54kT4yU/g3/4tahG8/e1w0026c1hEChZGAGgUTU3R5PLr10d3D3/60zB7Ntx/f2MEORGpqDACQKN9OU6bFrUGfvADOHgwepDcO98JDz7YeHUVkbIJIwBAfY8B5GIGH/5w1Bq45ZZoTOD974+mnLztNti/v9olFJEaF0YAaORfxUOGRFcKdXbCt78ddRNdeSVMmsTIv/kbePLJxq6/iBQsjAAAjdcCyDZsGFx+OaxdGw0SL1jA8LvvhnnzogHj666Dxx+HQ4eqXVIRqRFhBICQfgGbwfnnw113sWv9evjWt6I5B/71X6NxgkmT4LLL4NZbo1ZDSH8bEXmLMB4FAY3fAsjBx4yBj30sWnbvhp/+FFasiF7vvDNKNGlSFBjOOOPwMmFCdQsuIhURRgDQr1wYNw4WL44Wd3j+efjFL6LlySfhvvsOp21rg1NOgZNOipYZM6LXKVOiriYRaQiJAoCZLQD+HWgGbnH3L2XtH0Y0ofuZwE7go+6+Od53HXAlcAj4a3dflSTPUvMAWwD9MosmoJ85E666Ktq2e3c0fvDMM9HrCy/A3XcfOU/xhAnQ3n54aWuLto0fHy2trYfXW1oqXzcRSWzQAGBmzcDNwJ8QTf6+2sxWuPvvMpJdCex29xPN7BLgRuCjZjaTaA7hWcAk4CEzOyk+ZrA8S0ctgMGNGwfvfne0pLnDzp1RMHjhBdi6FbZtiy453bIFHnss2t+fUaOiJ5qml/T7zO0jR0atiuHDo9fMJde2IUOipbk5WrLW7fXXc26nqSnIbkCRgSRpAcwFOt19E501+sMAAAYKSURBVICZLQcWAZlf1ouA/xuv3wd83cws3r7c3Q8AL8VzBs+N0w2WZ2npP3/+zKJf9K2t0ThBLvv3R0Fg505Ipd76unMn7NkTLW+8Eb1u2/bW92++WdIijx9oZ1PT4UBQiiVpXkklTZuVbtyhQ1GgK2GeJUlbjjzjtGMHqnMDGnvoUDRn+PTpJc03SQBoA7ZmvO8C5vWXxt17zayb6P9iG/CbrGPb4vXB8gTAzJYASwDa29tJpVIJivxWR82ZQ097O4cKOLaedXd3V+ZEw4ZFg8mTJuV/rDv09GA9PXDgwICvduBAdBlrvNihQ9GzkA4dgr4+rLeXfW+8wVFDh2J9fUfsS6/jfrhVmF53xzLW+12yj8nML1fafP4OBaY72NPD0FzdbUXk2R8rQ56FpD148CBDhw5NflydO3jwID1vvklfib/DkgSAXKE5+xPrL01/23NdfprzX4G7LwOWAXR0dHhra2v/Je3PDTewL5WioGPrXGh13p9KMSqwOr+RSjE6sDrvCexz3lOm768k9wF0AZMz3rcD2/tLY2ZDgKOBXQMcmyRPEREpoyQBYDUww8ymmVkL0aDuiqw0K4DL4/WLgIfd3ePtl5jZMDObBswAnkyYp4iIlNGgXUBxn/7VwCqiSzZvc/f1ZrYUWOPuK4BbgTvjQd5dRF/oxOnuJRrc7QU+6e6HAHLlWfrqiYhIfxLdB+DuK4GVWduuz1jfD1zcz7FfAL6QJE8REamcMJ4FJCIiR1AAEBEJlAKAiEigFABERAJlXkfPyTGz14CXCzy8FQjrVmDVORSqc+Mrtr7Hu/sRz3mvqwBQDDNb4+4d1S5HJanOYVCdG1+56qsuIBGRQCkAiIgEKqQAsKzaBagC1TkMqnPjK0t9gxkDEBGRtwqpBSAiIhkUAEREAtVwAcDMFpjZBjPrNLNrc+wfZmb3xPufMLOplS9laSWo83wze9rMes3somqUsZQS1PczZvY7M1tnZj83s+OrUc5SSlDnj5vZs2a21sx+Hc/HXdcGq3NGuovMzM2s7i8LTfA5X2Fmr8Wf81oz+1hRJ3T3hlmIHi29EZgOtAC/BWZmpfkE8J/x+iXAPdUudwXqPBU4HbgDuKjaZa5Afd8NjIjXrwrkMx6Tsb4Q+Gm1y13uOsfpRgO/JJp6tqPa5a7A53wF8PVSnbPRWgD/O4G9u/cA6cnmMy0CvhOv3wecH09gX68GrbO7b3b3dUBfNQpYYknq+4i7743f/oZoxrl6lqTOf8h4O5J+plitI0n+LwP8M/BlYH8lC1cmSetcMo0WAHJNYN/WXxp37wXSE9jXqyR1biT51vdK4IGylqj8EtXZzD5pZhuJvhD/ukJlK5dB62xmc4DJ7v7jShasjJL+2/5I3L15n5lNzrE/sUYLAMVMYF+vGq0+g0lcXzP7c6AD+EpZS1R+iers7je7+wnA3wL/UPZSldeAdTazJuCrwDUVK1H5Jfmc/x8w1d1PBx7icG9GQRotABQzgX29SlLnRpKovmb2XuDvgYXufqBCZSuXfD/j5cCFZS1R+Q1W59HAqcCjZrYZOBtYUecDwYN+zu6+M+Pf87eAM4s5YaMFgGImsK9XSercSAatb9w18E2iL//fV6GMpZakzjMy3v4f4MUKlq8cBqyzu3e7e6u7T3X3qURjPQvdfU11ilsSST7nt2W8XQg8X9QZqz3yXYaR9A8ALxCNpv99vG0p0T8OgOHA94FO4ElgerXLXIE6n0X06+JNYCewvtplLnN9HwJeBdbGy4pql7kCdf53YH1c30eAWdUuc7nrnJX2Uer8KqCEn/O/xJ/zb+PP+ZRizqdHQYiIBKrRuoBERCQhBQARkUApAIiIBEoBQEQkUAoAIiKBUgAQKZKZjTWzT1S7HCL5UgAQKd5YoqfMitQVBQCR4n0JOCF+Pnu9P3dIAqIbwUSKFE8q9GN3P7XKRRHJi1oAIiKBUgAQEQmUAoBI8fYQPZ5YpK4oAIgUyd13Av9tZs9pEFjqiQaBRUQCpRaAiEigFABERAKlACAiEigFABGRQCkAiIgESgFARCRQCgAiIoH6/4QUVYtptVNsAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXAAAAEGCAYAAAB8Ys7jAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAeKUlEQVR4nO3de3RV5bnv8e+TBIhcDEgQkSABCtQgCJqCPSq2BRW1QGsZVlsveIoZp3aL2u3Yvdiq7A7tZTuK1jJsc1qox22tFFuKF7ZIhWKqggYohFuAIJeAwQBJgAC5veePBAwxkJW15lozc83fZ4wM1mWuOZ+Hhb+8zttrzjlERCR4UvwuQEREoqMAFxEJKAW4iEhAKcBFRAJKAS4iElBpidxYZmamy87OjuqzdXV1pKUltFzfqedwUM/hEEvPhYWF5c65Pi1fT+jfYHZ2Nh988EFUny0vLyczM9Pjijo29RwO6jkcYunZzHa29rp2oYiIBFSbAW5mc81sv5kVtXj9PjPbYmYbzOwX8StRRERaE8kI/A/ApOYvmNkXganAKOfcCOBJ70sTEZGzaTPAnXMrgIMtXv428DPn3ImmZfbHoTYRETmLaA9iDgOuNrPHgePAQ86591tb0MzygDyArKwsysvLo9pgZWVllKUGl3oOB/UcDvHoOdoATwN6AVcAnwPmm9lg18qdsZxz+UA+QG5urovlyHPYjlqDeg4L9RwOXvcc7Vkoe4C/uEargAYgfN+GiIiPoh2BLwS+BCw3s2FAZyC6fSOSdBoaGqipqeHEiRMcP36cEydOcOLECerq6qirq6O+vv7Un80f19XVcfDgQbp3797qcvX19Tjn4v7T0NBAe26zHOmyZ1quurqarl27RrVOL7afqHU211rPye6mm27yfATeZoCb2YvAF4BMM9sDPArMBeY2nVpYA9zV2u4TCYba2loOHjzIgQMHKC8v58CBAxw4cICqqioOHz7MkSNHPvXnyZ+TAd08qGtra/1uSRLIzPwuIRBGjBjBuHHjPF1nmwHunLvtDG/d7mklEhdHjhxh69at7Ny5kz179lBaWsqePXtO/Xz88cdtHlw555xz6NGjB927dz/153nnncdFF11Eeno6Xbp0afWn5XudOnUiNTWV1NRU0tLSTvvz5OPDhw+TmZn5qfdO/phZXH9SUlJOex6pSJdtbbkzXaGXqO37IaxXYnotXDcjSGKHDh1i7dq1rFmzhs2bN1NcXExxcTH79u07bblOnTrRv39/srKyuPzyyzn//PPJzMykd+/e9O7d+7THGRkZdOvWLaH3rAjjf9gnfzmJtJcCPIDq6upYs2YN//jHP3jvvfdYvXo1O3bsOPV+ZmYmw4YN4/rrr2fYsGEMHTqU7OxsBgwYQJ8+fUhJ0R0URJKBAjwgSkpKeOWVV3jjjTcoKCjg8OHDAAwZMoTc3Fzy8vK47LLLGDNmDH36fOqmZSKShBTgHdiWLVuYPXs2ixYtoqio8VY0n/3sZ7n99tu55pprGD9+PP369fO5ShHxiwK8gzl06BAvvPACzz33HB988AGpqamMHz+e2bNnM3nyZIYMGeJ3iSLSQSjAO4ht27bx9NNPM2/ePI4ePcro0aP5yU9+wj333EPfvn39Lk9EOiAFuM+Ki4t55JFHmD9/PmlpaXzjG9/g/vvvZ8yYMaE8I0NEIqcA90lpaSmPPfYY8+bNIz09ne9973vMnDlT+7RFJGIK8ASrq6tjzpw5/OhHP6KmpobvfOc7/PCHP9RuEhFpNwV4Aq1bt467776b1atXM2nSJObMmcPgwYP9LktEAkpXdCSAc45nnnmGsWPHUlpayvz583n99dcV3iISE43A46yiooI777yTV155hRtvvJF58+Zx/vnn+12WiCQBBXgcbdu2jcmTJ7Nt2zaeeuopZs6c2WFuJiQiwacAj5OCggKmTp0KwNKlS7nmmmt8rkhEko32gcfBkiVLuO666+jTpw+rVq1SeItIXCjAPbZo0SImT57MsGHDePvtt3Xpu4jETZsBbmZzzWx/0+w7Ld97yMycmelyQRpH3tOmTWP06NEsW7ZMdwUUkbiKZAT+B2BSyxfNbABwLbDL45oCadWqVdx8883k5OSwZMkSevXq5XdJIpLk2gxw59wK4GArb80G/gMI/VyYW7du5aabbuL8889n8eLFZGRk+F2SiIRAVPvAzWwKUOqc+5fH9QROVVUVU6ZMARp3oeheJiKSKO0+jdDMugIPA9dFuHwekAeQlZUV9cSebU2864eGhgamT5/O1q1befnll+nZs6enE5d2xJ7jTT2Hg3r2RjTngQ8BBgH/arooJQtYbWZjnXMftVzYOZcP5APk5ua6WG6P2tFurfrEE0+wePFinnrqqVPnfHuto/WcCOo5HNRz7Nod4M659cCpa8HN7EMg1znn3dAzAFatWsUjjzzC17/+dWbOnOl3OSISQpGcRvgi8C4w3Mz2mNm34l9Wx3bkyBG++c1v0r9/f37zm9/o8ngR8UWbI3Dn3G1tvJ/tWTUB8dBDD7F9+3aWLVtGz549/S5HREJKV2K209tvv81vf/tbvvvd7+oSeRHxlQK8HU6cOEFeXh4DBw5k1qxZfpcjIiGnuxG2w89//nM2b97M66+/Trdu3fwuR0RCTiPwCO3evZuf/vSn3HLLLdxwww1+lyMiogCP1MMPP4xzjl/84hd+lyIiAijAI1JYWMjzzz/Pgw8+yMCBA/0uR0QEUIC3yTnHQw89RJ8+ffjBD37gdzkiIqfoIGYbli1bxvLly3nmmWc499xz/S5HROQUjcDbMGvWLC688EJmzJjhdykiIqfRCPwsli9fzooVK/jVr35Fenq63+WIiJxGI/CzmDVrFv369eOee+7xuxQRkU/RCPwMVq5cyfLly5k9e7ZG3yLSIWkEfgazZ88mIyND+75FpMNSgLdi9+7dLFiwgBkzZtC9e3e/yxERaZUCvBW//vWvcc5x3333+V2KiMgZKcBbOHr0KPn5+Xzta1/TVZci0qEpwFt46aWXqKio0DRpItLhRTKl2lwz229mRc1e+y8z22xm68zsr2aWNNPS/O53v+Piiy/myiuv9LsUEZGzimQE/gdgUovX3gQucc6NAoqBpLhJyIYNG3j33XeZMWOG5rkUkQ6vzQB3zq0ADrZ4bYlzrq7p6XtAVhxqS7jf//73dOrUiTvuuMPvUkRE2uTFhTz/G3jpTG+aWR6QB5CVlUV5eXlUG6msrIzqc5E6ceIEzz33HDfccANmFnWdXop3zx2Reg4H9eyNmALczB4G6oAXzrSMcy4fyAfIzc11mZmZUW8vls+2ZcGCBRw8eJB77703rttpr45US6Ko53BQz7GLOsDN7C7gy8AE55zzriR/vPjii1xwwQVMnDjR71JERCIS1WmEZjYJ+B4wxTlX7W1JiVdZWclrr73GLbfcQmpqqt/liIhEJJLTCF8E3gWGm9keM/sW8GugB/Cmma01s9/Euc64WrhwISdOnOC2227zuxQRkYi1uQvFOddaqv0+DrX45k9/+hPZ2dmMGzfO71JERCIW+isxP/74Y958801uvfVWnfstIoES+gB/+eWXqa+v1+4TEQmc0Af4woUL+cxnPsPIkSP9LkVEpF1CHeBVVVW89dZbTJ06VbtPRCRwQh3gb7zxBrW1tUydOtXvUkRE2i3UAb5o0SJ69+7N5z//eb9LERFpt9AGeG1tLa+99hpf/vKXSUvT3M4iEjyhDfCCggIOHTrElClT/C5FRCQqoQ3wV155hS5dunDdddf5XYqISFRCG+BvvPEG48eP16zzIhJYoQzwPXv2sHHjRo2+RSTQQhngS5cuBVCAi0ighTLAlyxZQt++fXX1pYgEWugCvKGhgTfffJPrrrtOV1+KSKCFLsDXrl1LeXk51157rd+liIjEJHQBvmTJEgBNnSYigRfJjDxzzWy/mRU1e+08M3vTzLY2/dkrvmV6Z+nSpYwcOZJ+/fr5XYqISEwiGYH/AZjU4rXvA393zg0F/t70vMOrqanhnXfe4Qtf+ILfpYiIxKzNAHfOrQAOtnh5KvBc0+PngK94XFdcrF69mmPHjjF+/Hi/SxERiVm0d3Hq65zbB+Cc22dm559pQTPLA/IAsrKyKC8vj2qDlZWVUX2uucWLFwOQk5MTdR2J5EXPQaOew0E9eyPut+FzzuUD+QC5ubkuMzMz6nXF8lmAwsJChg8fTk5OTkzrSaRYew4i9RwO6jl20Z6FUmZm/QCa/tzvXUnxUV9fT0FBgXafiEjSiDbAFwF3NT2+C/ibN+XEz/r166msrFSAi0jSiOQ0wheBd4HhZrbHzL4F/Ay41sy2Atc2Pe/QVqxYAaAAF5Gk0eY+cOfcbWd4a4LHtcTVihUrGDhwIBdddJHfpYiIeCIUV2I65/jnP//JVVdd5XcpIiKeCUWA7969m48++ogrrrjC71JERDwTigBfuXIlAOPGjfO5EhER74QmwLt06cKll17qdykiIp4JTYCPGTOGzp07+12KiIhnkj7Aa2trKSws1O4TEUk6SR/gRUVFHDt2TAEuIkkn6QNcBzBFJFklfYCvWrWKzMxMBg0a5HcpIiKeSvoAX7lyJePGjdMExiKSdJI6wI8cOcKmTZsYO3as36WIiHguqQN83bp1OOcYM2aM36WIiHguqQN8zZo1AApwEUlKSR/gmZmZ9O/f3+9SREQ8l/QBPnr0aB3AFJGklLQBXltbS1FRkXafiEjSiinAzexBM9tgZkVm9qKZpXtVWKw2bdpETU2NAlxEklbUAW5m/YGZQK5z7hIgFbjVq8JidfIA5ujRo32uREQkPmLdhZIGnGNmaUBXYG/sJXlj7dq1dO3alWHDhvldiohIXLQ5J+aZOOdKzexJYBdwDFjinFvScjkzywPyALKysigvL49qe5WVle1aftWqVeTk5HDo0KGottcRtLfnZKCew0E9eyPqADezXsBUYBBQAfzZzG53zv138+Wcc/lAPkBubq7LzMyMuthIP+ucY8OGDdx2220Rf6ajCnr90VDP4aCeYxfLLpSJwA7n3MfOuVrgL8D/8qas2OzcuZPKykrt/xaRpBZLgO8CrjCzrtZ4ovUEYJM3ZcWmqKgIgJEjR/pciYhI/EQd4M65lcACYDWwvmld+R7VFZMNGzYAkJOT43MlIiLxE/U+cADn3KPAox7V4pkNGzbQv39/evbs6XcpIiJxk5RXYm7YsIERI0b4XYaISFwlXYDX19ezceNGLrnkEr9LERGJq6QL8B07dnD8+HGNwEUk6SVdgJ88gKkAF5Fkl7QBrjNQRCTZJV2AFxUVMXDgQHr06OF3KSIicZV0Aa4zUEQkLJIqwOvq6ti8ebMCXERCIakCfPv27dTU1CjARSQUkirAN21qvBWLDmCKSBgkVYAXFxcDaBIHEQmFpAvwvn37kpGR4XcpIiJxl1QBvmXLFo2+RSQ0kirAi4uLGT58uN9liIgkRNIEeEVFBfv379cIXERCI2kCXAcwRSRsYgpwM+tpZgvMbLOZbTKzz3tVWHudDHDtQhGRsIhpRh7gaeB/nHPTzKwz0NWDmqKyZcsWUlJSGDx4sF8liIgkVNQBbmbnAuOB6QDOuRqgxpuy2q+4uJhBgwbRuXNnv0oQEUmoWEbgg4GPgXlmdilQCNzvnDvafCEzywPyALKysigvL49qY5WVlWd9f+PGjQwaNCjq9XdEbfWcjNRzOKhnb8QS4GnAZcB9zrmVZvY08H3gx80Xcs7l0zRbfW5ursvMzIx6g2f6bENDAyUlJUycOPGMywRVsvUTCfUcDuo5drEcxNwD7HHOrWx6voDGQE+4vXv3Ul1drQOYIhIqUQe4c+4jYLeZnUzNCcBGT6pqJ51CKCJhFOtZKPcBLzSdgVIC3B17Se23detWAIYOHerH5kVEfBFTgDvn1gK5HtUStZKSEjp37kz//v39LkVEJGGS4krMkpISBg0aREpKUrQjIhKRpEi8kpISXcAjIqGjABcRCajAB/ihQ4eoqKhQgItI6AQ+wEtKSgAU4CISOgpwEZGASpoAHzRokM+ViIgkVlIEeJ8+fejRo4ffpYiIJFTgA3z79u3afSIioRT4ANcphCISVoEO8NraWnbt2qUAF5FQCnSA7969m/r6egW4iIRSoANcpxCKSJgpwEVEAirQAb5jxw7S0tJ0G1kRCaVAB/iuXbsYMGAAqampfpciIpJwMQe4maWa2Roze9WLgtpj165dXHTRRYnerIhIh+DFCPx+YJMH62k3BbiIhFlMAW5mWcBNwO+8KSdydXV1lJaWMnDgwERvWkSkQ4h1BP4U8B9Agwe1tMvevXupr6/XCFxEQivqSY3N7MvAfudcoZl94SzL5QF5AFlZWZSXl0e1vcrKytOer1u3DoCMjIyo19nRtew5DNRzOKhnb8QyK/2VwBQzuxFIB841s/92zt3efCHnXD6QD5Cbm+syMzOj3mDzz1ZVVQEwcuRIYllnR5fMvZ2Jeg4H9Ry7qHehOOd+4JzLcs5lA7cCb7UM73jatWsXgHahiEhoBfY88J07d9K7d2+6devmdykiIr6IZRfKKc655cByL9YVKZ1CKCJhF9gRuAJcRMJOAS4iElCBDPDKykqqqqoU4CISaoEM8J07dwI6A0VEwi2QAX7yFEJdRi8iYRboANcIXETCLLAB3qlTJ/r27et3KSIivglsgGdlZZGSEsjyRUQ8EcgELC0t1TRqIhJ6gQzwvXv3KsBFJPQCF+DOOY3ARUQIYIBXVFRw7NgxBbiIhF7gAry0tBRAAS4ioRe4AN+7dy8AF154oc+ViIj4K3ABrhG4iEijwAa4RuAiEnaBDPDevXuTnp7udykiIr6KOsDNbICZLTOzTWa2wczu97KwM9m7d69G3yIixDalWh3w78651WbWAyg0szedcxs9qq1VOgdcRKRRLLPS73POrW56fBjYBMQ9WRXgIiKNPJnU2MyygTHAylbeywPyALKysigvL49qG5WVldTW1lJWVkavXr2iXk+QVFZW+l1CwqnncFDP3og5wM2sO/Ay8IBzrqrl+865fCAfIDc312VmZka9rePHj+OcY+jQocSyniAJS5/NqedwUM+xi+ksFDPrRGN4v+Cc+4s3JZ2ZTiEUEflELGehGPB7YJNz7pfelXRmuohHROQTsYzArwTuAL5kZmubfm70qK5WKcBFRD4R9T5w51wBYB7W0qbS0lI6deoUyn1nIiItBepKzL1799KvXz9NpSYiQsACXOeAi4h8IlABvm/fPp2BIiLSJHAB3q9fP7/LEBHpEAIT4MeOHaOiokIBLiLSJDABvn//fgAFuIhIE0/uhZIIZWVlgAJcJFnU1NSwfft2qqur/S4lYXbu3HnW97t27cqQIUPo3LlzROtTgIuIL7Zv307Pnj0ZPny4Tg0GGhoa+OijjygqKiInJyeiSWsC87emABdJLtXV1fTt21fh3SQlJYULLriAhoYG/va3v1FXV9f2ZxJQlyfKyspIS0vTVZgiSUThfbqUlBTMjLKysohuPxuYv72ysjL9thaRUDCz5BuBa/eJiMgnFOAiIgEVmLNQ9u/fz5VXXul3GSISBw888ABr1671dJ2jR4/mqaeeOusyv/zlL5k7dy4AM2bM4Ctf+QqTJk3iqquu4r333uPSSy/l7rvv5tFHH2X//v288MILjB07lqNHj3Lfffexfv166urqeOyxx5g6dSrV1dVMnz6dzZs3c/HFF/Phhx8yZ84ccnNz+fa3v83777/PsWPHmDZtGrNmzYq5x0AEeF1dHeXl5RqBi4hnCgsLmTdvHitXrsQ5x7hx47jmmmvYtm0bf/7zn8nPz+dzn/scf/zjHykoKGDRokU88cQTLFy4kMcff5wvfelLzJ07l4qKCsaOHcvEiRN59tln6dWrF+vWraOoqIjRo0ef2t7jjz/OeeedR319PRMmTGDdunWMGjUqph4CEeBlZWU45xTgIkmqrZFyPBQUFPDVr36Vbt26AXDzzTfz9ttvM2jQIEaOHAnAiBEjmDBhAmbGyJEj+fDDDwFYsmQJixYt4sknnwQa5+vdtWsXBQUF3H///QBccsklpwX0/Pnzyc/Pp66ujn379rFx40Z/A9zMJgFPA6nA75xzP4upmjPYt28foHPARcQ7zrlWX+/SpcupxykpKaeep6SknDozxDnHyy+/zPDhwyNa544dO3jyySd5//336dWrF9OnT+f48eMx9xDLnJipwBzgBiAHuM3McmKuqBU7duwAICsrKx6rF5EQGj9+PAsXLqS6upqjR4/y17/+lauvvjqiz15//fU888wzpwJ7zZo1AFx11VXMnz8fgI0bN7J+/XoAqqqq6NatGxkZGZSVlbF48WJPeohlBD4W2OacKwEwsz8BU4GNXhTW3Nq1a0lNTSUnJy6/H0QkhC677DKmT5/O2LFjgcaDmL169Yrosz/+8Y954IEHGDVqFM45srOzefXVV7n33nu56667GDVqFGPGjGHUqFFkZGQwdOhQxowZw4gRIxg8eLBnJ2TYmYb8bX7QbBowyTk3o+n5HcA459y/tVguD8gDyMrKuvzkb6r2eP7553nnnXd49tlno6o1qCorK8nIyPC7jIRSz+FQWVlJRUUFl19+ud+leKq+vp7a2lrS09PZvn07EyZMoLi4OOKbU0HjwdW33nqLyZMnn7ryvE+fPoXOudyWy8YyAm9tQuNP/TZwzuUD+QC5ubkumkvhH3zwQe64445QXkavnsMhjD1XVFT4XYLnqqur+eIXv0htbS3OOZ599tl2hfdJ6enp9OrVq81/F7EE+B5gQLPnWcDeGNYnIhJoPXr04IMPPkjY9mK5EvN9YKiZDTKzzsCtwCJvyhKRMGhoaPC7hA6lvX8fUQe4c64O+DfgDWATMN85tyHa9YlIuHTt2pV9+/YpxJucvB/4yd0vkYjpPHDn3OvA67GsQ0TCaciQIRQWFrJv3z7MWjukFj61tbWUlJTgnOOcc85pc/lAXIkpIsmnc+fODBw4kIULF0Z069SgO378eESz7DjnuPrqqzn33HPbXFYBLiK+ufDCC7nzzjupqqpK+l0phw4diug8827dutGzZ8+I1qkAFxFfde/ene7du/tdRtylp6d7frpoYO4HLiIip4v6SsyoNmb2MbAzyo9nAuUelhME6jkc1HM4xNLzQOdcn5YvJjTAY2FmH7R2KWkyU8/hoJ7DIR49axeKiEhAKcBFRAIqSAGe73cBPlDP4aCew8HzngOzD1xERE4XpBG4iIg0owAXEQmoDhfgZjbJzLaY2TYz+34r73cxs5ea3l9pZtmJr9JbEfQ83sxWm1ld00xIgRdBz981s41mts7M/m5mA/2o00sR9Px/zGy9ma01s4J4zTGbKG3122y5aWbmzCzwpxVG8B1PN7OPm77jtWY2I6YNOuc6zA+Ns9tvBwYDnYF/ATktlrkX+E3T41uBl/yuOwE9ZwOjgP8HTPO75gT1/EWga9Pjb4fkez632eMpwP/4XXc8+21argewAngPyPW77gR8x9OBX3u1zY42Aj81UbJzrgY4OVFyc1OB55oeLwAmWLDvRdlmz865D51z64BkudtPJD0vc85VNz19j8YZn4Iskp6rmj3tRitTFAZIJP8tA/wE+AVwPJHFxUmkPXumowV4f2B3s+d7ml5rdRnXOKlEJdA7IdXFRyQ9J5v29vwtYHFcK4q/iHo2s++Y2XYaQ21mgmqLhzb7NbMxwADn3KuJLCyOIv13/bWmXYMLzGxAK+9HrKMFeCQTJUc0mXKAJFs/kYi4ZzO7HcgF/iuuFcVfpJOAz3HODQG+B/wo7lXFz1n7NbMUYDbw7wmrKP4i+Y5fAbKdc6OApXyyNyEqHS3AI5ko+dQyZpYGZAAHE1JdfIRxcuiIejazicDDwBTn3IkE1RYv7f2e/wR8Ja4VxVdb/fYALgGWm9mHwBXAooAfyGzzO3bOHWj2b/n/ApfHssGOFuCRTJS8CLir6fE04C3XdHQgoMI4OXSbPTf97/VvaQzv/T7U6LVIeh7a7OlNwNYE1ue1s/brnKt0zmU657Kdc9k0HueY4pxL3JTu3ovkO+7X7OkUGucTjp7fR25bOZJ7I1BM49Hch5te+08av1yAdODPwDZgFTDY75oT0PPnaPztfhQ4AGzwu+YE9LwUKAPWNv0s8rvmBPT8NLChqd9lwAi/a45nvy2WXU7Az0KJ8Dv+adN3/K+m7/izsWxPl9KLiARUR9uFIiIiEVKAi4gElAJcRCSgFOAiIgGlABcRCSgFuISemfU0s3v9rkOkvRTgItCTxrtcigSKAlwEfgYMabo/c9DvuSIhogt5JPSaJgV51Tl3ic+liLSLRuAiIgGlABcRCSgFuAgcpvH2piKBogCX0HPOHQD+aWZFOogpQaKDmCIiAaURuIhIQCnARUQCSgEuIhJQCnARkYBSgIuIBJQCXEQkoBTgIiIB9f8BON3O4/HsuOkAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "#%matplotlib notebook \n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "from scipy.signal import StateSpace,lsim,step\n", "import control as ctl\n", "\n", "La = 1600e-6\n", "Ra = 4.8\n", "T = 0.0637\n", "J = 2.5e-5\n", "c = 1e-8\n", " \n", "I = np.array([[J, 0], [0, La]]) #[J 0;0 La];\n", "K = np.array([[c, -T], [T, Ra]])#[c -T;T Ra];\n", "invI = np.linalg.inv(I)\n", "\n", "A = - np.dot(invI,K)\n", "B = invI\n", "C = np.eye(2)\n", "D = np.zeros((2,2))\n", "\n", "sys1 = StateSpace(A,B,C,D)\n", "#print(sys1)\n", "#print(A[1][0])\n", "\n", "tfinal = 1001\n", "t = np.linspace(0, 0.5, tfinal)\n", "Tl = np.linspace(0, 0, tfinal) \n", "va = np.linspace(1, 1, tfinal) \n", "ut = Tl, va\n", "u = np.transpose(ut)\n", "\n", "tout, yout, xout = lsim(sys1, U=u, T=t)\n", "plt.figure(1)\n", "\n", "plt.plot(tout, yout[:,1], 'r', linewidth=1.5, label='ia')\n", "plt.legend(loc='best', shadow=True, framealpha=1)\n", "plt.grid(alpha=0.3)\n", "plt.xlabel('t')\n", "plt.show()\n", "\n", "plt.figure(2)\n", "plt.plot(tout, yout[:,0], 'k', linewidth=1.5, label='omega')\n", "plt.legend(loc='best', shadow=True, framealpha=1)\n", "plt.grid(alpha=0.3)\n", "plt.xlabel('t')\n", "plt.show()\n", "\n", "\n", "\n" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "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.7.6" } }, "nbformat": 4, "nbformat_minor": 4 }