{ "cells": [ { "cell_type": "code", "execution_count": 50, "metadata": { "inputHidden": false, "outputHidden": false }, "outputs": [], "source": [ "import numpy as np\n", "from numpy import exp\n", "from cmath import sqrt\n", "import sympy as sp\n", "import matplotlib.pyplot as plt\n", "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 51, "metadata": { "inputHidden": false, "outputHidden": false }, "outputs": [], "source": [ "t = np.linspace(0,1,100)" ] }, { "cell_type": "code", "execution_count": 52, "metadata": { "inputHidden": false, "outputHidden": false }, "outputs": [], "source": [ "omega0 = 10\n", "Q = 0.3\n", "gamma = omega0/Q\n", "omegav = omega0 * sqrt(1 - 1/(4*Q**2))\n", "A = 1.0\n", "phi = 0" ] }, { "cell_type": "code", "execution_count": 53, "metadata": { "inputHidden": false, "outputHidden": false }, "outputs": [], "source": [ "def z(gamma, A, omegav, phi):\n", " return A * np.exp(-gamma/2*t) * np.exp(1j*(omegav*t + phi))" ] }, { "cell_type": "code", "execution_count": 54, "metadata": { "inputHidden": false, "outputHidden": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 54, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX8AAAD8CAYAAACfF6SlAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3Xl4HMWZ+PFvdffco2N0W7dsy5d838bc4TAmQAgEMJuDJeBNFpJs2B8k2Q2EZZOFAEsCgSxLAmFzQSCnuQPhCgaMjS18S5Ytybrve+7p+v0xkixZsiXLIx9SfZ5nnpnpqa6qEebtmurqt4WUEkVRFGVy0U52BxRFUZQTTwV/RVGUSUgFf0VRlElIBX9FUZRJSAV/RVGUSUgFf0VRlElIBX9FUZRJSAV/RVGUSUgFf0VRlEnIONkdOJKUlBSZn59/sruhKIpyWvn444+bpZSpI5U7ZYN/fn4+W7ZsOdndUBRFOa0IISpHU05N+yiKokxCKvgriqJMQir4K4qiTEKn7Jy/oihKrIVCIaqrq/H7/Se7K8fNbreTnZ2NxWIZ0/4q+CuKMmlUV1cTFxdHfn4+QoiT3Z0xk1LS0tJCdXU1BQUFY6ojJtM+QoinhBCNQoidR/hcCCEeEUKUCSG2CyEWx6JdRVGUY+H3+0lOTj6tAz+AEILk5OTj+gUTqzn/p4E1R/n8EqCw97Ee+J8YtasoinJMTvfA3+d4v0dMpn2klO8KIfKPUuQK4Jcyes/ID4UQiUKIKVLKuli0P1DJ+x9w8NebCeZEuPQ734x19YqiKBPCiVrtkwVUDXhf3bttECHEeiHEFiHElqampjE1JDWNGXHzcdU4xtZTRVGUcfbqq68yc+ZMpk+fzn333Tfk84ceeog5c+Ywf/58PvWpT1FZOarrto7JiQr+w/0+GXLneCnlE1LKpVLKpampI16dPKyCxYtp8FWRZc9Hmurm9IqinFoikQi33HILr7zyCrt37+aZZ55h9+7dg8osWrSILVu2sH37dq6++mruuOOOmPfjRAX/aiBnwPtsoHY8GrJZLdT0VOIwXPTsaRyPJhRFUcbso48+Yvr06UydOhWr1cp1113HX/7yl0FlzjvvPJxOJwArV66kuro65v04UUs9NwC3CiGeBVYAHeMx39+nNlJKyFxF/et7mF6UPl7NKIpyGvuPF3axu7YzpnXOyYzne5cVHbVMTU0NOTmHxsLZ2dls2rTpiOWffPJJLrnkkpj1sU9Mgr8Q4hngXCBFCFENfA+wAEgpHwdeBtYCZYAX+MdYtHsknUlV1HjLyK4vRIZNhKEuZFYU5dQQXfcy2JFW7vz6179my5YtvPPOOzHvR6xW+6wb4XMJ3BKLtkajK62T6pJa8t1z8Je24ZiTfKKaVhTlNDHSCH28ZGdnU1V1aP1LdXU1mZmZQ8q98cYb/OAHP+Cdd97BZrPFvB8TckisO+No97YRiHjp2jIupxYURVHGZNmyZezbt4/y8nKCwSDPPvssl19++aAy27Zt45/+6Z/YsGEDaWlp49KPCRn8ncKNED6qekoI7G3HDEROdpcURVEAMAyDRx99lIsvvpjZs2dzzTXXUFRUxF133cWGDRsAuP322+nu7uZzn/scCxcuHHJwiEk/Yl7jKcAhEml11VHZbWV6/CL8u1twLhqfo6eiKMqxWrt2LWvXrh207Z577ul//cYbb4x7HybmyF9PoSy9mtaIxBfpwbuj+WR3SVEU5ZQyIYO/w5bOgYwwjkgcVd178Je2qqkfRVGUASZk8Lfa0/C6TRICEaq9pRCW+EtaT3a3FEVRThkTMvhj95AWMbEabTT7GwmJEL6daupHURSlz4QM/sLlITUSoS2+CmHJosZXhn9vKzKkpn4URVFgggZ/zZlMaiTC/tRWrKRQ0bYdGTTxl7ad7K4piqKcEiZk8Le4E0kLR9iRFsTjDdPoP4hpSHw7W0521xRFUUZM6fz4448zb948Fi5cyJlnnjkk62csTMjg77LZiItoHPSYpHpbkMJOq96Ib08LMmye7O4pijKJjSal8/XXX8+OHTsoLi7mjjvu4Lbbbot5PyZk8HfadBwhK1ITeFLCaEYO++o/Rvoj+Pe3n+zuKYoyiY0mpXN8fHz/656ennG59eSEvMLXbTOwRhxABG2aE0t5OtXtGyH9Uvw7W3DMTDrZXVQU5WR75dtQvyO2dWbMg0uGTuMMNNqUzo899hgPPfQQwWCQN998M7b9ZKKO/K0GRih6G8e2aSkkek1MIvgTAtGpH3WHL0VRTpLRpnS+5ZZb2L9/Pz/84Q/5/ve/H/N+TMiRv8umY0bigW5qcp2kdRygMSWemu59TPMVEazqwpYXP2I9iqJMYCOM0MfLaFM697nuuuv46le/GvN+TMiRv8tm4DXjcJuSaoePZK0Vzchl94H3QBP4d6tVP4qinByjSem8b9++/tcvvfQShYWFMe/HhBz5Oy067bhJC4dp8jaSXpiM3p2F17cTkWHBt7uFhEsKTnY3FUWZhAamdI5EItx44439KZ2XLl3K5ZdfzqOPPsobb7yBxWLB4/Hwf//3f7HvR8xrPAUYukaPiCM1EqGxpx7n/E+R8DcvTUCbtZnE2nhCTV4sqc6T3VVFUSahkVI6P/zww+Pehwk57QPgN+JJC0do8jbimD+f5I5qhJ7M/vqt0c93q0RviqJMXhM2+Aes8aRGwjT5W7HNLSKhcz+akceB0q0YGU58e9S8v6Iok9eEDf5BSyJpkQhhGaHLCJGWApqRgxkJEUwNE6zsJNIdPNndVBRFOSkmbPAP2zykhqNZPBu9jbjnzSYhpIHQqPPtBwn+vWrqR1GUyWnCBv+ILYG0SDT4N/macMybT1JLGZo+hdKyj9ATbPjUvL+iKJPUhA3+whHN6Q/Q5G3CsWA+ie370IxcGiv2Y0xzEyhrQ4ZUojdFUSafmAR/IcQaIUSJEKJMCPHtYT7PFUK8JYTYJoTYLoRYO1w9sWS32XCEbQA0eBuwFRbiCdagWXJBSjqsLcigSaC8Y7y7oiiKMshIKZ37/P73v0cIwZYtW2Leh+MO/kIIHXgMuASYA6wTQsw5rNh3geeklIuA64CfHm+7I3HZdHy4SBYW6nrqEIZBXNFMEk0NodmoaNiOsGhq1Y+iKCfUaFI6A3R1dfHII4+wYsWKcelHLEb+y4EyKeUBKWUQeBa44rAyEuhLppMA1Mag3aNy2QzapJssqVPbHW3OuWQxiQ070YwsKnZuwzY9MXp7x2ESLSmKooyH0aR0Brjzzju54447sNvt49KPWFzhmwVUDXhfDRx+qLob+KsQ4muAC7hguIqEEOuB9QC5ubnH1SmXVafNdDHFDLCnN/g7Fi/B8+xDiNRVdDS+jcwyiOwJEG70Ykl3HVd7iqKcXn740Q/Z27o3pnXOSprFt5Z/66hlRpPSedu2bVRVVfHpT3+aBx98MKZ97BOLkf9wdxk4fCi9DnhaSpkNrAV+JYQY0raU8gkp5VIp5dLU1NTj6pTTatCOmymhELU9tZjSxLFwIQldFeiW6IGlPlAJqCWfiqKcOCOldDZNk29+85v893//97j2IxYj/2ogZ8D7bIZO63wZWAMgpfxACGEHUoDGGLQ/LLfNoF26yfT7CNtsNHmbSHen4y7MJxGTZkscFfu2kZZ5Ib49rcSdkzNypYqiTBgjjdDHy0gpnbu6uti5cyfnnnsuAPX19Vx++eVs2LCBpUuXxqwfsRj5bwYKhRAFQggr0RO6Gw4rcxD4FIAQYjZgB5pi0PYROW3RzJ5Z/i4A6nrqAHAsXUJiw07Qcjm4oxjbTA/Byk5Mb2g8u6MoigKMnNI5ISGB5uZmKioqqKioYOXKlTEP/BCD4C+lDAO3Aq8Be4iu6tklhLhHCNH3jf4VuFkI8QnwDHCDHOezrC6rQbt0kR2KpnCo6a4BwLl4CYnNe9CMXPzdXfjivdGrfUvbxrM7iqIowOCUzrNnz+aaa67pT+m8YcPh4+Zx7EcsKpFSvgy8fNi2uwa83g2sjkVbo+WyGXTgJqM3xcPAFT8JHfvRjRsIAQfrdpHjzsW3pxXnwrQT2UVFUSapkVI6D/T222+PSx8m7BW+TqtOu3TjlJIkSxy1PdHgb6Sm4shOx6P5sToyqNixFfvMJPylbciIWvKpKMrkMGGDv8sWnfYByLR5+kf+0Dv1U78dU+ZQW7IXo8CF9IUJVnWerO4qiqKcUBM4+Ou0EQdApsU9OPgvXUJiXTHCyMOMhGkKV4MG/r1q3l9RlMlh4gZ/q0GHdAOQqdmp7Y6u9QdwLF5MQmcFFusUNN1C5Z5irHkJ+EvUen9FUSaHCRv8HRadDnqnfaRB0AzS6o8Gd2t+PpakBFL1diz2XCq2b8Mxy0OorodIR+BkdltRFOWEmLDBX9MEVqsNnx5HVjgMHFruKYTAuXwZCdVbiZjZtNVWE+5d6ONTo39FUSaBCRv8AZw2g04jmUx/D8CgeX/XylUkVm1Bs+QBUF27Bz3Rpub9FUUZdyOldH766adJTU1l4cKFLFy4kJ///Ocx78OEDv5um0G7nkSmtx04LPivWomrpw67LRHDFk/F9q3YZ3oIlLUjw+oGL4qijI/RpnS+9tprKS4upri4mJtuuinm/ZjQwd9p1WkVHpzdTSTaEgcFf0tODpbMDFIjdWhGHpU7tmErTEQGIwQq1A1eFEUZH6NN6TzeYnKF76nKZTVoDiRCVwOZ02ZR01PT/5kQAtfKVcRv+wg5ZR7Bnh200gSGwL+3Dft0z0nsuaIo463+v/6LwJ7YpnS2zZ5Fxr/921HLjCalM8Af/vAH3n33XWbMmMGPfvSjQfvEwoQe+btsOo0yEcI+shxpg0b+EJ368dRuQzNyQQgqd32MbWqiWvKpKMq4GSmlM8Bll11GRUUF27dv54ILLuBLX/pSzPsxoUf+TptBvRm9gdgUSxx/765DStn/h3auWIE90Ea8U6NL5lJevJX5l55Px4sHCLf4MJIdJ7P7iqKMo5FG6ONlpJTOAMnJyf2vb775Zr71rdinn57YI3+rTk04AYhe6OWP+PvX+gNY0tKwTp9GirccM5JNw4F9kBU9Hqosn4qijIeRUjoD1NXV9b/esGEDs2fPjnk/JnbwtxnUhKIj/yypAwyd+lmxkvjSv4OWD0BV9W70ZDv+EhX8FUWJvdGkdH7kkUcoKipiwYIFPPLIIzz99NOx70fMazyFuKwGlUE32CAz0pvauaeWeanzDpVZtZLEZ57DmPNlTJubiuKPyZx5JT2b65EhE2GZ0MdHRVFOgpFSOt97773ce++949qHCR3ZnDaddulC6jYyAz5g6MjfuXw5OhFSHT1YbPlUfLI1uuQzZBIoV0s+FUWZmCZ08HfbDEBgutJw97SSYEuguqt6UBk9Ph57URHJzTsIR3LwdXXSLprB0NSN3RVFmbAmdPB3WqOzWmFHKnTXkxeXR2Vn5ZBy7rPOIm7nG2hGPiCo2LkV+7QEddJXUZQJa0IHf5c1epI34EiF7kbyE/Ip7ywfUs59ztk4e+pxu2zY47IoL96CfWYS4WYf4Wbfie62oijKuJvYwd/Wu2zTlgpd9eTH59PobcQb8g4qZ587F8PjIS1cjSlzqSsrRfYt+VQXfCmKMgFN8OAfHfn3WJPA10q+OxuAis6KQeWEruM660ziS98FkQ9SUlW5EyPFgU8t+VQUZQKa0MG/b86/25ICQL4lelvHio6KIWXdZ59DQvVWdEs6FnscB7Zujmb5PNCBDEVOWJ8VRZn4RkrpDPDcc88xZ84cioqKuP7662Pehwkd/N290z4dehIAuVJHIIaM/AFcq8/AkCFSnT0Y9mlUfLIV6/QECJv4D6gln4qixMZoUjrv27ePe++9l40bN7Jr1y5+/OMfx7wfEzr4O3tP+LaJaIZOm7eNTHfmsCN/w+PBsWAByQ2fEA7lEPD20BKuRVjUkk9FUWJnNCmdf/azn3HLLbfg8URjV1paWsz7EZMrfIUQa4CHAR34uZRyyO8YIcQ1wN2ABD6RUsb+d8xh+k74tmi96Zm7oyd9hxv5Q3TVT8Ljv0Fb8e8ITefAJ1somrYCf0nboIRwiqKc/v7+XCnNVd0xrTMlx81Z18w4apnRpHQuLS0FYPXq1UQiEe6++27WrFkT074e98hfCKEDjwGXAHOAdUKIOYeVKQS+A6yWUhYB/3K87Y6GzdCwGhoN4Wh+n77lnhWdFcOmVXWffTZOfzOJbh17XH7/vH+k1a+WfCqKEhOjSekcDofZt28fb7/9Ns888ww33XQT7e3tMe1HLEb+y4EyKeUBACHEs8AVwMBJrJuBx6SUbQBSysYYtDsiIQRJTivNPhOcydHlnhkr8YV9NHobSXelDypvmz0bIzWVNP8BWkM5tNa8TTA5erLXX9KGJdV5IrqtKMoJMNIIfbyMJqVzdnY2K1euxGKxUFBQwMyZM9m3bx/Lli2LWT9iMeefBVQNeF/du22gGcAMIcRGIcSHvdNEJ4THZaXNGwR3Rv/IH4Yu94TowcJ97jkk7HgNYRREy+0vxkh1qPX+iqLExGhSOn/mM5/hrbfeAqC5uZnS0lKmTp0a037EIvgPNxF++O8aAygEzgXWAT8XQiQOqUiI9UKILUKILU1NTTHoGiS5LLT2BMGd1j/nD8Mv9wSIu/BC4hr34nR6sDpTKN+2GfvMJALlHZhBteRTUZTjM5qUzhdffDHJycnMmTOH8847jwceeGDQDV5i0o8Y1FENDLy5ZDZQO0yZD6WUIaBcCFFC9GCweWAhKeUTwBMAS5cuHToxNgYep5Xa9k6YkgEtZaQ703EYjiOe9HWuXInucpIu6ijX8qnaVYxxhRPekwT2t+OYHdv/AIqiTD4jpXQWQvDQQw/x0EMPjVsfYjHy3wwUCiEKhBBW4Dpgw2Fl/gycByCESCE6DXQgBm2PKMll7R35p0N3AwLIi88bNscPgGa14j73XBL3vAkUEAmHqWnbh7Bq6gYviqJMGMcd/KWUYeBW4DVgD/CclHKXEOIeIUTfRNZrQIsQYjfwFnC7lLLleNsejSSXlQ5fiIgrDSJB8LVFl3seYdoHolM/CdUfo9uyMKxO9hd/hG1a9Mbuw52pVxRFOd3E5CIvKeXLUsoZUsppUsof9G67S0q5ofe1lFLeJqWcI6WcJ6V8NhbtjkaSywpAjzWa4oHuBvIT8qntriUYCQ67j/usM7HokjRrB4ZtGge2bsZamEikLUC4SS35VBTl9Dehr/CF6Jw/QKfee6FXb3ZPieRg58Fh99FcLlxnnklyxftEzHwCPd2069ET0OpqX0VRJoIJH/z7Rv4tou8q36Mv9+wTd+GFJO1/F92ah9AM9u/5CCPdqZZ8KooyIUz44N838m+UvStLBy73PFrwP+9crKaPVEcAqzOfsi2bolk+Kzox/eFx7rWiKMr4mvDBv2/k3xiwgMUJXfW4LC7SHGmUdwy/4gdAT0zEuXwZKVUfEjHz6WxqIOAJQkQSKIvtZdaKokwuI6V0/uY3v8nChQtZuHAhM2bMIDFxyGVRx23CB/9EpwWANm8IEnKgPTrPP90zndK20qPuG3/JJSSVvoVumQYIDlRvQ9h1fGreX1GUMRpNSucf/ehHFBcXU1xczNe+9jU++9nPxrwfEz742y06LqtOa08IPHnQHr2B+8ykmexv308oEjrivvEXXYRNeklxmFgcmezfugl7oac/y6eiKMqxGk1K54GeeeYZ1q1bF/N+xCSl86muP79PYh4c/BCkZKZnJiEzxIGOA8xMmjnsfnpiIu6zzya58n0aEgtoOPAe8iIDc0eQUF0P1kz3Cf4miqLEyltPP0FjZWyvNU3Lm8p5N6w/apnRpHTuU1lZSXl5Oeeff35M+wmTYOQPA67y9eRBoBN8bcxKmgUw4tRPwqcvJbniPTRLIQBV7dGfZ2rVj6IoYzGalM59nn32Wa6++mp0XY95PybHyN/ZO/L35Ec3tFeSlzEPq2Zlb+teLpt22RH3dZ93Hi7tuyRZdVrs6ezdtpHMrHX497YRf17uifkCiqLE3Egj9PEympTOfZ599lkee+yxcenH5Br5J+ZFN7RVYmgGhZ5CStpKjrqvZrcTd+GFpFRuRDKV2tI9aHkOggc7ifQc+XyBoijKcEaT0hmgpKSEtrY2Vq1aNS79mBTB3+O00tY37QPQVgFET/qWtJaMePI2/rLLSK16H806A6SkwV8BEgKlKtGboijHZjQpnSF6ove6664bt9vHToppnySXhZ5gBL/uxm5PPLTixzOTP+7747B39RrItXIFbrcgxYBGazJ79r7Latfl+Pa24lwU+xsrK4oysY2U0hng7rvvHtc+TI6Rf++FXu3eUHTev+3Qck9gxKkfYRjEX7qW1LI3kdp0qvbsxJjmji75jKgln4qinH4mRfBP6k3x0L/ip3fkP8MTvYdnSevRgz9A4lVXkVq/BcMyHWmatFCP9IcJHuwcv44riqKMk8kR/HtH/v1r/dsPgmkSZ40jy53F3ta9I9ZhnzmT+FkFpIS8aEYiJeUbQRPqal9FUU5Lkyr494/8I0HoqgNgVtKsEdf690m8+mrSKt5F6NMp37UVS65LpXhWFOW0NCmCv2fgyH/AWn+InvSt7KzEG/KOWE/8pWtJ6ynFsE7HjETosLUSbvASbvWPV9cVRVHGxaQI/omOaHK36Fr//OjGASd9JZJ97ftGrEd3u0m6+HzSOhrR9ERKqqKXZKurfRVFOd1MiuBv6BoJDkt0rX9iDiAGrfWH0Z30hejUT0bthwjLTEp3f4DmseLbo4K/oiijd+ONN5KWlsbcuXNPWh8mRfCH3qt8vSEwbBCf2T/tk+nKJM4SN+rg71i4kCkeP3YjD6Skx9VF4EA7ZjAynt1XFGUCueGGG3j11VdPah8mTfD3OHtH/hBd8dM77SOEYFbyLHa27BxVPUIIkj93FZn1JQg9mX21H0NYEtinrvZVFGV0zj77bJKSkk5qHybFFb4QHfnXtveemPXkQfm7/Z8tTF3IUzufwhvy4rQ4R6wr4coryXxiHeUpZ7J3/0bmzzoT3+5WHEUp49V9RVFirP2F/QRre2JapzXTReJl02Ja53iZRCP/3syeEB35d9ZCOADAwrSFRGSEXS27RlWXHh9P9sUriQ/HITHxxXnx721FmupqX0VRTg+TauTf2hNESonw5AMS2qsgZTrzU+YD8EnTJyzLWDaq+jyf/zzZN91HW0Y6ZfXFzDNWEqzqwpYXP35fQlGUmDldRujjJSYjfyHEGiFEiRCiTAjx7aOUu1oIIYUQS2PR7rHwuKwEwia+UORQds/2CgAS7YkUJBRQ3Fg86vpsUwsomG5BtxSyr+ZD0MC/p2Uceq4oihJ7xx38hRA68BhwCTAHWCeEmDNMuTjg68Dw9ysbZ4Py+wzI699nUdoiipuKj+nevBlfuJa0rgghM4jP6VNLPhVFGZV169axatUqSkpKyM7O5sknnzzhfYjFyH85UCalPCClDALPAlcMU+4/gfuBk3I5rGdgioe4KaBb+9f6Q/Skb0egg/LO8lHX6Vq9mgJZg2bkUtbwSfRq3xZfrLuuKMoE88wzz1BXV0coFKK6upovf/nLJ7wPsQj+WUDVgPfVvdv6CSEWATlSyhdj0N6YJLkGXOWradHRf+uhmzcvSFsAwCeNn4y6TqFpzPjcWdhFDpWt0SkjNfpXFOV0EIvgP9xtZvrnToQQGvAj4F9HrEiI9UKILUKILU1NTTHo2iEe54D8PgBps6Fxd//n+fH5JNgSKG4a/bw/QOJnryS/q5WesBef7sW/W837K4py6otF8K8Gcga8zwZqB7yPA+YCbwshKoCVwIbhTvpKKZ+QUi6VUi5NTU2NQdcOSXbZAGjtu+9uehG0lkMwus5XExoLUhcc00lfiN7jd97FM9Et06lo3UmgogPTq+7tqyinqmM5r3cqO97vEYvgvxkoFEIUCCGswHVA/40opZQdUsoUKWW+lDIf+BC4XEq5JQZtj1qc3UDXxKGrfNPmABKaDuXyX5S2iAMdB+gIdBxT3VO+eC3JQRs1PSVggr9EXe2rKKciu91OS0vLaX8AkFLS0tKC3W4fcx3Hvc5fShkWQtwKvAbowFNSyl1CiHuALVLKDUev4cTQNEGyy0pDZ+/55vSi6HPDbshaAsCC1N55/6ZPODv77FHXrbvdLF6Yx6vbywhIH75dzerevopyCsrOzqa6uppYTyufDHa7nezs7DHvH5OLvKSULwMvH7btriOUPTcWbY5FtsdBdVvvahxPPhiOQfP+c1Pmogud4sbiYwr+ADO+chVv37qPqq692EpcyFAEYdFj2HtFUY6XxWKhoKDgZHfjlDBp0jsA5CQ5qW7vvWmLpkPaLGg4lNLBYTiYlTSLbY3bjrluS1ISMzOzqPaWQsjEX9Yeq24riqLE3KQK/tkeB7XtfsIRM7ohrWjQyB9gecZyipuKR3Vnr8Od8S/X0hwMEjQD+HapVT+Kopy6JlXwz/E4iZiS+v55/znQ0wTdh+b/zsg6g7AZZnP95mOu35mdwZSkPOq8ZXQX16lEb4qinLImVfDP9kTTNVe19s77p/VmoWg8NPWzOG0xDsPBxtqNY2rjU9/4EtU95WhhjWBF53H1V1EUZbxMquCfk+QAoLqtd0pn4IqfXlbdyrKMZWysGVvwTynMxmtoRMwwLa9tP67+KoqijJdJFfynJDgQAqr6Vvy408CZMmjkD3BG5hkc7DpIVVfVMLWMbPm6z9Lgr8C7vxvTNI+324qiKDE3qYK/1dCYEm+nunXAydz0OYNG/gCrM1cD8H7N+2Nqp+j8ZdQF6rAbbtr//NaY+6soijJeJlXwB8hOch5a6w/RFT9Ne2HACD0vPo8sd9aY5/0B4hZOxZQmFS/vxQwEjqfLiqIoMTf5gr/HQVXbYSP/kBfaDqVyFkKwOnM1m+o2EYqMLU/Pqi9cRpO/GkdCPs1P/uJ4u60oihJTky7453ic1Hf6CYYHrPWHIev9z8g6A2/Ye8xZPvs43C56EiPEWeLZ99J2QnV1x9NtRVGUmJp0wT/b40BKqG3vW+45CxBD5v1XZKzAEAbv145t3h9g5nWrkVLiz1lCw/0PHkevFUXPnajqAAAgAElEQVRRYmvSBf+cpOha//55f6sLkgqgYeegcm6rm/mp83mv5r0xt5U+exrtsoV0ZzIVH5TRs+mjMdelKIoSS5Mu+Gd7omv9B837Zy6G6s1wWJrX83PPZ2/rXg52Hhxze475KSRak6medQ51d92J6T8pd7FUFEUZZNIF/4x4O4YmqBq43DNvFXTVDbqnL8DF+RcD8FrFa2NuL29tNF10gt1GU7uF5kcfHXNdiqIosTLpgr+ha0xJtA9e7pl7RvT54AeDyma4MliYupBXK14dc3uWJCcBV4BsRzpli66i+alf4Nuxc+QdFUVRxtGkC/4QXfEzaNondRbYE6Fy6MndNQVrKG0r5UDHgSGfjZZnZT4eWzoBfw3t+Sup++53kSF1q0dFUU6eSRv8B438NQ1yVw0Z+QNcmHchAsFr5WOf+olfFr3bTqZNY3/RlfhLSmn+n/8Zc32KoijHa1IG/2yPg6auAP5Q5NDGvFXQUgbdjYPKpjnTWJK+hFcqXhnzfT+NRDvaFBu5rlm0NBbTvXY9zY//L94tJ/Q2xoqiKP0mZfAfstwTjjjvD7Amfw3lHeWUtpWOuc34ZdkkWlOJo5q9toXo2bnU3HEHkU6V9llRlBNvUgb/YZd7TlkQvadv5dDgf0HeBWhCO65VP455KSAg21FAW/02um+4m3BjE3Xf+96Yf1EoiqKM1aQM/sOO/A0rZC+Fg0NP+iY7klmRsYJXyl/BlGNL0azHWbFNS2Rq4nyEuY2Pt/iI/8rX6XrlVdp///sx1akoijJWkzL4p7ptWA1tcGpngLwzoH4H+IdOxVw+/XKqu6vZVLdpzO06F6TiwI3btOLv2Mv+5LNxnbGKhnv+E9+OHWOuV1EU5VhNyuCvaYLcJCf7m3oGf5C7EqQJ1UPTMFyUdxEem4dn9z475nYdRcmgC2akL0PXtrLj3Rrst38fIzWV6q99nXCLuum7oignxqQM/gBzM+PZUdM+eGP2chD6sPP+Vt3KZws/y9vVb1PXPbYMnZrTgn2Gh1zXbHydtWiU88FrDWQ+8jCRtjZqvnkbMhweU92KoijHIibBXwixRghRIoQoE0J8e5jPbxNC7BZCbBdC/E0IkReLdo/H/OxEGjoDNHQOyLVjc8OU+VAxfDK3a2ZeA8Dzpc+PuV3nglS0gKBgygJ0bQtVe1qpD6WT8R934/3oIxruvU+dAFYUZdwdd/AXQujAY8AlwBxgnRBizmHFtgFLpZTzgd8D9x9vu8drfnYCANurOwZ/UHhRdNqnu2nIPpnuTM7OPps/7PsDwUhwTO3aZycjLBoLpl9AV3MV7oQ63v1dKc6LP03SP/4jbb/5Da2/eHpMdSuKooxWLEb+y4EyKeUBKWUQeBa4YmABKeVbUsq+s6sfAtkxaPe4FGUmoGuC7dWHTf3Mviw671/y0rD7XTfzOlr9rbxe+fqY2tVsOo65KdibbXjSMiHyEd3tfj74837Sbv9/xK1ZQ+P999P5yitjql9RFGU0YhH8s4CqAe+re7cdyZeBkx7ZHFadwjQ3nxw+8k+fC54C2L1h2P1WZa4iNy73uE78OhenIf0RVq26hpbqA+TN7mbnuzXUl3eR+cP7cCxZQu0d38K7efOY21AURTmaWAR/Mcy2YSethRCfB5YCDxzh8/VCiC1CiC1NTUOnXWJtQXYi26vbB8+xCxEd/Ze/A772IftoQmPdrHUUNxXzccPHY2rXNi0RLd5KSjCDhLR0WqvfwJVg5a1f70VqFnIeexRLTg5VX/kqvu3bx/r1FEVRjigWwb8ayBnwPhuoPbyQEOIC4N+By6WUgeEqklI+IaVcKqVcmpqaGoOuHd38nATavSGqWn2DP5hzBZhhKB0+lfNVM64ixZHCY8WPjaldoQmci9IIlLZz5hWfp/lgOVMXtNFW18OWVyvQExPJ/cVT6ElJHLzpZvy7d49cqaIoyjGIRfDfDBQKIQqEEFbgOmDQnIkQYhHwv0QDf+MwdZwU87MSAdh++JLPzMUQnwV7Xhh2P4fh4KZ5N7G5fvOYL/pyLUoDU5JlLyQ1r4DS9/9M4bIUPn6lkvoDHVjS08l7+hdobhcHb/wy/pKx5xVSFEU53HEHfyllGLgVeA3YAzwnpdwlhLhHCHF5b7EHADfwvBCiWAgx/IT6CTYzIw6rrg1d8aNpMOvTUPYGBLqH3ffqGVeT5kzjseLHxrQ005LhwpLpwlvcxNnX30BHYwOJqQdwJ9p4/Re7CfrDWLKyyHv6aYTNxsEvflFdBawoSszEZJ2/lPJlKeUMKeU0KeUPerfdJaXc0Pv6AillupRyYe/j8qPXeGJYDY3ZmfF8UjV0bp85l0PYHz0ADMOm21g/bz3bGrfxfu3QfECj4VycTqi6m8wps8idO58tG57jnHX5dDb7eO/5fdE+5uaS96tfosXFcfBLN9Dz4djTSyiKovSZtFf49lmQncDOmg4i5mGj99xV4EyBPUf+kXJl4ZVMcU0Z8+jfuSAVNPBubeSsdTfg6+qkatcbLL44jz0b69i/LTpDZs3NJe83v8GSlUnV+vV0vj62ZaaKoih9Jn3wn5eVQE8wQnnzYdM7mh4d/e99Gbytw+5r1a18dcFX2dG8gxcPvHjMbetxVuyzkvF+3EB6wXRmrT6HLS/8kRnL7KTmxvHmL/fS3hi9PMKSnkbuL3+JbfYsar7+DVqefEpdCawoyphN+uC/ICd60veTqo6hHy67CcI+2Pp/R9z/iulXMD91Pg9sfoB2/zDTRyNwLUvH7A7h39PKOZ+/Ec2w8O5vnmTN+rkIDV79352EgtE7jhkeD3lPP03cRRfR+MAD1N15JzI4tiuNFUWZ3CZ98J+W6sZp1Yde6QuQXgQF58BHP4PI8Ddc14TGXSvvojPYyUMfP3TM7dtnJKHFW+nZXI87KZlVV13HgY8/orlqJxfeWERLbTfv/Kakf5SvORxk/eghkr/6FTp+/wcqb7yRUOMps4BKUZTTxKQP/romWJSbyMb9R0invPKr0Flz1Ln/mUkz+WLRF/lT2Z/YUn9s9+UVusC1JB1/aRvh9gCL116OJzObt55+gqwZ8Sy7tICSTfXsfKfm0D6aRto3vkHmgw/i37Wb8s9epU4EK4pyTCZ98AdYU5RBWWM3pQ1dQz8svDia7uHD/zlqHV+Z/xWy3Fnc8+E9x5z0zbUsAyR4P25ANyycf8N62uvr+PjFP7FsbT7585L5+3P7qNw1+ACV8OlLKXjud+jx8Ry88UaafvpTlRJaUZRRUcEfuHhuBpqAF7cPk6df06Kj/+rNUHXkXDtOi5Pvrvwu5R3l3L/52JKWGkl2bNMT6dlcjzQl+QsWU7jiDD74wzO01ddw4ZeLSM5y8doTO2muHnyAshUWUvD8c8SvXUvzIz+h8vNfIFhZeUztK4oy+ajgD6TF2VlekMRL22uHX0Gz8HqwxcOHPz1qPWdmnckNRTfwu5Lf8dKB4bOCHolrWQaR9gCB/dFzD5+68atYrDZee/wRLFaNS/95AVaHwYuPbqe7zT9oX83lIvOB+8l88EECBw5w4DNX0vrb3yLNsd1vWFGUiU8F/16Xzs9kf1MPpQ3DXNFri4MlX4Ldf47e4/covr746yxOW8x/fPAfHGg/MOr2HUXJaE6D7g+jvz5ciR7O/dLN1JbspvivL+H22Pj0rQsI+sO88JNP8HUPnloSQpDw6UuZ+sIGnEuW0HDPf1J5/T+otBCKogxLBf9ea4qiUz8vbR+Sky7qrH8FhwdevgOOsr7eolm4/+z7cRgObnv7NnpCPUcsO5AwNFzLMvDvbiHcHh3Zzzn7fPIXLuHvv/0/OhobSMl2s/Yr8+ho9LHh4WIC3qErkCzp6eT87Akyf3gfwcpKyj/7WRoeeIBI9/BpKhRFmZxU8O+VGmdjRUEyL+6oG37qx+GBT90FB9+HHb8/al3prnR+ePYPqeis4BtvfoNAZNgkpkO4VkwBoGdTPRAdzV940y0gBK89/jCmGSF7VhJr/mkurbU9vPCTTwj6h57gFUKQcMUVTH35JRKu/AytTz7F/ovX0Pbs79QJYUVRABX8B7l0/hQONPVQMtyqH4BFX4ApC+H1OyFwhDK9Vk5ZyX+u/k821W/i9nduJ2yOHHSNJDv22cn0fFSPDEfn6+NT0zjvSzdTtWs7mzf8EYD8eSlcfNNcGiu7ePHRTwj4hq/b8HjI/P73yX/+eawF+dTffTflV15J51//qs4HKMokp4L/AGvm9k39DLPqB6IpH9Y+CF118O6w96MZ5LJpl/Gd5d/hraq3uGvjXZhy5IDrXjUFsyeEd0dz/7a5513IjFVn8f5zv6ZuXwkAUxelcuGNc2g40MmfH9qKt/PIy0sd8+aS96tfkfXIw8hQmJqvf4Pyq66m629/UwcBRZmkVPAfIMVt44xpKTy/pRp/KDJ8oZxlsPAf4IPHoPKDEeu8fvb13LrwVl448ALffvfbI14DYJuWiJHqoOf9Q+cehBBcePMtuJOSeemR+wl4o/l+Cpems/aW+bTXe/njgx/T2ew7UrUIIYi/6CKmvvgCmT+8D7Onh+pbbqX8iito/9OfVZoIRZlkVPA/zK3nT6e+08/T71ccudDF/wWJefDcF6HzCCeIB1g/fz3fXPJNXql4hZv/ejMdgWHyCPUSmsC9cgrBqi6CVYemluwuN2u/djudzU28/rNH+89L5BUlc/m/LMLfHeIPD3xMQ3nnUfsiDIOEK65g2ksvkvnD+wBB3Xe+Q9mFF9H8+OOEm5uPur+iKBODCv6HWTk1mfNnpfHTt8po9x5hNOxIhOt+CyEv/O4LED76CV0hBDfOvZH7z76fHc07+PzLn6eio+KI5Z1L0hFWne73Bx9YsmbOZvW1X6Dk/XfZ8uKf+rdPmZbAlf+6GN3Q+NN/b6X0o/oRv6ewWEi44goKNvyFnJ89gW3aVJp+/DD7zjufmtv+lZ4PPlBTQooygangP4xvrZlFdyDMY2+VHblQ2iy48nGo2QIv3XbU5Z99Lim4hCcufIK2QBvXvHgNf9z3x2FXFml2A9eydLyfNBFuH3xgWX7F1cxYeSZ//83TVBQfuoF8cpabz317KekF8bz+1G4++FMZZmTk4C2EwH3WWeQ+9RRTX36ZpOvX0f3eexz8xxspu+ACGh9+mMCB0V+voCjK6UGcqjnhly5dKrdsObYkabF0+/Of8JfiWt78f+eQ7XEeueCbP4B374elX4a1D0RPCo+gvqee7773XTbVb+LCvAu5c+WdeOyeQWXCbX7qH9iMe3UWiZdOHfRZyO/nmTv/H50tTfzDf/0IT0Zm/2eRsMnff1fKrr/XMmVaAhfcOIf4ZMcxfXfT76frb3+j409/pmfjRpAS28yZxF+yhriLLsI2derIlSiKclIIIT6WUi4dsZwK/sOr6/Bx7gNvc1FRBo9ctxAhxPAFpYQ37oaNP4aiz8KV/wuGdcT6TWny9K6n+cnWn+CwOLhl4S1cO/NaDM3oL9PyzF78e1uZ8p3laHZj0P4djfX8+t9uwxEXz7p77scRFz/o89KP6nn7tyVomuDcf5jF9CVpx/w3AAg1NNL12mt0vvwyvuJiAKz5+bjPPx/3OefgXLQQYR35+yqKcmKo4B8DP36jlB+/sY9/WzuL9WdPO3rhjQ/D63fBtPPhqifBmTSqNsrayrhv831sqtvE9MTpfGPxNzgn+xyEEARrumn8yTYSLskn7pycIftW797J7//rTtLypvK5O3+AxW4f9HlHk5e//nwXjZVdTF+SxlnXzsAZP/ZAHaqro+utt+j+25v0fPQRhEJoTifOFStwnXEGzuXLsRVOR2hqNlFRThYV/GPANCVfe3YbL22v4yfrFnHZgsyj77D1V/Div0Tv/Xv5T2DGRaNqR0rJmwff5IEtD1DTXcMMzwxumncTF+VdROtTuwk1eplyxzKEMTSo7tv8AS/8973kL1jEFbffiW4M/oUQCZts+2slm1+uwGLVWX11IbNWZRz5l8woRbq78W7aRPd779Hz3kZCVVUA6ImJOJYuwbloEY5Fi7DPmYN22EFJUZTxo4J/jPhDEb745EcUV7Xzqy8vZ8XU5KPvULcd/vRP0Lg7ekXwp74H7tRRtRUyQ7xS/go/3/FzyjvKSXems979Bc54Zyqeq2fgWpo+7H7b//Yarz/xE2afeS5rbvkm2jDnHVrrenj713up299BekE8q68uZMq0hFH1a1R9r6mhZ/NmvJs+wrd166G00oaBrbAQx9wi7EVF2GfNwjZjBprzKOdRFEUZMxX8Y6jdG+Sq/3mfug4/d356Dtctyzn6yDkcgLfvjU4FGXZYfjOc8Q1wjXDg6GVKk7cOvsXz+57n/er3+Un5d4jX3Oy+upPz885ninvKkH02/ek53nv2l8xafQ6X3HIbmj70ACBNyd4P6/nwL/vxdgSZviSN5ZcV4MlwjfpvMVrh1lZ827bh274D/86d+HbuxOzovb5BCCy5OdimF2KbPh3b9GlY8wuwFhSgu2Pfl4lOhk0i3SHM7iCRnhDSG8b0hTH9YWQwghmIIEMmhE1kRCIjMnquShJ9FgI0gRCAoSF0gTC06MOqISw6wqaj2QY8Oww0u4HmiD6G+1WqnBwq+MdYQ6ef254rZmNZCxfMTue+q+aR4rYdfafmffDO/bDjebA4Ye6VsOB6yDsj+j/cKNR117Hp7TdZ+l4u92f+grcSNjPTM5MVU1awYsoKlqQvwWWJBsyP/vJ7/v7bpylccQaXfv12dMMybJ1Bf5htrx+k+PWDhEMmhUvSWLI2n+RM9zH9TY6FlJJwbS3+vXvx791LoKSUQFlZ9BdC5NDV1EZqKpbcXKw5OVhysrFkZWHNysKSmYmRno44bFprMpBSEukIEm7yEm72EW7xE271E+kIEGkPYPYMf39pAHSBsOpoFg0svYFd06KLvPv+DUoJZrQdIhIZNpERExmKPoiMHCOEVUdzGmguS/TZaUF3WdDcFjSXBd1tRXNb0N0WNLcVzTbyqjhlbFTwHwemKXlqYzn3v1aCw6Jz/YpcvrAyj8zEEZZSNpXC+w/Drj9DsDt6dfCsS6Mnh/NWg/XoUyDSlDQ+so1gMMBfL97N+/Xvs61xGyEzhEAwLXEa81LmUZRchL24ibI/vcrUxcv49L98C4vtyPPt3s4gn/ztINvfriEciJBblMT883LInZOE0I7vnMBomcEgwYqK3kclwYoKQgcPEqyqItzQMLiwpmGkpGBMycCSloaRmoaRloaRmoKenIyRkoqR5EFPSjptzzPIsEmowUuwpotQTTehuh5CDV5k4NABUlg0dI8dw2NDT7Chx1vR4qzocdZDwbd3ZB6LEbmMmMhABDMYiT77I9FfFb7eXxjeMKY3FH3dE8L0hon0hDB7QoP6PZCwamhua/RgEBd91uOs/a+jz1b0OAvCog4Ux+KEBn8hxBrgYUAHfi6lvO+wz23AL4ElQAtwrZSy4mh1norBv09JfRc/er2Uv+6uRwjBp2al8anZaayennL0awKCPbDnBdj+HFRuhLAfdCtkzIfMRdFH2mxIng72wUs3fTubafn1HjzXzsS1KA1/2M+2xm0UNxazo3kHO5p30B6I3gVsZqWblbuS8KVYiHxmDtkZU8lwZTDFNYUMVwYpjhRs+qFfLb7uIDvermHXuzV4O4MkpDqYvXoKM5ZnEJd08oKoGQgQqq2NPmpqCNfXE6pvIFxfR6ixkXBT86GppMMIpxMjMRG9/5GAFh+PnpCIHh+HFheHHh+P5o5Dd7vQ3G40lyv6cDpP2C8M0x8mUNFJsLyDQGUnwepu6M3oKuw6likuLOkuLBlOjFQnlhQHWrz1uE/YnygyZPYfCCLdQcyuEGZPkEhX7/vuEJGuIGZXENM7fHZaYdOjB4a+A4TLcugA4Rrw7LYgbPpp87cZLycs+AshdKAUuBCoBjYD66SUuweU+WdgvpTyK0KI64ArpZTXHq3eUzn496lu8/KrDyv549YamrqiV+LmJDmYnRFPYbqbaalupiQ4SI+3kRZvx2Ud8A8z5IPK9+HAW1CzFeo+if4q6ONOj/5CiM+E+CykK43G9+YhpUH6P9gRjgSwJ0TvMmZxIoWgvqeesvYyytrLOLhtK45XywlaTP66pJ62+MFTA3HWOJLtySTZk0i0JZJoTyReT8BVnYHYlUS41gpI4gp00uc5yZwbR6LHjV23YzNs2HU7Fs1y0v9HM/1+ws0tRFqaCbe0EG5pIdLaRqS1lUh7G+H2diLt7ZjtHUQ6O4l0dg6aZjoSYbOhORwIpwPN6USzO6LvHXY0m/3Qs92OZrMibHaEzYawWtBsNoTVirBYos99ry1WMAwiHYJQg0moLky4ORSde9fAkm7Hku3CmhOHNcuNkeKM7jdJls7KiNl/MOg/h9EVPWBEuqMHDLM7+pk8QhpzdBE9EAx89E1D9T73vRd2vfechQWhT5wDxokM/quAu6WUF/e+/w6AlPLeAWVe6y3zgRDCAOqBVHmUxk+H4N9HSsm+xm7e29fM5opW9jV2U9HcQ9gc/PUMTRDvsBBnN3BaDewWDYdFx2poWDVJtllLVriK9GA1aaEqkkL1JISaSAg1YpUBfJFVtIT+HY/lQVz624PqDmk2gpqDiLAS0myENSstPjtvlyUSjAjypgaIZIRp1CXtmqRNM2nTInQJky4RoUtE8IoIQRHtc7w/mRlNyyhsWkpCIBWJSUNcBQcT91CdUEKTuwopTAwpsKBhIDCkho4Y8Ii+1xDoUiAAHQ0BaL3b+14LQMhD26IPMeiZAdv7Xh9pGwO2DdouJZaQxBYwsfpNrEGJJSixBkwsQRNrQGIJmlhCEiNkYglKjJDEEjYxQtHXelhi9D70sEQPmxhHO57oVoy0IozMRehpc9FsbqQ0MdvKCTftJdK0l0hbORwh46sETA2kJjAFSA2kEEgNzN5nKXq3CTAFRDSNsC6IaFr0tSYwNY2IJogIDVMITE0QIbrdFH0PgSk0JAIptGh7iGjdiOg5YgSI3vPFA/rY94eWh/X9SN9p7KL/tmy6E5vuxK47sekubFr0tVV3YNOc2HQHNt2OVXNi04/+CzZsBgmaAUJmkJAZJCyDA16HCPc/h3qfw5gyTFiGiMgIETMcfZYRIkQwZRiz91sO/Zd56Hsc/rpvi241+OennxzbX+cEBv+rgTVSypt6338BWCGlvHVAmZ29Zap73+/vLdN8WF3rgfUAubm5Syr7lguehoJhk4OtXho7/TR0+WnsDNDhC9HpD9HpC+MLRfCHIviCEYIRk2DYJBgxMU1J2JRETIkpJaYEM2Jilz48spO7Auk4EPzG+h5OenDhw44fJ37sBLARwkYQGyEshJHhCGU1Tjp8FqbEe5mb0YpNi6BhYhB91pBomAgkYSReHXxC4tPAJ6ArnEWnbzHenkVEgvkASM2L6dhP2LGfgPMAAUclQcNHSAjCAiJAuPe1CUQERBC9z2ASDWQm9AcYk+hnfQHEhP5gY/admzzsAb2BrzcwwdDAMnC7PMoA71gD1cDtQkqMCFjC0Yc9bLDIN4/lvmUUheZhxUoP3ZRqO9nHDirNPYTMHowIGKYFI2LFkFaMiBXNtGBIK3rEgmYa0YfUEKaGkH3RtW/FjgQiSGn2PkeAsdytTQdh9D5HD919r0Hr3ab1vu597n/fd/jV6DsUg+jffugBwwW80Tv8sNL3/Qe8Rw77uUBi0SxYNStWzcCmWbFoBlZhwdL72iIsWDULhmbBolmwCCP6uvd5LCIygtn/MDExo8/SRHLoWUoZfUYipaQ71MU5P10/pjZHG/xjMbE53H/Bw/9/GU0ZpJRPAE9AdOR//F07eayGxvQ0N9PTYruCxl/aRvNTO/n3C9cTd2bWqPY524zw4R+e5cM//A6/q5CLvvJ1smcVjal9X3eQ6r1tVJe0UVeWSlvdPPrOHriTbKRkx5GU6SIxzUliupOEVAeOuJM/PXS8pJREQiahQISgP0IoECHkD0efA4fe0+jFWtONo9mHZkrCuqDdbaHFotFq2gkFz8TjnYfb10HI30k42IWUPWD2IKUX2fuM6QVCvY/ePjBwhG1DCDtCt6PpdnTdhmbY0Qw7hsWObrGhW3pfW20YFhsWqw3dasewWrHYbBhWG4bVimG1ohsamqah6WLAI/peaPR/JjSB1vsQGtFtIrpdDNjet037/+3de3BcV33A8e/v7nu1D63e1sOvRDF2HENMEpJA0rxJHSDDK+OGBNIy0ELL0Jbp0Jl2hg4dZhhapikzLTSFUJqBxA3TQoZH0iRgEgJmcOImxI7tOLZs62FJ1q60q8e+7v31j7tRZEeW1pa00mrPZ2Znd7VXu+cn3f3dc8859xzLPUuwSq8D0z9zt2F635j589ePHYLMOI6UzuWEiu9P6ihasNGc43Z8593hs9P3hdJ90UELihZttKju8NrSENvpkVSOglMacuvo9HMcN/Gj4GtZ+utgFiP59wIz5x7oBM6e5P71bXpLzT5xILkIn11zgpckCHTXk376BHXbW7DC89dILMvDtR/+CF1bLuPxr9/Pri98nm233M51d99HsO78Dk6hiJ/uK1rpLl1wlpssMHgszfDJDCO945zuHefEyyM4M5q8PD6LaEOQuvoA4ZifcMxPKOojEPYRrPMRCHnxBT34Au7N47PweC08PmtGonnzl11LZ0dqK07pZtsOdsHBLrq3YsF9Xiw4FPM2xfwb94WcTSFvUyzdv5HE7TOTeummzuz1kYDAWr/FWr9FxCMUHIfeQoa+4iAjhQEcJ41TTFMsjGHnx3lTvUcEfyhKIBwnGGkiFI0TisQJxeOEojHq6usJx+OE4zEiiXrC8Qhen7fqD6jVRCxBAl4IuKNaVoPFaPbx4nb43gz04Xb43q2q+2ds86fAZTM6fD+gqnfN9b7V1OZfafmBCYa+9gKRd715xs/5FLJZnnv0u7zw4x8Sjng9js0AAA7PSURBVMe59q6PsPWGW2e9KOxC2bZD5nSW0cFJ0iNTZEayZJI5JkazTKbzTGYKFM8xBHBO8sbdYo1QtizBF/Tg9b9x8PEF3Of+oMc9KPk90wcnr8/CLoyRmzgNQ2miQxaxqRgWFqfzfbw2to+TE4ex1a2xe/0Boo1NpVsz0cZGIg1NRBoaiSQaqEs0EI7HZ70q2zAuRKWHeu4A7sc9KD6oql8SkS8Ce1X1MREJAg8Bl+PW+Heq6pyTxJvkP7fk9w8zuW+Its9dgfcChmMOHj3C09/+BgOHD5Jo7+S6nR/l4iuvrtjIkkLeJjdRJDdZIDdZdGveWZtCrohd1Onau2O7p8LOjAuNVNWtiYlgWW6t7PUmCo/XwuOV6TOH1+99fg9ev4XX93qit/D6PXhmGQevqkykkqQG+kj295Ea6CN1qp/RgX7GBofoDHWzKXYliUAreSfLgHOMsfoU/jUx4i2txJtbiDW3EmtuIRSNmRq6UVHmIq9Vzh7Lceof9xLc3EDj3Zsv6D1UlSN79/DL732HZH8vDR1dvP2OO9ly3U14a2Ca5mI+T+pUP8m+XlL9vSSnb30Usm+sh+z1+WlsW0t3fDutuS68RS8atwi9o4XEtevxBC+sM9AwloJJ/jUg/dRx0k+doOmPthK8JDH/L5yDY9scLC0NOdxzlFAszpbrb+LS62+ied2GRSzx8pjKpEn2uYl9pO+km+j7ehkbGiyNknFFG5tp6Oikob2TxJp2Eu2d1Mdb4UCOiT0DaM4mcHE90es7CXTXmxq9sSKZ5F8DtOAw+M8voKq0/fn2BV8Gr6qc3P879j3+GEdf+C2ObdO8fiPdV13DxsuvpGXDRSs24dnFAqODp0gN9Jdq8X2kBtwkP5V5Y1F7j89Hoq2dho4u99be4d6v6ThjPQQ7kyfzi14mfjOAFh1CW5uI3tCFv2Pp5j8yjMVgkn+NyB4Z5fQ3f0f0pi7it61ftPedTI9x6FfP8Mqzuxl47TCoUpdooGvLZbRv2kzHpi00dq590/oBSyk7MU56eIix4UHGTg0wOniK0cEBRk/1kx4ePqMWH4rFaWjvdGvyazqmk32suXnOzlV7PE9mdy/jewbAdghf3kL0xi58zWYKaqM6mORfQ5K7DjH50jCtn92+JOODJ0ZT9Lz4Akf37aXv4H4mUu4oXcvjpaG9g8bOtdS3tbsjWpqaqIsnCEYiBOoi+EOhWZOtqlIs5CnmcuSnpshOjJMdz5AdzzAxOsrk2CgTo0nGkyOMJ0fIjJwmNzlxxnsE6uqIt7SRaGsn0d7h3q/pILGmg2Dk/GrozmSBzLN9jD/XhxbcpB+7aS3epvNb/9gwlptJ/jXEHs9z6qvP42sL0/yJbUs6I6eqkh4eov/wKwyf6GHk5HFG+k6SHh5CHWfW3xHLwuN156hRx0EdG9u25xyvKZZFOF5PJNE4PTzSHUXjjqSJt7YRikQXHk/BZvxX/aR396JTRULbmojdus7U9I2qVckrfI1l5on4qb9jA6nvv8r4L/uIXt+5ZJ8lIu5wxpZWZo4xchybiVSKzMgwk2NjZCfGyU2Mk89OYReK2MUC6tiI5XEPBh4P3kAQXyCILxggGIkSrIsQjESpq08QikSXdNipOsrk84OknzyOnc4T3JQg9u71+JdwTQPDWElM8l8lwm9vZeqVJGNP9BC4uL7iScyyPNMXM61kqkr2cIqxnxyjODiJvytKw85NBDbWL3fRDKOiamOu2BogIiQ+0I0V9pLcdchdgck4Q75/nNPfepmRb+9Hiw4NH3kLzZ9+q0n8Rk0yNf9VxFPno+HDmzj94MuMPX6M+vdetNxFWhHsdJ6x/+1h8vlBrJCX+Hs2Erl6jVl31qhpJvmvMsFLEkSubWf8uX7862KEtzUvd5GWjRZsMs/2kdl9ErWVyDs7iN3UVdZkeIax2pnkvwrFd2wg3zdO6tHDeJtCNdeJqapMvTTM2E97sEdzBC9tpP73N5hhm4YxgznvXYXEa9F4z2Yk5GXkoQPYE4X5f2mVyJ1IM/z1F0k+fAgr5KXpE5fRdO8Wk/gN4ywm+a9Snqifpnu3YGfyJL/7Clpc3R3AxdEsyUcOMvyvL1JMZUl8sJuWz1xO8CLTmWsYszHNPquYvytK4v3dpB49THLXIRp2vmVVLVQN4ORsMr84SeaZPkCJ3thF9IZOrIDZtQ1jLuYbssrVvb0VZ7LA2I+PkfK/SuKD3Ut6BXClqK1MPH+K9JPHcTIFQm9rJv7u9XgT57+2gWHUIpP8a0D0uk6crE3m6RNYAQ/x925csbNzzkdVyR5MMvbTHopDk/jXxYjfu4XA2thyF80wqopJ/jUidstaNGcz/ss+nLxN4v0XI57q6vLJ9Ywx9ngP+Z403qYQjfdsJnhpY9UeyAxjOZnkXyNEhPgdGxC/ReZnJ3EyeRru3owVWPlrx+b7x0k/eZzsK0msqI/6Oy+i7qq2qjt4GcZKYpJ/DRER4retxxMPMPqDIww/8BKN927GW78y28nzAxOknzpOdv8IEvQQe/d6Iu9sx/Kv/AOWYax0JvnXoMg71uCJ+Uk+fIjB+/eR+ODFhC9bOVcC506kyezuJXtgBAl4iN68lui7OrBCZnc1jMVivk01KrS5kdbPXs7II4dIfvcguatGie/YgBVcnl1CHbcjN/NsH/ljY1hhr5v039lupmMwjCVgkn8N8zaGaPmTbaSfPE7mF71MHRghdts66q5oq9hwUHuiwOTeQcb39GOncnjiAeLv2UjdlW1V0R9hGNXKJP8aJx6L+O0bCG1tYvRHRxn97yNM/GqA6I1dhLY2LkmnqhYdsgeTTLwwRPZgEhzFvyFOfMdGQlsaV92FaIaxEpnkbwDg74zS/MfbmPrdadJP9JB8+CBW1E/kHW2E3tqMtym0oCGV9nie3KujTB0YIXs4heZsrKiPyLvaqdveiq+tbhGjMQxjPgtK/iLSAOwC1gM9wF2qmjprm7cBXwdigA18SVV3LeRzjaUhIoS3NRPa2kT2cIqJX/eTfuoE6adO4GkIEtyUwN8VxdcSxtscnrVZRlVxMnkKw1MUh6fI92bIH09THJ4CwIr6CL+1mdCljQQuTphavmEskwUt4C4iXwGSqvplEflrIKGqnz9rm0sAVdVXRaQdeB7YrKqjc723WcB9ZSimsmQPJskeSpF7bfSMFcLE70ECFuL3gK04ORvNFWHGHHIS8hJYF8O/PkZgYxx/Z3RVTC9hGCtVpRZwvxO4ofT4O8Bu4Izkr6qHZzzuF5EhoBmYM/kbK4M3ESRyTTuRa9pR26E4kqU4NElhaBJnsojmbZycjViCBD1YQS+emB9vUwhvcxhPzG+SvWGsQAtN/q2qOgCgqgMi0jLXxiJyFeAHXjvH658EPgmwdu3aBRbNWGzisfC1hPG1hDGz4xtGdZs3+YvIU0DbLC/9zfl8kIisAR4CPqaqs04ur6oPAA+A2+xzPu9vGIZhlG/e5K+qt5zrNREZFJE1pVr/GmDoHNvFgB8Df6uqey64tIZhGMaiWOgg7seAj5Uefwz44dkbiIgf+B/gP1X10QV+nmEYhrEIFpr8vwzcKiKvAreWniMiV4jIN0vb3AVcD9wnIv9Xur1tgZ9rGIZhLMCChnouJTPU0zAM4/yVO9TTTIhuGIZRg0zyNwzDqEEm+RuGYdSgFdvmLyLDwPEFvEUTcHqRilMtai3mWosXTMy1YiExr1PVeVdnWrHJf6FEZG85nR6rSa3FXGvxgom5VlQiZtPsYxiGUYNM8jcMw6hBqzn5P7DcBVgGtRZzrcULJuZaseQxr9o2f8MwDOPcVnPN3zAMwziHqk7+InK7iBwSkSOllcTOfj0gIrtKr/9GRNZXvpSLq4yY/1JEDojISyLytIisW45yLqb5Yp6x3YdEREWk6keGlBOziNxV+l/vF5HvVbqMi62MfXutiPxcRPaV9u8dy1HOxSIiD4rIkIi8fI7XRUS+Vvp7vCQi2xe1AKpalTfAg7sozEbcBWJeBLactc2ngW+UHu8Edi13uSsQ841AuPT4U7UQc2m7KPAMsAe4YrnLXYH/czewD3fpVICW5S53BWJ+APhU6fEWoGe5y73AmK8HtgMvn+P1HcBPAQGuBn6zmJ9fzTX/q4AjqnpUVfPAI7jLSs50J+7ykgDfB24WkWpeU3DemFX156o6WXq6B+iscBkXWzn/Z4C/B74CZCtZuCVSTsyfAP5FVVMAqjrrWhpVpJyYFYiVHseB/gqWb9Gp6jNAco5N7sSdCl/VXQelvrRuyqKo5uTfAZyc8by39LNZt1HVIjAGNFakdEujnJhn+jhuzaGazRuziFwOdKnqjypZsCVUzv/5EuASEXlORPaIyO0VK93SKCfmvwPuEZFe4CfAZypTtGVzvt/387LQNXyX02w1+LOHLpWzTTUpOx4RuQe4Avi9JS3R0pszZhGxgH8C7qtUgSqgnP+zF7fp5wbcs7tnRWSrqo4ucdmWSjkx/wHwH6r6VRG5BnioFPOsy8KuAkuav6q55t8LdM143smbTwOntxERL+6p4lynWStdOTEjIrfgrrH8PlXNVahsS2W+mKPAVmC3iPTgto0+VuWdvuXu2z9U1YKqHgMO4R4MqlU5MX8c+C8AVf01EMSdA2e1Kuv7fqGqOfn/FugWkQ2lpSJ34i4rOdPMZSY/BPxMSz0pVWremEtNIP+Gm/irvR0Y5olZVcdUtUlV16vqetx+jvepajWvBFTOvv0D3M59RKQJtxnoaEVLubjKifkEcDOAiGzGTf7DFS1lZT0GfLQ06udqYExVBxbrzau22UdViyLyZ8ATuCMFHlTV/SLyRWCvqj4GfAv31PAIbo1/5/KVeOHKjPkfgAjwaKlv+4Sqvm/ZCr1AZca8qpQZ8xPAbSJyALCBv1LVkeUr9cKUGfPngH8Xkb/Abf64r5orcyLyMG6zXVOpH+MLgA9AVb+B26+xAzgCTAJ/uKifX8V/O8MwDOMCVXOzj2EYhnGBTPI3DMOoQSb5G4Zh1CCT/A3DMGqQSf6GYRg1yCR/wzCMGmSSv2EYRg0yyd8wDKMG/T8dRnX+BFSnQAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "for Q in [0.2,0.3,0.4, 0.5, 0.6, 0.7, 1]:\n", " gamma = omega0/Q\n", " omegav = omega0 * sqrt(1 - 1/(4*Q**2))\n", " phi = 0\n", " A = 1.0\n", " plt.plot(t,np.real(z(gamma, A, omegav, phi)) ,label=Q)\n", "\n", "plt.legend()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Trying a numerical approach." ] }, { "cell_type": "code", "execution_count": 57, "metadata": {}, "outputs": [], "source": [ "from scipy.integrate import odeint" ] }, { "cell_type": "code", "execution_count": 99, "metadata": {}, "outputs": [], "source": [ "def sho(y, t, gamma, w2):\n", " x, xdot = y\n", " dydt = [xdot, -gamma*xdot - w2*x]\n", " return dydt" ] }, { "cell_type": "code", "execution_count": 130, "metadata": {}, "outputs": [], "source": [ "gamma = 6\n", "w2 = 9\n", "y0 = [np.pi/2, 0.0]\n", "t = np.linspace(0,10,100)\n", "sol = odeint(sho, y0, t, args=(gamma,w2))" ] }, { "cell_type": "code", "execution_count": 131, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 131, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD8CAYAAACMwORRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAG79JREFUeJzt3XuUnHWd5/H3p28J3UlIp7sDIZfuBMIluCp0i63sjCi6BmYXZmcdJeiMesCss+I4o2dHnNmjLp45M6O7M87shmEiy6IOwmEQNTooOwMoXojSIcgl4RI6t86FdC6QkJB0Ov3dP6o6FJ1OutJdlafqeT6vc3J46qlfqr4Vkk/9+vf8nt9PEYGZmaVLTdIFmJlZ6TnczcxSyOFuZpZCDnczsxRyuJuZpZDD3cwshRzuZmYp5HA3M0shh7uZWQrVJfXGra2t0dHRkdTbm5lVpVWrVu2MiLax2iUW7h0dHfT09CT19mZmVUnSxmLaeVjGzCyFxgx3SbdJ2iHpqRO0uUzS45KelvST0pZoZmYnq5ie++3A4uM9KWk6cDNwVURcCPxuaUozM7PxGjPcI+JhYPcJmlwL3BsRm/Ltd5SoNjMzG6dSjLmfCzRL+rGkVZJ+vwSvaWZmE1CK2TJ1QCdwOXAa8IiklRHx3MiGkpYCSwHmzZtXgrc2M7PRlKLn3gf8KCL2R8RO4GHgTaM1jIjlEdEVEV1tbWNO0xzVqo17WPbQOlZt3DP+is3MUq4UPffvAf9bUh3QALwV+JsSvO4xVm3cw7VfW8nA4BCT6mu44/puOtuby/FWZmZVbcxwl3QncBnQKqkP+AJQDxARt0TEWkk/Ap4AhoBbI+K40yYnYmXvLgYGhwhgYHCIlb27HO5mZqMYM9wjYkkRbb4CfKUkFZ1A94IWJtXXcPDw0NHHZmZ2rKq6Q7WzvZk7ru/mHee2MRQwbXJiqyeYmVW0qgp3yAX8X7//TTTU1vCPK4taYsHMLHOqLtwBWqZM4rfeOItvP7aF/YcGky7HzKziVGW4A3you51XDg3y3ce3JF2KmVnFqdpwv3jedDpaGvnqvzzHqg0nWh3BzCx7qjbcH9v0ElteepX+VwZYcusvfVOTmVmBqg33lb27ODIUABzOz3k3M7Ocqg337gUtNNTlyq+pkee8m5kVqNpwH57zPnv6ZM6dOcV3qpqZFajacIdcwL/r/DPo2/MqEZF0OWZmFaOqwx3g/FlT2XdokC0vvZp0KWZmFaP6w/3MaQA8s21fwpWYmVWOqg/3886cCsDabXsTrsTMrHJUfbhPmVRHe0sjz2x3z93MbFjVhzvA+WdOZe1299zNzIalJNynsWHnfl4dOJJ0KWZmFWHMcJd0m6Qdkk64u5Kkt0g6Iul9pSuvOBfMmspQwPM7PDRjZgbF9dxvBxafqIGkWuCvgPtLUNNJ84wZM7PXGzPcI+JhYKxlFz8JfBvYUYqiTta8GY00NtSyxjNmzMyAEoy5S5oN/EfglomXMz41NeK8M6fyjC+qmpkBpbmg+lXgsxEx5tVMSUsl9Ujq6e/vL8Fbv+b8M6fxzPZ9XobAzIzShHsXcJekDcD7gJsl/fZoDSNieUR0RURXW1tbCd76NRfMmspLBw7z5R8967XdzSzzJhzuETE/IjoiogO4B/gvEfHdCVd2kmoQALc8/AIfvHWlA97MMq2YqZB3Ao8A50nqk3SdpI9L+nj5yyvetr25hcMivHmHmVndWA0iYkmxLxYRH5lQNRPwrvNmsuyhFxBQX1fjzTvMLNPGDPdq0dkxg7PbmjgyFPzP97/Zm3eYWaalYvmBYW+YfTqDQ+FgN7PMS1W4d7Q0seWlVzl42GvMmFm2pSrc57c2EQGbdx9IuhQzs0SlKtw7WpsAWL9zf8KVmJklK1XhPr8lF+4bdjnczSzbUhXupzfWM6OpgfU7PSxjZtmWqnAH6GhpZP3OV5Iuw8wsUekL99YmNrjnbmYZl7pwn9/SxPa9B73lnpllWvrCvc0XVc3MUhfuHcMzZjwd0swyLH3hnp/r3utwN7MMS124T5lUR9vUSe65m1mmpS7cIXdR1WPuZpZl6Qz31ibfyGRmmVbMTky3Sdoh6anjPP9BSU/kf/1C0ptKX+bJ6WhtYucrh9h38HDSpZiZJaKYnvvtwOITPL8eeEdEvBH4ErC8BHVNyPzWRgC+4s2yzSyjxgz3iHgY2H2C538REcMJuhKYU6Laxu1A/gamb67c6M2yzSyTSj3mfh3wwxK/5kkbXs898GbZZpZNJQt3Se8kF+6fPUGbpZJ6JPX09/eX6q2P8W8Xth099mbZZpZFJQl3SW8EbgWujojjdpMjYnlEdEVEV1tb2/GaTVhnezOLZk1l9vTJ3HF9t/dUNbPMmXC4S5oH3Av8XkQ8N/GSSuPfzJ7OwBFvlm1m2VQ3VgNJdwKXAa2S+oAvAPUAEXEL8HmgBbhZEsBgRHSVq+Bitbc20t9ziP2HBmmaNObHNDNLlTFTLyKWjPH89cD1JauoRIYXENu46wCLzpqWcDVmZqdWKu9QBZg3IzfXfaOXITCzDEptuLe35MJ9wy4vQ2Bm2ZPacJ86uZ7WKQ3uuZtZJqU23AHavTqkmWVUysO9kY0eljGzDEp1uHe0NLHt5YMcPOzNss0sW1Id7sMXVTftdu/dzLIl1eFeONfdzCxLMhLuvqhqZtmS6nA/vbGe6Y31njFjZpmT6nAHaJ/hGTNmlj3pD3fPdTezDEp9uHe0NLJlz6sMDA4lXYqZ2SmT+nBvb2liKKBvj4dmzCw7Uh/uHa25ue7LHlrnjbLNLDNSH+4vv3oYgHsf28IHb13pgDezTBgz3CXdJmmHpKeO87wk/Z2kdZKekHRx6cscv7Xb9gIQwOHBIVb2HneLVzOz1Cim5347sPgEz18BLMz/Wgr8/cTLKp3uBa3kdv+D+roauhe0JFuQmdkpMGa4R8TDwO4TNLka+EbkrASmS5pVqgInqrO9mUvPbmXa5DruuL7bG2abWSaUYsx9NrC54HFf/lzF6GxvZt+hQd4w23upmlk2lCLcNcq5GLWhtFRSj6Se/v7+Erx1cRa0NREBm3ynqpllRCnCvQ+YW/B4DrB1tIYRsTwiuiKiq62trQRvXZz5rbkFxHp3+k5VM8uGUoT7CuD387NmuoGXI2JbCV63ZDry4b7e4W5mGVE3VgNJdwKXAa2S+oAvAPUAEXELcB9wJbAOOAB8tFzFjte0/GbZGxzuZpYRY4Z7RCwZ4/kAPlGyispkfmuTh2XMLDNSf4fqsPmtTR6WMbPMyFC4T6F/3yH2HTycdClmZmWXoXDPLSDmjTvMLAsyFO5TAE+HNLNsyEy4t7c0IsH6foe7maVfZsJ9cn0tZ51+Gut3vpJ0KWZmZZeZcAfPmDGz7MhkuOem5puZpVfmwn3vwUF27x9IuhQzs7LKVri3eY0ZM8uGTIX7gvwCYrf+bL33UjWzVMtUuL+49yAA9z+13Ztlm1mqZSrcH92QC3Nvlm1maZepcO9e0EKNN8s2swzIVLh3tjfzny6ejYDbP/oWb5ZtZqmVqXAH+I1zZxLA9MaGpEsxMyubzIX7wpm5BcSef9HLEJhZehUV7pIWS3pW0jpJN47y/DxJD0laLekJSVeWvtTSmN/aRI3g+R0OdzNLrzHDXVItsAy4AlgELJG0aESz/wbcHREXAdcAN5e60FKZXF/LvBmNvOBwN7MUK6bnfgmwLiJ6I2IAuAu4ekSbAKblj08HtpauxNI7Z+ZUnt+xL+kyzMzKpphwnw1sLnjclz9X6IvAhyT1AfcBnxzthSQtldQjqae/v38c5ZbGwjOmsH7nfg4fGUqsBjOzciom3DXKuZHLKi4Bbo+IOcCVwDclHfPaEbE8Iroioqutre3kqy2RhTOncPhIeMs9M0utYsK9D5hb8HgOxw67XAfcDRARjwCTgdZSFFgO5+RnzKzzuLuZpVQx4f4osFDSfEkN5C6YrhjRZhNwOYCkC8iFe3LjLmM4u2043D3ubmbpNGa4R8QgcANwP7CW3KyYpyXdJOmqfLPPAB+T9GvgTuAjUcE7YjRNqmP29NM8HdLMUquumEYRcR+5C6WF5z5fcLwGuLS0pZXXwjOm+EYmM0utzN2hOuyctim80P8KR4Yq9gcMM7Nxy2y4LzxjCocGh9iy59WkSzEzK7nMhvs5M6cC8NUHnvOmHWaWOpkN9/2HBgH4zmNbvCuTmaVOZsP9yS0vA96VyczSKbPh7l2ZzCzNMhvune3NXP3m2dQIvv7RS7wrk5mlSmbDHeBd589kKGDK5KKm+5uZVY1Mh/uis3KrFK/ZujfhSszMSivT4d7R0sRp9bWs3eY1ZswsXTId7rU14vxZU1mz7eWkSzEzK6lMhzvABbOmsWbrXip4nTMzs5OW+XBfNGsaew8OsuUlL0NgZunhcPdFVTNLocyH+/lnTkWCNdsc7maWHpkP98aGOua3NrnnbmapUlS4S1os6VlJ6yTdeJw275e0RtLTkr5V2jLLa9Gsae65m1mqjBnukmqBZcAVwCJgiaRFI9osBD4HXBoRFwJ/VIZay+aCWdPo2/MqL796OOlSzMxKopie+yXAuojojYgB4C7g6hFtPgYsi4g9ABGxo7RlltfwRdVn3Hs3s5QoJtxnA5sLHvflzxU6FzhX0s8lrZS0eLQXkrRUUo+knv7+/vFVXAYXzsqF+z883Ot13c0sFYoJd41ybuQdP3XAQuAyYAlwq6Tpx/ymiOUR0RURXW1tbSdba9ls3n0AgAef2eGNO8wsFYoJ9z5gbsHjOcDWUdp8LyIOR8R64FlyYV8VVq7fffTYG3eYWRoUE+6PAgslzZfUAFwDrBjR5rvAOwEktZIbpuktZaHl1L2ghbr8zh31td64w8yq35jhHhGDwA3A/cBa4O6IeFrSTZKuyje7H9glaQ3wEPBfI6Jqur+d7c187srzAfjs4vO9cYeZVb2idqmIiPuA+0ac+3zBcQCfzv+qSu/rnMuXfrCWV/IbZ5uZVbPM36E67PTT6lk4cwqrN7+UdClmZhPmcC9w0bzprN60x8v/mlnVc7gXuGheM3sOHGbjrgNJl2JmNiEO9wIXzctNzV+92fPczay6OdwLLJw5laaGWlZv8ri7mVU3h3uB2hrxprnTHe5mVvUc7iNcNG86a7ft5dWBI0mXYmY2bg73ES6a28zgUPDU1peTLsXMbNwc7iO8OX9RddmD67yAmJlVLYf7CBt3HUDAj5/r9wqRZla1HO4jrOzddXQ9Y68QaWbVyuE+QveCFuprcytE1nqFSDOrUg73ETrbm7n52k4Arr1krleINLOq5HAfxXsuPIMFbU1ehsDMqpbD/TguPbuVX63fzeEjQ0mXYmZ20hzux3HpOS3sHzjCr70EsJlVoaLCXdJiSc9KWifpxhO0e5+kkNRVuhKT0b2gBQl+vs6zZcys+owZ7pJqgWXAFcAiYImkRaO0mwr8IfDLUheZhOmNDbzhrNP5+Qs7ky7FzOykFdNzvwRYFxG9ETEA3AVcPUq7LwFfBg6WsL5Evf2cFlZv2sOBAW+9Z2bVpZhwnw1sLnjclz93lKSLgLkR8YMS1pa4S89u5fCR4NENvkvVzKpLMeGuUc4d3YdOUg3wN8BnxnwhaamkHkk9/f39xVeZkLd0zKC2BpY9+LyXITCzqlJMuPcBcwsezwG2FjyeCrwB+LGkDUA3sGK0i6oRsTwiuiKiq62tbfxVnyJrtu0lAn61YY/XmTGzqlJMuD8KLJQ0X1IDcA2wYvjJiHg5IlojoiMiOoCVwFUR0VOWik+hlb27GN4re8DrzJhZFRkz3CNiELgBuB9YC9wdEU9LuknSVeUuMEndC1poqMv9EdVIXmfGzKqGImLsVmXQ1dUVPT2V37lftXEP//mbPcycOon7PvWbSZdjZhknaVVEjHkvke9QHUNnezPXXjKPZ7bvY9crh5Iux8ysKA73Ivy7C89kKOBf176YdClmZkVxuBfhwrOmMXv6adz/tMPdzKqDw70IknjvhWfys+d38soh361qZpXP4V6k9154BgNHhvjcvU94vruZVTyHe5FqanI36n7/19t8Q5OZVTyHe5F+tX730WNvnG1mlc7hXqTuBS001OZvaKrxDU1mVtkc7kXqbG/mzo+9leameha0TfHG2WZW0RzuJ6GzYwYf/82zeXb7Ptbt2Jd0OWZmx+VwP0m/c/Ec6mrE3T19SZdiZnZcDveT1DZ1EpdfMJNvr+pjYHAo6XLMzEblcB+HD7xlLrv2D/An9/zaUyLNrCI53MdhyqQ6AL77+FbPeTeziuRwH4dHN+w5uvegN/Ews0rkcB+Hwk08hOe8m1nlKSrcJS2W9KykdZJuHOX5T0taI+kJSQ9Iai99qZWjs72Zb32sm4vnTQeCM0+fnHRJZmavM2a4S6oFlgFXAIuAJZIWjWi2GuiKiDcC9wBfLnWhlaazvZn/de3FSGL5T15Iuhwzs9cppud+CbAuInojYgC4C7i6sEFEPBQRB/IPVwJzSltmZZo9/TR+5+LZ3PXoZvr3eZcmM6scxYT7bGBzweO+/LnjuQ744USKqiZ/cNk5DAwO8YlvrfKsGTOrGMWEu0Y5N+qu2pI+BHQBXznO80sl9Ujq6e/vL77KCrZ7/wAS/Gr9Hq79mqdFmlllKCbc+4C5BY/nAFtHNpL0buDPgKsiYtQxiohYHhFdEdHV1tY2nnorTuE0yEOeFmlmFaKYcH8UWChpvqQG4BpgRWEDSRcB/0Au2HeUvszKNTwtcvjHm+bGhkTrMTODIsI9IgaBG4D7gbXA3RHxtKSbJF2Vb/YVYArwT5Iel7TiOC+XOp3tzdxxfTefunwhLU0N3PHLjRwZGnXUyszslFFEMkHU1dUVPT09ibx3ufzgia3c8K3VvPuCmfzBZed4zXczKzlJqyKia6x2vkO1hGZNm0yN4F/X7vDFVTNLlMO9hFYW7LN6aHCIX7ywM8FqzCzLHO4lNPLi6qZd+xOtx8yyqy7pAtJk+OLqyt5d9GzYzT+t2oIkPvCWeR5/N7NTyj33Eutsb+YT7zyHj/3GAgTc3dPn8XczO+Uc7mWyevNLKD8+c2hwiAfWvphsQWaWKR6WKZPh8feBwSGGAv75iW28df4Mntq6l+4FLR6mMbOy8jz3Mlq1cQ8re3fRWF/Ln9+3liMRCGioq+GO67sd8GZ20oqd5+6eexl1tjcfDfAn+l7iO49vJXhtaz6Hu5mVi8fcT5EPva2DhtrcH/dQwEsHBlj20DpfaDWzsnDP/RTpbG/mzqXd/PCpbXx39Ra+9tP1CJhU7yEaMys9h/spNDxM01hfy989uI4ADh4e4raf9bKy93RfaDWzknG4J+Ad581k+U97X5tJ8+R27ntyOw11NXzhP1zIngMDDnozmxDPlknI8Eya3v5X+PZjW46er8nPjfeMGjMbjWfLVLjhIZpVG/fwz09uY2BwiIjcxVbIDdfc/ON1XDxvOt0LWh3yZnZS3HOvAMO9+ObGBv7795/OBX3B87U14iNvb2d6YwNvP7sVyG3v56Ebs+wptudeVLhLWgz8LVAL3BoRfzni+UnAN4BOYBfwgYjYcKLXdLiPbjjoN+3az909fcfsRC6gpkZEBPW1rx+jB4e+WdqVbFhGUi2wDHgPuc2yH5W0IiLWFDS7DtgTEedIugb4K+AD4ys92wqHa773660cHhxCEkeGggACjm7jd2hwiD/9zpNAbqxeEkMRNNTW8NnF5/Hq4aFjQv94x/4yMEuXMXvukt4GfDEi3pt//DmAiPiLgjb359s8IqkO2A60xQle3D33sRUO19z0g6c5PDhEbf6K62A+4IvZrlXkvhRq8g8i8hduJYaGgrpa8WdXXkBDXS3PbN9L57xm6mrFk1te5pL5M6hTDY9t2kP32TOoranhV+t30b3g9cNDp/J4+MsvifdOQ33VVGul1zfRWsfTqSrlBdXZwOaCx33AW4/XJiIGJb0MtADeimgCCpcvOO/Mqcf85RgZ+oW9+sLQHz4cKnhwJMilPHD4SPDF77/2g9g3Htl49PiWn/S+9kIPFFb33IQ/33gIaJnSwK5XBo5+ruEvr0o4FtA6tYGd+yqzvmqqtdLrm0itUP4bGIsJd41ybrSh4LHaIGkpsBRg3rx5Rby1DSsM+uHHUFzoI3HkyLHHw18GdbU1vGNhGw888yJD8fq/iJVm+B9SjDhXKcfD35lJ15GGWiu9vonWerjMa0wVE+59wNyCx3OArcdp05cfljkd2D2iDRGxHFgOuWGZ8RRsr3cyoT/W8U/X9Y/5hZD0cX1dDX/8nvOK+vJKqr5PV3B91VRrpdc30Vrr62qO/tsrh2LC/VFgoaT5wBbgGuDaEW1WAB8GHgHeBzx4ovF2K7/jhf6Jjoe3CKyUMczjHXe2N4/ry8v1VV+tlV7fRGst50SGYqdCXgl8ldxUyNsi4s8l3QT0RMQKSZOBbwIXkeuxXxMRvcd/RV9QNTMbj5LeoRoR9wH3jTj3+YLjg8DvnmyRZmZWHl7P3cwshRzuZmYp5HA3M0shh7uZWQo53M3MUiixJX8l9QMbx/nbW8ne0gb+zNngz5wNE/nM7RHRNlajxMJ9IiT1FDPPM038mbPBnzkbTsVn9rCMmVkKOdzNzFKoWsN9edIFJMCfORv8mbOh7J+5KsfczczsxKq1525mZidQdeEuabGkZyWtk3Rj0vWUm6S5kh6StFbS05I+lXRNp4KkWkmrJf0g6VpOFUnTJd0j6Zn8/++3JV1TOUn64/zf6ack3ZlfXTZ1JN0maYekpwrOzZD0L5Kez/+35Gv/VlW4F2zWfQWwCFgiaVGyVZXdIPCZiLgA6AY+kYHPDPApYG3SRZxifwv8KCLOB95Eij+/pNnAHwJdEfEGcsuJX5NsVWVzO7B4xLkbgQciYiG5DSxL3lGtqnAHLgHWRURvRAwAdwFXJ1xTWUXEtoh4LH+8j9w/+NnJVlVekuYAvwXcmnQtp4qkacBvAv8HICIGIuKlZKsquzrgtPzubY0cu8NbKkTEwxy7M93VwNfzx18HfrvU71tt4T7aZt2pDrpCkjrIbYjyy2QrKbuvAn9Cfk/vjFgA9AP/Nz8cdaukpqSLKpeI2AL8D2ATsA14OSL+X7JVnVJnRMQ2yHXggJmlfoNqC/eiNuJOI0lTgG8DfxQRe5Oup1wk/XtgR0SsSrqWU6wOuBj4+4i4CNhPGX5UrxT5MeargfnAWUCTpA8lW1W6VFu4F7NZd+pIqicX7HdExL1J11NmlwJXSdpAbtjtXZL+MdmSTok+oC8ihn8qu4dc2KfVu4H1EdEfEYeBe4G3J1zTqfSipFkA+f/uKPUbVFu4H92sW1IDuQswKxKuqawkidw47NqI+Ouk6ym3iPhcRMyJiA5y/38fjIjU9+giYjuwWdJ5+VOXA2sSLKncNgHdkhrzf8cvJ8UXkEexAvhw/vjDwPdK/QZF7aFaKSJiUNINwP28tln30wmXVW6XAr8HPCnp8fy5P83va2vp8kngjnzHpRf4aML1lE1E/FLSPcBj5GaErSald6pKuhO4DGiV1Ad8AfhL4G5J15H7oiv5HtS+Q9XMLIWqbVjGzMyK4HA3M0shh7uZWQo53M3MUsjhbmaWQg53M7MUcribmaWQw93MLIX+PyBbQCWFq86yAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.plot(t, sol[:,0], \".-\", label='x(t)')" ] }, { "cell_type": "code", "execution_count": 63, "metadata": {}, "outputs": [], "source": [ "odeint?" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernel_info": { "name": "python3" }, "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.6.5" }, "nteract": { "version": "0.11.6" } }, "nbformat": 4, "nbformat_minor": 4 }