Skip to content

Instantly share code, notes, and snippets.

@gotoken
Created March 2, 2020 05:27
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save gotoken/e769347d122dc74da14468ebffa270c7 to your computer and use it in GitHub Desktop.
Save gotoken/e769347d122dc74da14468ebffa270c7 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{"cells":[{"metadata":{"trusted":true},"cell_type":"code","source":"# original: http://epirecip.es/epicookbook/chapters/kr08/2_1/python_original\nimport scipy.integrate as spi\nimport numpy as np\nimport pylab as pl\n%matplotlib inline\n\nbeta=1.4247\ngamma=0.14286\nTS=1.0\nND=70.0\nS0=1-1e-6\nI0=1e-6\nINPUT = (S0, I0, 0.0)\n\ndef diff_eqs(INP,t): \n\t'''The main set of equations'''\n\tY=np.zeros((3))\n\tV = INP \n\tY[0] = - beta * V[0] * V[1]\n\tY[1] = beta * V[0] * V[1] - gamma * V[1]\n\tY[2] = gamma * V[1]\n\treturn Y # For odeint\n\nt_start = 0.0; t_end = ND; t_inc = TS\nt_range = np.arange(t_start, t_end+t_inc, t_inc)\nRES = spi.odeint(diff_eqs,INPUT,t_range)\n\n#Ploting\npl.subplot(211)\npl.plot(RES[:,0], '-g', label='Susceptibles')\npl.plot(RES[:,2], '-k', label='Recovereds')\npl.legend(loc=0)\npl.title('Program_2_1.py')\npl.xlabel('Time')\npl.ylabel('Susceptibles and Recovereds')\npl.subplot(212)\npl.plot(RES[:,1], '-r', label='Infectious')\npl.xlabel('Time')\npl.ylabel('Infectious')\npl.show()","execution_count":1,"outputs":[{"output_type":"display_data","data":{"text/plain":"<Figure size 432x288 with 2 Axes>","image/png":"iVBORw0KGgoAAAANSUhEUgAAAYIAAAEeCAYAAACHXhKxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nOzdd3xUZdbA8d9JAiT0koD0gFIEDAGCiBUsgIrYBdZG2bUtbV0VUEFfldVddXVRFFlF7FhXAcUCCIgFSBCQqjSlSQfpkOS8fzyTMIRJcgOZRs7Xz/OZuXfuvXMyhjxzn3IeUVWMMcaUXDHhDsAYY0x4WUVgjDElnFUExhhTwllFYIwxJZxVBMYYU8JZRWCMMSWcVQTGGFPCWUVgjDElXKEVgYj8S0QqikgpEZkqIltF5KZQBGeMMSb4vNwRdFLVP4CuwDqgMXBvUKMyxhgTMl4qglK+x8uAd1R1exDjMcYYE2JxHo6ZKCLLgP3AXSKSBBwIbljGGGNCRbwknRORKsAfqpolIuWACqr6e9CjM8YYE3T53hGIyDUB9vlvfhSMgIwxxoRWQU1DV/geqwNnA9N82x2B6VhFYIwxJ4V8O4tVtbeq9gYUaKaq16rqtUDzkEVnShwRWSMi+0Vkj4hsEpFXRaR8uOMqLiJyr4gsEpHdIrJaRDyNwBORR0XkJxHJFJGHgxymKWG8jBpKVtWNftubcENIjQmWK1S1PNAaaAs86P+iOMU2GVJEvAyaKLa3A24BqgBdgH4i0sPDeSuA+4BPgxibKaG8/GOaLiJfiEgvEbkV94v4dZDjMgZVXQ9MBlqIyHQRGSEi3wL7gIYiUktEJojIdhFZISJ/yTlXRBJE5DUR2SEiS0XkPhFZ5/f6GhEZLCILgb0iEiciQ0Rkpe/b+hIRudrv+F4i8q2IPCMiO0VklYic7du/VkQ2+/59FPYz/UtV56lqpqouBz4BzvFw3muqOhnYXdixIjJOREaLyFe+n2WGiNT3vTZKRJ7Oc/xEERlU2HXNyavQikBV+wGjgZZAKjBGVfsHOzBjRKQubv7Kj75dNwO3ARWAX4F3cJMcawHXAf8QkYt8xz4EJAMNgUuAQLPhewKXA5VVNRNYCZwHVAL+D3hTRGr6Hd8OWAhUA94GxuPuWE7zXf/5ojRjiRt9cR6w2Os5RXAj8CiQCMwH3vLtfw3omXNHJSKJwEW4z9KUVKpaaAHqAxf7npfFDR/1dK4VK0UpwBpgD7AT98f+BSABN0DhEb/j6gJZ/r+LwOPAON/zVUBnv9f+DKzL8z59CollPnCl73kv4Be/187A9Z/V8Nu3DUgtws/6f8ACoEwRznkTeLiQY8YB4/22y/s+q7q+7aXAJb7n/YDPwv3/3Up4i5dcQ38BPgBe8u2qDXxc2HnGnICrVLWyqtZX1btUdb9v/1q/Y2oB21XVv6nkV9zvZ87r/sf7Pw+4T0RuEZH5vqafnUAL3DfqHJv8nu8HUNW8+zzdEYhIP1xfweWqetDLOUWU+7Op6h5gO+4zAXdXkHOHdBPwRhDe30QRL30Ef8W1Yf4BoKq/4IaUGhNq/rMfNwBVRaSC3756wHrf841AHb/X6hZ0PV8b+n9x35CrqWplYBGuc7dYiUgfYAhwkaquK+z445T78/qaq6riPjNwdxVXikhL4HTsi12J56UiOKiqh3I2fCMsCp+ObEwQqepa4DvgcRGJF5EUoC9H2sLfA4aKSBURqY37A1+Qcrjf6y0AItIbd0dQrETkRuAfuKaZVUU4r5SIxOP+zcb5fubYAk65TETOFZHSuL6C2b7PDF/lMxd3J/Ch3x2XKaG8VAQzROR+IEFELgHeByYGNyxjPOmJ6xDeAPwPeEhVv/K99giuI3k1MAXXvJlvE4yqLgGeBr7HNQGdAXwbhJgfw3U2z/XNldgjIqM9nPdfXNNTT+AB3/ObAUTkPBHZk+f4t3Ed5tuBNrjOY3+v4X5GaxYyheca8o0u6At0wt0mfwG8rIWdaEwEEZE7gR6qekG4Ywk2ERmH6xh/sIBjzsc1ESWranaoYjORqcCJNL5bz9dU9SbcNxJjooJv2GdD3Df8RsDfgefDGlSEEJFSwEDcFzqrBEzBTUOqmgUk+doZjYkmpXEj3Xbj8mR9ghuKGnQistiv2ce/5G2e8T/nvHzOydvkc6KxnY4bmlsTeLY4r22il5emoZdwU/0nAHtz9qvqv4MbmjHGmFDwkmNlg6/E4GZ0GmOMOYl4WpgGQETKqerewo8MrsTERE1OTg53GMYYE1UyMjK2qmpSoNcKvSMQkfbAK7gZk/V8k1BuV9W7ijdMb5KTk0lPTw/HWxtjTNQSkV/ze83LPIJngc64PCqo6gLgfA9vOtaXkXFRPq+LiIz0ZY1cKCKtPcRijDGmmHnK6Z4zI9FPlofTxuHyrefnUtywvka4jJIveonFGGNM8fLSWbxWRM4G1DeMdAAue2GBVHWmiCQXcMiVwOu+iWk/iEhlEampRy+CU2wmLp/I7ZNu93SsyLHpZQRBRBCEGIkhRmKIjYklVmKJi4mjbKmyJJRKICEugfKly1OvUj0aVmlIg8oNaFilIY2qNSIuJpTrn5iSTlU5ePBgbjl06FDuo385fPhw7qN/yczMzLdkZWUd9TxQyc7OPuYx73P/oqrHPM+7L+/+ombZzPlcvDz3fwy0L9AxhR2b9/WC9gU65s4772Tw4MH5/08/Tl7+Mt0B/AeX1XEd8CUuEd2Jqs3R2R/X+fYdUxGIyG24uwbq1at3XG9Ws0JNujbuWuhxAf8H4ftF8v2XrdlkazZZ2VlkazaHsw+z//B+9h3ex44DO/ht129MXjGZfYf35V4juXIyfzvrb/Rt1Zdypcsd189gTj7Z2dns3buXXbt25ZY//viD3bt355Y9e/awd+/e3Me9e/eyb98+9u/fz759+3KfHzhwILfs37+fQ4cOFR5AMYiJiSE2NpZq1aoxbNgwGjZsSEzMkcYG/y9W+T0v6r6C9nt9/XiPDcb5XiUkJLB0acHfw+Pj46lTpw6lSpXyfF0v8wiSVHWL5ysefW4yMElVj0neJSKfAo+r6izf9lTgPlXNKOiaaWlpGg2dxarK5r2bWbVjFcu2LuOVH1/h27XfUiW+Cn9t+1f6t+tP9XKWxPVkkp2dzY4dO9i0aRObN29m06ZNbNmyhW3btrF161a2bt3Ktm3b2LFjR27ZuXMn2dmFT+6Ni4ujXLlyR5WyZcuSkJBA2bJliY+PJyEhgYSEBOLj4ylTpgzx8fG5z3NK6dKlKV26NGXKlKFUqVK526VKlTqmxMXF5T7mlNjY2KO2Y2Jicv8Irl69mgoVKlCtWrWQ/WE0R1NVtm3bxu7du2nQoMFRr4lIhqqmBTrPyx3BdyKyGngXl6lw54mHC7g7AP/UwHU4kiY36okINcrXoEb5GrSv257erXrz3drvePK7JxnxzQjGzBvDgjsWcEr5U8IdqvEgOzubjRs38ttvv+WWtWvXsn79ejZs2MD69evZuHEjmZmZAc+vUqUKiYmJVK1alaSkJBo3bkyVKlWoXLkylStXplKlSrmlYsWKVKhQgfLly+c+lilTJsQ/cdEdOHCA5ORkqwTCSESoVq0aW7YU7bt7oRWBqjYSkTOBHsADIrIEt/rRm8cXaq4JuIW7x+OWANwVrP6BSHF23bP5X/f/kbEhg3PGnsMdk+7gf93/Z/9wIkR2djZr165l6dKlLF++nBUrVrBy5UpWrVrF6tWrj2lqqVixIrVr16ZWrVp07NiRmjVrcsopp1CjRg1q1KhB9erVSUpKolq1asTFlYz+IftdDr/j+X/g6bdTVecAc0TkH8C/cSlsC6wIROQdoAOQKG7R8IeAUr7rjQY+w61HuwK3GHnvIkcfpdrUasNjFz7GvV/dy9s/vc2NKfmmoDFBsmnTJhYsWMDChQtZsGABixYtYvny5ezffyQ1f4UKFTj11FNp0aIFV155JcnJydSvX5969epRt25dKlWqFMafwJji42VCWUXgatwdwam4vO9nFnaeqvYs5HWleDqdo9LfzvobHy39iP6T+9OxQUdqVahV+EnmuOzYsYO5c+cyd+5c5syZw5w5c/j9999zX69duzYtWrSgQ4cOnH766TRt2pSmTZuSlJRk33Cj0IgRI3j77beJjY0lJiaGl156iXbt2oU8jvnz57NhwwYuu+wyACZMmMCSJUsYMmQIvXr1omvXrlx33XVHnTN9+nSeeuopJk2aFNJYvdwRLMAtZfeIqn4f5HhKjNiYWMZdNY6Wo1ty+6TbmdBjgv3RKSYbNmxg5syZzJw5k2+++YZFi47MaWzatCmXXHIJrVu3pmXLlpxxxhkkJiYWcDUTTb7//nsmTZrEvHnzKFOmDFu3bg3Z6Km85s+fT3p6em5F0K1bN7p16xaWWArjZUJZQ1X9G7DIt/apKSaNqzXm8YseZ9LPk3h9wevhDidq7du3j8mTJzNw4ECaNm1K7dq16dmzJ2+88QZ16tRhxIgRTJkyhZ07d7J06VJef/11Bg0aRMeOHa0SOMls3LiRxMTE3M71xMREatWqRXJyMlu3bgUgPT2dDh06ADBjxgxSU1NJTU2lVatW7N69G4B//etfnHHGGbRs2ZIhQ4YAsHLlSrp06UKbNm0477zzWLZsGQC9evXijjvu4LzzzqNx48ZMmjSJQ4cOMXz4cN59911SU1N59913GTduHP36HVkxdcqUKUedk9fevXvp06cPbdu2pVWrVnzyyScALF68mDPPPJPU1FRSUlL45ZdfTvhz83JH0FxE3sAtfi0isgW4VVUDpo4wRTOg3QA+WvoRAz8fSKdTO1GzQs1whxQVfv/9dz755BM++ugjZsyYwcGDB4mPj+eCCy7gL3/5Cx06dKBly5YlppM20gz6fBDzf59frNdMPSWVZ7sUvIRCp06deOSRR2jcuDEXX3wx3bt354IL8l+U7qmnnmLUqFGcc8457Nmzh/j4eCZPnszHH3/M7NmzKVu2LNu3bwfgtttuY/To0TRq1IjZs2dz1113MW3aNADWrFnDjBkzWLlyJR07dmTFihU88sgjpKen8/zzbj2kcePGHfXegc7xN2LECC688ELGjh3Lzp07OfPMM7n44osZPXo0AwcO5MYbb+TQoUNkZXlJ9FAwL/9KxgB3q+rXACLSwbfv7BN+d0OMxPBS15do9kIzPljyAf3b9Q93SBFr48aNvPvuu3z44Yd8++23qCqnnXYad911F126dOG8884jISEh3GGaMCpfvjwZGRl88803fP3113Tv3p0nnngi3+PPOecc7r77bm688UauueYa6tSpw5QpU+jduzdly5YFoGrVquzZs4fvvvuO66+/PvfcgwePLIF9ww03EBMTQ6NGjWjYsGHu3UJBCjvnyy+/ZMKECTz11FOAG57722+/0b59e0aMGMG6deu45ppraNSoUZE+o0C8VATlcioBAFWdLiI2NbYYnZ50Og0qN2Dq6qlWEeSxf/9+PvnkE1577TW+/PJLsrOzSUlJ4eGHH+aaa66hefPm1rcSgQr75h5MsbGxdOjQgQ4dOnDGGWfw2muvERcXlztx78CBA7nHDhkyhMsvv5zPPvuMs846iylTpqCqx/xOZWdnU7lyZebPD3yXk/d4L7+ThZ2jqnz44Yc0adLkqP2nn3467dq149NPP6Vz5868/PLLXHjhhYW+X0G89BGsEpFhIpLsKw8Cq0/oXc0xLmpwEdPXTCcr+8Rv804GS5cu5a9//Ss1a9akZ8+eLF68mCFDhrBs2TIWLFjA8OHDadGihVUC5ijLly8/qs18/vz51K9fn+TkZDIyXNKCDz/8MPf1lStXcsYZZzB48GDS0tJYtmwZnTp1YuzYsezb51LEbN++nYoVK9KgQQPef/99wP2RXrBgQe513n//fbKzs3PnnTRp0oQKFSrk9jkEEugcf507d+a5557LTXvz448/ArBq1SoaNmzIgAED6NatGwsXLjyRjwzwVhH0AZKAj3wlkRI05j9ULmxwIbsO7mLexnnhDiVssrOzmTRpEp06daJZs2a8/PLLdO3alSlTprBmzRpGjBhxzD8WY/zt2bOHW2+9lWbNmpGSksKSJUt4+OGHeeihhxg4cCDnnXcesbGxucc/++yztGjRgpYtW5KQkMCll15Kly5d6NatG2lpaaSmpuY2zbz11lu88sortGzZkubNm+d23gI0adKECy64gEsvvZTRo0cTHx9Px44dWbJkSW5ncV6BzvE3bNgwDh8+TEpKCi1atGDYsGEAvPvuu7Ro0YLU1FSWLVvGLbfccuIfXFGz94W7tGnTRk9Gv+/+XXkYfeKbJ8IdSsgdOnRIx44dq6eddpoCWqtWLX3sscd006ZN4Q7NFMGSJUvCHUJY3Hrrrfr++++HO4yjBPp/AaRrPn9XC70jEJGvRKSy33YVEfnixKsg469G+Ro0T2rO1NVTwx1KyBw6dIj//ve/NG7cmD59+lCxYkXGjx/PmjVreOCBB6he3ZLyGRMKXjqLE9Uv0Zyq7hAR+xcaBBc2uJCX573MwcyDlImL/CRjxys7O5s333yTYcOG8dtvv9G2bVuef/55LrvsMmvzN1En77DQaOSljyBbRHIXARCR+oC3Fe9NkVzU4CL2Z+5n9vrZ4Q4laL799lvatWvHrbfeSvXq1Zk8eTKzZ8/m8ssvt0rAmDDxUhE8AMwSkTd8E8tmAkODG1bJdEHyBcRIDFNXnXzNQ7/++ivdu3fn3HPPZePGjbz++uvMnj2bLl26WAVgTJgVWhGo6udAa9x6BO8BbVTV+giCoHJ8ZdrUbMO0NdPCHUqxyc7OZtSoUTRv3pyJEyfy0EMPsXz5cm6++eajVrEyxoSPl85iwS1C31pVJwJlfesTmCC4sMGF/LDuB/Yc2hPuUE7YL7/8QocOHejXrx/nnnsuS5cu5eGHH6ZcOZuPaEwk8fKV7AWgPZCTVno3MCpoEZVwFzW4iMzsTGb9NivcoRy37Oxsnn76aVJSUvjpp5949dVXmTx5MvXr1w93aOYkFxsbS2pqKi1atOCKK65g587iWlAxtDp06EAol+T1UhG0U9W/AgfAjRoCSgc1qhLsnHrnUDq2dNT2E2zdupXLL7+ce+65h06dOrF48WJ69epl/QAmJBISEpg/fz6LFi2iatWqjBoVOd9Z81vGNBJ4qQgOi0gsvpFCIpIEFL7atjkuZUuVpX2d9lHZT/Dtt9/SqlUrpk2bxgsvvMDHH39MrVq24I4Jj/bt27N+/frc7SeffJK2bduSkpLCQw89lLv/9ddfJyUlhZYtW3LzzTcDbnDDRRddREpKChdddBG//fYbu3btIjk5OTdn0b59+6hbty6HDx8uMEX13XffTceOHRk8eHC+qaX3799Pjx49SElJoXv37rkr5WVlZdGrVy9atGjBGWecwTPPPBOUz8rLPIKRuFXJqovICOA64MGgRGMA10/w8PSH2b5/O1UTqoY7nEKpKk8//TRDhgyhfv36fP/997Ru3TrcYZkwGjRoUL4J2o5Xamoqzz7rLZldVlYWU6dOpW/fvoDL5PnLL78wZ84cVJVu3boxc+ZMqlWrxogRI/j2229JTEzMTTndr18/brnlFm699VbGjh3LgAED+Pjjj2nZsiUzZsygY8eOTJw4kc6dO1OqVKkCU1T//PPPTJkyhdjYWO6///6AqaVfeuklypYty8KFC1m4cGHuv5/58+ezfv363MWVgtXU5WXU0FvAfcDjwEbgKlV9PyjRGMD1EyjK9DXTwx1KoQ4ePMhNN93Evffey5VXXsm8efOsEjBhs3//flJTU6lWrRrbt2/nkksuAVxF8OWXX9KqVStat27NsmXL+OWXX5g2bRrXXXdd7gJFVau6L17ff/89f/rTnwC4+eabmTXL9dl17949N2/Q+PHj6d69+1EpqlNTU7n99tvZuHFjbkzXX399bn6jL7/8kieeeILU1FQ6dOiQm1p65syZ3HTTTQCkpKSQkpICQMOGDVm1ahX9+/fn888/p2LFikH53LwuXr8MWAYgIpVF5AFVHRGUiAxta7elXKlyTFs9jWtOvybc4eRr+/btXH311cycOZMRI0YwdOhQ6wswAJ6/uRe3nD6CXbt20bVrV0aNGsWAAQNQVYYOHcrtt99+1PEjR44sUsrobt26MXToULZv305GRgYXXnghe/fuLTBFtf8oOc0ntbT/e/irUqUKCxYs4IsvvmDUqFG89957jB07ttB4iyrfOwIRqSsiY0Rkkoj8WUTKisjTwM+ApZgIotKxpWlds3Wxr/BUnFavXs0555zDDz/8wFtvvcX9999vlYCJGJUqVWLkyJE89dRTHD58mM6dOzN27Fj27HHDstevX8/mzZu56KKLeO+999i2bRtAbtPQ2Wefzfjx4wGXdfTcc88F3MI3Z555JgMHDqRr167ExsYWmqLaX36ppc8//3zeeustABYtWpSbWnrr1q1kZ2dz7bXX8uijjzJvXnCyExd0R/A6MAP4EDeP4AdgMZCiqr8HJRqTq1lSM95b/F7ARTLCbe7cuXTt2pXDhw/z1Vdfcf7554c7JGOO0apVK1q2bMn48eO5+eabWbp0Ke3btwfcH/Q333yT5s2b88ADD3DBBRcQGxtLq1atGDduHCNHjqRPnz48+eSTJCUl8eqrr+Zet3v37lx//fVMnz49d99bb73FnXfeyWOPPcbhw4fp0aMHLVu2PCamYcOGMWjQIFJSUlBVkpOTmTRpEnfeeSe9e/cmJSWF1NRUzjzTTdVav349vXv3zu2gfvzxx4PyWUlOzXTMCyILVLWl3/YmoJ6qHgx4QoikpaVpKMfXhsvI2SMZ+PlANv59I6eUPyXc4eT64Ycf6NSpE9WqVWPy5Mk0bdo03CGZCLF06VJOP/30cIdhCPz/QkQyVDUt0PEF9hGISBUg5+vo77hZxeUAVHX7iYdr8tMsqRkAS7YsiZiKIKcSqF69OtOnT6dOnTrhDskYUwwKqggqARkcqQgAchqoFGgYrKDM0RXBhQ1ObD3S4jB79mw6d+5slYAxJ6F8KwJVTQ5hHCaPmuVrUqlMJZZsWRLuUJg9ezadOnUiKSnJKgFToEjs0ypp8mvuL4ilf4xQIkKzpGZhrwh++uknOnfuTGJiIl9//bVVAiZf8fHxbNu27bj+EJnioaps27btmPWPC+NpHoEJj2ZJzZiwfELY3n/Dhg1cfvnllC1blq+//pq6deuGLRYT+erUqcO6devYsmVLuEMp0eLj44v8hc0qggjWPKk5r/z4Clv2biGpXFJI33vPnj107dqV7du3880331CvXr3CTzIlWqlSpWjQoEG4wzDHId+KQEQKTHJjo4aCL6fDeOnWpSGtCDIzM+nRowcLFixg4sSJtGrVKmTvbYwJvYL6CDKAdN/jFtyM4l98zzOCH5rxHzkUKqrKgAED+PTTTxk1ahSXXXZZyN7bGBMe+VYEqtpAVRsCXwBXqGqiqlYDugIfhSrAkqxOxTqUL10+pBXBqFGjePHFF7n33nu54447Qva+xpjw8TJqqK2qfpazoaqTgQuCF5LJEeqRQ3PnzuXuu+/m8ssv54knngjJexpjws9LRbBVRB4UkWQRqS8iDwDbgh2YcZolNWPxlsVBf58dO3Zwww03ULNmTV577TVbWN6YEsTLv/aeQBJucZqPcZlHexZ4hik2zRKb8fue39m+P3h986pK7969WbduHe+++y7VqlUL2nsZYyKPl4VptqvqQFVt5SsDvY4YEpEuIrJcRFaIyJAAr3cQkV0iMt9Xhh/PD3Eyyx05tGVp0N7j2Wef5ZNPPuFf//oXZ511VtDexxgTmQqdRyAijYF7gGT/41W1wAQ4vnWORwGXAOuAuSIyQVXzNnh/o6pdixh3ieE/cuiceucU+/V/+OEH7rvvPq666ioGDRpU7Nc3xkQ+LxPK3gdGAy8DWUW49pnAClVdBSAi44ErgfAnz4ki9SvXJyEuISgdxnv27OFPf/oTderUYezYsZYjxpgSyktFkKmqLx7HtWsDa/221wHtAhzXXkQWABuAe1Q1+D2jUSRGYjg96XSWbC3+imDo0KGsWbOGGTNmUKVKlWK/vjEmOnjpLJ4oIneJSE0RqZpTPJwX6Otl3mxU84D6vgVwnsN1Rh97IZHbRCRdRNJLYh6TYAwhnTlzJs8//zz9+vXjvPPOK9ZrG2Oii5eK4FbgXuA73IzinBnHhVkH+Gcpq4P71p9LVf9Q1T2+558BpUQkMe+FVHWMqqapalpSUmhz7kSC5knNWffHOv44+EexXG/fvn306dOHBg0aBG3pO2NM9Ci0aUhVjzeL1FygkYg0ANYDPYA/+R8gIqcAm1RVReRMXMVkcxTy8B851K5OoNa1onnwwQdZuXIl06ZNo1y5cid8PWNMdPOUfVREWgDNgNwk16r6ekHnqGqmiPTDpaiIBcaq6mIRucP3+mjgOuBOEckE9gM91JKZH8N/5NCJVgTfffcdzz77LHfeeScdO3YsjvCMMVHOy/DRh4AOuIrgM+BSYBZQYEUAuc09n+XZN9rv+fPA80WKuARqULkBZWLLnHA/wYEDB+jTpw/16tXjn//8ZzFFZ4yJdl7uCK4DWgI/qmpvEamBG0pqQiQ2JpamiU1PeOTQM888w/Lly/n888+pUKFCMUVnjIl2XjqL96tqNpApIhWBzdjC9SF3oiOH1q9fz4gRI7jqqqvo3LlzMUZmjIl2XiqCdBGpDPwXN2JoHjAnqFGZY7So3oI1O9ew68Cu4zp/yJAhZGZm8vTTTxdzZMaYaOcl19BdqrrT17Z/CXCrqvYOfmjGX5uabQDI2Fj0NYG+++473nzzTe655x4aNrSbOWPM0YqUa1hV16jqwmAFY/KXVisNgLnr5xbpvOzsbAYMGEDt2rUZOnRoMEIzxkQ5W7w+SlQrW40GlRuQvtHLXL4jXn31VTIyMnjrrbdszoAxJiBbfSSKtK3dtkh3BDt37mTo0KGcc8459OxpS0gYYwIrtCIQkVNFpIzveQcRGeDrPDYh1rZWW37d9Stb9nrLt/T444+zdetWRo4caZlFjTH58nJH8CGQJSKnAa8ADYC3gxqVCSinnwM/Pg0AACAASURBVCB9Q+HNQ+vXr2fkyJHcdNNNtG7dOtihGWOimJeKIFtVM4GrgWdV9W9AzeCGZQJpU7MNgjB3Q+HNQ48++ihZWVn83//9XwgiM8ZEMy8VwWER6YnLQjrJt69U8EIy+alQpgJNE5sWWhGsWLGCV155hdtuu40GDY43Z6AxpqTwUhH0BtoDI1R1tS+b6JvBDcvkJ61WGukb0ikoN9/w4cMpXbo0Dz74YAgjM8ZEKy8TypYAg3EzilHV1ar6RLADM4G1rdWW3/f8zvrd6wO+vmDBAt555x0GDhzIKaecEuLojDHRyMuooSuA+cDnvu1UEZkQ7MBMYG1rtwXyn1j2wAMPULlyZe69995QhmWMiWJemoYexi1EvxNAVefjRg6ZMGhZoyVxMXEBRw7NmjWLTz/9lMGDB9saxMYYz7xUBJmqmjfTmS0eEyYJpRJoUb3FMR3Gqsr9999PzZo1GTBgQJiiM8ZEIy8VwSIR+RMQKyKNROQ53PrFJkza1mp7TIfxlClT+Oabb3jggQcoW7ZsGKMzxkQbLxVBf6A5cBB4B/gDGBTMoEzB0mqlsePADlbtWAW4u4Hhw4dTt25d/vznP4c5OmNMtPGyeP0+4AFfMRGgbS1fh/GGuZxa9VQmT57MDz/8wEsvvUSZMmXCHJ0xJtrkWxGIyEQK6AtQ1W5BicgUqkX1FpSJLUP6hnS6N+/O8OHDadCgAb172zIRxpiiK+iO4KmQRWGKpFRsKVJPSWXuhrlMmDCBjIwMxo4dS6lSNuHbGFN0+VYEqjoj57mIlAaa4u4QlqvqoRDEZgrQtlZbxs4by/BnhtOoUSNuvvnmcIdkjIlSXiaUXQ6sBEYCzwMrROTSYAdmCta2dlv2LdzHwoULeeihh4iLszWGjDHHx8tfj6eBjqq6Atz6BMCnwORgBmYK1rZmW5gONZJr0KNHj3CHY4yJYl6Gj27OqQR8VgGbgxSP8ei7id/BFsi+IBtszRljzAnItyIQkWtE5BpgsYh8JiK9RORWYCJQtBXUTbHatWsX999/P01aNWFL8hYm/jwx3CEZY6JYQXcEV/hKPLAJuADoAGwBLJFNGD366KNs2bKF1156jXqV6/HsD8+GOyRjTBQraNSQDUqPQMuXL+c///kPffr0oV3bdvQ/3J97v7qX+b/PJ/WU1HCHZ4yJQgU1Dd3ne3xOREbmLaEL0fi7++67KVu2LCNGjACgb6u+lCtVjv/M/k+YIzPGRKuCmoaW+h7TgYwAxYTYZ599xmeffcbw4cOpUaMGAFUSqtArtRdv//Q2m/ZsCnOExpholG9FoKo5PZD7VPU1/wLsC014JsehQ4f429/+RuPGjenfv/9Rrw1oN4BDWYcYnT46TNEZY6KZl+GjQz3uM0H06KOP8vPPP/PMM89QunTpo15rXK0xlzW6jBfSX+Bg5sEwRWiMiVYF9RFc6lt7oHae/oFxQGbIIjSMGTOGxx57jF69enHZZZcFPGZQu0Fs3ruZdxa9E+LojDHRrqA7gg24/oEDHN03MAHoHPzQDMDHH3/MnXfeyaWXXsqYMWPyPe7ihheTUiOFAZMH8OXKL0MYoTEm2hXUR7DA1x9wGm5Bmh+BecAkVd0RovhKtFmzZtGzZ0/S0tJ4//33C8wuKiJM6jmJ5MrJXP725Yz9cWwIIzXGRDMvfQSXYEnnQm7RokVcccUV1KtXj08//ZRy5coVek7dSnWZ1WcWHZM70ndCX4ZNG3bUcpbGGBOIl4rg37ikcx1U9QKgI/CMl4uLSBcRWS4iK0RkSIDXxdfvsEJEFopI66KFf/LZtGkTgwcPpn379iQkJPDFF1+QmJjo+fyKZSry6Z8+pU9qHx775jFu+OAGvl/7PdmaHcSojTHRLGhJ50QkFhgFXAo0A3qKSLM8h10KNPKV24AXvQR9Mlq7di0DBgwgOTmZp556iiuuuIKZM2eSnJxc5GuVii3Fy91eZsSFI/hk2SecPfZs6j1Tj4GTBzLz15nsPri7+H8AY0zUksKaDkTkRaA+8B5uYZrrgeXAtwCq+lE+57UHHlbVzr7tob7jH/c75iVguqq+49teDnRQ1Y35xZOWlqbp6elef75cmZmZHDhwoMjn5RXo81JVVJXs7Ozc51lZWWRlZZGZmUlWVhb79+9nz549uWXTpk0sXbo0t/z666/ExsZyyy23MGTIEBo1anTCsQLsOrCLiT9P5MOlHzL5l8kczHLDSxPLJtKwSkMaVmlIzfI1qVimIpXKVKJSfCUqlqlImdgylIkrQ5nYMsTHxRMXE5dbYmNiiZVYYiQmt4iIe0SOeg6u/8L/OXDMdg4JkErVyzGFnVMcvLyvMcFUOrY0ZeKOb11yEclQ1bRAr3lZj8A/6Ry4pHNVcQnpFAhYEQC1gbV+2+uAdh6OqQ3kWxEcr//973/ccMMNxX3ZExIfH0+TJk0466yz6Nu3LzfffDP169cv1veoFF+Jm1Ju4qaUm9h9cDdTV09l+dblrNqxilU7VzFn/Rw27dnE3sN7i/V9jTHFb/A5g3ni4ieK/bqFVgQnkHwu0NenvF+nvRyDiNyGazqiXr16xxVMSkoKTz755HGdGyCeY/bFxLhvxTklNjaW2NhY4uLiiI2NJSEhgfLly+eWatWqUa9ePWJjY4slJi8qlKnAVU2vCvhaZnYmuw/uZtfBXew+uJuDWQc5kHmAg5kHOZh1kMzsTLKys8jMznTPNcvdBWl2blHcds5+AEVz76LU978273aOgHdbHo4p7JziYJ3uJhK0rd02KNcttCIQkca4tvsaqtpCRFKAbqr6WCGnrgPq+m3Xwc1NKOoxqOoYYAy4pqHCYg6kSZMmNGnS5HhOLRHiYuKoklCFKgmWYdyYksZLZ/F/cSklDgOo6kLAy9qIc4FGItJAREr7zpmQ55gJwC2+0UNnAbsK6h8wxhhT/Lz0EZRV1Tl5mkMKTTGhqpki0g/4AogFxqrqYhG5w/f6aOAz4DJgBS6Rna2BYIwxIeZl1NBkoB/wvqq2FpHrgL6qGpZJZSKyBfj1OE9PBLYWYzjBFk3xRlOsEF3xRlOsEF3xRlOscGLx1lfVpEAveKkIGuLa588GdgCrgZtUdc1xBhM2IpKe3/CpSBRN8UZTrBBd8UZTrBBd8UZTrBC8eL2MGloFXCwi5YAYVbXZSMYYcxIptLNYRP4hIpVVda+q7haRKiJS2IghY4wxUcLLqKFLVXVnzoYv82jgpPiRL/88zpEpmuKNplghuuKNplghuuKNplghSPF66SNYCLRV1YO+7QQgXVWbByMgY4wxoeVl+OibwFQReRU367cP8FpQozLGGBMyhd4RgEsnDVyMSwnxpap+EezAjDHGhIaXzuJyuD/+9+Dap8qISP5LZUWowtZGCDcRGSsim0Vkkd++qiLylYj84nuMiPwPIlJXRL4WkaUislhEBvr2R1y8IhIvInNEZIEv1v+L1FhziEisiPwoIpN825Ec6xoR+UlE5otIum9fJMdbWUQ+EJFlvt/f9pEYr4g08X2mOeUPERkUrFi9dBbPBOJFpDYwBTf7d1xxvHmoeFwbIdzGAV3y7BsCTFXVRsBU33YkyAT+rqqnA2cBf/V9npEY70HgQlVtCaQCXXzpTCIx1hwDgaV+25EcK7iFq1L9xrdHcrz/AT5X1aZAS9znHHHxqupy32eaCrTBZV74H8GKNSd/fn4FmOd77A/c53v+Y2HnRVIB2gNf+G0PBYaGO64AcSYDi/y2lwM1fc9rAsvDHWM+cX+CW9I0ouMFyuLW3W4XqbHiEi9OBS7ErQ8e0b8HwBogMc++iIwXqIibECvREK9ffJ2Ab4MZq5c7AhG3yMyNwKe+fV46mSNJfuseRLoa6kvC53usHuZ4jiEiyUArYDYRGq+vqWU+bmW9r1Q1YmMFngXuA/zXFo3UWMENIPlSRDJ86eIhcuNtiFtP5VVf09vLvqbvSI03Rw/gHd/zoMTqpSIYhPsG/T91SeMaAl8Xx5uHkKd1D0zRiEh54ENgkKr+Ee548qOqWepusesAZ4pIi3DHFIiIdMUtDZsR7liK4BxVbY1rdv2riJwf7oAKEAe0Bl5U1VbAXiKgGaggvszN3YD3g/k+hVYEqjpDVbup6j9926tUdUAwgwoCT+seRKBNIlITwPdY6FrRoeIbMPAh8JYeWa40YuMFUDcxcjquLyYSYz0H6CYia4DxwIUi8iaRGSsAqrrB97gZ14Z9JpEb7zpgne+OEOADXMUQqfGCq2Dnqeom33ZQYvUyauhrEZmWtxTHm4eQl7URItEE4Fbf81txbfFhJyICvAIsVdV/+70UcfGKSJKIVPY9T8ANg15GBMaqqkNVtY6qJuN+R6ep6k1EYKzgRhSKSIWc57i27EVEaLyq+juwVkRyVqi6CFhChMbr05MjzUIQpFi9zCxu47cZD1wLZKrqfcURQKiIyGW49tectRFGhDmko4jIO0AHXJrZTcBDwMfAe0A94DfgelXdHq4Yc4jIucA3wE8cacu+H9dPEFHxiltR7zXc//cY4D1VfUREqhFhsfoTkQ7AParaNVJj9TUT/8+3GQe8raojIjVeABFJBV4GSgOrcKMgY4jAeEWkLK5vs6Gq7vLtC8pn62lCWYAAZ6jqBYUfaYwxJtJ5WbO4qt9mDG5M6ylBi8gYY0xIeRkGmoEbYSO4iUSrgb7BDMoYY0zoHFfTkDHGmJNHvncEItIWWOvraUdEbsF1FP8KPByuzpTExERNTk4Ox1sbY0zUysjI2Kr5rFlcUNPQS7ihdvgmiTyBSzORiks+d10xx+lJcnIy6enp4XhrY4yJWiLya36vFVQRxPp96+8OjFHVD4EPfdP1jTHGnAQKmlAWKyI5FcVFgP8ksmjLNVQyrFgBmyNpUqQxJhoUVBG8A8wQkU+A/bgJRIjIacCuEMRmvFq4EK67Dho1ggsugAMHwh2RMSaK5FsR+Gbe/h2XJ/9cPTK8KAbXV2DCbf58uPpqaNkSvvoKevWCZcvgscfCHZkxJooU2MSjqj8E2Pdz8MIxnq1ZA+3aQUICPPQQDBwIVapAdjb8859w/fWugjDGmEJ4SUNtItE338ChQzB9Ojz8sKsEAP79b6haFfr2hczMcEZojIkSVhFEq9mzoVw5OOOMo/dXqwbPPw8ZGfDMM+GJzRgTVawiiFZz5kBaGsTGHvvaddfBVVfB8OHwyy+hj80YE1WCWhGISBcRWS4iK0Qk4EpAItJBROaLyGIRmRHMeE4aBw+6juJ27QK/LgKjRkGZMvCXv4ClETHGFCBoFYGIxAKjcCvsNAN6ikizPMdUBl4Auqlqc+D6YMVzUlmwAA4fhjPPzP+YWrXgiSdgxgz44ovQxWaMiTrBvCM4E1jhW9ryEG7pvSvzHPMn4CNV/Q1yl7szhZkzxz0WVBEA9OkDdevCiIhag8cYE2GCWRHUxq2uk2Odb5+/xkAVEZkuIhm+xHbHEJHbRCRdRNK3bNkSpHCjyJw5cMopUKdOwceVLg333guzZrlRRsYYE0AwKwIJsC9vY3UcbqGby4HOwDARaXzMSapjVDVNVdOSkgImzytZ5sxxdwMS6CPOo29fSEqyuwJjTL6CWRGsA+r6bdcBNgQ45nNV3auqW4GZgM2CKsiOHbB8eeHNQjnKloW773b9BBkZwY3NGBOVglkRzAUaiUgDESkN9AAm5DnmE+A8EYnzLdTcDlgaxJiiX04K7vxGDAVy551QqRL84x/BickYE9WCVhGoaibQD/gC98f9PVVdLCJ3iMgdvmOWAp8DC4E5wMuquihYMZ0UcjqK09K8n1OpEvTrBx99BEuWBCcuY0zUirqlKtPS0rREL0xz5ZWuaWjZsqKdt3Ur1K8P114Lr78enNiMMRFLRDJUNeA3SJtZHE1UXWoJr/0D/hIT4bbb4O23YfXq4o/NGBO1rCKIJuvWwaZNx1cRANxzD8TEwL/+VbxxGWOimlUE0cTrRLL81K7t1iwYOxY25B3AZYwpqawiiCZz5rhJYieyzsDgwS499b//XXxxGWOimlUE0WT2bEhNdcnkjtepp0LPnjB6NGzbVnyxGWOillUE0SIry80hON5mIX9Dh8LevfCf/5z4tYwxUc8qgmixdKn7410cFUHz5m6t4+eegz/+OPHrGWOimlUE0eJEO4rzeuAB2LkTXnyxeK5njIlanioCESknIjG+541FpJuIlApuaOYoP/3k8gY1alQ812vTBjp3dp3G+/cXzzWNMVHJ6x3BTCBeRGoDU4HewLhgBWUCWLUKGjZ08wCKy/33w+bN8PLLxXdNY0zU8fpXRVR1H3AN8JyqXo1bdcyESk5FUJzOPx/OO8+tZGZ3BcaUWJ4rAhFpD9wIfOrbFxeckMwxVINTEQA89pibXPbCC8V/bWNMVPBaEQwChgL/82UQbQh8HbywzFE2b4Z9+9wcgOJ2/vmur+Dxx20EkTEllKeKQFVnqGo3Vf2nb3uVqg4Ibmgm16pV7jEYdwTgVi/bts1mGxtTQnlq3hGRrzl2mUlU9cJij8gca+VK9xisiqBNG7juOnj6abduQWJicN7HGBORvDYN3QPc6yvDgPlACV4UIMRy7giSk4P3Ho884pqfnngieO9hjIlIXpuGMvzKt6p6N25ZSRMKq1a5zKHx8cF7j9NPh1tugeefd+mujTElhtcJZVX9SqKIdAZOCXJsJkewRgzl9dBDkJ3tRhIZY0oMr01DGbimoAzge+DvQN/CThKRLiKyXERWiMiQAo5rKyJZInKdx3hKllBVBMnJcPvtboLZ0qXBfz9jTETw2jTUQFUb+h4bqWonVZ1V0DkiEguMAi7FTT7rKSLHTELzHfdP3CL3Jq8DB2D9+uAMHQ1k2DAoXx4GDnTzF4wxJz2vTUOlRGSAiHzgK/085Bo6E1jhG2p6CBgPXBnguP7Ah8DmIkVeUqxZ4x5DcUcAUL266zj+6iv45JPQvKcxJqy8Ng29CLQBXvCVNr59BakNrPXbXufbl8uXu+hqYHRBFxKR20QkXUTSt2zZ4jHkk0Swh44Gctdd0KIF/O1vlnrCmBLAa0XQVlVvVdVpvtIbaFvIORJgX962hmeBwaqaVdCFVHWMqqapalpSUpLHkE8SwZ5MFkhcHIwc6e5GnnwydO9rjAkLrxVBlojkNlL7UkwU+McbdwdQ12+7DpB3xfQ0YLyIrAGuA14Qkas8xlQyrFrl0k9Xrx7a9+3YEa6/3qWe+PXX0L63MSakvFYE9wJfi8h0EZkBTMONHCrIXKCRiDQQkdJAD2CC/wG+zudkVU0GPgDuUtWPi/QTnOxyRgxJoBusIHvqKfe+99wT+vc2xoSMpxQTqjpVRBoBTXBNPstU9WAh52SKSD/caKBYYKwvYd0dvtcL7BcwPqtWhW7EUF716rk1C4YNc53Hl1wSnjiMMUElWsAQQRG5UFWnicg1gV5X1Y+CFlk+0tLSND29hGS3UHVDOW+/PXwJ4Q4cgJQUOHTIrZJWoUJ44jDGnBARyVDVtECvFdY0dIHv8YoApWuxRWgCy0k/HcqO4rzi42HcOPjtN7j33vDFYYwJmgKbhlT1Id/TR1R1tf9rItIgaFEZJxxDRwM5+2z4+99dn8G111oTkTEnGa+dxR8G2PdBcQZiAgjH0NH8PPIING0KffvaAjbGnGQKvCMQkaZAc6BSnn6CikAQU2EaIDTpp71KSHBNRDl3B//9b7gjMsYUk8LuCJrg+gIqc3T/QGvgL8ENzYQk/XRRtGvn+glefhk+/zzc0RhjikmBo4ZyDxJpr6rfhyCeQpWoUUPnn+/G8c+YEe5IjjhwwK1otn07zJsHNWuGOyJjjAcnMmooxx0iUtnvglVEZGyxRGfyF6r000URHw/jx8OuXdCjB2RmhjsiY8wJ8loRpKjqzpwNVd0BtApOSAY4kn460ioCgDPOgJdegpkz4cEHwx2NMeYEea0IYkSkSs6GiFTF46xkc5xW+0brRmJFAHDzzW6i2z//CRMmFH68MSZief1j/jTwnYh8gMsgegMwImhRmcgaOpqfZ5+F9HS31nFGRvhSYRhjTojXFcpeB64FNgFbgGtU9Y1gBlbiRUNFEB8PH3wAMTFw3XWwd2+4IzLGHAevTUMAVYG9qvocsMVmFgfZqlVQrlzo008XVXIyvPUWLFwIPXta57ExUcjrUpUPAYOBob5dpYA3gxWUIbzpp4vq0kvh+edh4kTo39/WOjYmynjtI7gaN0poHoCqbhARS0MZTCtXRleb+513usR0Tzzh0lcPHVr4OcaYiOC1aeiQuplnCiAi5YIXkmH/fli2zA3TjCYjRsCNN7o1DN6wLiRjooXXO4L3ROQloLKI/AXoA1iymWD56SfIynIzeKNJTAyMHQsbN0KfPpCUBF26hDsqY0whCrwjEJEyAKr6FC7b6Ie4/EPDfZ3GJhgyMtxj69bhjeN4lC4NH30ELVrAVVfBZ5+FOyJjTCEKaxr6HkBE3lDVr1T1XlW9R1W/CkFsJde8eVC1qmtrj0aVKsHUqdC8OVx9tetENsZErMIqgtIicitwtohck7cUdnER6SIiy0VkhYgMCfD6jSKy0Fe+E5GWx/uDnFTmzXPNQtEwYig/VavClClumctrr4WPPw53RMaYfBRWEdwBnMWxaagLXapSRGKBUcClQDOgp4g0y3PYauACVU0BHgXGFPUHOOkcPOj6CKKxWSivKlXcovetW8P118P774c7ImNMAIUtVTkLmCUi6ar6ShGvfSawQlVXAYjIeOBKYInf9b/zO/4HoE4R3+Pks3gxHD4cfR3F+alcGb780s016N4d1q2DQYOi+27HmJOMp1FDqvqKiJwNJPuf40s9kZ/awFq/7XVAuwKO7wtMDvSCiNwG3AZQL1rbzb2K5o7i/FSs6O4MbrkF7r4bfv4ZnnsO4ixvoTGRwOvM4jeAp4Bzgba+EnCBA//TAuwLOOVURDriKoLBgV5X1TGqmqaqaUlJSV5Cjl7z5rnO1kjOMXQ8ypaF996DwYNh9Gi4/HK3poExJuy8fiVLA5qpl+XMjlgH1PXbrgNsyHuQiKQALwOXquq2Ilz/5DRvnrsbOBmbTmJi3Mzj005zM5HPOccNNW3cONyRGVOieZ1ZvAg4pYjXngs0EpEGIlIa6AEclbheROoBHwE3q+rPRbz+yefwYViw4ORqFgrkz392ax5v3Oj6Qt5+O9wRGVOieb0jSASWiMgc4GDOTlXtlt8JqpopIv2AL4BYYKyqLhaRO3yvjwaGA9WAF8R9A87Mb03NEmHpUjdq6GTpKC7IRRe5Sq9nT5eWYupUGDnSZVw1xoSU14rg4eO5uKp+BnyWZ99ov+d/Bv58PNc+Kc2b5x5P9juCHHXqwNdfw8MPwz/+AT/84O4OWtp0EmNCyevCNDMClWAHV+JkZED58tCoUbgjCZ24OHjsMTfEdNs2SEtz6yAfOBDuyIwpMQrLNbRbRP4IUHaLyB+hCrLEmDcPWrVynaolzcUXuzkUN97ospimpsI334Q7KmNKhAL/4qhqBVWtGKBUUNWKoQqyRMjKgvnzS0b/QH6qVYNx4+CLL1xfyfnnw+23w+bN4Y7MmJNaCfzqGaGWL4d9+0pO/0BBOnWCRYvc5LOxY91w08cfd+s0GGOKnVUEkaKkdRQXplw5ePppVyF07OgWu2na1K2PnJUV7uiMOalYRRApMjIgIcH9sTNHNGkCn3ziRhclJsJNN7m1Dt54AzIzwx2dMScFqwgixbx5roM0NjbckUSmDh1g7lx49123+M0tt7gZyWPGuP4EY8xxs4ogEmRnw48/luyOYi9iYuCGG1yn+oQJbinM22+HunXhgQdg7drCr2GMOYZVBJFg8WLYvdsNHTWFE4ErrnAT0KZMgbPPdjmMkpPhmmtcplPrRzDGM6sIIsG//w3x8XDZZeGOJLqIuFQVH38Mq1bBfffBzJlu1FH9+m77p5/CHaUxEc8qgnBbscJ1fN5xB5xS1Lx+Jlf9+m6I6bp1rh+hVSt45hm3VGZqqpu9vHgxFCmBrjElg1UE4TZiBJQq5b69mhMXH+/6ESZOhA0bXCK7hAQYNsyNNmrSxK2JMGuWy/ZqjEGKtsRA+KWlpWl6enq4wygeK1a44aL9+7tvryZ4Nmxww1A//himTXNDTytUcHMUOnVyKS4aNz4514EwBhCRjPyyO1tFEE69e8P48a59u2bNcEdTcuza5SqDL75wZc0at796dTj3XDjvPFdSUtzdmjEngYIqAls0Nlxy+gb697dKINQqVYKrr3ZFFVaudBXDrFku0d1HH7njypRx/Qtt27qsqK1bu6al0qXDG78xxczuCMLF7gYi17p18O23bgLb3Llu1vfeve61UqVcc94ZZ7jStKkrp55qdw8molnTUKSxvoHokpXlkgIuWAALF7ohqQsXHj2BLTbWVQannQYNGx4pyclQrx5Urmz9DyasrCKIJJ9+Cn/9q0utvHKl3Q1Esz/+cBXEsmVHHletcv9f/8izXEe5cq5CqFsXatU6utSo4fonatRwCxNZhWGCIGx9BCLSBfgPbs3il1X1iTyvi+/1y4B9QC9VnRfMmMJm3ToYONC1Pzdr5mbEWiUQ3SpWdP0HbdsevV8VduxwFcKvv8Jvv7m7h99+c2XxYvj998Czn+PjXeqMatVckr1q1VypUsWVypXdY6VKrlSseOQxPt4qEXNcglYRiEgsMAq4BFgHzBWRCaq6xO+wS4FGvtIOeNH3GP1U3R//n35yqRCeecYNWfzHP+Dvf7cOx5OZCFSt6kreSiJHVhZs3eqGtW7a5O4QN21yZetWt2zntm2uItm2DXbudDmpChIb6+4oKlRwj+XKHV3KlnUlIeHIY3z8saVMmaNL6dJHl1KljjzmlLg4q4SiWDDvCM4EVqjqKgARGQ9cCfhXBFcCr6trn/pBRCqLSE1V3Vjs0fz8s2uW8ReoWSzvPtVjS3a2obYPYwAAByBJREFU+6Oemen+QR886HIF7d4Ne/bA9u3uW9/OnUeuc9ll8Nxzrt3YmNhY1xRUo4a341Xd79eOHa788YcbBptT/H//cp7v2+c6ubdscUNk9+93+3JKYRXL8fxMOZVCTilVyu2Pizv6MVCJicn/0b/k7BM5+jHQc/8SaJ/XAt6f58i7r7DHvM8Dbbdv7+a+FLNgVgS1Af90kOs49tt+oGNqA0dVBCJyG3AbQL169Y4vmgUL3IpXxS0uzn07yvkWVqGCu1Xv0cONKklJcTNaK1cu/vc2JYeIa/6pWNGl0zhRqm5m9YEDR8r+/e5LjX85dMgV/+eHDx95zCmZmUeeZ2Ud2ef/hSnvo3/Jzj768dAh9zynZGW5mP2P8f9ilvOa/76c/XmPK2rJ+by8PPf/fINh8OCoqwgC3Sfm/XS8HIOqjgHGgOssPq5orrrq6G/ouREECCFQrexfcr7ZlMRF5s3JQeRIU09FW3486PJWFl4qj0CVSZD+5gSzIlgH1PXbrgNsOI5jikepUu6bujHGhFqgZqAIEsyvtHOBRiLSQERKAz2ACXmOmQDcIs5ZwK6g9A8YY4zJV9DuCFQ1U0T6AV/gho+OVdXFInKH7/XRwGe4oaMrcMNHewcrHmOMMYFF3YQyEdkC/HqcpycCW4sxnGCLpnijKVaIrnijKVaIrnijKVY4sXjrq2pSoBeiriI4ESKSnt/MukgUTfFGU6wQXfFGU6wQXfFGU6wQvHht2IsxxpRwVhEYY0wJV9IqgjHhDqCIoineaIoVoiveaIoVoiveaIoVghRvieojMMYYc6ySdkdgjDEmjxJTEYhIFxFZLiIrRGRIuOPJS0TGishmEVnkt6+qiHwlIr/4HquEM8YcIlJXRL4WkaUislhEBvr2R1y8IhIvInNEZIEv1v+L1FhziEisiPwoIpN825Ec6xoR+UlE5otIum9fJMdbWUQ+EJFlvt/f9pEYr4g08X2mOeUPERkUrFhLREXglxL7UqAZ0FNEmoU3qmOMA7rk2TcEmKqqjYCpvu1IkAn8Xf+/vbsJraMKwzj+f6QutLZ+1KqBLoKCX6i5rRqsEWljESnSjRtFpQs3Qv1aCUEQly5E7cpNxY1SUesXLjSi6MKF1NbaRmsqxYJF2oggWpUi7evinGvGmtTU5N45eJ4fDHNm7oU8Geby3jln7pmIK4AbgE35eJaY9ygwGhFDQAe4Lf+KvcSsXQ8DexvbJWcFWBsRncZtjSXn3Qy8GxGXA0Ok41xc3oiYzMe0A1xL+sHtG/Qqa0T87xdgNfBeY3sMGGs71ww5B4GJxvYkMJDbA8Bk2xlnyf0W6bkTRecFzgR2kmbBLTIrab6tD4BR4J3SzwPgAHD+CfuKzAssBb4lj42WnreR71bgk15mreKKgNmnuy7dhZHnXsrrC1rO8w+SBoGVwKcUmjd3tewCpoD3I6LYrMCzwKNA82EBpWaFNFvwuKQdebp4KDfvxcAPwAu5622LpMWUm7frTmBrbvckay2FYE7TXdupkXQWsA14JCJ+/rf3tyUijkW6xF4BDEu6qu1MM5F0OzAVETvaznIKRiJiFanbdZOkm9sOdBKLgFXAcxGxEviVArqBTiZP2LkBeLWXf6eWQtC/6a4X1mFJAwB5PdVynr9IOp1UBF6KiNfz7mLzAkTET8BHpLGYErOOABskHQBeBkYlvUiZWQGIiO/zeorUhz1MuXkPAgfzFSHAa6TCUGpeSAV2Z0Qczts9yVpLIZjLlNglehvYmNsbSX3xrZMk4Hlgb0Q83XipuLySlks6J7fPANYBX1Ng1ogYi4gVETFIOkc/jIh7KDArgKTFkpZ026S+7AkKzRsRh4DvJF2Wd91CenRukXmzu5juFoJeZW17IKSPAy7rgX3AfuCxtvPMkG8r6RGdf5C+udwHLCMNHH6T1+e1nTNnvYnUtbYb2JWX9SXmBa4BPs9ZJ4DH8/7isp6Qew3Tg8VFZiX1uX+Rly+7n6tS8+ZsHeCzfD68CZxbal7SzQ0/Amc39vUkq39ZbGZWuVq6hszMbBYuBGZmlXMhMDOrnAuBmVnlXAjMzCq3qO0AZqWS1L1VD+Ai4BhpigKA3yLixlaCmS0w3z5qNgeSngCORMRTbWcxW2juGjL7DyQdyes1kj6W9IqkfZKelHR3fgbCHkmX5Pctl7RN0va8jLT7H5hNcyEwm78h0jMErgbuBS6NiGFgC/Bgfs9m4JmIuB64I79mVgSPEZjN3/bIUwNL2g+M5/17gLW5vQ64Mk3TBMBSSUsi4pe+JjWbgQuB2fwdbbSPN7aPM/0ZOw1YHRG/9zOY2Vy4a8isP8aBB7obkjotZjH7GxcCs/54CLhO0m5JXwH3tx3IrMu3j5qZVc5XBGZmlXMhMDOrnAuBmVnlXAjMzCrnQmBmVjkXAjOzyrkQmJlVzoXAzKxyfwKMi8OOpUftiAAAAABJRU5ErkJggg==\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.6.7","mimetype":"text/x-python","codemirror_mode":{"name":"ipython","version":3},"pygments_lexer":"ipython3","nbconvert_exporter":"python","file_extension":".py"}},"nbformat":4,"nbformat_minor":4}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment