{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# The puzzle about \"tangential to the slow direction\" - sheep and rabbits model" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [], "source": [ "# import python libraries\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "from typing import Tuple \n", "from scipy.integrate import solve_ivp # used in numerical integration\n", "#from scipy import optimize" ] }, { "cell_type": "code", "execution_count": 18, "metadata": {}, "outputs": [], "source": [ "def simulate(n_timesteps: int, state0: list=[0.0, 1.0], h: float=0.02) -> Tuple[np.ndarray, np.ndarray]:\n", "\n", " def diff(t, state):\n", " x, y = state \n", " return y*(3-x-2*y), y*(2-x-y) # dx/dt, dy/dt\n", " #return y-x, x*x\n", "\n", " t = np.arange(0, n_timesteps * h, h)\n", " \n", " sol = solve_ivp(diff, \n", " y0=state0, \n", " t_span=(0.0, n_timesteps*h),\n", " dense_output=True)\n", " \n", " return sol.sol(t).T, t" ] }, { "cell_type": "code", "execution_count": 49, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 49, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAD4CAYAAAAAczaOAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAAA8AUlEQVR4nO3dd3xUVf7/8ddJryQhDQiE0JuCIFWQKlWQRdSl2YEVAcXuuiqKouIiiihVUUERVwGVDlJEeijSWwgEkgDpvc+c3x8T95cvG8gAk9xk5vN8PHiQmXtv5nO55J0zZ849R2mtEUII4RicjC5ACCFExZHQF0IIByKhL4QQDkRCXwghHIiEvhBCOBAXowsoTVBQkI6IiDC6DCGEqDL279+fpLUOLmu/Shn6ERER7Nu3z+gyhBCiylBKxVizn3TvCCGEA5HQF0IIByKhL4QQDkRCXwghHIiEvhBCOJAyQ18pVUcptUUpdUIpdUwp9Wwp+yil1KdKqSil1GGlVJsS2/oppU4Vb3vV1icghBDCeta09IuAF7TWzYCOwHilVPOr9ukPNCr+MxaYA6CUcgY+L97eHBheyrFCCCEqSJmhr7W+pLU+UPx1JnACCLtqt8HAIm2xG/BXStUE2gNRWutorXUBsLR4XyGEuGXabCZjzRrSlq8wupQq44b69JVSEUBrYM9Vm8KAiyUexxY/d63nS/veY5VS+5RS+xITE2+kLCGEA8revZvzDz5E3PMvkL58ObI2iHWsviNXKeUDLAMmaa0zrt5cyiH6Os//75NazwfmA7Rt21aunhCiVHknT5Iw/SOyt2/HpVZNan7wPn6DBqFUaXEjrmZV6CulXLEE/nda6+Wl7BIL1CnxuDYQD7hd43khhLghBbFxJH46k4yVq3CqVo2Ql18mYOQInNzdjS6tSikz9JXl1+eXwAmt9Yxr7PYrMEEptRToAKRrrS8ppRKBRkqpekAcMAwYYZvShRCOoCg1leS5c0ld8j04ORE4ejSBY0bjXK2a0aVVSda09DsDDwNHlFJ/Fj/3GhAOoLWeC6wBBgBRQA7wePG2IqXUBGA94Aws1Fofs+UJCCHskzk3l5RvFpH8xReYc3Lwu38IwRMm4FqjhtGlVWllhr7Wejul982X3EcD46+xbQ2WXwpCCFEmXVRE2vLlJM36jKLERHx69iTkuUm4N2pkdGl2oVJOrSyEcDxaa7I2bSJhxscUREfj2bo1YZ98jNeddxpdml2R0BdCGC5n/34Spn9E7sGDuNWvT+3PZuHTq5eMyCkHEvpCCMPkR0WRMONjsjZvxiUkhBrvTMF/yBCUi0RTeZF/WSFEhSu8fJnEWbNIX/EzTl5eBD/3HNUfeRgnT0+jS7N7EvpCiApjysggecECUhYtBrOZ6g8/TOBT/8AlIMDo0hyGhL4QotyZ8/NJ/W4JSfPmYc7IwO++QQRNfAa32qXOyiLKkYS+EKLcaJOJ9JUrSfz0U4riL+F9992EvPA8Hk2bGl2aw5LQF0LYnNaarN9/J/HjT8g/dQqP226j1nvv4d2xo9GlOTwJfSGETeXs30/CjI/J3b8f1/Bwwj6egW/fvignWaivMpDQF0LYRN6pUyTO+Jis33/HJTiYGm9Nxn/oUJSrq9GliRIk9IUQt6Tg4kUSP51FxqpVOPn6EvzC81QfNUqGX1ZSEvpCiJtSmJBgmf3yPz+iXFwIHDOGwCefwNnPz+jSxHVI6AshbogpI4PkLxeSsmgRurAQ/wcfIGjcOFxDQowuTVhBQl8IYRVzbi6p331H0oIvMKenU23gQIKfmYhbeLjRpYkbIKEvhLguXVhI2vIVJH3+OUUJCfh060bwc5NkrH0VZVehnzR3Ht53dcKzZUujSxGiytNmM5nr1pEwcyaFMRfwbNOGsBkf4dW2rdGliVtgzXKJC4GBQILW+rZStr8EjCzx/ZoBwVrrFKXUeSATMAFFWuty+99iysgg9bvvSPzkE/z6dCP4+Um4RkhLRIgbpbUme/t2Ej7+mPzjJ3Bv3Jjac+fg062bTHVsB5Rl0avr7KBUVyALWFRa6F+17yDgOa11z+LH54G2WuukGymqbdu2et++fTdyCACmrGyS580lZeEXoDRBgztR/bVPcfL2veHvJYQjyjl4kMQZH5MTGYlr7doEP/sM1QYMQDk7G12aKINSar81Desyb5HTWm8DUqx83eHA91bua3POPt6EvPAC9b+bg0+DaiQu20109w5kfP0hZf1yE8KR5Z0+zcWnxxMzfAT5584R+sbrNFizGr9BgyTw7UyZLX0ApVQEsOp6LX2llBcQCzTUWqcUP3cOSAU0ME9rPf86x48FxgKEh4ffGRMTcwOnUQqtyf7pM658PIf8FI1XuBehb7+HR6e+t/Z9hbAjBbFxJM2aRfqvv+Lk7U3g6NGWee29vIwuTdwga1v6tgz9vwOjtNaDSjxXS2sdr5QKATYCE4vfOVzXzXbvlEbn5ZD20fMk/mcrpgLw71iP4Hdn4xJWzybfX4iqqCg5maS580hduhTl5ET1h0cROHo0zv7+RpcmbpK1oW/L0TvDuKprR2sdX/x3glJqBdAeKDP0bUl5eBHwr7lUezyaxDfHkbrjHBn9+xP0QA+qv/wxysOjIssRwlCmzExSvvqK5K+/Qefn4z90KEHjn8Y1NNTo0kQFscm0d0opP6Ab8EuJ57yVUr5/fQ30AY7a4vVuhnOt+tT4Yj31v5qBZ5gXCUu2Et39TrJ+mG1USUJUGHNeHslffc3Z3n1Imj0H3+7dqL9qJTWnvC2B72CsGb3zPdAdCAKuAJMBVwCt9dzifR4D+mmth5U4rj6wovihC7BEaz3VmqJs2b1TGm02k/XtNBI+X0RBOng39CP0nX/j3vrucntNIYxgLigg7aefSJ47j6KEBLy7dCH4uUl4tmhhdGnCxmzap1/Ryjv0/6JzMkh5fyJJP+/BXAQB3ZsSPGU2zsG1yv21hShPuqiI9F9+Ienz2RTGx+PZ9k6Cn3kG7/btjS5NlBMJ/RtQdP4YiW9MIC3yEs7uEDxqIP7Pvodyc6uwGoSwBW0ykbFmDYmffUZhzAU8Wra0hH3nu+TGKjtns3H6jsAlogU1F2+h3pwpuAe7cfnL1ZzreSfZKxcZXZoQVtFmMxnrNxA9eDDxL72Mk6cXtWfPJuKHpfh06SyBL/5LWvpX0SYTmQsmk/DFTxRmKXxbBBHy7ie4NbvTkHqEuB6tNVlbt5L46SzyT5zArX59gp+ZiG+fPrI8oYOR7p1bZM5IIuWdp0laexjMUL3vHQRO/hxn/0BD6xICLGGfs2sXCTNnknfoMK7h4QSPf5pqAwfKHbQOSkLfRgpPRZL45iTSD6Xg7AUhTz6E31Nvyg+WMEzOvn0kzvyUnMhIXGrWJOjpcfj/7W+yFq2Dk9C3sdzVC7ny7xnkXjbhUcOd0NffwOueoUaXJRxI7uHDJM78lOwdO3AODiLoH0/h/9CDOMmAA4GEfrnQhQVkfP4qCYvWUJSjqNa6JiHvfoZrg+ZGlybsWN6JEyR+OousLVtwDgggcMwYAoYPk4XHxf8hoV+OzEmxJL39NCmbTwMQOKgjga9/ipNPNYMrE/Yk/+xZEmd9Rua6dTj5+hL45BMEjHoYZx9vo0sTlZCEfgUoOPQ7CZNfIvNkJi4+ipBxj1LtiZdleJy4JfnnzpE8dy7pK1fh5OFBwKOPEPjYYzj7+RldmqjEJPQritbkrPicyzNmk5+k8azjRejkqXh26Wd0ZaKKyY+OJmnOXDJWr0a5uREwYgSBo5/EpXp1o0sTVYCEfgXT+TmkzXiexB+2YMpT+HWIIGTqbFxq1ze6NFHJ5Z89S9LsOWSsWYPy8CBg2DACn3wCl6Ago0sTVYiEvkFMl6NJenMcKdtjcHKGwKE9qP7KRzh5yqIU4v/KP3OGpDlzyFi7DuXpSfURw6n++OO4BMq9IOLGSegbrGDvGq68/TpZZ3NxreZEyMSx+I56Rvr7BXmnTpM0Zw6Z69fj5OlJwMiRVH/8MenGEbdEQr8y0Jrs7z/kymdfk58CnnV9CX3rPTw73WN0ZcIAeSdPkvT5bDI3bsTJ25uAUaOo/tijuAQEGF2asAMS+pWIzs0kbfqzJC7bgSnPCb8O9Ql+9zNc68iSjY4g7/hxEmfPJuu3TTj5+BDw8CgCH31UliYUNiWhXwmZ4k6S/NZ4UnbEgZMicEgPAv85XRahtlO5R4+RNHs2WZs34+TrS/VHHqH6Iw/L0EtRLmw2tbJSaqFSKkEpVepSh0qp7kqpdKXUn8V/3iyxrZ9S6pRSKkop9eqNnYL9cQ5rSsiCTdT/6kN86rmT9OMWznZtT/rXM9Fms9HlCRvJOXiQi0+N4/wDD5Czbx9BEyfQcNNvBE+cIIEvDGfNcoldgSxgkdb6tlK2dwde1FoPvOp5Z+A00BuIBSKB4Vrr42UVZa8t/f/DbCZn6Qdcmf0NeUlOeNS29Pd7dZH+/qpIa032jp0kz59Pzt69OPv5Uf2xRwkYNQpnX1+jy7NruUW5pOWlUdOnptGlGMpmLX2t9TYg5SZqaA9Eaa2jtdYFwFJg8E18H/vk5ITXiNeI2LCHWqPaUZScRszoicSOupeCmGijqxNW0mYzGRs2cP6BB7k4ejQF588T8uorNNy8iaBx4yTwy9m22G0M+WUIL257kcrYVV0Z2WqVhU5KqUNKqbVKqb9WXA4DLpbYJ7b4uVIppcYqpfYppfYlJibaqKzKT3n54/f6IhqsXEFQtxpkHYwiesC9JPxrPKbMTKPLE9egCwtJ+/lnogfdR9wzz2LKzKTGlLdp8NtGAh97DCdvmR+nPF3OvszzW59n/KbxuDu7M6nNJBkObSUXG3yPA0BdrXWWUmoA8DPQCCjtClzzV7HWej4wHyzdOzaoq0pxqt2C4Hlb8N+9jMQP3iZ52WbS1nYi+Okx+D8+QebvryTMeXmkLVtGypcLKYyPx71xY2p9NJ1qffuiXGzx4ySup8hcxNKTS5l1cBYmbeLZNs/yaPNHcXWWtQSsdcstfa11htY6q/jrNYCrUioIS8u+ToldawPxt/p69s6141BqLdtPxORhuPkUcHn6XM7d04nsLRuMLs2hmbKySFqwgKhe93DlnXdxCQmh9pzZ1PvlZ/zuvVcCvwIcSTzCiNUjmBY5jTahbVgxeAWjbx8tgX+Dbvl/qlKqBnBFa62VUu2x/CJJBtKARkqpekAcMAwYcauv5xCcXfEcPpm6900gc+ZEEpZFcmHcs/i0rk/IuzNxb9DQ6AodRuGVBFK//ZbUpUsxZ2bi3bkzgf8Yi1e7dtKdUEEyCzL59MCn/HDqB4I9g/mo20f0rttb/v1vUpmhr5T6HugOBCmlYoHJgCuA1nou8AAwTilVBOQCw7TlE5UipdQEYD3gDCzUWh8rl7OwU8o7kGqvLcFn5AFS3p9I8vYoogcNImBQT4L/OVVu7ilHeadOk/LVV6SvXg1FRfjecw+BY8fiefv/DGAT5URrzbrz6/gw8kNS8lIY2Wwk4+8Yj4+bj9GlVWlyc1ZVoTVFu74j8d/vkXbSjJO7C0FPPkbA2Ik4ubsbXZ1d+GvYZcpXX5G9YwfK0xP/+++n+qOP4BYebnR5DuVCxgWm7pnKzvidtAhswRud3qBFYIuyD3RgckeuvSrMJe/Hd0j88j9kxbniUt2LkBdfpdrfhqKcbDUYy7HoggLSV68h5auvyD99GufgIKqPHEXAsL/Lu6kKVmAqYOHRhSw4vAA3ZzeeafMMDzV+CGcnGchQFgl9e5ceR/b8Z0lYfpC8VDfcI2oQ+sa7eHfubHRlVYYpLY3U//xI6rffUpSQgHujhlR/7HGqDRooi40bYO+lvbyz+x3OZ5ynX0Q/Xmr3EiFeIUaXVWVI6DsIfX4nGTMnkbgtlcJsF7zbtSTkX2/j0bSp0aVVWnnHj5OyZAkZq1aj8/LwvqsT1R9/Au8uneXDQQMk5ybz0b6PWBm9kto+tXm94+t0DpPGy42S0HckZhPmyG9Inf0BSQcV5kInqvXuRdCkF3CvLzN5gqULJ2P9BlKXLCH34EGUpyd+AwcSMHKE/II0iFmbWX5mOR/v/5icohyevO1JRt8+Gg8XD6NLq5Ik9B1RXjqmte+S/O0PpJzyQpuc8LtvEEETJuBWp07Zx9uhwsuXSf3hB9J+/AlTUhKudcMJGD4c/yFDZPIzA51KOcU7u9/hUOIh2tVox+sdX6e+nywteisk9B1Z0hmKlr1E8pr9pEb5oHHC//6hBI17CtdatYyurtzpoiKytm8nfdkyMjdvAbMZn65dCRg1Eu/OneUDbwPlFOYw59AcFh9fjJ+7Hy+2fZGB9QdKt5oNSOgLOL2BwmWvkLwzibRoX3Byxv+hvxM4+klca9rfjIQFMTGkLV9B+ooVFCUk4Fy9On5/+xsBw4c57DudymTzhc28v/d9LmdfZmijoTx353P4ucu7LVuR0BcWRQWwezaFa6aTdNiVtGhPcHLGb/B9BI4ejXu9qt3nb87JIXPjRtJ+WkZOZCQ4OeF9dxf8hw7Ft3t3lIzCMdylrEu8v/d9tlzcQqOARrzZ8U3uCLnD6LLsjoS++L8yLsFvkync9RPJ0aGknXFBF5rw7duXoLFj8Gje3OgKraYLCsjavoOM1avJ3LIFnZODa3g4/vffj9+Qv+EaGmp0iQIwmU18f/J7Pj34KQDjWo1jVPNRuDrJXDnlQUJflO7CbljzIkXnj5FypSmphwswZ+fg1bYtAaNG4turF8q18v1Q6qIicvbtJ2P1ajI2bMCcno6znx++/frhN/BePNu2lX7hSuRkykne3vk2R5OP0iWsC290fINaPvb/eZKRJPTFtZlNcOAb2DQFU0YGafndST2YTmFcPC4hIfg/MJRqAwcZPtzTlJVF9vbtZG7eTPbv2zClp6O8vPDt1Qu/gffifdddlfIXlCPLLcplzqE5LDq2CD93P/7Z/p/0jegrv5ArgIS+KFtOCmyZCvsWot0DyAocQeqey2Tv3AVmMx63347foIH49OhRIR+E6qIi8o4dI3vvXnJ27SY7MhIKC3H298enWzd8evbEp+vdOHl6lnst4sbtjN/JO7veITYrVj6oNYCEvrDe5SOw5mW4sBNq3kFhh9fJOBhH+spfyT9+AgC3unXx7toV704d8WzVCpfAwFt+2cKEBPKOHyfv2DFyDx8md99+zNnZltdr2ACfbt3w7dkTzzvukEVkKrGUvBSmR05nZfRKIqpF8GanN2lXo53RZTkcCX1xY7SGo8tgw+uQeRnaPgG93qTgShpZ2/4ga/sf5OzZi87PB8A1LAz3Zk1xq1sXt/C6uNYIxcnXFydvH5y8PMFkQhcWYs7NpSglBVNKKkUJCRRcuEDBhRgKzsdgSkqyvLZSuNWvj1e7tnh36IBXu3a4BAUZ+I8hrKG1ZmX0Sv4d+W+yCrN48rYnGdNyDO7OMuurEST0xc3Jz4TNU2HvPPAKgn7vw21DQSnMeXnkHT1K7qHD5B4+TH5UFIUXLqALC63+9i7BwbjWDcctvC4eTZvg0aIF7k2a4uwja8pWJRcyLjBl9xT2XNrDHcF3MLnTZBoGyOI+RrJZ6CulFgIDgQSt9f+sIKGUGgm8UvwwCxintT5UvO08kAmYgCJrCgIJ/Uoh/k9YNQniD0L9HnDvRxDY4H9202YzRVeuUHTlCqbsbMxZ2ZhzclAuzihXV5S7Oy7Vq+NcvTougYE4eXlV+KkI2yk0F/LNsW+Ye2gurk6uPHfnczzQ+AGclNzlbDRbhn5XLGG+6BqhfxdwQmudqpTqD7ylte5QvO080FZrnXQjxUvoVxJmE0R+CZumgLkQer4OHZ8GmdvcIR1JPMJbu97idOpp7gm/h392+KdMfVyJWBv6ZS6XqLXeppSKuM72nSUe7sayALqwB07O0GEsNBsEq5+39PefWAmDZ0OQvJV3FNmF2cw6OIslJ5YQ7BXMJz0+oVd4L6PLEjfJ1u/JngTWlnisgQ1Kqf1KqbHXO1ApNVYptU8ptS8xMdHGZYlbUq0mDFsCQ+ZD4imY2xl2fW55JyDs2rbYbQz+eTBLTixhWNNh/DL4Fwn8Kq7Mlr61lFI9sIR+lxJPd9ZaxyulQoCNSqmTWuttpR2vtZ4PzAdL946t6hI2ohS0+jvU7wYrJ8H61+DUWrh/geWXgrAraXlpTIucxqroVTT0b8j0btNlvhw7YZOWvlKqJfAFMFhrnfzX81rr+OK/E4AVQHtbvJ4wkG8NGP69pYsnbj/M7QJRm4yuStjQxpiNDP5lMOvOreOpVk/xw8AfJPDtyC2HvlIqHFgOPKy1Pl3ieW+llO9fXwN9gKO3+nqiElAKWo+EMVvAOwi+HQp75htdlbhFSblJPL/1eZ7f+jyhXqEsHbiU8XeMx81ZZiq1J2V27yilvge6A0FKqVhgMuAKoLWeC7wJBAKzi+fX+GtoZiiwovg5F2CJ1npdOZyDMEpIUxizGZaNhrUvQXYi9PyX0VWJG6S1ZlX0KqZFTiO3MJdn2zzLYy0ew8XJZr2/ohKRm7PErTMVWcb0H1wMfd+DTuONrkhY6XL2ZabsmsIfcX/QKrgVUzpPkWULqyibDdkUokzOLjBoJuRnwPp/Qc07IKKz0VWJ69Ba89OZn5ixbwYmbeKVdq8wvOlwnOUeDLsnoS9sw8nZ8uHupUPwy3iYEAnOMu1xZXQx8yJv73ybPZf30L5Ge9666y3q+Mpyko5C7p0WtuPuA/2mQeo5OLrc6GrEVbTW/HDyB4b+OpSjyUd5s9ObfNHnCwl8ByMtfWFbjfuCby04vc4yrl9UCpezLzN552R2xu+kU81OTOk8hRreNYwuSxhAQl/YllIQ2gJSoo2uRGBp3a8+t5r39rxHkbmI1zu8zkNNHpKVrByYhL6wPVOBTMpWCaTkpfDu7nfZGLORO4LvYGqXqYRXCze6LGEwCX1hW6ZCuPQntBhidCUObcuFLby16y0yCzKZ1GYSj7V4TEbmCEBCX9jayVWQlw6N+xldiUPKLMhk2t5p/HL2F5oENGFBnwU0DmhsdFmiEpHQF7ZjKoTf/w0BEdCoj9HVOJw9l/bwxo43uJJzhTG3j2Fcq3G4yrBZcRUJfWE7uz6HhGPw9++kT78CFZgK+OTAJyw+vpiIahEs7r+YlsEtjS5LVFIS+sI2LuyGze9C04HQbKDR1TiM6LRoXt72MqdSTzGsyTCeb/s8ni6eRpclKjEJfXHr0uPgP4+Afx0Y/LnR1TiEv6ZR+HDvh3i6ePJZz8/oVqeb0WWJKkBCX9yarARYdB8U5sLDP4Onv9EV2b20vDTe2vUWmy5solPNTkztMpVgr2CjyxJVhIS+uHk5KbDob5ARD6OWQ2hzoyuye5GXI3n1j1dJyUvhxbYv8nDzh3FSMpuKsJ6Evrg5GfGw+H7LnbcjfoC6nYyuyK6ZtZmFRxcy6+As6vjW4dsB39IisIXRZYkqSEJf3LikKFg8BHJTYOSPlnVzRblJz0/nte2vsS12G/0j+jP5rsl4u3obXZaoosp8X6iUWqiUSlBKlbrUobL4VCkVpZQ6rJRqU2JbP6XUqeJtr9qycGGQi3thYV8ozIHHVkngl7OjSUd5aOVD7IzfyWsdXmNa12kS+OKWWNMZ+DVwvdsr+wONiv+MBeYAKKWcgc+LtzcHhiulpNO3KvtzCXx9r2UK5SfWQ63WRldk15afWc4jax9Bo1nUbxHDmw6XidLELSuze0drvU0pFXGdXQYDi7Rl3cXdSil/pVRNIAKI0lpHAyillhbve/yWqxYVy2yCjW/Crs+gXld48Bvwqm50VXbLZDbx0f6PWHx8MZ1qduLDrh/i7+FvdFnCTtiiTz8MuFjicWzxc6U93+Fa30QpNRbLOwXCw2UmwEojNw2WPQlRv0H7sZY1cOXW/nKTVZDFy9te5o+4PxjRdAQvtXtJFigXNmWL/02lvd/U13m+VFrr+cB8sCyMboO6xK26fAT+8yikxcDAT6Dt40ZXZNcuZl5k4qaJnM84zxsd3+ChJg8ZXZKwQ7YI/Vig5HprtYF4wO0az4vKTms48A2seRk8A+DRlVD3LqOrsmunUk7xj43/oNBcyLze8+hQ85pvioW4Jba4q+NX4JHiUTwdgXSt9SUgEmiklKqnlHIDhhXvKyqz/CxYPhZWPmsJ+qe2S+CXs0OJh3h8/eM4OzmzuP9iCXxRrsps6Sulvge6A0FKqVhgMuAKoLWeC6wBBgBRQA7wePG2IqXUBGA94Aws1FofK4dzELZy5Tj8+CgkR0GP1+HuF8BJ7vYsT7sv7eaZzc8Q5BnEgj4LCPMJM7okYeesGb0zvIztGhh/jW1rsPxSEJXZX905a18Fj2rwyC+WUTqiXP2Z8CcTN02kTrU6zO89nyDPIKNLEg5AhgU4uuxkWPmMZcWr+t1hyHzwDTW6KrsXlRrF+E3jCfEKkcAXFUpC35FFbYKfx0FuKvSZCh2flu6cCpCal8q4TeNwd3ZnXu95EviiQknoO6LCPNj0NuyeDcFNYdQyqHG70VU5BLM289r210jOTWbxgMXU9q1tdEnCwUjoO5orx2DZGMuyhu3/Ab3fBldZaamirDy7ku1x23mtw2syS6YwhIS+ozCbYPcc2DQFPPxg5E/QqLfRVTmUv9aybRnckr83+bvR5QgHJaHvCJLPwi/j4cIuaDIABn0KPrLSUkXbELOBpNwkpnaeKgufCMNI6Nszsxkiv4DfJoOTK/xtLrQaBjJToyE2nt9ILe9adKzV0ehShAOT0LdXqTGW1v35P6DhPZbWvZ/c+GOko8lHaVejnbTyhaEk9O2N1rD/K9jwBqDgvlnQ+mFp3VcCOYU5+Ln5GV2GcHAS+vYkPRZ+mQDRWyw3Wt33GfjXKfMwUTFCvEKIy4ozugzh4CT07YHWcPBbWP+aZZTOvTOg7RPSuq9kOtbsyI+nfyQpN0luyBKGkc7Fqi49FpY8BL9OgBotYdwOaPekBH4lNLzpcDSa6fumG12KcGAS+lWV2QyRX8LnHeH8dug3zTLvffV6RlcmriHCL4Kxt49ldfRqlp5canQ5wkFJ905VlHwWfn0GYrZb+u4HzYSACKOrElYY23Isx5KP8cHeD/B396dfvX5GlyQcjLT0qxJTEeyYCXPusixleN9n8PDPEvhViLOTM9O6TuOOkDt45Y9X+PH0j0aXJByMhH5VcfkIfNELNr5pGXc/fg+0kaGYVZG3qzeze82mU61OTNk1hXd3v0uhqdDosoSDsCr0lVL9lFKnlFJRSqlXS9n+klLqz+I/R5VSJqVU9eJt55VSR4q37bP1Cdi9onzY/C7M7w4ZcfDg1/D3b6FaTaMrE7fAy9WLz3t+zuMtHueHUz/w2LrHiMmIMbos4QCUZeGr6+yglDNwGuiNZRH0SGC41vr4NfYfBDynte5Z/Pg80FZrnWRtUW3bttX79snvBy7utYy7TzoFrYZD3/fAq7rRVQkbW39+PVN2TaHQXMiLbV/kwcYPouQdnLhBSqn9Wuu2Ze1nTUu/PRCltY7WWhcAS4HB19l/OPC9dWWKUuVnWZYu/LIPFObAyGUwZK4Evp3qG9GX5fctp1VwK97Z/Q5Pb3qaS1mXjC5L2ClrQj8MuFjicWzxc/9DKeUF9AOWlXhaAxuUUvuVUmOv9SJKqbFKqX1KqX2JiYlWlGWnon6DOZ1gzxxoNxqe3gWN7jG6KlHOQr1Dmdd7Hq+2f5X9V/Yz+JfBfHfiO0xmk9GlCTtjTeiX9j7zWn1Cg4AdWuuUEs911lq3AfoD45VSpa64rbWer7Vuq7VuGxzsgNP+ZiXCstHw7VBwdofH18K908Hd1+jKxE0ymzXpOYUkZORxOT2P+LRc0nIKMJlL//FxUk6MbDaSFYNX0CakDR/s/YBH1z1KVGpUBVcu7Jk14/RjgZITuNQG4q+x7zCu6trRWscX/52glFqBpbto242Xaqf+mkJhw+uWrpzu/4Quz4GLu9GVCSvkFZo4fimDM1cyiU7K5lxiNjHJOSRn55OaU1hqwCsFvu4u1PDzoF6QN/WCfGhW05c76wYQ5u9JmE8Yc+6Zw6roVXwY+SEPrnqQMbePYfTto3FzdjPgLIU9seaDXBcsH+T2AuKwfJA7Qmt97Kr9/IBzQB2tdXbxc96Ak9Y6s/jrjcAUrfW6672mw3yQm3QGVj1nmf44/C4Y9AkENzG6KnEdVzLy2BGVxIELqRy6mM6JSxkUFQe7m7MT4YFeRAR6EezrQaC3G/5erni6OeOkFArILjCRnltIek4BcWl5nE/O5kJyDgUmMwAhvu50rB/IPc1D6dY4GJPKZNreaaw5t4Zw33D+2eGfdAnrYuC/gKisrP0gt8yWvta6SCk1AVgPOAMLtdbHlFJPFW+fW7zrEGDDX4FfLBRYUTwSwQVYUlbgO4SiAtjxCWz7t2V92kGfWqY/dpLbJiqb3AITO88msT0qiR1RSZy+kgVYWuq31/ZjTNf6tKrtT/Oa1QgL8MTZ6cZH3RSZzJy8nMnBC6nsi0ll+5kkfj0Uj4uTomP9QIa0nkjf7gP5+MCHjPttHL3Ce/Fyu5ep5VPL1qcrHECZLX0j2HVLP2YXrHzWMgyzxf3Q7wPwDTW6KlFCek4hm05eYf2xy/x+OpG8QjPuLk60r1edLg2D6NwwiOY1q+F0EwFvDZNZ8+fFVDYeT2DNkUtcSMnB282ZfrcFUa3GLlZdXAQaxrQcw2MtHpMuHwFY39KX0K8ouWnw21uWBU78wuHej6BxH6OrEsXScwtZd/QSqw5fYtfZZIrMmtBq7vRpXoM+LUJpF1EdD1fnCq9La03k+VR+2n+R1YcvkV1gokUdE75hazmWvl26fMR/SehXFlrD8Z9h7SuQnQgdn7Z8WOvuY3RlDq+gyMzWUwn8/Gccv51IoKDITESgF/1uq0nfFqG0qu1fbq35m5GdX8TyA7F8tfM80YnZVA86h1fNlaQXxdOjTg9eavcSdXxl0RxHJaFfGaRdhDUvwul1ULOVZTbMWq2NrsrhHYlN54d9F1h1+BJpOYUEersxqFUthrQOo2Vtv0p/N6zZrPkjKomF28/x+5lL+IbsxCVwM07KzKMtHmX07aPxcvUyukxRwST0jWQ2wZ55ljlz0NDjX9DhKXCWmayNkpVfxK9/xvP93gsciUvHw9WJPs1rMKR1GF0aBeHqXDU/RD8al87nW6JYd/I0XjXW4eR7gCCPEF5u/yL9IvpV+l9gwnYk9I0SfxBWToJLf0LD3pa++4C6RlflsI7EprNkbwy//BlPToGJpjV8GdEhnMF3hOHn6Wp0eTZz+koms7dEser0TjxCV6I84mgV1JrXO71G0+pNjS5PVAAJ/YqWlw6bp0LkAvAKgn7vw21DZepjA+QWmFhxMI4le2M4GpeBh6sTg1rWYniHcFrX8bfr1m90YhYzNp5i/YVf8QhZj3LOZUjDoUy6cyIBHgFGlyfKkYR+RfnvB7WvQtYVy/q0Pd8AT3+jK3M4cWm5LNp1nqV7L5KeW2i3rXprHItP54P1B9ib9gNuAbvwcPbimTYTGN7s77g4STejPZLQrwgp0bDmJcskaTVawsBPoPadRlflULTW7ItJ5asd51h/7Apaa/q2qMHjnevRLiLArlv11th7LoV3N2zmjOlbXLzPEuoRwdSur9OhZgejSxM2JqFfnoryYeensG06OLlAz9eh3Rj5oLYC5ReZWHXoEl/tPMfRuAyqebgwvH04D3eqS+0AGblSktaaLacSeHvTDyS5/YiTWxov3D6dx9r0Nbo0YUMS+uXl3B+w+nlIOg3NB1vuqK0mt8NXlNTsAhbvjmHRrhiSsvJpGOLDY3dFcH+bMLzc5Jfu9ZjMmqX7zvLRrm9JuXQnfZrX5NX+TakfLPeM2AMJfVvLTrLMhHnoe/APhwFyR21Fik3N4cvt51i69yK5hSa6Nwnmic71uLtRkMN34dyo3AITX26PZs7Ws+QXmRnRIZxnezUi0Edmdq3KJPRtxWyGg4stC5IXZEPnZ+DuF8FNuhAqwvH4DOZvO8vKw5dQwOA7whjbtT5Nasg6A7cqMTOfmZtO8/3ei3i6OjOuewOe6FwPT7eKn25C3DoJfVu4cswy9fHFPVC3M9w7A0JkzHN501qz62wyc7dFs+10It5uzozoEM7jnetRy9/T6PLsTlRCFtPWnWTj8SvU9PPghT5NuL91WKWagkKUTUL/VhRkw9YPYNfn4OEHfadaFiaXboRyZTJr1h29zLxtZzkcm06QjzuPd45gVMe6Djfk0gh7opN5b80JDsWmc1tYNd64tzkd6gcaXZawkoT+zTq5Bta+DOkXLXPc954iC5KXs/wiEz/tj2Xe79FcSMmhXpA3Y7vWZ0jrMENmtnRkZrNm5eF4pq09SXx6HgNur8E/+zejTnXpzqzsbLaISvE36wfMxLKIyhda6w+u2t4d+AXLylkAy7XWU6w5ttJIu2iZCfPUaghuBo+vg7qdjK7KruUWmPh+7wXmbTvLlYx8WtXx57UBTendvMZNLUYibp2Tk2LwHWH0aV6DBX9YPuz97XgCT3Spx/geDfD1kHdcVZ01yyU6Y1kusTeW9XIjgeFa6+Ml9ukOvKi1Hnijx5amQlv6pkLYPQe2vm+5u7b7q9BpPDjLf+7ykplXyOLdMXz5xzmSswvoWL86E3s24q4GgTISp5K5nJ7Hh+tPsvxAHEE+brzYpwkPtq0jv5QrIVu29NsDUVrr6OJvvBQYDFw3uG1wbPm7uNfyQe2Vo9C4H/T/UCZHK0dpOQV8teM8X+04R0ZeEd0aBzOhZ0PaRUj3WWVVw8+DGQ/dwaOdInhn1XFeXX6Eb3bF8MbAZtzVIMjo8sRNsCb0w4CLJR7HAqXdw91JKXUIiMfS6j92A8dWvOxk+GYQeAXC37+DpvfKB7XlJCkrny/+OMe3u2PIyi+iT/NQJvRsSMva/kaXJqzUqo4/Pz7VidVHLvH+mpOMWLCHPs1DeW1AMyKCvI0uT9wAa0K/tCS8uk/oAFBXa52llBoA/Aw0svJYy4soNRYYCxAeHm5FWbfIOxCGLYE6HWQVq3JyOT2PedvO8v3eCxQUmbm3ZS3G92hA0xrVjC5N3ASlFANb1uKeZqF8uf0cn2+JovfHv/N453pM6NmQatLfXyVYE/qxQMk12Gpjac3/l9Y6o8TXa5RSs5VSQdYcW+K4+cB8sPTpW1X9rWrYq0JextHEp+Uye2sU/4mMxaw1Q1qHMa57A7nd3054uDozvkdDHryzNtM3nGLBH9Es2x/Lc70bM6xdHVyq6II0jsKaD3JdsHwY2wuIw/Jh7Iji7pu/9qkBXNFaa6VUe+AnoC6WETvXPbY0ho/TFzflUnous7ec5YfIi2g0D7atw7huDWS4n507GpfOlFXH2XsuhSahvrwxsDldGkl/f0Wz2Qe5WusipdQEYD2WEF+otT6mlHqqePtc4AFgnFKqCMgFhmnLb5NSj73psxKV0uX0POZsjeL7vf8/7J/u3kBmu3QQt4X58cPYjqw7epn31p5g1Jd7uKdZCK8NaCbv7iohuTlL3LQrGXnM2XqWJXsvYDZrHmxbm/E9GkrYO7C8QhNf7zzPZ5ujyCs08dhdEUzs1UjuqK4AckeuKDcJGXnMLhH2D9xpCXvpxhF/SczMZ8bGUyyNvEiAlxsv9GnMsHbhMr6/HEnoC5tLyMxj7tZovtsTQ5FZM7RNGBN6NCI8UMJelK5kf3+zmtV4c2BzOjWQ+XzKg4S+sJnEzHzm/n6Wb3dbwv7+1mFM6NmQuoEyPluUTWvNmiOXeW/NCeLScul/Ww1eGyDz+diaTefeEY4pMTOf+dvOsnh3DAVFZoa0rs3Eng3lZhxxQ5RS3NuyJr2ahbBgWzSzt55l08kERnepx9M9GuLjLjFUkaSlL/5HWk4B87ZF8/WO8+QXmfhb6zAm9mxEPQl7YQOX0/OYtu4kKw7GEeLrzsv9msr8/TYg3TvihmXmFbJw+3m++COarIIi7mtVi2d7NZJhd6JcHLiQytsrj3PoYhqtavvx5qAW3Fk3wOiyqiwJfWG13AITi3efZ87Ws6TmFNKneSjP92ks0yWIcmc2a37+M45p605yJSOfv91Ri1f6N6Wmn6yQdqOkT1+UKb/IxA+RF/lscxQJmfl0bRzMC70b06qOv9GlCQfh5KS4v01t+raowZytZ5n/RzTrj13hqW4NGNu1vqzXWw6kpe+Aikxmlh+IY+amM8Sl5dI+ojov9GksS+MJw11MyeH9tSdYc+QyYf6evNq/KQNb1pR1Fqwg3Tvif5jNmlVHLvHJxtNEJ2XTsrYfL/Zpwt2NguSHSlQqu6OTeXvlcU5cyqBdRACTB7XgtjA/o8uq1CT0xX9prdl4/AozNp7m5OVMmoT68nyfxvRpHiphLyotk1nzn30Xmb7+FCk5BTx0Zx1e7NuEYF93o0urlCT0BVprtkclMX3DaQ5dTKNekDeT7mnEoJa1ZHicqDIy8gqZtekMX+88j7uLMxN6NuTxzhG4u0h/f0kS+g4u8nwK/15/ir3nUgjz9+TZXo24v02YzHUuqqzoxCymrj7BppMJ1A304l8DmtFb3q3+l4S+gzocm8ZHG07z++lEgn3dmdCjIcPa15FWkbAbv59O5J1Vx4lKyKJzw0DeHNiCJjV8jS7LcBL6DiYqIYuPNpxi7dHL+Hu5Mq5bAx7pFCFD3oRdKjSZ+W53DB//dobMvEJGdqjL870bE+DtZnRphpHQdxDxabnM/O0MP+6/iKerM6Pvrs/ou+vhK+uVCgeQml3Ax7+d5rs9F/Bxd2HSPY0Y1bEurg7YjWnT0FdK9QNmYln96gut9QdXbR8JvFL8MAsYp7U+VLztPJAJmIAia4qS0C9banYBs7dG8c2uGNDwcKe6PN29AYE+MrJBOJ5TlzN5Z9Vxtkcl0TDEhzcGNqdb42Cjy6pQNgt9pdRf69z2xrLQeSQwXGt9vMQ+dwEntNapSqn+wFta6w7F284DbbXWSdYWL6F/bdn5RXy5/RwLtkWTXVDE0Da1efaeRrJalXB4Wmt+O5HAu6uPE5OcQ6+mIfzrXsdZstGW0zC0B6K01tHF33gpMBj4b+hrrXeW2H83UPvGyhVlyS8y8f2eC3y2JYqkrAL6tgjlxT5NaBQqH2AJAZYpnHs3D6Vr4yC+3nGeWZuj6PvJNh7tJEs2lmRN6IcBF0s8jgU6XGf/J4G1JR5rYINSSgPztNbzSztIKTUWGAsQHh5uRVmOwWTW/PJnHDM2niY2NZeO9asz/5GmtAmX2QiFKI27izP/6NaA+9vUZvr6U3y54xwrDsbxQp8m/L1dHYdfstGa7p0Hgb5a69HFjx8G2mutJ5aybw9gNtBFa51c/FwtrXW8UioE2AhM1Fpvu95rSveO5a3qphMJ/Hv9KU5dyaRFrWq80q+pTJkgxA06GpfO2yuPEXk+lWY1qzF5UHM62uE8U7bs3okF6pR4XBuIL+UFWwJfAP3/CnwArXV88d8JSqkVWLqLrhv6jm7vuRSmrTvJ/phU6gV589mI1gy4rabcRSvETbgtzI///KMTq49c4v01Jxk2f7dDL9loTehHAo2UUvWAOGAYMKLkDkqpcGA58LDW+nSJ570BJ611ZvHXfYAptire3hyPz+Df60+y5VQiodXceW/I7TzYtrZDDj8TwpaUUgxsWYt7moUyf1s0c4qXbBxzdz2e7t4QbwdastHaIZsDgE+wDNlcqLWeqpR6CkBrPVcp9QUwFIgpPqRIa91WKVUfWFH8nAuwRGs9tazXc7TunZjkbGZsPM0vf8ZTzcOFp3s05FG5sUqIcnMpPZdpa0/y85/xhPi680q/pgyp4ks2ys1ZVUBCRh6fbj7D0r0XcXFWPNG5Hv/o2gA/LxllIERF2B+TypSVxzgUm06rOv5MHtS8yg6SkNCvxNJzC5n3+1kW7jhHkUkzvH04E3s2JKSah9GlCeFwzGbNioOWJRsTMqvuko2yXGIllFdo4pud55m99SzpuYUMvqMWz/duTN1Ab6NLE8JhOTkpht5Zm3631WD21igW/HGO9ceuMK67ZclGD1f76maVln4FMJk1yw7E8vHG01xKz6N7k2Be6tuEFrVkJSAhKpuLKTlMXX2CdccsSza+NqAZA26vUemHSkv3TiWgtWbLqQSmrbWMtW9V249X+jflrgZBRpcmhCjDrrPJvL3yGCcvZ9K+XnUmD2peqRtqEvoGO3AhlQ/WnmTvuRQiAr14qW/TKtFaEEL8fyazZmnkBaavP0VabiHD2tXhhT5NCKqEExtK6BvkbGIW09db5rUP8nHn2XsaMaxdHRlrL0QVlp5TyMxNZ1i06zyers48e08jHukUgZtL5fm5ltCvYAkZeczcdIalkRfxcHFibNcGjL67nkPd9CGEvYtKyOLd1cfZeiqR+kHevD6wGT2ahFSKd/AS+hUkM6+Q+dui+eKPcxSazIzsEM6Eno0I9q18b/+EELax5WQC76w6TnRSNt0aB/PGwGY0DDF2xlsJ/XJWUGTmuz0xzNocRUp2AQNb1uTFPk2ICJLhl0I4goIiM4t2nWfmpjPkFJh4pFNdJvVqbNjNlRL65cRs1qw8HM9HG05zISWHuxoE8mr/prSs7W90aUIIAyRn5fPRxtN8v/cC/p6uvNCnCcPa1cGlgj/Hk9AvB9vPJPHBuhMcjcugWc1qvNq/KV1lqmMhBHAsPp0pK4+z51wKTWv48ubA5tzVsOKGZ0vo29DRuHSmrTvJH2eSCPP35MW+jRncqmpPziSEsD2tNWuPXmbq6hPEpeXSr4VlCufwwPKfwlmmYbCBiyk5fLThFD//GY+/lyuv39uMUR3r2t1t2UII21BKMeD2mvRsGsIXf0Tz+ZazbD6VwOgu9Xi6R0N8KsFoPmnplyIlu4BZm8/w7e4YnJ0ss18+1b0B1Txk9kshhPUup+fx4bqTLD8YR4ivOy/3a8r95TSFs3Tv3IScgiIWbj/HvN+jyS4o4qG2dZh0T2Nq+Mnsl0KIm3fgQipvrzzOoYtptKrtx+T7Wth8CmcJ/RtQZDLzn32xfPLbaRIy8+ndPJRX+jUxfNytEMJ+XD2F85DWYbzSr6nNGpXWhr5VY4qUUv2UUqeUUlFKqVdL2a6UUp8Wbz+slGpj7bFG0lqz7uhl+nyyjddWHKFOdS9+eqoTCx5pK4EvhLCpv6Zw3vJid8b3aMDqI5foMX0rszadIa/QVGF1lNnSV0o5A6eB3lgWSY8Ehmutj5fYZwAwERgAdABmaq07WHNsaSqipZ+aXcCT30Ry4EIaDYK9eaVfU3o3D5Xhl0KICnEhOYf31vz/KZz/dW8z+t9285My2rKl3x6I0lpHa60LgKXA4Kv2GQws0ha7AX+lVE0rjzWEv5crodU8+OD+21k/qSt9WsgMmEKIihMe6MXch+9kyZgO+Hq48PR3Bxg2fze5BeXb6rdm/FAYcLHE41gsrfmy9gmz8lgAlFJjgbEA4eHhVpR1a5RSzBl1Z7m/jhBCXM9dDYJY/czdLI28wJHYdDzdyndIuDWhX1rz9+o+oWvtY82xlie1ng/MB0v3jhV1CSGEXXB2UozsUPcaTWLbsib0Y4E6JR7XBuKt3MfNimOFEEJUEGv69COBRkqpekopN2AY8OtV+/wKPFI8iqcjkK61vmTlsUIIISpImS19rXWRUmoCsB5wBhZqrY8ppZ4q3j4XWINl5E4UkAM8fr1jy+VMhBBClEluzhJCCDtg05uzhBBC2AcJfSGEcCAS+kII4UAk9IUQwoFUyg9ylVKJQEwFvVwQkFRBr1UZONr5guOds6OdL8g5A9TVWgeXdVClDP2KpJTaZ80n3vbC0c4XHO+cHe18Qc75Rkj3jhBCOBAJfSGEcCAS+sWTvDkQRztfcLxzdrTzBTlnqzl8n74QQjgSaekLIYQDkdAXQggH4hChfysLu1dVVpxzd6VUulLqz+I/bxpRp60opRYqpRKUUkevsd0er3FZ52xv17iOUmqLUuqEUuqYUurZUvaxm+ts5fne+DXWWtv1HyxTOp8F6mNZ1OUQ0PyqfQYAa7Gs9NUR2GN03RVwzt2BVUbXasNz7gq0AY5eY7tdXWMrz9nernFNoE3x177AaXv+WbbyfG/4GjtCS/9WFnavqirtgvTlRWu9DUi5zi72do2tOWe7orW+pLU+UPx1JnACyzrcJdnNdbbyfG+YI4T+tRZtv9F9qhJrz6eTUuqQUmqtUqpFxZRmGHu7xtayy2uslIoAWgN7rtpkl9f5OucLN3iNrVkjt6q7lYXdqyprzucAlrk6spRSA4CfgUblXZiB7O0aW8Mur7FSygdYBkzSWmdcvbmUQ6r0dS7jfG/4GjtCS/9WFnavqso8H611htY6q/jrNYCrUiqo4kqscPZ2jctkj9dYKeWKJQC/01ovL2UXu7rOZZ3vzVxjRwj9W1nYvaoq85yVUjWUUqr46/ZY/i8kV3ilFcfernGZ7O0aF5/Ll8AJrfWMa+xmN9fZmvO9mWts9907+hYWdq+qrDznB4BxSqkiIBcYpouHA1RFSqnvsYxkCFJKxQKTAVewz2sMVp2zXV1joDPwMHBEKfVn8XOvAeFgl9fZmvO94Wss0zAIIYQDcYTuHSGEEMUk9IUQwoFI6AshhAOR0BdCCAcioS+EEA5EQl8IIRyIhL4QQjiQ/we4C/KMXmNsIQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig = plt.figure()\n", "ax = fig.add_subplot(1, 1, 1)\n", "data1, data_t = simulate(1500, state0=[0,0.01], h= 0.02)\n", "ax.plot(data1[:, 0], data1[:, 1])\n", "data2, data_t = simulate(200, state0=[0,0.5], h= 0.02)\n", "ax.plot(data2[:, 0], data2[:, 1])\n", "data3, data_t = simulate(200, state0=[2,1.5], h= 0.02)\n", "ax.plot(data3[:, 0], data3[:, 1])\n", "data4, data_t = simulate(140, state0=[2,2], h= 0.02)\n", "ax.plot(data4[:, 0], data4[:, 1])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### The text and the drawing in the book claim that trajectories leave the unstable node (0, 0) tangental to to the slow direction, here the y-axis. Below are trajectories for a very short time starting from (0.01, 0.01) and (0.0, 0.01)." ] }, { "cell_type": "code", "execution_count": 57, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 57, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAD4CAYAAADiry33AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAAAv3klEQVR4nO3dd3yN5//H8dclQ+zYK4m9VxB7dOGLqt3SraWqLVJare5WJx1Kh1bRb6cYsapWhw5qhSRICBErEiLIkMg81++P+/T7S1U5leRKzsnn+XicB8d9H9fnPpJ3bvd1n+ujtNYIIYRwXaWKugAhhBCFS4JeCCFcnAS9EEK4OAl6IYRwcRL0Qgjh4tyLuoArqVatmq5fv35RlyGEEE5j9+7diVrr6lfaViyDvn79+oSEhBR1GUII4TSUUsf/aZtcuhFCCBcnQS+EEC5Ogl4IIVycBL0QQrg4CXohhHBxEvRCCOHiJOiFEMLFSdALIURel5Jg/dOQkVzUlRQYCXohhPjT8W3wSU/YtQCObS3qagqMBL0QQuTmwOY34L8DoZQ7PLgJmg80WkJc0iU2RpwulL+7WC6BIIQQxlw4BsEPQexOaHcXDJwFpSsYLWHdvnimB+/F3a0UPRtXo1zpgo1mCXohRMm1dxl8P9X6/YiF0Gak0eHTMnOY8V0kS0JO0s6nEu+Pbl/gIQ8S9EKIkigjBdZNg71B4NsVhs+HyvWMlrA3NonAoDCOnUvj0RsbMaVvUzzcCudqugS9EKJkiQ2B4LGQdAJufAZ6PQlu5qLQZtPM/z2GdzZGUa18ab4d15VujaoW6pgS9EKIksGWC1veg81vQsW68MB68OtqtITTyRlMXRrGH0fOMaB1Ld4c3gbvsp6FPq4EvRDC9SWdhJUPw/Gt0HoE3PoelPE2WsLGiNM8HbyXzGwbM0e04Y4AX5RSRsaWoBdCuLaIVfDdZOuMfugn0G40GApYgPSsHF5de4DFO0/Qum5F5oxuT6Pq5Y2NDxL0QghXlZVmfcI19Cuo0wFGLICqjYyWsP9UMoFBoRw5m8bDNzTkib7N8HQ3//ElCXohhOuJC4XgcXDuCPR6wpp0dfMwNrzNplm09SizNkRRuZwH34zrQo/G1YyNfzkJeiGE67DZYNuH8NMMKFcd7v8OGvQyWkJCSgZPLAvn98OJ9G1Zk5kj2lKlXOFPuF6NBL0QwjWkxMOqCRDzC7S4DW6bC2WrGC3hpwNnmLZ8L+lZObw+rDV3dfYzNuF6NRL0Qgjnd3AdrH4McjKsgO9wn9EJ14zsXN5Yd4Avtx2nRe2KfHCnP41rmF1G4Wok6IUQziv7Emx63lptslZbaxmD6k2NlnAgPoXAoFAOnbnI2J4NeKp/M0q7uxmt4Vok6IUQzun0fmvC9ewB6DYRbnkR3EsbG15rzX//OMab6w9S0cuDLx/sTO+m1Y2N/29I0AshnIvWsHM+bHrB+tDTvSuh0c1GS0i8mMm0ZeFsjjrLzc1rMGtkW6qVN/dD5t+SoBdCOI+LZ2H1o3B4EzTtD0M+gnJmb1vcHJXAtGXhpGTkMGNIK+7tWq9YTLhejQS9EMI5HP4RVj1itfgb+A50Gmd8wnXmhoN8vvUYzWpW4JtxXWlWq/hMuF6NBL0QonjLyYQfX4HtH0GNlnDfaqjZ0mgJh86kMnlxKAdPpzKme32mD2iOl0fxmnC9Ggl6IUTxdTYKlo+FM/ug88PQ9xXwKGNseK01X+84wWtrIylf2p3Px3TipuY1jI1fUBxadEEp1V8pFaWUilZKTb/CdqWUmmvfvlcp1SHPNm+l1HKl1EGl1AGlVLeCPAAhhAvSGkIWwac3QGoc3LnEavFnMOTPp2Xx0Je7eWHVfro2rMr6x3s5ZciDA2f0Sik34COgLxAL7FJKrdFaR+bZbQDQxP7oAsyz/wowB9igtR6plPIEyhZg/UIIV5N+HtZMgoNrrbtphs6DCrWMlrDlcCJTl4aRlJ7NC4Na8kD3+pQqVbwnXK/GkUs3nYForXUMgFIqCBgC5A36IcCXWmsNbLefxdcG0oDewBgArXUWkFVw5QshXErMr9a68WmJ0O916PoolDK32mNWjo13NkUx/7cYGtcoz38f6EzLOhWNjV9YHAn6usDJPM9j+f+z9avtUxfIAc4Cnyul2gG7gUCtddrlgyilxgPjAfz8/BytXwjhCnKyYPPrsHUOVG0Mdy2B2u2MlnDk7EUmLw4lIi6Fe7r68dzAlpTxdJ4J16tx5Efllf6/oh3cxx3oAMzTWrfHOsP/2zV+AK31fK11gNY6oHr14vnpMiFEITh3BBb1g63vW2vUPPyr0ZDXWhO08wSD5m4hLukS8+/tyGtD27hMyINjZ/SxgG+e5z5AnIP7aCBWa73D/ufL+YegF0KUMFpD2Lewbpq1VvwdX0HLwUZLSErPYnrwPjZEnKZH46q8d4c/NSt6Ga3BBEeCfhfQRCnVADgFjAbuumyfNcBE+/X7LkCy1joeQCl1UinVTGsdBdzCX6/tCyFKoktJsHYKRKyA+r1g2KdQqa7RErYdOceUJWGcS8vkmQHNeahXQ6eecL2aawa91jpHKTUR2Ai4AYu01hFKqQn27Z8A64CBQDSQDjyQ56+YBHxjv+Mm5rJtQoiS5vg2WPEQpMZbC5H1eBxKmbtMkp1rY/YPh5j36xEaVC3HZ/f1oI1PJWPjFwVl3ShTvAQEBOiQkJCiLkMIUZByc+C3t+G3WeBdz1pS2Kej0RKOJaYRGBRKeGwyozv58uJtLSnr6RqfG1VK7dZaB1xpm2scoRCieLtw3DqLP7kD2t1lffiptLl1YrTWBO85xUur9+PuVop5d3dgQJvaxsYvahL0QojCtW+5dT0erLP4NiONDp98KZvnVu5j7d54ujaswnt3+FPH29wnbIsDCXohROHITLXuqAlfDL5dYPhnULme0RJ2HTvP40FhnEnJYNp/mjHhhka4ueiE69VI0AshCl7sbggeC0nH4Ybp0HsauJmLm5xcG3N/jubDnw/jW6Usyx/pjr+vt7HxixsJeiFEwbHlWh982vwGVKgDY9ZBPbPrGJ48n05gUCh7TiQxooMPrwxpRfnSJTvqSvbRCyEKTnIsrHgYjm+BVsNh0Gyr1Z9Bq0JP8fyq/SgFc+9sz+B2dYyOX1xJ0Ash8i9yNayZDLnZ1mqT7e402v0pNSObF1dHsDL0FAH1KvP+aH98KstCuX+SoBdCXL+sNNjwDOz5Auq0t+6qqdrIaAm7j1/g8SWhxCVlMKVPUx67qRHubuZWvHQGEvRCiOsTFwbB4+BcNPScAjc+C+6exobPtWk+2hzNnJ8OU7uSF0sf7krHelWMje9MJOiFEP+OzWb1b/3xFShXHe5fAw16Gy3hVNIlpgSFsfPYeYb41+HVoa2p6OVhtAZnIkEvhHBc6mlYOQFiNkPzQTD4Ayhr9ix67d44nlmxD61h9qh2DGvvY3R8ZyRBL4RwTNQGWP0oZKXDoPeh4xijE65pmTm8vCaCZbtj8ff1Zs5of+pVLWdsfGcmQS+EuLrsS/DDi7BzPtRsAyMXQvVmRksIP5lEYFAox8+nM+nmxky+pQkeMuHqMAl6IcQ/OxNpfcI1IRK6PgZ9XgL30saGz7VpPv3tCO9tOkSNCqUJeqgrXRpWNTa+q5CgF0L8ndawawFsfA68KsHdwdCkj9ES4pMvMXVJONtiznFrm9q8MawNlcrKhOv1kKAXQvxVWiKsfgwObYAm/WDIx1DebB/nDfvjeTp4H9m5NmaNbMvtHX1QBucDXI0EvRDi/x352bqr5lISDJgFnccbnXBNz8rh1bUHWLzzBG19KjFndHsaVJMJ1/ySoBdCQE4W/DwD/vgAqjeHe1ZArdZGS9h/KpnJQaEcTUxjwg2NmNq3KZ7uMuFaECTohSjpEg9bE67x4dBpHPR7DTzMNeaw2TQLtxxl1saDVCnnyTdju9C9cTVj45cEEvRClFRaw54vYcN0cPeC0Yuh+UCjJSSkZPDEsnB+P5xIv5Y1mTmiLZXLmVtGoaSQoBeiJEo/D98FwoE10OAGGPYpVDTbQ/XHyDM8FbyX9Kwc3hjWhjs7+8qEayGRoBeipDm2BVaMh4tnoO8M6DYJSpm7Fp6Rncsb6w7w5bbjtKxdkbl3+tO4hrlG4SWRBL0QJUVuNvzyFvz+LlRpCGN/gLodjJZw8HQKkxeHcujMRcb1bMC0/s0o7e5mtIaSSIJeiJLg/FFrSeFTIdD+Hug/E0qXNza81pov/jjGG+sPUtHLgy8e7MwNTc3em1+SSdAL4erCl8D3T4AqBSM/h9bDjQ6feDGTp5bv5eeDCdzcvAazRralWnlzyygICXohXFdGihXw+5aCXzcYPh+8/YyW8EtUAk8u20tKRjavDG7Ffd3qyYRrEXAo6JVS/YE5gBuwQGv91mXblX37QCAdGKO13mPfdgxIBXKBHK11QIFVL4S4spO7rHvjk2Phpueg51RwM3del5mTy6wNUSzccpRmNSvw9bjONK9V0dj44q+u+S+vlHIDPgL6ArHALqXUGq11ZJ7dBgBN7I8uwDz7r3+6SWudWGBVCyGuzJYLv78Hv7wJlerCA+vBr8u1X1eAohNSmbQ4jAPxKdzfrR7PDGyBl4dMuBYlR37EdwaitdYxAEqpIGAIkDfohwBfaq01sF0p5a2Uqq21ji/wioUQV5Yca902eXwrtB4Jg96zVp40RGvNtztP8OraSMp6urPw/gBuaVHT2PjinzkS9HWBk3mex/LXs/V/2qcuEA9oYJNSSgOfaq3nX2kQpdR4YDyAn5/Z64hCOL3I1bBmknVGP+xTaDvK6GJk59OyeDp4Lz9EnqFXk2q8e0c7alTwMja+uDpHgv5KXy36X+zTQ2sdp5SqAfyglDqotf7tbztbPwDmAwQEBFz+9wshriQrzVrCYM+XULcjjFhg3SNv0NboRKYuDeNCWjbP39qCB3s0oFQpmXAtThwJ+ljAN89zHyDO0X201n/+mqCUWol1KehvQS+E+JfiwqwJ13NHoNcTcOMz4GauMUdWjo33fjjEp78doWG1ciwa04lWdcxdKhKOcyTodwFNlFINgFPAaOCuy/ZZA0y0X7/vAiRrreOVUuWAUlrrVPvv+wEzCq58IUogmw22fQg/zYBy1eH+76BBL6MlxJy9SGBQGPtOJXNnZz9eHNSSMp4y4VpcXTPotdY5SqmJwEas2ysXaa0jlFIT7Ns/AdZh3VoZjXV75QP2l9cEVtrvm3UHvtVabyjwoxCipEg9bTUGidkMzQfB4A+gbBVjw2utWRYSy8vfReDpXopP7ulI/9a1jI0vro+ybpQpXgICAnRISEhRlyFE8RK1AVY/Clnp0P9N6DjG6IRrcno2z67cx/f74unWsCqzR/lTq5JMuBYXSqnd//Q5JflkrBDFXfYl2PQC7PoMaraBkQuhejOjJew8ep7Hg0JJSM3k6f7NGd+7IW4y4eo0JOiFKM7ORFoTrgmR0PUx6PMSuJtbJyY718bcnw7z0eZo/KqUJfiR7rTz9TY2vigYEvRCFEdaw87PYNPz4FUR7g6GJn2MlnDiXDqBS0IJPZHE7R19eHlwK8qVlshwRvKvJkRxk5YIqx+DQxugcV8Y+jGUr2G0hJWhsbywKgKlYO6d7Rncro7R8UXBkqAXojg58rN1V82lC9aa8V0eNjrhmpqRzQur9rMqLI5O9Ssze5Q/PpXLGhtfFA4JeiGKg5ws+HkG/PEBVGsG9wRDrTZGS9h9/AKPLwklLimDqX2b8uiNjXB3M9diUBQeCXohilriYWvCNT4cAsZCv9fA09xZdK5N89HmaOb8dJjalbxY+nBXOtYzd2++KHwS9EIUFa0h9CtY/7R1J83ob6H5rUZLOJV0iSlBYew8dp4h/nV4dWhrKnqZW0ZBmCFBL0RRuHQBvgu0Vp1s0NtacbKi2QnP7/fGM33FXrSG2aPaMay9j9HxhTkS9EKYdmyrtW78xdPQ5xXoPhlKmbsWnpaZw8trIli2OxZ/X2/mjPanXtVyxsYX5knQC2FKbjb8OhN+fxcq14exm6ylhQ3aG5tEYFAYx86lMfGmxgT2aYKHTLi6PAl6IUw4fxRWPASxu8D/HhgwE0qXNza8zaaZ/3sM72yMonqF0ix+qCtdG1Y1Nr4oWhL0QhS2vctg7RTrfvgRC6HNSKPDn07OYOrSMP44co6BbWrxxrA2eJf1NFqDKFoS9EIUlowUWDcN9gaBbxcY/hlUrme0hE0Rp3k6eC8Z2TZmjmjDHQG+KIMfwBLFgwS9EIUhNsS6Nz7pBNwwHXpPAzdz326XsnJ57ftIvtlxgtZ1KzJndHsaVTd3qUgULxL0QhQkWy5smQ2b37BulxyzDup1M1pCZFwKk4NCiU64yPjeDXmyXzM83WXCtSSToBeioCTHwoqH4fgWaDUcBs2GMt7GhrfZNJ//cYyZ6w/iXdaDr8Z2pleT6sbGF8WXBL0QBSFyDayZZN1COXQetLvT6GJkZ1MzeXJZOL8eOkufFjWYOaItVcubW7deFG8S9ELkR1YabHgG9nwBddpbd9VUbWS0hM0HE5i2PJzUjBxeHdqae7r4yYSr+AsJeiGuV3w4BI+zFiXrOQVufBbczd22mJGdy1vrD/LfP47RvFYFvn2oK01rVjA2vnAeEvRC/Fs2G+yYBz++DGWrwn2roOGNRks4dCaVyYtDOXg6lTHd6zN9QHO8PNyM1iCchwS9EP9G6hlY9Qgc+Qma3QqDP4By5j5hqrXm6x0neG1tJOVLu/P5mE7c1Nxs9ynhfCTohXDUoU2w+lHITIVb34OAB41OuJ5Py+Kp5Xv58cAZbmhanbdvb0uNCl7GxhfOS4JeiGvJzrAu0+yYBzVawf3fQY0WRkvYcjiRqUvDSErP5oVBLXmge31KlZIJV+EYCXohribhoPUJ1zP7ocsEa1lhD3Nn0Vk5Nt7dFMWnv8XQuEZ5Pn+gE63qVDI2vnANEvRCXInWELIINj4LnuXhrqXQ9D9GS4g5e5HAoDD2nUrm7i5+PH9rS8p4yoSr+PccCnqlVH9gDuAGLNBav3XZdmXfPhBIB8Zorffk2e4GhACntNaDCqh2IQpH2jnrw09R30Ojm2HoJ1ChprHhtdYsDTnJy2siKe1Rik/v7ch/WtUyNr5wPdcMentIfwT0BWKBXUqpNVrryDy7DQCa2B9dgHn2X/8UCBwAKhZQ3UIUjphfYeXDkJYI/3kDujxitPtTcno2z6zcy7p9p+neqCrv3eFPrUoy4Sryx5Ez+s5AtNY6BkApFQQMAfIG/RDgS621BrYrpbyVUrW11vFKKR/gVuB1YGrBli9EAcnJgs2vw9Y5ULUx3LUEarczWsKOmHNMWRJGQmom0wc0Z3yvhjLhKgqEI0FfFziZ53ksfz1b/6d96gLxwPvAU8BVP7KnlBoPjAfw8/NzoCwhCsi5I9aEa1wodBxjncl7muuhmp1rY+5Ph/lwczT1qpQl+JHutPP1Nja+cH2OBP2VTim0I/sopQYBCVrr3UqpG682iNZ6PjAfICAg4PK/X4iCpzWEL4bvnwQ3D7jjS2g5xGgJJ86lE7gklNATSdze0YeXB7eiXGm5R0IULEe+omIB3zzPfYA4B/cZCQxWSg0EvICKSqmvtdb3XH/JQhSAS0nw/VTYHwz1esLwT6GSj9ESVobG8sKqCJSCD+9qz6C2dYyOL0oOR4J+F9BEKdUAOAWMBu66bJ81wET79fsuQLLWOh54xv7Afkb/pIS8KHIndliLkaWcgpufh55ToZS52xZTM7J5YdV+VoXF0al+ZWaP8sencllj44uS55pBr7XOUUpNBDZi3V65SGsdoZSaYN/+CbAO69bKaKzbKx8ovJKFuE65OfD7u/DrW1DJFx7cCL6djJaw58QFAoNCiUvKYGrfpjx6YyPc3aT7kyhcyrpRpngJCAjQISEhRV2GcCVJJ2HFQ3BiG7QdBQPfAS9zd/vm2jQfb47m/Z8OU7uSF3NG+9OxXhVj4wvXp5TarbUOuNI2mfURrm//Clj7uLW88LD50G6U0eFPJV1iSlAYO4+dZ4h/HV4d2pqKXh5GaxAlmwS9cF2ZF2HD0xD6NdQNgBGfQZWGRktYty+e6cF7ybVpZo9qx7D2Zid8hQAJeuGq4kJh+Vg4HwO9noQbp1u3UBqSlpnDjO8iWRJyEn9fb+aM9qdeVXP35guRlwS9cC02G2z7EH6aAeVrwJi1UL+n0RL2xSYTGBTK0XNpTLypMYF9muAhE66iCEnQC9eRetpapybmF2hxG9w2F8qam/C02TQLtsTw9sYoqpUvzeKHutK1obnuU0L8Ewl64RqiNljdn7LS4bY50OF+o92fzqRk8MTScLZEJzKgdS3eHN4G77LmGoULcTUS9MK5ZV+CH16EnfOhVhsYsRCqNzNawg+RZ3hqeTgZ2TbeGt6GUZ18UQZ/yAhxLRL0wnmdibQWI0uIhK6PQZ+XwL20seEzsnN5/fsDfLX9OK3qVGTune1pVL28sfGFcJQEvXA+WsOuBbDxOetDT3cHQ5M+Rks4EJ/C5MWhHE64yPjeDXmiX1NKu0v3J1E8SdAL55KWCKsnwqH10LgvDP3YurvGEK01//3jGG+uP0ilMh58NbYzvZpUNza+ENdDgl44jyObYeUEuHQe+s+ELg8bnXBNvJjJtGXhbI46yy3NazBrZFuqljd3qUiI6yVBL4q/nCz4+VX4Yy5Uawb3LLcmXg36JSqBJ5ftJSUjm1cGt+K+bvVkwlU4DQl6UbwlRlsTrvFhEPAg9HsdPM0t6ZuZk8usDVEs3HKUZjUr8PW4zjSvJa2PhXORoBfFk9YQ9g2sewrcPWHU19aHoAyKTkhl0uIwDsSncH+3ejwzsAVeHjLhKpyPBL0ofi4lWatNRqyE+r1g2KdQqa6x4bXWfLvzBK+ujaSspzsL7w/glhY1jY0vREGToBfFy/Ft1rrxKXFwy4vQ43Gj3Z8upGUxfcVeNkacoVeTarx7eztqVPQyNr4QhUGCXhQPuTnw29vw2yzw9oOxm8Dnij0UCs0f0YlMXRrOubRMnh3YnHE9G1KqlEy4CucnQS+KXtIJCH4ITm6HdnfCgFlGuz9l59p474dDfPLrERpUK8eC+3vQum4lY+MLUdgk6EXR2h8M300BNAxfAG1vNzr8scQ0AoNCCY9N5s7OvrwwqCVlPeXbQrgW+YoWRSPzIqx/yrqzxqcTjFgAlesbG15rTfCeU7y0ej/ubqWYd3cHBrSpbWx8IUySoBfmndoDwePgwlHoPQ1umA5u5r4Uky9l8/yq/XwXHkeXBlWYPcqfOt5ljI0vhGkS9MIcm836dOvPr0L5WnD/Wqjfw2gJu4+fZ/LiME6nZDDtP82YcEMj3GTCVbg4CXphRkq81f3p6K/QYjAMngtlKhsbPifXxkebjzDnp0P4VC7L8gndaO9nbnwhipIEvSh8B9fB6scgJ8Nq79fhPqOLkcVeSGfKkjB2HbvA8PZ1eWVIKyp4mWsULkRRk6AXhSf7Emx63lo7vlYbGLEIqjc1WsJ34XE8u3IfWsP7o/wZ2t7cJ2yFKC4cCnqlVH9gDuAGLNBav3XZdmXfPhBIB8ZorfcopbyA34DS9rGWa61fKsD6RXF1JgKWj4WzB6DbROtTrga7P13MzOHlNREs3x1Lez9v5oxqj19Vc4uhCVGcXDPolVJuwEdAXyAW2KWUWqO1jsyz2wCgif3RBZhn/zUTuFlrfVEp5QFsUUqt11pvL+DjEMWF1rDzM+tM3qsS3BMMjc12fwo/mURgUCgnzqcz+ebGTLqlCR5upYzWIERx4sgZfWcgWmsdA6CUCgKGAHmDfgjwpdZaA9uVUt5Kqdpa63jgon0fD/tDF1j1onhJS7SuxR/aAE36wZCPoby57ks2m2b+7zG8szGK6hVKs/ihrnRpWNXY+EIUV44EfV3gZJ7nsVhn69fapy4Qb/8fwW6gMfCR1nrHlQZRSo0HxgP4+fk5VLwoRo78bO/+dAH6vwVdJhidcD2dnMHUpWH8ceQcA9vU4s1hbalUViZchQDHgv5K362Xn5X/4z5a61zAXynlDaxUSrXWWu//285azwfmAwQEBMhZv7PIyYKfZ8AfH9i7P62AWq2NlrAx4jRPB+8lM9vGzBFtuCPAV7o/CZGHI0EfC/jmee4DxP3bfbTWSUqpX4D+wN+CXjihxGgIfhDiw4uk+9OlrFxe+z6Sb3acoHXdiswZ3Z5G1csbG18IZ+FI0O8CmiilGgCngNHAXZftswaYaL9+3wVI1lrHK6WqA9n2kC8D9AFmFlz5okhoDaFfw/qn7d2fvoEWg4yWcCA+hcmLQzmccJGHezfkiX7N8HSXCVchruSaQa+1zlFKTQQ2Yt1euUhrHaGUmmDf/gmwDuvWymis2ysfsL+8NvCF/Tp9KWCp1nptwR+GMOby7k/D50PFOsaG11rz+dZjvLX+IJXKevDV2M70amJuwlcIZ6SsG2WKl4CAAB0SElLUZYjL/dn9KTUebnoOegQa7f6UeDGTJ5eF80vUWW5pXoNZI9tStby5e/OFKM6UUru11lfs1iOfjBXX9pfuT/XgwU3g09FoCb9EJfDksnBSMnKYMaQV93atJxOuQjhIgl5c3eXdnwa+DaUrGBs+MyeXmeujWLT1KM1qVuCbcV1pVsvc+EK4Agl68c+KuPtTdEIqkxaHcSA+hfu61ePZgS3w8jB3qUgIVyFBL/4uM9W6o6YIuz8t3nmSGWsjKOvpzsL7A7ilRU1j4wvhaiToxV+d2gPBY+HCMej9FNzwtNHuTxfSspi+Yi8bI87Qq0k13r29HTUqehkbXwhXJEEvLH/p/lSzSLo/bTtyjilLwjiXlslzA1swtmcDSkn3JyHyTYJe/L37021zoGwVY8Nn59qY/cMh5v16hAZVy7Hg/h60rlvJ2PhCuDoJ+pIuaj2serTIuj8dS0wjcEkY4SeTGN3Jlxdva0lZT/myFKIgyXdUSZV9CTa9ALs+K5LuT1prVuw5xYur9+NWSvHx3R0Y2Ka2sfGFKEkk6EuiM5HWhGtCJHR9DPq8ZLT7U0pGNs+v3M+a8Dg6N6jC+6P8qeNdxtj4QpQ0EvQlidZW/9aNzxVZ96fdx88TGBRGfHIGT/ZryiM3NsZNJlyFKFQS9CVF2jl796f10LgvDJ1ntPtTTq6NjzYfYe7Ph6nj7cWyCd3o4FfZ2PhClGQS9CXBkc327k/ni6T706mkS0wJCmPnsfMM8a/Dq0NbU9FLuj8JYYoEvSvLyYLNr8HWuVCtCdy9DGq3NVrC2r1xPLNiHzab5r072jG8g4/R8YUQEvSu69wRWP4gxIdBxwfgP28Y7f6UlpnDK99FsDQklna+3swd7U+9quWMjS+E+H8S9K5Gawj7FtZNAzcPuOMraDnYaAn7YpMJDArl6Lk0HrupEY/3aYqHm3R/EqKoSNC7kktJ8P1Ua9XJej2t7k+V6hob3mbTLNgSw9sbo6harjTfjutKt0ZVjY0vhLgyCXpXcWK7tW58yim4+QXoOcVo96czKRk8sTScLdGJ9G9Vi7dGtMG7rKex8YUQ/0yC3tnl5sDv78Kvb4G3H4zdBD5X7CZWaH6MPMO05eFkZNt4c3gbRnfyle5PQhQjEvTOLOmk1cP1xDZoOwoGvgNeFY0Nn5Gdy+vfH+Cr7cdpWbsic+9sT+Ma5Y2NL4RwjAS9s4pYCWsCQdtg2HxoN8ro8AdPpzB5cSiHzlzkoV4NePI/zSjtLt2fhCiOJOidTVYarH8KQr+Guh2t7k9VGhobXmvNV9uP89r3B6jo5cEXD3bmhqbmPmErhPj3JOidSVyYtRjZuSPQcyrc9Kx1C6Uh5y5m8tTyvfx0MIEbm1XnndvbUa28ucXQhBDXR4LeGdhssP0j+PEVKFcd7l8DDXobLeH3w2eZujSc5PRsXrqtJWO615cJVyGchAR9cZd6BlZNgCM/Q/NBMPgDo92fsnJsvLMpivm/xdCkRnm+fLAzLWqbm/AVQuSfBH1xdmij1f0pKw0GzbaWMjB4Fh1z9iKTg0LZfyqFu7v48fytLSnjKROuQjgbhz6XrpTqr5SKUkpFK6WmX2G7UkrNtW/fq5TqYP9zX6XUZqXUAaVUhFIqsKAPwCVlZ8C6p+DbO6BCLRj/CwQ8aCzktdYs3XWSW+duIfbCJebf25HXh7WRkBfCSV3zjF4p5QZ8BPQFYoFdSqk1WuvIPLsNAJrYH12AefZfc4AntNZ7lFIVgN1KqR8ue63IK+GgNeF6Zj90fRRueQk8vIwNn5yezbMr9/H9vni6N6rKe3f4U6uSufGFEAXPkUs3nYForXUMgFIqCBgC5A3rIcCXWmsNbFdKeSulamut44F4AK11qlLqAFD3stcKsBYjC1kEG58Fz/Jw1zJo2s9oCTuPnufxoFASUjN5un9zHu7dkFLS/UkIp+dI0NcFTuZ5Hot1tn6tfepiD3kApVR9oD2w40qDKKXGA+MB/Pz8HCjLhaSfhzWT4OBaaHQzDP0EKtQ0NnxOro25Px3mw83R+FUpS/Aj3Wnn621sfCFE4XIk6K90Sqf/zT5KqfJAMPC41jrlSoNorecD8wECAgIu//td19HfYMV4SEuEfq9bl2tKmVvS9+T5dAKDQtlzIomRHX14eXArypeWOXohXIkj39GxgG+e5z5AnKP7KKU8sEL+G631iusv1cXkZsPmN2DLbKjaGO4Mgjr+RktYHXaK51fuB2DOaH+G+Jtb0lgIYY4jQb8LaKKUagCcAkYDd122zxpgov36fRcgWWsdr6xP1CwEDmit3yvAup3b+RgIHgendkOH+6w+rp7mui9dzMzhxdX7WbHnFB3rVeb9Uf74VjHXfUoIYdY1g15rnaOUmghsBNyARVrrCKXUBPv2T4B1wEAgGkgHHrC/vAdwL7BPKRVm/7NntdbrCvQonEn4Evj+CevyzO1fQKuhRocPO5lEYFCodcnmliZMurkx7tL9SQiXpqwbZYqXgIAAHRISUtRlFKyMFCvg9y0Fv+5W9ydv32u/roDk2jSf/HqE2T8comZFL94f7U+n+uY+YSuEKFxKqd1a6ys2o5BZNxNO7rLujU8+CTc+C72fNNr9KT75ElOXhLMt5hy3tq3NG8PaUKmMucXQhBBFS4K+MNlyYct7sPlNqFgXHlgPfl2NlrBh/2mmr9hLVo6NWSPbcntHH1mMTIgSRoK+sCSfgpUPw7HfodVwa62aMt7Ghr+Ulcur30fy7Y4TtKlbibl3tqdBNXMTvkKI4kOCvjAc+A5WT7RuoRzyMfjfZXQxsoi4ZCYvDuXI2TQe7t2QJ/o1w9NdJlyFKKkk6AtSVrq1hMHuz6FOexixEKo2Mja81ppFW48xc/1BvMt68PXYLvRsUs3Y+EKI4kmCvqCc3gfLx0LiIejxONz0HLh7Ghv+bGom05aH80vUWfq0qMmskW2pUs7c+EKI4kuCPr+0hh2fwA8vQpkqcO9KaHST0RJ+iUrgyWV7Sc3IZsaQVtzbtZ5MuAoh/keCPj8unoVVj0D0D9B0AAz5EMqZu1SSmZPLrA1RLNxylGY1K/DNuC40q1XB2PhCCOcgQX+9on+ElY9ARjIMeBs6P2R0wjU6IZVJi8M4EJ/C/d3q8czAFnh5SGMQIcTfSdD/WzmZ8NMM2PYhVG8B962Cmq2MDa+1ZvHOk8xYG0FZT3cW3h/ALS3MLWkshHA+EvT/RuJhWP4gnN4LncZBv9fAo4yx4ZPSs5gevI8NEafp2bga793RjhoVpfuTEOLqJOgdoTWEfgXrnwb30jD6W2h+q9EStsecY8qSMBIvZvLcwBaM7dlAuj8JIRwiQX8tly7Ad4EQuRrq97IWI6tYx9jw2bk23v/xEB//coQGVcux4pEetPGpZGx8IYTzk6C/muN/QPBDcPG01aS7R6DRxciOn0sjMCiMsJNJjArw5cXbWlJOuj8JIf4lSY0ryc2B32bBb2+Ddz14cBP4dDRawoo9sbywaj+lSik+vKs9g9qa+1+EEMK1SNBfLumEdRZ/cju0uxMGvg2lzd2bnpqRzQur9rMqLI7O9aswe7Q/db3NTfgKIVyPBH1e+4PhuymgbTB8AbS93ejwe05cIDAolLikDKb0acrEmxvjJhOuQoh8kqAHyLxo3VET9jX4dIIRC6ByfWPD59o0836JZvaPh6ldyYulD3elYz3p/iSEKBgS9HGh1mJk52Og15Nw43RwM9d9KS7pElOWhLHj6Hlua1eH14e1pqKXdH8SQhSckhv0Nhts+wB+ehXK14Axa6F+T6MlbNgfz9PB+8jJtfHu7e0Y3qGuLEYmhChwJTPoU09b3Z9ifoEWt8Ftc6GsuUsl6Vk5vLo2ksU7T9LWpxJzRkv3JyFE4Sl5QR+1AVY/ajUJGfQ+dBxjdDGy/aeSmRwUytHENCbc0IipfZtK9ychRKEqOUGfnQE/vmStHV+zDYxcCNWbGRveZtMs2nqUWRui/tf9qUdj6f4khCh8JSPoEw5YE64JEdDlEejzMniYWwzsbGomTy4L59dD0v1JCGGeawe91hCyEDY+Z33o6e7l0KSv0RI2RyUwbVk4qRk5vDq0Nfd08ZMJVyGEUa4b9OnnYc0kOLgWGt0MQz+BCubWbc/IzmXmhoN8vvUYzWtV4NuHutK0pnR/EkKY59AsoFKqv1IqSikVrZSafoXtSik11759r1KqQ55ti5RSCUqp/QVZ+FUd/Q3mdYdDG6014+8ONhry0QmpDPv4Dz7feowx3euz6rEeEvJCiCJzzTN6pZQb8BHQF4gFdiml1mitI/PsNgBoYn90AebZfwX4L/Ah8GXBlf0PcrNh8xuwZTZUbQx3LYHa7Qp92D9d3v1p0ZgAbm4u3Z+EEEXLkUs3nYForXUMgFIqCBgC5A36IcCXWmsNbFdKeSulamut47XWvyml6hd04X9z6QJ8PQJO7Yb298KAmeBp7t70XJtm4rd7WL//NL2aVOPd26X7kxCieHAk6OsCJ/M8j+X/z9avtk9dIN7RQpRS44HxAH5+fo6+7P95eUOVhtB9ErQa9u9fn09upRQNqpXj2YHNGdezoXR/EkIUG44E/ZUSS1/HPleltZ4PzAcICAj4V6+1KlDWYmRF6Kn+zYt0fCGEuBJHJmNjAd88z32AuOvYRwghRBFwJOh3AU2UUg2UUp7AaGDNZfusAe6z333TFUjWWjt82UYIIUThuWbQa61zgInARuAAsFRrHaGUmqCUmmDfbR0QA0QDnwGP/vl6pdRiYBvQTCkVq5QaW8DHIIQQ4iqUdaNM8RIQEKBDQkKKugwhhHAaSqndWuuAK22TZROFEMLFSdALIYSLk6AXQggXJ0EvhBAurlhOxiqlzgLHr/Pl1YDEAizH2cjxl+zjB3kPSurx19NaV7/ShmIZ9PmhlAr5p5nnkkCOv2QfP8h7UNKP/0rk0o0QQrg4CXohhHBxrhj084u6gCImxy9K+ntQ0o//b1zuGr0QQoi/csUzeiGEEHlI0AshhItzmqDPZ4Pyq77WWVzve6CU8lVKbVZKHVBKRSilAs1Xn3/5+Rqwb3dTSoUqpdaaq7rg5PN7wFsptVwpddD+ddDNbPUFI5/vwRT71/9+pdRipVTJ6fWptS72D8ANOAI0BDyBcKDlZfsMBNZjdbvqCuxw9LXO8Mjne1Ab6GD/fQXgkLO9B/k5/jzbpwLfAmuL+nhMHz/wBTDO/ntPwLuoj8nke4DV2vQoUMb+fCkwpqiPydTDWc7o/9egXGudBfzZoDyv/zUo11pvB7yVUrUdfK0zuO73QFtN2vcAaK1TsfoK1DVZfAHIz9cASikf4FagaPtNXr/rPn6lVEWgN7AQQGudpbVOMlh7QcnX1wBW69QySil3oCwlqAueswT9PzUfd2QfR17rDPLzHvyPUqo+0B7YUfAlFqr8Hv/7wFOArZDqK2z5Of6GwFngc/ulqwVKqXKFWWwhue73QGt9CngHOAHEY3XB21SItRYrzhL0+WlQnu/G5cVEvpu0K6XKA8HA41rrlAKszYTrPn6l1CAgQWu9u+DLMiY///7uQAdgnta6PZAGOONcVX6+Bipjne03AOoA5ZRS9xRwfcWWswR9fhqUu0rj8nw1aVdKeWCF/Dda6xWFWGdhyc/x9wAGK6WOYf13/2al1NeFV2qhyO/3QKzW+s//xS3HCn5nk5/3oA9wVGt9VmudDawAuhdircVLUU8SOPLAOiOJwfpp/OckTKvL9rmVv07C7HT0tc7wyOd7oIAvgfeL+jiK4vgv2+dGnHMyNl/HD/wONLP//mXg7aI+JpPvAdAFiMC6Nq+wJqcnFfUxmXq4O/LDoKhprXOUUn82KHcDFml7g3L79k+wGpQPxGpQng48cLXXFsFh5Et+3gOsM9p7gX1KqTD7nz2rtV5n8BDyJZ/H7/QK4PgnAd8opTyxwtLp3pt85sAOpdRyYA+QA4RSgpZKkCUQhBDCxTnLNXohhBDXSYJeCCFcnAS9EEK4OAl6IYRwcRL0Qgjh4iTohRDCxUnQCyGEi/s/XTR6hn8HeFkAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig = plt.figure()\n", "ax = fig.add_subplot(1, 1, 1)\n", "data1, data_t = simulate(50, state0=[0.01,0.01], h= 0.02)\n", "ax.plot(data1[:, 0], data1[:, 1])\n", "data2, data_t = simulate(50, state0=[0.0,0.01], h= 0.02)\n", "ax.plot(data2[:, 0], data2[:, 1])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### No way thes trajectories leave (0, 0) tangential to the y-axis. In fact, they are more aligned with the fast direction, i.e. the x-axis. Trajectories are just \"more tangential\" to the y-axis at the start than later on, since with increasing time trajectories get more and more aligned with the fast direction." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### When time is reversed, trajectories follow the same trajectories for all dynamical systems, just the direction is reversed. Reversibility is not required for that. Accordingly, the same thing applies for stable nodes: trajectories are \"more tangential\" to the slow direction when approaching a stable node than they were when further away from the node, but for example when time is reversed for this sheep-rabbit model, trajectories approaching (0, 0) are not tangential to the y-axis, but in fact more aligned with the x-axis - just less so than far awau from (0, 0)." ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "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.9" } }, "nbformat": 4, "nbformat_minor": 4 }