Skip to content

Instantly share code, notes, and snippets.

@mimakaev
Created December 4, 2020 03:13
Show Gist options
  • Save mimakaev/0327bf6ffe7057ee0e0625092ec8e318 to your computer and use it in GitHub Desktop.
Save mimakaev/0327bf6ffe7057ee0e0625092ec8e318 to your computer and use it in GitHub Desktop.
polychrom_repulsive_potential
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"metadata": {},
"cell_type": "markdown",
"source": "### Polychrom polymer repulsive potential \n\n$$ U(x) = U_0 \\left(1 + \\left(\\frac{r \\sqrt{6/7}}{R_{max}}\\right) ^{12} \\left(\\left(\\frac{r \\sqrt{6/7}}{R_{max}}\\right)^2 - 1\\right) \\right)$$ \n\nIn the example below $R_{max} = 1.05$ and $U_0 = 4$"
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "import numpy as np \nimport matplotlib.pyplot as plt \n%matplotlib inline \n\ndef f(x):\n term1 = ( x* np.sqrt(6/7)) / (1.05 )\n term2 = 1 + (term1**12) * (term1**2 - 1) * (823543/46656)\n return 4 * term2\n\nx = np.linspace(0,1.05,500)\nplt.figure()\nplt.plot(x,f(x))\nplt.xlabel(\"radius\")\nplt.ylabel(\"Value of potential\")\n\nplt.figure()\nplt.plot(x[:-1],np.diff(f(x)))\nplt.xlabel(\"radius\")\nplt.ylabel(\"Derivative of potential\")\n",
"execution_count": 2,
"outputs": [
{
"output_type": "execute_result",
"execution_count": 2,
"data": {
"text/plain": "Text(0, 0.5, 'Derivative of potential')"
},
"metadata": {}
},
{
"output_type": "display_data",
"data": {
"text/plain": "<Figure size 432x288 with 1 Axes>",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEGCAYAAABo25JHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAAAlD0lEQVR4nO3deZgV5Zn38e/dO0uzCE2zNi3KjoDYIuAukiBGTNQ4kHGUxISJSSZ7ZuLMXDFxJu+8Jpn4xpiJwSVq4jrqRNxiFDW4gTSrsknLvnaz79DLPX+c0rdtm+7TTdeps/w+11VXn1P1nDq/Yrupp6qex9wdERHJXFlRBxARkWipEIiIZDgVAhGRDKdCICKS4VQIREQyXE7UAVqqe/fuXlpaGnUMEZGUsnDhwp3uXtTYtpQrBKWlpZSXl0cdQ0QkpZjZhhNtU9eQiEiGUyEQEclwKgQiIhlOhUBEJMOpEIiIZLjQC4GZZZvZYjN7tpFt+Wb2mJlVmNl8MysNO4+IiHxcIs4IvgWsPMG2G4E97n46cDtwWwLyiIhIPaE+R2BmfYHLgZ8C322kyZXAj4PXTwB3mpl5CGNjr95+gOeWbW0s5CdXNfL5RpphjbRsrF2L9nmiHbRhnjCO75PtTubXNb79teTzrf19NoycLCMn28jOMnKysoKfRna2kfvh+4+2x9p8+L4gN5uCnCza5WVTkJNNVlacv4AiCRT2A2X/D/hHoPAE2/sAmwDcvcbM9gHdgJ31G5nZTGAmQElJSauCVFQe5NevVnxsnaZikETLy8miICeLgtzsj4pDQV42HfOz6VSQS+d2uXRqF/wsyKFT8L5bhzx6FBbQvWMeOdm6tCdtK7RCYGafASrdfaGZXXQy+3L3WcAsgLKyslb98335yF5cPvLy1n5/I+saadeSzzfarrF9xvfdjX9v6/fXeL7WH0djDU8my8nmaey7G1tV51DrTk1tHTV1Tm2dU1Mb/Kyro7bOqW7w/sN21bV1HKuu42hNLUeO13K0uo4j1bUcrbccqa7lSHUdh47VULn/IPuOVLPvSDXHauoaPWYz6NYhj6LCAnoU5lPcKZ++XdvTv1t7Sk5pT/9uHejaPjfuM0sRCPeM4FxgqplNAQqATmb2R3e/rl6bLUA/YLOZ5QCdgV0hZmqVRrs5WvT3TH8ppWWOVtey/2g1+4/UsO9INbsOHqPyQGypOnCUyv2x1yu27afqwLGPfbYwP4dTizowpGchQ3p2Ykiv2M9TOuRFdDSS7EIrBO5+M3AzQHBG8P0GRQBgNnAD8DZwDfBKGNcHRFJNQW42BbnZ9DhRp2o9R47XsmnPYTbsOsyGXYfYuPswH1Qd5OWVlTxevvmjdn26tGNM/66MKenCWf27MrRXJ3LVzSREMOicmd0KlLv7bOBe4A9mVgHsBqYlOo9IqmuXl82g4kIGFX+8arg7VQePsXr7AVZu28/SzftYsG43zyyN3TTRIS+b8ad144JBRVwwsIjS7h2iiC9JwFLtP+BlZWWu0UdFWm/r3iMs2riHtz/YxetrdrJx92EA+ndrz+ThPZlyRi9G9u2s6wxpxswWuntZo9tUCEQy2/qdh5i7popXVlXyxpqd1NQ5fbu244pRvfmbsn46U0gTKgQiEpd9h6v5y4rtPLtsG29U7KS2zhk/oBvTzynh08OLyc/JjjqitJIKgYi02I79R/nv8k08umATm/ccobhTPjeedyrTx5ZQWJAbdTxpIRUCEWm1ujpn7poqZs1dy1sf7KKwIIfrx/fnK+cPoEt73ZKaKlQIRKRNLN20l9/N/YAX3ttOYX4OX7v4dGZMKKUgV11GyU6FQETa1Krt+/nZn1fzyqpKenUu4J+nDOUzI3vpTqMk1lQh0NMkItJiQ3p24r4ZZ/PozHF065jHPzyymOvve4d1Ow9FHU1aQYVARFpt3IBuPP318/jJ1OEs2biXT98+l9/99QNq61KrpyHTqRCIyEnJzjJumFDKnO9dyMVDiviPF1Yx/e55bAoeVJPkp0IgIm2iR6cC7rruLH7x+VGs2Lqfy371Os8t2xZ1LImDCoGItBkz45qz+vLCt85nUHFHvv7wIn763ApqahsfVluSgwqBiLS5fqe059GZ47l+fH/ufn0d1907n10HjzX/QYmECoGIhCIvJ4tbrxzBL68dxeKNe7nqt2+xXncVJSUVAhEJ1VVj+vLIzHHsP1LNVb99i8Ub90QdSRpQIRCR0I0p6cqTN02gY34O0++ex6urK6OOJPWoEIhIQgwo6shTX5vAaUUd+fsHF/LKqh1RR5JAaIXAzArM7B0zW2pmy83sJ420mWFmVWa2JFi+HFYeEYle9475PPTlcxjcs5C//8NCXl6hYpAMwjwjOAZc4u6jgNHAZDMb10i7x9x9dLDcE2IeEUkCXdrn8ccbz2For07c9NBC5qxUMYhaaIXAYw4Gb3ODRc+diwid2+fyhxvPYUjPTnztoUUs3LA76kgZLdRrBGaWbWZLgErgJXef30izq81smZk9YWb9TrCfmWZWbmblVVVVYUYWkQTp3C6X33/xbHp1LuBL95ezZseBqCNlrFALgbvXuvtooC8w1sxGNGjyDFDq7iOBl4AHTrCfWe5e5u5lRUVFYUYWkQTq3jGfB790DrnZWVx/3zts23ck6kgZKSF3Dbn7XuBVYHKD9bvc/cPHDe8BzkpEHhFJHiXd2vPAl87mwNEavvxAOUeO10YdKeOEeddQkZl1CV63AyYBqxq06VXv7VRgZVh5RCR5De/dmTumj2bFtv3805PLSLUJs1JdmGcEvYBXzWwZsIDYNYJnzexWM5satPlmcGvpUuCbwIwQ84hIErtkSDHf/9RgZi/dyu/mro06TkbRVJUikjTcnW88vJjn39vG72eczUWDe0QdKW1oqkoRSQlmxs8/P5LBxYV89/Gl7Nh/NOpIGUGFQESSSvu8HO78whiOHK/l248u0bSXCaBCICJJ5/QeHfnJ1OG8vXYXv32tIuo4aU+FQESS0ufL+nLFqN7c/vIaPXkcMhUCEUlKZsZPPzeC3l0K+O7jSzl8vCbqSGlLhUBEklanglx+fs0oNuw6zM9fXB11nLSlQiAiSW3cgG5cP74/97+1nvL16iIKgwqBiCS9f5o8hD5d2vGPTyzjaLWGoGhrKgQikvQ65Odw29UjWbvzELe/9H7UcdKOCoGIpIRzT+/OtLP7cc8b63hfQ1a3KRUCEUkZ/zh5CIUFOfzo6fc0MF0bUiEQkZRxSoc8fvDpwcxbu5vZS7dGHSdtqBCISEqZdnYJZ/TpzE+fW8mBo9VRx0kLKgQiklKys4x/++wIqg4e4445a6KOkxZUCEQk5Yzu14XPn9WX+99az8Zdh6OOk/JUCEQkJX3vU4PJzjJ+/hc9cXyywpyqssDM3jGzpcEsZD9ppE2+mT1mZhVmNt/MSsPKIyLppbhTAV85fwDPLN3K0k17o46T0sI8IzgGXOLuo4DRwGQzG9egzY3AHnc/HbgduC3EPCKSZmZeMIBuHfL4P8+v1O2kJyG0QuAxB4O3ucHS8HfqSuCB4PUTwEQzs7AyiUh6KSzI5duXDmT+ut3MWVkZdZyUFeo1AjPLNrMlQCWxyevnN2jSB9gE4O41wD6gWyP7mWlm5WZWXlVVFWZkEUkx08aWcGr3Dtz251XUaTazVgm1ELh7rbuPBvoCY81sRCv3M8vdy9y9rKioqE0zikhqy83O4juTBrGm8iDPvbst6jgpKSF3Dbn7XuBVYHKDTVuAfgBmlgN0BnYlIpOIpI/Lz+jF6T06csecNZrjuBXCvGuoyMy6BK/bAZOAVQ2azQZuCF5fA7ziuuIjIi2UnWV8a+JAnRW0UphnBL2AV81sGbCA2DWCZ83sVjObGrS5F+hmZhXAd4EfhphHRNLYlDN6MVBnBa2SE9aO3X0ZcGYj639U7/VR4PNhZRCRzJGdZXzr0oF84+HFPPfuNqaO6h11pJShJ4tFJG1MGdGLQcWxswLdQRQ/FQIRSRtZWcbXLz6disqDzFml5wridcJCYGYHzGx/I8sBM9ufyJAiIvG6/Ixe9OnSjt/99YOoo6SMExYCdy90906NLIXu3imRIUVE4pWTncVXzj+V8g17KF+/O+o4KSHuriEz62FmJR8uYYYSETkZ157dj67tc7nrr2ujjpISmi0EZjbVzNYA64C/AuuBF0LOJSLSau3zcrh+fCkvr9xBRaUmum9OPGcE/waMA95391OBicC8UFOJiJykGyaUUpCbxe90VtCseApBtbvvArLMLMvdXwXKQs4lInJSTumQx7Vl/Xh6yVaqDhyLOk5Si6cQ7DWzjsBc4CEz+xVwKNxYIiInb8aEUo7X1vHIOxujjpLU4ikEVwJHgO8AfwY+AK4IM5SISFsYUNSRCwcV8cd5G6iurYs6TtJqthC4+6FgOOkad3/A3e8IuopERJLejAmlVB44xgvvbY86StJq6oGyN4KfDR8s0wNlIpIyLhxURGm39tz/5rqooyStph4oOy/42fDBMj1QJiIpIyvLuH58KYs27mXZ5r1Rx0lK8TxH8Id41omIJKtryvrSIS+b+99aH3WUpBTPxeLh9d8EM4mdFU4cEZG216kgl6vP6suzS7ex+9DxqOMknaauEdxsZgeAkfWvDwA7gKeb27GZ9TOzV81shZktN7NvNdLmIjPbZ2ZLguVHje1LRORkfeGcEo7X1vHUos1RR0k6TV0j+A93LwR+3uD6QDd3vzmOfdcA33P3YcSeTP66mQ1rpN3r7j46WG5t3WGIiDRtSM9OjCnpwsPvbEQz4n5cPLeP3mxmfcxsgpld8OESx+e2ufui4PUBYCXQ5+Qji4i0zvSxJaytOsQ76zQqaX3xXCz+v8CbwL8CPwiW77fkS8yslNi0lfMb2TzezJaa2QtmNryR7ZjZTDMrN7Pyqqqqlny1iMhHPjOyN4UFOXrSuIF45iz+HDDY3Vs1WEcwPMWTwLfdveHzB4uA/u5+0MymAH8CBjbch7vPAmYBlJWV6ZxORFqlXV42nzuzD48u2MQth47TtUNe1JGSQjx3Da0FcluzczPLJVYEHnL3pxpud/f97n4weP08kGtm3VvzXSIi8Zh2dgnHa+p4avGWqKMkjXjOCA4DS8xsDvDRWYG7f7OpD5mZAfcCK939lydo0xPY4e5uZmOJFSYNXyEioRnWuxOj+3XhkXc28qVzS4n9U5XZ4ikEs4Olpc4F/g5418yWBOv+GSgBcPe7gGuAm8yshtjAdtNcl/NFJGTTzu7HD596lyWb9nJmSdeo40Su2ULg7g+YWTugxN1Xx7tjd38DaLLUuvudwJ3x7lNEpC1MGdmLW2Yv58lFm1UIiO+uoSuAJcSGoMbMRptZa84QRESSQqeCXD49vCfPLN3GsZraqONELp6LxT8GxgJ7Adx9CTAgtEQiIglw9Vl92XekmjkrK6OOErl4p6rc12CdZngQkZR23undKe6UryEniK8QLDezLwDZZjbQzH4NvBVyLhGRUGVnGZ89sw+vra5i58HMntM4nkLwD8RGID0GPAzsAz4xgJyISKq5Zkxfauqcp5dsjTpKpOIpBJe7+7+4+9nB8q/A1LCDiYiEbWBxISP7dubJhZndPRRPIWhspNF4Rh8VEUl6V4/py4pt+1m5LXNn4G1qPoLLgusBfczsjnrL/cSGmBYRSXlTR/UmN9sy+qygqTOCrUA5cBRYWG+ZDXw6/GgiIuHr2iGPiwb34Nll26iry8yBDZqamGapuz8AnA48wv8vBM+4+54E5RMRCd0Vo3qzff9RFqzPzHkK4rlGMAFYA/wG+C/g/XgmphERSRWXDu1Bu9xsZi/NzLuH4ikEvwQ+5e4XuvsFxLqFbg83lohI4rTPy+HSYcW88N52qmsz73nZeApBbv3B5tz9fVo5P4GISLKaOqo3uw8d582KnVFHSbh4CkG5md1jZhcFy93ELiKLiKSNCwZ1p7Agh2eWbos6SsLFUwhuAlYA3wyWFcE6EZG0kZ+TzeThPfnL8u0crc6sEUmbLQTBXMV3Aj8BbgF+09r5i0VEktnU0b05cKyG11ZXRR0loeKZj+By4APgV8QKQoWZXRbH5/qZ2atmtsLMlpvZJ8Ynspg7zKzCzJaZ2ZjWHISISFsYP6Ab3Tvm8UyG3T0Uz1SV/wlc7O4VAGZ2GvAc8EIzn6sBvufui8ysEFhoZi+5+4p6bS4DBgbLOcBvg58iIgmXk53FlDN68Xj5Jg4eq6Fjfjz/RKa+eK4RHPiwCATWAgea+5C7b3P3RcHrA8BKoE+DZlcCD3rMPKCLmfWKL7qISNu7YlRvjlbXMWfljqijJEy8dw09b2YzzOwG4BlggZldZWZXxfMlZlYKnAnMb7CpD7Cp3vvNfLJYYGYzzazczMqrqjKr705EEuuskq70KMznhXe3Rx0lYeIpBAXADuBC4CKgCmgHXAF8prkPm1lH4Eng2+7equH93H2Wu5e5e1lRUVFrdiEiEpesLGPyiJ689n4lh49nxviazXaAufsXW7tzM8slVgQecvenGmmyBehX733fYJ2ISGQmj+jJg29v4K+rq7jsjPTvrY7njKBVzMyAe4GV7v7LEzSbDVwf3D00Dtjn7pn3NIeIJJWxpadwSoc8XngvM7qHwrwkfi7wd8C7ZrYkWPfPQAmAu98FPA9MASqAw0Crzz5ERNpKTnYWnxpWzLPLtnG0upaC3OyoI4XqhIXAzL7l7r8ys3Pd/c2W7tjd3wCsmTYOfL2l+xYRCdvkET15dMEm3lizk0uHFUcdJ1RNdQ19+L/zXyciiIhIMplwWmzsoUzoHmqqa2ilma0BepvZsnrrjdh/5keGG01EJDp5OVlMGlrMyyt3UF1bR252aJdUI3fCQuDu082sJ/AiMDVxkUREksPkET15avEW3v5gFxcMSt9b15ssce6+3d1HAduAwmDZ6u4bEhFORCRKFwwqokNedtp3D8Uz6NyFaKpKEclABbnZXDykBy+t2E5tGk9sr6kqRUSacNmIXuw8eDytJ7bXVJUiIk24aHAR+TlZ/DmNu4c0VaWISBM65Odw/sAiXlqxg9ijT+lHU1WKiDRj0rAebNl7hBXbWjVuZtKLZ9C5Y8SuE5xovCARkbR2yZBizN7lpRU7GN67c9Rx2lz6PiEhItJGigrzGVPSlZdWpOdkNSoEIiJxmDSsmOVb97Nl75Goo7S5uAuBmbUPM4iISDKbFAw8l45TWMbzQNkEM1sBrArejzKz/wo9mYhIEjmtqCMDijqkZfdQPGcEtxN7iGwXgLsvBfRksYhknEnDipm3dhf7j1ZHHaVNxdU15O6bGqyqDSGLiEhSmzS0mOpa57XVVVFHaVPxFIJNZjYBcDPLNbPvAyub+5CZ3WdmlWb23gm2X2Rm+8xsSbD8qIXZRUQS6sySrnTrkJd23UPxFIKvEptFrA+xieVHE9+sYvcDk5tp87q7jw6WW+PYp4hIZLKzjIlDe/DaqkqO19RFHafNNFsI3H2nu/+tuxe7ew93v87dd8XxublA+o7SJCIZadKwnhw4VsP8dc3+M5gymn2y2Mx+D3xigA13/1IbfP94M1sKbAW+7+7LT5BhJjAToKSkpA2+VkSkdc47vTsFuVm8vGIH5w9Mj8lq4ukaehZ4LljmAJ2Ag23w3YuA/sHEN78G/nSihu4+y93L3L2sqCg9fuFFJDW1y8tOu0Ho4ukaerLe8hBwLVB2sl/s7vvd/WDw+nkg18y6n+x+RUTCNmlYMVv3HWX51vQYhK41Q0wMBHqc7BebWU8zs+D12CBL+nS6iUjamjikB2akzd1D8VwjOEDsGoEFP7cD/xTH5x4BLgK6m9lm4BaCCW3c/S7gGuAmM6sBjgDTPF3Os0QkrXXrmM9ZJV2Zs2oH35k0KOo4Jy2eYagLW7Njd5/ezPY7gTtbs28RkahNHFrMbX9exbZ9R+jVuV3UcU7KCbuGzGxMU0siQ4qIJJtJw2I95C+vrIw4yclr6ozgP5vY5sAlbZxFRCRlnFbUkdJu7Zmzcgd/N65/1HFOygkLgbtfnMggIiKpxMyYOLSYP7y9gUPHauiQ32xPe9KK664hMxthZtea2fUfLmEHExFJdpcOLeZ4bR2vr0ntQejimY/gFmIPfP0auBj4GTA15FwiIkmvrLQrndvlpvx1gnjOCK4BJgLb3f2LwCgg/WZvFhFpodzsLC4eXMQrqyqprUvdu9/jKQRH3L0OqDGzTkAl0C/cWCIiqWHi0GJ2HzrO4o17oo7SavEUgnIz6wLcDSwkNkbQ22GGEhFJFRcOLiIny3gphecybuo5gt+Y2bnu/jV33xs8DTwJuCHoIhIRyXidCnIZN6Abc1L4OkFTZwTvA78ws/Vm9jMzO9Pd17v7skSFExFJBROH9qCi8iDrdh6KOkqrnLAQuPuv3H08cCGxweDuM7NVZnaLmaX+4BoiIm3k0qHFAMxJ0e6heIah3uDut7n7mcB04LPEMWexiEim6HdKe4b0LOTldC0EZpZjZleY2UPAC8Bq4KrQk4mIpJBLhxazYP0e9h4+HnWUFmvqYvEkM7sP2Ax8hdgMZae5+zR3fzpRAUVEUsHEoT2orXNeW516Txk3dUZwM/AWMNTdp7r7w+6emldCRERCNqpvF7p3zE/J20ibulh8ibvf4+6tekrCzO4zs0oze+8E283M7jCzCjNbpqGtRSSVZWUZlw7twdzVVRyvqYs6Tou0ZqrKeN0PTG5i+2XEpr0cCMwEfhtiFhGR0E0cWsyBYzW8s2531FFaJLRC4O5zgaZ+Na4EHvSYeUAXM+sVVh4RkbCdd3p38nOyUu7uoTDPCJrTB9hU7/3mYN0nmNlMMys3s/KqqtS7ECMimaFdXjbnD+zOyyt3kEpTsEdZCOLm7rPcvczdy4qKiqKOIyJyQpcOLWbzniOs3nEg6ihxi7IQbOHjo5j2DdaJiKSsS4YEcxmvSJ3uoSgLwWzg+uDuoXHAPnffFmEeEZGT1qNTAaP6deGlFBqELrRCYGaPEBuuerCZbTazG83sq2b21aDJ88BaoILYENdfCyuLiEgiTRrag6Wb9lJ54GjUUeIS2mzL7j69me0OfD2s7xcRicrEocX84i/v88rKSqaNLYk6TrNS4mKxiEgqGdKzkD5d2qXMbaQqBCIibczMmDSsmDcqdnLkeG3UcZqlQiAiEoJLhxZztLqONyt2Rh2lWSoEIiIhGHvqKRTm56RE95AKgYhICPJysrhgcBEvr6ykri65nzJWIRARCcmkocXsPHiMZVv2RR2lSSoEIiIhuWhwEdlZlvRPGasQiIiEpEv7PMr6d0366wQqBCIiIZo0rJhV2w+waffhqKOckAqBiEiIJg4tBmBOEp8VqBCIiITo1O4dOK2oAy8n8SB0KgQiIiG7dFgx89buYt+R6qijNEqFQEQkZJOH96SmzpP27iEVAhGRkI3u14U+Xdrx3LvJOeWKCoGISMjMjCln9OT1NVVJ2T2kQiAikgBTzuhFdW1ydg+FWgjMbLKZrTazCjP7YSPbZ5hZlZktCZYvh5lHRCQqydw9FOZUldnAb4DLgGHAdDMb1kjTx9x9dLDcE1YeEZEoJXP3UJhnBGOBCndf6+7HgUeBK0P8PhGRpJas3UNhFoI+wKZ67zcH6xq62syWmdkTZtavsR2Z2UwzKzez8qqqqjCyioiELlm7h6K+WPwMUOruI4GXgAcaa+Tus9y9zN3LioqKEhpQRKStfKx76HDydA+FWQi2APX/h983WPcRd9/l7seCt/cAZ4WYR0QkclNH9aG61pPqrCDMQrAAGGhmp5pZHjANmF2/gZn1qvd2KrAyxDwiIpEb0acTA3t05KlFm6OO8pHQCoG71wDfAF4k9g/84+6+3MxuNbOpQbNvmtlyM1sKfBOYEVYeEZFkYGZ8bkwfyjfsYcOuQ1HHAUK+RuDuz7v7IHc/zd1/Gqz7kbvPDl7f7O7D3X2Uu1/s7qvCzCMikgw+O7oPZvDUoi3NN06AqC8Wi4hknN5d2jHhtG78z+ItuEc/sb0KgYhIBD53Zl827j7Mwg17oo6iQiAiEoXJI3rSLjebx8s3Nd84ZCoEIiIR6Jifw9RRvXlm6bbIh5xQIRARich14/pzpLqW/4n4VlIVAhGRiJzRtzOj+nbmofkbI71orEIgIhKhvx3XnzWVB5m/bndkGVQIREQidMXI3nRtn8s9r6+NLIMKgYhIhNrlZXPDhFJeXlnJ+zsORJJBhUBEJGI3jC+lXW42d/31g0i+X4VARCRiXTvk8Tdn92P2kq1s3nM44d+vQiAikgRmXjCA7CzjFy+uTvh3qxCIiCSB3l3aceN5p/KnJVtZtnlvQr9bhUBEJEncdNFpdO+Yx78/tzKhzxWoEIiIJInCgly+96nBvLNuNw+/szFh36tCICKSRP6mrB/nD+zOvz+7korKxNxOGmohMLPJZrbazCrM7IeNbM83s8eC7fPNrDTMPCIiyS4ry/jZNSPpkJ/DjN8voHL/0fC/M6wdm1k28BvgMmAYMN3MhjVodiOwx91PB24Hbgsrj4hIqujVuR33zShj96HjTL3zTd5YszPUawZhnhGMBSrcfa27HwceBa5s0OZK4IHg9RPARDOzEDOJiKSEkX278N9fHU9OtnHdvfMZ/x+vhDYMRU4oe43pA9SfcWEzcM6J2rh7jZntA7oBO+s3MrOZwEyAkpKSsPKKiCSV4b0789J3LuSZpVt564OdFBXmh/I9YRaCNuPus4BZAGVlZdFP8CkikiDt8rK59ux+XHt2v9C+I8yuoS1A/eR9g3WNtjGzHKAzsCvETCIi0kCYhWABMNDMTjWzPGAaMLtBm9nADcHra4BXPMrZGUREMlBoXUNBn/83gBeBbOA+d19uZrcC5e4+G7gX+IOZVQC7iRULERFJoFCvEbj788DzDdb9qN7ro8Dnw8wgIiJN05PFIiIZToVARCTDqRCIiGQ4FQIRkQxnqXa3pplVARta+fHuNHhqOU3pONOLjjO9RHWc/d29qLENKVcIToaZlbt7WdQ5wqbjTC86zvSSjMepriERkQynQiAikuEyrRDMijpAgug404uOM70k3XFm1DUCERH5pEw7IxARkQZUCEREMlxaFgIzm2xmq82swsx+2Mj2fDN7LNg+38xKI4h50uI4zu+a2QozW2Zmc8ysfxQ5T1Zzx1mv3dVm5maWVLfmxSue4zSza4Pf0+Vm9nCiM7aFOP7clpjZq2a2OPizOyWKnCfDzO4zs0oze+8E283M7gh+DZaZ2ZhEZ/wYd0+rhdiQ1x8AA4A8YCkwrEGbrwF3Ba+nAY9FnTuk47wYaB+8vildjzNoVwjMBeYBZVHnDun3cyCwGOgavO8Rde6QjnMWcFPwehiwPurcrTjOC4AxwHsn2D4FeAEwYBwwP8q86XhGMBaocPe17n4ceBS4skGbK4EHgtdPABPNzBKYsS00e5zu/qq7Hw7eziM2S1yqief3E+DfgNuAo4kM14biOc6vAL9x9z0A7l6Z4IxtIZ7jdKBT8LozsDWB+dqEu88lNsfKiVwJPOgx84AuZtYrMek+KR0LQR9gU733m4N1jbZx9xpgH9AtIenaTjzHWd+NxP4HkmqaPc7gtLqfuz+XyGBtLJ7fz0HAIDN708zmmdnkhKVrO/Ec54+B68xsM7H5TP4hMdESqqV/f0OVEpPXy8kxs+uAMuDCqLO0NTPLAn4JzIg4SiLkEOseuojY2d1cMzvD3fdGGSoE04H73f0/zWw8sVkMR7h7XdTB0lU6nhFsAfrVe983WNdoGzPLIXb6uSsh6dpOPMeJmV0K/Asw1d2PJShbW2ruOAuBEcBrZraeWH/r7BS8YBzP7+dmYLa7V7v7OuB9YoUhlcRznDcCjwO4+9tAAbGB2tJJXH9/EyUdC8ECYKCZnWpmecQuBs9u0GY2cEPw+hrgFQ+u4KSQZo/TzM4EfkesCKRifzI0c5zuvs/du7t7qbuXErsWMtXdy6OJ22rx/Ln9E7GzAcysO7GuorUJzNgW4jnOjcBEADMbSqwQVCU0ZfhmA9cHdw+NA/a5+7aowqRd15C715jZN4AXid2hcJ+7LzezW4Fyd58N3EvsdLOC2AWdadElbp04j/PnQEfgv4Nr4RvdfWpkoVshzuNMeXEe54vAp8xsBVAL/MDdU+pMNs7j/B5wt5l9h9iF4xmp9h81M3uEWNHuHlzruAXIBXD3u4hd+5gCVACHgS9GkzRGQ0yIiGS4dOwaEhGRFlAhEBHJcCoEIiIZToVARCTDqRCIiGQ4FQKRNmBmF5nZs8HrqU2NkiqSbNLuOQKRthQMRmgtGd4guBc+LZ5vkMygMwKRBsysNBgv/0HgPeBeMysP5gD4Sb12k81slZktAq6qt36Gmd0ZvL7fzK6pt+1g8LOXmc01syVm9p6ZnZ+wAxRpQGcEIo0bCNzg7vPM7BR3321m2cAcMxtJbJyfu4FLiD0d+lgL9/8F4EV3/2mw3/ZtGV6kJXRGINK4DcE48QDXBv/rXwwMJzZZyhBgnbuvCYY/+GML978A+KKZ/Rg4w90PtFFukRZTIRBp3CEAMzsV+D4w0d1HAs8RGwQtXjUEf8+CIbPz4KOJSy4gNuLk/WZ2fdtFF2kZFQKRpnUiVhT2mVkxcFmwfhVQamanBe+nn+Dz64GzgtdTCQYes9j80Tvc/W7gHmLTGopEQtcIRJrg7kvNbDGxf/g3AW8G64+a2UzgOTM7DLxObG6Ehu4GnjazpcCfCc40iI1M+QMzqwYOAjojkMho9FERkQynriERkQynQiAikuFUCEREMpwKgYhIhlMhEBHJcCoEIiIZToVARCTD/S82K9DEUi5FJQAAAABJRU5ErkJggg==\n"
},
"metadata": {
"needs_background": "light"
}
},
{
"output_type": "display_data",
"data": {
"text/plain": "<Figure size 432x288 with 1 Axes>",
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZAAAAEGCAYAAABLgMOSAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAAAqmUlEQVR4nO3deXxV9Z3/8dcnNzskEJYECEFAkVVFjfsuCGpbsZvV1ha7Oe20M63tdMaOnbFjp49xOtM6ttZaa63WX622zmMU64KIW9WiRgUFWWWRnUAg+57P7497QmMMcHNz7z03yfv5eJzHvefcc8/5nLB88t3N3REREemtjLADEBGR/kkJRERE4qIEIiIicVECERGRuCiBiIhIXDLDDiCVRo0a5RMnTgw7DBGRfuX111/f6+6jux8fVAlk4sSJVFRUhB2GiEi/YmZbejquKiwREYmLEoiIiMRFCUREROKiBCIiInFRAhERkbiEmkDM7GIzW2tmG8zs+h4+zzGzB4PPXzGziV0++25wfK2ZzU9p4CIiEl4CMbMI8HPgEmAGcJWZzeh22heB/e5+DHAL8J/Bd2cAVwIzgYuB24PriYhIioQ5DuRUYIO7bwQwsweABcA7Xc5ZAHw/eP8QcJuZWXD8AXdvBjaZ2Ybgen9JRqD3vLSJqvqW9x80e/9uD9/rdgrW7azun/d0nQ9co6cvHeE7sdw7UfF/8Brx/JwOf41YvtNTcEf62UbMiGQYmREjMyODzAwjMxJ97Xo8kmFkRYJjGRlkRozcrAh5wZaTmUFGRgw/HJF+LswEUgps7bK/DTjtUOe4e5uZVQMjg+PLun23tKebmNm1wLUAEyZMiCvQ+199j/V76g7uawkVOZLcrIyDCSU3O/qanx1hWF4WhXlZDOu2Dc/PYtTQHIoLchk1NJvMiJonJf0N+JHo7n4ncCdAeXl5XP/1P3XdefHct9t+t89j+c4HPu/hOxz+Pj3Hdvhr9HxO9897+M4RrtHTQx8p/mT9nHra7XCnrd1p63DaOzpobXfaO6L7be0dwavT1tFBe4d3+byDptZ2GlvaaWztoLG1nebWdhoPHmunqbWd+uZ2dhxoYvXOWmoaW6ltbuvh6aIlo5FDciguyKGkMIfSojwmjhzCUSOHcNTIfCaMyCc3SzW2Er4wE8h2oKzL/vjgWE/nbDOzTGAYsC/G74bqA1U3MdVoqNpjMGlr76C2qY3qxlYONLZSWdvMntomdtc0Uxm87qpu4vUt+6lpen+yGTcsl2ljC5k+toBpYwqZPraQSaOGEFHVmaRQmAnkNWCKmU0i+p//lcCnu52zCFhItG3jE8Az7u5mtgi438x+AowDpgCvpixykQTIjGRQNCSboiHZRzz3QEMLm/c1sGVfPVv2NbCxso7VO2t5YV0lbR3R4tSQ7AgnTiiifGIRp0wcweyy4QzJGfCVDBKi0P52BW0aXwcWAxHgbndfZWY3ARXuvgj4NXBf0EheRTTJEJz3B6IN7m3A19y9PZQHEUmB4fnZzM7PZnbZ8Pcdb25rZ/3uOlbvrOGtbdVUbNnPrUvX4w6ZGcbJRxVxwbRiLpxWzJTioTF1whCJlfVUlz1QlZeXu2bjlYGupqmVN7bsZ9nGKp5bu4c1u2oBKB2ex/yZY1gwexzHjx+mZDJIVNW3UN3YylEj8uPuHWhmr7t7+QeOK4GIDGw7qxt5bm0lS1fv4YV1lbS0dzBxZD6XzS7lU6eUUTo8L+wQJYluf24DP3pyLWt+cHHcnS8OlUBUQSoywI0dlsdVp07gqlMnUN3YyuKVu1i0Yge3PbOe255Zz5zpJXz29KM4+5hRGr8yADU0txPJMHIyE981XAlEZBAZlpfFFaeUccUpZWw/0Mj9r2zhgVe3suSd3Rw9egh/d+EUPnz8WI1DGUDqmtvIz44kpcpSf0tEBqnS4Xl8Z/40Xv7uhfzPp2aTFcngmw8u56JbXuCh17fR1t4RdoiSAPXNbQxNUm88JRCRQS4nM8LlJ5by+N+fwx1Xn0xeVoR/+OMKLrn1zzy/rjLs8KSPGlrak9adWwlERADIyDAunjWGx/7+bO64+mRa2jtYePerfP43r7Kxsu7IF5C0VNfcxpDs5MxcoAQiIu9jFk0kT113LjdcOp2Kzfu5+NY/c8fz76paqx+qb25TCUREUisnM8KXz53M0m+fxwVTR3PzE2v46O0vszYYVyL9Q72qsEQkLMWFudxx9cnc/pmT2FndyGW3vcjvXtnS44Sakn7UiC4ioTIzLj1uLE9+81xOmzySG/5vJV+//02qG1vDDk2OoD7oxpsMSiAiErNRQ3O455pTuP6SaTy5aheX//wlNbCnuTqVQEQkXWRkGF8572geuPZ0ahpbufznL/HShr1hhyU9aGvvoLmtQ20gIpJeTpk4goe/dhZjh+Xxubtf5f5X3gs7JOmmviU6SbmqsEQk7ZSNyOehr57BuVNG8c//9za3P7ch7JCki/pg1cuCXJVARCQNFeRmcefnyrnshHH86Mm1/OeTa9RDK03UBQlkaE5WUq6vyRRFpM+yIhnc8qnZDM3N5BfPvUtjSzs3fmSG1hwJWW2wFPLQJJVAlEBEJCEiGcYPL59FflaEu17cRG5WhH+6eKqSSIj+WgJRAhGRNGdm3PCh6TS1tXPH8+9SkJvJ1y44JuywBq26puS2gSiBiEhCmRk3XTaLhuZ2/mvxWoZkR7jmrElhhzUo1TVHB3qqBCIi/UZGhvGjTxxPXXMb//andygtyueiGSVhhzXoJLsNRL2wRCQpMiMZ3HrliRxfOoxvPPAmq3ZUhx3SoNPZBjIkWwlERPqZvOwIv/pcOcPysvjiPRXsrmkKO6RBpa4puhZIJElr3SuBiEhSFRfm8uuFp1DT1Mq1v62gua097JAGjbokrgUCSiAikgIzxhVyy6dms2JbNf/+p9VhhzNo1Da3Ja39A5RARCRF5s8cw7XnTua+ZVtYtGJH2OEMCnVNbRSoBCIiA8F35k+l/Kgivvu/b/GupoFPunqVQERkoMiKZPCzT59ITlaEr/3uDZpa1R6STMlcCwSUQEQkxcYOy+PHV5zAml213LJkXdjhDGi1TW1Jm0gRlEBEJAQXTC3m06dN4M4/b+TVTVVhhzNg1TW3JW0aE1ACEZGQ3HDpdMqK8vn2H5cfHPAmiePuqsISkYFpSE4mP77iBLbtb+SHj6lrb6I1tXbQ3uFqRBeRgemUiSO49pzJ/P7V9/jLu/vCDmdAqU3yRIqgBCIiIfvm3GMpG5HHDQ+/rVHqCZTsqdxBCUREQpaXHeEHC2axsbKeO5/fGHY4A0ayF5MCJRARSQPnTy3mQ8eP5WfPbmDz3vqwwxkQOksgSiAiMuD964dnkBPJ4F8eWYm7hx1Ov1fbnNy1QEAJRETSRElhLt+edyx/Xr+Xp97ZHXY4/d7BNhANJBSRweDq04/imOKh/Mfjq2lp6wg7nH6tTiUQERlMMiMZ3PCh6Wze18Bv/7I57HD6tYOrEeZEknaPUBKImY0wsyVmtj54LTrEeQuDc9ab2cIux39oZlvNTNN5igwwF0wt5txjR/PTpevZX98Sdjj9Vm1TG9mRDHIyQ0ggZlZrZjU9bLVmVtPH+14PLHX3KcDSYL/7/UcANwKnAacCN3ZJNI8Gx0RkALrh0unUNbdx69L1YYfSb9U1tya1+goOk0DcvcDdC3vYCty9sI/3XQDcG7y/F7i8h3PmA0vcvcrd9wNLgIuD2Ja5+84+xiAiaWrqmAKuPHUC9y3bwiZ1641LXVNy58GCXlRhmVmxmU3o3Pp435IuCWAXUNLDOaXA1i7724JjvWJm15pZhZlVVFZW9j5SEQnFN+dOISti3Pq0pnyPR11ze/gJxMwuM7P1wCbgeWAz8EQM33vazFb2sC3oep5HO3wnrdO3u9/p7uXuXj569Ohk3UZEEqy4IJeFZ07kkRU7WLe7Nuxw+p1Qq7C6+AFwOrDO3ScBc4BlR/qSu89191k9bI8Au81sLEDwuqeHS2wHyrrsjw+Oicgg8ZVzj2ZIdqYWnopDXXNy10OH2BJIq7vvAzLMLMPdnwXK+3jfRUBnr6qFwCM9nLMYmGdmRUHj+bzgmIgMEkVDsvnC2ZN4YuUuVm6vDjucfqWuKbnroUNsCeSAmQ0FXgB+Z2a3An1t1boZuCioGpsb7GNm5WZ2F4C7VxEt/bwWbDcFxzCzH5nZNiDfzLaZ2ff7GI+IpKkvnTOJYXlZKoX0Um1TG0OSXAKJ5eoLgCbgOuAzwDDgpr7cNCjRzOnheAXwpS77dwN393DePwL/2JcYRKR/KMzN4m/Om8yPnlzL8q0HmF02POyQ0p67U9PUyrC85E1jAjGUQNy93t3b3b3N3e91958GCUBEJCU+d8ZEhuVlcfuzG8IOpV9oau2gtd0pzA0pgZjZi8Fr9wGFiRhIKCISs6E5mSw84yieemc369Uj64hqmqKrERbmhTeQ8OzgtfuAwkQMJBQR6ZVrzppEXlaEXzz/btihpL2axiCBhFUC6WRm98VyTEQkmUYMyebKU8tYtHwH2/Y3hB1OWussgYTeBgLM7LpjZpnAyckJR0Tk0L58zmTM4FcvaOnbw6nuLIGElUDM7LtmVgsc37X9A9hNz+M2RESSatzwPC6fXcoDr21lb11z2OGkrZrG6FTuhSFOpvgf7l4A/Fe39o+R7v7dpEYlInIIf3Pe0TS3dfC7Ze+FHUra+msjevjdeL9rZqVmdqaZndu5JTUqEZFDOKZ4KOdPHc19y7bQ3NYedjhpKZ0a0W8GXgK+B3wn2P4hqVGJiBzGF86axN66Zh57S6s69KSmqY28rAjZmcldMzCWCrKPAlPdXRWOIpIWzpkyimOKh/LrFzfx0RNLMbOwQ0or1Q2tSR8DArH1wtoIJLccJCLSC2bGF86axKodNby2eX/Y4aSdmqbWpFdfQWwJpAFYbma/NLOfdm7JDkxE5HA+emIpw/OzuPvFTWGHknZqmlqT3oAOsVVhLQo2EZG0kZcd4apTJ/DL599la1UDZSPyww4pbdQ0tjFqaHbS7xNLL6x7gT8Ay4LJFO8NjomIhOqzpx8FwP2vqktvV6mYiRdi64X1EWA58GSwP9vMVCIRkdCNG57HnOkl/LFiKy1tHWGHkzaqG1NThRVLG8j3gVOBAwDuvhyYnLSIRER64dOnTWBvXQtPvbMr7FDSgrtT05g+jeit7t59LUmlehFJC+dOGc34ojyNTA/Ut7TT4cmfyh1iSyCrzOzTQMTMppjZz4CXkxyXiEhMIhnGVadO4C8b9/FuZV3Y4YQuVaPQIbYE8ndEZ+RtBu4HqoFvJDMoEZHeuKK8jMwM4/evqBTSORNvWjSiAx9y9xvc/ZRg+x5wWbIDExGJ1eiCHObPHMNDb2yjqXVwz49Vk6Kp3CG2BNLTzLuajVdE0spnTpvAgYZWnlg5uOfHqmnqnMo9xIGEZnYJcClQ2m3keSHQluzARER64/TJI5kwIp+HXt/GR08cH3Y4oflrCSTcRvQdQAXQBLzeZVsEzE96ZCIivZCRYXz8pPG8/O6+Qb3k7cG1QMJsRHf3FcGI82OA3/PXBPKou2v2MhFJOx8/uRSA/319e8iRhKezEb0gyasRQmxtIGcC64GfA7cD67SglIiko/FF+Zx59EgeemMrHR0edjihqGlsY2hOJpmR5K4FArElkJ8A89z9PHc/l2j11S3JDUtEJD6fPLmMrVWNvLKpKuxQQhGdyj35pQ+ILYFkufvazh13X4fWBxGRNDV/5hgKcjL54+tbww4lFDUpmgcLYksgFWZ2l5mdH2y/Itq4LiKSdvKyI3z4hHE88fYu6poHX4fR6hTNgwWxJZCvAu8Afx9s7wTHRETS0ifLx9PY2s7jg3DN9AMNrQzPT00COWJFmbs3m9ltwFKikyiudfeWpEcmIhKnE8uGM3FkPo+s2M4Vp5SFHU5K7W9oYXbZ8JTcK5b1QD4EvAvcCtwGbAgGGYqIpCUzY8HsUl5+dx+7a5rCDidl3J0Dja0MH5I+VVg/Bi5w9/Pd/TzgAtQLS0TS3ILZ43CHR1fsCDuUlGlsbaelrYPheclfzhZiSyC17r6hy/5GoDZJ8YiIJMTk0UM5fvwwHlk+eBLI/oboIMKiFLWBxNoL63Ezu8bMFgKPAq+Z2cfM7GNJjk9EJG6XnTCOt7dXD5p1QvbXR5unh+enTwkkF9gNnAecD1QCecBHgA8nLTIRkT76yAnjMGPQlEIOpLgEEksvrM+nIhARkUQrKczlzKNHsmj5dq6bOwUzCzukpNrfEC2BFA1JnxKIiEi/teCEUjbva+CtbdVhh5J0Bxo6q7DSpw1ERKTfuvi4MWRnZvDw8oE/Q29nI3rovbDM7BvB61kpiUREJAkKc7O4cGoxj67YSVt7R9jhJNX+hhaG5mSSnZmassHh7tLZ9vGzRN/UzEaY2RIzWx+8Fh3ivIXBOeuDHmCYWb6ZPWZma8xslZndnOj4RGRguWz2OPbWNfPq5oE9Q28qpzGBwyeQ1Wa2HphqZm912d42s7f6eN/rgaXuPoXoFCnXdz/BzEYANwKnAacCN3ZJNP/t7tOAE4GzNDJeRA7ngqnF5GVFePztgT031v6GFopS1IUXDr8i4VXAOcAGol12O7cPB699sQC4N3h/L3B5D+fMB5a4e1WwAuIS4GJ3b3D3Z4MYW4A3gMG7ALKIHFFedoQLpxXz5MrdtA/ghab2p1EJBHff5e4nADuBgmDb4e5b+njfEnfv/FVgF1DSwzmlQNcJ/bcFxw4ys+FEk9nSQ93IzK41swozq6isrOxT0CLSf1163NhoNdYAXmjqQENLygYRQgzjQMzsPOC3wGbAgDIzW+juLxzhe08DY3r46IauO+7uZtbrXwnMLJPoWu0/dfeNhzrP3e8E7gQoLy8fuL96iMhhXTBtNLlZGTz+9k7OOHpk2OEkxYGG1pQNIoQYEgh/XdJ2LYCZHUv0P+6TD/cld597qM/MbLeZjXX3nWY2FtjTw2nbiY587zQeeK7L/p3Aenf/nxieQUQGufzsTC6cVswTK3fx/ctmEskYWIMK2zucmqbWlJZAwlrSdhGwMHi/EHikh3MWA/PMrChoPJ8XHMPM/h0YBnyzj3GIyCByyaxoNVbFAOyNVd3YinvqpjGB8Ja0vRm4KOjlNTfYx8zKzewuAHevAn4AvBZsN7l7lZmNJ1oNNgN4w8yWm9mX+hiPiAwCF04rJiczY0D2xjo4jUk6tYEQXb72a0SXswX4M3B7X27q7vuAOT0crwC+1GX/buDubudsI9oWIyLSK0NyMrlgarQa68aPzCRjAFVjpXoaE4ihBOLuze7+E3f/WLDd4u7NqQhORCTRLj1+LHtqm6nYsj/sUBJqf33nTLzp1QYiIjJgzBmg1VhhVGEpgYjIoDIkJ5Pzjh3N4lW7cB84Pfs71wJJ1Xro0IsEYmb5yQxERCRV5s8cw87qJt7ePnCmeN/f0EJmhlGQE0vTdmIcMYGY2Zlm9g6wJtg/wcz61IguIhKmC6cVE8kwnlq1O+xQEqZzGpNULpoVSwnkFqLzUu0DcPcVwLnJDEpEJJmKhmRz6sQRLF61K+xQEmZfXTMjh+Sk9J4xVWG5+9Zuh9qTEIuISMrMm1nC+j11bKysCzuUhNhb18yogtQ1oENsCWSrmZ0JuJllmdk/AKuTHJeISFLNmxmdqm/JOwOjGmtffUtalkC+QnQgYSnR+almB/siIv1W6fA8ZpUW8tRASSB1LYwamtoEEktzvbn7Z5IeiYhIis2bMYZbnl7Hntomigtyww4nbk2t7dQ1tzFyaPpVYb1kZk+Z2ReD9TdERAaEeTNLcIen3+lpQvD+Y29ddHKQUemWQNz9WOB7wEyikxf+ycyuTnpkIiJJNrWkgKNG5vPUO/27N9beuugo9FRXYcXaC+tVd/8W0bXJq/jrcrQiIv2WmTFvRgkvb9hHbVNr2OHEbV9QAhmZbgnEzArNbKGZPQG8THR521OTHpmISArMmzmGlvYOnlvbf5e83heUQEYOSbMqLGAF0Z5XN7n7se7+T+7+enLDEhFJjZMmFDFySHa/7o1VebANJP16YU32gTTjmIhIF5EM46IZJfzprZ20tHWQndn/5pjdV9fCkOwIedmRlN73kD8pM/uf4O0iM/vAlprwRESSb+70Euqa23hl076wQ4nLvvrmlLd/wOFLIPcFr/+dikBERMJy1jGjyMnMYOnqPZwzZXTY4fTa3rrmlHfhhcOUQLq0c8x29+e7bkTbREREBoS87AhnHzOKpWt298s1QvbVtYRSAomlsm9hD8euSXAcIiKhunB6MVurGlm/p/9Nrrg3hGlM4PBtIFeZ2aPApG7tH88SHQsiIjJgzJlWAsDTq/tXb6z2DqeqPpwqrMO1gXSO+RgF/LjL8VrgrWQGJSKSamOG5TKrtJBnVu/hb88/JuxwYnagoYUOT/0YEDhMAnH3LcAW4IzUhSMiEp4500r42TPrqapvYUQI/yHHY199MI1JQRpVYXUys9PN7DUzqzOzFjNrN7OaVAQnIpJKc6YX0+Hw7Jr+M7ni3tpgGpMUrwUCsTWi3wZcBawH8oAvAT9PZlAiImGYNW4YxQU5LF3Tf9pB9naWQNKpG29X7r4BiLh7u7v/Brg4uWGJiKReRoYxZ3oxL6zbS0tbR9jhxKSzBJJWvbC6aDCzbGC5mf3IzK6L8XsiIv3OnGnRUemvbuofnU131zaRHclgeH5Wyu8dSyL4bHDe14F6oAz4eDKDEhEJS+eo9P7SnXdPTTPFhTmYWcrvHUsCORlwd69x939z928FVVoiIgNOXnaEs/rRqPTdNU2UFIazHG8sCeQjwDozu8/MPmxmsczgKyLSb83pR6PSowkk9e0fENuStp8HjgH+SLQ31rtmdleyAxMRCUvnqPSlq9O/O++emmaKC9K3BIK7twJPAA8ArwOXJzEmEZFQjRmWy8xxhSxN83aQ+uY2apvb0rcKy8wuMbN7iI4D+ThwFzAmyXGJiIRqzvQS3nhvP1XBOIt0tCfowpu2VVjA54CHganufo27P+7ubckNS0QkXHP7waj03TVNAOlbheXuV7n7w+7enIqARETSQeeo9Gf6QQJJuxKImb0YvNaaWU2XrVZzYYnIQJeRYVw4rZjn11Wm7aj0yqAKqzjd2kDc/ezgtcDdC7tsBe5emLoQRUTCMWd6eo9K313TRG5WBoW54YyuOGwVlplFzGxNqoIREUknZ6f5qPRdNc2UFOaGMgodjpBA3L0dWGtmE1IUj4hI2kj3Uek7DjQydlg41VcQWy+sImCVmS3turRtX25qZiPMbImZrQ9eiw5x3sLgnPVmtrDL8SfNbIWZrTKzO8ws0pd4REQOJZ1Hpe880Mi4YXmh3T+WirN/ScJ9rweWuvvNZnZ9sP9PXU8wsxHAjUA54MDrZrbI3fcDV7h7jUXLbQ8BnyQ6yFFEJKHmTCvhBlby9OrdHFtSEHY4B7V3OLtrmxk3PLwEEks33ueBzUBW8P414I0+3ncBcG/w/l56Htk+H1ji7lVB0lhCsA6Ju3f2AssEsokmGBGRhOtcKz3dpjXZU9tEe4czdngaV2GZ2ZeJ/pb/y+BQKdGBhX1R4u47g/e7gJIezikFtnbZ3xYc64xrMbAHqA3i65GZXWtmFWZWUVlZ2cewRWQwmjMtOip9X136DIfbcaARIL1LIMDXgLOAGgB3Xw8UH+lLZva0ma3sYVvQ9TyPtkz1ugTh7vOBsUAOcOFhzrvT3cvdvXz06NG9vY2ICHOnl+AOz65Nn19CdxyIDiIMsw0klgTS7O4HJ4MJpnM/4n/47j7X3Wf1sD0C7DazscH1xhItSXS3nejiVZ3GB8e63qMJeIRolZiISFLMHFcYjEpPn+68O6ujJZC0rsICnjezfwbyzOwiotO6P9rH+y4COntVLSSaBLpbDMwzs6Kgl9Y8YLGZDe2SfDKBDwEaqyIiSZOOa6XvONBEQU4mhbmpX8q2UywJ5HqgEngb+BvgceB7fbzvzcBFZrYemBvsY2blnWuNuHsV8AOijfavATcFx4YAi8zsLWA50dLLHX2MR0TksDrXSn9l076wQwGCMSAhlj4ghm687t5hZg8DD7t7QioA3X0fMKeH4xXAl7rs3w3c3e2c3cApiYhDRCRWnWulL129h3OmhN+eurO6ibEhtn/A4SdTNDP7vpntBdYSHZFeaWb/mrrwRETSQ152hLOPGcXTq9NjVPq2/Q2UFqVpAgGuI9r76hR3H+HuI4DTgLPM7LqURCcikkbmTC9h2/5G1u0Od1R6bVMr+xtamTAiP9Q4DpdAPgtc5e6bOg+4+0bgaqKLTImIDCpzpkdHMIQ9ueLWqmgPrLKi9E0gWe6+t/vBoB0kvGZ/EZGQlBTmclzpsNDXSt+6vwEgrUsgh1sIOH0XCRYRSaI504t5c+sB9oY4Kn1rVTSBlI1I3zaQE7qtRHhwRULguFQFKCKSTg6OSg9xqdv3qhooyM1kWF64lUGHW5Ew0m0lwq4rEqoKS0QGpZnjCikpzAl1csWtVQ2UFeWHtpBUp1gGEoqISMDMuHBaCX9eX0lzW3soMbxX1RB6+wcogYiI9Nrc6cXUt7TzysbUr5Xe0eFs3d8YevsHKIGIiPTaWceMIjcrI5TeWHtqm2lp61AJRESkP8rN6hyVviflo9I3VkYHMU4aNTSl9+2JEoiISBzmzRzD9gONrNxec+STE2jj3noAJo8ektL79kQJREQkDhdNLyGSYTy5aueRT06gjZX15GZlMKYw3Jl4QQlERCQuRUOyOX3yCJ5cuSul9924t45Jo4aSkRFuF15QAhERidvFM8fwbmU963fXpuyem/bWp0X1FSiBiIjEbd7MMQApK4U0t7WztaqBo0cpgYiI9GslhbmcNGE4T65KTQLZWtVAh8MklUBERPq/S2aNZdWOmoMTHCbThj1BD6w06MILSiAiIn0yP4XVWGt31WIGU0qUQERE+r0JI/OZMbYwJdVYa3fXMGFEPvnZmUm/VyyUQERE+ujiWWN4fct+9tQ0JfU+a3bVMrWkIKn36A0lEBGRPrpkVrQaa3ESSyFNre1s3lvPtDFKICIiA8YxxUM5pngoj76VvFHpG/bU0eEwdUxh0u7RW0ogIiJ9ZGYsOGEcr26qYseBxqTcY82u6GDFqSqBiIgMLJfNHgfAoyt2JOX6a3bWkJOZwcSR4U/j3kkJREQkAY4aOYTZZcN5ZHlyEshb26qZMa6QzEj6/LedPpGIiPRzC2aP452dNWzYk9i5sdraO3h7ezUnjB+e0Ov2lRKIiEiCfOj4sWQYLEpwKWRDZR2Nre0cP35YQq/bV0ogIiIJUlyQy5lHj+Lh5TsSulLhW1urATheJRARkYHroyeW8l5VA69sqkrYNVdsO0BBTiaT02QW3k5KICIiCXTpcWMpyM3kwde2Juyab753gOPLhqXFIlJdKYGIiCRQXnaEy2eX8vjbO6luaO3z9aobWlm9q4ZTJ45MQHSJpQQiIpJgnzqljOa2Dh5evr3P13p1cxXucPrkEQmILLGUQEREEmxW6TBmlRby+1ff63Nj+rKN+8jJzOCEsuGJCS6BlEBERJLgylMmsGZXLW9tq+7TdV7ZtI+TJhSRmxVJUGSJowQiIpIEl80eR352hHtf3hz3NfbWNbNqRw1nHJ1+7R+gBCIikhSFuVl86pQyFq3Ywc7q+CZYfGb1HtxhzvTiBEeXGEogIiJJ8oWzJtHhzj1xlkKWrN5N6fA8ZoxNnyncuwolgZjZCDNbYmbrg9eiQ5y3MDhnvZkt7OHzRWa2MvkRi4j0XtmIfC6ZNZb7X3mPuua2Xn23saWdP6+vZO70YszSa/xHp7BKINcDS919CrA02H8fMxsB3AicBpwK3Ng10ZjZx4C61IQrIhKfa8+dTG1TG795cVOvvvf06t00tXYwb+aYJEXWd2ElkAXAvcH7e4HLezhnPrDE3avcfT+wBLgYwMyGAt8C/j35oYqIxO+EsuHMm1HCL1/YSFV9S8zf+0PFVkqH53HG5PRsQIfwEkiJu3eu/bgLKOnhnFKg61wA24JjAD8Afgw0JC1CEZEE+c78qTS0tHH7sxtiOn/b/gZe3LCXT5w8Pu2mL+kqaQnEzJ42s5U9bAu6nufRUTYxj7Qxs9nA0e7+fzGef62ZVZhZRWVlZa+eQUQkEaaUFPDxk8bz279sYWPlkWve73lpMwZ8snx88oPrg6QlEHef6+6zetgeAXab2ViA4HVPD5fYDpR12R8fHDsDKDezzcCLwLFm9txh4rjT3cvdvXz06NGJeTgRkV76zvyp5GZlcP3/vk1Hx6F/Z95d08R9y7Zw+YmljC9Kn+VrexJWFdYioLNX1ULgkR7OWQzMM7OioPF8HrDY3X/h7uPcfSJwNrDO3c9PQcwiInErLszlex+ewaubq7j9uUNXZf3kqXW0dTjfmDMlhdHFJ6wEcjNwkZmtB+YG+5hZuZndBeDuVUTbOl4LtpuCYyIi/dInTx7P5bPH8eMl63jsrZ0f+PyJt3fyYMVWvnT2JI4amV5rf/TEErlqVrorLy/3ioqKsMMQkUGsqbWdz9z1Cm++t59vz5vKF8+eRE5mBotX7ea6B5dzbMlQ/viVM8nOTJ9x3mb2uruXf+C4EoiISGo1trTzrT8s54mVuyjMzWRoTiY7qpuYOa6Q33z+FIoLcsMO8X0OlUAywwhGRGQwy8uO8IurT+bld/eyaPkOmlrbOW3ySD52Uik5mek36+6hKIGIiITkzKNHcebRo8IOI27pU8kmIiL9ihKIiIjERQlERETiogQiIiJxUQIREZG4KIGIiEhclEBERCQuSiAiIhKXQTWViZlVAlvi/PooYG8Cw0lXes6BRc85sIT1nEe5+wfWwxhUCaQvzKyip7lgBho958Ci5xxY0u05VYUlIiJxUQIREZG4KIHE7s6wA0gRPefAouccWNLqOdUGIiIicVEJRERE4qIEIiIicVEC6cbMLjaztWa2wcyu7+HzHDN7MPj8FTObGEKYfRbDc37LzN4xs7fMbKmZHRVGnH11pOfsct7HzczNLG26SPZGLM9pZlcEf6arzOz+VMeYCDH8vZ1gZs+a2ZvB391Lw4izr8zsbjPbY2YrD/G5mdlPg5/DW2Z2UqpjBMDdtQUbEAHeBSYD2cAKYEa3c/4WuCN4fyXwYNhxJ+k5LwDyg/dfHajPGZxXALwALAPKw447SX+eU4A3gaJgvzjsuJP0nHcCXw3ezwA2hx13nM96LnASsPIQn18KPAEYcDrwShhxqgTyfqcCG9x9o7u3AA8AC7qdswC4N3j/EDDHzCyFMSbCEZ/T3Z9194ZgdxkwPsUxJkIsf54APwD+E2hKZXAJFMtzfhn4ubvvB3D3PSmOMRFieU4HCoP3w4AdKYwvYdz9BaDqMKcsAH7rUcuA4WY2NjXR/ZUSyPuVAlu77G8LjvV4jru3AdXAyJRElzixPGdXXyT6205/c8TnDIr+Ze7+WCoDS7BY/jyPBY41s5fMbJmZXZyy6BInluf8PnC1mW0DHgf+LjWhpVxv/w0nRWaqbyj9i5ldDZQD54UdS6KZWQbwE+CakENJhUyi1VjnEy1NvmBmx7n7gTCDSoKrgHvc/cdmdgZwn5nNcveOsAMbiFQCeb/tQFmX/fHBsR7PMbNMosXkfSmJLnFieU7MbC5wA3CZuzenKLZEOtJzFgCzgOfMbDPRuuRF/bAhPZY/z23AIndvdfdNwDqiCaU/ieU5vwj8AcDd/wLkEp2AcKCJ6d9wsimBvN9rwBQzm2Rm2UQbyRd1O2cRsDB4/wngGQ9atfqRIz6nmZ0I/JJo8uiP9eVwhOd092p3H+XuE919ItG2nsvcvSKccOMWy9/bh4mWPjCzUUSrtDamMMZEiOU53wPmAJjZdKIJpDKlUabGIuBzQW+s04Fqd9+Z6iBUhdWFu7eZ2deBxUR7fNzt7qvM7Cagwt0XAb8mWizeQLSR68rwIo5PjM/5X8BQ4I9BH4H33P2y0IKOQ4zP2e/F+JyLgXlm9g7QDnzH3ftVyTnG5/w28Cszu45og/o1/fAXPMzs90QT/qigPedGIAvA3e8g2r5zKbABaAA+H0qc/fBnKyIiaUBVWCIiEhclEBERiYsSiIiIxEUJRERE4qIEIiIicVECEQmRmZ1vZn8K3l92uBmDRdKNxoGIJEEwwab1ZgqNYBzDgBibIoODSiAiCWJmE4O1Kn4LrAR+bWYVwfob/9blvIvNbI2ZvQF8rMvxa8zstuD9PWb2iS6f1QWvY83sBTNbbmYrzeyclD2gSDcqgYgk1hRgobsvM7MR7l5lZhFgqZkdT3QOql8BFxIdRfxgL6//aWCxu/8wuG5+IoMX6Q2VQEQSa0uwPgPAFUEp401gJtEFjqYBm9x9fTDFxv/r5fVfAz5vZt8HjnP32gTFLdJrSiAiiVUPYGaTgH8A5rj78cBjRCf2i1Ubwb/PYNr5bDi40NC5RGdevcfMPpe40EV6RwlEJDkKiSaTajMrAS4Jjq8BJprZ0cH+VYf4/mbg5OD9ZQQT6Vl0bfrd7v4r4C6iy56KhEJtICJJ4O4rzOxNogljK/BScLzJzK4FHjOzBuDPRNcl6e5XwCNmtgJ4kqBkQ3SG1u+YWStQB6gEIqHRbLwiIhIXVWGJiEhclEBERCQuSiAiIhIXJRAREYmLEoiIiMRFCUREROKiBCIiInH5/4VHUd4UeTHeAAAAAElFTkSuQmCC\n"
},
"metadata": {
"needs_background": "light"
}
}
]
},
{
"metadata": {
"trusted": true
},
"cell_type": "code",
"source": "",
"execution_count": null,
"outputs": []
}
],
"metadata": {
"kernelspec": {
"name": "python3",
"display_name": "Python 3",
"language": "python"
},
"language_info": {
"name": "python",
"version": "3.7.6",
"mimetype": "text/x-python",
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"pygments_lexer": "ipython3",
"nbconvert_exporter": "python",
"file_extension": ".py"
},
"gist": {
"id": "",
"data": {
"description": "polychrom_repulsive_potential",
"public": true
}
}
},
"nbformat": 4,
"nbformat_minor": 4
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment