diff --git a/examples/jupyter/pincell.ipynb b/examples/jupyter/pincell.ipynb index 9e8775805c..0f11ca6559 100644 --- a/examples/jupyter/pincell.ipynb +++ b/examples/jupyter/pincell.ipynb @@ -166,15 +166,15 @@ }, { "cell_type": "code", - "execution_count": 9, + "execution_count": 7, "metadata": {}, "outputs": [], "source": [ - "zirconium = openmc.Material(2, \"zirconium\")\n", + "zirconium = openmc.Material(name=\"zirconium\")\n", "zirconium.add_element('Zr', 1.0)\n", "zirconium.set_density('g/cm3', 6.6)\n", "\n", - "water = openmc.Material(3, \"h2o\")\n", + "water = openmc.Material(name=\"h2o\")\n", "water.add_nuclide('H1', 2.0)\n", "water.add_nuclide('O16', 1.0)\n", "water.set_density('g/cm3', 1.0)" @@ -189,7 +189,7 @@ }, { "cell_type": "code", - "execution_count": 10, + "execution_count": 8, "metadata": {}, "outputs": [], "source": [ @@ -205,7 +205,7 @@ }, { "cell_type": "code", - "execution_count": 11, + "execution_count": 9, "metadata": {}, "outputs": [], "source": [ @@ -221,7 +221,7 @@ }, { "cell_type": "code", - "execution_count": 12, + "execution_count": 10, "metadata": {}, "outputs": [ { @@ -230,7 +230,7 @@ "True" ] }, - "execution_count": 12, + "execution_count": 10, "metadata": {}, "output_type": "execute_result" } @@ -251,7 +251,7 @@ }, { "cell_type": "code", - "execution_count": 13, + "execution_count": 11, "metadata": {}, "outputs": [ { @@ -266,7 +266,7 @@ " \r\n", " \r\n", " \r\n", - " \r\n", + " \r\n", " \r\n", " \r\n", " \r\n", @@ -274,7 +274,7 @@ " \r\n", " \r\n", " \r\n", - " \r\n", + " \r\n", " \r\n", " \r\n", " \r\n", @@ -306,7 +306,7 @@ }, { "cell_type": "code", - "execution_count": 14, + "execution_count": 12, "metadata": {}, "outputs": [ { @@ -321,7 +321,7 @@ " \r\n", " \r\n", " \r\n", - " \r\n", + " \r\n", " \r\n", " \r\n", " \r\n", @@ -329,7 +329,7 @@ " \r\n", " \r\n", " \r\n", - " \r\n", + " \r\n", " \r\n", " \r\n", " \r\n", @@ -368,7 +368,7 @@ }, { "cell_type": "code", - "execution_count": 15, + "execution_count": 13, "metadata": {}, "outputs": [ { @@ -416,7 +416,7 @@ }, { "cell_type": "code", - "execution_count": 16, + "execution_count": 14, "metadata": {}, "outputs": [], "source": [ @@ -437,7 +437,7 @@ }, { "cell_type": "code", - "execution_count": 17, + "execution_count": 15, "metadata": {}, "outputs": [], "source": [ @@ -480,7 +480,7 @@ }, { "cell_type": "code", - "execution_count": 18, + "execution_count": 16, "metadata": {}, "outputs": [], "source": [ @@ -498,7 +498,7 @@ }, { "cell_type": "code", - "execution_count": 19, + "execution_count": 17, "metadata": {}, "outputs": [], "source": [ @@ -515,7 +515,7 @@ }, { "cell_type": "code", - "execution_count": 20, + "execution_count": 18, "metadata": {}, "outputs": [ { @@ -541,7 +541,7 @@ }, { "cell_type": "code", - "execution_count": 21, + "execution_count": 19, "metadata": {}, "outputs": [], "source": [ @@ -558,7 +558,7 @@ }, { "cell_type": "code", - "execution_count": 22, + "execution_count": 20, "metadata": {}, "outputs": [ { @@ -567,7 +567,7 @@ "(array([-1., -1., 0.]), array([1., 1., 1.]))" ] }, - "execution_count": 22, + "execution_count": 20, "metadata": {}, "output_type": "execute_result" } @@ -585,7 +585,7 @@ }, { "cell_type": "code", - "execution_count": 23, + "execution_count": 21, "metadata": {}, "outputs": [], "source": [ @@ -605,7 +605,7 @@ }, { "cell_type": "code", - "execution_count": 24, + "execution_count": 22, "metadata": {}, "outputs": [], "source": [ @@ -628,7 +628,7 @@ }, { "cell_type": "code", - "execution_count": 25, + "execution_count": 23, "metadata": {}, "outputs": [], "source": [ @@ -648,22 +648,22 @@ }, { "cell_type": "code", - "execution_count": 26, + "execution_count": 24, "metadata": {}, "outputs": [ { "data": { "text/plain": [ - "" + "" ] }, - "execution_count": 26, + "execution_count": 24, "metadata": {}, "output_type": "execute_result" }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAARUAAAD8CAYAAABZ0jAcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy86wFpkAAAACXBIWXMAAAsTAAALEwEAmpwYAAASOElEQVR4nO3dfawc1X3G8e9TB6vBRSEE8xJCkzQgCCjFwJWBEiXQFAQWrUPVVNCKRFEiiyqWmqiqSoVE+TOlSqNSJVAnRQ1SC00VnFhg3tWW0ogEm+A3DMFQRzG2sCHUNIVCTX/9Y2fN+nr37uzu2Zkzs89HWt3debl75s6Z556ZffkpIjAzS+UX6m6AmbWLQ8XMknKomFlSDhUzS8qhYmZJOVTMLKkkoSLpNkl7JW0dMF+Sbpa0Q9JmSef0zLtM0jPFvOtStMfM6pNqpPJ3wGULzL8cOLW4rQJuAZC0CPhaMf8M4GpJZyRqk5nVIEmoRMQjwM8WWGQlcHt0PAYcLelEYDmwIyKej4g3gTuLZc2sod5R0fOcBPy05/GuYlq/6ef1+wWSVtEZ5bBkyZJzTz/99Om01Erbse2pqfzeU870YLVuGzdufCkilo6zblWhoj7TYoHph0+MWAOsAZibm4sNGzaka531dcVHzl5w/gmnnDadJ37jrQVn373lR9N5XjtI0k/GXbeqUNkFnNzz+H3AbmDxgOlWsWEBkpN+bXXQ5KOqUFkHrJZ0J53Tm/0RsUfSPuBUSR8EXgCuAn6vojbNtCaFSBnzt8chU58koSLpDuAi4FhJu4A/A44AiIhbgfXACmAH8Brw2WLeAUmrgfuBRcBtEbEtRZvsbW0LkDI8mqmPmvjVB76mMtwsBskoHDALk7QxIubGWdfvqDWzpKq6pmIV8QilnO7fySOW9BwqLeAgGV/v384Bk4ZDpaEcJOk5YNJwqDSIg6Q6DpjxOVQawGFSL19/GY1DJWMOk7w4XMrxS8qZcqDky/tmYR6pZMSdtTk8ahnMoZIBh0lz+YLu4Xz6UzMHSnt4X3Z4pFITd8B28mmRRypmlphHKhXzCGU2zPKIxSOVCjlQZs8s7nOHSkVmsXNZx6zte5/+TNmsdSjrb5ZOhxwqU+IwsX5mIVxSlT1dsHSppD+W9GRx2yrpLUnHFPN2StpSzGvFd0Q6UGyYNveRiUOlTOnSiPiLiFgWEcuAPwX+NSJ6KxpeXMwf6zsxc9LmzmJptbWvpDj9OVi6FKAow7ESGFS+7mrgjgTPm5W2dhCbrjaeDqU4/RlU0vQwko6kU8j9Oz2TA3hA0saitGnjOFBsUm3qQylCpXTpUuA3gX+fd+pzYUScQ+f06QuSPtb3SaRVkjZI2rBv377JWpxQmzqD1astfSlFqAwqadrPVcw79YmI3cXPvcBaOqdTh4mINRExFxFzS5eOVTc6ubZ0AstHG/pUilB5nKJ0qaTFdIJj3fyFJL0L+DjwvZ5pSyQd1b0PXApsTdAmM6vJxKESEQeAbunS7cC3I2KbpGslXduz6JXAAxHx3z3TjgcelbQJ+CFwT0TcN2mbqtCG/yiWp6b3LZc9HUPTd7o1Q52vCLnsaYUcKFaVpvY1v02/pKbuYGu2Jr6PxSOVEhwoVrcm9UGHyhBN2pnWbk3piw6VBTRlJ9rsaEKfdKiYWVIOlQGa8B/BZlPufdOh0kfuO80s5z7qUJkn551l1ivXvupQMbOkHCo9ck1+s0Fy7LMOFTNLyqFSyDHxzcrIre86VMhvp5iNKqc+PPOhktPOMJtELn155kPFzNKa2a8+yCXVzVLK4asSPFIxs6RmMlQ8SrG2q7OPV1VL+SJJ+3vqKd9Qdt3UHCg2K+rq6xNfU+mppXwJnRpAj0taFxHzy57+W0RcMea6ZtYQKUYqB2spR8SbQLeW8rTXHZlHKTZr6ujzVdZSvkDSJkn3SjpzxHWzLXtqZoeqqpbyE8D7I+Is4K+B746wbmdihmVPzexwldRSjohXI+Lnxf31wBGSji2zbio+9bFZVXXfr6SWsqQTJKm4v7x43pfLrJuCA8VmXZXHwMSv/kTEAUndWsqLgNu6tZSL+bcCvwP8gaQDwOvAVdGpt9p33UnbZGb1aX0tZY9SzN5W9u37rqVsZtlodah4lGJ2qCqOiVaHiplVz6FiZkm1NlR86mPW37SPjdaGipnVw6FiZkm1MlR86mO2sGkeI60MFTOrj0PFzJJqXaj41MesnGkdK60LFTOrV2vq/niEYja6adQJ8kjFzJJyqJhZUq0IFZ/6mE0m5THUilAxs3w4VMwsqarKnv6+pM3F7fuSzuqZt1PSlqIcarnviDSzbFVV9vQ/gI9HxCuSLgfWAOf1zL84Il6atC1mVr9Kyp5GxPcj4pXi4WN06vsk4Yu0ZmmkOpaqLHva9Tng3p7HATwgaaOkVYNWctlTs2ZI8Y7a0qVLJV1MJ1Q+2jP5wojYLek44EFJT0fEI4f9wog1dE6bmJuba15dEbMZUUnZUwBJvwp8E1gZES93p0fE7uLnXmAtndMpM2uoqsqe/jJwF3BNRPy4Z/oSSUd17wOXAlvLPrGvp5illeKYqqrs6Q3Ae4CvFyWVDxTVz44H1hbT3gH8Q0TcN2mbzKw+ST6lHBHrgfXzpt3ac//zwOf7rPc8cNb86WbWXH5HrZkl5VAxs6QaGyq+SGs2HZMeW40NFTPLk0PFzJJyqJhZUg4VM0vKoWJmSTlUzCwph4qZJdXIUNmx7anhC5nZ2N71i+88d9x1GxkqZpYvh4qZJeVQMbOkHCpmlpRDxcyScqiYWVIOFTNLqqqyp5J0czF/s6Rzyq5rZs0ycaj0lD29HDgDuFrSGfMWuxw4tbitAm4ZYV0za5BKyp4Wj2+PjseAoyWdWHJdM2uQqsqeDlqmdMnU3rKnb751YOJGm9l0pAiVMmVPBy1TumRqRKyJiLmImFu8KEllETObghRHZ5myp4OWWVxiXTNrkErKnhaPP128CnQ+sD8i9pRc18wapKqyp+uBFcAO4DXgswutO2mbzKw+iuh7CSNrR7/zyPjoKafV3Qyz1np0xzP85+uv9bvmOZTfUWtmSTlUzCwph4qZJeVQMbOkHCpmllQjQ+WUM/2ZQ7Np2v8/r28cd91GhoqZ5cuhYmZJOVTMLCmHipkl5VAxs6QcKmaWlEPFzJJqbKjcveVHdTfBrJUmPbYaGypmlieHipkl5VAxs6QcKmaW1EShIukYSQ9Kerb4+e4+y5ws6Z8lbZe0TdIf9sy7UdILkp4sbitGeX5frDVLK8UxNelI5Trg4Yg4FXi4eDzfAeCPIuLDwPnAF+aVNv1qRCwrbusnbI+Z1WzSUFkJfKu4/y3gk/MXiIg9EfFEcf+/gO0MqEJoZs03aagcX9Tvofh53EILS/oAcDbwg57JqyVtlnRbv9OnnnUPlj3dt2/fhM02s2kZGiqSHpK0tc9tpELqkn4J+A7wxYh4tZh8C/AhYBmwB/jKoPV7y54uXbr04HRfVzFLI9WxNLSYWET8xqB5kl6UdGJE7JF0IrB3wHJH0AmUv4+Iu3p+94s9y3wDuHuUxptZfiY9/VkHfKa4/xnge/MXkCTgb4HtEfGX8+ad2PPwSmDrhO0xs5pNGipfBi6R9CxwSfEYSe+V1H0l50LgGuDX+7x0fJOkLZI2AxcDX5qwPWZWs4lqKUfEy8An+kzfTad2MhHxKNC3fGJEXDPJ85tZflrxjlpfrDWbTMpjqBWhYmb5cKiYWVITXVPJSXf4dsVHzq65JWbNMY1LBx6pmFlSrQsVX7Q1K2dax0rrQsXM6uVQMbOkWhkqPgUyW9g0j5FWhoqZ1cehYmZJtTZUfApk1t+0j43WhoqZ1cOhYmZJtTpUfApkdqgqjolWh4qZVa/1oeLRillHVcdC60PFzKo19bKnxXI7i++ifVLShlHXn5RHKzbrqjwGqih72nVxUdp0bsz1J+JgsVlVdd+fetnTKa9vZpmpquxpAA9I2ihp1RjrJyl76tGKzZo6+vzQr5OU9BBwQp9Z14/wPBdGxG5JxwEPSno6Ih4ZYX0iYg2wBmBubi5GWdfMqlNJ2dOiDhARsVfSWmA58AhQan0za44qyp4ukXRU9z5wKW+XNx26fmo+BbJZUVdfr6Ls6fHAo5I2AT8E7omI+xZaf9ocLNZ2dfbxKsqePg+cNcr6ZtZcran7MyrXCbI2ymEU7rfpm1lSMx8qOSS7WQq59OWZDxXIZ2eYjSunPuxQKeS0U8xGkVvfdaiYWVIOlR65Jb7ZMDn2WYeKmSXlUJknx+Q36yfXvupQ6SPXnWXWlXMfdagMkPNOs9mWe990qJhZUg6VBeT+H8FmTxP6pENliCbsRJsNTemLDpUSmrIzrb2a1Adn9qsPRuWvSrA6NClMujxSGVETd7I1U1P7mkNlDE3d2dYcTe5jUy97Kum0otxp9/aqpC8W826U9ELPvBWTtKdKTd7plrem962plz2NiGeKcqfLgHOB14C1PYt8tTs/ItbPX9/MmqXqsqefAJ6LiJ9M+LxZaPp/FMtPG/pUVWVPu64C7pg3bbWkzZJu63f6lLs2dALLQ1v60tBQkfSQpK19bitHeSJJi4HfAv6pZ/ItwIeAZcAe4CsLrD9xLeVpaUtnsPq0qQ9VUva0cDnwRES82PO7D96X9A3g7gXakXUtZb+PxcbRpjDpmnrZ0x5XM+/Upwiirit5uxxqY7Wxk9h0tLWvVFH2FElHFvPvmrf+TZK2SNoMXAx8acL2ZKGtncXSaXMfUUR2ZxJDzc3NxYYNG+puRik+HbJeTQkTSRsjYm6cdf3ZnynztRaD5oRJCn6bfkVmqVPZoWZt3ztUKjRrnctmc5/79KdiPh2aDbMYJl0eqZhZUh6p1MQjlnaa5RFKl0cqNXMnbA/vyw6PVDLgUUtzOUgO51DJSG8HdcDkzWEymE9/MuVOmy/vm4V5pJIxnxblxWFSjkOlARwu9XKYjMah0iC+5lIdB8n4HCoN5YBJz0GShkOlBRww43OQpOdQaRlffynHYTI9fknZzJLySKWlfEp0OI9OquFQmQH9Dqa2B40DpD4ThYqkTwE3Ah8GlkdE3y+OlXQZ8FfAIuCbEdH9guxjgH8EPgDsBH43Il6ZpE1WzvyDrukh4xDJx6Qjla3AbwN/M2gBSYuAr9H5Nv1dwOOS1kXEU7xdi/nLkq4rHv/JhG2yMTRpNOMAydtEoRIR2wEkLbTYcmBHRDxfLHsnnRrMTxU/LyqW+xbwLzhUsjHs4J1W6Dg0mq2KayonAT/tebwLOK+4f0gtZkkDazFLWgWsKh6+Ianxhcf6OBZ4qe5GTEnpbRvyTyo3bd1np4274tBQkfQQcEKfWddHxEIVCQ/+ij7TRi421Fv2VNKGcWuS5Kyt2wXt3bY2b9e4605US7mkXcDJPY/fB+wu7o9Si9nMGqCKN789Dpwq6YOSFgNX0anBDKPVYjazBpgoVCRdKWkXcAFwj6T7i+kHaylHxAFgNXA/sB34dkRsK35F31rMJayZpN0Za+t2QXu3zds1TyNrKZtZvvzZHzNLyqFiZkk1IlQkfUrSNkn/J2ngy3eSLpP0jKQdxTt0sybpGEkPSnq2+PnuAcvtlLRF0pOTvNQ3bcP+/uq4uZi/WdI5dbRzHCW27SJJ+4t99KSkG+po56gk3SZp76D3fY21zyIi+xudzxadRucdt3MDllkEPAf8CrAY2AScUXfbh2zXTcB1xf3rgD8fsNxO4Ni62ztkW4b+/YEVwL103rt0PvCDutudcNsuAu6uu61jbNvHgHOArQPmj7zPGjFSiYjtEfHMkMUOfhwgIt4Euh8HyNlKOh9PoPj5yfqaMrEyf/+VwO3R8RhwdPH+pNw1sW+VEhGPAD9bYJGR91kjQqWkfh8HOKmmtpR1yMcUgEEfUwjgAUkbi48r5KjM37+J+wjKt/sCSZsk3SvpzGqaNnUj77Nsvk8ll48DpLbQdo3way6MiN3FZ6MelPR08R8mJ2X+/lnuoxLKtPsJ4P0R8XNJK4DvAqdOu2EVGHmfZRMqMd2PA9Rmoe2SVOpjChGxu/i5V9JaOsPx3EKlzN8/y31UwtB2R8SrPffXS/q6pGMjoukfNhx5n7Xp9GehjwPkaujHFCQtkXRU9z5wKZ3vsclNmb//OuDTxSsK5wP7u6d/mRu6bZJOUPHxaknL6RxbL1fe0vRG32d1X30ueYX6SjqJ+QbwInB/Mf29wPp5V6p/TOdK/fV1t7vEdr0HeBh4tvh5zPztovOKw6biti3n7er39weuBa4t7ovOF3Y9B2xhwCt5Od5KbNvqYv9sAh4Dfq3uNpfcrjuAPcD/FsfY5ybdZ36bvpkl1abTHzPLgEPFzJJyqJhZUg4VM0vKoWJmSTlUzCwph4qZJfX//7lTdZ+xkhAAAAAASUVORK5CYII=\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAARUAAAD8CAYAAABZ0jAcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy86wFpkAAAACXBIWXMAAAsTAAALEwEAmpwYAAASQ0lEQVR4nO3de6wc5X3G8e9TB0uNi0II5hICSYgohDSFwhGXUgVoAgKrrUPVVNCKRFEiiwpHTVRVpUKi/JlSpVGhCdRNUUFqoakCiQXmrgZKIxJshC9gri5JjC1sCDKkRKFOf/1jZ836ePec2Z13Z96ZfT7S6uzO5ew7Z955zjuzl58iAjOzVH6p6QaYWbc4VMwsKYeKmSXlUDGzpBwqZpaUQ8XMkkoSKpJukrRL0pYR8yXpOknPS9ok6dSBeRdKeqaYd2WK9phZc1KNVP4ZuHCB+RcBxxe3VcANAJKWAF8r5p8EXCrppERtMrMGJAmViHgY+MkCi6wEbomeR4FDJB0FnA48HxHbIuIt4LZiWTNrqXfU9DxHAz8eeLy9mDZs+hnDfoGkVfRGOSxbtuy0E088cTottdLe2Pyjqfzegz967FR+r5W3YcOGVyJi+STr1hUqGjItFph+4MSINcAagLm5uVi/fn261tlQDx33hYUXOPqs6TzxQmNe4Jxt10/neW0fST+cdN26QmU7cMzA4/cBO4ClI6ZbzRYNkIwMa6uDJh91hcpaYLWk2+id3uyJiJ2SdgPHS/og8BJwCfBHNbVpprUpRMqYvz0OmeYkCRVJtwLnAodJ2g78FXAQQETcCKwDVgDPA28Cny3m7ZW0GrgXWALcFBFPpmiTva1rAVKGRzPNURu/+sDXVBY3i0EyDgfMwiRtiIi5Sdb1O2rNLKm6rqlYTTxCKaf/d/KIJT2HSgc4SCY3+LdzwKThUGkpB0l6Dpg0HCot4iCpjwNmcg6VFnCYNMvXX8bjUMmYwyQvDpdy/JJyphwo+fK+WZhHKhlxZ20Pj1pGc6hkwGHSXr6geyCf/jTMgdId3pc9Hqk0xB2wm3xa5JGKmSXmkUrNPEKZDbM8YvFIpUYOlNkzi/vcoVKTWexc1jNr+96nP1M2ax3Khpul0yGHypQ4TGyYWQiXVGVPFyxdKunPJT1R3LZI+oWkQ4t5L0raXMzrxHdEOlBsMV3uI5VDpUzp0oj4m4g4JSJOAf4SeCgiBqu7nFfMn+g7MXPS5c5iaXW1r6Q4/dlXuhSgKMOxEnhqxPKXArcmeN6sdLWD2HR18XQoxenPqJKmB5D0TnqF3L81MDmA+yRtKEqbto4DxarqUh9KESqlS5cCvwv817xTn7Mj4lR6p09XSPrY0CeRVklaL2n97t27q7U4oS51BmtWV/pSilAZVdJ0mEuYd+oTETuKn7uAO+idTh0gItZExFxEzC1fPlHd6OS60gksH13oUylC5TGK0qWSltILjrXzF5L0LuAc4DsD05ZJOrh/H7gA2JKgTWbWkMqhEhF7gX7p0q3ANyPiSUmXS7p8YNGLgfsi4n8Gph0BPCJpI/AD4K6IuKdqm+rQhf8olqe29y2XPZ1A23e6tUOTrwi57GmNHChWl7b2Nb9Nv6S27mBrtza+j8UjlRIcKNa0NvVBh8oi2rQzrdva0hcdKgtoy0602dGGPulQMbOkHCojtOE/gs2m3PumQ2WI3HeaWc591KEyT847y2xQrn3VoWJmSTlUBuSa/Gaj5NhnHSpmlpRDpZBj4puVkVvfdaiQ304xG1dOfXjmQyWnnWFWRS59eeZDxczSmtmvPsgl1c1SyuGrEjxSMbOkZjJUPEqxrmuyj9dVS/lcSXsG6ilfXXbd1BwoNiua6uuVr6kM1FI+n14NoMckrY2I+WVP/zMifmfCdc2sJVKMVPbVUo6It4B+LeVprzs2j1Js1jTR5+uspXyWpI2S7pb0kTHXzbbsqZntr65ayo8D74+Ik4HrgW+PsW5vYoZlT83sQLXUUo6I1yPip8X9dcBBkg4rs24qPvWxWVV336+llrKkIyWpuH968byvllk3BQeKzbo6j4HKr/5ExF5J/VrKS4Cb+rWUi/k3An8A/ImkvcDPgEuiV2916LpV22Rmzel8LWWPUszeVvbt+66lbGbZ6HSoeJRitr86jolOh4qZ1c+hYmZJdTZUfOpjNty0j43OhoqZNcOhYmZJdTJUfOpjtrBpHiOdDBUza45DxcyS6lyo+NTHrJxpHSudCxUza1Zn6v54hGI2vmnUCfJIxcyScqiYWVKdCBWf+phVk/IY6kSomFk+HCpmllRdZU//WNKm4vY9SScPzHtR0uaiHGq574g0s2zVVfb0v4FzIuI1SRcBa4AzBuafFxGvVG2LmTWvlrKnEfG9iHitePgovfo+SfgirVkaqY6lOsue9n0OuHvgcQD3SdogadWolVz21KwdUryjtnTpUknn0QuV3xqYfHZE7JB0OHC/pKcj4uEDfmHEGnqnTczNzbWvrojZjKil7CmApF8HvgGsjIhX+9MjYkfxcxdwB73TKTNrqbrKnh4L3A5cFhHPDkxfJung/n3gAmBL2Sf29RSztFIcU3WVPb0aeA/w9aKk8t6i+tkRwB3FtHcA/xoR91Rtk5k1J8mnlCNiHbBu3rQbB+5/Hvj8kPW2ASfPn25m7eV31JpZUg4VM0uqtaHii7Rm01H12GptqJhZnhwqZpaUQ8XMknKomFlSDhUzS8qhYmZJOVTMLKlWhsobm3/UdBPMOu1Xly4/bdJ1WxkqZpYvh4qZJeVQMbOkHCpmlpRDxcyScqiYWVIOFTNLqq6yp5J0XTF/k6RTy65rZu1SOVQGyp5eBJwEXCrppHmLXQQcX9xWATeMsa6ZtUgtZU+Lx7dEz6PAIZKOKrmumbVIXWVPRy1TumTqYNnTPb/4WeVGm9l0pAiVMmVPRy1TumRqRKyJiLmImHvXkl8es4lmVpcUdX/KlD0dtczSEuuaWYvUUva0ePzp4lWgM4E9EbGz5Lpm1iJ1lT1dB6wAngfeBD670LpV22Rmzamr7GkAV5Rd18zay++oNbOkHCpmlpRDxcyScqiYWVIOFTNLqpWhcvBHj226CWad9uxbuzdMum4rQ8XM8uVQMbOkHCpmlpRDxcyScqiYWVIOFTNLyqFiZkm1NlTO2XZ9000w66Sqx1ZrQ8XM8uRQMbOkHCpmlpRDxcySqhQqkg6VdL+k54qf7x6yzDGS/kPSVklPSvrTgXnXSHpJ0hPFbcU4z++LtWZppTimqo5UrgQejIjjgQeLx/PtBf4sIj4MnAlcMa+06Vcj4pTi5u+qNWu5qqGyEri5uH8z8Mn5C0TEzoh4vLj/BrCVEVUIzaz9qobKEUX9Hoqfhy+0sKQPAL8BfH9g8mpJmyTdNOz0aWDdfWVPd+/eXbHZZjYti4aKpAckbRlyG6uQuqRfAb4FfDEiXi8m3wB8CDgF2Al8ZdT6g2VPly9fvm+6r6uYpZHqWFq07k9EfGLUPEkvSzoqInZKOgrYNWK5g+gFyr9ExO0Dv/vlgWX+EbhznMabWX6qnv6sBT5T3P8M8J35C0gS8E/A1oj423nzjhp4eDGwpWJ7zKxhVUPly8D5kp4Dzi8eI+m9kvqv5JwNXAb89pCXjq+VtFnSJuA84EsV22NmDatU9jQiXgU+PmT6Dnq1k4mIRwCNWP+yKs9vZvnpxDtqfbHWrJqUx1AnQsXM8uFQMbOkKl1TyUl/+PbQcV9ouCVm7TGNSwceqZhZUp0LFV+0NStnWsdK50LFzJrlUDGzpDoZKj4FMlvYNI+RToaKmTXHoWJmSXU2VHwKZDbctI+NzoaKmTXDoWJmSXU6VHwKZLa/Oo6JToeKmdWv86Hi0YpZT13HQudDxczqNfWyp8VyLxbfRfuEpPXjrl+VRys26+o8Buooe9p3XlHadG7C9StxsNisqrvvT73s6ZTXN7PM1FX2NID7JG2QtGqC9ZOUPfVoxWZNE31+0a+TlPQAcOSQWVeN8TxnR8QOSYcD90t6OiIeHmN9ImINsAZgbm4uxlnXzOpTS9nTog4QEbFL0h3A6cDDQKn1zaw96ih7ukzSwf37wAW8Xd500fVT8ymQzYqm+nodZU+PAB6RtBH4AXBXRNyz0PrT5mCxrmuyj9dR9nQbcPI465tZe3Wm7s+4XCfIuiiHUbjfpm9mSc18qOSQ7GYp5NKXZz5UIJ+dYTapnPqwQ6WQ004xG0dufdehYmZJOVQG5Jb4ZovJsc86VMwsKYfKPDkmv9kwufZVh8oQue4ss76c+6hDZYScd5rNttz7pkPFzJJyqCwg9/8INnva0CcdKotow0602dCWvuhQKaEtO9O6q019cGa/+mBc/qoEa0KbwqTPI5UxtXEnWzu1ta85VCbQ1p1t7dHmPjb1sqeSTijKnfZvr0v6YjHvGkkvDcxbUaU9dWrzTre8tb1vTb3saUQ8U5Q7PQU4DXgTuGNgka/250fEuvnrm1m71F329OPACxHxw4rPm4W2/0ex/HShT9VV9rTvEuDWedNWS9ok6aZhp0+560InsDx0pS8tGiqSHpC0Zcht5ThPJGkp8HvAvw9MvgH4EHAKsBP4ygLrV66lPC1d6QzWnC71IUVMXpZY0jPAuQNlS78bESeMWHYlcEVEXDBi/geAOyPi1xZ73rm5uVi/fv3E7Z4mv4/FxpFrmEjaEBFzk6w79bKnAy5l3qlPEUR9F/N2OdTWyrWTWH662lfqKHuKpHcW82+ft/61kjZL2gScB3ypYnuy0NXOYul0uY9UOv1pSs6nP/P5dMgGtSVMqpz++LM/U+bPDBm0J0xS8Nv0azJLncr2N2v73qFSo1nrXDab+9ynPzXz6dBsmMUw6fNIxcyS8kilIR6xdNMsj1D6PFJpmDthd3hf9nikkgGPWtrLQXIgh0pGBjuoAyZvDpPRfPqTKXfafHnfLMwjlYz5tCgvDpNyHCot4HBplsNkPA6VFvE1l/o4SCbnUGkpB0x6DpI0HCod4ICZnIMkPYdKx/j6SzkOk+nxS8pmlpRHKh3lU6IDeXRSD4fKDBh2MHU9aBwgzakUKpI+BVwDfBg4PSKGfnGspAuBvwOWAN+IiP4XZB8K/BvwAeBF4A8j4rUqbbJy5h90bQ8Zh0g+qo5UtgC/D/zDqAUkLQG+Ru/b9LcDj0laGxFP8XYt5i9LurJ4/BcV22QTaNNoxgGSt0qhEhFbASQttNjpwPMRsa1Y9jZ6NZifKn6eWyx3M/BdHCrZWOzgnVboODTarY5rKkcDPx54vB04o7i/Xy1mSSNrMUtaBawqHv5cUusLjw1xGPBK042YkvLbpr+fbkvS6uo+G1pptIxFQ0XSA8CRQ2ZdFRELVSTc9yuGTBu72FBErAHWFG1aP2lNkpx1dbugu9vW5e2adN1FQyUiPjHpLy9sB44ZePw+YEdx/2VJRw3UYt5V8bnMrGF1vPntMeB4SR+UtBS4hF4NZhivFrOZtUClUJF0saTtwFnAXZLuLabvq6UcEXuB1cC9wFbgmxHxZPErhtZiLmFNlXZnrKvbBd3dNm/XPK2spWxm+fJnf8wsKYeKmSXVilCR9ClJT0r6P0kjX76TdKGkZyQ9X7xDN2uSDpV0v6Tnip/vHrHci5I2S3qiykt907bY31891xXzN0k6tYl2TqLEtp0raU+xj56QdHUT7RyXpJsk7Rr1vq+J9llEZH+j99miE+i943ZuxDJLgBeA44ClwEbgpKbbvsh2XQtcWdy/EvjrEcu9CBzWdHsX2ZZF//7ACuBueu9dOhP4ftPtTrht5wJ3Nt3WCbbtY8CpwJYR88feZ60YqUTE1oh4ZpHF9n0cICLeAvofB8jZSnofT6D4+cnmmlJZmb//SuCW6HkUOKR4f1Lu2ti3SomIh4GfLLDI2PusFaFS0rCPAxzdUFvK2u9jCsCojykEcJ+kDcXHFXJU5u/fxn0E5dt9lqSNku6W9JF6mjZ1Y++zbL5PJZePA6S20HaN8WvOjogdxWej7pf0dPEfJidl/v5Z7qMSyrT7ceD9EfFTSSuAbwPHT7thNRh7n2UTKjHdjwM0ZqHtklTqYwoRsaP4uUvSHfSG47mFSpm/f5b7qIRF2x0Rrw/cXyfp65IOi4i2f9hw7H3WpdOfhT4OkKtFP6YgaZmkg/v3gQvofY9Nbsr8/dcCny5eUTgT2NM//cvcotsm6UgV3wEi6XR6x9artbc0vfH3WdNXn0teob6YXmL+HHgZuLeY/l5g3bwr1c/Su1J/VdPtLrFd7wEeBJ4rfh46f7voveKwsbg9mfN2Dfv7A5cDlxf3Re8Lu14ANjPilbwcbyW2bXWxfzYCjwK/2XSbS27XrcBO4H+LY+xzVfeZ36ZvZkl16fTHzDLgUDGzpBwqZpaUQ8XMknKomFlSDhUzS8qhYmZJ/T9i5lecQhCVeQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] @@ -687,22 +687,22 @@ }, { "cell_type": "code", - "execution_count": 27, + "execution_count": 25, "metadata": {}, "outputs": [ { "data": { "text/plain": [ - "" + "" ] }, - "execution_count": 27, + "execution_count": 25, "metadata": {}, "output_type": "execute_result" }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAARUAAAD8CAYAAABZ0jAcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy86wFpkAAAACXBIWXMAAAsTAAALEwEAmpwYAAAQOUlEQVR4nO3de6wc5X3G8ecpsSPFRaVggx0uaZpYIU4VXLPiUqoENwVhq63jqpXsViSKEllUWGqiqhIVEs2fLVUaiSoBOS1qkBpQq+AEkcM1auumkVPOQb4cFwiGGsXYwYakpilVgPbXP3aOWda7Z2fPvDu3/X6k1ZmdmXf3nZ3Z58zszuzPESEASOVnqu4AgHYhVAAkRagASIpQAZAUoQIgKUIFQFJJQsX23bZP2J4fMt2277B92PYB2xt6pt1g+5ls2i0p+gOgOqn2VP5W0g2LTN8kaW122yHpTkmyfZakL2XT10nabntdoj4BqECSUImIPZJ+tMgsWyTdE117JZ1je42kKyQdjojnI+J1Sfdl8wJoqHeU9DwXSvpBz/2j2bhB468c9AC2d6i7l6MVK1Zcfumll06mp8jtjfm5iTzusl+6fCKPi/zm5uZejohVS2lbVqh4wLhYZPyZIyN2SdolSZ1OJ2ZnZ9P1DgP98P3LFp/hogltPv+5f9HJqw+/MZnnxWm2X1hq27JC5aiki3vuXyTpmKTlQ8ajZCMDpEYG9ZWgqY+yQuUBSTtt36fu4c2piDhu+6SktbbfK+lFSdsk/V5JfZpqTQqRPPqXh5CpTpJQsX2vpGslrbR9VNKfSlomSRFxl6QZSZslHZb0mqRPZdPetL1T0iOSzpJ0d0QcStEnvKVtAZIHezPVSRIqEbF9xPSQdPOQaTPqhg4SmsYgGaX3NSFgJoczagEkVdZnKigJeyj5LLxO7LGkR6i0AEGydBwSpUeoNBRBkh4Bkwah0iAESXkImKUjVBqAMKkWn7+Mh1CpMcKkXgiXfPhKuaYIlPpi3SyOPZUaYWNtDvZahiNUaoAwaS4+0D0Thz8VI1Dag3XZxZ5KRdgA24nDIvZUACTGnkrJ2EOZDtO8x8KeSokIlOkzjeucUCnJNG5c6Jq2dc/hz4RN2waFwabpcIhQmRDCBINMQ7ikKnu6aOlS239se192m7f9v7bPzaYdsX0wm9aKuhsECkZp8zZSOFTylC6NiL+IiPURsV7Sn0j654jorWi4MZveKdqfqrV5Y0Fabd1WUhz+nC5dKklZGY4tkv59yPzbJd2b4Hlrpa0bCCarjYdDKQ5/hpU0PYPtd6lbyP3rPaND0qO257LSpo1DoKCoNm1DKUIld+lSSb8p6V/7Dn2uiYgN6h4+3Wz7IwOfxN5he9b27MmTJ4v1OKE2bQyoVlu2pRShMqyk6SDb1HfoExHHsr8nJO1W93DqDBGxKyI6EdFZtWpJdaOTa8tGgPpowzaVIlSeUFa61PZydYPjgf6ZbP+cpI9K+mbPuBW2z14YlnS9pPkEfQJQkcIf1A4rXWr7pmz6XdmsWyU9GhH/3dP8Akm7bS/05WsR8XDRPpWhDf9RUE8/fP+yRn9w625F0mbpdDoxO1vdKS0ECspQZbDYnlvqKR5c+zMmAgVlaeq2xmn6OTV1BaPZmngeC3sqORAoqFqTtkFCZYQmrUy0W1O2RUJlEU1ZiZgeTdgmCRUASREqQzThPwKmU923TUJlgLqvNKDO2yih0qfOKwvoVddtlVABkBSh0qOuyQ8MU8dtllABkBShkqlj4gN51G3bJVRUv5UCjKtO2/DUh0qdVgZQRF225akPFQBpTe1PH9Ql1YGU6vBTCeypAEhqKkOFvRS0XZXbeFm1lK+1faqnnvJtedumRqBgWlS1rRf+TKWnlvJ16tYAesL2AxHRX/b0XyLiN5bYFkBDpNhTOV1LOSJel7RQS3nSbcfGXgqmTRXbfJm1lK+2vd/2Q7Y/NGbb2pY9BfB2ZdVSflLSeyLiMkl/JekbY7Ttjqxh2VMAZyqllnJEvBoRP8mGZyQts70yT9tUOPTBtCp72y+llrLt1c5qm9q+InveV/K0TYFAwbQr8z1QVi3l35H0B7bflPQ/krZFt97qwLZF+wSgOq2vpcxeCvCWvKfvU0sZQG20OlTYSwHeroz3RKtDBUD5CBUASbU2VDj0AQab9HujtaECoBqECoCkWhkqHPoAi5vke6SVoQKgOoQKgKRaFyoc+gD5TOq90rpQAVCt1tT9YQ8FGN8k6gSxpwIgKUIFQFKtCBUOfYBiUr6HWhEqAOqDUAGQVFllT3/f9oHs9l3bl/VMO2L7YFYONd9vRAKorbLKnv6HpI9GxI9tb5K0S9KVPdM3RsTLRfsCoHqllD2NiO9GxI+zu3vVre+TBB/SAmmkei+VWfZ0waclPdRzPyQ9anvO9o5hjSh7CjRDijNqc5cutb1R3VD51Z7R10TEMdvnS3rM9tMRseeMB4zYpe5hkzqdTvPqigBTopSyp5Jk+8OS/lrSloh4ZWF8RBzL/p6QtFvdwykADVVW2dNLJN0v6caI+H7P+BW2z14YlnS9pPm8T8znKUBaKd5TZZU9vU3SeZK+nJVUfjOrfnaBpN3ZuHdI+lpEPFy0TwCqk+Qq5YiYkTTTN+6unuHPSPrMgHbPS7qsfzyA5uKMWgBJESoAkmpsqPAhLTAZRd9bjQ0VAPVEqABIilABkBShAiApQgVAUoQKgKQIFQBJNTJU3pifq7oLQKt9+J26fKltGxkqAOqLUAGQFKECIClCBUBShAqApAgVAEkRKgCSKqvsqW3fkU0/YHtD3rYAmqVwqPSUPd0kaZ2k7bbX9c22SdLa7LZD0p1jtAXQIKWUPc3u3xNdeyWdY3tNzrYAGiTFr+kPKnt6ZY55LszZVlK37Km6ezm65JJLtPrwC8V6DWCoA/aSr4VJsaeSp+zpsHlyl0yNiF0R0YmIzqpVq8bsIoCypNhTyVP2dNg8y3O0BdAgpZQ9ze5/IvsW6CpJpyLieM62ABqkrLKnM5I2Szos6TVJn1qsbdE+AaiOIwZ+hFFrnU4nZmdnq+4G0Fq257J652PjjFoASREqAJIiVAAkRagASIpQAZAUoQIgKUIFQFKECoCkCBUASREqAJIiVAAkRagASIpQAZAUoQIgKUIFQFKECoCkCBUASREqAJIqFCq2z7X9mO1ns78/P2Cei23/o+2nbB+y/Yc90z5v+0Xb+7Lb5iL9AVC9onsqt0j6dkSslfTt7H6/NyX9UUR8UNJVkm7uK236xYhYn91mCvYHQMWKhsoWSV/Nhr8q6eP9M0TE8Yh4Mhv+L0lPqVuZEEALFQ2VC7L6Pcr+nr/YzLZ/QdIvS/pez+idtg/YvnvQ4VNP2x22Z23Pnjx5smC3AUzKyFCx/bjt+QG3sQqp2/5ZSV+X9NmIeDUbfaek90laL+m4pC8Ma0/ZU6AZRhYTi4hfHzbN9ku210TEcdtrJJ0YMt8ydQPl7yLi/p7Hfqlnnq9IenCczgOon6KHPw9I+mQ2/ElJ3+yfwbYl/Y2kpyLiL/umrem5u1XSfMH+AKhY0VD5M0nX2X5W0nXZfdl+t+2Fb3KukXSjpF8b8NXx7bYP2j4gaaOkzxXsD4CKFaqlHBGvSPrYgPHH1K2drIj4jiQPaX9jkecHUD+cUQsgKUIFQFKECoCkCBUASREqAJIiVAAkRagASIpQAZAUoQIgKUIFQFKECoCkCBUASREqAJIiVAAkRagASIpQAZAUoQIgKUIFQFITL3uazXck+y3afbZnx20PoDnKKHu6YGNW2rSzxPYAGmDiZU8n3B5AzZRV9jQkPWp7zvaOJbSn7CnQECNLdNh+XNLqAZNuHeN5romIY7bPl/SY7acjYs8Y7RURuyTtkqROpxPjtAVQnlLKnmZ1gBQRJ2zvlnSFpD2ScrUH0BxllD1dYfvshWFJ1+ut8qYj2wNoljLKnl4g6Tu290v6N0nfioiHF2sPoLnKKHv6vKTLxmkPoLk4oxZAUoQKgKQIFQBJESoAkiJUACRFqABIilABkBShAiApQgVAUoQKgKQIFQBJESoAkiJUACRFqABIilABkBShAiApQgVAUoQKgKQmXvbU9geycqcLt1dtfzab9nnbL/ZM21ykPwCqN/GypxHxTFbudL2kyyW9Jml3zyxfXJgeETP97QE0S9llTz8m6bmIeKHg8wKoqbLKni7YJunevnE7bR+wffegwycAzTIyVGw/bnt+wG3LOE9ke7mk35L0Dz2j75T0PknrJR2X9IVF2lNLGWiAUsqeZjZJejIiXup57NPDtr8i6cFF+kEtZaABJl72tMd29R36ZEG0YKveKocKoKHKKHsq2+/Kpt/f1/522wdtH5C0UdLnCvYHQMUmXvY0u/+apPMGzHdjkecHUD+cUQsgKUIFQFKECoCkCBUASREqAJIiVAAkRagASIpQAZAUoQIgKUIFQFKECoCkCBUASREqAJIiVAAkRagASIpQAZAUoQIgKUIFQFKECoCkitZS/l3bh2z/n+3OIvPdYPsZ24dt39IzfmQtZgDNUnRPZV7Sb0vaM2wG22dJ+pK6dX/WSdpue102eWQtZgDNUihUIuKpiHhmxGxXSDocEc9HxOuS7lO3BrM0fi1mADVXqERHThdK+kHP/aOSrsyG31aL2fbQWsy2d0jakd39qe02Fh5bKenlqjsxIW1dtrYu1weW2nBkqNh+XNLqAZNujYjFKhKefogB48YuW9pb9tT2bEQM/Qynqdq6XFJ7l63Ny7XUtoVqKed0VNLFPfcvknQsGx6nFjOABijjK+UnJK21/V7byyVtU7cGszReLWYADVD0K+Wtto9KulrSt2w/ko0/XUs5It6UtFPSI5KekvT3EXEoe4iBtZhz2FWk3zXW1uWS2rtsLFcfR4z98QYADMUZtQCSIlQAJNWIUCl6OUBd5b1MwfYR2wdt7yvyVd+kjXr93XVHNv2A7Q1V9HMpcizbtbZPZeton+3bqujnuGzfbfvEsPO+lrTOIqL2N0kfVPdknH+S1Bkyz1mSnpP0i5KWS9ovaV3VfR+xXLdLuiUbvkXSnw+Z74iklVX3d8SyjHz9JW2W9JC65y5dJel7Vfc74bJdK+nBqvu6hGX7iKQNkuaHTB97nTViTyWKXw5QV226TCHP679F0j3RtVfSOdn5SXXXxG0rl4jYI+lHi8wy9jprRKjkNOhygAsr6kteb7tMQdKwyxRC0qO257LLFeooz+vfxHUk5e/31bb3237I9ofK6drEjb3Oyrj2J5e6XA6Q2mLLNcbDXBMRx7Jrox6z/XT2H6ZO8rz+tVxHOeTp95OS3hMRP7G9WdI3JK2ddMdKMPY6q02oxGQvB6jMYstlO9dlChFxLPt7wvZudXfH6xYqeV7/Wq6jHEb2OyJe7Rmesf1l2ysjoukXG469ztp0+LPY5QB1NfIyBdsrbJ+9MCzpenV/x6Zu8rz+D0j6RPaNwlWSTi0c/tXcyGWzvdq2s+Er1H1vvVJ6T9Mbf51V/elzzk+ot6qbmD+V9JKkR7Lx75Y00/dJ9ffV/aT+1qr7nWO5zlP3x6mezf6e279c6n7jsD+7Harzcg16/SXdJOmmbNjq/mDXc5IOasg3eXW85Vi2ndn62S9pr6RfqbrPOZfrXknHJb2Rvcc+XXSdcZo+gKTadPgDoAYIFQBJESoAkiJUACRFqABIilABkBShAiCp/wcE0F5YC2dk8AAAAABJRU5ErkJggg==\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAARUAAAD8CAYAAABZ0jAcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy86wFpkAAAACXBIWXMAAAsTAAALEwEAmpwYAAAQQUlEQVR4nO3dbawc5XnG8esqwVLiolKwMQ7BaRpZIY4aXLMyUFcJbgrCVlvXVSuZViSKEllUsdREVSVXSDQfW6o0ElUCclqrQWpBqYITC8xrlNalkROOLXxsYwiGmmJsYUMiaEoV6ubuh50Dy3p3z+yZZ+dt/z9pdWZn5tl9Zmf2OjO7M3s7IgQAqfxc1R0A0C6ECoCkCBUASREqAJIiVAAkRagASCpJqNjeafu07cNDptv2HbaP2Z61vaZn2o22n8mmbU/RHwDVSbWn8g+SbhwxfYOkldltq6Q7Jcn2eZK+kk1fJekm26sS9QlABZKESkTslfSjEbNsknR3dO2TdKHt5ZLWSjoWEc9HxJuS7s3mBdBQ7yrpeS6T9GLP/RPZuEHjrx70ALa3qruXo8WLF191xRVXTKanyO2pN2cn8rirFn10Io+L/Pbv3/9KRCxdSNuyQsUDxsWI8eeOjNghaYckdTqdmJmZSdc7DLTmPy8fOf0KLZvI8/5ML4+cfmDFiyOnozjbLyy0bVmhckJS7xb6PkknJS0aMh4lmy9A6mRQXwma+igrVHZL2mb7XnUPb16LiFO2z0haafsDkl6StEXSH5bUp6nWpBDJo395CJnqJAkV2/dIuk7SEtsnJP2FpPMlKSLukrRH0kZJxyS9IenT2bSztrdJeljSeZJ2RsSRFH3C29oWIHmwN1MdN/GnD/hMZX7TGCTjIGBGs70/IjoLacsZtQCSKuszFZSEPZR85l4n9ljSI1RagCBZuN7XjoBJg1BpKIIkPQImDUKlQQiS8hAwC0eoNABhUi0+fxkPoVJjhEm9EC758JVyTREo9cW6GY09lRphY20O9lqGI1RqgDBpLj7QPReHPxUjUNqDddnFnkpF2ADbicMi9lQAJMaeSsnYQ5kO07zHwp5KiQiU6TON65xQKck0blzomrZ1z+HPhE3bBoXBpulwiFCZEMIEg0xDuKQqezqydKntP7P9ZHY7bPv/bF+UTTtu+1A2rRW/EUmgYD5t3kYKh0qe0qUR8dcRsToiVkv6c0n/GhG9FQ3XZ9MX9JuYddLmjQVptXVbSXH481bpUknKynBskvTUkPlvknRPguetlbZuIJisNh4OpTj8GVbS9By236NuIfdv9owOSY/Y3p+VNm0cAgVFtWkbShEquUuXSvptSf/ed+izLiLWqHv49DnbHxv4JPZW2zO2Z86cOVOsxwm1aWNAtdqyLaUIlWElTQfZor5Dn4g4mf09LWmXuodT54iIHRHRiYjO0qULqhudXFs2AtRHG7apFKHyhLLSpbYXqRscu/tnsv0Lkj4u6ds94xbbvmBuWNINkg4n6BOAihQOlYg4K2mudOlRSd+IiCO2b7F9S8+smyU9EhH/3TNumaTHbR+U9ANJD0TEQ0X7VIY2/EdBPTV926Ls6QI0faWjGar8RoiypyUiUFCWpm5rnKafU1NXMJqtieexsKeSA4GCqjVpGyRU5tGklYl2a8q2SKiM0JSViOnRhG2SUAGQFKEyRBP+I2A61X3bJFQGqPtKA+q8jRIqfeq8soBedd1WCRUASREqPeqa/MAwddxmCRUASREqmTomPpBH3bZdQkX1WynAuOq0DU99qNRpZQBF1GVbnvpQAZDW1P70QV1SHUipDj+VwJ4KgKSmMlTYS0HbVbmNl1VL+Trbr/XUU74tb9vUCBRMi6q29cKfqfTUUr5e3RpAT9jeHRH9ZU//LSJ+a4FtATREij2Vt2opR8SbkuZqKU+67djYS8G0qWKbL7OW8rW2D9p+0PZHxmxb27KnAN6prFrKByS9PyKulPS3kr41RtvuyBqWPQVwrlJqKUfE6xHxk2x4j6TzbS/J0zYVDn0wrcre9kuppWz7UtvOhtdmz/tqnrYpECiYdmW+Bwp/+xMRZ23P1VI+T9LOuVrK2fS7JP2+pD+2fVbS/0jaEt16qwPbFu0TgOq0vpYyeynA2/Kevk8tZQC10epQYS8FeKcy3hOtDhUA5SNUACTV2lDh0AcYbNLvjdaGCoBqECoAkmplqHDoA4w2yfdIK0MFQHUIFQBJtS5UOPQB8pnUe6V1oQKgWq2p+8MeCjC+SdQJYk8FQFKECoCkWhEqHPoAxaR8D7UiVADUB6ECIKmyyp7+ke3Z7PY921f2TDtu+1BWDjXfb0QCqK2yyp7+h6SPR8SPbW+QtEPS1T3T10fEK0X7AqB6pZQ9jYjvRcSPs7v71K3vkwQf0gJppHovlVn2dM5nJD3Ycz8kPWJ7v+2twxpR9hRohhRn1OYuXWp7vbqh8us9o9dFxEnbl0h61PbTEbH3nAeM2KHuYZM6nU7z6ooAU6KUsqeSZPujkv5O0qaIeHVufESczP6elrRL3cMpAA1VVtnTFZLuk3RzRPywZ/xi2xfMDUu6QdLhvE/M5ylAWineU2WVPb1N0sWSvpqVVD6bVT9bJmlXNu5dkv4pIh4q2icA1UlylXJE7JG0p2/cXT3Dn5X02QHtnpd0Zf94AM3FGbUAkiJUACTV2FDhQ1pgMoq+txobKgDqiVABkBShAiApQgVAUoQKgKQIFQBJESoAkmpkqDz15mzVXQBa7d2/cv5VC23byFABUF+ECoCkCBUASREqAJIiVAAkRagASIpQAZBUWWVPbfuObPqs7TV52wJolsKh0lP2dIOkVZJusr2qb7YNklZmt62S7hyjLYAGKaXsaXb/7ujaJ+lC28tztgXQICl+TX9Q2dOrc8xzWc62krplT9Xdy9GKFSt0YMULxXoNYCgf8v6Ftk2xp5Kn7OmweXKXTI2IHRHRiYjO0qVLx+wigLKk2FPJU/Z02DyLcrQF0CCllD3N7n8y+xboGkmvRcSpnG0BNEhZZU/3SNoo6ZikNyR9elTbon0CUB1HDPwIo9Y6nU7MzMxU3Q2gtWzvz+qdj40zagEkRagASIpQAZAUoQIgKUIFQFKECoCkCBUASREqAJIiVAAkRagASIpQAZAUoQIgKUIFQFKECoCkCBUASREqAJIiVAAkRagASKpQqNi+yPajtp/N/v7igHkut/1d20dtH7H9Jz3Tvmj7JdtPZreNRfoDoHpF91S2S/pORKyU9J3sfr+zkv40Ij4s6RpJn+srbfrliFid3fYU7A+AihUNlU2Svp4Nf13S7/bPEBGnIuJANvxfko6qW5kQQAsVDZVlWf0eZX8vGTWz7V+S9KuSvt8zepvtWds7Bx0+9bTdanvG9syZM2cKdhvApMwbKrYfs314wG2sQuq2f17SNyV9PiJez0bfKemDklZLOiXpS8PaU/YUaIZ5i4lFxG8Om2b7ZdvLI+KU7eWSTg+Z73x1A+UfI+K+nsd+uWeer0m6f5zOA6ifooc/uyV9Khv+lKRv989g25L+XtLRiPibvmnLe+5ulnS4YH8AVKxoqPylpOttPyvp+uy+bL/X9tw3Oesk3SzpNwZ8dXy77UO2ZyWtl/SFgv0BULFCtZQj4lVJnxgw/qS6tZMVEY9L8pD2Nxd5fgD1wxm1AJIiVAAkRagASIpQAZAUoQIgKUIFQFKECoCkCBUASREqAJIiVAAkRagASIpQAZAUoQIgKUIFQFKECoCkCBUASREqAJIiVAAkNfGyp9l8x7Pfon3S9sy47QE0RxllT+esz0qbdhbYHkADTLzs6YTbA6iZssqehqRHbO+3vXUB7Sl7CjTEvCU6bD8m6dIBk24d43nWRcRJ25dIetT20xGxd4z2iogdknZIUqfTiXHaAihPKWVPszpAiojTtndJWitpr6Rc7QE0RxllTxfbvmBuWNINeru86bztATRLGWVPl0l63PZBST+Q9EBEPDSqPYDmKqPs6fOSrhynPYDm4oxaAEkRKgCSIlQAJEWoAEiKUAGQFKECIClCBUBShAqApAgVAEkRKgCSIlQAJEWoAEiKUAGQFKECIClCBUBShAqApAgVAEkRKgCSmnjZU9sfysqdzt1et/35bNoXbb/UM21jkf4AqN7Ey55GxDNZudPVkq6S9IakXT2zfHluekTs6W8PoFnKLnv6CUnPRcQLBZ8XQE2VVfZ0zhZJ9/SN22Z71vbOQYdPAJpl3lCx/ZjtwwNum8Z5ItuLJP2OpH/uGX2npA9KWi3plKQvjWhPLWWgAUope5rZIOlARLzc89hvDdv+mqT7R/SDWspAA0y87GmPm9R36JMF0ZzNerscKoCGKqPsqWy/J5t+X1/7220fsj0rab2kLxTsD4CKTbzsaXb/DUkXD5jv5iLPD6B+OKMWQFKECoCkCBUASREqAJIiVAAkRagASIpQAZAUoQIgKUIFQFKECoCkCBUASREqAJIiVAAkRagASIpQAZAUoQIgKUIFQFKECoCkCBUASRWtpfwHto/Y/pntzoj5brT9jO1jtrf3jJ+3FjOAZim6p3JY0u9J2jtsBtvnSfqKunV/Vkm6yfaqbPK8tZgBNEuhUImIoxHxzDyzrZV0LCKej4g3Jd2rbg1mafxazABqrlCJjpwuk/Riz/0Tkq7Oht9Ri9n20FrMtrdK2prd/antNhYeWyLplao7MSFtXba2LteHFtpw3lCx/ZikSwdMujUiRlUkfOshBowbu2xpb9lT2zMRMfQznKZq63JJ7V22Ni/XQtsWqqWc0wlJl/fcf5+kk9nwOLWYATRAGV8pPyFppe0P2F4kaYu6NZil8WoxA2iAol8pb7Z9QtK1kh6w/XA2/q1ayhFxVtI2SQ9LOirpGxFxJHuIgbWYc9hRpN811tblktq7bCxXH0eM/fEGAAzFGbUAkiJUACTViFApejlAXeW9TMH2cduHbD9Z5Ku+SZvv9XfXHdn0WdtrqujnQuRYtutsv5atoydt31ZFP8dle6ft08PO+1rQOouI2t8kfVjdk3H+RVJnyDznSXpO0i9LWiTpoKRVVfd9nuW6XdL2bHi7pL8aMt9xSUuq7u88yzLv6y9po6QH1T136RpJ36+63wmX7TpJ91fd1wUs28ckrZF0eMj0sddZI/ZUovjlAHXVpssU8rz+myTdHV37JF2YnZ9Ud03ctnKJiL2SfjRilrHXWSNCJadBlwNcVlFf8nrHZQqShl2mEJIesb0/u1yhjvK8/k1cR1L+fl9r+6DtB21/pJyuTdzY66yMa39yqcvlAKmNWq4xHmZdRJzMro161PbT2X+YOsnz+tdyHeWQp98HJL0/In5ie6Okb0laOemOlWDsdVabUInJXg5QmVHLZTvXZQoRcTL7e9r2LnV3x+sWKnle/1quoxzm7XdEvN4zvMf2V20viYimX2w49jpr0+HPqMsB6mreyxRsL7Z9wdywpBvU/R2busnz+u+W9MnsG4VrJL02d/hXc/Mum+1LbTsbXqvue+vV0nua3vjrrOpPn3N+Qr1Z3cT8qaSXJT2cjX+vpD19n1T/UN1P6m+tut85lutidX+c6tns70X9y6XuNw4Hs9uROi/XoNdf0i2SbsmGre4Pdj0n6ZCGfJNXx1uOZduWrZ+DkvZJ+rWq+5xzue6RdErS/2bvsc8UXWecpg8gqTYd/gCoAUIFQFKECoCkCBUASREqAJIiVAAkRagASOr/AWK2ZMJSs3TpAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] @@ -726,16 +726,16 @@ }, { "cell_type": "code", - "execution_count": 28, + "execution_count": 26, "metadata": {}, "outputs": [ { "data": { "text/plain": [ - "" + "" ] }, - "execution_count": 28, + "execution_count": 26, "metadata": {}, "output_type": "execute_result" }, @@ -774,7 +774,7 @@ }, { "cell_type": "code", - "execution_count": 29, + "execution_count": 27, "metadata": {}, "outputs": [], "source": [ @@ -792,7 +792,7 @@ }, { "cell_type": "code", - "execution_count": 30, + "execution_count": 28, "metadata": {}, "outputs": [], "source": [ @@ -810,18 +810,18 @@ }, { "cell_type": "code", - "execution_count": 33, + "execution_count": 29, "metadata": {}, "outputs": [], "source": [ - "fuel = openmc.Cell(1, 'fuel')\n", + "fuel = openmc.Cell(name='fuel')\n", "fuel.fill = uo2\n", "fuel.region = fuel_region\n", "\n", - "gap = openmc.Cell(2, 'air gap')\n", + "gap = openmc.Cell(name='air gap')\n", "gap.region = gap_region\n", "\n", - "clad = openmc.Cell(3, 'clad')\n", + "clad = openmc.Cell(name='clad')\n", "clad.fill = zirconium\n", "clad.region = clad_region" ] @@ -835,7 +835,7 @@ }, { "cell_type": "code", - "execution_count": 34, + "execution_count": 30, "metadata": {}, "outputs": [], "source": [ @@ -855,13 +855,13 @@ }, { "cell_type": "code", - "execution_count": 35, + "execution_count": 31, "metadata": {}, "outputs": [], "source": [ "water_region = +left & -right & +bottom & -top & +clad_outer_radius\n", "\n", - "moderator = openmc.Cell(4, 'moderator')\n", + "moderator = openmc.Cell(name='moderator')\n", "moderator.fill = water\n", "moderator.region = water_region" ] @@ -875,7 +875,7 @@ }, { "cell_type": "code", - "execution_count": 36, + "execution_count": 32, "metadata": {}, "outputs": [ { @@ -884,7 +884,7 @@ "openmc.region.Intersection" ] }, - "execution_count": 36, + "execution_count": 32, "metadata": {}, "output_type": "execute_result" } @@ -904,7 +904,7 @@ }, { "cell_type": "code", - "execution_count": 37, + "execution_count": 33, "metadata": {}, "outputs": [], "source": [ @@ -920,7 +920,7 @@ }, { "cell_type": "code", - "execution_count": 38, + "execution_count": 34, "metadata": {}, "outputs": [ { @@ -929,10 +929,10 @@ "text": [ "\r\n", "\r\n", - " \r\n", - " \r\n", - " \r\n", - " \r\n", + " \r\n", + " \r\n", + " \r\n", + " \r\n", " \r\n", " \r\n", " \r\n", @@ -967,7 +967,7 @@ }, { "cell_type": "code", - "execution_count": 39, + "execution_count": 35, "metadata": {}, "outputs": [], "source": [ @@ -985,25 +985,20 @@ }, { "cell_type": "code", - "execution_count": 40, + "execution_count": 36, "metadata": {}, "outputs": [], "source": [ - "# OpenMC simulation parameters\n", - "batches = 100\n", - "inactive = 10\n", - "particles = 1000\n", - "\n", - "# Instantiate a Settings object\n", "settings = openmc.Settings()\n", - "settings.batches = batches\n", - "settings.inactive = inactive\n", - "settings.particles = particles" + "settings.source = source\n", + "settings.batches = 100\n", + "settings.inactive = 10\n", + "settings.particles = 1000" ] }, { "cell_type": "code", - "execution_count": 41, + "execution_count": 37, "metadata": {}, "outputs": [ { @@ -1016,6 +1011,11 @@ " 1000\r\n", " 100\r\n", " 10\r\n", + " \r\n", + " \r\n", + " 0 0 0\r\n", + " \r\n", + " \r\n", "\r\n" ] } @@ -1040,7 +1040,7 @@ }, { "cell_type": "code", - "execution_count": 42, + "execution_count": 38, "metadata": {}, "outputs": [], "source": [ @@ -1059,7 +1059,7 @@ }, { "cell_type": "code", - "execution_count": 43, + "execution_count": 39, "metadata": {}, "outputs": [], "source": [ @@ -1076,7 +1076,7 @@ }, { "cell_type": "code", - "execution_count": 44, + "execution_count": 40, "metadata": {}, "outputs": [ { @@ -1086,7 +1086,7 @@ "\r\n", "\r\n", " \r\n", - " 1\r\n", + " 3\r\n", " \r\n", " \r\n", " 1\r\n", @@ -1114,7 +1114,7 @@ }, { "cell_type": "code", - "execution_count": 45, + "execution_count": 41, "metadata": { "scrolled": true }, @@ -1152,7 +1152,7 @@ " License | https://docs.openmc.org/en/latest/license.html\n", " Version | 0.12.0\n", " Git SHA1 | 3d90a9f857ec72eae897e054d4225180f1fa4d93\n", - " Date/Time | 2020-08-15 07:03:19\n", + " Date/Time | 2020-08-25 14:58:51\n", " OpenMP Threads | 4\n", "\n", " Reading settings XML file...\n", @@ -1286,20 +1286,20 @@ "\n", " =======================> TIMING STATISTICS <=======================\n", "\n", - " Total time for initialization = 6.9749e-01 seconds\n", - " Reading cross sections = 6.8627e-01 seconds\n", - " Total time in simulation = 1.9684e+00 seconds\n", - " Time in transport only = 1.9468e+00 seconds\n", - " Time in inactive batches = 1.5675e-01 seconds\n", - " Time in active batches = 1.8117e+00 seconds\n", - " Time synchronizing fission bank = 4.5360e-03 seconds\n", - " Sampling source sites = 3.6973e-03 seconds\n", - " SEND/RECV source sites = 6.8224e-04 seconds\n", - " Time accumulating tallies = 1.0140e-04 seconds\n", - " Total time for finalization = 5.5400e-05 seconds\n", - " Total time elapsed = 2.6701e+00 seconds\n", - " Calculation Rate (inactive) = 63796.4 particles/second\n", - " Calculation Rate (active) = 49677.1 particles/second\n", + " Total time for initialization = 6.9022e-01 seconds\n", + " Reading cross sections = 6.7913e-01 seconds\n", + " Total time in simulation = 1.7892e+00 seconds\n", + " Time in transport only = 1.7650e+00 seconds\n", + " Time in inactive batches = 1.5005e-01 seconds\n", + " Time in active batches = 1.6391e+00 seconds\n", + " Time synchronizing fission bank = 4.2308e-03 seconds\n", + " Sampling source sites = 3.4593e-03 seconds\n", + " SEND/RECV source sites = 6.2601e-04 seconds\n", + " Time accumulating tallies = 9.5555e-05 seconds\n", + " Total time for finalization = 7.4948e-05 seconds\n", + " Total time elapsed = 2.4836e+00 seconds\n", + " Calculation Rate (inactive) = 66645.8 particles/second\n", + " Calculation Rate (active) = 54907.5 particles/second\n", "\n", " ============================> RESULTS <============================\n", "\n", @@ -1325,7 +1325,7 @@ }, { "cell_type": "code", - "execution_count": 46, + "execution_count": 42, "metadata": {}, "outputs": [ { @@ -1334,7 +1334,7 @@ "text": [ " ============================> TALLY 1 <============================\r\n", "\r\n", - " Cell 1\r\n", + " Cell 3\r\n", " U235\r\n", " Total Reaction Rate 0.726151 +/- 0.00251702\r\n", " Fission Rate 0.543836 +/- 0.00205084\r\n", @@ -1358,7 +1358,7 @@ }, { "cell_type": "code", - "execution_count": 47, + "execution_count": 43, "metadata": {}, "outputs": [], "source": [ @@ -1379,7 +1379,7 @@ }, { "cell_type": "code", - "execution_count": 48, + "execution_count": 44, "metadata": {}, "outputs": [ { @@ -1393,7 +1393,7 @@ " 1.26 1.26\r\n", " 200 200\r\n", " \r\n", - " \r\n", + " \r\n", " \r\n", "\r\n" ] @@ -1414,7 +1414,7 @@ }, { "cell_type": "code", - "execution_count": 49, + "execution_count": 45, "metadata": {}, "outputs": [ { @@ -1450,7 +1450,7 @@ " License | https://docs.openmc.org/en/latest/license.html\n", " Version | 0.12.0\n", " Git SHA1 | 3d90a9f857ec72eae897e054d4225180f1fa4d93\n", - " Date/Time | 2020-08-15 07:03:32\n", + " Date/Time | 2020-08-25 14:58:54\n", " OpenMP Threads | 4\n", "\n", " Reading settings XML file...\n", @@ -1490,7 +1490,7 @@ }, { "cell_type": "code", - "execution_count": 50, + "execution_count": 46, "metadata": {}, "outputs": [], "source": [ @@ -1506,19 +1506,19 @@ }, { "cell_type": "code", - "execution_count": 51, + "execution_count": 47, "metadata": { "scrolled": false }, "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAMgAAADIAgMAAADQNkYNAAAABGdBTUEAALGPC/xhBQAAACBjSFJNAAB6JgAAgIQAAPoAAACA6AAAdTAAAOpgAAA6mAAAF3CculE8AAAADFBMVEUAAP9yEhL//////wAZPRNOAAAAAWJLR0QCZgt8ZAAAAAd0SU1FB+QIDwgDIDdCihAAAAI7SURBVGje7ZmxkcIwEEUhcAnuhxIIEIGvAwioxiU4QASU4GougA4gsM6SfNxh766sfx7N3Ix+zJv9X7KxtLtaZWVlZWWl1Fr12sQQpXLaRRMRTKFe2kbk+Na8POpNkbbmWlMqtkwxRsJllIotMxQ56F7neWX8LmrdGtNpPWc//Z4cLsbpfp6xN85XNRA904Sd2V/stXlJ16EFcL4O7Q/SnUPOrK/9xfzSvQ44K0dFhjLCmq0nRYYyvDPrq2rfka4RnVlfoyJ9GdGZ9XUdI89aWGYbpTITNUKYgvLlnXFhSsqXd7bjEcKXc8YhvYEjhdzY/Db9lUKebP4+/b6lkK7m8nNRhDAlE8WFoREuig/DpKej+DAbGmGiuDAUUsgItWR8ejZ/yaZ3+SmET+/zL4JIC8Ys2VpI7/JPkSKEbCnkyiNPCilDyI5A+AVzSzZFVAhRBCKssVvlBRB5W8iNwZCThDymiLyT5F5iSCsh3RSRN5/cfgwxoiikkpFmgqgwopZAPmTkcxnkKCO3BZDQs088/RhykpFHRjLyP5AE70uadz/Z/1iK/2TgYwF9klJ8K4GPeJrTBXDsSXMeS3O2RA69wNEaOMAD1wTgMgJceYCLVZobH3IVBS68wLUauLwDLQKgEYG0O4CmCtC6ARpEQBsKaXYBLTWgcQe0B4EmJNLqBBqqQNsWaA4jLWig0Q2005GmPTAaAAYQyJgDGKYgIxtgMISMn4AhFzJKAwZ2yFgQGT4iI86srKysrD/rC4LWcCSWwIp+AAAAJXRFWHRkYXRlOmNyZWF0ZQAyMDIwLTA4LTE1VDA3OjAzOjMyKzAxOjAwzyqh1gAAACV0RVh0ZGF0ZTptb2RpZnkAMjAyMC0wOC0xNVQwNzowMzozMiswMTowML53GWoAAAAASUVORK5CYII=\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAMgAAADIAgMAAADQNkYNAAAABGdBTUEAALGPC/xhBQAAACBjSFJNAAB6JgAAgIQAAPoAAACA6AAAdTAAAOpgAAA6mAAAF3CculE8AAAADFBMVEUAAP9yEhL//////wAZPRNOAAAAAWJLR0QCZgt8ZAAAAAd0SU1FB+QIGQ86Nr0vrD0AAAI7SURBVGje7ZmxkcIwEEUhcAnuhxIIEIGvAwioxiU4QASU4GougA4gsM6SfNxh766sfx7N3Ix+zJv9X7KxtLtaZWVlZWWl1Fr12sQQpXLaRRMRTKFe2kbk+Na8POpNkbbmWlMqtkwxRsJllIotMxQ56F7neWX8LmrdGtNpPWc//Z4cLsbpfp6xN85XNRA904Sd2V/stXlJ16EFcL4O7Q/SnUPOrK/9xfzSvQ44K0dFhjLCmq0nRYYyvDPrq2rfka4RnVlfoyJ9GdGZ9XUdI89aWGYbpTITNUKYgvLlnXFhSsqXd7bjEcKXc8YhvYEjhdzY/Db9lUKebP4+/b6lkK7m8nNRhDAlE8WFoREuig/DpKej+DAbGmGiuDAUUsgItWR8ejZ/yaZ3+SmET+/zL4JIC8Ys2VpI7/JPkSKEbCnkyiNPCilDyI5A+AVzSzZFVAhRBCKssVvlBRB5W8iNwZCThDymiLyT5F5iSCsh3RSRN5/cfgwxoiikkpFmgqgwopZAPmTkcxnkKCO3BZDQs088/RhykpFHRjLyP5AE70uadz/Z/1iK/2TgYwF9klJ8K4GPeJrTBXDsSXMeS3O2RA69wNEaOMAD1wTgMgJceYCLVZobH3IVBS68wLUauLwDLQKgEYG0O4CmCtC6ARpEQBsKaXYBLTWgcQe0B4EmJNLqBBqqQNsWaA4jLWig0Q2005GmPTAaAAYQyJgDGKYgIxtgMISMn4AhFzJKAwZ2yFgQGT4iI86srKysrD/rC4LWcCSWwIp+AAAAJXRFWHRkYXRlOmNyZWF0ZQAyMDIwLTA4LTI1VDE0OjU4OjU0KzAxOjAwkH0c4gAAACV0RVh0ZGF0ZTptb2RpZnkAMjAyMC0wOC0yNVQxNDo1ODo1NCswMTowMOEgpF4AAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, - "execution_count": 51, + "execution_count": 47, "metadata": {}, "output_type": "execute_result" } @@ -1537,17 +1537,17 @@ }, { "cell_type": "code", - "execution_count": 52, + "execution_count": 48, "metadata": {}, "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAMgAAADIAgMAAADQNkYNAAAABGdBTUEAALGPC/xhBQAAACBjSFJNAAB6JgAAgIQAAPoAAACA6AAAdTAAAOpgAAA6mAAAF3CculE8AAAADFBMVEUAAP9yEhL//////wAZPRNOAAAAAWJLR0QCZgt8ZAAAAAd0SU1FB+QIDwgDItlM6zwAAAI7SURBVGje7ZmxkcIwEEUhcAnuhxIIEIGvAwioxiU4QASU4GougA4gsM6SfNxh766sfx7N3Ix+zJv9X7KxtLtaZWVlZWWl1Fr12sQQpXLaRRMRTKFe2kbk+Na8POpNkbbmWlMqtkwxRsJllIotMxQ56F7neWX8LmrdGtNpPWc//Z4cLsbpfp6xN85XNRA904Sd2V/stXlJ16EFcL4O7Q/SnUPOrK/9xfzSvQ44K0dFhjLCmq0nRYYyvDPrq2rfka4RnVlfoyJ9GdGZ9XUdI89aWGYbpTITNUKYgvLlnXFhSsqXd7bjEcKXc8YhvYEjhdzY/Db9lUKebP4+/b6lkK7m8nNRhDAlE8WFoREuig/DpKej+DAbGmGiuDAUUsgItWR8ejZ/yaZ3+SmET+/zL4JIC8Ys2VpI7/JPkSKEbCnkyiNPCilDyI5A+AVzSzZFVAhRBCKssVvlBRB5W8iNwZCThDymiLyT5F5iSCsh3RSRN5/cfgwxoiikkpFmgqgwopZAPmTkcxnkKCO3BZDQs088/RhykpFHRjLyP5AE70uadz/Z/1iK/2TgYwF9klJ8K4GPeJrTBXDsSXMeS3O2RA69wNEaOMAD1wTgMgJceYCLVZobH3IVBS68wLUauLwDLQKgEYG0O4CmCtC6ARpEQBsKaXYBLTWgcQe0B4EmJNLqBBqqQNsWaA4jLWig0Q2005GmPTAaAAYQyJgDGKYgIxtgMISMn4AhFzJKAwZ2yFgQGT4iI86srKysrD/rC4LWcCSWwIp+AAAAJXRFWHRkYXRlOmNyZWF0ZQAyMDIwLTA4LTE1VDA3OjAzOjM0KzAxOjAwrPqU7AAAACV0RVh0ZGF0ZTptb2RpZnkAMjAyMC0wOC0xNVQwNzowMzozNCswMTowMN2nLFAAAAAASUVORK5CYII=\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAMgAAADIAgMAAADQNkYNAAAABGdBTUEAALGPC/xhBQAAACBjSFJNAAB6JgAAgIQAAPoAAACA6AAAdTAAAOpgAAA6mAAAF3CculE8AAAADFBMVEUAAP9yEhL//////wAZPRNOAAAAAWJLR0QCZgt8ZAAAAAd0SU1FB+QIGQ86Nr0vrD0AAAI7SURBVGje7ZmxkcIwEEUhcAnuhxIIEIGvAwioxiU4QASU4GougA4gsM6SfNxh766sfx7N3Ix+zJv9X7KxtLtaZWVlZWWl1Fr12sQQpXLaRRMRTKFe2kbk+Na8POpNkbbmWlMqtkwxRsJllIotMxQ56F7neWX8LmrdGtNpPWc//Z4cLsbpfp6xN85XNRA904Sd2V/stXlJ16EFcL4O7Q/SnUPOrK/9xfzSvQ44K0dFhjLCmq0nRYYyvDPrq2rfka4RnVlfoyJ9GdGZ9XUdI89aWGYbpTITNUKYgvLlnXFhSsqXd7bjEcKXc8YhvYEjhdzY/Db9lUKebP4+/b6lkK7m8nNRhDAlE8WFoREuig/DpKej+DAbGmGiuDAUUsgItWR8ejZ/yaZ3+SmET+/zL4JIC8Ys2VpI7/JPkSKEbCnkyiNPCilDyI5A+AVzSzZFVAhRBCKssVvlBRB5W8iNwZCThDymiLyT5F5iSCsh3RSRN5/cfgwxoiikkpFmgqgwopZAPmTkcxnkKCO3BZDQs088/RhykpFHRjLyP5AE70uadz/Z/1iK/2TgYwF9klJ8K4GPeJrTBXDsSXMeS3O2RA69wNEaOMAD1wTgMgJceYCLVZobH3IVBS68wLUauLwDLQKgEYG0O4CmCtC6ARpEQBsKaXYBLTWgcQe0B4EmJNLqBBqqQNsWaA4jLWig0Q2005GmPTAaAAYQyJgDGKYgIxtgMISMn4AhFzJKAwZ2yFgQGT4iI86srKysrD/rC4LWcCSWwIp+AAAAJXRFWHRkYXRlOmNyZWF0ZQAyMDIwLTA4LTI1VDE0OjU4OjU0KzAxOjAwkH0c4gAAACV0RVh0ZGF0ZTptb2RpZnkAMjAyMC0wOC0yNVQxNDo1ODo1NCswMTowMOEgpF4AAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, - "execution_count": 52, + "execution_count": 48, "metadata": {}, "output_type": "execute_result" }