{ "cells": [ { "cell_type": "code", "execution_count": 38, "metadata": {}, "outputs": [], "source": [ "import sympy\n", "import matplotlib.pyplot as plt\n", "import numpy\n", "import control \n", "from control.matlab import *\n", "sympy.init_printing()\n", "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 39, "metadata": {}, "outputs": [], "source": [ "from ipywidgets import interact" ] }, { "cell_type": "code", "execution_count": 90, "metadata": {}, "outputs": [], "source": [ "ts = numpy.linspace(0, 40,1000)\n", "\n", "s = tf([1,0],1)\n", "\n", "def secondorder(wn_in, zeta_in,Tp_in):\n", " G = wn_in**2/((s**2+2*zeta_in*wn_in*s + wn_in**2)*(1+Tp_in*s))\n", " \n", " #G = wn_in**2*(1+Tp_in*s)/((s**2+2*zeta_in*wn_in*s + wn_in**2))\n", " \n", " plt.figure(figsize=(12, 6))\n", " ax_impulse = plt.subplot2grid((2, 2), (0, 0))\n", " ax_step = plt.subplot2grid((2, 2), (1, 0))\n", " ax_complex = plt.subplot2grid((2, 2), (0, 1), rowspan=2)\n", "\n", " y,t = impulse(G,ts)\n", " ax_impulse.plot(t , y)\n", " ax_impulse.set_title('Impulse response')\n", " ax_impulse.set_ylim(-1, 2)\n", "\n", " y,t = step(G,ts)\n", " ax_step.set_title('Step response')\n", " ax_step.plot(t ,y)\n", " ax_step.set_ylim(0, 2)\n", "\n", "\n", " ax_complex.set_title('Poles plot')\n", " pole,zero = pzmap(G, plot = False)\n", " ax_complex.scatter(pole.real, pole.imag, marker='x', s=30)\n", " ax_complex.scatter(zero.real, zero.imag, marker='o', s=30)\n", " #ax_complex.axhline(0, color='black')\n", " #ax_complex.axvline(0, color='black')\n", " ax_complex.axis([-10,0.10, -2, 2])" ] }, { "cell_type": "code", "execution_count": 82, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtYAAAF1CAYAAADMcK0bAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nOzdd3yV9d3/8dcni0AYYYRN2EMcKERAcUDrLnW0taKtW6kt3ta2d++u+2d722VbW+u21K2to1YtrbgVN0pARDZhSQgjjLASRpLP74/rCh5jwsp1RpL38/E4j5xzXd/r+n7OyTlXPvme7zB3R0REREREGiYt2QGIiIiIiDQFSqxFRERERCKgxFpEREREJAJKrEVEREREIqDEWkREREQkAkqsRUREREQioMRaUpqZXWZmbyc7DhERiYaZ/cLMHk1wnSvM7JRE1inNkxJrOSC6KImISKzw70KFmW03s3Vm9oCZtU52XA1hZmPNrDjZcUjjpcRamj0zy0h2DCIijdSX3b01MBw4FvjfJMcjklRKrOWghd0z3jGzW8yszMyWmdnx4fZVZrbezC6NKf+gmd1jZi+b2TYze8PMeof7+piZxya3ZjbNzK6qo14L61xvZlvMbI6ZHRHua2FmN5vZJ2HLyT1m1vIA4t8E/CLcfoWZLTCzzWb2YkyM+6q33ucW7j/ezGaEx80ws+NrPc9fhrFsM7OXzKxTuC/bzB41s43hazzDzLqE+9qZ2X1mtsbMVpvZr8ws/dB/oyIiDePuq4HngZprY3czm2Jmm8ysyMyuru9YMxttZu+G17qPzGxszL7Lwr8x28xsuZl9o55z/MLMnjKzJ8Kys8xsWD1lW5jZn82sJLz9OdyWEz6H7mEr/HYz696Al0WaISXWcqhGAXOAjsDfgccJWisGAN8E7qj1leA3gF8CnYDZwN8Ooc7TgJOAQUAucAGwMdz3u3D70WEMPYAb9hP/MqAz8GszOxf4KfAVIA94C3jsAOqt97mZWQfgOeA2gtfpT8BzZtYx5tiLgMvDOLKA/w63Xwq0A3qFx14DVIT7HgIqw+d5TBjf5/4RERFJFDPrBZwFfBhuegwoBroDXwN+Y2ZfrOO4HgTXyV8BHQiugf80s7ww0b0NONPd2wDHE1xj63MO8I/wPH8HnjWzzDrK/QwYTfD3YhgwEvhfd98BnAmUuHvr8FZyEC+DiBJrOWTL3f0Bd68CniBIAG90913u/hKwmyDxq/Gcu7/p7rsILmrHhRfig7EHaAMMAczdF7j7GjMz4Grge+6+yd23Ab8BJuzjXCXufru7V7p7BfAt4LfhOSvD448OW5/rrPcAntuXgCXu/khYz2PAQuDLMcc+4O6LwxieJLjQ1zzXjsAAd69y95nuvjVstT4TuN7dd7j7euCW/TxXEZF4edbMyoC3gTcIEuhewAnAj9x9p7vPBu4FLq7j+G8CU919qrtXu/vLQCFBkg5QDRxhZi3dfY27z9tHLDPd/Sl330PQkJFNkEDX9g2Cv1fr3b0U+L96YhM5aEqs5VCti7lfAeDutbfFtlivqrnj7tuBTQQtGQfM3V8D7gDuBNaZ2WQza0vQwtwKmBl+lVgGvBBur8+qWo97A7fGHL8JMKDHPurd33PrDqysVc9Kgtb0Gmtj7pfz6Wv2CPAi8Hj4VeXvw5aX3kAmsCYm1r8QtHiLiCTaue6e6+693f07YSNBd6CmkaNG7Wtfjd7A+TXXs/CadgLQLWxBvoDgG7s1ZvacmQ3ZRyyx1+JqPm0xr632tXllPeVEDpoSa0mUva3TYReRDkAJsCPc3CqmbNf6TuLut7n7COBwgq4ZPwQ2ECTyh4cX+Fx3bxcOqKn3VLUerwK+FXN8rru3dPd391Hv/p5bCcEfjVj5wOp9xFXzPPe4+/+5+1CCrz/HA5eEce4COsXE2dbdD9/fOUVEEqQE6GBmbWK21XftWwU8Uuvam+PuNwG4+4vufirQjeAbv7/uo97Ya3Ea0DOMpa74Yq/N+THlav9tEDkoSqwlUc4ysxPMLIugP/L77r4q/BpuNfBNM0s3syuA/nWdwMyONbNRYcvtDmAnUBW2TPwVuMXMOodle5jZ6QcR3z3AT8zs8PD4dmZ2/r7q3d9zA6YCg8zsIjPLMLMLgKHAf/YXjJmNM7Mjw0GJWwm6hlSFXVBeAv5oZm3NLM3M+pvZyQfxXEVE4ia8/r0L/DYciH0UcCV1j615FPiymZ0e/g3ItmDKu55m1sXMzg77Wu8CtvPZa29tI8zsKxYMhr8+PGZ6HeUeA/437MfdiWA8Ts282uuAjmbW7hCeuogSa0mYvwM/J+gmMYKgj1uNqwlagDcStAi/W8852hIk0JsJvrrbCNwc7vsRUARMN7OtwCvA4AMNzt2fIRgA+Xh4/FyCvsz7q7fe5+buGwlamn8QHvM/wHh333AAIXUFniJIqhcQ9F2sufBfQjDQcX4Y01MErTkiIqniQqAPQUvwM8DPw/7TnxEm4ecQDB4vJWjB/iFBfpJGcP0sIbi+ngx8Zx91/oug68hmgj7TXwn7W9f2K4J+3HOAj4FZ4TbcfSFB4r0s7JqiLiJyUMxd33pIfJnZg0Cxuze5+U2b8nMTEWkszOwXBIO9v5nsWKR5U4u1iIiIiEgEGpxYm1kvM3vdgoU15pnZd+soY2Z2mwWTxM8xs+ENrVdERA6NrtsiIvHR4K4gZtaNYFqcWeEI4JkE0+/MjylzFvBfBPNSjgJudfdRDapYREQOia7bIiLx0eAW63DC9lnh/W0EA61qz1V5DvCwB6YDueGFXUREEkzXbRGR+Ii0j7WZ9SFYYvn9Wrt68NkFOYqpe6J4ERFJIF23RUSikxHVicKFMf5JsNTy1tq76zjkc31QzGwiMBEgJydnxJAh+1pgSUQkdc2cOXODu+9r9c+k03VbRCQQ1TU7ksQ6XDjjn8Df3P3pOooUE7MiEvWshuTuk4HJAAUFBV5YWBhFeCIiCWdmtZezTym6bouIfCqqa3YUs4IYcB+wwN3/VE+xKcAl4Sjz0cCWcAU5ERFJMF23RUTiI4oW6zEEKxx9bGazw20/BfIB3P0egqWdzyJYGa8cuDyCekVE5NDoui0iEgcNTqzd/W3q7osXW8aBSQ2tS0REGk7XbRGR+NDKiyIiIiIiEVBiLSIiIiISASXWIiIiIiIRUGItIiIiIhIBJdYiIiIiIhFQYi0iIiIiEgEl1iIiIiIiEVBiLSIiIiISASXWIiIiIiIRUGItIiIiIhIBJdYiIiIiIhFQYi0iIiIiEgEl1iIiIiIiEVBiLSIiIiISASXWIiIiIiIRUGItIiIiIhKBSBJrM7vfzNab2dx69o81sy1mNju83RBFvSIiIiIiqSIjovM8CNwBPLyPMm+5+/iI6hMRERERSSmRtFi7+5vApijOJSIiIiLSGCWyj/VxZvaRmT1vZofXVcDMJppZoZkVlpaWJjA0EREREZGGSVRiPQvo7e7DgNuBZ+sq5O6T3b3A3Qvy8vISFJqIiIiISMMlJLF2963uvj28PxXINLNOiahbRERERCQREpJYm1lXM7Pw/siw3o2JqFtEREREJBEimRXEzB4DxgKdzKwY+DmQCeDu9wBfA75tZpVABTDB3T2KukVEREREUkEkibW7X7if/XcQTMcnIiIiItIkaeVFEREREZEIKLEWEREREYmAEmsRERERkQgosRYRERERiYASaxERERGRCCixFhERERGJgBJrEREREZEIKLEWEWlmzOx+M1tvZnPr2T/WzLaY2ezwdkOiYxQRaYwiWSBGREQalQcJFu16eB9l3nL38YkJR0SkaVCLtYhIM+PubwKbkh2HiEhTo8RaRETqcpyZfWRmz5vZ4fUVMrOJZlZoZoWlpaWJjE9EJOUosRYRkdpmAb3dfRhwO/BsfQXdfbK7F7h7QV5eXsICFBFJRUqsRUTkM9x9q7tvD+9PBTLNrFOSwxIRSXlKrEVE5DPMrKuZWXh/JMHfio3JjUpEJPVpVhARkWbGzB4DxgKdzKwY+DmQCeDu9wBfA75tZpVABTDB3T1J4YqINBpKrEVEmhl3v3A/++8gmI5PREQOQiRdQQ5gsQEzs9vMrMjM5pjZ8CjqFRERERFJFVH1sX4QOGMf+88EBoa3icDdEdUrIiIiIpISIkmsD2CxgXOAhz0wHcg1s25R1C0iIiIikgoSNStID2BVzOPicNtnaKEBEREREWmsEpVYWx3bPjfCXAsNiIiIiEhjlajEuhjoFfO4J1CSoLpFREREROIuUYn1FOCScHaQ0cAWd1+ToLpFREREROIuknmsD2CxganAWUARUA5cHkW9IiIiIiKpIpLE+gAWG3BgUhR1iYiIiIikokR1BRERERERadKUWIuIiIiIRECJtYiIiIhIBJRYi4iIiIhEQIm1iIiIiEgElFiLiIiIiEQgkun2mqpVm8r52/ufMH3ZRrZU7KFjThaj+3XkgmN70atDq2SHJyIiIiIpRC3Wdaiqdv78ymLG3TyNe99aRlZGGkf0aIcDd00r4gt/nMbvXljI7srqZIcqIiIiIilCLda17NxTxfWPz+aFeWs59+ju/PjMw+jaLnvv/pKyCv708mLunraUwhWbuPubI+jUukUSIxYRERGRVKAW6xjV1c51j33IC/PWcsP4ofx5wjGfSaoBuue25Obzh3Hbhcfw8eotXDh5Ohu270pSxCIiIiKSKpRYx7j5pUW8NH8dN4wfyhUn9N1n2bOHdeeBy0ayanM5F9/3ATt2VSYoShERERFJRUqsQ+8UbeCuaUu5cGQ+l4/pc0DHHNe/I/d8cwSL1m7l+0/Oprra4xukiIiIiKQsJdbAtp17+J+n5tAvL4eff3koZnbAx44d3JmfnnUYL85bx91vLI1jlCIiIiKSypRYA7e/VkTJlgpuPn8Y2ZnpB338lSf05UtHdeOWlxczd/WWOEQoIiIiIqmu2SfWKzfu4IF3lnP+iJ4Mz29/SOcwM3597hF0bJ3F9U/MZueeqoijFBEREZFUF0libWZnmNkiMysysx/Xsf8yMys1s9nh7aoo6o3CTc8vJDM9jR+cNrhB58ltlcXN5w+jaP12bnllcUTRiYiIiEhj0eDE2szSgTuBM4GhwIVmNrSOok+4+9Hh7d6G1huFuau38PzctUw8qR9d2mbv/4D9OHFgHueP6Ml9by1nybptEUQoIiIiIo1FFC3WI4Eid1/m7ruBx4FzIjhv3N35ehFtsjP2O7XewfjxmUPIaZHB//vXXNw1S4iIiIhIcxFFYt0DWBXzuDjcVttXzWyOmT1lZr0iqLdBitZv44V5a7n0uD60zc6M7LwdW7fgf84YzPRlm5jyUUlk5xURERGR1BZFYl3X3HS1m2r/DfRx96OAV4CH6jyR2UQzKzSzwtLS0ghCq989bywjOyM90tbqGhcem89RPdtx0/MLNZBRREREpJmIIrEuBmJboHsCn2mqdfeN7l6z7vdfgRF1ncjdJ7t7gbsX5OXlRRBa3bZU7OE/c0o4b3gPOuRkRX7+tDTjp2cdxpotO7n/neWRn19EREREUk8UifUMYKCZ9TWzLGACMCW2gJl1i3l4NrAggnoP2ZSPSti5p5oLj82PWx2j+3XklMM6c/frS9m4fdf+DxARERGRRq3BibW7VwLXAi8SJMxPuvs8M7vRzM4Oi11nZvPM7CPgOuCyhtbbEI9/8AlDu7XliB5t41rPj88cwo7dldz+WlFc6xERERGR5ItkHmt3n+rug9y9v7v/Otx2g7tPCe//xN0Pd/dh7j7O3RdGUe+hmLt6C/NKtjJhZK+DWrr8UAzo3IYLjs3n0ekr+WRjeVzrEhE5UGZ2v5mtN7O59ew3M7stXJtgjpkNT3SMIiKNUbNbefGxDz6hRUYa5xxd18Ql0bv+lIGkpxm3vbYkIfWJiByAB4Ez9rH/TGBgeJsI3J2AmEREGr1mlViX765kyuwSvnRkN9q1jG6KvX3p0jabb47uzdOzillWuj0hdYqI7Iu7vwls2keRc4CHPTAdyK01VkYkblaXVXDHa0v2rgXh7tzx2hJWl1UkOTKR/WtWifVzc9awbVclFxyb2Gm0rzm5Py0y0rntVbVai0ijcKDrEyR0mlRpHp6ZVczNLy3mxv/Mx9258T/zufmlxTwzqzjZoYnsV0ayA0ikJ2asol+nHEb27ZDQevPatOCS43sz+c1lTBo3gIFd2iS0fhGRg3Qg6xMEG90nA5MBCgoKtNysNNikcQPYuGM3D7yzggfeWQHA5WP6MGncgOQGJnIAmk2LddH6bRSu3MwFx8Z/0GJdvnVSf1plpvNntVqLSOrb7/oEIvFiZtwwfuhntt0wfmhS/naLHKxmk1g//sEqMtKMr47omZT6O+RkcfmYvjw3Zw0L1mxNSgwiIgdoCnBJODvIaGCLu69JdlDSPNR0/4hV0y1EJNU1i8R6V2UVT3+4mlOHdqFT6xZJi+PqE/vRpkUGf35lcdJiEBExs8eA94DBZlZsZlea2TVmdk1YZCqwDCgiWC33O0kKVZqhO18v4oF3VnD5mD4s/+1ZXD6mDw+8s4I7X9eaEJL6mkUf65fnr2PTjt0JH7RYW7tWmVx5Yl/+/MoSPi7ewpE92yU1HhFpntz9wv3sd2BSgsIR+YzzhgffLE8aN2Bvt5COOVl7t4uksmaRWD8xYxU9clty4sC8ZIfCFSf05YF3VnDLK4u5/7Jjkx1OUm3buYfF67azdstOyip2A5CdkU6Xttn0aN+S3h1akZamPnUiIs1Jj9yWXPuFgXsfm9lnHouksiafWK/aVM5bSzbsXagl2dpmZzLxpH784cVFzPpkM8Pz2yc7pIRxd2Z9UsYLc9cwbVEpRaXb2VeXuTYtMjiqVztOHpTHqUO70rdTTuKCFRERETlITT6xfrJwFWZwfkFyu4HEuuz4Ptz39nJueXkxj1w5KtnhxN2uyiqmzC7hvreXs3DtNjLTjdH9OvLlYd05vHtbuue2pENOFgaU765i7dadfLKxnI+Ky5i5cjO/mbqQ30xdyJCubbhwZD7nHtMjYQv8iIiIiByoJp1YV1ZV84/CYk4elEeP3JbJDmevnBYZXHNyP34zdSEzVmzi2D6JnVc7UdydF+et5TdTF/LJpnKGdG3D7756JGce2Y222fUnxn065TC6X0e+HvaJL95czivz1/H0h6v5+ZR5/Pb5BXxleE+uOak/+R1bJerpiIiIiOxTk06sX124nrVbd/KLsw9Pdiifc/HoPkx+czl/fGkRj088LtnhRG51WQU/emoObxdtYFCX1jxw+bGMHZR3SPOQ9mzfisvG9OWyMX2Zu3oLj05fyVOFxTz+wSecPay7Ft0RERGRlNCkp9t7dPpKurbN5pTDOic7lM9pmZXOd8b2Z/qyTbxbtCHZ4UTqnzOLOeOWN/nwk83ceM7hTL3uRMYN7hzJ5P5H9GjHTV89ird+NI4rT+jLS/PXcdqf3+QHT37E6rKKCKIXEREROTRNNrFevmEHby3ZwEWj8slIT82nedGofLq2zeZPLy9uEhPf76qs4idPf8wP/vERh3VvywvXn8Qlx/WJy+vfpW02P/vSUN7+0Re4+sR+/HtOCeNunsavn5vP5h27I69PREREZH9SM+OMwL1vLSMz3ZiQ5Lmr9yU7M51JXxhA4crNvLmkcbdar9+6kwsnT+exDz7hO2P789jVo+nVIf79nzvkZPHTsw7j9f8ey9nDunPv28s56Q+vc9e0Iip2V8W9/sZgT1U123buYeP2XWwp38OOXZXsqqyiurrx/zMnIiKSSiLpY21mZwC3AunAve5+U639LYCHgRHARuACd18RRd11KSmr4MnCVZxf0IvObbPjVU0kvl7Qk3umLeVPLy3ipIGdIukukWjLSrdzyf0fsGnHbu76xnDOOrJbwmPokduSm88fxlUn9uUPLyzi9y8s4qF3V/C9UwbxtRE9U/Zbi4bavGM3RaXbWbp+O8s27GDNlp2s37qT0m27KN2+i4rdVVTWk0CbQbuWmbRvlUVuq0w6tMqic9tserZvSY/clnTPbUmP9i3p2jY7JaaqFBERSXUNTqzNLB24EzgVKAZmmNkUd58fU+xKYLO7DzCzCcDvgAsaWnd9fv/CQgC+M7Z/vKqITIuMdP7rCwP48dMf89rC9XzxsC7JDumgfLSqjMsfnIEBj08czVE9c5Maz5CubbnvsmP5YPkmbnp+AT9++mP++tYyfnj6EE4/vEuj/MelxvZdlcxZVcaHq8r48JMyPiouo3Tbrr37szLS6NYum85tWnBYt7ac1KYFOS3Syc5IJzsznayMNKqqnT1V1VRWOzv3VFFWvofN5bspK9/Dmi07mb2qjI21utJkpaeR37EVfTvl0C8vh36dcujbqTX98nLomJPVqF9TERGRKEXRYj0SKHL3ZQBm9jhwDhCbWJ8D/CK8/xRwh5mZR9yxePOO3TwyfSXPzi7h+lMG0rN945iK7asjenLXtKXc/NJixg7u3GhaB98t2sBVDxfSISeLR64clVILuIzs24F/fvt4Xpq/jt+/sJBrHp3JMfm5/PiMIYzq1zHZ4R2QyqpqPiou443FG3hzcSlzisuoaXzul5fDiQM7MbRbW/p3bs2AvNb0yG0ZyUqVFburWF1WQUlZBavLKli5sZzlG7azrHQHbywqZXdV9d6ybbIzwkQ7h/yOOfTu0Ir8jq3o3aEVeW1apGzSXVlVzZaKPWwu38OWit1s3rGHsoo9lIX/ZGwu37338eAubbnhy0OTHbKIiDQCUSTWPYBVMY+Lgdqrnuwt4+6VZrYF6Ah8pmOxmU0EJgLk5+cfdCAPvruCW19dwmlDuzBp3ICDPj5ZMtPT+MFpg/ju47N5elZxSi1mU5/pyzZyxUMz6N0hh0euHJmSXW7MjNMP78oXh3TmqZnF3PLKYi6YPJ0TB3bi2yf357j+HVMu8Sspq+DNxaW8uaSUt5dsYOvOStIMhvXK5dpxAxjRpwNH98ylXav4LZDTMiudAZ1bM6Bz68/tq6p2SsoqWFq6neUbdrB8ww6Wle5gxorNTPmohNheJ9mZaeR3aEV+hxy652aT17oFeW2CW+c22XRsnUXr7AxysjIO6Z9Jd6diTxVbKyrZUrGHLRV72Br+rLmVle9mc/nnk+ZtOyvrPW+aQW7YPSa3ZSaZGan1HhERkdQVRWJd11+d2i3RB1IGd58MTAYoKCg46NbsCSN7cerQLhzevW3KJUz7c/aw7jzwzgpufmkRXzqqG62yUneK8ZkrN3HFgzPo2b4Vf7t6FJ1at0h2SPuUkZ7GhJH5nHN0Dx56bwX3vrWci+59n2E923HNyf05dWiXpPXB3rmniveXbwqS6cWlLFm/HYCubbM544iunDQojxMGdCK3VVZS4qstPc3o1aEVvTq0Yuzgz+7bXVlN8eZyPtkU3jaWszL8WbhyE2Xle+o9b8vMdHJaZNC6RToZ6WmkGaSZYWakWZDQ79xTxa7K6uC2p4qdldVU7WcAZruWmUGC3CqL9q2y6NcpZ2/SXNO3PLdVFrlhX/N2rTJp0yIjkpZ/ERFpfqLI3oqB2CbWnkBJPWWKzSwDaAdsiqDuz+jWriXd2qXOCosHw8z4f+MP46t3v8fkN5dx/SmDkh1SnT5aVcZl98+gS9ts/n5V6ifVsVpmpXPNyf257Pg+/HNWMZPfXMa3/zaLzm1acH5BTy4oyI/7So7uzsK123hrSSlvLdnAB8s3sauymqyMNEb17cDXC3px8uA8BnZu3ej+OczKSKNfXmv65X2+pRuC6Rg3bN8dDKzctosN23exY1cl23dVhj+rKN9dSWW14+5UV0O1O9XupKcZLTLSyc5Mo0VGOi0y0miRmUab7EzaZmfSrmVwa9syI/iZnUnblpmNpluViIg0DVEk1jOAgWbWF1gNTAAuqlVmCnAp8B7wNeC1qPtXNwUjenfgS0d24y9vLOPCkfl0SbHuFXNXb+Hi+94nNyeTv189KiW7fxyI7Mx0vjGqNxOOzeeVBet4YsYq7p62lDtfX8pRPdtx+uFdOeWwLgzq0vDktrKqmkXrtjHrkzJmrtjE20Ub2bA9GHA4qEtrvjGqNycN6sSovh1pmZUexdNLWS0y0umRG8w4IiIi0hQ1OLEO+0xfC7xIMN3e/e4+z8xuBArdfQpwH/CImRURtFRPaGi9TdWPzhjCy/PX8bvnF/KnC45Odjh7LVy7lYvve5822Zn8/arRjfabgVjpaUEf7NMP70pJWQXPzl7Ni/PW8YcXF/GHFxfRvlUmI3p34Kie7eiXFwzQ69I2m3YtM8mM6TpSXe2UVQTzRK/dupOl67dTVLqdJeu28/HqLZSH82l3at2C4/t35MSBnThxYB5d2zXOf0xERESkbpaqDccFBQVeWFiY7DCS4o8vLeL214r4+1WjOH5Ap2SHQ9H67UyY/B7pacaT3zqO3h1TZ/aPeFi7ZSdvLimlcMUmZqzYzPINOz5XpmVm0Lpc5U5lVTW1u/q2yc5gQOfWHNWjHcN7t2d4fnt6tm/Z6Lp3yKEzs5nuXpDsOBKpOV+3RaRxi+qanboj5JqxSeMG8O+PSvjZs3N5/rsnkp2ZvC4Cyzfs4KK/TgeMv189uskn1QBd22Xz9YJefD2cnaVidxUrNgYzYGwIVy/cujMYiJeWZmSlp9ExJ4uO4awX/fJyyGudulPNiYiISHwosU5B2Znp/OrcI/nmfe9z1+tFfP+0wfs/KA5WbSrnor9Op7LaeXziaPrXMyitqWuZlc5h3dpyWLe2yQ5FREREUljTXOe5CThhYCfOPbo7d7+xlHklWxJe/6pN5UyYPJ3y3VU8euUoBnVpk/AYRERERBoTJdYp7IYvH05uqyy++/hsKsIBcImwuqyCi+6dzrade/jbVaMY2l0ttSIiIiL7o8Q6hXXIyeJPXx9G0frt/GbqgoTUuWZLBRf9dTpl5Xt49KpRHNGjXULqFREREWnslFinuBMH5nHVCX15ZPpKnpuzJq511XT/2Lh9Nw9fMZKjeubGtT4RERGRpkSJdSPwwzMGMzw/lx/8YzYfF8env/Witdv46t3vUla+h4evHMkx+e3jUo+IiIhIU6XEuhFokZHOXy4uoGNOC5oD8T4AACAASURBVK56eAYlZRWRnn/mys18/S/vAfDkt45juJJqERERkYOmxLqRyGvTgr9eUkD5riou/Ot01m7ZGcl5nyxcxYWTp9O+VSb//PbxDO6q2T9EREREDoUS60ZkaPe2PHTlSDZu380Fk99jaen2Qz7Xzj1V/L9n5/I/T83h2L7teeY7Y+jVoVWE0YqIiIg0L0qsG5nh+e155MqRbN9ZyXl3vsOrC9Yd9Dk+/GQzX7rtLR6ZvpKrTujLQ5ePpH1OVhyiFREREWk+lFg3Qsfkt+fZSWPontuSKx8q5PtPzD6gftcrNuzgusc+5Ly73t278Mv/jh9KRrreBiIiIiINpSXNG6leHVrxr2vHcOdrRdw1bSn/nlPCl47sxhlHdOWIHu3o1LoFVdVO8eYKZqzYxIvz1vLWkg1kZ6Zx7bgBfOvkfrTJzkz20xARERFpMpRYN2ItMtL5/mmD+fqxvbj3reU88+Fqnp1dUmfZnu1b8r1TBnHhyF50bpud4EhFREREmj4l1k1Az/at+MXZh/OzLx3GnOIyFq3dzuby3aSZ0a1dNsN65dKnYyvMLNmhioiIiDRZSqybkMz0NEb07sCI3h2SHYqIpDAzOwO4FUgH7nX3m2rtvwz4A7A63HSHu9+b0CBFRBqhBiXWZtYBeALoA6wAvu7um+soVwV8HD78xN3Pbki9IiJyaMwsHbgTOBUoBmaY2RR3n1+r6BPufm3CAxQRacQaOh3Ej4FX3X0g8Gr4uC4V7n50eFNSLSKSPCOBIndf5u67gceBc5Ick4hIk9DQxPoc4KHw/kPAuQ08n4iIxFcPYFXM4+JwW21fNbM5ZvaUmfWq72RmNtHMCs2ssLS0NOpYRUQalYYm1l3cfQ1A+LNzPeWywwvvdDOrN/nWBVpEJO7qGsXstR7/G+jj7kcBr/BpA8rnD3Sf7O4F7l6Ql5cXYZgiIo3PfvtYm9krQNc6dv3sIOrJd/cSM+sHvGZmH7v70tqF3H0yMBmgoKCg9oVeREQarhiIbYHuCXxmnk533xjz8K/A7xIQl4hIo7ffxNrdT6lvn5mtM7Nu7r7GzLoB6+s5R0n4c5mZTQOOAT6XWIuISNzNAAaaWV+CWT8mABfFFqi5rocPzwYWJDZEEZHGqaFdQaYAl4b3LwX+VbuAmbU3sxbh/U7AGKD26HMREUkAd68ErgVeJEiYn3T3eWZ2o5nVDC6/zszmmdlHwHXAZcmJVkSkcWnoPNY3AU+a2ZXAJ8D5AGZWAFzj7lcBhwF/MbNqgkT+pjqmdRIRkQRx96nA1Frbboi5/xPgJ4mOS0SksWtQYh32w/tiHdsLgavC++8CRzakHhERERGRVNfQriAiIiIiIoISaxERERGRSCixFhERERGJgBJrEREREZEIKLEWEREREYmAEmsRERERkQgosRYRERERiYASaxERERGRCCixFhERERGJgBJrEREREZEIKLEWEREREYmAEmsRERERkQgosRYRERERiYASaxERERGRCCixFhERERGJQIMSazM738zmmVm1mRXso9wZZrbIzIrM7McNqVNEREREJBU1tMV6LvAV4M36CphZOnAncCYwFLjQzIY2sF4RERERkZSS0ZCD3X0BgJntq9hIoMjdl4VlHwfOAeY3pG4RERERkVSSiD7WPYBVMY+Lw20iIiIiIk3GfluszewVoGsdu37m7v86gDrqas72euqaCEwMH243s0UHcP7aOgEbDuG4eEiVWFIlDlAsdUmVOCB1YkmVOODQY+kddSAiIpLa9ptYu/spDayjGOgV87gnUFJPXZOByQ2pzMwK3b3egZSJlCqxpEocoFhSOQ5InVhSJQ5IrVhERCS1JaIryAxgoJn1NbMsYAIwJQH1ioiIiIgkTEOn2zvPzIqB44DnzOzFcHt3M5sK4O6VwLXAi8AC4El3n9ewsEVEREREUktDZwV5Bnimju0lwFkxj6cCUxtS10FoUFeSiKVKLKkSByiWuqRKHJA6saRKHJBasYiISAoz9zrHEYqIiByUgoICLywsTHYYIiIHzcxmRjGeRkuai4iIiIhEoNEm1vtbJt3MWpjZE+H+982sT5zi6GVmr5vZgnB59+/WUWasmW0xs9nh7YY4xbLCzD4O6/hcs5EFbgtfkzlmNjxOcQyOea6zzWyrmV1fq0zcXhMzu9/M1pvZ3JhtHczsZTNbEv5sX8+xl4ZllpjZpXGI4w9mtjB8/Z8xs9x6jt3n7zKiWH5hZqtjfgdn1XPsPj9rEcTxREwMK8xsdj3HRv2a1PnZTcZ7RUREmoZGmVjbgS2TfiWw2d0HALcAv4tTOJXAD9z9MGA0MKmOWADecvejw9uNcYoFYFxYR11fZ5wJDAxvE4G74xGAuy+qea7ACKCcOvriE7/X5EHgjFrbfgy86u4DgVfDx59hZh2AnwOjCFYM/Xl9SVUD4ngZOMLdjwIWAz/Zx/H7+l1GEQvALTG/g8+NgzjAz1qD4nD3C2LeL/8Ent7H8VG+JvV9dpPxXhERkSagUSbWxCyT7u67gZpl0mOdAzwU3n8K+KLZvtdePxTuvsbdZ4X3txHMfJKqK0ueAzzsgelArpl1i3OdXwSWuvvKONezl7u/CWyqtTn2/fAQcG4dh54OvOzum9x9M0ESXFcyeshxuPtL4Uw5ANMJ5nWPu3pekwNxIJ+1SOIIP59fBx471PMfZCz1fXYT/l4REZGmobEm1geyTPreMmEiswXoGM+gwu4mxwDv17H7ODP7yMyeN7PD4xSCAy+Z2UwLVrGsLRnLy0+g/kQpEa9JjS7uvgaChAroXEeZRL8+VwDP17Nvf7/LqFwbdku5v54W10S+JicC69x9ST374/aa1PrspuJ7RUREGoHGmlgfyDLpB7yUehTMrDXB19jXu/vWWrtnAb3dfRhwO/BsnMIY4+7DCb62n2RmJ9UOs45j4vmaZAFnA/+oY3eiXpODkbDXx8x+RtAV4W/1FNnf7zIKdwP9gaOBNcAf6wq1jm3xes9cyL5bq+Pymuzns1vvYXVsa1RTLO2v77wlaJyKiEhT0lgT6wNZJn1vGTPLANpxaF+F75eZZRL8Yf6bu3+uf6i7b3X37eH9qUCmmXWKOo5w/nDcfT1Bn+aRtYoc8PLyETkTmOXu62rvSNRrEmNdTbeX8Of6Osok5PUJB7qNB77h9cx3eQC/ywZz93XuXuXu1cBf66kjUa9JBvAV4In6ysTjNanns5sy75V4SbFxKs3K6rIK7nhtCTUffXfnjteWsLqsIsmRiUgUGmtifSDLpE8Bakbqfw14rb4kpiHCfqH3AQvc/U/1lOla07/bzEYSvO4bI44jx8za1NwHTgPm1io2BbjEAqOBLTVfecdJvS2QiXhNaol9P1wK/KuOMi8Cp5lZ+7BbxGnhtsiY2RnAj4Cz3b28njIH8ruMIpbY/vXn1VPHgXzWonAKsNDdi+vaGY/XZB+f3ZR4r8RZyoxTaW6emVXMzS8t5sb/zMfdufE/87n5pcU8M6vOt76INDINWnkxWdy90sxqlklPB+5393lmdiNQ6O5TCP5gPmJmRQQt1RPiFM4Y4GLg45hpwn4K5Iex3kOQ2H/bzCqBCmBCHJL8LsAz4d+9DODv7v6CmV0TE8dUghUxiwhm6rg84hj2MrNWwKnAt2K2xcYSt9fEzB4DxgKdzKyYYPaGm4AnzexK4BPg/LBsAXCNu1/l7pvM7JcEySTAje5+yN9y1BPHT4AWwMvh72q6u19jZt2Be939LOr5XR5qHPuIZayZHU3QhWEF4e8qNpb6PmtRxuHu91FHX/x4vybU/9lN+HslCerqIz6qvjLh+6BmnMqGhETYRE0aN4CNO3bzwDsreOCdFQBcPqYPk8YNSG5gIhIJrbwokTOzE4DfA4cDVQSzLVzv7jPM7DLgKnc/IYkhijRrZnY+cLq7XxU+vhgY6e7/FVNmXlimOHy8NCyzsda5JhJM30l+fv6IlSsTNgFQo+Xu9P3Jp7NbLv/tWejLAJHkMq28KKnIzNoC/yEYkNiBoNXr/4BdSYjFzEzvcZHPi2ycirtPdvcCdy/Iy8uLU7hNR033j1g13UJEpPFT0iFRGwTg7o+FA+Mqwrmb55jZYcA9BNPsbTezMtg7+8DNZvaJma0zs3vMrGW4b6yZFZvZT81sgwWr732jvsrNbJqZ/drM3iHo7tLPzNqZ2X1mtsaClQZ/FQ7ewswGmNkbFqwCucHMnog5l5vZdWa2LNz3h5pE3czSzOx/zWylBSsJPmxm7cJ9fcJjLw2f0wYLZgGpOe9IMyu0YEXKdWb2p5h9o83sXTMrs2AqwrGR/WZEPpUy41SamztfL+KBd1Zw+Zg+LP/tWVw+pg8PvLOCO18vSnZoIhKBRtnHWlLaYqDKzB4iGBA1PVxAA3dfEPazrt0V5HdAP4Jp3/YAfwdu4NNVCbsCnQhav0cDU82s0N0X1RPDxQSzHSwimBbtH8A6YACQQ9Civgr4C/BL4CVgHJAF1P4a6LxwW2vglfCc9wKXhbdxBLNGPAzcEdZd4wRgMME/Gx+Y2dPuvgC4FbjV3R+xYKq3IwDMrAfwXHiOFwgW1/mnmQ1x99J6nqvIQUuxcSrNynnDgzWhJo0bgJlxw/ihdMzJ2rtdRBo39bGWyIUt0z8imOmhK8GgyavdfV3tPtbhLAPbgaPcfWm47TiCwWl9wxbbV4B27r4j3P8k8LG7/7KOuqcBb7r7DeHjLgQD0HLdvSLcdiEw0d3HmdnDwE6CwWfFtc7lwJk1g+TM7DvAV939i2b2KvBPd78r3DeYYJaKlgRfqy8HesX0T/0A+JO7P25mbwKvA7e7+4aY+n5EsNz5xTHbXgxfi5rZGURSVkFBgRcWFiY7DBGRg6Y+1pKy3H2Bu1/m7j0JWmO7A3+up3ge0AqYGXZ/KCNorY3trLm5JqkOrQzPWZ/Y2Q56A5nAmpjz/4VPV9P7H4JW7Q/MbJ6ZXbGPc8XW2z18HLsvg2D2ihprY+6XE7R6QzA/8CBgoZnNMLPxMbGeXxNnGOsJQLyXnRcREZEIqCuIxJW7LzSzB/l02r3aX5FsIJhu73B3X13PadqbWU5Mcp3Pvucwjq1jFcHAyU4eLG1fO761wNWwdzaTV8zsTXev6fDYC6iZXi6fTwd4lRAkwsTsqyTocrLP73Q9WLL7wrC/9leAp8ysYxjrI+5+9b6OFxERkdSkFmuJlJkNMbMfmFnP8HEvgoVipodF1gE9wwFTxKz6d4uZdQ6P6WFmp9c69f+ZWZaZnUiwamFdy6R/TrgAzkvAH82sbTjosL+ZnRzWdX5NrMBmgqS8KuYUP7RgEZBewHf5dGXAx4DvhYO/WgO/AZ6oK3mv4zX6ppnlhc+9LNxcBTwKfNnMTjezdDPLtmDwpjpfioiINAJKrCVq2wgWmnjfzHYQJNRzgR+E+18jaAFea2Y1/Yt/RLBozXQz20rQp3pwzDnXEiS9JcDfCBbpWHgQMV1CMDBxfniep/i0e8WxYazbCWZB+K67L4859l/ATGA2wcDC+8Lt9wOPAG8S9KfeCfwXB+YMYF5Y560Ei+PsdPdVBKvd/RQoJWjB/iH6nIqIiDQKGrwoKS0cvPho2F870XU7MDCmW4iI7IMGL4pIY6XBiyIiIiIiKaTBibWZ9TKz181sQTirwnfrKGNmdpuZFZnZHDMb3tB6RURERERSSRSzglQCP3D3WWbWhmDatJfdPXbN1jOBgeFtFHB3+FNkn9x9GvuZZSOOdVsy6hUREZHGqcEt1u6+xt1nhfe3AQsIVsiLdQ7wsAemA7lmprl5RURERKTJiLSPtZn1AY4B3q+1qwefXWijmM8n3yIiIiIijVZkC8SEc/n+E7je3bfW3l3HIZ+bjsTMJgITAXJyckYMGTIkqvBERBJq5syZG9w9b/8lRUSkqYgksTazTIKk+m/u/nQdRYoJVrCr0ZNPV7Dby90nA5NB0zaJSONmZiv3X0pERJqSKGYFMYJFMxa4+5/qKTYFuCScHWQ0sCVcEU9EREREpEmIosV6DHAx8LGZzQ63/RTIB3D3e4CpwFkEq+uVA5dHUK+IiIiISMpocGLt7m9Tdx/q2DIOTGpoXSIiIiIiqUorL4qIiIiIRECJtYiIiIhIBJRYi4iIiIhEQIm1iIiIiEgElFiLiIiIiERAibWIiIiISASUWIuIiIiIRECJtYiIiIhIBJRYi4iIiIhEQIm1iIiIiEgElFiLiIiIiERAibWIiIiISASUWIuIiIiIRECJtYiIiIhIBJRYi4iIiIhEQIm1iIiIiEgElFiLiIiIiERAibWISDNiZh3M7GUzWxL+bF9PuSozmx3epiQ6ThGRxkiJtYhI8/Jj4FV3Hwi8Gj6uS4W7Hx3ezk5ceCIijZcSaxGR5uUc4KHw/kPAuUmMRUSkSVFiLSLSvHRx9zUA4c/O9ZTLNrNCM5tuZvUm32Y2MSxXWFpaGo94RUQajYxkByAiItEys1eArnXs+tlBnCbf3UvMrB/wmpl97O5Laxdy98nAZICCggI/pIBFRJoIJdYiIk2Mu59S3z4zW2dm3dx9jZl1A9bXc46S8OcyM5sGHAN8LrEWEZFPqSuIiEjzMgW4NLx/KfCv2gXMrL2ZtQjvdwLGAPMTFqGISCOlxFpEpHm5CTjVzJYAp4aPMbMCM7s3LHMYUGhmHwGvAze5uxJrEZH9UFcQEZFmxN03Al+sY3shcFV4/13gyASHJiLS6KnFWkREREQkAkqsRUREREQiEElibWb3m9l6M5tbz/6xZrYlZnncG6KoV0REREQkVUTVx/pB4A7g4X2Uecvdx0dUX5Oyu7KahWu3snH7bjrkZDG0e1sy0/VlgoiIiEhjEkli7e5vmlmfKM7VnOzcU8U9byzl4fdWsmnH7r3b27XM5Juj8/n22AG0bqHxpSIiIiKNQSKztuPCqZtKgP9293kJrDvlFG8u58oHC1m0bhunHNaFc4/pTvfclpSUVfDcnDXcNW0pz35Ywt3fHM5RPXOTHa6IiIiI7EeiEutZQG93325mZwHPAgNrFzKzicBEgPz8/ASFlnirNpXztXvepWJ3FQ9cfizjBnfeu294fnvGH9WdmSs3cd1js7ngL9O56xvDGTek8z7OKCIiIiLJlpCOvO6+1d23h/enApnhal61y0129wJ3L8jLy0tEaAm3pWIPVzw4g4rdVTx5zXGfSapjjejdgWcmHU//zjl869GZvLd0Y4IjFREREZGDkZDE2sy6mpmF90eG9Ta7TNHd+enTH7N8ww7uuXgEQ7q23Wf5zm2yeeSKUfTu0IqrHy5kaen2BEUqIiIiIgcrqun2HgPeAwabWbGZXWlm15jZNWGRrwFzwz7WtwET3N2jqLsxmfJRCc99vIbvnTqI4/t/rsG+Tu1zsnjoipFkZaTxnUdnUb67Ms5RioiIiMihiGpWkAv3s/8Ogun4mq0tFXv4xZR5HN0rl2+d1O+gju2e25JbJxzNJfd/wC+mzOP3XxsWpyhFRERE5FBpsuQEuf3VJZRV7OHX5x1BxiHMUX3iwDy+M7Y/TxYWM23R+jhEKCIiIiINocQ6AT7ZWM5D763g6yN6cXj3dod8nuu+OJABnVvzk6c/ZtvOPdEFKCIiIiINpsQ6Ae5+owgz4/unDWrQeVpkpPOHrx3Fuq07+f0LiyKKTkRERESioMQ6zkrKKnhqZjEXFPSiS9vsBp/vmPz2XDy6N397fyUL126NIEIRERERiYIS6zib/OYy3OFbJx/cgMV9+d6pg2jbMpMb/z2fZji5ioiIiEhKUmIdR1sq9vDEjFWcc3QPerZvFdl5c1tl8f1TB/Hu0o28NH9dZOcVERERkUOnxDqOnp5VTMWeKi4f0yfyc180Mp9BXVrzm6kL2FNVHfn5RUREROTgKLGOE3fnkekrObpXLkf0OPSZQOqTkZ7Gj84YwsqN5fyjsDjy84uIiIjIwVFiHSfvLd3IstIdXDy6d9zq+MKQzhyTn8vtry1h556quNUjIiIiIvunxDpOHpuxitxWmXzpqG5xq8PM+OFpg1mzZSd/e/+TuNUjIiIiIvunxDoOtu3cw0vz1vLlo7qTnZke17qOH9CJ4/t35K7Xi9ixqzKudYmIiIhI/ZRYx8Hzc9eyq7Ka84b3SEh9/336YDbu2M2D765ISH0iIiIi8nlKrOPgmVmr6dsph2N65SakvuH57Rk3OI9731qmVmsRERGRJFFiHbGSsgqmL9/IuUf3wMwSVu+1XxjI5vI9/F19rUVERESSQol1xKZ8VII7nHtM94TWO6J3e8YM6Mjkt5ZphhARERGRJFBiHbHn565lWM929O6Yk/C6rx03kNJtu/hH4aqE1y0iIiLS3CmxjlBJWQUfrSrj9CO6JqX+0f06UNC7Pfe8sYzdlVqNUURERCSRlFhH6KV5awE44/DkJNZmxrVfGMDqsgqe/XB1UmIQERERaa6UWEfo+blrGdylDf3yWicthpMH5XFkj3bcNa2Iyiq1WouIiIgkihLriGzYvosZKzYlrRtIjZpW6xUby3nu4zVJjUVERESkOclIdgBNxSvz11HtyesGEuvUw7owuEsb7nitiC8f1Z20tMRN+5dKysp38+GqMhau2cbqsnK2VFRSXe20zEqnc5sW9OmUw5CubRjarS0Z6fofU0QkFawuq+CZWcVMGjcAM8PdufP1Is4b3pMeuS2THZ7IPimxjsiL89aS36EVh3Vrk+xQSEszJn1hANc99iEvzFvLWUd2S3ZICbNt5x6e+XA1//6ohJkrN1PtwfbcVpnktswkPc2o2F3F+m27qAx3tmmRwah+HRl/VDdOHdqFnBb6WEjTZWbnA78ADgNGunthPeXOAG4F0oF73f2mhAUpzdozs4q5+aXFbNyxmxvGD+XG/8zngXdWAMGaDSKpTBlEBCp2V/Hu0o1cNCo/oYvC7MuXjuzGn19ZzG2vLuGMw7s2+VbrsvLd3D1tKY9MX0n57iqGdG3DpHEDGDOgE4d1bUu7VpmfKV9ZVU3x5grmrN7Ce0s3Mm3Rel5ZsI7szDTGH9WdK8b0ZWj3tkl6NiJxNRf4CvCX+gqYWTpwJ3AqUAzMMLMp7j4/MSFKczZp3AA27tjNA++s2JtQXz6mD5PGDUhuYCIHQIl1BKYv28iuymrGDe6c7FD2Sk8zrh03gO8/+RGvLFjHaSnQRSUeqqqdB95Zzq2vLmH7rkrOGdady8f0Zdh+lpPPSE+jT6cc+nTK4exh3amudgpXbuaZD1fz7IereWpmMcf378i14wZw/IBOCXo2IvHn7guA/TUCjASK3H1ZWPZx4BxAibXEnZlxw/ihe5NqgBvGD02ZhiuRfVHH0gi8vmg9LTPTGdm3Q7JD+Yyzh3Wnd8dW3P5aEe6e7HAit2TdNr5y97v86rkFFPRuzwvfPYk/Tzhmv0l1XdLSjJF9O/DbrxzJez/5Aj86YwhLS7dz0b3vc9FfpzNz5aY4PAORlNUDiF1pqjjc9jlmNtHMCs2ssLS0NCHBSdPm7tz4n8/+D3fjf+Y3yb9j0vQosW4gd+e1hesZM6Aj2ZnpyQ7nMzLS05g0dgAfr97CtEVN6w/eP2cWM/72t/lk4w5uu/AY7r/sWAZ3jaZ/e26rLL49tj9v/HAcN4wfyuJ12/jq3e9x1UMzWFq6PZI6ROLJzF4xs7l13M450FPUsa3OrMbdJ7t7gbsX5OXlHXrQIqE7Xy/igXdWcPmYPiz/7VlcPqYPD7yzgjtfL0p2aCL7pa4gDbS0dAfFmyu45uT+yQ6lTucN78Gtry7h1leXMHZwXqP/Km13ZTW//M98Hpm+ktH9OnD7hcPJa9MiLnVlZ6ZzxQl9mTCyFw+8s4K7py3l9Fve5NLj+3DdFwZ+rt+2SKpw91MaeIpioFfM455ASQPPKXJAzhveE2DvrCA3jB9Kx5ysvdtFUplarBto2qL1AIwdnJotNZnpaXxnXH9mryrj7aINyQ6nQbbt3MPlD37AI9NX8q2T+vHolaPillTHapWVwaRxA3j9v8dyfkFP7n9nOWNvfp1H3luhRXhi6GvaJmUGMNDM+ppZFjABmJLkmKSZ6JHbkmu/MHBvQ1CwPsNATbUnjUIkLdZmdj8wHljv7kfUsd8Ipm06CygHLnP3WVHUnWyvL1rPoC6t6dm+VbJDqdf/b+/O46uq7r2Pf36ZB5IQEghhDmESUBRTJmnF4TqgVWvVl0O9Kg7XqrW9t316HV7a+9Q+z320L2/tvdpaq1atvQ61tXIrFaFFUBAEFFDmJICGMQlDGDJnPX+cAz0eTsIJ2WdKvu/XK6+cnL3OXj/W2cOPvdde66ozB/Hk38r5+fzNTB9RmJBXrfccbOCW3yxn466DPH71BL55ZvSvXPTNSeffrzyNb00ZyiN/XsdDb63lpQ+38cAlpzBjVOLfDWiPc47t++vZtPsgldWH2bG/gV119ew80MD+I80camzhUEML9c2tACQZJJmRnGTkZKSSm5lCXmYqeZmpFOVkUNw7gwF5mfTPy2BgfiaD87NIS9H/8aPFzL4B/BfQF3jbzFY55y40swH4htWb6ZxrMbN7gLn4htt73jm3NoZhi4gkBK+6grwAPAm81M7yi4GR/p/JwC/9vxPaocYWPtqyl1lnlcQ6lA6lpyTz7RmlPPzWWj6srGVaaWKNcrF9fz3XPbOU6oONPHtTGTNiPPrKuAF5vHL7FN5dt5v/O2c9t/xmOdNHFPLAzFO6xRB9uw40sHLbPlZs28uqL/azefchDjW2HFuemZpMce8MivMyGJSfRa/0ZLLTUsjyj//tnKPNOVpaHXUNLdQ1NFNX30ztoSbW7qij+mDjl+pLTjIG52dSUpjN8L69fL/9r4ty07vtf1hixTn3JvBmiPd34Lv4cfTvOcCcKIYmIpLwPEmsnXOLzGxYB0UuB15yvnvFS82st5kVO+cSes7tpRW1NLc6zh4Vn91AAl1TNpgntHWpSAAAGOdJREFU/1bOE/M3M3V4QcIkK9UHG/nWs8vYd6SJ/759MmcMyY91SIDv1uSF4/pzzuh+vLx0G//5t81c8l/vc9XEQfzgwtEU5WbEOsSwHWps4YPNNSzYsIcPymvYvr8egIzUJE4b2JtvThzIqP45jC7KYUS/XuRlpnZp+2lqaWN3XQM79tdTta+eLTWH2VJzmMqaw3xYWUtD89+712SlJVNSmP2lZLukMJuSvtnkZqiPu4iIxJdoPbzY3tBNCZ1YL6moJT0liYlD4yPZ60hGajJ3nzOCH81ey3ubquNqzO32HDjSzI3PLWPXgQZevm1S3CTVgdJSkpg1vYRvThzEkws28+KSbfx5zU5mTR/GbdOHk5+dFusQj+Oco7LmMAs27GHBxj18tGUvza2OnPQUzhpRyKzpJZQNzWfsgFxSIzDVe1pKEoP7ZDG4T9Zxt63a2hy76hqorD7MlppDVPqT7jVVB5jz6c5jM2kCFPZK9yfb2Qzuk0VRbgb9czPon5dOUW4GOUq8RUQkyqKVWIc1dJOZ3QHcATBkyJBIx9RlSypqKBuWH3fD7LXnuklDeO6DLTz2zkbOHtk3rmdjPNzYws0vfERl9WGeu7mMM4fG1xjhwfKyUnnwkrHcOGUYj83dwFMLKvjN4q3cOGUot361hH45sb2C3dDcyrIte48l09tqjwAwsl8vZp1Vwjlj+nHm0PyIJNKdkZRkDOidyYDemUwf+eUuS40trXxee+RYsl1ZfYgtNYeZt243tYebjltXdloyvbPSyMlIITczldyMVHIyUkhLTiIl2UhNTiIlyUhONlpbHc2tbTT5fze3ttHQ3HpsFs8HLxkbrSYQEZEEFq3EOqyhm5xzzwDPAJSVlcX1EAO1hxrZsOsg/+vC0bEOJWxpKUl8/4JRfPfVVcxevYMrzgg530PMNTS3csdvV7Cm6gBPXT+Rr46M/642Rw0pyOLJ6ydy73kH+cWCcn79fiUvLNnK5acP4IbJQzltUF7UuuF8sfcICzdV897GahaX11Df3Ep6ShLTSgu4bXoJM0b3Y3Cf+H3oNlh6SjIji3IYWXT8eOX1Ta3srmtg54EGdtc1sKvO97uu/u99vLfvr+dgQzPNrW20+BPoljZHS5sjNclITUkiNTnp2Ov0lCSy01PQYCciIhKuaCXWs4F7/NPiTgYOJHz/6krfTHxTSwtiHEnnfP20AfxqYSWPz9vIzFOL4240hpbWNr7zyicsLq/l8asncNH4xJyKfVRRDk9cewbfO38Uv1pUyVurtvP6iirGD8zlitMHctH4/p6PJHO4sYUV2/axcGM1CzftoaL6MACD8jO56sxBnDumH1NL428iIy9kpiUfm6JeREQkVrwabu8VYAZQaGZVwI+AVADn3NP4niyfCZTjG27vFi/qjaUlFTX0Sk/htIF5sQ6lU5KSjH+9eAw3Pf8RLy/dxqzp8TOiSVub44dvrGHeut3878vGxWRIPa8NK8zm3688lQdmjuFPq3bw6kef85O31/OTt9czfmAu00oL+cqwPkwc0puCXuGPyd3Y0kpl9WE27T7Iqi/2s2LrPtbtrKO1zZGWksSU4QVcP3koM0b3ZXhhdsI8rCoiIpLIvBoV5LoTLHfA3V7UFS8+rKhlckkfUmLcJ/VkfG1kIV8dWcgT8zdx2ekDKOxEQhcpzjl+NHstf/xkOz+4YBQ3TRsW65A8lZORyo1ThnLjlKFsrTnMO2t38df1u3lh8VaeWVQJQH5WKsP79qI4L4O8zFRyM1NJNqO5rY2mljb2Hm5iT10ju+sa2Lb3CK3+J/kyUpM4Y3A+d80opWxYHyYN60NmWve7Ki0iIhLvNKX5Sdh5oJ7KmsNcPzn+H7AMxcz40dfHcdETi3jsnQ08dtWEWIfET+duPDaj4t3njIh1OBE1rDCbO88u5c6zS2lobmX1F/v5dPsBKqp9D+St3VHHgXpfv2CHb5zn1CSjoFc6fXPSGVOcw8xTixnVP4dRRb0o7dsr5g8dioiIiBLrk/JhRS1Awk20EmhEv17cOr2EXy2q5LpJQ2I6lN0v36vgF+9VcP3kIdx38Zge1W0hIzWZycMLmDw8sfrqi4iIyPF0meskLKmoJT8rlTH9jx+dIJF857yRFOWm8/Bba2lpbTvxByLg5aXbePSdDVw2YQCPXD6+RyXVIiIi0r0ose4k5xwfVtQytbQgrseBDkev9BQeunQsn24/wDPvV0a9/j99sp2H3vqM88b04/FrJpCc4O0pIiIiPZsS6076fO8Rtu+vZ2oCdwMJdMmpxVxyajE/m7eJDbvqolbvnE938i+vr2JKSQFP3TBRfYRFREQk4Smb6aQlx/pXd48+sWbGI1eMJy8zjX9+bTVNLZHvEjJ/3W7ufeUTJg7J57mby7rluMoiIiLS8yix7qQlFbUU5aYzvBtNRNEnO43/d+WprN9Zx0/eXhfRuhZuquau333MuAG5/OaWr5CVpudnRUREpHtQYt0Jvv7VNUwrLex2D9mdP7aIO742nJc+3MYbK6siUsc7n+3i9hdXUNqvFy/OmkRORmpE6hERERGJBSXWnbBp9yFqDjUl3DTm4frhhaOZVlrAA29+yvKtez1d9xsrq7jrdysZNzCXV26fTO+sNE/XLyIiIhJrSqw7YUlFDdB9+lcHS0lO4snrJzKodyazXljO2h0HurzOtjbHz+dv5ge/X8200kJevlVJtYiIiHRPSqw7YUlFLUP6ZDEoPyvWoURMn+w0fnvbZHLSU7jh2WWs6MKV64MNzdz58kp+Nn8TV04cyHM3l5Gdrj7VIiIi0j0psQ5Ta5tjaWVtt71aHWhg70xevWMq+Vlp3PDsMv6wsgrnXKfWsbi8houeeJ+/btjDw5eO5fGrJ5CeotE/REREpPtSYh2mtTsOcLChpdv2rw42pCCLP3x7GhMG9+b7v1/N3f/9Mdv315/wc5/XHuG7r37CDc8uIz0lidf/aQqzppd0u4c9RURERILpvnyYjo5f3VMSa/B1C3nl9ik8vbCCn8/fzPx1e7h0QjFfnzCAMwb3Ji8zFedgx4F6Vmzdx/+s3sGCjXtIS0ni7nNK+c65IzVGtYiIiPQYSqzDtKSilpH9etEvJyPWoURVcpJx9zkjuOKMgfxiQTmzV+3gjx9vByAtOYlW52ht83UT6ZuTzj+dXcrN04ZRlNuz2klEREREiXUYmlraWL5lL9eUDYp1KDEzsHcm/+cbp/LQpWNZsXUf63YeoPZwE6lJSRTlZXDG4N6cUpxLcpK6fIiIiEjPpMQ6DKur9lPf3MrU0sJYhxJzGanJTB9ZyPSRagsRERGRQHp4MQxLymsxgynD+8Q6FBERERGJU0qsw7CkooZxA3I1sYmIiIiItEuJ9QnUN7Xyyef7maZuICIiIiLSASXWJ7By2z6aWtt61DB7IiIiItJ5SqxPYElFDSlJxqRh6l8tIiIiIu1TYn0CSypqOX1wb7LTNYCKiIiIiLRPiXUH6hqaWVO1n2nqBiIiIiIiJ6DEugPLt+ylzaHxq0VERETkhJRYd2BxeS3pKUmcMaR3rEMREekyM7vazNaaWZuZlXVQbquZfWpmq8xsRTRjFBFJZOo43IElFTWUDcsnIzU51qGIiHjhM+BK4FdhlD3HOVcT4XhERLoVXbFuR/XBRjbsOshZI9QNRES6B+fceufcxljHISLSXXmSWJvZRWa20czKzey+EMtvNrNq/23FVWZ2mxf1RtKSCt+FmrPUv1pEeh4HvGtmK83sjo4KmtkdZrbCzFZUV1dHKTwRkfjU5a4gZpYMPAX8A1AFLDez2c65dUFFX3PO3dPV+qJlSXktuRkpjB+YF+tQRETCZmbzgf4hFj3onHsrzNWc5ZzbYWb9gHlmtsE5tyhUQefcM8AzAGVlZe6kghYR6Sa86GM9CSh3zlUCmNmrwOVAcGKdUBZX1DBleAHJSRbrUEREwuacO9+Ddezw/95jZm/iO86HTKxFROTvvOgKMhD4IuDvKv97wb5pZmvM7A0zG+xBvRHzee0RqvbVq3+1iPQ4ZpZtZjlHXwMX4HvoUURETsCLxDrUJd3g24H/Awxzzp0GzAdeDLmiOOmrt/ho/2ol1iLSjZjZN8ysCpgKvG1mc/3vDzCzOf5iRcAHZrYa+Ah42zn3TmwiFhFJLF50BakCAq9ADwJ2BBZwztUG/Plr4NFQK4qXvnqLy2soyk2ntG92rEIQEfGcc+5N4M0Q7+8AZvpfVwITohyaiEi34MUV6+XASDMrMbM04FpgdmABMysO+PMyYL0H9UZEW5tjSUUtZ5UWYqb+1SIiIiISni5fsXbOtZjZPcBcIBl43jm31sx+DKxwzs0G7jWzy4AWYC9wc1frjZQNuw6y93AT09QNREREREQ6wZOZF51zc4A5Qe89HPD6fuB+L+qKtGPjV48oiHEkIiIiIpJINPNikEWbayjtm01xXmasQxERERGRBKLEOkB9UytLK2uZMbpfrEMRERERkQSjxDrA0spamlraOHtU31iHIiIiIiIJRol1gIWbqslITWJSSZ9YhyIiIiIiCUaJdYD3Nu5h6vACMlKTYx2KiIiIiCQYJdZ+W2sOs7X2iPpXi4iIiMhJUWLtt3CTbwp19a8WERERkZOhxNrvvY17GFaQxbBCTWMuIiIiIp2nxBo41NjC4opazhmjbiAiIiIicnKUWAOLNlXT1NLGheP6xzoUEREREUlQSqyBd9fuIj8rlbKh+bEORUREREQSVI9PrJtb2/jrhj2cd0oRKck9vjlERERE5CT1+ExyaWUtBxta1A1ERERERLqkxyfW767dTWZqMl8dWRjrUEREREQkgfXoxLqltY2/fLaTGaP7arZFEREREemSHp1YL66opeZQE5efPiDWoYiIiIhIguvRifVbn2wnJyNF05iLiIiISJf12MS6vqmVuWt3MXN8sbqBiIiIiEiX9djEet763RxualU3EBERERHxRI9NrF9Z9jmD8jOZPLwg1qGIiIiISDfQIxPriupDfFhZy3WThpCcZLEOR0RERES6gR6ZWL+y7HNSkoxrygbHOhQRERER6SZ6XGJd19DMa8u/4KLx/embkx7rcERERESkm+hxifXLS7dxsLGFO88ujXUoIiIiItKN9KjEur6plec/2MLXRvVl/MC8WIcjIiIiIt1Ij0qsf/1+JTWHmvjOuSNiHYqIiIiIdDM9JrHedaCBpxdWcOG4Ir4yrE+swxERERGRbqZHJNbOOX74hzW0OccDM0+JdTgiIiIi0g15klib2UVmttHMys3svhDL083sNf/yZWY2zIt6w/X0wkoWbarmgZmnMLQgO5pVi4jEFTP7qZltMLM1ZvammfVup1yHx3URETlelxNrM0sGngIuBsYC15nZ2KBitwL7nHMjgJ8Bj3a13nD9duk2Hn1nA5ecVsy3Jg+NVrUiIvFqHjDeOXcasAm4P7hAmMd1EREJ4sUV60lAuXOu0jnXBLwKXB5U5nLgRf/rN4DzzCxiUx42tbSxcts+vv3ySh7602ecO6Yf/3HNBJI0y6KI9HDOuXedcy3+P5cCg0IUC+e4LiIiQVI8WMdA4IuAv6uAye2Vcc61mNkBoACo8aD+Y559v5Kfzt1IY0sbANlpyXz/H0Zx1zkjNHW5iMjxZgGvhXg/nOO6iIgE8SKxDpWxupMog5ndAdzh//OQmW08iXgKCUjY730E7j2JlXjkS7HEULzEAYollHiJA+InlniJA04+lpj1PTOz+UD/EIsedM695S/zINAC/C7UKkK8d9wx27+ewON2o5l91vmIIy6etqdAiqtzFFfnKK7OGe3FSrxIrKuAwQF/DwJ2tFOmysxSgDxgb/CKnHPPAM90JRgzW+GcK+vKOrwSL7HESxygWOI5DoifWOIlDoivWMLlnDu/o+VmdhNwKXCecy5UwhzOcf1oXceO2/HaVoqrcxRX5yiuzonnuLxYjxd9rJcDI82sxMzSgGuB2UFlZgM3+V9fBfytnYO5iIhEkJldBPwrcJlz7kg7xcI5rouISJAuJ9b+h2DuAeYC64HXnXNrzezHZnaZv9hzQIGZlQP/AmjoJhGR2HgSyAHmmdkqM3sawMwGmNkcaP+4HquARUQShRddQXDOzQHmBL33cMDrBuBqL+oKQ5e6kngsXmKJlzhAsYQSL3FA/MQSL3FAfMXSZf5hT0O9vwOYGfD3ccf1MMRrWymuzlFcnaO4Oqdbx2XqkSEiIiIi0nU9YkpzEREREZFIS9jEOl6mUTezwWa2wMzWm9laM/tuiDIzzOyAvz/jKjN7ONS6PIhlq5l96q/juKdbzec//W2yxswmRiiO0QH/1lVmVmdm3wsqE7E2MbPnzWxP4LBfZtbHzOaZ2Wb/7/x2PnuTv8xm/8gJXscR7nTSHX6XHsXyb2a2PeA7mNnOZz2b2rqdOF4LiGGrma1q57Net0nIfTcW20oiMbOr/e3VZmZlQcvu928nG83swnY+X+I/Jm/2f/dpEYgxJttUGHFFfZ8LM66YHJc6iCcuzu9BdcbNuT5EvXFx7g+qM6Z5QFA90c0JnHMJ9wMkAxXAcCANWA2MDSpzF/C0//W1wGsRiqUYmOh/nYNviuDgWGYAf45Cu2wFCjtYPhP4C74xaqcAy6L0Xe0ChkarTYCvAROBzwLeewy4z//6PuDREJ/rA1T6f+f7X+d7HMcFQIr/9aOh4gjnu/Qoln8DfhDG99fhvtbVOIKWPw48HKU2CbnvxmJbSaQf4BR8472+B5QFvD/Wv32kAyX+7SY5xOdfB671v34a+HaE443aNhVGLFHf58KMKybHpZP99xOl83tQnXFzru/s90IMzv0hvtOo5gFB9UQ1J0jUK9ZxM426c26nc+5j/+uD+J6gH+h1PR65HHjJ+SwFeptZcYTrPA+ocM5ti3A9xzjnFnH8OOmB28OLwBUhPnohMM85t9c5tw+YB1zkZRwuvOmkPddOm4TD06mtO4rDv39eA7xysuvvZCzt7btR31YSiXNuvXMu1ORdlwOvOucanXNbgHJ8288x/u/4XHzHZGi/fT0R7W3KI1GfTj5Wx6V2xM35PVCCneuDxeLcHyjqeUCgaOcEiZpYh5puN3gD/9I06sDRadQjxn876gxgWYjFU81stZn9xczGRSgEB7xrZivNNxtasHDazWvX0v5JLRptclSRc24n+A6QQL8QZaLdPrPwXUUI5UTfpVfu8d8afL6dW2HRbJOvArudc5vbWR6xNgnad+NxW0kE4bRJAbA/IImLdLvFbJvqQDztc6HE+rgUl+f3QHFwrg8Wj+f+QPGSBwSK2HHek+H2YsCzadS9Yma9gD8A33PO1QUt/hjfLZBD/j51fwJGRiCMs5xzO8ysH74xajf4/6d2LMwQn4lkm6QBlwH3h1gcrTbpjKi1j3U8nTSc+Lv0wi+BR/D9Gx/Bd8t8VnCoIT4XqW3mOjq+shiRNgned8O88BXVfSnaLIwp0UN9LMR7ETsuhxlj1LepjuIihvtcOO0VJ8eluDu/f6ni+DjXB4urc3+gBMwDAp1UuyVqYu3ZNOpeMLNUfDva75xzfwxeHrjzOefmmNkvzKzQOVfjZRzONw4tzrk9ZvYmvltqgTtX2NMUe+Ri4GPn3O4QsUalTQLsNrNi59xO/y2wPSHKVOHr83XUIHz9SD1lJ55OOpzvsssCvxcz+zXw5xDForLN+PfRK4Ez2ysTiTZpZ9+Nm20lVtwJpkRvRzjbSg2+29Ap/iuNJ709nSjGWG1T4bZdtPe5MNorLo5LxNn5PVC8nOtD1Btv5/5A8ZQHBIrYcT5Ru4LEzTTq/n5dzwHrnXP/0U6Z/kf7f5nZJHztXutxHNlmlnP0Nb6HUT4LKjYb+EfzmQIcOHorJELavVoUjTYJErg93ASEuuo2F7jAzPL9t2gv8L/nGQtjOukwv0svYgnsY/eNduqI1tTW5wMbnHNVoRZGok062HfjYltJQLOBa803YkMJvitPHwUW8B+DF+A7JkP77euFqG9TJxJn+1xgXHFzXCKOzu+B4uVcH6LOeDz3B4qnPCBQ5I7zLspPr3r1g+8p1034nh5+0P/ej/EdGAAygN/je4DmI2B4hOKYju/WwBpglf9nJnAncKe/zD3AWnxPNy8FpkUgjuH+9a/213W0TQLjMOApf5t9SsAT/RGIJwvfDpIX8F5U2gTfTrwTaMb3P85b8fW/+yuw2f+7j79sGfBswGdn+beZcuCWCMRRjq/P1tFt5eiT7QOAOR19lxGI5bf+7WANvoNMcXAs7e1rXsbhf/+Fo9tGQNlIt0l7+27Ut5VE+sGXEFYBjcBuYG7Asgf928lG4OKA9+cAAwK+y4/87fZ7ID1CcUZ9mwojpqjvc2HGFZPjUgfxxMX5PSimuDjXh4grrs79QbHFLA8IiiOqOYFmXhQRERER8UCidgUREREREYkrSqxFRERERDygxFpERERExANKrEVEREREPKDEWkRERETEA0qsRUREREQ8oMRaRERERMQDSqxFRERERDzw/wF8k1G8wAgyXgAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "secondorder(1.0, 0.2, 0.2)" ] }, { "cell_type": "code", "execution_count": 92, "metadata": {}, "outputs": [ { "data": { "application/vnd.jupyter.widget-view+json": { "model_id": "e9a06aa45f6a4aa6884a7803c298ccd4", "version_major": 2, "version_minor": 0 }, "text/plain": [ "interactive(children=(IntSlider(value=5, description='wn_in', max=10, min=1), FloatSlider(value=0.55, descript…" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "interact(secondorder, wn_in=(1, 10), zeta_in=(0.1, 1.0), Tp_in=(0.01, 20));" ] } ], "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.3" } }, "nbformat": 4, "nbformat_minor": 4 }