{ "cells": [ { "cell_type": "code", "execution_count": 49, "id": "d063afc4", "metadata": {}, "outputs": [], "source": [ "# importing standard libraries\n", "import numpy as np\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 50, "id": "1159aacf", "metadata": {}, "outputs": [], "source": [ "# setting up quantities\n", "dt=0.1 # the time interval in days\n", "t=np.arange(0,200+dt/2,dt)\n", "nNi=1000 # the initial number of Ni atoms" ] }, { "cell_type": "code", "execution_count": 51, "id": "f3b66e55", "metadata": {}, "outputs": [], "source": [ "# a python function that computes the derivative\n", "def nuclDecay(t,nCo,nNi0):\n", " return(0.113*nNi0*np.exp(-t/8.8)-0.009*nCo)" ] }, { "cell_type": "code", "execution_count": 52, "id": "8c60e8a6", "metadata": {}, "outputs": [], "source": [ "# first doing Euler (RK1) method\n", "nCoEu=np.zeros(t.size)\n", "nCoEu[0]=0 # the initial condition\n", "for i in range(1,t.size):\n", " nCoEu[i] = nCoEu[i-1]+dt*nuclDecay(t[i-1],nCoEu[i-1],nNi)" ] }, { "cell_type": "code", "execution_count": 53, "id": "735c903b", "metadata": {}, "outputs": [], "source": [ "# now trying an RK2 method\n", "nCoRK2=np.zeros(t.size)\n", "nCoRK2[0]=0 # the initial condition\n", "for i in range(1,t.size):\n", " k1=nuclDecay(t[i-1],nCoRK2[i-1],nNi)\n", " k2=nuclDecay(t[i],nCoRK2[i-1]+k1*dt,nNi)\n", " nCoRK2[i] = nCoRK2[i-1]+dt/2*(k1+k2)" ] }, { "cell_type": "code", "execution_count": 54, "id": "851d8b0b", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0, 0.5, 'Number of $^{56}$Co atoms')" ] }, "execution_count": 54, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAEGCAYAAACQO2mwAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAA/DUlEQVR4nO3dd3hUZfbA8e+ZmfRCCISOAopilGpABMSKHbAvrgUU5aeCbW1Y1rbriq4NUFSsWAEVBRtFFBULGKQ3ASmGGmoIhNTz+2Nu3AETyITM3CE5n+eZZ+a+c8vJzWRObnnPK6qKMcYYUxaP2wEYY4yJXJYkjDHGlMuShDHGmHJZkjDGGFMuSxLGGGPK5XM7gKpUt25dbdasmdthGGPMIWXWrFmbVTWtrPeqVZJo1qwZmZmZbodhjDGHFBFZXd57drrJGGNMuSxJGGOMKZclCWOMMeWqVtckjDGRobCwkKysLPbs2eN2KCZAbGwsTZo0ISoqqsLLWJIwxlS5rKwskpKSaNasGSLidjgGUFW2bNlCVlYWzZs3r/BydrrJGFPl9uzZQ506dSxBRBARoU6dOkEf3YU1SYjI7SKyUEQWiMj7IhIrIqkiMkVEljnPtQPmv1dElovIUhE5K5yxGmMOjiWIyFOZ30nYTjeJSGPgFiBdVfNEZCzQB0gHpqrqEBEZDAwG7hGRdOf9Y4FGwFcicpSqFocr5nBTVVZk5zIvawebduZTUFRCSnwUTVPjadskhdSEaLdDNMbUMOG+JuED4kSkEIgH1gH3Aqc4748CpgH3AL2B0aqaD6wUkeVAJ+CnMMcccrsLihj142pG/7KGrC07OULWUU+2E00h2zWRPzSNbGqT3jCZ89o05ML2jWmUEud22MZENK/XS+vWrf+c7tOnD4MHDy53/jfffJPMzEyef/75cIT3F9OmTSM6OpouXboA0K9fP84//3wuueSSSq3vYJcvFbYkoaprReQpYA2QB0xW1ckiUl9V1zvzrBeRes4ijYGfA1aR5bTtRUQGAAMADjvssFD+CCHxzZJN3PPhXFrt/oUnk3/g+IRZ+Irz/jJfbkx9ZuS15p0p7Rk+pQ1ntWnKjaccQasGyS5EbUzki4uLY86cOSFbf1FRET5f1X2FTps2jcTExD+TRKQI2zUJ51pDb6A5/tNHCSJy5f4WKaPtL8PoqepIVc1Q1Yy0tDJLj0QkVeWZyUt5dNQEXtWHeCv6CU7wLsXX4Qq4cCRcMxGu+xqu+AjOHkLikV05XX/hjej/MjPhHzRbPJK/Df2Suz6Yy4YddpuhMRXVrFkzNm/eDEBmZiannHLKX+bJzs7m4osvpmPHjnTs2JEffvgBgIcffpgBAwZw5plncvXVV++1zLRp0zj55JO57LLLOOqooxg8eDDvvvsunTp1onXr1qxYsaLcda9atYqXXnqJZ599lnbt2vH9998D8N1339GlSxdatGjBhx9+CPi/O+666y6OO+44WrduzZgxY/5sHzRoEOnp6Zx33nls2rSpSvZXOE83nQGsVNVsABEZB3QBNopIQ+cooiFQ+pNlAU0Dlm+C//TUIU9VeWjCQrbOGM2k2FeIio6DU5+CDn3BV8Z1h5ZnQOcboSgfln9F8oyXuX3le9wQ/xlD5/Wix7yzuPXM47ima3O8HrtYaCLLI58uZNG6nCpdZ3qjZB7qeex+58nLy6Ndu3Z/Tt9777387W9/q9D6b731Vm6//Xa6devGmjVrOOuss1i8eDEAs2bNYvr06cTF/fWU79y5c1m8eDGpqam0aNGC6667jpkzZzJ06FCGDx/Oc889V+66b7jhBhITE7nzzjsBeO2111i/fj3Tp09nyZIl9OrVi0suuYRx48YxZ84c5s6dy+bNm+nYsSPdu3fnp59+YunSpcyfP5+NGzeSnp7OtddeW8E9Wr5wJok1QGcRicd/uul0IBPYBfQFhjjP4535JwDvicgz+I88WgIzwxhvyLzwzXJKZr7K8Og3oWln5JLXIbnRgRf0xUCr8/yP9XOJ++oRBq94h/6+r7jzy758Ou9knry4DUc3SAr5z2BMpDuY001fffUVixYt+nM6JyeHnTt3AtCrV68yEwRAx44dadiwIQBHHHEEZ555JgCtW7fmm2++OeC693XBBRfg8XhIT09n48aNAEyfPp3LL78cr9dL/fr1Ofnkk/nll1/47rvv/mxv1KgRp512WqV+9n2F85rEDBH5EPgVKAJmAyOBRGCsiPTHn0gudeZf6NwBtciZf2B1uLPpq0UbWfrVmwyPfgM96mzk0jchqhIXoRu2havGwYqvqfvlPYwqeoIvNv/I1c9fxaDzTuDKzofbLYgmIhzoP/5w8/l8lJSUAJTbZ6CkpISffvqpzGSQkJBQ7rpjYmL+fO3xeP6c9ng8FBUVHXDd+1ufqu71XJZQ/M2HtZ+Eqj6kqq1U9ThVvUpV81V1i6qerqotneetAfM/pqpHqOrRqvplOGMNhU079/DGh+N4OvplSpp2Ri57q3IJItARpyE3TIeT7+Ec+ZnJ0Xcz+dP3ueGdWWzfXVA1gRtTjTRr1oxZs2YB8NFHH5U5z5lnnrnXXU5VeQG8vHUnJSWVe0QRqHv37owZM4bi4mKys7P57rvv6NSpE927d2f06NEUFxezfv36P49cDpb1uA6jf4/7hX8XP4ckpuG5/H3/6aOq4IuBU+9D/u9bkus05O3oIWT89iy9h37DgrU7qmYbxhxiSq9JlD5Kb3996KGHuPXWWznppJPwer1lLjts2DAyMzNp06YN6enpvPTSS1UWV3nr7tmzJx9//PFeF67LcuGFF9KmTRvatm3LaaedxpNPPkmDBg248MILadmyJa1bt+bGG2/k5JNPrpJ4ZX+HLoeajIwMjdRBh6Yv28zvo27gKt9XSL/PoFm30GyoYDdMug9mvcFCaclNhbfyj0tOo3e7v9w9bEzILF68mGOOOcbtMEwZyvrdiMgsVc0oa347kgiDouIS3h0/gSt9X1Hc8frQJQiA6Hjo+Rxc9hbHRK1nfPQDvDNmNI9/sZjikurzD4ExJjwsSYTB+NlruSbnJYqiU/Cddn94NpreG8/131ArpS6jYx5j1w8vc+0bM9m5pzA82zfGVAuWJEKspESZPfU9OnmWEtXjQYhLCd/G045CBnyNt+UZ/DvqDc5d9TiXvzTdOt8ZYyrMkkSITV64nkt2jWFXfBOkw9UHXqCqxdaCy0fDSXfwN+833Lftn1z5wmSWbKjazk3GmOrJkkSI/TR1HO08K4g79Q7wujTGk8cDpz8IvUdwomcxLxfex6AXP+X7ZdnuxGOMOWRYkgihxetzOH3LaHbFpOFpf4Xb4UD7K5CrPqJ51HbGeu7nqTfH8EHmH25HZYyJYJYkQmjidz/S3TsfT8f+Vdcn4mC1OAXPdVNISUpgTPSjfDnuTV7+doXbURlT5bxeL+3ateO4446jZ8+ebN++HYBVq1Zx3HHH/TnfK6+8QocOHdi2bRt33XUXrVq1ok2bNlx44YV/LlOVpk2bxo8//vjndL9+/f4s3lcZB7v8gViSCJHdBUXUWvQuxXiJO+Eat8PZW71j8Fw/legGx/BK9DMsm/wyj3+5eL/d/Y051JTWblqwYAGpqam88MILf5nn7bffZvjw4UyePJnatWvTo0cPFixYwLx58zjqqKN4/PHHqzyufZNEpLMkESJfzfuD3nzDjsPOgKQGbofzV0kN8PT7DE/z7jwV9TI6fSj3fDSPouIStyMzpsqdeOKJrF27dq+2sWPHMmTIECZPnkzdunUBf8mM0jEiOnfuTFZW1l/WVV1LgpfHpSup1V/WzHHUkZ2UdBvgdijli0lCrhiLfnwD9y18n5fn5DBw120M/fvxxEaVXa7AmKB9ORg2zK/adTZoDecMqdCsxcXFTJ06lf79+//Ztnr1agYNGsTs2bNp0KDsf+Jef/31ckuLV8eS4OWxI4kQ2JFXSPMNk8j11cZzxCluh7N/vhjk4leh4/X8n+9zeiz/F/1f/9E63ZlDXmntpjp16rB161Z69Ojx53tpaWkcdthhjB07tsxlH3vsMXw+H1dcUfYNJ6UlwWNiYv5SEnzVqlWAvyT4oEGDaNeuHb169Yr4kuDlsSOJEJg693fOkdnkHvk3Et267TUYHi+c+19ISOOSaf8hde1O+r58PyP7d6duYoRccDeHrgr+x1/VSq9J7Nixg/PPP58XXniBW265BYD4+Hi+/PJLunXrRr169fZKBqNGjeKzzz5j6tSp5Zbero4lwctjRxIhsCnzY+KkgLqd+7gdSsWJwCn3wHnPcKpnDg9svY9+Iybzx9bdbkdmzEGpVasWw4YN46mnnqKw8H9HyGlpaUycOJH77ruPSZMmATBx4kSeeOIJJkyYQHx8/EFt91ArCV4eSxJVLK+gmCM3TSEnqi5yWGQNaF4hHfsjl75JO+8qnt19Lze++DnLNh74A21MJGvfvj1t27Zl9OjRe7U3b96cCRMmcO211zJjxgwGDRrEzp076dGjB+3ateOGG26o9DYPtZLg5QlbqXARORoYE9DUAngQeMtpbwasAi5T1W3OMvcC/YFi4BZVnbS/bURCqfBpC9fQaWwG24+6hEZXjHA1loPy+7cUv385GwoTuEEe4F/X9qZd0xS3ozKHCCsVHrkitlS4qi5V1Xaq2g44HtgNfAwMBqaqaktgqjONiKQDfYBjgbOBESIS8bfcrJ41hXjJJ+34Xm6HcnBanIy332c0iC1ilP6Tf70ymh+Xb3Y7KmNMmLl1uul0YIWqrgZ6A6Oc9lHABc7r3sBoZ4jTlcByoFO4Aw2GqhK/eioFEk3UEaE9BAyLxh3w9p9EraRE3vI8wvNvvsWkhRvcjsoYE0ZuJYk+wPvO6/qquh7Aea7ntDcGAgsLZTltexGRASKSKSKZ2dnuFqxbkZ1LRmEm2XVPOPixqyNF2lF4r5tMbGpj3vA9zgfvjWSs1XsyFWA9+CNPZX4nYU8SIhIN9AI+ONCsZbT95SdU1ZGqmqGqGWlpaVURYqXNmzOL5p6NxB97rqtxVLlaTfD2n4Sv4bG8HPUsP497nle//93tqEwEi42NZcuWLZYoIoiqsmXLFmJjY4Nazo2b+M8BflXVjc70RhFpqKrrRaQhUNrHPAtoGrBcE2BdGOMMWv5vXwNQu83ZLkcSAgl18Pb7lJL3r+CZVS/xr4m5PLX7Fu4486iw3rNtDg1NmjQhKysLt4/uzd5iY2Np0qRJUMu4kSQu53+nmgAmAH2BIc7z+ID290TkGaAR0BKYGcY4g6Kq1Nkyk22+etSu3dztcEIjJgnPlR+gH13HPxe/w/Dvc3lg1x08ekFrvB5LFOZ/oqKiaN68mv4d1DBhPd0kIvFAD2BcQPMQoIeILHPeGwKgqguBscAiYCIwUFWLwxlvMJZt3EmHkoXsqH+Cv2NadeWLQS59E+3Ql5t9n5D+68Pc9n4mBUVWGNCY6iisRxKquhuos0/bFvx3O5U1/2PAY2EI7aAtnv8LvSUHT6tT3Q4l9DxepOdQiK/DFdOfodaSXG4c9SDPX3UicdERf5eyMSYI1uO6iuT99i0AtY8NbbGtiCECZzwEZ/6b870z6Ld6MNe9Mo0du60woDHViSWJKqCq1N3svx4htZu5HU54dbkZLniRrt5F3L3xLq57aRKbdu5xOypjTBWxJFEF1m/Po3XJYrbWPb56X48oT7u/4/nbO7T2/cGQHXdz04hPrTCgMdWEJYkqsPS3JdSX7cQ0O8HtUNzT6lw8V39Ms+gdDMsbzO0jPmTpBisMaMyhzpJEFdi27CcA6h3T1eVIXNasG95rv6BeHIwsup+HXnqHX9dsczsqY8xBsCRRBXzrZ1GIj+jGbd0OxX0N2+K7bjLJScm8xiMMffV1vl9mHaqMOVRZkjhIRcUlNMpdyMb4o8Fno7gBUOcIfNdPIabOYbzieZz3Ro3gi/nr3Y7KGFMJliQO0rIN20lnJfkN2rsdSmRJboSv/0Q8jdryvO85po1+htEz17gdlTEmSJYkDtIfi2cRL/kkHtHZ7VAiT3wqvn4ToMUpPBk1khXjH+elb1e4HZUxJgiWJA7SnjX+kfDSjj4EhyoNh+gEvH8fQ3H6hdwf9R465SGGfLHYqoMac4iwJHGQojcvJE/i8KRaMbNy+aLxXvIaJRn9udH3Kc1+HMy9H/5KYbHVezIm0lmSOAiqStruZWyMOxI8tiv3y+PFc97TaPe76eObxmnz7uKmN6aTm1/kdmTGmP2wb7aDsHbbLo7S1eTXTXc7lEODCHLa/XDuU5zhnc3ANbcx4MUv2ZRjZTyMiVRBJQkRuVREkpzXD4jIOBHpEJrQIt+q5YtJkjxim1j/iKB0uh5Pn3dpHbWWIdvu4NbnP2D5JuudbUwkCvZI4p+qulNEugFnAaOAF6s+rEPD9pWzAUhrmeFyJIegVufiveYLGsUV8WLBYP414g1mrtzqdlTGmH0EmyRKB/05D3hRVccD0RVdWERSRORDEVkiIotF5EQRSRWRKSKyzHmuHTD/vSKyXESWishZQcYaehvmU4yH+Mat3Y7k0NTkeHwDppKQUp9XeJR3XhvK5/Os050xkSTYJLFWRF4GLgO+EJGYINcxFJioqq2AtsBiYDAwVVVbAlOdaUQkHegDHAucDYwQkYga0SYlZwmboppAdLzboRy6UpsTNeArPI3bM8z3HLPH/ItXv1tht8gaEyGCTRKXAZOAs1V1O5AK3FWRBUUkGegOvAagqgXOOnrjP22F83yB87o3MFpV81V1JbAc6BRkvCGTV1DM4UUryal1tNuhHPriU/H1+5TiY3rzQNS7+CYP5tEJ8ykusURhjNuCShKqultVx6nqMmd6vapOruDiLYBs4A0RmS0ir4pIAlBfVdeXrg+o58zfGPgjYPksp20vIjJARDJFJDM7O3yF5Fat20BTyYb6x4Ztm9VaVCzeS99ET7yZfr7JnJh5G4Pe/N5ukTXGZcHe3ZQhIh+LyK8iMk9E5ovIvAou7gM64L+W0R7YhXNqqbzNldH2l38tVXWkqmaoakZaWloFQzl4m1bOByCxyXFh22a15/EgZ/37f7fIrrqZG16YwNrteW5HZkyNFezppneBN4CLgZ7A+c5zRWQBWao6w5n+EH/S2CgiDQGc500B8zcNWL4JsC7IeENm97olAKQ1tyRR5Tpdj+fvY2gVnc0zO/7BvcNHMeeP7W5HZUyNFGySyFbVCaq6UlVXlz4qsqCqbgD+EJHSk/inA4uACUBfp60vMN55PQHoIyIxItIcaAnMDDLekPFsWUoRXqLTjnQ7lOrpqDPxXf8VtZMTeaX4n7wx8lk+mxcx/yMYU2P4gpz/IRF5Ff9dSPmljao6roLL3wy8KyLRwO/ANfgT1VgR6Q+sAS511rlQRMbiTyRFwEBVLS57teGXuHMlm6Ia08gb5XYo1Vf9dKL+7xsK3/87Q9c+x9Njsli56R4Gnd4SqYljiRvjgmCTxDVAKyAKKK3OpkCFkoSqzgHK6nl2ejnzPwY8FmSMIVdcojQoWM3O1JZuh1L9JaYRdc1nFI+/mTvmj+GTb9dxd/bD/PvSDGJ8EXVHtDHVUrBJoq2q1vieY2s353AYG1mWeq7bodQMvhi8F72Mph3NBV8/yuGLb+Kmlx7hiX49qJtoowEaE0rBXpP42enkVqOtW7WQKCkmrtExbodSc4gg3e+Ay96mddRa/pV9C3cMfZsFa3e4HZkx1VqwSaIbMMcpkxHsLbDVxs4/FgFQ5/Aaf1AVfum98PWfSFpiFC8X3strLz3F+Dlr3Y7KmGor2CRxNv67jM4k+Ftgq43iTUsBSG5a4w+q3NGoHVE3fo+3cXue9Q5j/Yd38/jnC6yHtjEhEGyP69VACv7E0BNIqegtsNVJXM4KtnjqQEyS26HUXIn1/Be0M/pzg+8zuv58AwNf/YoduwvdjsyYaiXYHte34u9QV895vCMiN4cisEiWuucPtsYe5nYYxheN9/xnoNdwukYt4b61N3Hb8Hf4baONTWFMVQn2dFN/4ARVfVBVHwQ6A9dXfViRa09hMQ1LNpCfdLjboZhSHa7Ge82XNEgQRuTdw4svPMWkhRvcjsqYaiHYJCH8b0wJnNc1qlfT2o2bqCs5SGpzt0MxgZp2JPrG7/A2bMOznudY8f5dPDNpkV2nMOYgBZsk3gBmiMjDIvIw8DPwepVHFcE2r/HXbIqtb+U4Ik5SA6L7f0FR+37c5JtAxvQBDHxlClty8w+8rDGmTMFeuH4Gf6/rrcA24BpVfTYUgUWq3A0rAEhtauNIRCRfNL7eQ6HnMLpGLeWhdf/HPUNf49c129yOzJhDUrAXrp9Q1V9VdZiqDlXV2SLyRKiCi0TFW34HIKWRleSIaMf3xXvdFOokJ/BS4QN8/spDjPphpY14Z0yQgj3d1KOMtnOqIpBDRdSO1eyQZCQuxe1QzIE0akf0Td+jR/bgn95R1Jl4A3e/9yO7bCAjYyqsQklCRG4UkfnA0U5P69LHSmB+aEOMLMl5f7A1upHbYZiKiqtN1N/fp+T0RzjXO5MblvbntmHvsXxTrtuRGXNIqOiRxHv4O89N4H8d6XoCx6vqFSGKLeKUlCj1itazO8H6SBxSPB48J92Gp9+nNI0vZNiuO3nl+cetnIcxFVChJKGqO1R1lapeDuQA9YHDgeNEpHsoA4wkG7fvpCGbKandzO1QTGU060b0wB/wNG7PE57n2fXhQP75wUzyCiJmmBJjIk6wF66vA74DJgGPOM8PV31YkWnDmmX4pITotCPcDsVUVlIDYq79nOIut/F33zdcOf8abh32rvXSNqYcwV64vhXoCKxW1VOB9kB2RRcWkVVO5dg5IpLptKWKyBQRWeY81w6Y/14RWe5UnT0ryFir3M51ywBIbmh3Nh3SvD68Zz4CV46jeXwew3Lv4N0XHub9Gavt7idj9hFsktijqnsARCRGVZcAwXYYOFVV26lq6Qh1g4GpqtoS/7Cog531pwN9gGPxV58dISKuDkVWtGUlAKlNLUlUC0eeTvSgn5HDu/CI51WSP7uOu9/5np17rEigMaWCTRJZIpICfAJMEZHxwMGOTt8bGOW8HgVcENA+WlXzVXUlsBzodJDbOjg7sijCS0xKY1fDMFUosR4xfT+m5PRHONs3i1uXX8Ndz77GvKztbkdmTEQItsf1haq6XVUfBv4JvMb/vtQrtApgsojMEpEBTlt9VV3vrH89/uqyAI2BPwKWzXLa9iIiA0QkU0Qys7MrfOarUmJ2r2OLpy54bGzlasW5+8nbfzJpSbE8n38vk14ezMhvl1FitZ9MDRfskcSfVPVbVZ2gqgVBLNZVVTvg74A38AB3RpVVOPAvf7GqOlJVM1Q1Iy0tLYhQgpeUv4GcmPoh3YZxUZMMYgb+QMnRPbnL+z7pX/Xl5pc/Y932PLcjM8Y1lU4SlaGq65znTcDH+E8fbRSRhgDO8yZn9iygacDiTTj4U1uVpqrUKdpEXrx1pKvW4lKI7jMK7TmME6J/5z8b/o/nnnucz+a59tEzxlVhSxIikiAiSaWv8Q+BugB/B72+zmx9gfHO6wlAHxGJEZHm+IdNnRmuePe1LTeP+mxFk5u4FYIJFxHk+L5E3fQDMQ1a8SRDKfqgP/e/9z05dlHb1DBBJwkRaSsig5xH2yAWrQ9MF5G5+L/sP1fVicAQoIeILMNfG2oIgKouBMYCi4CJwEBVda3XU/a61fikBF+q9bauMeocQeyAyRSfcj+9vDMYtPRq7n/6BWau3Op2ZMaETdiGL1XV31W1rfM4VlUfc9q3qOrpqtrSed4asMxjqnqEqh6tql8GE2tV277BX/01Ia2Zm2GYcPP68J5yN57rp1A7JYXhhQ8x77WBPPPlXAqKStyOzpiQs+FLKygvezUAKY2st3WN1Ph4Ygf+QEGHa7nO9wXn/nQF/xj2NovX57gdmTEhZcOXVlDJdv/duLXq29jWNVZ0PNG9noUrPqR5/B6ezfkHk0b8gxFfLaao2I4qTPVUFcOXvlblUUUg78617CAJiUlyOxTjtpY9iLl5BiWtenKb9wO6f9eH263+k6mmqmL40udCEFfEic9bx7aoegee0dQMCXWI6fMmXPY2LeN28syO2/ji+dt56Ws7qjDVS0UHHTpSRLoCBA5fCiSKSI04SZ9SsJFdsQ3dDsNEmvRexNzyC8WtenGb9wO6TevDP55/j+Wb7KjCVA8VPZJ4DijrU5/nvFetFRSVUF83U5BoNZtMGRLqENvnDdQ5qnh62y18Nvx2Rn6zxI4qzCGvokmimarO27dRVTOBZlUaUQTK3rKZZNmNp5YlCVM+Se9FzK2ZFLfqzW3eD+j6zWXcMXQUC9bucDs0Yyqtokkidj/vxVVFIJFs2wb/7a++2pYkzAHEp/55VHFE/G6eyfkHP790A099+quNgGcOSRVNEr+IyF/6Q4hIf2BW1YYUeXZtzgIgoU7TA8xpjJ+k9yL21kyK2l7Fdd4v6JN5Kf98+jmmL9vsdmjGBKWiSeI24BoRmSYiTzuPb4Hr8I9WV63lb1sLQK16liRMEOJSiLlwGFzzJam1avFU/qNsfusqHnpvGtt2BVM82Rj3VChJqOpGVe2Cf1zrVc7jEVU9UVU3hC68yFCcsx6wJGEq6fAuxN/yE4Un3UNP30xuX/p3hj71EONnZ9lwqSbiBdtP4htVHe48vg5VUJHGm7uBXOLxxFpHOlNJvhiiTr8P700/Et0wnYd1BGnjLuXulz/i9+xct6MzplwHTBIicpWInCoiH4jI+yJyYzgCiyQxeZvY7q3jdhimOkg7mvgBkyk57zkyYv7gP+sHMGnYQIZ+OZc9hXZh20SeihxJdATOU9VLVfVyoFWIY4o4CQXZ5EbXdTsMU114PHg6XkP0bb9SdOxF3Oj9hIt/vohH//skUxdV+7O35hBTkSSRA9QRketF5BIgIcQxRZzaxZvJj7WSHKaKJdYj7rJXod8X1K6dyn8KhiDv/437X5tA1rbdbkdnDFCxJPFP/KPF1QaigQqNH1EWEfGKyGwR+cyZThWRKSKyzHmuHTDvvSKyXESWishZld3mwcrLL6KubqM4sYFbIZjqrllXEm7+kaIej9Et+jceXHMtnzw7iJenLrAxK4zrKpIkTgLmAW2AC/CffqqsW4HFAdODgamq2hKY6kwjIulAH+BY4GxghIh4D2K7lbY5ez0xUoQn2eo2mRDyRuHrOojoW3+luNX5DPJ8xNnfXsjDTz3NN0s2HXh5Y0KkIknicuAB4HbgSqBSF65FpAlwHvBqQHNvYJTzehT+JFTaPlpV81V1JbAc6FSZ7R6s7Zv840jE1G7kxuZNTZPckPjL34SrJ1A3JZH/7Pk38u4l3PfyByzfZHdBmfCrSJI4FqivqtmqWgBUthDNc8DdQODxc31VXQ/gPJee+G8M/BEwX5bTFnZ/9raua30kTBi1OJmEW2ZQ1OMxToxZwaPrBvDD8Gv57yc/sSOv0O3oTA1S0WsSTwRMTwp2IyJyPrBJVStawqOs0e7K7HUkIgNEJFNEMrOzs4MN7YBKe1un1DusytdtzH75ovF1HUTM7XMpbHc1V3mncN3sS3jpyXt4/6cVFJdYRzwTegdMEqr6rap+FzD9cSW20xXoJSKrgNHAaSLyDrBRRBoCOM+lJ1+zgMB/3ZsA68qJb6SqZqhqRlpaWiVC2z/d6e9tnZTWpMrXbUyFJNQl7sKheG6cTlTjttyjr5HxZU8eemYoM37f4nZ0ppoLdvjSSlHVe1W1iao2w39B+mtVvRKYAPR1ZuuL/y4qnPY+IhIjIs2BlsDMcMS6L1/uBraThETtrxCuMWFQ/1gSr/8c7fMejZK8/Dv3IXa9cREPvjbOem2bkAlLktiPIUAPEVkG9HCmUdWFwFhgETARGKiqrnRHjdmziR0+621tIoQI0uo8Em7PpPC0R+gavYwH1/Tn52FX88QH37I5N9/tCE01IxUpMCYib6vqVSJyqzNsaUTKyMjQzMzMKl3non91gugE0u/5pkrXa0yV2LWZ3V89TszsN9ijPkZxPp4ut9D31NbERbty17g5BInILFXNKOu9ih5JHC8ihwPXikhtpxPcn4+qCzXyJBdvoyDGSnKYCJVQl/jeT+O9+Re05ZncJB9x8Y89Gf7EPYydYRe3zcGraJJ4Cf9pn1b4BxkKfFTtv+4RpKi4hFTdQVG8JQkT4eocQeKV78B1XxPT4BjuLn6FjM/P5bGnhvDN4o1WktxUWkXHkximqscAr6tqC1VtHvBoEeIYXbN1xzbiJR9PQtXfNWVMSDQ5nuQbJqGXjyEtJYkHdw+h1vvn8vDwl5m5cqvb0ZlDULDjSdwoIm1FZJDzaBOqwCJBTrb/rltfcn2XIzEmCCLI0WeTdNsMis4fTqu4HTyy9R7yX+/Jwy+OYl7WdrcjNIeQoJKEiNwCvIu/Z3Q94F0RqXTBv0iXu8WfJGJSrG6TOQR5vPgyrib+zvkUnPFvMmLX8vDGW9j48oX8+9Ux/LZxp9sRmkNAsLfAXgecoKoPquqDQGfg+qoPKzLs2e6v7Z+QahVgzSEsKo7objcTd+cC8rvfz0nRy3ggawDLnr+YIW99wpotVpbclC/YJCFAYH+FYsouoVEtFOb4O4An13WlbJQxVSsmkZjT7ib2zgXknfgPzoiaz10r+pH53GX8d/RE1m7PcztCE4F8Qc7/BjBDREpLc1wAvFalEUUQzfUniaQ6diRhqpG4FOLOegi6DWTX10/T89dXYfGPfLyoO6uPuYE+Z3WnaWq821GaCBHshetngGuArcA24BpVfS4EcUUE7+5sckhAfDFuh2JM1UuoS0LPx4m6fR757fpyofcHbl9yOb88exn/ffdTVm/Z5XaEJgIEeySBqv4K/BqCWCJO1J4t5Hhqk+x2IMaEUnJDEi98Fs64h9xvnqHn7Dfx/jadz5d0ZkzLG7j03B40r1vjRi02DrdrN0W0+IKt5EbVPvCMxlQHSQ1I7PUkUXcsJK/TQM6KmsPdv/dlydAL+O8oG/SopqpwkhC/GjXyTlLxNvJjrLifqWES00g47zGi71zErs63c2r0Qu5aeR2rhvfkidfesX4WNUyFk4T6+/V/ErpQIouqUqtkO0VxVpLD1FDxqSSc/TCxdy5iV5e76Rq9nHv+GMiul8/hieeHM/23bCv3UQMEe7rpZxHpGJJIIkzOrt2kyC7USnKYmi4uhYQz7yfu7kXsOe1RWsdv5p7ND1DnndN4+ul/88WcNVZIsBoLNkmcCvwkIitEZJ6IzBeReaEIzG3bsv3DlnqT6h1gTmNqiJgkYrvfSuJdCyk8/3kaJfu4M/cp2nx8Ki8MuZMxPyxhT6Erw76YEAr27qZzQhJFBMrd4u9tHVPL+kgYsxdfNFEZV1GrwxUUL51I/Ff/5ZYtr7J18mjemnouvs4DuKhbW1Lio92O1FSBoJKEqq6u7IZEJBb4Dohxtvuhqj7kjEcxBmgGrAIuU9VtzjL3Av3x9+y+RVUnVXb7wcrb5h/bOq621W0ypkweD95jziX1mHPR1T+hU55iQNZY8n74hPHTu7MxvR89Tz+VFmmJbkdqDkKwBf5ERK4UkQed6cNEpFMFF88HTlPVtkA74GwR6QwMBqaqaktgqjONiKTjHw/7WOBsYISIhG2orYIdGwFIrmNJwpgDkcNPpM51H8FNM8g/5iIu9n7HrUuuZM2wc3n6xRH8uMwuch+qgr0mMQI4Ebjcmd4JvFCRBdWv9EbrKOehQG9glNM+Cn+pD5z20aqar6orgeVARRPSQSvZtRmwkhzGBKVeK1L6vEzUnYvZ1fVeMmLXcsfGe0l7+2SGP3kf42b8Rn6RXbc4lASbJE5Q1YHAHgDntFCFTzyKiFdE5gCbgCmqOgOor6rrnfWtx1+CHKAx8EfA4llO277rHCAimSKSmZ2dHeSPs59Y87ZSqF6i42tV2TqNqTES6pLQYzCJ9yymoPfL1Kldi1vyRnDaF6fw3mP9ee3z6WTvzHc7SlMBwSaJQueUjwKISBpQUtGFVbVYVdsBTYBOInLcfmYvq7rsX45XVXWkqmaoakZaWtXdrurN28oOTzJItS1ya0zo+aKJbt+H1Nt+RK+ZSOHhJ9FXx9N3Zk9m/rcXw157g8yVW+xUVAQLNkkMAz4G6ovIY8B04D/BblRVtwPT8F9r2CgiDQGc503ObFlAYA/vJsC6YLdVWVEF28n1WNUmY6qECHL4iaRdOwbPrXPIbX89p0Ut5JY/biP5jZMY8eRgxk5fxO6CIrcjNfsItgrsu8Dd+BPDOuACVf2gIsuKSJqIpDiv44AzgCXABKCvM1tfYLzzegLQR0RiRKQ50BKYGUy8ByOucDt5vpRwbc6YmqP24aRc8CRx9/xG/nnDqFM7hYF5L3HelFP49D99eHHMeH7PtjpRkSKoW2Cd21jPBU7Cf5opWkRWquqeCizeEBjlnK7yAGNV9TMR+QkYKyL9gTXApQCqulBExgKLgCJgoKqG7YpXfPEOdsQdEa7NGVPzRMcT07EvMR37olmzyJv2IheuGE/04in8svAoJta9iCNPuYLTjm2Cz2u1SN0iwZwLdL60dwLvOE2XA7VV9dIQxBa0jIwMzczMrJJ1bXmoKSvTTidj0JtVsj5jTAXs3krujLcomvEqKXv+YLMm84X3NPLbXMFZJ3XjsDo2GFIoiMgsVc0o671ge1wf7fRzKPWNiMytfGiRKb+wkBR2onFWJtyYsIpPJfHU2+DkWyhe/jUl00ZwxboJeOd8woxfW/F5nZ4c3u1yTm/bjBhf2LpN1WjBJonZItJZVX8GEJETgB+qPix37di6mXqiSIKVCTfGFR4P3qPOoN5RZ8DODeT8/BZHzxrFCdv+S86EF/jk05PIOeZyTjnlDFrWT3I72mqtQklCRObjv/00CrhaRNY404fjv2ZQrezctol6gC/JKsAa47qkBiT3uBvOuIuSld+T992rXLj6S6IXT2LBwmaMrHUe9bpcQY8OR5MQE/Rgm+YAKrpHzw9pFBFm13Z/SY7oJBtLwpiIIYKnRXfqt+gOedvIzXyfBjPeZMDOF9gzcSRTJp7A2sMvoE23XnQ+sh4ej/VxqgoVShIHU9jvUJS/w99zOz7FjiSMiUhxtUk86SYST7oJXTebnG9f4Yzl44lbM50N7z7Ge76TKTzuMrp3O5kjrMDgQQn2FtgM4H78p5l8+HtFq6q2CUFsrinc6a/blFi7vsuRGGMORBq1p97lI6DwGQoWf4H++DZ9NnyKb+4nLJjdjNeSepDc6e/06HiclS+vhGBP4L0L3AXMJ4hyHIeakl1bAEhOteJ+xhwyomKJbnMRDdtcBLs2szNzNPUy36X/zlco+vo1pk9ty4pG53N4l0s4Kb2p3R1VQcEmiWxVnRCSSCKI7N5KgfqIjreyHMYckhLqknTyIJJOHoRuWsK2H9+i/cIPOGXDv9j50VNMGteRTYf35JguPencsgFeu35RrmCTxEMi8ir+cR/+LOGoquOqNCqXefdsJceTRF0r7mfMIU/qtSLtgv9Ar39T9Pt35P7wFj1WTyJu9XdsXfUIn3i6kHNkb9p0OZsOzeog9ne/l2CTxDVAK/y3wpaeblKgWiWJqIJt5HpqYfc2GVONeDz4jjyFhkeeAkX5FCyZzJ6f3+P8tV8Ts2wyG36rzVhfN/KPuZCOXc6gVcNkSxgEnyTaqmrrkEQSQWILd7DbivsZU335Yog+rieNjusJBbvYvfBzimeM5qINE4la8Clr5qUxOvZkitIvpNMJ3TmqQVKNTRjBJomfRSRdVatdB7pACcXb2RbX0u0wjDHhEJ1AfPvLiG9/GeRtJ3fueLy/jOayLePwzvmQVb/WZ3RMV4qPPo92nc/g2Ma1alTCCDZJdAP6icjv+K9JVMtbYJNLctgcY3WbjKlx4lJI7NyXxM59YddmcuZ8TPTscVy6eTy+BeNYNz+VD6NOJP/I80jvfBbtDqtT7TvtBZskzg5JFBGkoLCIWuRSEpfqdijGGDcl1CW56/Ukd73e38N7/ueUzPqQ3hsnE73kc7IXJzPBewI7W5zD0Z3P4/gW9avlXVLBJom+5bQ/erCBRIod27JJE0XiLUkYYxxxtUnsdCWJna6E/Fx2L/ySPZkfcvb6acQun8KOZQ/wpSeDrY1Pp+Hx59MlvVm1qSMV7E+xK+B1LP6aTourLhz37dqxmTTAm2BJwhhThphE4jtcSnyHS6Ewj7ylX5Ez4wNOXfsNCVnfUfDHo8z8OJ1VdbqT2LYXJ3ZoR/3kWLejrrSgkoSqPh04LSJP4R9m9IBEpCnwFtAA/+2zI1V1qIikAmOAZsAq4DJV3eYscy/QHygGblHVScHEWxl5O/wlOaIS7ZqEMeYAouKIO64nTY/rCcVFFK35meyZH9Pq98l02/YCTHuBxV835duELnDUObQ+4VRaNTy0Lnwf7PFQPNCigvMWAXeo6q8ikgTMEpEpQD9gqqoOEZHBwGDgHhFJB/oAxwKNgK9E5KhQD2Gav3MrADGJNpaEMSYIXh++5t1o3LwbALp5GdmzxpO86HMu3vEB3rljyJ5Ti898Hdl52Bk07nAWnY4+nLjoyC4PEmyBv9JxJQC8QBrwr4osq6rrgfXO650ishhoDPQGTnFmGwVMA+5x2keraj6wUkSWA52An4KJOVgFu/xJIi7ZkoQxpvKkbkvqnXUnnHUn7N7KjvlfsHvOBE7f8D3xK7+i4Pf7maWtWFW7CzHHnEXbDp1pkZYYcUcZwR5JBI4rUQRsVNWiYDcqIs2A9sAMoL6TQFDV9SJSz5mtMfBzwGJZTtu+6xoADAA47LDDgg3lL4pztwGQkGL9rY0xVSQ+lVonXEmtE66EogIKfp/Oxtmf02Ll15y44xX4+RXW/ZTKZ77j2XXYyTRsfw6dWjWPiKOMio5M9+B+3lNVrdDRhDN/IvARcJuq5uwna5b1hv6lQXUkMBIgIyPjL+8HqyRvOwDJliSMMaHgiyb6qNNoetRp/ukda9k69wsKF07k9E0/EL9yCkW/P8hsjmJlyolEH30W6R260LK+O2VCKnoksauMtnjgOqAOFTzlJCJR+BPEuwFFATeKSEPnKKIhsMlpzwKaBizeBFhXwXgrTfZsZ49GERubEOpNGWMM1GpMavfrSe1+PRQXkr96JhtnfUqTlV/TcccbMPMNsmckM8nblm0NupCSfgbHt2lDvTDdMVXRken+vKvJueh8K3AtMBp4urzlAok/Bb4GLFbVZwLemoC//8UQ53l8QPt7IvIM/gvXLYGZFdnWwfDmbydXEjl0b1gzxhyyvFHEtOjKYS26+qdzN7F13hfsWjCZLht/JHnd97DuCVZMbsjPse3Ja3ISjdr1oMPRzUPWL6PCa3VuVf0HcAX+C8wdSm9VraCuwFXAfBGZ47Tdhz85jBWR/sAa4FIAVV0oImOBRfivfwwM9Z1NAL6CHHZ5kqwCrDHGfYn1SO3Sj9Qu/UCVko2L2DhnIr5lX3Pm1q+J/f0LilfcxwJtwea2N3L6xddXeQgVvSbxX+Ai/Of+W6tqbrAbUtXplH2dAeD0cpZ5DHgs2G0djJiiHPJ8NiauMSbCiOBpcCwNzz4Wzr4DigrIX/MLG2dPJHXlNBJSokKy2YoeSdyBv6DfA8D9ARdPSgv8VZsh3GKLctgVa8OWGmMinC9671NTodpMRWZSVU9Io4gg8SW5bI+uNjnPGGMOSo358q+oJM2lJKaW22EYY0xEsCQRIL8gnyTJQ2NT3A7FGGMigiWJADnb/MX9PHEp7gZijDERwpJEgN1OBViPjSVhjDGAJYm97M7ZAkBUoiUJY4wBSxJ7KS0THptkFWCNMQYsSeylINfKhBtjTCBLEgGKnbEkElIsSRhjDFiS2Iv+WSY8zd1AjDEmQliSCCB7tpOn0fhi4twOxRhjIoIliQCe/B3slCS3wzDGmIhhSSJAVEEOuzxWAdYYY0pZkggQXZTDHq8lCWOMKRW2JCEir4vIJhFZENCWKiJTRGSZ81w74L17RWS5iCwVkbPCEWNs8S4KfHa6yRhjSoXzSOJN4Ox92gYDU1W1JTDVmUZE0oE+wLHOMiNExBvqAONKcimKsiRhjDGlwpYkVPU7YOs+zb3xD4WK83xBQPtoVc1X1ZXAcqBTqGOM190UR1uSMMaYUm5fk6ivqusBnOd6Tntj4I+A+bKctr8QkQEikikimdnZ2ZUOREtKSNTdlMTYgEPGGFPK7SRRnrLGwtayZlTVkaqaoaoZaWmV7wS3e3cuUVIMliSMMeZPbieJjSLSEMB53uS0ZwFNA+ZrAqwLZSC7dmwDwBNno9IZY0wpt5PEBKCv87ovMD6gvY+IxIhIc6AlMDOUgex2KsB6LUkYY8yffOHakIi8D5wC1BWRLOAhYAgwVkT6A2uASwFUdaGIjAUWAUXAQFUtDmV8ebn+I4moBEsSxhhTKmxJQlUvL+et08uZ/zHgsdBFtLcCJ0lEJ9Q+wJzGGFNzuH26KWIU7t4BQGxiiruBGGNMBLEk4Sh2kkRcsh1JGGNMKUsSjpI8f5JISLbxrY0xppQliVL5OZSokJBoF66NMaaUJQmH5OeQK3GIJ+Qloowx5pBhScLhKdjJbhLcDsMYYyKKJQlHVOFO8jyWJIwxJpAlCUdUUS57vJYkjDEmkCUJR2xxLgU+G5XOGGMCWZJwxJbsotAGHDLGmL1YknAk6C6KLUkYY8xeLEngH3AoQXejMZYkjDEmkCUJYE/ebqKlGGKsI50xxgSyJAHsyvGPJSFxNiqdMcYEsiTB/5KEDThkjDF7syQB5O/aDtiAQ8YYs6+IThIicraILBWR5SIyOFTbybcBh4wxpkwRmyRExAu8AJwDpAOXi0h6KLZV6BxJ2IBDxhizt4hNEkAnYLmq/q6qBcBooHcoNlRUOuBQko0lYYwxgSI5STQG/giYznLa9iIiA0QkU0Qys7OzK7WhuLqH8WvCSSTVTqtcpMYYU0353A5gP6SMNv1Lg+pIYCRARkbGX96viNYnXwQnX1SZRY0xplqL5COJLKBpwHQTYJ1LsRhjTI0UyUniF6CliDQXkWigDzDB5ZiMMaZGidjTTapaJCKDgEmAF3hdVRe6HJYxxtQoEZskAFT1C+ALt+MwxpiaKpJPNxljjHGZJQljjDHlsiRhjDGmXJYkjDHGlEtUK9X/LCKJSDawupKL1wU2V2E4VcXiCk6kxgWRG5vFFZzqGNfhqlpmyYlqlSQOhohkqmqG23Hsy+IKTqTGBZEbm8UVnJoWl51uMsYYUy5LEsYYY8plSeJ/RrodQDksruBEalwQubFZXMGpUXHZNQljjDHlsiMJY4wx5bIkYYwxplyWJAAROVtElorIchEZ7GIcTUXkGxFZLCILReRWp/1hEVkrInOcx7kuxLZKROY728902lJFZIqILHOea4c5pqMD9skcEckRkdvc2F8i8rqIbBKRBQFt5e4fEbnX+bwtFZGzwhzXf0VkiYjME5GPRSTFaW8mInkB++2lMMdV7u/N5f01JiCmVSIyx2kP5/4q77sh9J8xVa3RD/xlyFcALYBoYC6Q7lIsDYEOzusk4DcgHXgYuNPl/bQKqLtP25PAYOf1YOAJl3+PG4DD3dhfQHegA7DgQPvH+Z3OBWKA5s7nzxvGuM4EfM7rJwLiahY4nwv7q8zfm9v7a5/3nwYedGF/lffdEPLPmB1JQCdguar+rqoFwGigtxuBqOp6Vf3Veb0TWEwZ43pHkN7AKOf1KOAC90LhdGCFqla2x/1BUdXvgK37NJe3f3oDo1U1X1VXAsvxfw7DEpeqTlbVImfyZ/yjPoZVOfurPK7ur1IiIsBlwPuh2Pb+7Oe7IeSfMUsS/h39R8B0FhHwxSwizYD2wAynaZBzeuD1cJ/WcSgwWURmicgAp62+qq4H/4cYqOdCXKX6sPcfr9v7C8rfP5H0mbsW+DJgurmIzBaRb0XkJBfiKev3Fin76yRgo6ouC2gL+/7a57sh5J8xSxIgZbS5el+wiCQCHwG3qWoO8CJwBNAOWI//kDfcuqpqB+AcYKCIdHchhjKJf3jbXsAHTlMk7K/9iYjPnIjcDxQB7zpN64HDVLU98A/gPRFJDmNI5f3eImJ/AZez9z8iYd9fZXw3lDtrGW2V2meWJPwZtmnAdBNgnUuxICJR+D8E76rqOABV3aiqxapaArxCiA6190dV1znPm4CPnRg2ikhDJ+6GwKZwx+U4B/hVVTc6Mbq+vxzl7R/XP3Mi0hc4H7hCnZPYzqmJLc7rWfjPYx8Vrpj283uLhP3lAy4CxpS2hXt/lfXdQBg+Y5Yk4BegpYg0d/4j7QNMcCMQ55zna8BiVX0moL1hwGwXAgv2XTbEcSWISFLpa/wXPhfg3099ndn6AuPDGVeAvf7Dc3t/BShv/0wA+ohIjIg0B1oCM8MVlIicDdwD9FLV3QHtaSLidV63cOL6PYxxlfd7c3V/Oc4AlqhqVmlDOPdXed8NhOMzFo4r85H+AM7Ff7fACuB+F+Pohv+QcB4wx3mcC7wNzHfaJwANwxxXC/x3SswFFpbuI6AOMBVY5jynurDP4oEtQK2AtrDvL/xJaj1QiP+/uP772z/A/c7nbSlwTpjjWo7/fHXpZ+wlZ96Lnd/vXOBXoGeY4yr39+bm/nLa3wRu2GfecO6v8r4bQv4Zs7IcxhhjymWnm4wxxpTLkoQxxphyWZIwxhhTLksSxhhjymVJwhhjTLksSZgaR0TqBFTu3BBQeTRXREaEaJu3icjVZbQ3C6w4WgXbiRaR75zOX8YcNPsgmRpH/b1k24G/PDWQq6pPhWp7zhf2tfiri4aUqhaIyFTgb/yv3IYxlWZHEsY4ROQUEfnMef2wiIwSkcnOGAIXiciT4h9TY6JTIgEROd4p7jZLRCbt02u41Gn4y4YUBSwzV0R+AgYGbL+ZiHwvIr86jy5O+9si0jtgvndFpJeIHCsiM52joHki0tKZ5RPgilDsI1PzWJIwpnxHAOfhL7v8DvCNqrYG8oDznEQxHLhEVY8HXgceK2M9XYFZAdNvALeo6on7zLcJ6KH+Qop/A4Y57a8C1wCISC2gC/AFcAMwVFXbARn4ewiDv5xFx0r+zMbsxU43GVO+L1W1UETm4x/UaKLTPh//gDNHA8cBU/yldfDiL+mwr4b46/+XfsmnqOq3zntv4y9QCBAFPC8i7YBinGJxqvqtiLwgIvXwF5n7SFWLnCOR+0WkCTBOnRLWqlosIgUikqT+sQeMqTRLEsaULx9AVUtEpFD/V8OmBP/fjgALyzgi2FceEOu8Fsov2Xw7sBFoi/8of0/Ae2/jP4XUB//1DVT1PRGZgf9oZ5KIXKeqXzvzx+yzvDGVYqebjKm8pUCaiJwI/lLOInJsGfMtBo4EUNXtwA4R6ea8F3jtoBawXv2lsq/Cf2RS6k3gNmcdC53ttQB+V9Vh+AvitXHa6wDZqlp48D+iqeksSRhTSeof7vYS4AkRmYu/MmeXMmb9Ev/YyaWuAV5wThflBbSPAPqKyM/4TzXtCtjWRvzJ5o2A+f8GLBCROUAr4C2n/VT81yyMOWhWBdaYMBCRj4G7de+hL4NZPh7/tZAOqrrjAPOOA+5V1aWV2ZYxgexIwpjwGIz/AnbQROQMYAkwvAIJIhr4xBKEqSp2JGGMMaZcdiRhjDGmXJYkjDHGlMuShDHGmHJZkjDGGFMuSxLGGGPK9f8CO9olj3jYQAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# plotting the result\n", "plt.plot(t,nCoEu,label='Euler method')\n", "plt.plot(t,nCoRK2,label='RK2 method')\n", "plt.legend()\n", "plt.xlabel('Time (days)')\n", "plt.ylabel(r'Number of $^{56}$Co atoms')" ] }, { "cell_type": "code", "execution_count": 55, "id": "aa06e471", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0, 0.5, 'Number of atoms')" ] }, "execution_count": 55, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAEGCAYAAACUzrmNAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABKAklEQVR4nO3dd3iUZfbw8e9J7w1CAoTeQXoRkI4o2FBRZN1VXEUEe/mpqOvqa1l7r4uiK4KCIgIqKIiAIL0joYcWSgihBFIg5X7/uCcQ4qRnMinnc11zzeSeZ+Y5PBly5u5ijEEppZQqiIe7A1BKKVXxabJQSilVKE0WSimlCqXJQimlVKE0WSillCqUl7sDcJWaNWuahg0bujsMpZSqVNasWXPUGBOZt7zKJouGDRuyevVqd4ehlFKViojsdVauzVBKKaUKpclCKaVUoTRZKKWUKpQmC6WUUoXSZKGUUqpQLksWIvKZiBwRkT9zlUWIyDwR2eG4D8/13BMislNEtonI5bnKO4vIJsdz74qIuCpmpZRSzrmyZvE/YHCesnHAfGNMM2C+42dEpDUwAmjjeM2HIuLpeM1HwGigmeOW9z2VUkq5mMuShTHmd+BYnuKhwBeOx18A1+Yqn2KMOWOM2Q3sBLqJSG0gxBizzNi11Cfmeo1LfLF0D7M2HHTlKZRSqtIp7z6LKGPMIQDHfS1HeV1gf67j4h1ldR2P85Y7JSKjRWS1iKxOTEwsUYBfr9zHD5oslFLqAhWlg9tZP4QpoNwpY8x4Y0wXY0yXyMi/zFYvkrAAb06kni3Ra5VSqqoq72SR4GhawnF/xFEeD9TLdVwMcNBRHuOk3GXCA3w4kZrhylMopVSlU97JYhYw0vF4JDAzV/kIEfEVkUbYjuyVjqaqUyLS3TEK6tZcr3GJsABvjmuyUEqpC7hsIUER+RroB9QUkXjgGeBl4BsRuQPYB9wIYIzZLCLfALFAJnCPMSbL8VZjsSOr/IE5jpvLhAX4cCL1LMYYdJSuUkpZLksWxpi/5fPUwHyOfxF40Un5auCiMgytQOEB3mRmG1LOZhHkW2UX5VVKqWKpKB3cFUaYvw8Ax1O0k1sppXJossgjLMAbgJNp2m+hlFI5NFnkERbgqFno8FmllDpHk0Ue4Y6ahQ6fVUqp8zRZ5BF6LllozUIppXJossgjp4NbaxZKKXWeJos8fLw8CPTx1Il5SimViyYLJ8ICfDiRps1QSimVQ5OFE3YxQa1ZKKVUDk0WToQH+OjQWaWUykWThROhAd6c1JqFUkqdo8nCifAAb61ZKKVULposnAjz9+FkWgbZ2fnus6SUUtWKJgsnwgK8yTZwKj3T3aEopVSFoMnCiXDH+lA6fFYppSxNFk7krDyrE/OUUsrSZOFEzsqzuj6UUkpZmiycCNOVZ5VS6gKaLJwI15qFUkpdQJOFEyF+XojAMa1ZKKUUoMnCKS9PD0L9vTmWcsbdoSilVIWgySIfNQJ9OJaizVBKKQWaLPJVI9CXpNOaLJRSCjRZ5CtCaxZKKXWOJot81AjyIUmThVJKAZos8lUj0O5pkaWLCSqllCaL/EQE+mCMzrVQSinQZJGviCBfAO23UEopNFnkq2agncV9VEdEKaUqiePpx1l1eJVL3tvLJe9aBUQE2WShNQulVEWUkZXB9uPb2ZC4gY1HN7IxcSP7T+1HEJbdvIxA78AyPZ8mi3xEBOYkC53FrZRyL2MMCakJNjEk2sSw5dgWzmTZv0+R/pG0i2zHsGbDaBfZDh8PnzKPwS3JQkQeAkYBBtgE/BMIAKYCDYE9wHBjzHHH8U8AdwBZwP3GmF9cHWPOYoI6fFYpVd7SMtOITYo9lxg2Jm7kSNoRAHw8fGhdozU3tbiJdpHtaB/ZnqiAKETEpTGVe7IQkbrA/UBrY0yaiHwDjABaA/ONMS+LyDhgHPC4iLR2PN8GqAP8KiLNjTFZrozT29ODsABvncWtlHK54+nHWXdkHeuOrGPtkbXEJsWSmW23da4XXI+utbvStmZb2ke2p0V4C7w9vcs9Rnc1Q3kB/iKSga1RHASeAPo5nv8CWAg8DgwFphhjzgC7RWQn0A1Y5uog3TaLO2EzbP8FDqyBk/sh7QR4+YJvCNRsBrVaQcPeULsDeOgYBaUqE2MMB04fYN2RdaxJWMO6I+uIOxkHgLeHNxfVvIiRrUfSoVYH2kW2I8Ivws0RW+WeLIwxB0TkdWAfkAbMNcbMFZEoY8whxzGHRKSW4yV1geW53iLeUfYXIjIaGA1Qv379UsdaI9CHpPLqszAGts2GRa/AoQ22LKIJ1GgCNVtAdgakJkHcQtjwtX3ePwLaXAvtb4aYLuDiaqhSqviysrPYcWIHaxPW2ppDwtpzTUrB3sF0qNWBq5tcTadanWhTsw2+nr5ujtg5dzRDhWNrC42AE8C3IvKPgl7ipMzptGpjzHhgPECXLl1KPfW6RqAvcUdPl/ZtCnfyAMy82yaCGk1hyGvQ5joIinR+/KkE2L3I1j7Wfw2rP4NabaDnfXDRMPAq+84tpVTRZJtsdhzfwarDq1h5eCVrEtaQfDYZgKiAKDpHd6ZTrU50iupE07CmeEjlaB1wRzPUpcBuY0wigIhMB3oCCSJS21GrqA0ccRwfD9TL9foYbLOVy0UE+bBqj4ubofYth69HQOYZuOJ16PxP8Czk1xIcBe2G21t6Mmz+HpZ/BDPGwG/PQ++HodNIcEO7plLVjTGGuJNxrDy8klWHV7H68GqOnzkOQExQDJc2uJQuUV3oHNWZOkF13BxtybkjWewDuotIALYZaiCwGkgBRgIvO+5nOo6fBXwlIm9iO7ibASvLI9Cc9aGysw0eHi5o4tk+F765BULqws3fQM2mxX8PvxDoPBI63Qo75sHiN+CnR2zyGPgMtLpam6eUKkPGGPYm7z2XHFYdXkVSehIA0YHR9I7pTdfornSL7lapk0Ne7uizWCEi04C1QCawDtt0FAR8IyJ3YBPKjY7jNztGTMU6jr/H1SOhckQE+pBt4ERaxrl5F2Vm/yr45laIbAH/mA6BNUv3fiLQ/DJoNgi2/wy/PmsTUf2ecNWbtlNcKVUix9OPs+LQCpYeXMrSg0tJSE0A7PyGi2tfTLfobnSL7kZMcIzLh7C6i1tGQxljngGeyVN8BlvLcHb8i8CLro4rrxrn1oc6U7bJ4mQ8fH0TBEfD378rfaLITQRaDIGmg2D9JJs0Pu4FPe6Bvo+DT9nO6lSqKsrIymB94nqWHVzG0oNLiU2KxWAI9gmme+3udK/dna7RXWkY0rDKJoe8dAZ3AWrkWh+qaa1CDi6q7CyYPtr2Udw+Lf9O7NLy9ILOt0HLq+HXf8Mf78Cf0+Hqt6Hppa45p1KVlDGG3cm7zyWHVYdXkZaZhqd40j6yPXd3uJuedXrSukZrvDyq55/N6vmvLqIajvWhynRi3h9vw94/4NqPS9ZHUVyBNWDoB9DhH/DDAzBpmO38vuwF29+hVDWVmpHKskPLWBy/mD8O/sHhlMMANAhpwNAmQ+lZpyddo7sS5BPk5kgrBk0WBYh0NEMdPV1Gcy2O74FFr0LrodB+RNm8Z1E16AF3/Q4LXoSl78Gu32Do+9C4X/nGoZQb7U3ey+/xv7M4fjGrE1aTkZ1BkHcQPer0YHS70fSo3YOY4Bh3h1khabIoQHiAD14ewpFT6WXzhj8/CeIJl7/knhFK3n5w2fPQ8iqYMRYmDoWud9oyb//yj0cpFzubdZY1CWtsgjiwmL3JewFoHNqYv7f6O31i+tChVge8PXSYeWE0WRTAw0OoGeRL4qkyqFnsXgzbfoJLn4VQpxPQy0/9i2HMEjsnY/mHsGcxDJsA0Re5Ny6lysDRtKP8Hv87i/YvYvmh5aRmpuLj4UPX2l25ueXN9Inpo7WHEtBkUYjIYF+OlEWyWPgyBEXDxWNL/15lwScABr8ETQfC92PhkwG2htFttM7LUJXOvuR9/LbvN37b/xvrj6zHYIgKiOLKxlfSJ6YP3aK7EeAd4O4wKzVNFoWIDPYlIbmUzVC7F8PeJTDkVdsUVJE0vRTGLrXLjcx5DHbOh2s/LNvhvEqVMWMMsUmxzN83nwX7F7DzxE4AWkW0YmyHsQyoN4Dm4c2rzbDW8qDJohC1gn3ZdOBk6d5k8eu2VtFpZNkEVdaCIu0M8pXjYe7T8FFPuPYjW+tQqoLIyM5g9eHV52oQR1KP4CmedI7qzLhu4+hfr3+VmjFd0WiyKERksC9Jp8+QlW3wLMmSH4nb7AKBA56ueLWK3ETg4rugwSXw3R0w6Xro9RD0/1fha1Up5SIZWRksO7SMX/b8woL9Czh19hR+nn5cUvcSBtQfQJ+6fQjzC3N3mNWC/hUoRGSwL9nG7sUdGVyCpYNXfgKePhW3VpFX9EVw5wL4eRwsecsudDhsgvs75VW1cUGC2LeAUxmnCPYJpn+9/lxa/1K61+mOv5eO3itvmiwKUcuRII6cSi9+skhPtntPXDTMdTO1XcEnAK55Fxr2gh8etMuFXD/erjullAs4q0EEewczoP4ALmt4GT1q93DL7nDqvGIlC8deFPWMMRtdFE+Fk5MgSjR8NnYGnD0NXe4o26DKS7vhUKcjfHsbTL4BLnkQBvxLlz5XZSIjO4NlB/+aIPrX78/lDS/XBFHBFJosRGQhcI3j2PVAoogsMsY87NrQKoZawbafoUTDZzd+Y3e7i+lSxlGVo5rNYNSv8PMTdqmSfcvghs8gVMepq+IzxrAhcQM/xv3I3D1zOX7muCaISqIoNYtQY0yyiIwCPjfGPCMi1aZmUTOohDWLE/vtZLf+T1X+eQve/nYBwoa97PpSH/eC6/4LzS93d2Sqkth1Yhc/xf3E7N2zOXD6AL6evvSv158rGl3BJXUvwcdTd3es6IqSLLwcO9cNB55ycTwVjr+PJ8G+XsVPFpu+tfdtbyz7oNyl7Q1Qu4NtlvpqOPS8Hwb+W5ullFMJKQnM2T2Hn3b/xNZjW/EQD7rX7s7dHe5mYP2BBHrrcvmVSVGSxXPAL8ASY8wqEWkM7HBtWBVLZEgJlvzYPB1iukFEI9cE5S41m9pmqV+egKXv2tFSN3wGYfUKf62q8lIyUpi7Zy4/xP3A6sOrMRja1mzL410fZ3CjwdT018melVWhycIY8y3wba6f44Bhrgyqooks7vpQx/fC4U0w6HnXBeVO3n5w1Vu2WWrWA/Df3nYSX4sh7o5MuUG2yWZNwhpm7JzBvL3zSMtMo0FIA8a2H8sVja+gQUgDd4eoykBROrgbAfcBDXMfb4y5xnVhVSy1QvzYFH+i6C/YNsfet7zSJfFUGBcNczRLjYSvR0CPe+2+317a/lwdHDx9kJm7ZjJr5yziT8cT6B3IFY2u4Nqm19I+sr0utVHFFKUZagYwAfgByHZpNBVUsWsW236CyJZQo4nrgqooajSBO36FX56EZe/D3qW2WaqqNb8pANIy0/h176/M3DWTlYdWYjBcHH3xuX4IXayv6ipKskg3xrzr8kgqsMhgX1LOZpFyJpNA30IuWeox2PMHXPJA+QRXEXj7wVVvQuO+MPM++G8fuPoduOh6d0emysjmpM1M2z6Nn3f/zOmM09QNqsvYDmO5psk11A3S2f3VQVGSxTsi8gwwFzj39doYs9ZlUVUw52dxn6FRYckibgGYLGhxRTlEVsG0HmqbpabdDtP+CbsXweCXdWOlSiolI4XZu2fz7bZv2XJsC36eflzW8DKubXotnaM64yEe7g5RlaOiJIu2wC3AAM43QxnHz9VCdKidmJeQnE6jmoUM94tbCL6hULeT6wOriMIbwO0/242V/ngH9q+EGz6HWi3dHZkqotikWL7d/i2z42aTmplKs/BmPHnxk1zZ+EpCfHTf9uqqKMniOqCxMeasq4OpqHKSxeGThexrYQzsWgiNeoOHp+sDq6g8vWHQc9CwD3x/F4zvB1e8Ch1vqfwTFKuo1IxUZu+ezbTt09ictBk/Tz8ub3g5NzS/QTurFVC0ZLEBCAOOuDaUiis6xCaLQ4Uli2NxcHIf9KpG/RUFaXYpjP0Dpt8Js+6DuEV2yK2ffjutKHYe38mUbVP4Me5HUjJSaBrWlHHdxnF1k6urZS0iIyOD+Ph40tNLueFZJeDn50dMTAze3kWbVFuUZBEFbBWRVVzYZ1Fths4G+noR4ufF4ZNpBR8Yt9DeN+7v8pgqjeBouGUGLH4TFv4HDqyBGz+3CxQqt8jMzmTR/kV8tfUrVh5eiY+HD4MbDebG5jdW+1pEfHw8wcHBNGzYsEpfB2MMSUlJxMfH06hR0UYuFiVZPFO6sKqG2qH+hdcs4hZCaD2IaFwuMVUaHp7Q91FoeAl8Nwo+HWSbqbqP1WapcnQi/QTf7fiOqdumcijlENGB0TzQ6QGGNRtGuF+4u8OrENLT06t8ogAQEWrUqEFiYmKRX1OUGdyLRCQK6OooWmmMqXZNUlGhfhwuaC9uY+wcg2aX6R/A/DToCWOWwMx77HIhcQth6PsQVMvdkVVpW5K28PXWr5m9ezZnss7QLbobj3V9jH71+uHloVva5FXVE0WO4v47izKDezjwGrAQEOA9EXnUGDOtJAFWVrVD/NhyKDn/A47FQepRqH9x+QVVGQVEwIivYMV/Yd6/7X7fQz/QFWzLWGZ2JvP3zeerLV+x9sha/L38uabJNfyt5d9oFt7M3eGpSqgoXyueArrm1CZEJBL4FahWySI61I+jp89wNjMbHy8n48v3Lbf39bqXb2CVkQh0HwON+tjO76+GQ9dRdi0tH50BXBqpGal8v/N7voz9kgOnDxATFMP/dfk/rm16LaG+oe4OT1ViRUkWHnmanZKAajcbp3aoH8bY7VVjwp38Qdu/HPzCoGbzco+t0opqDaPm2zkZy96H3Yth2CdQu727I6t0ElIS+GrrV3y7/VtOnT1Fx1odebTLo/Sr1w/P6jyMW5WZoiSLn0XkF+Brx883AXNKc1IRCQM+BS7CTvC7HdgGTMUuWLgHGG6MOe44/gngDiALuN8Y80tpzl8SuedaOE8WK6FeN/Codnm0dLz94PIXoelA+H4sfDIQBj4NPe7Ta1kE245tY2LsRGbvnk22yWZg/YGMbDOS9pGacCurhQsX8vTTT9OmTRtGjBhBnz59ePrpp0lOTqZLly6MHDnSLXEV+r/RGPMo8F+gHdAeGG+MeayU530H+NkY09LxnluAccB8Y0wzYL7jZ0SkNTACaAMMBj4UkXL/qlQ71C5Z4XREVOoxSNwK9bS/osSaDIC7l0GLwbYvY+I1cDLe3VFVSMYY/jjwB3fOvZMbfriBeXvnMbz5cH687kfe7PemJopKTkQICgoiPT2dmJgYZs6cyYEDB/D29iYmJuaC4x555JFzP7/++us8++yzAPTs2bPM4ypKB/crxpjHgelOyopNREKAPsBtAI6Z4WdFZCjQz3HYF9gO9ceBocAUY8wZYLeI7AS6ActKcv6SKnAWd/wqe19f+ytKJSAChn8J6ybBnMdt5/dVb+uChA5Z2VnM2zePCZsmsPXYViL9I3mg0wPc2PxG7Y+oQnr37k3fvn1JSEjg4Ycfpm3btvTo0YO77rqLG264gYEDBwLg6+vL9OnTeeKJJ6hZ88JNpZYuXVrmcRWlnj/ISVlpdrlpDCQCn4vIOhH5VEQCgShjzCEAx33OeMq6wP5cr493lP2FiIwWkdUisro444eLIsTPiwAfT+c1i/jVIB5Qp5quB1WWRKDTLTBmMdRoZhck/PaftvZWTWVkZfD9ju8ZOnMojy56lPTMdJ7r+Rw/D/uZUW1HaaKoYjwcza/h4eGcOXOGmJgYwsPtPBhPz/ONKl5eXowePZq33nrrL+8RFBRU5nHlW7MQkbHA3UBjEdmY66lg4I9SnrMTcJ8xZoWIvIOjySm/UJyUGWcHGmPGA+MBunTp4vSYkhIRokP9SHA21+LQert/hY7kKTs1msDtv8Afb8HCV2DPErvsecvqs5pvWmYa03dM5/M/PychNYFWEa14o+8bDKw/UDutq7Dp06fzyy+/cOLECe699166devGfffdx+LFi+nTp88Fx95zzz20a9eOxx4rbc9A4QpqhvoK25H9Ehf+MT9ljCnN17x4IN4Ys8Lx8zTH+yeISG1jzCERqc35tajigdwbPMcAB0tx/hKLDvHjkLMlPw5tgCYDyz+gqs7TC/o8Cs0Hw/djYMrfoP3NMPgl8A9zd3Quk3w2malbpzJpyySOpR+jU61OPNvzWS6pc0m1mTBWEfy/HzYTe7CAuVUl0LpOCM9c3abAY66//nquv/7CptcJEyY4PTYkJIRbb72Vd999F39/124FkG8zlDHmpDFmjzHmb8aYvUAa9ht9kIjUL+kJjTGHgf0i0sJRNBCIBWYBOd38I4GZjsezgBEi4uvY4rUZsLKk5y+N6FC/v/ZZJB+C0wk63NOVotvCnQts4tg4FT7sATt/dXdUZe7kmZO8t+49Lp92Oe+ue5fWNVrzv8H/44shX9Crbi9NFMqpBx98kAkTJpCSkuLS8xSlg/tq4E2gDvbbfgPs6KWC02PB7gMmi4gPEAf8E5u4vhGRO4B9wI0AxpjNIvINNqFkAvcYY7JKce4SqxPqT8KpM2RmZePl6cizhzY4nuzgjpCqDy8fGPAvaDHEDrGdNAw63waXvQC+we6OrlROnjnJxNiJTN4ymZSMFAY1GMSdbe+kVY1W7g6tWiusBlBRREREMHz4cCZMmMDtt9/usvMUZZ7FC0B34FdjTEcR6Q/8rTQnNcasB7o4ecppW44x5kXgxdKcsyzEhPuTlW04nJxrrsWh9YBA1EXuDK36qNsZ7vodFrwAS9+HXb/Z5UIa9Sn8tRVM8tlkvoz9kkmxkzidcZpBDQYxpv0YmofrxE5VPI888gjvv/++S89RlGSRYYxJEhEPEfEwxiwQkVdcGlUFlZMg4o+n5UoWG6BmM/At+9EHKh/efrZG0eJKmDEWvrgaOo20K9lWgr6M5LPJTIqdxKTYSZzKOMWl9S9lTPsxtIhoUfiLlXI4ffr0ucdRUVGkpqY6fa6sFCVZnBCRIOB3bNPREWxzULUTE247kOKP5+rkPrgeGvZyT0DVXYMeMHap3Sdj2QewYy5c+Qa0vNLdkTl16uwpJm2ZxJexX3Lq7CkG1h/I2PZjNUmoSqEoyWIotnP7IeDvQCjwnCuDqqhqh/khAvHHHRn8dCKcOqid2+7kE2BrGW2ug1n3w5Sb7eMhr1aYpc/TM9OZsnUKn/75KSfPnGRAvQGM7TCWlhG6L7mqPIqyn0VOF3s2dmZ1teXr5UlUsN/5msWRWHsfVTk6wqq0up1h9EL4421Y9CrsWmCH2Lb/m9v2F8nMzmTGzhl8tOEjjqQe4ZK6l3B/x/tpXaO1W+JRqjR055Niign3P1+zOLLF3tfS//wVgqe3HV7b6hpby5gxFjZ9a5cMCW9QbmFkm2zm7Z3H++veZ0/yHtpHtufl3i/TNbpr4S9WqoLSZFFMMeH+rN573P6QuAX8wytMc4dyiGwB/5wDqyfAr8/Ch92h/5Nw8RibUFzEGMOyg8t4e+3bbDm2haZhTXm3/7v0q9dP50ioSi/fSXkiMt9xXy1HPuUnJjyAQyfTyczKtjWLyFa6jWpF5OEB3e6Eu5dDw94w918wvp9dSt4FNiVuYtTcUdz1612cPHOSF3u9yLSrp9G/fn9NFKpKKKhmUVtE+gLXiMgU8qzRZIxZ69LIKqhzcy1OphFzZCu0vcHdIamChNWDm6fC1h/tSrYTBtlhtpc+a1e5LaWDpw/y9tq3mbN7DhF+EYzrNo4bm9+Ij6dP6WNXqgIpKFn8G7tmUwx2BnduBhjgqqAqspz5FQkH9xBz5iTU0lm2FZ4ItLoaGveDhS/D8o9s8rjshRJ3gJ86e4pPN33KpNhJiAh3tr2T2y+6nSAfnW+jqqZ8k4UxZhowTUSeNsY8X44xVWg5cy1S92+yBZosKg/fYLsrX/sR8ONDtgN83WQ7N6NW0YaxZmZn8t327/hww4ccSz/GVY2v4oFODxAdGO3i4JVyr6IMnX1eRK7BblgEsNAY86Nrw6q4cuZaZOeMhIrUZFHpRLeF2+fCuokw7xn4+BLoeR/0/r98Z+IbY1h8YDFvrH6DuJNxdI7qzIcDP6RNTR02rcpWpd1WVUReAh7ALuQXCzzgKKuWcuZa+B/fDoG1ILCGu0NSJeHhYRcivHc1tL0RlrwF73eFTdPAXLgVyvbj2xk9bzT3zL+HLJPF2/3f5vPLP9dEoVyiqNuqHj58mBEjRtCkSRNat27NFVdcwfbt210WV1GGzl4JdDDGZAOIyBfAOuAJl0VVwcWE+xN+LA5idAZupRcUCdd9bBPHnMfguztg1acw5FVORjTgg/UfMHXbVIJ9ghnXbRzDmw/H24XDb5Uqyraqxhiuu+46Ro4cyZQpUwBYv349CQkJNG/umoUoizrPIgzI2fCo2u/hWD/cn6iE/VDjEneHospK/e52z4x1X5I1/zm+/+oK3o2sxUmTxY0tbuS+jvfp9qWqXDjbVtXHx46uy9lWdcGCBXh7ezNmzJhzr+vQocO5x2+++SafffYZAKNGjeLBBx8sdVxFSRYvAetEZAF2+GwfqnGtAqBFaCYhpJAR3hj9jlmFeHiyvl4HXmrZidjj2+iUdponT2XQomUt8Ap0d3SqvM0ZB4c3le17RreFIS8XeEhRtlX9888/6dy5s9PXr1mzhs8//5wVK1ZgjOHiiy+mb9++dOzYsVShF6WD+2sRWQh0xSaLxx273VVbrXztjq8JXjHEFHKsqhwSUxN5e+3bzNo1i1oBtXil9ysMCWyI/DwO5jwKaz6HwS9D477uDlVVccXZVtWZJUuWcN111xEYGHju/RYvXuz6ZAFgjDmE3d5UAY04BEBcdrQmi0ouIzuDybGT+Xjjx5zNOsuotqO4s+2dBHg79iu5dSZs+QHmPgUTr7H7gQ96zi4poqq2QmoA7tSmTRumTZvm9DmTZ4BGWSl0NJT6q8iMA2QaD7alh7k7FFUK646sY/gPw3ljzRt0ierCjKEzeKDTA+cTBdgJe62vgXtW2Vnfe5faPcB/fNguUa+UGwwYMIAzZ87wySefnCtbtWoVixYtok+fPsyYMYPU1FRSUlL4/vvv6d27d6nPqQsJloDfyTj2ShRxx8+6OxRVAifPnOStNW/x3Y7viA6M5t3+79K/fv+CX+TtB70ego63wKJXYNUE2PgN9H4Iut8N3v7lE7xS2OG133//PQ8++CAvv/wyfn5+NGzYkLfffptmzZpx22230a1bN8B2cJe2CQpACqqyiIgHsNEYU+k2mO7SpYtZvXq1a978416sOebHazWfZ8roHq45hypzxhh+iPuB11e9TvLZZG5pfQtj24+9sCZRVEd32Al9236CkBgY+DS0HW7nb6hKa8uWLbRqVX0m2jr794rIGmNMl7zHFvjJdsyt2CAi9cs2xErMGEiK43RQQ/YmpRZ+vKoQ4k7GccfcO3hqyVPUD6nP1Kum8kiXR0qWKMDuu/63r2DkjxBYE76/Cz7pB3ELyzJspSqMojRD1QY2i8hKIGfXPIwx17gsqors1GHISCE7vDGHDqaTdjYLfx9Pd0el8pGemc6nmz5lwp8T8Pfy5989/s2wZsPwkDKqATTqbedn/DkN5j8HE4dCo75w6TN29z6lqoiiJIv/5/IoKpOknQD4RjeHzbD3WAoto0PcHJRyZvXh1Ty77Fn2Ju/lqsZX8UiXR6jpX7PsT+ThAe2G2x36Vn8Gi1+HTwbYlW4HPK0jp1SVUJR5FotEpAHQzBjzq4gEANX3q/SxXQCEx7QE9rLnaKomiwomJSOFt9a8xdRtU6kbVJfxg8bTo0459C15+0GPu6HTLbDsA1j6Pmz9CdrfDP3G2b01lKqkirKQ4J3ANOC/jqK6wAwXxlSxHdsNHt7UbdAUgD1JKYW8QJWnJQeWcO3Ma/lm2zf8o9U/mH7N9PJJFLn5Btvk8MAGO1Jq07fwXic7I1iH26pKqijNUPcA3YAVAMaYHSJSfTedPrEPwuoREuBHjUAfdidqsqgITqSf4LXVrzFr1ywahzZm4pCJdKjVwb1BBdaw+2d0H2uH264cD+u+tHuB97inTHbqU6q8FKWX74wx5tyEAhHxwu6UVz2d2AthDQBoEhnErsTTbg5Izd0zl6EzhzI7bjaj243m26u/dX+iyC00Bq55D+5ZAc0GweI34O128NsLkHqs8NcrVQEUJVksEpEnAX8RGQR8C/zg2rAqsBP7IMyOJG4aFcSOI6ddNr1eFSwpLYmHFjzEI4seISogiilXTeG+jvdV3P2vazaDG/8HY5dCs0vh99c1aahKoyjJYhyQCGwC7gJmA/9yZVAV1tkUSEmEcFuzaFYriJNpGSSePuPmwKqfX/f+ynUzr2NR/CIe6PQAX135FS0iKsmoo6jWeZLGa5o0VIVXlNFQ2Y4Nj1Zgm5+2mer6VfrEfnsflpMsggHYmXCaWsF+7oqqWkk+m8zLK17mh7gfaBXRiv/0+g9Nw5u6O6ySyUkafR6zfRq/vwbLP4buY2zHuPZpVEuVeVvVK4FdwLvA+8BOERni6sAqpBN77b0jWTSPsvs17zii/RblYenBpVw/83pm757NmPZjmHzl5MqbKHKLag3Dv7A1jaYDHTWNtvDLU5B8yN3RqXJW1G1VPT096dChw7nbnj17XBpXUUZDvQH0N8bsBBCRJsBPwJzSnFhEPIHVwAFjzFUiEgFMBRoCe4DhxpjjjmOfAO4AsoD7jTG/lObcJXZin7139FlEBvsS4ufFjiOn3BJOdZGakcqba95k6rapNAptxKT+k7ioZqVbrqxwUW1s0kiItXuCL//QjqDqcDNc8gBENHZ3hKocFGVbVQB/f3/Wr19fbnEVpc/iSE6icIgDjpTBuR8AtuT6eRww3xjTDJjv+BkRaQ2MANoAg4EPHYmm/B3fA15+EGRHDosIzaKC2Z6gNQtXWX9kPTf+cCPfbPuGW1vfyjdXfVM1E0VuUa1h2Cdw31ro+A9Y/zW81xmm3Q6H/3R3dMrFnG2rGh4eDpzfVrUgkyZNolu3bnTo0IG77rqLrKysMokr35qFiORs1bRZRGYD32D7LG4EVpXmpCISA1wJvAg87CgeCvRzPP4CWAg87iifYow5A+wWkZ3YeR/LShNDieSMhBI5V9SsVhBzYxPKPZSqLiM7g483fMynmz6ldmBtJlw+ga7RXd0dVvmKaARXvQV9H7czwld/Bn9+B80uh94P233Dlcu8svIVth7bWqbv2TKiJY93e7zAY4qyrSpAWlrauX23GzVqxPfff8+WLVuYOnUqf/zxB97e3tx9991MnjyZW2+9tdSxF9QMdXWuxwlAzn6SiUB4Kc/7NvAYEJyrLMqxIx/GmEO5Jv7VBZbnOi7eUfYXIjIaGA1Qv74LFso9sfdcE1SOprWCmLJqP0mnz1AjyLfsz1kN7U/ez7jF49h4dCPXNr2Wcd3GEehdjffADo6Gy563CWLlp7DiI/jscqjfE3reZ3fv06XRq4yibqvqrBlq/vz5rFmzhq5d7RertLQ0atUqmznU+SYLY8w/y+QMeYjIVdimrTUi0q8oL3FS5nQ0ljFmPDAe7H4WJY0xXyf2Qd0Ll3lvFmXz3Y4jpzVZlJIxhh/jfuSF5S/gKZ681vc1Bjcc7O6wKg7/cOj7qF1/au1EW9uY8jeIaGLL2t8MPiVccl39RWE1gIrIGMPIkSN56aWXyvy9izIaqpGIvCki00VkVs6tFOe8BLhGRPYAU4ABIjIJSBCR2o5z1uZ8v0g8kHsFthjgYCnOXzLpyZB2/C81i2a1HCOiErSTuzROnT3F44sf58klT9IyoiXfXfOdJor8+ATaJUTuXw/DJoBfCPz0CLzVxs7VOKXNotXVwIEDmTZtGkeO2D+fx44dY+/evWXy3kUZDTUDmICdtZ1d2hMaY54AngBw1Cz+zxjzDxF5DRgJvOy4n+l4ySzgKxF5E6gDNANWljaOYjsZb+/zrBxaO9SPYD8vthzWZFFS64+sZ9zicRxOOcy9He5lVNtReHpU34WNi8zTC9reABcNg33L7Cq3v78Of7xjd+3rcY/tLFfVRuvWrXnhhRe47LLLyM7Oxtvbmw8++IAGDRqU+r2LkizSjTHvlvpMhXsZ+EZE7gD2YTvSMcZsFpFvgFggE7jHGFM23fvFkeyozITEXFAsIrSuHULsweRyD6myy8zO5JONn/Dxxo+pHVibL4Z8QfvI9u4Oq/IRgQY97S1plx1yu24yrJ8ETQZA93vsvfZrVCmnTzsfhXnTTTdx0003lfn5ipIs3hGRZ4C5wLl1LYwxa0t7cmPMQuyoJ4wxScDAfI57ETtyyn2SHTWLkDp/eapNnVC+WrmXrGyDp4ezLhaV15HUIzz2+2OsSVjD1Y2v5smLnyTIJ8jdYVV+NZrAlW9A/6fs6KmV42HyMNuv0e1OO2fDL9TdUapKqCjJoi1wCzCA881QxvFz9ZF8EMTDjkzJo3WdENIzstl9NIWmtfQPXmGWHljKE0ueIC0zjRd7vcg1TarnDr0uFRABff7PjpaKnWWTxs/jYP7z0H6ETRy1Wrk7SlWJFCVZXAc0zr1MebWUfACCosDT+y9Pta5td8qLPZSsyaIAmdmZfLThIz7Z+AlNwprwRt83aByms5JdyssX2t1obwfXwcpPYN0kWD0BGvWBbqOh+RDb/6FUAYrSiLkBCHNxHBVf8kGnTVBg51p4e4r2WxQgMTWRO+feyfiN47m26bV8deVXmijKW52OcO2H8PAWuPRZu+vj1H/AO+3tHhu6ix9AtdlyoLj/zqJ8nYgCtorIKi7ss6hebQcnD0Bkc6dP+Xh50DwqmM0HT5ZzUJXD0oNLeWKxNjtVGIE1oNdD0OM+2P6zbaKa/xwseAlaXgmdb4NGfatlh7ifnx9JSUnUqFEDkarb/2iMISkpCT+/oq+WXZRk8UzJQ6pCkg9Ck/75Pt26dgi/bT2CMaZKf8iKIys7i482fMT4jeNpHNqYzy7/jCZhTdwdlsrh6QWtrrK3xG2w5gvY8BXEzoDwRtB5JHT4+7m10KqDmJgY4uPjSUys+rUsPz+/C1axLUxR9rNYVKqIqoL0ZDh7CkKcrjIC2E7ub9fEk3jqDLVCdG+L4+nHefT3R1lxaAXXNr2WJ7o9QYC3zi6usCJbwOD/wMB/w5YfYM3n8OuzdpLfudpGvypf2/D29qZRo0buDqNCKjRZiMgpzi+v4QN4AynGmBBXBlahnJtj4bzPAuzwWYA/D55kQDVPFpuPbuahhQ+RlJbEcz2f47pm17k7JFVU3n7nO8QTt8PaL2D9ZIidafdxyaltOBkVqKq2Qr8mGGOCjTEhjpsfMAy7CVL1kXzA3hdQs2hTJwQPgfX7q3e/xfQd07l1jl3hcuIVEzVRVGaRzeHyF+HhrXZZkdB6tm/jzVYw+UbY/D1kpLs7SlVOij1ezhgzQ0TGuSKYCutcssi/ZhHo60WL6BDW7TteTkFVLGezzvLSypeYtn0a3Wt359U+rxLuV9rFiVWF4O1nlxVpewMc3WlrGhumwLe3gV+YLe9wM9TpdMHy/apqKUozVO61cj2ALuSz6muVldMMFVy7wMM61Avjp40Hyc42eFSjmdyHUw7z8MKH2XR0E6PajuLeDvfq2k5VVc2mcOkzMOBfELcQ1n9l522s+hQiW9qk0e4mbaaqgopSs8i9r0UmdsvToS6JpqJKPgCBtcDLp8DDOtYL4+uV+9idlEKTyOoxOW/loZU8+vujpGem81a/t7i0waXuDkmVBw9Pu19404GQdsI2Sa3/Cub9G379f9D0UjtTvMUQ8PZ3d7SqDBRlNJRL9rWoVJIPQmj+/RU5OtQPA2DdvhNVPlkYY/gy9kveWPMGDUIa8Pngz2kcqpPsqiX/MOjyT3s7usMmjQ1TYMcv4BMMra62TVWN+upM8UqsoG1V/13A64wx5nkXxFMxJR+0484L0TQyiGBfL9bvP84NnYs+frmyOZt1lueWPcfMXTMZWH8gL/Z6sXrvZKfOq9nsfDPVnsWw6Vu7NtWGr2zt/KLr7fLpdbV/o7IpKM2nOCkLBO4AagDVJ1mcTijSfsceHkK7eqGs33/C9TG5ydG0ozy44EE2JG5gTPsxjG0/Fg+p2mPvVQl4eELjfvZ2xRuwY65NHKs/hxUf2y9fbW+EdsNtglEVXkHbqr6R81hEgoEHgH9id7d7I7/XVTmZZyE1CYKK1mHXoV4YHy+KI+1sFv4+VauTNzYplvt/u5+TZ07yet/Xubzh5e4OSVUG3n7Q+hp7SzsBW3+Ejd/A76/B769C7fbQ5npocy2EN3RzsCo/BX4lFJEIEXkB2IhNLJ2MMY8bY44U9LoqJcXxTy3ikgedG4STlW1Yt79qDaH9ec/PjJwzEhFh4pCJmihUyfiHQcd/wMhZ8MhWuPwlEE/49Rm7oOF/+8KSt+wih6pCKajP4jXgemA80NYY43xbpqrutGM/4yIOBezSMAIPgeVxx+jZpKYLAysf2SabD9d/yH83/pcOkR14q/9b1PSv/P8uVQEER0OPu+3t+F47Szx2hl1m5NdnIbqdrW20vtZu6qTcqqA+i0ewq8z+C3gq1+J4gu3grh7LfZxyJIsi1ixC/LxpUyeUFXFJLgyqfKRmpPLkkieZv28+1zW9jn91/xc+ngUPH1aqRMIbwCX329uJfTZxbJ5hZ4zPfw6i29qk0fpaO9dDlbuC+iy01xLO1yyK2GcBcHGjCCYu30t6RhZ+3pWz3+JwymHu++0+th/fzuNdH+fvrf6uq+mq8hFW3+7w1/M+OLEftsyyieO35+0tspVd3LDlFVC7Y5Vf3LCi0KtcmNPFq1kAXNy4BmczsyvtqKgtSVv4+09/Z/+p/Xww8AP+0fofmiiUe4TVgx73wKh58NBmGPwyBNa0/RqfDIC32sCPD8POX+1gFOUyOkOmMKcTIKCG0+1U89OtYQQisCLuGN0b13BhcGVv0f5FPPr7o4T6hjJxyESahzvf8EmpchcaA93H2lvqMTscd+uPsOFru02sb4idOd7ySmg2CPxC3R1xlaLJojCnEorVBAUQGuBNy+gQVuxOAirPGPLJWybz6qpXaRnRkvcHvE9kQKS7Q1LKuYAIu5xI+xGQkQZxi2DbT7BtDmyeDh7e0LAXtLjCJo4I3aOitDRZFOZ0Qol2CuvRuAaTV1SOfous7CxeW/0ak7dMpn+9/rzc+2XdqEhVHt7+0GKwvWVnQfxqmzi2/gRzHoU5QM3m0Owymzjq9yx0nTf1V5osCnM6oUQzTPs0r8lnf+xmeVwS/VpU3G0pUzNSeez3x1gUv4hbWt/CI50f0RVjVeXl4Qn1L7a3Qc9B0i7YMc+uU7VyPCx7H3yC7MzyZoNsAilg6wF1niaLghhT4ppF98Y18PXyYNH2xAqbLBJSErjvt/vYdnwbT138FCNajnB3SEqVrRpN7K37GDibArt/t30d2x39HQBRbc8njpguxeqfrE40WRQk7ThknS12nwWAn7cn3RvXYNG2xAsXea8gdh7fyZhfx3Dq7CneG/AefWL6uDskpVzLJ9Aumd5iiP0ieGSLTRw75sEf78CSN+0quY16Q5MB0Li/TTQ6EhDQZFGw08Vb6iOvfi0i+X8/xLIvKZX6NSpOH8CahDXc99t9+Hr68sWQL2gZ0dLdISlVvkQgqrW99XrQrlkVtxDiFsCu32DbbHtcaD3bZNVkgL0PiHBbyO6myaIgpw/b+xLu+tW3uR1NtGj7EW7p0bCMgiqdX/f+yuO/P06doDp8POhj6gYVvk+HUlWef5hdWqTNtbbWcSzOkTgW2CXW130JiF30sEl/W+uo3x28fN0bdznSZFGQU8WfvZ1bo5qB1I8IYMG2xAqRLL7e+jUvrXiJtpFt+WDAB4T5hbk7JKUqHpHzfR1dR0FWJhxcZ2sccQtg6Xt2UqCXv+1Ib9gLGvaBOh2r9CgrTRYFKcHs7dxEhIGtajF5xT5On8kkyNc9l9sYw3vr3uOTTZ/QL6Yfr/Z9FX8v3epSqSLx9IJ6Xe2t3+OQngx7/7C1jj1L4LcX7HHeAba2cS55dKhSneXl/tdLROoBE4FoIBsYb4x5R0QigKlAQ+w+38ONMccdr3kCu+lSFnC/MeaXcgk25Qh4+YFvcInf4oq2tfn8jz38tvUI17Qv/yF6GdkZPLv0WWbtmsWwZsP4V/d/4eWh3xGUKjG/kPMd5QApSbB3iU0ce5bYhQ/BDtHNnTxqt6/U28q6I/JM4BFjzFrHpkprRGQecBsw3xjzsoiMA8YBj4tIa2AE0AaoA/wqIs2NMVkujzQlCQJqlmo0ROf64UQG+/Lzn4fKPVmkZqTy8KKH+ePAH9zd/m7GtB+jazwpVdYCa0DrofYGcDrR1jz2LLbJ49dnbblPMNTrBvV72CRStzP4VJyBL4Up92RhjDkEHHI8PiUiW4C6wFCgn+OwL4CFwOOO8inGmDPAbhHZCXQDlrk82NSjdtGyUvDwEC5vE8V3aw6U6+55SWlJ3DP/HrYc28IzPZ7hhuY3lMt5lar2giLPd5aDHVWZU+vYtxwWvAgY8PCytY2c5FGvu31tBeXWOpGINAQ6AiuAKEciwRhzSERyOgrqAstzvSzeUebs/UYDowHq169f+gBTEkudLACGXFSbScv3sWh7IoMvKllneXEcOH2A0XNHcyT1CO/0f4d+9fq5/JxKqXwE1YKLrrc3sPO39q+Cfcts8lj5iZ1ZDhDR5HzyqN+jQs3zcFuyEJEg4DvgQWNMcgHNI86eMM4ONMaMx+7sR5cuXZweUywpSVCzRanf5uJGEUQE+vDjxoMuTxa7Tuxi9LzRpGWm8clln9ChVgeXnk8pVUz+4dD8MnsDyDwDB9fD/uU2eWz7CdZPss8F1ISYrhDT2d7X6WT7TNzALclCRLyxiWKyMWa6ozhBRGo7ahW1gZx9vuOBerleHgMcLJdAy6AZCsDL04Or29Xm61X7OZmWQai/a0ZIbErcxNj5Y/H28Obzyz+nRUTpE51SysW8fM+vZ3XJA5CdDUk7ztc84lfD9jmOgwUiW55PHnW7QK1Wdk0sV4fp8jPkIbYKMQHYYox5M9dTs4CRwMuO+5m5yr8SkTexHdzNgJUuD/RsCmSklkmyABjWOYYvlu3lx40H+fvFDcrkPXNbcWgF9/92P+F+4Xwy6BPqhdQr/EVKqYrHwwMiW9hb59tsWdpxOLAG4tfAgdV2Rd11jtqHT5Cd4xHT5XwCCY4q87DcUbO4BLgF2CQi6x1lT2KTxDcicgewD7gRwBizWUS+AWKxI6nuKZ+RUEftfUDZJIu2dUNpHhXEd2viyzxZzN83n0cXPUqDkAb8d9B/qRVQMRcuVEqVkH+43dip6aX255xZ5vGrbfKIX2UnC2Zn2ucf213mS5O4YzTUEpz3QwAMzOc1LwIvuiwoZ3KSRWDZjE4QEYZ1iuGlOVvZlXiaJpFBZfK+M3bO4Jmlz3BRzYv4cOCHhPrq7mBKVXm5Z5m3v8mWZaTBoY2QuMUla1jpHtz5Sc1JFmVTswC4rmNdPD2Eqav2l8n7Tdw8kaf/eJqLoy/mk0GfaKJQqjrzdiw/ktN0VcY0WeTnXDNU2e2hXSvEj8Ftopmych+pZzNL/D45y3e8tvo1BjUYxPsD39ed7ZRSLqXJIj8pifa+jJqhcozs2ZDk9ExmrCvZgK5sk82LK15k/MbxXN/sel7r8xo+nlV38TKlVMWgySI/qUftulA+gWX6tl0bhtOqdghfLN2DMcWbCpKRncETi59g6rap/LPNP3m2x7O6BapSqlxosshPSpKtVZTx7EkR4baeDdiWcIolO48W+XVpmWk8uOBBZu+ezQOdHuDhLg/rOk9KqXKjySI/KYll2l+R29AOdakV7Mt7v+0s0vGnzp5izLwxLI5fzNPdn2ZU21EuiUsppfKjySI/ZTR72xk/b0/u6tuElbuPsSIuqcBjj6Yd5fZfbmfj0Y282udVhrcY7pKYlFKqIJos8pNytMw7t3O7uVt9agb5FFi7OHT6ELf9fBt7Tu7hvQHvMbjRYJfFo5RSBdFkkZ+Uoy5rhgLw9/FkdJ/GLNl5lKVO+i7iTsZxy5xbOJZ2jPGXjadX3V4ui0UppQqjycKZsymQmeayZqgct/ZoSN0wf174aQtZ2edHRsUmxXLbnNvIyM7g88Gf07FWR5fGoZRShdFk4UzOHIsyWhcqP37enjw2uAWxh5KZvjYegNWHV3PHL3fg5+XHxCETdeVYpVSFoMnCmdRj9t7FNQuAa9rXoUO9MF79ZRs/xy1gzK9jiAyIZOKQiTQIKfvVaZVSqiQ0WTiT5kgW/uEuP5WI8PzQizghK3hs8UM0CWvC/wb/j+hA1++op5RSRaXJwpnU4/bev+xXbnRmS8ov+NWdSkZKfUY3e5UIv/I5r1JKFZUmC2fSHMnCBcv85vXppk95fvnz9KrTm8iUe3h25i6S0zNcfl6llCoOTRbO5DRD+YW57BTGGN5c8ybvrH2HKxpdwbsD3+at4V05eCKdcd9tLPa6UUop5UqaLJxJPQa+oeDpmr2hMrIzePqPp/n8z8+5qcVNvNT7Jbw9vOncIIJHL2/B7E2Hmbhsr0vOrZRSJeGObVUrvrTj4B/mkrdOzUjlkUWPsOTAEu5ufzdj2o+5YEHA0b0bs3L3MZ7/MZYmkUH0aub6EVlKKVUYrVk4k3bMJf0Vx9KPcccvd7D04FKe6fEMYzuM/cvKsR4ewtsjOtC0VhBjJ61h6+HkMo9DKaWKS5OFM2nHy3wk1P5T+7l1zq3sOLGDt/u9zQ3Nb8j32BA/bz67rSv+Pp7cOmEluxJPl2ksSilVXJosnEk9VqZzLGKTYrll9i2cOHOCTy/7lP71+xf6mjph/nx5x8VkZRv+Nn45O49owlBKuY8mC2fKsBlqwb4F3Pbzbfh4+jBxyEQ61OpQ5Ne2iA7m69HdyTaGEeOXsXbf8TKJSSmlikuTRV7ZWZB+stTNUMYYPvvzMx5Y8ABNw5oy+YrJNA5tXOz3aR4VzNS7ehDg48WI8cuZuf5AqeJSSqmS0GSRV9oJe1+KZqizWWd5+o+neWvNW1ze8HI+u/wzIgNKvjdGk8ggZtxzCR1iwnhgynqe+n4TaWezSvx+SilVXJos8sqZkFfCZqijaUe5c+6dzNw1k7Htx/Jqn1fx8/IrdVgRgT5MGnUxo/s0ZvKKfVz9/hJtllJKlRtNFnnlLPVRgprFmoQ1DP9hOLFJsbza51Xu7nD3X4bGloaPlwdPXtGKSXdczOn0TK7/cCmPTdvA0dNnyuwcSinljCaLvHKWJy9Gn4Uxhv/9+T/u+OUOArwDmHzlZIY0GuKiAKFXs5rMf6Qvd/VtzPS1B+j76gJe+Xkrx1LOuuycSqnqTWdw53VuEcGi1SyOph3l2aXPsih+EYMaDOK5ns8R5BPkwgCtQF8vnhjSihs71+Od+Tv4eNEuvli6h2GdYvh79/q0jA5xeQxKqepDk0VexdjLYt7eeTy37DnSMtMY120cN7e8uUybnYqiaa0g3vtbR+4f0JSPFu1i6ur9fLl8L50bhHN1u9oMvqg20aGl7zNRSlVvUlVXN+3SpYtZvXp18V84/3lY8iY8nQQezlvpDp0+xGurX2Pe3nm0qdGG//T6D43Dij8s1hWOp5zlu7XxfLs6nm0JpwDoVD+M3s0i6dmkBh3rh+Pjpa2PSinnRGSNMaZL3nKtWeSVdtzWKpwkitSMVL6M/ZJPN30KwP0d7+e2i27D28O7vKPMV3igD6N6N2ZU78bsPHKan/88xLzYBN77bQfvzN+Bn7cH7eqG0aZuCG3rhtK2bigNagRqAlFKFajSJAsRGQy8A3gCnxpjXnbJidL+utTHifQTfLv9WybGTuTEmRNcWv9SHu36KHWC6rgkhLLStFYQ9w5oxr0DmnEyLYMVcUksi0tiw/4TfL1yH59nZAPg6SHEhPvTsEYgjWoGEhPuT3SoH1EhfkQF+1ErxBc/b083/2uUUu5UKZKFiHgCHwCDgHhglYjMMsbElvnJUo9h/MM5mprIisMrWLBvAQv2LyAjO4PedXtzV/u7aB/ZvsxP62qh/t5c1iaay9rYvb0zs7LZlZjC5oMn2X00hbijKew5msLqPcdIcTLhL9jPi1B/b0L9vQnx8z73ODTAmwAfT/y9PfHzdtz7eOLn5YG/o9zXyxMvT8HbU/D08MDLQ/DyFLzyPvYUvDyk3Pt9lFKFqxTJAugG7DTGxAGIyBRgKFDmyeLerP2s884m+dsBAIT7hjO8xXCua3odLSJalPXp3MbL04MW0cG0iA6+oNwYQ3JaJoeT00lw3I6cOkPiqTOcTMvgZFoGyWkZ7Eo8fe7nM5nZZRqbp4fgISAICOceewiICGKL8fAQBFvmIcC5Y8BDzj8njjJbkr+CclRh6auwBFdo+nPnuVWV8+P9vfD1KtvWgMqSLOoC+3P9HA9cnPcgERkNjAaoX79+iU5UP7g+Ud5+NGh2JZ2jOtMyoiWeHtWnCUZECA2wNYa8iSQ/WdmG9Iws0jOySHPcp2dkX/A4MyubzGxDZnY2mVnG8diQmZVNVrYhI8tccIwxkG3AYMBAtjEXlBljE1vun7MNgCE725ZlG+xxjucLUtBAj8KGgBT63oW+vuTnLuwAU/g7qCqosC9GJVFZkoWzf/lf/hcYY8YD48GOhirJiR67cWZJXlateXoIgb5eBPpWlo+TUqq4KssQmHigXq6fY4CDbopFKaWqncqSLFYBzUSkkYj4ACOAWW6OSSmlqo1K0W5gjMkUkXuBX7BDZz8zxmx2c1hKKVVtVIpkAWCMmQ3MdnccSilVHVWWZiillFJupMlCKaVUoTRZKKWUKpQmC6WUUoWqskuUi0gisLeEL68JHC3DcMqKxlU8GlfxaFzFU1XjamCMicxbWGWTRWmIyGpn67m7m8ZVPBpX8WhcxVPd4tJmKKWUUoXSZKGUUqpQmiycG+/uAPKhcRWPxlU8GlfxVKu4tM9CKaVUobRmoZRSqlCaLJRSShVKk0UuIjJYRLaJyE4RGefGOOqJyAIR2SIim0XkAUf5syJyQETWO25XuCG2PSKyyXH+1Y6yCBGZJyI7HPfh5RxTi1zXZL2IJIvIg+66XiLymYgcEZE/c5Xle41E5AnHZ26biFxeznG9JiJbRWSjiHwvImGO8oYikpbr2n1cznHl+7tz8/WamiumPSKy3lFeLtergL8Nrv98GWP0ZvttPIFdQGPAB9gAtHZTLLWBTo7HwcB2oDXwLPB/br5Oe4CaecpeBcY5Ho8DXnHz7/Ew0MBd1wvoA3QC/izsGjl+rxsAX6CR4zPoWY5xXQZ4OR6/kiuuhrmPc8P1cvq7c/f1yvP8G8C/y/N6FfC3weWfL61ZnNcN2GmMiTPGnAWmAEPdEYgx5pAxZq3j8SlgC3Yf8opqKPCF4/EXwLXuC4WBwC5jTEln75eaMeZ34Fie4vyu0VBgijHmjDFmN7AT+1ksl7iMMXONMZmOH5djd6EsV/lcr/y49XrlEBEBhgNfu+LcBcSU398Gl3++NFmcVxfYn+vneCrAH2gRaQh0BFY4iu51NBl8Vt7NPQ4GmCsia0RktKMsyhhzCOyHGajlhrhyjODC/8Duvl458rtGFelzdzswJ9fPjURknYgsEpHebojH2e+uolyv3kCCMWZHrrJyvV55/ja4/POlyeI8cVLm1nHFIhIEfAc8aIxJBj4CmgAdgEPYanB5u8QY0wkYAtwjIn3cEINTYrfcvQb41lFUEa5XYSrE505EngIygcmOokNAfWNMR+Bh4CsRCSnHkPL73VWI6wX8jQu/lJTr9XLytyHfQ52Uleh6abI4Lx6ol+vnGOCgm2JBRLyxH4bJxpjpAMaYBGNMljEmG/gEF1W/C2KMOei4PwJ874ghQURqO+KuDRwp77gchgBrjTEJjhjdfr1yye8auf1zJyIjgauAvxtHQ7ej2SLJ8XgNtq27eXnFVMDvriJcLy/gemBqTll5Xi9nfxsoh8+XJovzVgHNRKSR4xvqCGCWOwJxtIdOALYYY97MVV4712HXAX/mfa2L4woUkeCcx9jO0T+x12mk47CRwMzyjCuXC77tuft65ZHfNZoFjBARXxFpBDQDVpZXUCIyGHgcuMYYk5qrPFJEPB2PGzviiivHuPL73bn1ejlcCmw1xsTnFJTX9crvbwPl8flyde99ZboBV2BHF+wCnnJjHL2wVcWNwHrH7QrgS2CTo3wWULuc42qMHVmxAdicc42AGsB8YIfjPsIN1ywASAJCc5W55XphE9YhIAP7ze6Ogq4R8JTjM7cNGFLOce3EtmnnfM4+dhw7zPE73gCsBa4u57jy/d2583o5yv8HjMlzbLlcrwL+Nrj886XLfSillCqUNkMppZQqlCYLpZRShdJkoZRSqlCaLJRSShVKk4VSSqlCabJQ1ZqI1Mi1UujhXCudnhaRD110zgdF5FYn5Q1zr3BaBufxEZHfHZPIlCoV/RCpas3YWbcdwC6LDZw2xrzuqvM5/nDfjl3N1KWMMWdFZD5wE+eX8VCqRLRmoZQTItJPRH50PH5WRL4QkbmOPQyuF5FXxe7r8bNj+QVEpLNjEbk1IvJLnlnIOQZglyTJzPWaDSKyDLgn1/kbishiEVnruPV0lH8pIkNzHTdZRK4RkTYistJRK9ooIs0ch8wA/u6Ka6SqF00WShVNE+BK7JLPk4AFxpi2QBpwpSNhvAfcYIzpDHwGvOjkfS4B1uT6+XPgfmNMjzzHHQEGGbto403Au47yT4F/AohIKNATmA2MAd4xxnQAumBnHINdJqNrCf/NSp2jzVBKFc0cY0yGiGzCbrD0s6N8E3bjmxbARcA8u3wPntilIvKqjd2DIOePfZgxZpHjuS+xiyECeAPvi0gHIAvHonTGmEUi8oGI1MIuZvedMSbTUTN5SkRigOnGsXS2MSZLRM6KSLCx+x8oVSKaLJQqmjMAxphsEckw59fJycb+PxJgs5MaQl5pgJ/jsZD/ctEPAQlAe2wLQHqu577ENi2NwPZ/YIz5SkRWYGs/v4jIKGPMb47jffO8Xqli02YopcrGNiBSRHqAXUZaRNo4OW4L0BTAGHMCOCkivRzP5e5bCAUOGbtE9y3YmkqO/wEPOt5js+N8jYE4Y8y72IX32jnKawCJxpiM0v8TVXWmyUKpMmDsVrw3AK+IyAbsaqA9nRw6B7u3c45/Ah84mpHScpV/CIwUkeXYJqiUXOdKwCadz3MdfxPwp4isB1oCEx3l/bF9GkqViq46q1Q5E5HvgcfMhVtyFuf1Adi+kk7GmJOFHDsdeMIYs60k51Iqh9YslCp/47Ad3cUmIpcCW4H3ipAofIAZmihUWdCahVJKqUJpzUIppVShNFkopZQqlCYLpZRShdJkoZRSqlCaLJRSShXq/wOw01zLCr8jtAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# just for fun, (not required) here is the number of all atoms involved as a function fo time\n", "nNi=1000*np.exp(-t/8.8)\n", "nFe=1000-nNi-nCoRK2\n", "plt.plot(t,nNi,label=r'$^{56}$Ni')\n", "plt.plot(t,nCoRK2,label=r'$^{56}$Co')\n", "plt.plot(t,nFe,label=r'$^{56}$Fe')\n", "plt.legend()\n", "plt.xlabel('Time (days)')\n", "plt.ylabel(r'Number of atoms')" ] }, { "cell_type": "code", "execution_count": 57, "id": "11e5113b", "metadata": {}, "outputs": [], "source": [ "# by knowing the masses of the three isotopes, one can compute the luminosity of the SN\n", "# (also not required and just for fun)\n", "mNi=55.942128*1.66e-24\n", "mCo=55.939838*1.66e-24\n", "mFe=55.934936*1.66e-24\n", "c=3e10 # the speed of light\n", "lNiCo=-(mNi-mCo)*c**2*(nNi[1:]-nNi[:-1])/dt\n", "lCoFe=(mCo-mFe)*c**2*(nFe[1:]-nFe[:-1])/dt" ] }, { "cell_type": "code", "execution_count": 48, "id": "901be3aa", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 48, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAEGCAYAAAB7DNKzAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABA70lEQVR4nO3dd3xT1/n48c/jPTGW2RgsscLeiAAhTDuQnbRJaJNCVjPbZrTfjl86krZpmjZt2rRN0ibNTpsNaRbYjDDD3iNM29gYzLAxGG/7/P6QLIyRbdnW8Hjer5delq6O7n10LevxOeeec8QYg1JKKVWXoEAHoJRSqmXTRKGUUqpemiiUUkrVSxOFUkqpemmiUEopVa+QQAfgTSJyDXBNbGzsdwcMGBDocJRSqlXZtGnTSWNM59rbpS1eHjt27FizcePGQIehlFKtiohsMsaMrb1dm56UUkrVSxOFUkqpemmiUEopVa821ZmtlGqfysvLyc7OpqSkJNChtAoREREkJiYSGhrqUXlNFEqpVi87O5vY2FisVisiEuhwWjRjDKdOnSI7OxubzebRa7TpSSnV6pWUlJCQkKBJwgMiQkJCQqNqX5oolFJtgiYJzzX2XGmiUEopVS9NFEoppeqliUIppVS9NFE4ZZ/N5q3db/Hd1O+y/uj6QIejlGqFRIQf/vCHrsfPPPMMjz/+OAATJ05s9v6PHTvGnDlz6Nu3L4MHD+bKK69k3759zd5vQ9pUohCRa0TkXwUFBY1+7Ru73+DpDU+z9uhalmUt80F0Sqm2Ljw8nI8++oiTJ09e9NyaNWuatW9jDDfccANTp07l4MGD7N69m9/97nfk5uY2a7+eaFOJwhjziTHmnri4uEa/dkriFNf95dnLaYuTJSqlfCskJIR77rmHZ5999qLnYmJiPNrH9u3bSU9Pv2j7smXLCA0N5b777nNtGzlyJJMnTwbgz3/+M0OHDmXo0KH85S9/adobqIMOuHMa120coRJBuSkh62wW6WfS6RPXJ9BhKaUa6YlPdrE754xX9zm4Rwd+dc0Qj8o++OCDDB8+nB//+MdNOlZxcTG33XYbH3/88QUD4nbu3MmYMWPcvmbTpk28+uqrrFu3DmMM48ePZ8qUKYwaNapJMdTWpmoUzREWHMaITnbX4xVZKwIYjVKqterQoQNz587lueeea7DsW2+95aoFVN/uuusujhw5wi233OLxMVetWsUNN9xAdHQ0MTEx3HjjjaxcubI5b+MCWqOo4aq+09l4wpEglmcv5/ahtwc2IKVUo3n6n78vPfzww4wePZo77rij3nK33XYbt9122wXbDh8+zLXXXntR89WQIUP44IMP3O7H103lWqOoYWrv8/0UW45voaC08Z3iSillsVi4+eab+fe//93o1+7du5fnn3+eSZMmXbB9+vTplJaW8tJLL7m2bdiwgeXLl3P55ZezYMECioqKOHfuHPPnz3f1XXiDJooaOkV2onNYPwAqTSVrcpp3lYJSqv364Q9/6Pbqp4YkJye7vZRWRJg/fz5paWn07duXIUOG8Pjjj9OjRw9Gjx7N7bffjt1uZ/z48dx9991e658AbXq6yMTuk/k48wDgaH6abZsd4IiUUq1FYWGh637Xrl0pKipy+1xT9ejRg/fee8/tc48++iiPPvpos4/hjtYoavnmoBTX/VVHVlFRVRHAaJRSKvA0UdQyossQQkxHAApKC9h+YntgA1JKqQDTRFGLiNAnapzr8fLs5QGMRimlAk8ThRvTa1z9tCJbx1Mopdo3TRRu3DJsOqbK0c9/4PQBss9mBzgipZQKHE0UbnSKjiWqaqDrsTY/KaXaM00UdRgcN951X5uflFLtmSaKOszuO911f8OxDRSWNf8aaKWUao3aVKJoznoUtaUMGEhlSXcAyqvKWXVkVbP3qZRq2/y5sFBwcDAjR4503TIyMnxyHGhjiaI561HUFh8dRlzV+SHwSw4vafY+lVJtl78XFoqMjGTr1q2um9Vq9clxoI0lCm8b2+Vy1/2VR1ZSVlkWwGiUUi1ZfQsLNXVRoboWMarPW2+9hd1uZ+TIkdx7771UVlY26vXu6FxP9ZjZdyRLN1gICsvjXPk51h1dx+RE783IqJTygS9+Csd2eHef3YbB7N/XW6SuhYWas6hQXYsYVT83cuRIAGw2G/Pnz2fPnj28++67rF69mtDQUB544AHefvtt5s6d6/l7dUMTRT3G2xKoWDqEsATHAiBLDi/RRKGUapSaiwoBrkWFaieKt956i9///uJkdPToUW655RbWr19/wfbqpqealixZwqZNmxg3zjG7RHFxMV26dGn2e9BEUY9ucREkyBjO4kgUy7KW8YuqXxAcFBzgyJRSdWrgP39fqWthIU8XFWrMIkZ1McYwb948nnrqKY/Ke0r7KBowoedoqIwFIK8kj20ntgU4IqVUS1TXwkKjRo1q8qJCdS1iVJcZM2bwwQcfcPz4cQDy8vLIzMxs/JupRRNFA8bbOlF2ZpDrsV79pJRyp76FhZq6qFBdixjVZfDgwfz2t78lJSWF4cOHk5yczNGjR5v6llzE12utBsLYsWPNxo0bvbKvjJPnmPHCP4nq/SoAiTGJfH7j54iIV/avlGq+PXv2MGjQoIYLKhd350xENhljxtYuqzWKBiQlRGEJGkwwEQBkF2azL983A2iUUqol0kTRABFhvK0LFJ3PvEsPLw1gREop5V+aKDxgt1kozNN+CqVU+6SJwgN2m4WKc5cQLKEA7M3fS+aZ5l9JoJRSrYEmCg8M6BJLXHgMCUHDXdtSM1IDGJFSSvmPJgoPBAUJ46zxFOUPdW1bmLEwgBEppZT/aKLwkN1m4ejRPoQFhQOwL38fhwoOBTgqpVSgnTp1yjXVd7du3ejZs6frcVnZhROJ/uUvf6GoqKjBfU6dOhVvXeLvDZooPGS3JUBVOP1jz19irM1PSqmEhATXVN/33XcfjzzyiOtxWFjYBWU9TRQtjSYKDw3p0YGosGCiy8/PDrkoY1EAI1JKtVRLlixh1KhRDBs2jDvvvJPS0lKee+45cnJymDZtGtOmTQPg/vvvZ+zYsQwZMoRf/epXAY66bjopoIdCg4MY3TuerCNBRHSKoKSyhAOnD3Ag/wD94vsFOjyllNOw14f5bN875jU8fXlJSQm33347S5YsYcCAAcydO5cXXniBhx9+mD//+c8sW7aMTp06AfDkk09isViorKxkxowZbN++neHDhzdwBP/TGkUj2G0W9ueWc2n38xN0pWZq85NS6rzKykpsNhsDBgwAYN68eaxYscJt2ffee4/Ro0czatQodu3axe7du/0Zqsc0UTSC3WbBGEgKP58oFmYs9HgaYaVU21e97kRD0tPTeeaZZ1iyZAnbt2/nqquuoqSkxMfRNY02PTXCyF4dCQ0Wigv6ExkSSXFFMekF6ew/vZ8B8QMCHZ5SCs+ah3yppKSEjIwMDhw4QL9+/XjzzTeZMmUKALGxsZw9e5ZOnTpx5swZoqOjiYuLIzc3ly+++IKpU6cGNPa6tKkahYhcIyL/Kigo8Mn+I0KDGZHYkU0Z55iSOMW1fWG6jqlQSjlERETw6quvctNNNzFs2DCCgoJc62jfc889zJ49m2nTpjFixAhGjRrFkCFDuPPOOz1ecyIQdJrxRvrDwq/514pDPHdXCD9Z9UMAesb05Isbv9Cpx5UKEJ1mvPF0mnEfstssVFQZoiqHEhvmWPnuSOERtp7YGtjAlFLKRzRRNNKYpHiCBDZnFnKF9QrX9k8OfhLAqJRSync0UTRSbEQog3t0YH36Ka7uc7Vr+6KMRZRVltXzSqWUL7XFZnRfaey50kTRBHZrAlsOn2awZTg9Y3oCcKbsDCuzVwY4MqXap4iICE6dOqXJwgPGGE6dOkVERITHr9HLY5vAbovnldXp7M45y5W2K3lpx0sAfHLoE2YkzQhwdEq1P4mJiWRnZ3PixIlAh9IqREREkJiY6HF5TRRNMM5qAWBdeh7XjLrGlSiWZy+noLSAuPC4QIanVLsTGhqKzWYLdBhtljY9NUFCTDj9usSwPj0PW5yNoQmOdSoqqip0okClVJtTb6IQkUQR+ZGIfCwiG0RkhYg8LyJXiUi7TjJ2m4VNGflUVhmu7nu+U/vTQ58GMCqllPK+Or/sReRV4BWgDHga+BbwALAYmAWsEpHL/RFkS2S3WjhbWsGeo2eYZZ1FsAQDsOX4FjIKMgIbnFJKeVF9tYI/GWNSjDHPGWPWGGMOGGN2GmM+MsZ8H5gK5PgnzJbHbnP0U6xPzyMhMoHLE8/nzI8OfBSosJRSyuvqTBTGmJ0AInK1u2YmY0yZMeaAL4NryXp0jCQxPpINGXkAfKP/N1zPfXzgY8qrygMVmlJKeZUn/QxzgP0i8gcR0clUarBbLaxPz8MYw6Sek+gS2QWAvJI8VmS5n39eKaVamwYThTHmNmAUcBB4VUS+EpF7RCTW59G1cHabhVPnyjh44hwhQSFc1+8613Mf7v8wgJEppZT3eHTlkjHmDPAh8A7QHbgB2Cwi3/dhbC1ezX4KgBv73+h6bnXOao6dOxaQuJRSypsaTBTONR7mA0uBUMBujJkNjAB+5OP4WjRbp2g6xYS7+ikSYxO5tPulAFSZKhYcWBDA6JRSyjs8qVHcBDxrjBlujPmjMeY4gDGmCLjTp9G1cCKC3RbvqlHAhZ3a8/fPp8pUBSI0pZTyGk/6KOYaY9z2zBpjlng/pNbFbrVw5HQx2flFAEzvPd01hUfOuRxWHVkVyPCUUqrZPGl6utQ5KrtQRMpEpFJEzvgjuNbAbksAcDU/hQWHcX3f613P/+fr/wQiLKWU8hpPmp7+jmNU9n4gErgb+Jsvg2pNLukWS2xEyAXNT7cMvAXBsSzq6iOryTyTGajwlFKq2Ty96ukAEGyMqTTGvApM821YrUdwkDDOamFdjUTRK7bXBSO13/n6nUCEppRSXuFJoigSkTBgq3PQ3SNAtI/jalXsNguHTpzjxNlS17ZvD/y26/6CAwsoKi8KRGhKKdVsniSK7wDBwPeAc0Av4Bv1vqKdqR5PsTHjfK3i0h6XYu1gBaCwvFDX1FZKtVqeXPWUaYwpNsacMcY8YYx5tD3P8eTO0B5xRIQGXdD8FCRBzBk4x/X4v1//V5dpVEq1SvVNM75DRLbXdfNnkC1dWEgQo3tfOJ4C4Lq+1xEVEgXAwYKDfJXzVSDCU0qpZqmvRnE1cA2w0Hm71Xn7HPjA96G1LnabhT3HznCm5PyssTFhMVzf73rX41d2vRKAyJRSqnnqm2Y80xiTCUwyxvzYGLPDefspcIX/Qmwd7FYLxsCmjPwLts8dMte1qNG6o+vYdWpXIMJTSqkm86QzO1pELqt+ICIT0aueLjKqdzwhQXJBPwVAz5iepFhTXI9f3fmqv0NTSqlm8SRR3AX8Q0QyRCQdeJ52PseTO5FhwQxPjGN9+qmLnrtz6PnTlZaZRtaZLH+GppRSzVJfZ/YEERFjzCZjzAhgODDSGDPSGLPZfyG2HuNsFnYcKaC4rPKC7QMtA5nYYyLgmFX29d2vByI8pZRqkvpqFPOATSLyjojcDkQZYwr8E1brNN5mobzSsCUr/6LnatYqFhxYwMnik/4MTSmlmqy+zuz7jDGjgceBeOA15+p2vxORy0WcPbTKZUySBREuukwWwN7NzpCEIQCUVpbyyk69Akop1Tp4MuDua2PMs8aYWcB0YBWONSrW+Tq41iYuMpRB3Tq4ZpKtSUS4b8R9rsfv7X2PE0Un/BmeUko1iSfTjFuqbzhmj10LPGqMGevz6Fohu83Cpsx8yiouXrBoSuKUC2oV/975b3+Hp5RSjebJVU+bgRPAPhxTjZ8A0kVks4iM8WVwrZHdZqGkvIqdORd354gID4x8wPX4/b3vk3su15/hKaVUo3mSKBYCVxpjOhljEoDZwHvAAzgulfU5EYkWkU0icrU/jtcc46yOCQLd9VMATO45meGdhgNQVlXGyzte9ltsSinVFJ4kirHGmEXVD4wxqcDlxpi1QHh9LxSRV0TkuIjsrLV9lojsFZEDIvJTD2L4CY7k1OJ1jg2nT6doNtSRKGrXKj7Y94EubKSUatE8SRR5IvITEUly3n4M5Duverq4If5CrwGzam5wvu4fOGomg4FvichgERkmIp/WunURkZnAbqDVtNHYbRbWZ+RRWeV+ttiJPSYypquj1a7CVPDXzX/1Z3hKKdUoniSKbwOJwALnrZdzWzBwc30vNMasAGr/a20HDhhjDhljyoB3gOuc80hdXet2HMdqepc6j/ldEXEbs4jcIyIbRWTjiROBvZrIbrNwtqSCvcfOun1eRPjR2B+5HqdlprHl+BZ/haeUUo1Sb6Jw/vf/F2PM940xo5y37xtjThhjypq4LkVPoOYcFtnObW4ZYx4zxjwM/Ad4yRjjthZjjPmXMWasMWZs586dmxCW91QvZOTuMtlqQzsNZbZttuvxMxue0fUqlFItUr2JwhhTCXR2LoXqLeLuUA29yBjzmjHmUy/G4TOJ8VH0iIuos0O72kOjHyI0KBSA7Se3szBjoT/CU0qpRvGk6SkDWC0ivxCRR6tvzThmNo7mq2qJQE4z9tci2W0W1qXn1VtL6BnTk9sG3eZ6/MzGZzhXfs4f4SmllMc8SRQ5wKfOsrE1bk21AegvIjZnTWUO8L9m7K9FstsSOFlYSsaponrL3T38bhIiEgA4XnScF7a+4I/wlFLKY55M4fGEMeYJ4Jnq+87HDRKR/wJfAZeISLaI3GWMqQC+BywC9gDvGWPa3Go+dls8gNtpx2vqENaBH40737H91p632Je/z6exKaVUY3gyhccEEdmN40sdERkhIh4NtDPGfMsY090YE2qMSTTG/Nu5/XNjzABjTF9jzJPNegcXxnqNiPyroCDwk9z27RyDJTrsooWM3LnKdhXjuo0DoNJU8tu1v6XKfZ+9Ukr5nSdNT3/BsfTpKQBjzDbgch/G1GTGmE+MMffExcUFOhREBLvV0mCHdnXZn4//OSESAsCW41v4YJ8uS66Uahk8SRQYY2ovyVbptqC6gN1mITu/mJzTxQ2W7dOxD7cPvd31+JmNz5B1VlfCU0oFnieJIsu5TrYRkTAR+RHOZihVP0/GU9R034j76BPXB4DiimJ+ufqX2gSllAo4TxLFfcCDOAbFZQMjnY9VAwZ170BMeIhH/RQA4cHhPHnZkwQ714TamLuR/+z5jy9DVEqpBnly1dNJY8ytxpiuxpguxpjbjDH1X8qjAAgOEsZa4+ucINCdoZ2GXrBs6l82/4W9eXt9EZ5SSnmkzkQhIj93LlZU1/PTW9q03y3pqqdq46wW9h8v5FRhqcevuX/E/QyIHwA4Fjj60fIfUVRe/3gMpZTylfpqFDuAT0RkiYj8UUR+LCK/FJE3RWQHcA0tbDnUlnTVU7Xxrn6KfI9fExocyh8v/yORIZEAZJzJ4Ddrf6NzQSmlAqLORGGM+dgYMwlHH8UuHLPFngHeAuzGmEeMMbrocwOGJcYRHhLk0WWyNfXp2IefX/pz1+NPD33K/APzvR2eUko1KKShAsaY/TiWQFVNEB4SzKjeHT2+8qmma/tey4ZjG1hwYAEAv137W/rE9WFkl5HeDVIpperh0TgK1Tx2q4VdOQWcLSlv9Gv/3/j/R//4/gCUV5Xz0LKHyClsc3MoKqVaME0UfmC3JVBlYFOm5/0U1SJDIvnb9L8RH+6YOyqvJI/vL/2+zjKrlPIbT+Z6qvPKJ+WZ0UkdCQmSJjU/gWM68menPUtIkKOlcF/+Ph5Z9ghllWXeDFMppdzypEaxTkTeF5ErRcTdokOqAVFhIQzpGdfoDu2axnQdwy8v/aXr8VdHv+JnK39GZZXOpqKU8i1PEsUA4F/Ad4ADIvI7ERng27CapiWOo6g23mZhW1YBJeVN/2K/of8NPDDyAdfj1MxUvWxWKeVznozMNsaYNGPMt4C7gXnAehFZLiITfB5hI7TEcRTV7FYLZZVVbM063az93Df8Pm4ddKvr8Yf7P+TJdU/qnFBKKZ/xpI8iQUQeEpGNwI+A7wOdgB8COhGRh8ZZLYjQqOk83BERfjzux1zd5/yg+Hf3vssTXz2hzVBKKZ9ocBwFjhXq3gSuN8Zk19i+UURe9E1YbU9cVCiXdI1lfRM7tGsKkiB+PenXVFZV8kXGFwB8tP8jSitL+c2k3xAaFNrsYyilVDVP+ih+boz5Tc0kISI3ARhjnvZZZG2Q3WZhU2Y+FZXNbyYKDQrlqclPcV3f61zbPjv0GQ8ufpCzZWebvX+llKrmSaL4qZttP/N2IO2B3WahqKySXTlnvLK/4KBgfj3p19w84GbXtq+OfsXcL+ZytPCoV46hlFL1zR47W0T+BvQUkedq3F4DKvwWYRtitzqGpDTnMtnagiSIn1/68wuuhjpw+gDf/vzbbD2+1WvHUUq1X/XVKHKAjUAJsKnG7X841tBWjdSlQwTWhCiPFzLylIhw/4j7efKyJ12D8k4Wn+SOhXfw9p639fJZpVSz1NmZbYzZBmwTkbeNMa2iBiEi1wDX9OvXL9Ch1Mlus7BoVy5VVYagIO+OX7y277V0i+rGo8sfpaC0gApTwe/X/56tx7fyqwm/IiYsxqvHU0q1D/U1Pb3nvLtFRLbXuO0Qke1+iq9RWvI4imp2WwIFxeXsP17om/13t/Pe1e8xJGGIa9vCjIV885NvsvHYRp8cUynVttXX9PSQ8+fVOBYpqr5VP1ZNcL6fwneryfaI6cEbs9+4oJP7SOER7lx0J3/a+CdKKz1fbU8ppepbuKj6spmTQJYxJhMIB0bg6L9QTdDLEkm3DhFe76eoLSw4jF9M+AV/vPyPdAjrAIDB8Nqu17j5k5vZcGyDT4+vlGo7PLk8dgUQISI9gSXAHcBrvgyqLRMR7DYLGzLy/NLJPMs2i4+u/YhJPSa5th0qOMSdi+7ksVWPkVfi24SllGr9PEkUYowpAm4E/maMuQEY7Nuw2rZxNgu5Z0o5nFfkl+N1je7KCzNf4LHxjxEVEuXa/r+D/+Oa+dfwztfvUF7V+EWVlFLtgydTeIhz8r9bgbsa8TpVh/E2Rz/FuvQ8khKioeQMnNwHpw7C6cNQchpKCqCiFIJCICgYQqMgKgGiEyC6M3TsDfE2iOzo0TFFhDkD5zCt1zSe3vA0aZlpAJwpO8OT657k7T1v8+iYR5naayo6m7xSqiZPvvAfxjESe74xZpeI9AGW+TSqNq5fp2gmR2bQ9auPYdPXcGwH1Jz9NTQaIjpASDhUVUFVBZSfcySP2qISwNIHugyC7iOg2wjoOgTCoi4ui6N28eepf2ZF9gp+t+53HCk8AkDGmQx+sOwHjOk6hkfGPMKIziN88daVUq2QeNpOLiKxOGYd9811nV40duxYs3FjC7wUtKwItrwFa5+H/HRKCSXceikkTYTuIyGhH8QnORKEO5XlUHQKCo87ah55hxy3UwcgdycUO5dalSDoNAB62aH3BOg13pFMatUUSipKeHvP27y842UKyy/8tU7qMYn7RtzHyC4jvX8elFItkohsMsaMvWh7Q4lCRIYBbwAWQIATwFxjzC5fBOoNLS5RGANb/wOLfwXnTkCineWxV/G9LT1J+9k1dIuL8M4xCrLh2HY4uh1ytkDWOkczFkB0F+g9HpIugz5TofMlrsSRV5LHP7f9k/f2vkdFrbGVE7pP4P6R9zOqy6jmx6iUatGakyjWAI8ZY5Y5H08FfmeMmeiDOJulxsjs7+7fvz/Q4TicyYGP7oGMlZBoh5mPQ9JEth8p4Nq/r+a5b43i2hE9fHPsqio4uRcOfwWH1zp+nj7seC62uyNh9JkKtinQoTuZZzJ5cduLfJ7++UULIY3uMpq5Q+YyNXEqwUHBvolXKRVQzUkU24wxIxra1pK0mBpF+kp4/3YoL4ZZv4NRcyHIcaFZRWUVI55I5cbRifzm+qH+iyk/Aw4th0NfOm7FzstjOw+CvtNhQAoZHXvyr12v8ln6ZxcljN6xvblt8G1c1/c6okLd94MopVqn5iSK+cBmHIsXAdwGjDXGXO/tIL2lRSSKrz9zJIl4G9zypqOpp5bv/Hsdx8+UsuiRy/0fHzhqHLk7HAnj4DLIXA2VZRAWC32nktl7HC+VZPFZ9hIqqi5skooLj+Mb/b/BTQNuIjE2MTDxK6W8qjmJIh54ArgMRx/FCuBxY0y+LwL1hoAniq8/g3e/Az1Gwq0fQJTFbbG/L93PM6n72PKLZOKjw/wboztl5xy1jf2LYF8qnHUMwM/tMZz/du7Be0UZnK24cOyHIEzqOYlbLrmFyT0na7OUUq1YkxNFaxTQRHFkE7x6FXQdDHM/hvDYOouuT8/j5n9+xUtzx5I8uKsfg/SAMY7LdquTRvYGigTmJ3TjrbhYsqtKLnpJ9+jufHPAN7mx/410iuwUgKCVUs3RnBrFAOBHgJUa4y6MMdO9HKPXBCxRnDsJL0xyXN569xKI6Vxv8ZLySoY/kcq8CUk8dlULH+x+7hQcWAx7P6fywGJWBlfybsc4VkeEUfsTFCIhTE6czPX9rmdy4mRdw1upVqKuROHJgLv3gReBl4FKbwfWZhgDH3/P0Tn83aUNJgmAiNBgRiZ29OqKdz4TnQAjboERtxBcUcrU9BVM/fpTsvZ9wfshpSyIjSY/2NHsVGEqWJa1jGVZy7BEWLi6z9Vc1+86BsQPCPCbUEo1hSc1ik3GmDF+iscrAlKj2PQ6fPIDmPU0XHqfxy97ZtFeXlh+kO2/SiE6vBXOjFJVBUc2UrbnY1IPfsb7QefYHOF+XMjghMFc3+96rrRdSVx4y10zRKn2qjlNT48Dx4H5gGshA2NMi/032O+JoigP/jbGMY3G7Z9dNAK6Piv2nWDuK+t58y47k/s3XAtp8U7sI3PH23ycsZCPzRmOh1yc/EKDQpncczJX9bmKKb2mEB5cx0h0pZRfNSdRpLvZbIwxfbwVnLf5PVF8+oijRnHvCujWuDERhaUVDH98EQ9O68cPUy6+hLY1qyzIZu3mf7Lg8GKWmrOUuUmgMaExzEyayVV9rmJc13F61ZRSAaRXPfnKqYPw93Ew7i648o9N2sW1f19FRGgw7907wcvBtRwFBYf5YsNzLMhZwS5T7LZMl8guzLLN4qo+VzHIMkhnsVXKzxrdmS0i040xS0XkRnfPG2M+8maArdbKP0FwKEz+UZN3YbdaeGNtJqUVlYSHtM3/qOPiejNn5jPMAdJP7OTzLS/y2bG1ZJnzy7IeLz7OG7vf4I3db2CLs3GV7Squ7HMlvWJ7BS5wpVTdNQoRecIY8ysRedXN08YYc6dvQ2s6v9Uo8tIdfRP2e2D275u8m9Rdx7jnzU28f98ExlndD85ri4wx7Di2kc+2vczC4+vJqzUhYbUhCUOYZZ3FFdYr6B7T3c9RKtV+tIumJ79PCrjoMVj3Ijy8Azo0fWK//HNljPpNGv93xSU8OK2fFwNsPSqqKlibvYrPdr7GkhNbKKbKbbkRnUcwyzqLFGsKXaK6+DlKpdq25nRmdwTmcvGAux94N0Tv8UuNorwY/jTQMfvqza83e3cpzy6nW1wkb9xpb35srVxReRFfZi7hsz3/YU3eLiouGtLnmDpkdNfRzLLOYmbSTB0JrpQXNGfA3efAWmAH1PFvXnu08yPHWg/j7vbK7uw2Cwu25FBRWUVIsCdLmbddUaFRXNnvGq7sdw0FpQUszVzCor3vszZvF5XOpGEwbMrdxKbcTTy17inGdR/HLOssZvSeQXxEfIDfgVJtiyc1is3GmNF+iscr/FKjeGU2FJ2EB9c3atxEXT7eeoSH3tnKJ9+7jGGJOhjNnfySfBZnLmbRvg/YkLeHKjc1jWAJ4tLuE7jCegUzkmbQIaxDACJVqnVqTtPTI0Ah8Ck64M6h4Ag8Oxim/Rym/J9Xdnm0oJgJTy3l51cN4u7JLXaISotxsvgkaZlpLNw3ny35e9ykDMecU5N6TiLFmsLUXlM1aSjVgOY0PZUBfwQeA9ffowHa77fZrvmOn0PdXjncJN3jIultiWJDRp4mCg90iuzEtwZ+i28N/Ba553IdSWP/Arad3usqU2EqWJ69nOXZywmRECb2nEhKkiNp6BQiSnnOk0TxKNDPGHPS18G0Gjs/hO4jIKGvV3drt1lYsicXY4wONmuErtFduW3wbdw2+DZyCnNIzUhl4YGP2VVwwFWmwlSwInsFK7JXECIhXNrjUlKSUpjee7omDaUa4Emv6S6gqMFS7UVBNuRshiE3eH3XdquF/KJyDhwv9Pq+24seMT24fejtvHP9fD6/4XMeHv0wg+MuvOS4wlSw6sgqfrnml0x9dwr3pd3LR/s/4nTJ6cAErVQL50mNohLYKiLLuLCPosVeHutTBxY7fg6Y5fVd222OwXbr0vPo37XuBY+UZ3p16MVdw+7irmF3kXU2i8WZi0k99Ck78/e5ylSYSlbnrGF1zhp+LU9g7zaeFKujpmGJaD+DH5Wqjyed2fPcbTfGNH/wgI/4tDP73dvgyGZ4ZJdXrnaqyRjD+N8tYULfBP46Z5RX963OO1J4xJE0Dn7K9vyv3ZYJliDGdh1HijWFGb1nkBCZ4OcolfK/djEyu5rPEkVlOfyhj6PZ6drnvL9/4MH/bGZTRj5f/Wy69lP4wdHCo6RmppJ26Au25e1yWyYIYWy3saQkOS651cF9qq1q8lVPzmnGL8omLXmacZ/J3gClZ6DfTJ8dYrzNwmfbj5KdX0wvS5TPjqMcusd0Z96QecwbMo9j546RlplG2qEv2HJqh6tMFYb1xzaw/tgGnlz3JGO6jiHFmsLM3jPpHNUG1hBRqgGe9FHUzC4RwE1A+2y8zVzj+Gm9zGeHqO6nWJ+ep4nCz7pFd+M7g7/DdwZ/h9xzuSw+vJjUQ1+w5eS2GteFGzbmbmRj7kaeWvcUo7qMciWNrtFdAxq/Ur7SpKYnEVlljPHdt2Uz+azp6a1vQkEWPLjO+/t2qqoyjPpNGrOGdOPpbw732XGU544XHXf0aaR/weYTW90O7gMY1XkkKdYrmJk0k27R3fwao1Le0Jymp5rTdwThqGG0v0tyqiohaz0M9f5lsTUFBQnjrPGsz2ixA9/bnS5RXfj2oG/z7UHf5mTxSVfS2HR8ywXTiGw5sZUtJ7by9IanGdFpOCnWK0hOStap0VWr50nT059q3K8A0nE0P7Uvx/dAaQH09v0qdHabhcV7jnP8bAldYiN8fjzluU6RnZgzcA5zBs7hZPFJlh5eSmr6F2zI3XRB0th2cjvbTm7njxv/yPCEoSRbryDZmkzPmJ4BjF6ppmkwURhjptXeJiIPA/suLh1YNdaj8P7OD3/l+NlrvPf3XYvd5rgUc0N6PlcN1/9GW6pOkZ24+ZKbufmSm8kryWPJ4SWkpi9kw7GNVNaYaHn7qZ1sP7WTP236E0Msg0i2XkFKUgq9OujKfap1aGofxWFjTG8fxOMVPumj+PhB2PsF/N9Br4+fqK28sorhj6dy89hEnrhuqE+PpbwvvyTfUdPIWMi6o+svSBo1Dew4gBTbLJKTkrHGWf0bpFJuNGdSQLf7a2Y8rc/R7dBtuM+TBEBocBBjkuJZl679FK1RfEQ83xjwDb4x4BucLjnN0qylpGWksvboWipMpavc16f38fWWfTy35Tn6x/Ul2eaoafTt6N05xJRqrqYmirY3Sq8+FWVw4mu49H6/HdJus/Ds4n0UFJUTFxXqt+Mq7+oY0ZEb+9/Ijf1vpKC0gOXZy0lLX8TqnDWU11gjfH/BQfZvfZ7ntz5Pnw5WR59GUjID4gfowEsVcHUmChE5i/uEIECkzyJqiU7uhcoyR43CT+w2C8bAxsw8ZgzS6/PbgrjwOK7tey3X9r2WwrJCvsz+krSMVFZlr6SsRtI4dCaDf27/J//c/k+SYhJJdjZPDbIM0qShAqLORGGMaX+XwNbl6HbHTz8mipG9OhIaLKxP10TRFsWExXB1n6u5us/VnCs/x8rslaRmpLIyezklVWWucpmF2by842Ve3vEyidE9XM1TQxKGaNJQftPUpqf25dh2CI3y+voT9YkIDWZEYkcdT9EORIdGM8s2i1m2WRSVF7E6ZzVpGal8mbWM4krXhM1kn8vh1Z2v8urOV+ke1YVkq6OmMbzzcIKkfa+zrnxLE4Unju+GLoMhKNivh7XbLPxrxSGKyiqICtNfVXsQFRpFclIyyUnJlFSUOJNGGl8eXsK5yhJXuaNFx3lj9xu8sfsNukR0ctU0RnYZqUlDeZ1++3ji5H7oc9FwEp8bZ7Pw/JcH2XL4NJP66Yyl7U1ESAQzes9gRu8ZlFWW8VXOV6RmprIscwlnK865yh0vOcnbe97m7T1v0znCwoykFFKsKYzuMppgP/9zo9omTRQNKT0LZ49CJx8M4mvAmKR4gsSxkJEmivYtLDiMKb2mMKXXFMonPM7ao2tJy0xjaWYaBeXnV0Q8UZLHO3vf4Z2972AJ78gMZ+1kXLdxhATpn7tqGv3kNOSUc93lhP5+P3SHiFAG9+jABh1PoWoIDQ5lcuJkJidO5hcTfsGGYxtIy0xjSUYq+WVnXOXySk/z/r73eX/f+3QMjWWGNYXkpGTs3e2EBukl18pzmigactKZKDoNCMjh7dYE3l6XSVlFFWEh2vasLhQaFMrEHhOZ2GMij41/jM25m0nNTGVxxiJOlZ52lTtdfpYP93/Ih/s/pENoDNN6zyDFmsKl3S8lLDgscG9AtQqaKBpych9IEFhsATm83RbPK6vT2XHkNGOS2ucyIMozIUEh2LvbsXe38zP7z9hyfAtpmWkszljE8ZJTrnJnygv5+ODHfHzwY2JCIpnaewYpSSlM7DmR8ODwAL4D1VJpomjIqf3QMQlCAvMHNM5avZBRviYK5bHgoGDGdhvL2G5j+Yn9J2w/sd2x5Gv6Qo4Vn3CVK6wo5tNDn/LpoU+JCo5gSq9ppFhTmNRzEpEh7WtcraqbrpndkH9eDtGd4bYPvbO/Jpj55+X0io/k1TvsAYtBtQ3GGHae3ElaZhqp6V9wpOiY23KRQeFM7jWFZGsyl/e8nKhQXW2xPfD2pIDtx+ks6DkmoCGMs1r4dFsOlVWG4CAdjauaTkQY1nkYwzoP45Exj7A7bzdpGWmkpS/k8LkjrnLFVaWkZqaSmplKRFAYlyVOJjkphcsTLycmLCaA70AFgiaK+pQWQnEexAV23YDxNgv/XX+YPUfPMLRnXEBjUW2HiDAkYQhDEobw0OiH2Je/z5Ec0heScfawq1xJVRmLDy9h8eElhAWFMLHHZaRYU5jSawodwjoE8B0of9FEUZ+CLMfPjoFdesNuc/RNbMjI00ShfEJEuMRyCZdYLuF7I7/HwdMHXc1TB86ku8qVVVXwZfaXfJn9JSESzIQeE0lOSmZ67+nEhetns63SRFGf085EEeAaRY+OkSTGR7I+PY87JgXm6ivVfogI/eL70S++H/ePvJ9Dpw+RlplGWsZC9p4+4CpXYSpZeWQlK4+s5NdfPc747peSnJTC9N7TiY+ID+A7UN6miaI+Bc7qd8fAL1lpt1pYvu8ExhidNVT5VZ+Ofbi3473cO+JeMs9kOpJG+iJ253/tKlNhqlids4bVOWv4zVe/Zmy3saRYr2B67+l0itRZBVo7TRT1OZ0FQaEQ0y3QkWC3WfhoyxEOnTxH387amagCI6lDEncPu5u7h91N9tlsV01jx6ndrjKVVLHu2HrWHVvPk2t/y5iuY0i2pjCj9wy6RHUJYPSqqTRR1Of0YYhLhKDAj4iu7qdYn56niUK1CImxidwx9A7uGHoHOYU5LM5cTFrGIrae3O4qU4VhQ+5GNuRu5Kl1TzGq8wiSrVcwM2km3aID/w+Y8kybGkchItcA1/Tr1++7+/fvb/4OX54JoZEw75Pm76uZjDGMe3Ixk/t35tlbRgY6HKXqlHsul8WHF5OWkcrm41swdaycPDxhGCk2x5KvPWJ6+DlK5U5d4yjaVKKo5rUBd88OBetlcMOLzd+XFzzw9ia2ZRWw+qfTAx2KUh45UXSCJYeXkJaxiI25m6iqI2kMtQwm2XYFyb2T6dUh8H2C7ZUOuGssY6AwF2JaTpuq3Wrh8x3HyM4vIjFeR8qqlq9zVGfmDJzDnIFzOFV8iqVZS0nLSGX90fVUUuUqtzNvNzvzdvPspmcZFD/AtXqfNc4auOCViyaKuhTnQ2VZi+jIrma3JQCO8RSaKFRrkxCZwE0DbuKmATdxuuQ0y7KWsShjIeuOrqXCnE8ae/L3sSd/H89teY7+cX1Jsc0iJSmFPh37BDD69k0TRV0Kjzt+tqAaxSXdYomNCGF9eh43jEoMdDhKNVnHiI7c0P8Gbuh/AwWlBXyZ9SVpGamsyVlDualwldtfcJD9W//BP7b+g74drCTbHDWN/h3762XifqSJoi6FuY6fMV0DG0cNwUHCOKuF9bqQkWpD4sLjuK7fdVzX7zrOlp1lefZy0jJSWXVkJWVV55PGwTMZHNz2Ii9uexFrTC9X0hhoGahJw8c0UdSlBSYKcFwmu/Tr45wsLKVTjK4doNqW2LBYru5zNVf3uZpz5edYkb2CtIxUVmavoKSqzFUuozCLl3a8xEs7XqJXdA+Snc1TgxMGa9LwAU0UdalOFLEtK1FUr0+xIT2P2cO6BzgapXwnOjSa2bbZzLbNpqi8iFVHVpGWmcbyrGUUV5a6ymWdy+GVna/wys5X6BHV1dERbk1mWKdhBEngx0C1BZoo6lKYCyEREN6yZscc1jOOiNAg1mdoolDtR1RoFCnWFFKsKZRUlLA6ZzWpGaksP7yMc5XFrnI5Rbm8vvt1Xt/9Ol0jO5FsnUWKNYURnUdo0mgGTRR1Oeu8NLaFVWPDQoIY3Tte+ylUuxUREsGM3jOY0XsGpZWlfJXzFWkZqSw7vISzFUWucrnFJ3lrz1u8tectOkdYmGl1DO4b3WU0wUHBAXwHrY8mirqcO97i+ieq2W0W/rpkP2dKyukQERrocJQKmPDgcKb2msrUXlMpryxn7dG1pGWmsTQzjYLyQle5EyV5/Pfr//Lfr/+LJawjM60pJFuTGdt1LCFB+jXYED1DdSnKgw4tc1oBu9WCMbApI59pA1vO5btKBVJocCiTEyczOXEyv5jwCzYc20BaZhpLMlLJLzvjKpdXdpr39r3He/veIz6sA9OTkklJSmFc93GEBuk/Xu5ooqhLcT50HRLoKNwa1TuekCBhfUaeJgql3AgNCmVij4lM7DGRx8Y/xqbcTaRlprE4YxGnSk+7yuWXneHD/R/y4f4P6RAaw/SkmSQnJTOh+wRCgzVpVNNEUZfifIi0BDoKtyLDghmeGKf9FEp5ICQohPHdxzO++3h+Zv8ZW45vcSWN4yWnXOXOlBey4MACFhxYQGxIFFN7zyDFmsKEHhMID27fl6JronCnogzKCiGy5a7SNc5m4ZVV6ZSUVxIRqh1zSnkiOCiYsd3GMrbbWH5i/wnbT2xnUcYi0jIWklt80lXubEURnxz6hE8OfUJ0cCRTek8jJSmFST0nERESEcB3EBiaKNwpznf8jGq5iWK8zcI/lx9iy+HTTOibEOhwlGp1giSIkV1GMrLLSP5v3P+x8+RO5+p9CzlSdMxV7lxlMZ+nf87n6Z8TGRzO5YlTSbYmM7nnZKJC28eca5oo3Cl2Num04BrFmCQLIo6FjDRRKNU8QRLE8M7DGd55OI+OeZTdebtJy0gjNf0Lss7luMoVV5ayKHMRizIXEREUxmWJk0lOSmFKrylEh0YH8B34liYKd6prFC20jwIgLjKUQd06sD7jFNA/0OEo1WaICEMShjAkYQgPjX6Ivfl7Sc1IJS19IRmFWa5yJVVlLD68hMWHlxAWFMKkHpNJtiYztddUYsNiA/gOvE8ThTtFLb9GAY7xFO9uyKK8sorQYB11qpS3iQgDLQMZaBnI90d9nwOnD7iapw6cSXeVK6uqYFn2MpZlLyNEgpnYYyLJ1hSm9ZpGXHhcAN+Bd2iicMfVR9FyaxTgSBSvrclg55ECRvVu2UlNqdZOROgf35/+8f15YOQDHDp9iLTMNFLTF7Kv4ICrXIWpZMWRlaw4spIQCWJ890tJsV7BtF7TiI9onX+nmijcaQV9FHB+gsD16XmaKJTysz4d+3Bvx3u5d8S9ZBRksPjwYlLTF7Inf6+rTIWpYnXOGlbnrOHXPMG4buNItqYwo/cMEiJbT9+irpntzuLHYc3f4RcnWtxcT7VNf+ZLbJ2i+fft4wIdilIKyDqbxeLMxaRlLGTHqd1uywQhjOk6hmRrCjN7z6RzVGc/R+merpndGEV5jmanFp4kwNH89PmOo1RVGYKCWn68SrV1vWJ7ccfQO7hj6B3kFOY4+jQyFrHt5A5XmSoMG3I3siF3I0+te4pRnUeSYruCGb1n0C265Sy/XE0ThTslBRDROjqg7DYL72zIYm/uWQZ1b1lToivV3vWI6cG8IfOYN2Qex84dY8nhJaSmL2LLia0YHK05BsPmE1vYfGILv1//e0Z0Gkayc6bbHjEtY745TRTulJ6F8NZxeZvddr6fQhOFUi1Xt+hu3DroVm4ddCsnik6w5PAS0jIWsTF3E1Wc7wLYdnIH207u4JmNzzDUMphkmyNp9IrtFbDYtY/CnX+nQGgkzP3Ye0H50MSnljCqdzz/uHV0oENRSjXSyeKTLD28lLSMVDYc20AlVW7LDYq/hBTnOuFJHZJ8Eov2UTRG6VmIbhmdS56w2yysPngKY4yuF6xUK9MpshM3X3IzN19yM/kl+SzLWkZqxiLWHV1LhTmfNPbk72VP/l7+uvmvDIjrR7LtClKSUujTsY/PY9RE4U7p2Ra3BGp97LYEFmzNIeNUEbZObXcaAaXauviIeG7sfyM39r+RgtIClmUtIy0jlTU5a6gwla5y+woOsG/rAf6x9R/07WAlxTab5KRk+nXs55N/FjVRuFN6FsJjAh2Fx+w2xxiK9emnNFEo1UbEhcdxfb/rub7f9ZwtO8uXWV+SlpHG6pyVlFVVuModPJPBC9te4IVtL2CN7U2y9QpuueQWukZ7b4VOnfehNmNaVWc2QN/OMViiw1ifnh/oUJRSPhAbFss1fa/huRnPsWLOKv5w+R9I7j2TiKCwC8plnD3MSzteorii2KvH1xpFbRUlYCohrPXUKEQEu9XinCBQKdWWRYdGM9s2m9m22RSVF7HqyCrSMtNYnrWM4spS+sf1xRpn9eoxNVHUVnrW8bMV1SjA0aG9cNcxck4X06NjZKDDUUr5QVRoFCnWFFKsKRRXFLPmyBqfHEcTRW2tOFEAbMjI47qRPQMcjVLK3yJDIpmRNMMn+9Y+itpaaaIY1L0DMeEhuo62UsrrNFHUVlbo+NmK+igAgoOEsdZ4TRRKKa/TRFFbK61RgGPa8f3HC8k7VxboUJRSbUiLTxQiMlVEVorIiyIy1ecHbMWJYnyNfgqllPIWnyYKEXlFRI6LyM5a22eJyF4ROSAiP21gNwYoBCKAbF/F6lKdKFpZ0xPAsMQ4wkOCtPlJKeVVvr7q6TXg78Ab1RtEJBj4B5CM44t/g4j8DwgGnqr1+juBlcaY5SLSFfgzcKtPIy53DlQJi/LpYXwhPCSYUb07aqJQSnmVTxOFMWaFiFhrbbYDB4wxhwBE5B3gOmPMU8DV9ewuHwiv60kRuQe4x/mwUET21lW2AZ2AkzzRsYkv9xlHXB6SH/gwkgs1Ki4/0rgap6XGBS03trYYl9tpaQMxjqInkFXjcTYwvq7CInIjcAXQEUftxC1jzL+AfzU3OBHZ6G6a3UDTuBpH42qclhoXtNzY2lNcgUgU7qY2rHNRDGPMR8BHvgtHKaVUfQJx1VM2UHOppkQgJwBxKKWU8kAgEsUGoL+I2EQkDJgD/C8AcdSl2c1XPqJxNY7G1TgtNS5oubG1m7h8uhSqiPwXmIqjcyUX+JUx5t8iciXwFxxXOr1ijHnSZ0EopZRqlja5ZrZSSinvafEjs5VSSgWWJooaGjli3Jdx9BKRZSKyR0R2ichDzu2Pi8gREdnqvF0ZgNgyRGSH8/gbndssIpImIvudP+P9HNMlNc7JVhE5IyIPB+J8uZuNoL7zIyI/c37e9orIFX6O648i8rWIbBeR+SLS0bndKiLFNc7bi36Oq87fW4DP17s1YsoQka3O7f48X3V9N/j2M2aM0Zuj+S0YOAj0AcKAbcDgAMXSHRjtvB8L7AMGA48DPwrwecoAOtXa9gfgp877PwWeDvDv8RiOgUN+P1/A5cBoYGdD58f5O92GYyCpzfn5C/ZjXClAiPP+0zXistYsF4Dz5fb3FujzVev5PwG/DMD5quu7waefMa1RnOcaMW6MKQPeAa4LRCDGmKPGmM3O+2eBPTgGKrZU1wGvO++/DlwfuFCYARw0xmQG4uDGmBVA7TlU6jo/1wHvGGNKjTHpwAEcn0O/xGWMSTXGVDgfrsVxqbpf1XG+6hLQ81VNRAS4GfivL45dn3q+G3z6GdNEcZ67EeMB/3J2ToEyCljn3PQ9Z1PBK/5u4nEyQKqIbHJOmwLQ1RhzFBwfZKBLAOKqNocL/4ADfb6g7vPTkj5zdwJf1HhsE5EtIrJcRCYHIB53v7eWcr4mA7nGmP01tvn9fNX6bvDpZ0wTxXmNGjHuDyISA3wIPGyMOQO8APQFRgJHcVR//W2SMWY0MBt4UEQuD0AMbjnH5VwLvO/c1BLOV31axGdORB4DKoC3nZuOAr2NMaOAR4H/iEgHP4ZU1++tRZwv4Ftc+M+I38+Xm++GOou62dboc6aJ4rwWNWJcREJxfBDeNo5pTDDG5BpjKo0xVcBL+KjaXR9jTI7z53FgvjOGXBHp7oy7O3Dc33E5zQY2G2NynTEG/Hw51XV+Av6ZE5F5OCbjvNU4G7WdzRSnnPc34WjXHuCvmOr5vbWE8xUC3Ai8W73N3+fL3XcDPv6MaaI4r8WMGHe2gf4b2GOM+XON7d1rFLsB2Fn7tT6OK1pEYqvv4+gM3YnjPM1zFpsHfOzPuGq44D+9QJ+vGuo6P/8D5ohIuIjYgP7Aen8FJSKzgJ8A1xpjimps7yyO5QAQkT7OuA75Ma66fm8BPV9OM4GvjTGutXH8eb7q+m7A158xf/TUt5YbcCWOqwgOAo8FMI7LcFQPtwNbnbcrgTeBHc7t/wO6+zmuPjiuoNgG7Ko+R0ACsATY7/xpCcA5iwJOAXE1tvn9fOFIVEeBchz/zd1V3/kBHnN+3vYCs/0c1wEc7dfVn7EXnWW/4fz9bgM2A9f4Oa46f2+BPF/O7a8B99Uq68/zVdd3g08/YzoyWymlVL206UkppVS9NFEopZSqlyYKpZRS9dJEoZRSql6aKJRSStVLE4Vqt0QkocaMn8dqzFhaKCLP++iYD4vIXDfbrTVnKvXCccJEZIVzgJhSzaIfItVuGcdo2pHgmNoaKDTGPOOr4zm/tO/EMSupTxljykRkCXAL56fmUKpJtEahVC0iMlVEPnXef1xEXheRVOcaBDeKyB/EsSbHQud0CojIGOeEcJtEZFGt0cXVpuOYYqSixmu2ichXwIM1jm8VkZUistl5m+jc/qaIXFej3Nsicq2IDBGR9c7a0HYR6e8ssgC41RfnSLUvmiiUalhf4CocUza/BSwzxgwDioGrnMnib8A3jTFjgFcAd+vATwI21Xj8KvADY8yEWuWOA8nGMfniLcBzzu0vA3cAiEgcMBH4HLgP+KsxZiQwFsdIYnBMfTGuie9ZKRdtelKqYV8YY8pFZAeOhZEWOrfvwLFozSXAUCDNMRUPwTimf6itO471A6q/6DsaY5Y7n3sTx6SGAKHA30VkJFCJc4I5Y8xyEfmHiHTBMTHdh8aYCmeN5DERSQQ+Ms7pr40xlSJSJiKxxrF2gVJNoolCqYaVAhhjqkSk3Jyf96YKx9+QALvc1AxqKwYinPeFuqd7fgTIBUbgqPWX1HjuTRzNSXNw9HdgjPmPiKzDUetZJCJ3G2OWOsuH13q9Uo2mTU9KNd9eoLOITADHNNAiMsRNuT1APwBjzGmgQEQucz5Xsy8hDjhqHNNsfwdHDaXaa8DDzn3sch6vD3DIGPMcjkn0hju3JwAnjDHlzX+Lqj3TRKFUMxnH0rnfBJ4WkW04ZvSc6KboFzjWYq52B/APZ9NRcY3tzwPzRGQtjmanczWOlYsj4bxao/wtwE4R2QoMBN5wbp+Gow9DqWbR2WOV8iMRmQ/82Fy4jGZjXh+Fo29ktDGmoIGyHwE/M8bsbcqxlKqmNQql/OunODq1G01EZgJfA3/zIEmEAQs0SShv0BqFUkqpemmNQimlVL00USillKqXJgqllFL10kShlFKqXpoolFJK1ev/A4eIE3BVe1Q9AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.semilogy(t[:-1],lNiCo,label=r'Ni $\\to$ Co')\n", "plt.plot(t[:-1],lCoFe,label=r'Co $\\to$ Fe')\n", "plt.plot(t[:-1],lCoFe+lNiCo,label='Total',lw=3)\n", "plt.ylim([1e-5,5e-4])\n", "plt.xlabel('Time (days)')\n", "plt.ylabel('Luminosity (erg/day)')\n", "plt.legend()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "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.9.12" } }, "nbformat": 4, "nbformat_minor": 5 }