diff --git a/examples/jupyter/pincell.ipynb b/examples/jupyter/pincell.ipynb index 7d0f4e65a..06ac8dc58 100644 --- a/examples/jupyter/pincell.ipynb +++ b/examples/jupyter/pincell.ipynb @@ -173,7 +173,7 @@ "name": "stderr", "output_type": "stream", "text": [ - "/home/sam/openmc/openmc/openmc/mixin.py:71: IDWarning: Another Material instance already exists with id=2.\n", + "/home/master/anaconda3/envs/OMC/lib/python3.8/site-packages/openmc/mixin.py:68: IDWarning: Another Material instance already exists with id=2.\n", " warn(msg, IDWarning)\n" ] } @@ -218,7 +218,7 @@ "metadata": {}, "outputs": [], "source": [ - "mats = openmc.Materials([uo2, zirconium, water])" + "materials = openmc.Materials([uo2, zirconium, water])" ] }, { @@ -245,10 +245,10 @@ } ], "source": [ - "mats = openmc.Materials()\n", - "mats.append(uo2)\n", - "mats += [zirconium, water]\n", - "isinstance(mats, list)" + "materials = openmc.Materials()\n", + "materials.append(uo2)\n", + "materials += [zirconium, water]\n", + "isinstance(materials, list)" ] }, { @@ -294,7 +294,7 @@ } ], "source": [ - "mats.export_to_xml()\n", + "materials.export_to_xml()\n", "!cat materials.xml" ] }, @@ -353,7 +353,7 @@ "water.remove_nuclide('O16')\n", "water.add_element('O', 1.0)\n", "\n", - "mats.export_to_xml()\n", + "materials.export_to_xml()\n", "!cat materials.xml" ] }, @@ -386,14 +386,14 @@ "text": [ "\n", "\n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", - " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", " ...\n", " \n", " \n", @@ -493,7 +493,7 @@ "metadata": {}, "outputs": [], "source": [ - "sph = openmc.Sphere(r=1.0)" + "sphere = openmc.Sphere(r=1.0)" ] }, { @@ -511,8 +511,8 @@ "metadata": {}, "outputs": [], "source": [ - "inside_sphere = -sph\n", - "outside_sphere = +sph" + "inside_sphere = -sphere\n", + "outside_sphere = +sphere" ] }, { @@ -555,7 +555,7 @@ "outputs": [], "source": [ "z_plane = openmc.ZPlane(z0=0)\n", - "northern_hemisphere = -sph & +z_plane" + "northern_hemisphere = -sphere & +z_plane" ] }, { @@ -663,7 +663,7 @@ { "data": { "text/plain": [ - "" + "" ] }, "execution_count": 24, @@ -672,7 +672,7 @@ }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAARUAAAD8CAYAAABZ0jAcAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAATXElEQVR4nO3dfbBcdX3H8fenUEJrp5IQBiMqkIoGqh3AK9oyY1tFRP9IsKJGpxUcHGpb2pk6OsIwUzuoLdg/sE6tmsEHfBiC0nEaRxzksf6hQW4skgCFhDjVpCiRIJ1OMJrw7R/nt/HkZnfv3t3fnqf9vGZ27u552Ps79/zO5/7O2YevIgIzs1x+re4GmFm3OFTMLCuHipll5VAxs6wcKmaWlUPFzLLKEiqSPiPpcUnbBsyXpI9J2iHpfklnl+ZdLGl7ul2coz1mVp9cI5XPARcMmf964LR0uwz4BICkFcAHgFcA5wAfkLQ8U5vMrAZZQiUivgXsHbLIOuDzUdgMHCdpFfA64LaI2BsRTwK3MTyczKzhjq7o95wE/Kj0eFeaNmj6ESRdRjHK4VnPetbL1qxZM52W2sh2bntwKs+7+iVnTOV5bXRbtmz5aUScMM66VYXKxCJiA7ABYG5uLubn52tuUfetf9HZQ+e/8AVTCvb/fWbo7I2PfG86v9cOkfTf465bVajsBp5fevy8NG038EcLpt9dUZusZLEAaZJ+bXXQNEdVobIJuFzSRoqLsk9FxGOSbgX+oXRx9nzgyoraNNPaFCKjWLg9Dpn6ZAkVSTdSjDhWStpF8YrOrwNExCeBW4A3ADuAfcA707y9kj4I3Jue6uqIGHbB18bQtQAZhUcz9VEbv/rA11QWN4tBshQOmOEkbYmIuXHW9TtqzSyr1rz6Y6PxCGU0vb+TRyz5OVQ6wEEyvvLfzgGTh0OlpRwk+Tlg8nCotIiDpDoOmPE5VFrAYVIvX39ZGodKgzlMmsXhMhq/pNxQDpTm8r4ZziOVBnFnbQ+PWgZzqDSAw6S9fEH3SD79qZkDpTu8LwseqdTEHbCbfFrkkYqZZeaRSsU8QpkNszxi8UilQg6U2TOL+9yhUpFZ7FxWmLV979OfKZu1DmX9zdLpkENlShwm1s8shEuusqcXSHo4lTW9os/86yTdl26PSPpZad7B0rxNOdpTNweKLabLfWTikYqko4CPA6+lKAZ2r6RNEXGo0lRE/G1p+b8Gzio9xdMRceak7WiKLncWy2v9i87u5Iglx+nPOcCOiNgJkMpwrAMGla97G8W37XeKw8TG0cXToRynP0spXXoycCpwZ2nysZLmJW2WdGGG9lTOgWKT6lIfqvol5fXAzRFxsDTt5FQK4O3ARyX9Tr8VJV2Wwmd+z549VbR1JF3qDFavrvSlHKEyqKRpP+uBG8sTImJ3+rmTouTpWUeuVtRSjoi5iJg74YSx6kZn15VOYM3RhT6VI1TuBU6TdKqkYyiC44hXcSStAZYD3ylNWy5pWbq/EjiXwddizKwFJg6ViDgAXA7cCjwEfDkiHpB0taS1pUXXAxvj8JKIpwPzkr4P3AVcU37VqMm68B/FmqntfctlT8fQ9p1u7VDnK0Iue1ohB4pVpa19zW/TH1Fbd7C1Wxvfx+KRyggcKFa3NvVBh8oi2rQzrdva0hcdKkO0ZSfa7GhDn3SomFlWDpUB2vAfwWZT0/umQ6WPpu80syb3UYfKAk3eWWZlTe2rDhUzy8qhUtLU5DcbpIl91qFiZlk5VJImJr7ZKJrWdx0qNG+nmC1Vk/rwzIdKk3aG2SSa0pdnPlTMLK+Z/eqDpqS6WU5N+KoEj1TMLKuZDBWPUqzr6uzjVdVSvkTSnlLN5HeV5l0saXu6XZyjPcM4UGxW1NXXK6mlnNwUEZcvWHcFRQnUOSCALWndJydtl5nVI8dI5VAt5Yj4BdCrpTyK1wG3RcTeFCS3ARdkaFNfHqXYrKmjz1dZS/lNku6XdLOkXkXDpdRhbmTZUzM7XFUXar8GnBIRv0cxGrlhqU/QxLKnZnakSmopR8QTEbE/PbweeNmo6+biUx+bVVX3/UpqKUtaVXq4lqI8KhSlUs9PNZWXA+enaVk5UGzWVXkMTPzqT0QckNSrpXwU8JleLWVgPiI2AX+T6iofAPYCl6R190r6IEUwAVwdEXsnbZOZ1afztZQ9SjH7lVHfvu9aymbWGJ0OFY9SzA5XxTHR6VAxs+o5VMwsq86Gik99zPqb9rHR2VAxs3o4VMwsq06Gik99zIab5jHSyVAxs/o4VMwsq86Fik99zEYzrWOlc6FiZvXqTN0fj1DMlm4adYI8UjGzrBwqZpZVJ0LFpz5mk8l5DHUiVMysORwqZpZVVWVP3yPpwVT35w5JJ5fmHSyVQ920cF0za5eqyp7+JzAXEfsk/QXwEeCtad7TEXHmpO0ws2aopOxpRNwVEfvSw80U9X2y8EVaszxyHUtVlj3tuRT4Runxsamc6WZJFw5ayWVPzdqh0nfUSvpTYA74w9LkkyNit6TVwJ2StkbEowvXjYgNwAYoSnRU0mAzW7JKyp4CSDoPuApYWyqBSkTsTj93AncDZ2Vok5nVpKqyp2cBn6IIlMdL05dLWpburwTOBcoXeIfy9RSzvHIcU1WVPf0n4LeAr0gC+GFErAVOBz4l6RmKgLtmwatGZtYyWa6pRMQtwC0Lpv1d6f55A9b7NvDSHG0ws2bwO2rNLCuHipll1dpQ8UVas+mY9NhqbaiYWTM5VMwsK4eKmWXlUDGzrBwqZpaVQ8XMsnKomFlWrQyVndv88SCzaVq+7DdeNu66rQwVM2suh4qZZeVQMbOsHCpmlpVDxcyycqiYWVYOFTPLqqqyp8sk3ZTm3yPplNK8K9P0hyW9Lkd7zKw+E4dKqezp64EzgLdJOmPBYpcCT0bEC4HrgGvTumdQfPv+7wIXAP+ans/MWqqSsqfp8Q3p/s3Aa1R8rf46YGNE7I+IHwA70vOZWUtVVfb00DIRcQB4Cjh+xHWBw8ue7j94IEOzzWwaWnOhNiI2RMRcRMwtO6rSaq1mtgRVlT09tIyko4FnA0+MuK6ZtUglZU/T44vT/YuAOyMi0vT16dWhU4HTgO9maJOZ1aSqsqefBr4gaQewlyJ4SMt9maJ+8gHgryLi4KRtMrP6qBgwtMuKY38zzn/BmrqbYdZZ3/zhf7H35/s0zrqtuVBrZu3gUDGzrBwqZpaVQ8XMsnKomFlWrQyV1S9Z+HlFM8vpyf1Pbxl33VaGipk1l0PFzLJyqJhZVg4VM8vKoWJmWTlUzCwrh4qZZdXaUNn4yPfqboJZJ016bLU2VMysmRwqZpaVQ8XMsnKomFlWE4WKpBWSbpO0Pf1c3meZMyV9R9IDku6X9NbSvM9J+oGk+9LtzKX8fl+sNcsrxzE16UjlCuCOiDgNuCM9Xmgf8I6I6JU2/aik40rz3xcRZ6bbfRO2x8xqNmmolMuZ3gBcuHCBiHgkIran+/8DPA6cMOHvNbOGmjRUToyIx9L9HwMnDltY0jnAMcCjpckfTqdF10laNmTdQ2VP9+zZM2GzzWxaFg0VSbdL2tbndlgR9lQcbGC9D0mrgC8A74yIZ9LkK4E1wMuBFcD7B61fLnt6wgm/Guj4uopZHrmOpUWLiUXEeYPmSfqJpFUR8VgKjccHLPfbwNeBqyJic+m5e6Oc/ZI+C7x3Sa03s8aZ9PSnXM70YuDfFy6QSqF+Ffh8RNy8YN6q9FMU12O2TdgeM6vZpKFyDfBaSduB89JjJM1Juj4t8xbgVcAlfV46/pKkrcBWYCXwoQnbY2Y1m6iWckQ8Abymz/R54F3p/heBLw5Y/9WT/H4za55OvKPWF2vNJpPzGOpEqJhZczhUzCyria6pNElv+Lb+RWfX3BKz9pjGpQOPVMwsq86Fii/amo1mWsdK50LFzOrlUDGzrDoZKj4FMhtumsdIJ0PFzOrjUDGzrDobKj4FMutv2sdGZ0PFzOrhUDGzrDodKj4FMjtcFcdEp0PFzKrX+VDxaMWsUNWx0PlQMbNqTb3saVruYOn7aTeVpp8q6R5JOyTdlL4kOzuPVmzWVXkMVFH2FODpUmnTtaXp1wLXRcQLgSeBSydsz0AOFptVVff9qZc9HSSV5Xg10CvbsaT1zayZqip7emwqWbpZUi84jgd+FhEH0uNdwEmDflGOsqcerdisqaPPL/p1kpJuB57TZ9ZV5QcREZIGlT09OSJ2S1oN3Jlq/Ty1lIZGxAZgA8Dc3NzA8qpmVq9Kyp5GxO70c6eku4GzgH8DjpN0dBqtPA/YPcY2mFmDVFH2dLmkZen+SuBc4MFU0P0u4KJh6+fmUyCbFXX19SrKnp4OzEv6PkWIXBMRD6Z57wfeI2kHxTWWT0/YnpE4WKzr6uzjVZQ9/Tbw0gHr7wTOmaQNZtYsnan7s1SuE2Rd1IRRuN+mb2ZZzXyoNCHZzXJoSl+e+VCB5uwMs3E1qQ87VJIm7RSzpWha33WomFlWDpWSpiW+2WKa2GcdKmaWlUNlgSYmv1k/Te2rDpU+mrqzzHqa3EcdKgM0eafZbGt633SomFlWDpUhmv4fwWZPG/qkQ2URbdiJNhva0hcdKiNoy8607mpTH5zZrz5YKn9VgtWhTWHS45HKErVxJ1s7tbWvOVTG0Nadbe3R5j429bKnkv64VPL0Pkk/79X+kfQ5ST8ozTtzkvZUqc073Zqt7X1r6mVPI+KuXslTioqE+4BvlhZ5X6kk6n0TtsfMalZ12dOLgG9ExL4Jf28jtP0/ijVPF/pUVWVPe9YDNy6Y9mFJ90u6rlcfqE260AmsGbrSlxYNFUm3S9rW57auvFwqDjawHGmqYPhS4NbS5CuBNcDLgRUUdYAGrT9xLeVp6UpnsPp0qQ9VUvY0eQvw1Yj4Zem5e6Oc/ZI+C7x3SDsaXUvZ72OxcXQpTHqmXva05G0sOPVJQYQkUVyP2TZhe2rXxU5i09HVvlJF2VMknQI8H/iPBet/SdJWYCuwEvjQhO1phK52Fsuny31ExaWQdpmbm4v5+fm6mzESnw5ZWVvCRNKWiJgbZ11/9mfKfK3FoD1hkoPfpl+RWepUdrhZ2/cOlQrNWuey2dznPv2pmE+HZsMshkmPRypmlpVHKjXxiKWbZnmE0uORSs3cCbvD+7LgkUoDeNTSXg6SIzlUGqTcQR0wzeYwGcynPw3lTttc3jfDeaTSYD4tahaHyWgcKi3gcKmXw2RpHCot4msu1XGQjM+h0lIOmPwcJHk4VDrAATM+B0l+DpWO8fWX0ThMpscvKZtZVh6pdJRPiY7k0Uk1HCozoN/B1PWgcYDUZ6JQkfRm4O+B04FzIqLvF8dKugD4Z+Ao4PqI6H1B9qnARuB4YAvwZxHxi0naZKNZeNC1PWQcIs0x6UhlG/AnwKcGLSDpKODjwGuBXcC9kjZFxIPAtcB1EbFR0ieBS4FPTNgmG0ObRjMOkGabKFQi4iGAomzPQOcAOyJiZ1p2I7BO0kMUBdvfnpa7gWLU41BpiMUO3mmFjkOj3aq4pnIS8KPS413AKyhOeX4WEQdK008a9CSSLgMuSw/3S2p94bE+VgI/rbsRUzLytt00/J9U03R1n7143BUXDRVJtwPP6TPrqogYVpEwq3LZU0nz49YkabKubhd0d9u6vF3jrjtRLeUR7aaoTtjzvDTtCeA4SUen0Upvupm1WBVvfrsXOE3SqZKOAdYDm6IojXgXcFFabrFazGbWAhOFiqQ3StoF/D7wdUm3punPlXQLQBqFXA7cCjwEfDkiHkhP8X7gPZJ2UFxj+fSIv3rDJO1usK5uF3R327xdC7SylrKZNZc/+2NmWTlUzCyrVoSKpDdLekDSM5IGvnwn6QJJD0vaIemKKts4DkkrJN0maXv6uXzAcgcl3Zdum6pu56gW+/tLWibppjT/HkmnVN/K8YywbZdI2lPaT++qo51LIekzkh4f9J4vFT6Wtvl+SaO92zEiGn+j+GzRi4G7gbkByxwFPAqsBo4Bvg+cUXfbF9mujwBXpPtXANcOWO7/6m7rCNuy6N8f+Evgk+n+euCmutudcdsuAf6l7rYucbteBZwNbBsw/w3ANwABrwTuGeV5WzFSiYiHIuLhRRY79HGAKD6UuBFYN/3WTWQdxccTSD8vrLEtkxrl71/e3puB12iRz3g0RBv71qIi4lvA3iGLrAM+H4XNFO8rW7XY87YiVEbU7+MAA9/23xAnRsRj6f6PgRMHLHespHlJmyU1NXhG+fsfWiaKtxo8RfFWgqYbtW+9KZ0m3Czp+X3mt81Yx1Rjvk+lKR8HyG3YdpUfRERIGvT6/skRsVvSauBOSVsj4tHcbbWJfA24MSL2S/pzihHZq2tuUy0aEyoxvY8D1GrYdkn6iaRVEfFYGlY+PuA5dqefOyXdDZxFcY7fJKP8/XvL7JJ0NPBsio9rNN2i2xYR5e24nuJ6WduNdUx16fSn78cBam7TYjZRfDwBBnxMQdJyScvS/ZXAucCDlbVwdKP8/cvbexFwZ6Qrgg236LYtuNawluLd4223CXhHehXolcBTpdP1weq+Aj3iVeo3UpzP7Qd+Atyapj8XuGXB1epHKP6LX1V3u0fYruOBO4DtwO3AijR9juIb8gD+ANhK8YrDVuDSuts9ZHuO+PsDVwNr0/1jga8AO4DvAqvrbnPGbftH4IG0n+4C1tTd5hG26UbgMeCX6fi6FHg38O40XxRfsPZo6nt9X3ldePPb9M0sqy6d/phZAzhUzCwrh4qZZeVQMbOsHCpmlpVDxcyycqiYWVb/D5MtR6Xki1ByAAAAAElFTkSuQmCC\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAARUAAAD8CAYAAABZ0jAcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy86wFpkAAAACXBIWXMAAAsTAAALEwEAmpwYAAASQElEQVR4nO3de6wc5X3G8e9TglXFJeV+CaEhiSwIaQuFIy6lSqAJCFBbh6qh0IpEKJFFhaUmqqpSIVH+TKnSqFQJ1E1JQGqhVMGJBeauppRGJNgI3zDEhjqKsYUNQbgpbajbX//YWbM+3j1ndufdmXdmn4+0OrtzOfvOmXee887s5aeIwMwslZ9pugFm1i0OFTNLyqFiZkk5VMwsKYeKmSXlUDGzpJKEiqQ7Je2RtHnEfEm6TdJ2SRslnT0w7zJJLxbzbkzRHjNrTqqRyjeAyxaYfzmwrLitAG4HkHQY8JVi/hnANZLOSNQmM2tAklCJiCeBHy+wyHLg7uh5GjhS0knAucD2iHg5It4G7i2WNbOWeldNz3My8KOBxzuLacOmnzfsF0haQW+Uw9KlS885/fTTp9NSK23bGzum8nuXHXXqVH6vlbd+/frXIuK4SdatK1Q0ZFosMP3QiRGrgFUAc3NzsW7dunSts6Euu++6Becfyy/V1JKDPXzV1xt53lki6YeTrltXqOwEThl4/D5gF7BkxHSr2WIBkpNhbXXQ5KOuUFkDrJR0L73TmzcjYrekvcAySR8AXgGuBn6vpjbNtDaFSBnzt8ch05wkoSLpHuAi4FhJO4E/Aw4HiIg7gLXAFcB24C3gumLefkkrgUeAw4A7I2JLijbZO7oWIGV4NNMctfGrD3xNZXGzGCTjcMAsTNL6iJibZF2/o9bMkqrrmorVxCOUcvp/J49Y0nOodICDZHKDfzsHTBoOlZZykKTngEnDodIiDpL6OGAm51BpAYdJs3z9ZTwOlYw5TPLicCnHLylnyoGSL++bhXmkkhF31vbwqGU0h0oGHCbt5Qu6h/LpT8McKN3hfdnjkUpD3AG7yadFHqmYWWIeqdTMI5TZMMsjFo9UauRAmT2zuM8dKjWZxc5lPbO27336M2Wz1qFsuFk6HXKoTInDxIaZhXBJVfZ0wdKlkv5Y0nPFbbOk/5V0dDFvh6RNxbxOfEekA8UW0+U+UjlUypQujYi/iIizIuIs4E+Bf4mIwYqGFxfzJ/pOzJx0ubNYWl3tKylOfw6ULgUoynAsB54fsfw1wD0JnjcrXe0gNl1dPB1KcfozqqTpISS9m14h928OTA7gUUnri9KmreNAsaq61IdShErp0qXAbwL/Nu/U58KIOJve6dMNkj469EmkFZLWSVq3d+/eai1OqEudwZrVlb6UIlRGlTQd5mrmnfpExK7i5x5gNb3TqUNExKqImIuIueOOm6hudHJd6QSWjy70qRSh8gxF6VJJS+gFx5r5C0n6eeBjwLcHpi2VdET/PnApsDlBm8ysIZVDJSL2A/3SpVuB+yJii6TrJV0/sOiVwKMR8Z8D004AnpK0Afg+8GBEPFy1TXXown8Uy1Pb+5bLnk6g7Tvd2qHJV4Rc9rRGDhSrS1v7mt+mX1Jbd7C1Wxvfx+KRSgkOFGtam/qgQ2URbdqZ1m1t6YsOlQW0ZSfa7GhDn3SomFlSDpUR2vAfwWZT7n3ToTJE7jvNLOc+6lCZJ+edZTYo177qUDGzpBwqA3JNfrNRcuyzDhUzS8qhUsgx8c3KyK3vOlTIb6eYjSunPjzzoZLTzjCrIpe+PPOhYmZpzexXH+SS6mYp5fBVCR6pmFlSMxkqHqVY1zXZx+uqpXyRpDcH6infXHbd1BwoNiua6uuVr6kM1FK+hF4NoGckrYmI+WVP/zUifmPCdc2sJVKMVA7UUo6It4F+LeVprzs2j1Js1jTR5+uspXyBpA2SHpL0kTHXzbbsqZkdrK5ays8C74+IM4G/Br41xrq9iRmWPTWzQ9VSSzki9kXET4r7a4HDJR1bZt1UfOpjs6ruvl9LLWVJJ0pScf/c4nlfL7NuCg4Um3V1HgOVX/2JiP2S+rWUDwPu7NdSLubfAfwO8AeS9gP/BVwdvXqrQ9et2iYza07nayl7lGL2jrJv33ctZTPLRqdDxaMUs4PVcUx0OlTMrH4OFTNLqrOh4lMfs+GmfWx0NlTMrBkOFTNLqpOh4lMfs4VN8xjpZKiYWXMcKmaWVOdCxac+ZuVM61jpXKiYWbM6U/fHIxSz8U2jTpBHKmaWlEPFzJLqRKj41MesmpTHUCdCxczy4VAxs6TqKnv6+5I2FrfvSjpzYN4OSZuKcqjlviPSzLJVV9nTfwc+FhFvSLocWAWcNzD/4oh4rWpbzKx5tZQ9jYjvRsQbxcOn6dX3ScIXac3SSHUs1Vn2tO+zwEMDjwN4VNJ6SStGreSyp2btkOIdtaVLl0q6mF6o/NrA5AsjYpek44HHJL0QEU8e8gsjVtE7bWJubq59dUXMZkQtZU8BJP0y8DVgeUS83p8eEbuKn3uA1fROp8yspeoqe/oLwP3AtRHxg4HpSyUd0b8PXApsLvvEvp5illaKY6qusqc3A8cAXy1KKu8vqp+dAKwupr0L+IeIeLhqm8ysOUk+pRwRa4G186bdMXD/c8Dnhqz3MnDm/Olm1l5+R62ZJeVQMbOkWhsqvkhrNh1Vj63WhoqZ5cmhYmZJOVTMLCmHipkl5VAxs6QcKmaWlEPFzJJqZahse2NH000w67T3fPCYcyZdt5WhYmb5cqiYWVIOFTNLyqFiZkk5VMwsKYeKmSXlUDGzpOoqeypJtxXzN0o6u+y6ZtYulUNloOzp5cAZwDWSzpi32OXAsuK2Arh9jHXNrEVqKXtaPL47ep4GjpR0Usl1zaxF6ip7OmqZ0iVTB8uevr3vvys32symI0WolCl7OmqZ0iVTI2JVRMxFxNyS9/zsmE00s7qkqPtTpuzpqGWWlFjXzFqklrKnxeNPF68CnQ+8GRG7S65rZi1SV9nTtcAVwHbgLeC6hdat2iYza05dZU8DuKHsumbWXn5HrZkl5VAxs6QcKmaWlEPFzJJyqJhZUq0MlWVHndp0E8w6bd/Lr6+fdN1WhoqZ5cuhYmZJOVTMLCmHipkl5VAxs6QcKmaWlEPFzJJqbag8fNXXm26CWSdVPbZaGypmlieHipkl5VAxs6QcKmaWVKVQkXS0pMckbSt+HjVkmVMk/bOkrZK2SPrDgXm3SHpF0nPF7Ypxnt8Xa83SSnFMVR2p3Ag8ERHLgCeKx/PtB/4oIj4MnA/cMK+06Zcj4qzi5u+qNWu5qqGyHLiruH8X8Mn5C0TE7oh4trj/H8BWRlQhNLP2qxoqJxT1eyh+Hr/QwpJOBX4F+N7A5JWSNkq6c9jp08C6B8qe7t27t2KzzWxaFg0VSY9L2jzkNlYhdUk/B3wT+HxE7Csm3w58CDgL2A18adT6g2VPjzvuuAPTfV3FLI1Ux9KidX8i4hOj5kl6VdJJEbFb0knAnhHLHU4vUP4+Iu4f+N2vDizzt8AD4zTezPJT9fRnDfCZ4v5ngG/PX0CSgL8DtkbEX86bd9LAwyuBzRXbY2YNqxoqXwQukbQNuKR4jKT3Suq/knMhcC3w60NeOr5V0iZJG4GLgS9UbI+ZNaxS2dOIeB34+JDpu+jVTiYingI0Yv1rqzy/meWnE++o9cVas2pSHkOdCBUzy4dDxcySqnRNJSf94dtl913XcEvM2mMalw48UjGzpDoXKr5oa1bOtI6VzoWKmTXLoWJmSXUyVHwKZLawaR4jnQwVM2uOQ8XMkupsqPgUyGy4aR8bnQ0VM2uGQ8XMkup0qPgUyOxgdRwTnQ4VM6tf50PFoxWznrqOhc6HipnVa+plT4vldhTfRfucpHXjrl+VRys26+o8Buooe9p3cVHadG7C9StxsNisqrvvT73s6ZTXN7PM1FX2NIBHJa2XtGKC9ZOUPfVoxWZNE31+0a+TlPQ4cOKQWTeN8TwXRsQuSccDj0l6ISKeHGN9ImIVsApgbm4uxlnXzOpTS9nTog4QEbFH0mrgXOBJoNT6ZtYedZQ9XSrpiP594FLeKW+66Pqp+RTIZkVTfb2OsqcnAE9J2gB8H3gwIh5eaP1pc7BY1zXZx+soe/oycOY465tZe3Wm7s+4XCfIuiiHUbjfpm9mSc18qOSQ7GYp5NKXZz5UIJ+dYTapnPqwQ6WQ004xG0dufdehYmZJOVQG5Jb4ZovJsc86VMwsKYfKPDkmv9kwufZVh8oQue4ss76c+6hDZYScd5rNttz7pkPFzJJyqCwg9/8INnva0CcdKotow0602dCWvuhQKaEtO9O6q019cGa/+mBc/qoEa0KbwqTPI5UxtXEnWzu1ta85VCbQ1p1t7dHmPjb1sqeSTivKnfZv+yR9vph3i6RXBuZdUaU9dWrzTre8tb1vTb3saUS8WJQ7PQs4B3gLWD2wyJf78yNi7fz1zaxd6i57+nHgpYj4YcXnzULb/6NYfrrQp+oqe9p3NXDPvGkrJW2UdOew06fcdaETWB660pcWDRVJj0vaPOS2fJwnkrQE+C3gnwYm3w58CDgL2A18aYH1K9dSnpaudAZrTpf6kCImL0ss6UXgooGypd+JiNNGLLscuCEiLh0x/1TggYj4xcWed25uLtatWzdxu6fJ72OxceQaJpLWR8TcJOtOvezpgGuYd+pTBFHflbxTDrW1cu0klp+u9pU6yp4i6d3F/PvnrX+rpE2SNgIXA1+o2J4sdLWzWDpd7iOVTn+akvPpz3w+HbJBbQmTKqc//uzPlPkzQwbtCZMU/Db9msxSp7KDzdq+d6jUaNY6l83mPvfpT818OjQbZjFM+jxSMbOkPFJpiEcs3TTLI5Q+j1Qa5k7YHd6XPR6pZMCjlvZykBzKoZKRwQ7qgMmbw2Q0n/5kyp02X943C/NIJWM+LcqLw6Qch0oLOFya5TAZj0OlRXzNpT4Oksk5VFrKAZOegyQNh0oHOGAm5yBJz6HSMb7+Uo7DZHr8krKZJeWRSkf5lOhQHp3Uw6EyA4YdTF0PGgdIcyqFiqRPAbcAHwbOjYihXxwr6TLgr4DDgK9FRP8Lso8G/hE4FdgBXBURb1Rpk5Uz/6Bre8g4RPJRdaSyGfht4G9GLSDpMOAr9L5NfyfwjKQ1EfE879Ri/qKkG4vHf1KxTTaBNo1mHCB5qxQqEbEVQNJCi50LbI+Il4tl76VXg/n54udFxXJ3Ad/BoZKNxQ7eaYWOQ6Pd6rimcjLwo4HHO4HzivsH1WKWNLIWs6QVwIri4U8ltb7w2BDHAq813YgpKb1t+t1vTLclaXV1nw2tNFrGoqEi6XHgxCGzboqIhSoSHvgVQ6aNXWwoIlYBq4o2rZu0JknOurpd0N1t6/J2TbruoqESEZ+Y9JcXdgKnDDx+H7CruP+qpJMGajHvqfhcZtawOt789gywTNIHJC0BrqZXgxnGq8VsZi1QKVQkXSlpJ3AB8KCkR4rpB2opR8R+YCXwCLAVuC8ithS/Ymgt5hJWVWl3xrq6XdDdbfN2zdPKWspmli9/9sfMknKomFlSrQgVSZ+StEXS/0ka+fKdpMskvShpe/EO3axJOlrSY5K2FT+PGrHcDkmbJD1X5aW+aVvs76+e24r5GyWd3UQ7J1Fi2y6S9Gaxj56TdHMT7RyXpDsl7Rn1vq+J9llEZH+j99mi0+i943ZuxDKHAS8BHwSWABuAM5pu+yLbdStwY3H/RuDPRyy3Azi26fYusi2L/v2BK4CH6L136Xzge023O+G2XQQ80HRbJ9i2jwJnA5tHzB97n7VipBIRWyPixUUWO/BxgIh4G+h/HCBny+l9PIHi5yeba0plZf7+y4G7o+dp4Mji/Um5a2PfKiUingR+vMAiY++zVoRKScM+DnByQ20p66CPKQCjPqYQwKOS1hcfV8hRmb9/G/cRlG/3BZI2SHpI0kfqadrUjb3Psvk+lVw+DpDaQts1xq+5MCJ2FZ+NekzSC8V/mJyU+ftnuY9KKNPuZ4H3R8RPJF0BfAtYNu2G1WDsfZZNqMR0Pw7QmIW2S1KpjylExK7i5x5Jq+kNx3MLlTJ//yz3UQmLtjsi9g3cXyvpq5KOjYi2f9hw7H3WpdOfhT4OkKtFP6YgaamkI/r3gUvpfY9Nbsr8/dcAny5eUTgfeLN/+pe5RbdN0okqvgNE0rn0jq3Xa29peuPvs6avPpe8Qn0lvcT8KfAq8Egx/b3A2nlXqn9A70r9TU23u8R2HQM8AWwrfh49f7voveKwobhtyXm7hv39geuB64v7oveFXS8BmxjxSl6OtxLbtrLYPxuAp4FfbbrNJbfrHmA38D/FMfbZqvvMb9M3s6S6dPpjZhlwqJhZUg4VM0vKoWJmSTlUzCwph4qZJeVQMbOk/h954lbyY/ECLgAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] @@ -702,7 +702,7 @@ { "data": { "text/plain": [ - "" + "" ] }, "execution_count": 25, @@ -711,7 +711,7 @@ }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAARUAAAD8CAYAAABZ0jAcAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAARZ0lEQVR4nO3de4xc5X3G8e9TXBuJqsHGDnUAGxxQjau0NmwhLVLaAuGSP2zS0MRUVWxk5KYtrVREBAipkUjSQprKUdS0iQUEiCIucRXVUUCuwab80ZiwVg02RuDFNIldBzs2IFWmTmx+/eO8Sw/rnd3ZnXfPZfb5SKOZc5t9z8yZZ99zmfkpIjAzy+WX6m6AmfUXh4qZZeVQMbOsHCpmlpVDxcyycqiYWVZZQkXS/ZIOStrVYbokfVXSkKQXJF1UmrZK0p50W5WjPWZWn1w9lQeAa8aYfi1wQbqtBf4ZQNIc4HPApcAlwOckzc7UJjOrQZZQiYhngCNjzLICeCgK24DTJc0HrgY2R8SRiHgD2MzY4WRmDTejor9zFvCT0vC+NK7T+JNIWkvRy+G00067ePHixVPTUuva4ef3TcnznvFbZ0/J81r3tm/f/rOImDeZZasKlZ5FxHpgPcDAwEAMDg7W3KL+98D7bx17hqnaUf3x2JNXH/zyFP1hGybpR5NdtqpQ2Q+cUxo+O43bD/z+iPFPV9QmKxk3QBpktLY6aJqjqlDZCNws6RGKg7JvRcQBSZuAvy0dnL0KuKOiNk1rbQqRboxcH4dMfbKEiqSHKXoccyXtozij88sAEfF14HHgY8AQcBS4MU07IunzwHPpqe6KiLEO+Nok9FuAdMO9mfqojT994GMq45uOQTIRDpixSdoeEQOTWdZX1JpZVq05+2PdcQ+lO8Ovk3ss+TlU+oCDZPLKr50DJg+HSks5SPJzwOThUGkRB0l1HDCT51BpAYdJvXz8ZWIcKg3mMGkWh0t3fEq5oRwozeX3ZmzuqTSIN9b2cK+lM4dKAzhM2ssHdE/m3Z+aOVD6h9/LgnsqNfEG2J+8W+Seipll5p5KxdxDmR6mc4/FPZUKOVCmn+n4njtUKjIdNy4rTLf33rs/U2y6bVA2uum0O+RQmSIOExvNdAiXXGVPr5H0ciprevso09dJ2pFur0h6szTtRGnaxhztqZsDxcbTz9tIzz0VSacAXwM+SlEM7DlJGyNi9/A8EfHXpfn/ElhWeoq3I2Jpr+1oin7eWCyvB95/a1/2WHLs/lwCDEXEXoBUhmMFsLvD/DdQ/Np+X3GY2GT04+5Qjt2fiZQuXQicB2wpjT5V0qCkbZKuy9CeyjlQrFf9tA1VfUp5JbAhIk6Uxi1MpQD+GPiKpA+OtqCktSl8Bg8dOlRFW7vSTxuD1atftqUcodKppOloVgIPl0dExP50v5ei5OmykxcrailHxEBEDMybN6m60dn1y0ZgzdEP21SOUHkOuEDSeZJmUgTHSWdxJC2mKOn9g9K42ZJmpcdzgcvofCzGzFqg51CJiOPAzcAm4CXgsYh4UdJdkpaXZl0JPBLvLYl4ITAo6XlgK3B3+axRk/XDfxRrprZvWy57Ogltf9OtHeo8I+SypxVyoFhV2rqt+TL9LrX1DbZ2a+N1LO6pdMGBYnVr0zboUBlHm95M629t2RYdKmNoy5to00cbtkmHipll5VDpoA3/EWx6avq26VAZRdPfNLMmb6MOlRGa/GaZlTV1W3WomFlWDpWSpia/WSdN3GYdKmaWlUMlaWLim3WjaduuQ4XmvSlmE9WkbXjah0qT3gyzXjRlW572oWJmeU3bnz5oSqqb5dSEn0pwT8XMspqWoeJeivW7Orfxqmopr5Z0qFQz+abStFWS9qTbqhztGYsDxaaLurb1SmopJ49GxM0jlp1DUQJ1AAhge1r2jV7bZWb1yNFTebeWckT8HBiupdyNq4HNEXEkBclm4JoMbRqVeyk23dSxzVdZS/kTkl6QtEHScEXDidRhbmTZUzN7r6oO1H4PODcifpOiN/LgRJ+giWVPzexkldRSjojDEXEsDd4LXNztsrl418emq6q3/UpqKUuaXxpcTlEeFYpSqVelmsqzgavSuKwcKDbdVfkZ6PnsT0QclzRcS/kU4P7hWsrAYERsBP4q1VU+DhwBVqdlj0j6PEUwAdwVEUd6bZOZ1afvaym7l2L2/7q9fN+1lM2sMfo6VNxLMXuvKj4TfR0qZlY9h4qZZdW3oeJdH7PRTfVno29Dxczq4VAxs6z6MlS862M2tqn8jPRlqJhZfRwqZpZV34WKd33MujNVn5W+CxUzq1ff1P1xD8Vs4qaiTpB7KmaWlUPFzLLqi1Dxro9Zb3J+hvoiVMysORwqZpZVVWVPb5G0O9X9eUrSwtK0E6VyqBtHLmtm7VJV2dP/BAYi4qikPwO+BHwqTXs7Ipb22g4za4ZKyp5GxNaIOJoGt1HU98nCB2nN8sj1Waqy7OmwNcATpeFTUznTbZKu67SQy56atUOlV9RK+hNgAPi90uiFEbFf0iJgi6SdEfHqyGUjYj2wHooSHZU02MwmrJKypwCSrgTuBJaXSqASEfvT/V7gaWBZhjaZWU2qKnu6DPgGRaAcLI2fLWlWejwXuAwoH+Adk4+nmOWV4zNVVdnTvwd+BfiOJIAfR8Ry4ELgG5LeoQi4u0ecNTKzlslyTCUiHgceHzHub0qPr+yw3H8AH8rRBjNrBl9Ra2ZZOVTMLKvWhooP0ppNjV4/W60NFTNrJoeKmWXlUDGzrBwqZpaVQ8XMsnKomFlWDhUzy6qVoXL4+X11N8Gsr50748yLJ7tsK0PFzJrLoWJmWTlUzCwrh4qZZeVQMbOsHCpmlpVDxcyyqqrs6SxJj6bpz0o6tzTtjjT+ZUlX52iPmdWn51AplT29FlgC3CBpyYjZ1gBvRMT5wDrgnrTsEopf3/8N4Brgn9LzmVlLVVL2NA0/mB5vAK5Q8bP6K4BHIuJYRLwGDKXnM7OWyvFr+qOVPb200zyppMdbwBlp/LYRy45aMlXSWmAtwIIFC1j9oy9naLqZjeZG/cP2yS7bmgO1EbE+IgYiYmDevHl1N8fMOqiq7Om780iaAbwPONzlsmbWIpWUPU3Dq9Lj64EtERFp/Mp0dug84ALghxnaZGY1qars6X3AtyQNAUcogoc032MU9ZOPA38RESd6bZOZ1UdFh6FdBgYGYnBwsO5mmPUtSdsjYmAyy7bmQK2ZtYNDxcyycqiYWVYOFTPLyqFiZlk5VMwsK4eKmWXlUDGzrBwqZpaVQ8XMsnKomFlWDhUzy8qhYmZZOVTMLCuHipll5VAxs6wcKmaWlUPFzLLqKVQkzZG0WdKedD97lHmWSvqBpBclvSDpU6VpD0h6TdKOdFvaS3vMrH699lRuB56KiAuAp9LwSEeBT0fEcGnTr0g6vTT9sxGxNN129NgeM6tZr6FSLmf6IHDdyBki4pWI2JMe/zdwEHA1MLM+1WuonBkRB9LjnwJnjjWzpEuAmcCrpdFfTLtF6yTNGmPZtZIGJQ0eOnSox2ab2VQZN1QkPSlp1yi39xRhT8XBOtb7kDQf+BZwY0S8k0bfASwGfhuYA9zWaXmXPTVrh3GLiUXElZ2mSXpd0vyIOJBC42CH+X4V+D5wZ0S8W5C91Ms5JumbwK0Tar2ZNU6vuz/lcqargH8dOUMqhfpd4KGI2DBi2vx0L4rjMbt6bI+Z1azXULkb+KikPcCVaRhJA5LuTfN8EvgIsHqUU8fflrQT2AnMBb7QY3vMrGYue2pmJ3HZUzNrDIeKmWXlUDGzrBwqZpaVQ8XMsnKomFlWDhUzy8qhYmZZOVTMLCuHipll5VAxs6wcKmaWlUPFzLJyqJhZVg4VM8vKoWJmWTlUzCwrh4qZZTXlZU/TfCdKv0+7sTT+PEnPShqS9Gj6kWwza7Eqyp4CvF0qbbq8NP4eYF1EnA+8AazpsT1mVrMpL3vaSSrLcTkwXLZjQsubWTNVVfb01FSydJuk4eA4A3gzIo6n4X3AWZ3+kMuemrXDuBUKJT0J/Nook+4sD0RESOpU72NhROyXtAjYkmr9vDWRhkbEemA9FCU6JrKsmVWnkrKnEbE/3e+V9DSwDPgX4HRJM1Jv5Wxg/yTWwcwapIqyp7MlzUqP5wKXAbtTQfetwPVjLW9m7VJF2dMLgUFJz1OEyN0RsTtNuw24RdIQxTGW+3psj5nVzGVPzewkLntqZo3hUDGzrBwqZpaVQ8XMsnKomFlWDhUzy8qhYmZZOVTMLCuHipll5VAxs6wcKmaWlUPFzLJyqJhZVg4VM8vKoWJmWTlUzCwrh4qZZeVQMbOsprzsqaQ/KJU83SHpf4dr/0h6QNJrpWlLe2mPmdVvysueRsTW4ZKnFBUJjwL/Vprls6WSqDt6bI+Z1azqsqfXA09ExNEe/66ZNVRVZU+HrQQeHjHui5JekLRuuD6QmbVXVWVPSRUMPwRsKo2+gyKMZlKUNL0NuKvD8muBtQALFiwYr9lmVpNKyp4mnwS+GxG/KD33cC/nmKRvAreO0Q7XUjZrgSkve1pyAyN2fVIQIUkUx2N29dgeM6tZFWVPkXQucA7w7yOW/7akncBOYC7whR7bY2Y1G3f3ZywRcRi4YpTxg8BNpeH/As4aZb7Le/n7ZtY8vqLWzLJyqJhZVg4VM8vKoWJmWTlUzCwrh4qZZeVQMbOsHCpmlpVDxcyycqiYWVYOFTPLyqFiZlk5VMwsK4eKmWXlUDGzrBwqZpaVQ8XMsnKomFlWDhUzy6rXWsp/JOlFSe9IGhhjvmskvSxpSNLtpfHnSXo2jX9U0sxe2mNm9eu1p7IL+EPgmU4zSDoF+BpwLbAEuEHSkjT5HmBdRJwPvAGs6bE9ZlaznkIlIl6KiJfHme0SYCgi9kbEz4FHgBWp1s/lwIY0Xze1mM2s4Xoq0dGls4CflIb3AZcCZwBvRsTx0viTyngMK5c9paho2I+Fx+YCP6u7EVOkX9etX9fr1ye7YE+1lCNirIqEWZXLnkoajIiOx3Daql/XC/p33fp5vSa7bE+1lLu0n6I64bCz07jDwOmSZqTeyvB4M2uxKk4pPwdckM70zARWAhsjIoCtwPVpvvFqMZtZC/R6SvnjkvYBvwN8X9KmNP4Dkh4HSL2Qm4FNwEvAYxHxYnqK24BbJA1RHGO5r8s/vb6XdjdYv64X9O+6eb1GUNFhMDPLw1fUmllWDhUzy6oVodLr1wGaStIcSZsl7Un3szvMd0LSjnTbWHU7uzXe6y9pVvo6xlD6esa51bdycrpYt9WSDpXep5vqaOdESLpf0sFO13yp8NW0zi9IuqirJ46Ixt+ACykuxnkaGOgwzynAq8AiYCbwPLCk7raPs15fAm5Pj28H7ukw3//U3dYu1mXc1x/4c+Dr6fFK4NG6251x3VYD/1h3Wye4Xh8BLgJ2dZj+MeAJQMCHgWe7ed5W9FSih68DTH3rerKC4usJ0P6vKXTz+pfXdwNwRfq6RtO1cdsaV0Q8AxwZY5YVwENR2EZxXdn88Z63FaHSpdG+DtDxsv+GODMiDqTHPwXO7DDfqZIGJW2T1NTg6eb1f3eeKC41eIviUoKm63bb+kTaTdgg6ZxRprfNpD5TVXz3pytN+TpAbmOtV3kgIkJSp/P7CyNiv6RFwBZJOyPi1dxttZ58D3g4Io5J+lOKHtnlNbepFo0JlZi6rwPUaqz1kvS6pPkRcSB1Kw92eI796X6vpKeBZRT7+E3Szes/PM8+STOA91F8XaPpxl23iCivx70Ux8vablKfqX7a/Rn16wA1t2k8Gym+ngAdvqYgabakWenxXOAyYHdlLexeN69/eX2vB7ZEOiLYcOOu24hjDcsprh5vu43Ap9NZoA8Db5V21zur+wh0l0epP06xP3cMeB3YlMZ/AHh8xNHqVyj+i99Zd7u7WK8zgKeAPcCTwJw0fgC4Nz3+XWAnxRmHncCauts9xvqc9PoDdwHL0+NTge8AQ8APgUV1tznjuv0d8GJ6n7YCi+tucxfr9DBwAPhF+nytAT4DfCZNF8UPrL2atr1Rz7yOvPkyfTPLqp92f8ysARwqZpaVQ8XMsnKomFlWDhUzy8qhYmZZOVTMLKv/A+YXQ8OS8hxTAAAAAElFTkSuQmCC\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAARUAAAD8CAYAAABZ0jAcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy86wFpkAAAACXBIWXMAAAsTAAALEwEAmpwYAAAQPklEQVR4nO3dbawc5XnG8euqg6XGpaVgAw6BJI2sEKcprjkyUFcJboqFrbauq1ayW5E0SmRRxVITVZWokGg+tlRpJKoE5LQWQWpBqYITCwwYorYujZxwDvJrDMFQRzG2sCGR3RQU6vbuh51jlvXuntkzz87b/n/S6szOzLP7zM7sdWZ2Z/Z2RAgAUvmZqjsAoF0IFQBJESoAkiJUACRFqABIilABkFSSULG93fYp24cGTLfte2wftX3A9squabfafj6bdkeK/gCoTqo9lfsl3Tpk+jpJy7LbFkn3SpLtBZK+lE1fLmmz7eWJ+gSgAklCJSL2SPrRkFk2SHogOvZKusT2UkmrJB2NiJci4k1JD2XzAmiod5T0PFdJ+mHX/ePZuH7jb+j3ALa3qLOXo0WLFl1/7bXXjqenyO3Y2TfG8rjv/fmfHcvjIr+ZmZlXI2LJfNqWFSruMy6GjL9wZMQ2SdskaWpqKqanp9P1Dn398e6DQ6f/ckn96HX/2g9X9MyTw/YP5tu2rFA5LunqrvvvlnRC0sIB41GyuQKkTvr1laCpj7JCZaekrbYfUufw5kxEnLR9WtIy2++T9LKkTZL+sKQ+TbQmhUgevctDyFQnSajYflDSzZIW2z4u6S8lXSRJEXGfpF2S1ks6Kul1SZ/Mpp2zvVXSE5IWSNoeEYdT9AlvaVuA5MHeTHXcxJ8+4DOVuU1ikIyCgBnO9kxETM2nLWfUAkiqrM9UUBL2UPKZfZ3YY0mPUGkBgmT+ul87AiYNQqWhCJL0CJg0CJUGIUjKQ8DMH6HSAIRJtfj8ZTSESo0RJvVCuOTDV8o1RaDUF+tmOPZUaoSNtTnYaxmMUKkBwqS5+ED3Qhz+VIxAaQ/WZQd7KhVhA2wnDovYUwGQGHsqJWMPZTJM8h4LeyolIlAmzySuc0KlJJO4caFj0tY9hz9jNmkbFPqbpMMhQmVMCBP0Mwnhkqrs6dDSpbb/3Pa+7HbI9v/avjSbdsz2wWxaK34jkkDBXNq8jRQOlTylSyPibyJiRUSskPQXkv4tIrorGq7Jps/rNzHrpM0bC9Jq67aS4vDnfOlSScrKcGyQ9L0B82+W9GCC562Vtm4gGK82Hg6lOPwZVNL0ArbfqU4h9693jQ5Ju23PZKVNG4dAQVFt2oZShEru0qWSflvSf/Qc+qyOiJXqHD59xvZH+j6JvcX2tO3p06dPF+txQm3aGFCttmxLKUJlUEnTfjap59AnIk5kf09J2qHO4dQFImJbRExFxNSSJfOqG51cWzYC1EcbtqkUofKMstKltheqExw7e2ey/QuSPirpm13jFtm+eHZY0lpJhxL0CUBFCodKRJyTNFu69Iikr0XEYdu32769a9aNknZHxH93jbtC0tO290v6rqRHI+Lxon0qQxv+o6Cemr5tUfZ0Hpq+0tEMVX4jRNnTEhEoKEtTtzVO08+pqSsYzdbE81jYU8mBQEHVmrQNEipzaNLKRLs1ZVskVIZoykrE5GjCNkmoAEiKUBmgCf8RMJnqvm0SKn3UfaUBdd5GCZUedV5ZQLe6bquECoCkCJUudU1+YJA6brOECoCkCJVMHRMfyKNu2y6hovqtFGBUddqGJz5U6rQygCLqsi1PfKgASGtif/qgLqkOpFSHn0pgTwVAUhMZKuyloO2q3MbLqqV8s+0zXfWU78rbNjUCBZOiqm298GcqXbWUb1GnBtAztndGRG/Z03+PiN+aZ1sADZFiT+V8LeWIeFPSbC3lcbcdGXspmDRVbPNl1lK+yfZ+24/Z/tCIbWtb9hTA25VVS/lZSe+JiOsk/Z2kb4zQtjOyhmVPAVyolFrKEXE2In6SDe+SdJHtxXnapsKhDyZV2dt+KbWUbV9p29nwqux5X8vTNgUCBZOuzPdA4W9/IuKc7dlaygskbZ+tpZxNv0/S70v6E9vnJL0haVN06q32bVu0TwCq0/payuylAG/Je/o+tZQB1EarQ4W9FODtynhPtDpUAJSPUAGQVGtDhUMfoL9xvzdaGyoAqkGoAEiqlaHCoQ8w3DjfI60MFQDVIVQAJNW6UOHQB8hnXO+V1oUKgGq1pu4PeyjA6MZRJ4g9FQBJESoAkmpFqHDoAxST8j3UilABUB+ECoCkyip7+ke2D2S3b9u+rmvaMdsHs3Ko+X4jEkBtlVX29D8lfTQifmx7naRtkm7omr4mIl4t2hcA1Sul7GlEfDsifpzd3atOfZ8k+JAWSCPVe6nMsqezPiXpsa77IWm37RnbWwY1ouwp0AwpzqjNXbrU9hp1QuXXu0avjogTti+X9KTt5yJizwUPGLFNncMmTU1NNa+uCDAhSil7Kkm2f0XS30vaEBGvzY6PiBPZ31OSdqhzOAWgocoqe3qNpIcl3RYR3+8av8j2xbPDktZKOpT3ifk8BUgrxXuqrLKnd0m6TNKXs5LK57LqZ1dI2pGNe4ekf4qIx4v2CUB1klylHBG7JO3qGXdf1/CnJX26T7uXJF3XOx5Ac3FGLYCkCBUASTU2VPiQFhiPou+txoYKgHoiVAAkRagASIpQAZAUoQIgKUIFQFKECoCkGhkqx86+UXUXgFa7bNny6+fbtpGhAqC+CBUASREqAJIiVAAkRagASIpQAZAUoQIgqbLKntr2Pdn0A7ZX5m0LoFkKh0pX2dN1kpZL2mx7ec9s6yQty25bJN07QlsADVJK2dPs/gPRsVfSJbaX5mwLoEFS/Jp+v7KnN+SY56qcbSV1yp6qs5eja665Rvev/XCxXgMYyC98b2a+bVPsqeQpezpontwlUyNiW0RMRcTUkiVLRuwigLKk2FPJU/Z00DwLc7QF0CCllD3N7n88+xboRklnIuJkzrYAGqSssqe7JK2XdFTS65I+Oaxt0T4BqI4j+n6EUWtTU1MxPT1ddTeA1rI9k9U7Hxln1AJIilABkBShAiApQgVAUoQKgKQIFQBJESoAkiJUACRFqABIilABkBShAiApQgVAUoQKgKQIFQBJESoAkiJUACRFqABIilABkFShULF9qe0nbb+Q/f3FPvNcbftfbB+xfdj2n3ZN+7ztl23vy27ri/QHQPWK7qncIelbEbFM0rey+73OSfqziPigpBslfaantOkXI2JFdttVsD8AKlY0VDZI+mo2/FVJv9s7Q0ScjIhns+H/knREncqEAFqoaKhckdXvUfb38mEz236vpF+V9J2u0VttH7C9vd/hU1fbLbanbU+fPn26YLcBjMucoWL7KduH+txGKqRu++ckfV3SZyPibDb6Xknvl7RC0klJXxjUnrKnQDPMWUwsIn5z0DTbr9heGhEnbS+VdGrAfBepEyj/GBEPdz32K13zfEXSI6N0HkD9FD382SnpE9nwJyR9s3cG25b0D5KORMTf9kxb2nV3o6RDBfsDoGJFQ+WvJN1i+wVJt2T3Zftdtme/yVkt6TZJv9Hnq+O7bR+0fUDSGkmfK9gfABUrVEs5Il6T9LE+40+oUztZEfG0JA9of1uR5wdQP5xRCyApQgVAUoQKgKQIFQBJESoAkiJUACRFqABIilABkBShAiApQgVAUoQKgKQIFQBJESoAkiJUACRFqABIilABkBShAiApQgVAUmMve5rNdyz7Ldp9tqdHbQ+gOcooezprTVbadGqe7QE0wNjLno65PYCaKavsaUjabXvG9pZ5tKfsKdAQc5bosP2UpCv7TLpzhOdZHREnbF8u6Unbz0XEnhHaKyK2SdomSVNTUzFKWwDlKaXsaVYHSBFxyvYOSask7ZGUqz2A5iij7Oki2xfPDktaq7fKm87ZHkCzlFH29ApJT9veL+m7kh6NiMeHtQfQXGWUPX1J0nWjtAfQXJxRCyApQgVAUoQKgKQIFQBJESoAkiJUACRFqABIilABkBShAiApQgVAUoQKgKQIFQBJESoAkiJUACRFqABIilABkBShAiApQgVAUmMve2r7A1m509nbWdufzaZ93vbLXdPWF+kPgOqNvexpRDyflTtdIel6Sa9L2tE1yxdnp0fErt72AJql7LKnH5P0YkT8oODzAqipssqeztok6cGecVttH7C9vd/hE4BmmTNUbD9l+1Cf24ZRnsj2Qkm/I+mfu0bfK+n9klZIOinpC0PaU0sZaIBSyp5m1kl6NiJe6Xrs88O2vyLpkSH9oJYy0ABjL3vaZbN6Dn2yIJq1UW+VQwXQUGWUPZXtd2bTH+5pf7ftg7YPSFoj6XMF+wOgYmMve5rdf13SZX3mu63I8wOoH86oBZAUoQIgKUIFQFKECoCkCBUASREqAJIiVAAkRagASIpQAZAUoQIgKUIFQFKECoCkCBUASREqAJIiVAAkRagASIpQAZAUoQIgKUIFQFJFayn/ge3Dtv/P9tSQ+W61/bzto7bv6Bo/Zy1mAM1SdE/lkKTfk7Rn0Ay2F0j6kjp1f5ZL2mx7eTZ5zlrMAJqlUKhExJGIeH6O2VZJOhoRL0XEm5IeUqcGszR6LWYANVeoREdOV0n6Ydf945JuyIbfVovZ9sBazLa3SNqS3f2p7TYWHlss6dWqOzEmbV22ti7XB+bbcM5Qsf2UpCv7TLozIoZVJDz/EH3GjVy2tLvsqe3piBj4GU5TtXW5pPYuW5uXa75tC9VSzum4pKu77r9b0olseJRazAAaoIyvlJ+RtMz2+2wvlLRJnRrM0mi1mAE0QNGvlDfaPi7pJkmP2n4iG3++lnJEnJO0VdITko5I+lpEHM4eom8t5hy2Fel3jbV1uaT2LhvL1cMRI3+8AQADcUYtgKQIFQBJNSJUil4OUFd5L1Owfcz2Qdv7inzVN25zvf7uuCebfsD2yir6OR85lu1m22eydbTP9l1V9HNUtrfbPjXovK95rbOIqP1N0gfVORnnXyVNDZhngaQXJf2SpIWS9ktaXnXf51iuuyXdkQ3fIemvB8x3TNLiqvs7x7LM+fpLWi/pMXXOXbpR0neq7nfCZbtZ0iNV93Uey/YRSSslHRowfeR11og9lSh+OUBdtekyhTyv/wZJD0THXkmXZOcn1V0Tt61cImKPpB8NmWXkddaIUMmp3+UAV1XUl7zedpmCpEGXKYSk3bZnsssV6ijP69/EdSTl7/dNtvfbfsz2h8rp2tiNvM7KuPYnl7pcDpDasOUa4WFWR8SJ7NqoJ20/l/2HqZM8r38t11EOefr9rKT3RMRPbK+X9A1Jy8bdsRKMvM5qEyox3ssBKjNsuWznukwhIk5kf0/Z3qHO7njdQiXP61/LdZTDnP2OiLNdw7tsf9n24oho+sWGI6+zNh3+DLscoK7mvEzB9iLbF88OS1qrzu/Y1E2e13+npI9n3yjcKOnM7OFfzc25bLavtO1seJU6763XSu9peqOvs6o/fc75CfVGdRLzp5JekfRENv5dknb1fFL9fXU+qb+z6n7nWK7L1Plxqheyv5f2Lpc63zjsz26H67xc/V5/SbdLuj0btjo/2PWipIMa8E1eHW85lm1rtn72S9or6deq7nPO5XpQ0klJ/5O9xz5VdJ1xmj6ApNp0+AOgBggVAEkRKgCSIlQAJEWoAEiKUAGQFKECIKn/B8T3ox0mtTgDAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] @@ -741,7 +741,7 @@ { "data": { "text/plain": [ - "" + "" ] }, "execution_count": 26, @@ -750,7 +750,7 @@ }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAARUAAAD8CAYAAABZ0jAcAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAARIElEQVR4nO3dXYxc5X3H8e+vUBuJqsHGFnUAGyxQjatWBm+BFiltgfCSC0MamhipwkRGbtrQSkWJAHGBRBIV0gtXUdMGCwghioDEVVRHAbkGQ7mJCWvVwcYI/ELT4DrYwYBUmTrB/HtxnqWH9c7u7J5nzsvs7yONdua8zDxn5pzfPuecOfNXRGBmlsuvNd0AMxsuDhUzy8qhYmZZOVTMLCuHipll5VAxs6yyhIqkhyQdkrSrx3hJ+pqkvZJelHRRadwaSXvSbU2O9phZc3L1VB4Grplk/LXA+em2DvhnAEnzgbuBS4CLgbslzcvUJjNrQJZQiYjngCOTTHId8EgUtgGnSVoEXA1siYgjEfEWsIXJw8nMWu7kml7nTOBnpcevp2G9hp9A0jqKXg6nnnrqymXLlg2mpda/7QN63pUDel7r2/bt238REQtnMm9doVJZRGwANgCMjIzE6Ohowy2aBdTQ604VVr6yZOAk/XSm89Z19ucAcHbp8VlpWK/hVjdNcGurLrV1FqorVDYBN6WzQJcC70TEQWAzcJWkeekA7VVpmA3asG2Uw7Y8HZZl90fSo8AfAwskvU5xRufXASLiG8ATwCeAvcBR4LNp3BFJXwJeSE91T0RMdsDXZmI2bmQTLbN3m2qRJVQi4sYpxgfw+R7jHgIeytEOK5mNQTKV8nvigBkYf6PWzLLqzNkf65N7KP0Ze5/cY8nOoTIMHCQz512i7BwqXeUgyc8Bk4VDpUscJPVxwMyYQ6ULHCbN8vGXaXGotJnDpF0cLn3xKeW2cqC0lz+bSbmn0iZeWbvDvZaeHCpt4DDpLh/QPYF3f5rmQBke/iwB91Sa4xVwOHm3yD0VM8vLPZW6uYcyO8ziHot7KnVyoMw+s/Azd6jUZRauXJbMss/euz+DNstWKOthFu0OOVQGxWFiE5kF4ZKr7Ok1kl5JZU3vmGD8ekk70u1VSW+Xxh0vjduUoz2Nc6DYVIZ4HancU5F0EvB14OMUxcBekLQpInaPTRMRf1ua/q+BC0tP8W5ErKjajtYY4pXFMhND2WPJ0VO5GNgbEfsj4pfAYxRlTnu5EXg0w+u2i0tD2EwM4XqTI1SmU7p0CXAusLU0+BRJo5K2Sbo+Q3vqN2QrhTVgiNahug/UrgY2RsTx0rAlEXFA0lJgq6SdEbFv/IzlWsqLFy+up7X9GKKVwRo2JLtDOXoq0ylduppxuz4RcSD93Q88y4ePt5Sn2xARIxExsnDhjOpG5+dAsdyGYJ3KESovAOdLOlfSHIrgOOEsjqRlwDzgR6Vh8yTNTfcXAJcBu8fPa2bdUXn3JyLek3QrRQ3kk4CHIuIlSfcAoxExFjCrgcdStcIxFwD3S3qfIuDuLZ81arUh+I9iLdXx3SB9eBvvhpGRkRgdHW2uAQ4Uq0ODm6ak7RExMpN5fe3PdDlQrC4dXdf8Nf1+dfQDto7r4Nf63VPphwPFmtahddChMpUOfZg25DqyLjpUJtORD9FmkQ6skw4VM8vKodJLB/4j2CzV8nXToTKRln9oZm1eRx0q47X4wzL7kJauqw4VM8vKoVLW0uQ366mF66xDxcyycqiMaWHim/WlZeuuQwVa96GYTVuL1mGHSos+DLNKWrIuO1TMLKvZ+9MHLUl1s6xa8FMJ7qmYWVazM1TcS7Fh1+A6Xlct5ZslHS7VTL6lNG6NpD3ptiZHeyZv7MBfwawdGlrXa6mlnDweEbeOm3c+cDcwQrEXuD3N+1bVdplZM5qopVx2NbAlIo6kINkCXJOhTRNzL8VmmwbW+TprKX9K0ouSNkoaq2g4nTrM61LN5dHDhw9naLaZDUJdB2p/AJwTEb9H0Rv51nSfoJVlT83sBLXUUo6INyPiWHr4ALCy33mz8a6PzVY1r/u11FKWtKj0cBXwcrq/Gbgq1VSeB1yVhuXlQLHZrsZtoK5ayn8jaRXwHnAEuDnNe0TSlyiCCeCeiDhStU1m1pzhr6XsXorZ/+tzc3ctZTNrjeEOFfdSzD6shm1iuEPFzGrnUDGzrIY3VLzrYzaxAW8bwxsqZtYIh4qZZTWcoeJdH7PJDXAbGc5QMbPGOFTMLKvhCxXv+pj1Z0DbyvCFipk1anjq/riHYjZ9A6gT5J6KmWXlUDGzrIYjVLzrY1ZNxm1oOELFzFrDoWJmWdVV9vQ2SbtT3Z+nJS0pjTteKoe6afy8ZtYtdZU9/Q9gJCKOSvpL4KvAZ9K4dyNiRdV2mFk71FL2NCKeiYij6eE2ivo+efggrVkembalOsuejlkLPFl6fEoqZ7pN0vW9ZnLZU7NuqPUbtZL+HBgB/qg0eElEHJC0FNgqaWdE7Bs/b0RsADZAUaKjlgab2bTVUvYUQNKVwF3AqlIJVCLiQPq7H3gWuDBDm8ysIXWVPb0QuJ8iUA6Vhs+TNDfdXwBcBpQP8E7Ox1PM8sqwTdVV9vTvgd8AvicJ4L8iYhVwAXC/pPcpAu7ecWeNzKxjul321D0Vs/zCZU/NrEUcKmaWVXdDxbs+ZoNRcdvqbqiYWSs5VMwsK4eKmWXlUDGzrBwqZpaVQ8XMsnKomFlW3QyV7U03wGy4rWTlypnO281QMbPWcqiYWVYOFTPLyqFiZlk5VMwsK4eKmWXlUDGzrOoqezpX0uNp/POSzimNuzMNf0XS1TnaY2bNqRwqpbKn1wLLgRslLR832VrgrYg4D1gP3JfmXU7x6/u/A1wD/FN6PjPrqFrKnqbH30r3NwJXqPhZ/euAxyLiWES8BuxNz2dmHZWjQuFEZU8v6TVNKunxDnB6Gr5t3LwTlkyVtA5YB7B48WL4aYaWm9mEtmv7jC+G6cyB2ojYEBEjETGycOHCpptjZj3UVfb0g2kknQx8BHizz3nNrENqKXuaHq9J928AtkZRxWwTsDqdHToXOB/4cYY2mVlD6ip7+iDwbUl7gSMUwUOa7rsU9ZPfAz4fEcertsnMmtPtsqdmNhAue2pmreFQMbOsHCpmlpVDxcyycqiYWVYOFTPLyqFiZlk5VMwsK4eKmWXlUDGzrBwqZpaVQ8XMsnKomFlWDhUzy8qhYmZZOVTMLCuHipll5VAxs6wqhYqk+ZK2SNqT/s6bYJoVkn4k6SVJL0r6TGncw5Jek7Qj3VZUaY+ZNa9qT+UO4OmIOB94Oj0e7yhwU0SMlTb9B0mnlcZ/MSJWpNuOiu0xs4ZVDZVyOdNvAdePnyAiXo2IPen+fwOHAFcDMxtSVUPljIg4mO7/HDhjsoklXQzMAfaVBn8l7RatlzR3knnXSRqVNHr48OGKzTazQZkyVCQ9JWnXBLcPFWFPxcF61vuQtAj4NvDZiHg/Db4TWAb8PjAfuL3X/C57atYNUxYTi4gre42T9IakRRFxMIXGoR7T/SbwQ+CuiPigIHupl3NM0jeBL0yr9WbWOlV3f8rlTNcA/zp+glQK9fvAIxGxcdy4RemvKI7H7KrYHjNrWNVQuRf4uKQ9wJXpMZJGJD2Qpvk08DHg5glOHX9H0k5gJ7AA+HLF9phZw1z21MxO4LKnZtYaDhUzy8qhYmZZOVTMLCuHipll5VAxs6wcKmaWlUPFzLJyqJhZVg4VM8vKoWJmWTlUzCwrh4qZZeVQMbOsHCpmlpVDxcyycqiYWVYOFTPLauBlT9N0x0u/T7upNPxcSc9L2ivp8fQj2WbWYXWUPQV4t1TadFVp+H3A+og4D3gLWFuxPWbWsIGXPe0lleW4HBgr2zGt+c2sneoqe3pKKlm6TdJYcJwOvB0R76XHrwNn9nohlz0164YpKxRKegr4rQlG3VV+EBEhqVe9jyURcUDSUmBrqvXzznQaGhEbgA1QlOiYzrxmVp9ayp5GxIH0d7+kZ4ELgX8BTpN0cuqtnAUcmMEymFmL1FH2dJ6kuen+AuAyYHcq6P4McMNk85tZt9RR9vQCYFTSTyhC5N6I2J3G3Q7cJmkvxTGWByu2x8wa5rKnZnYClz01s9ZwqJhZVg4VM8vKoWJmWTlUzCwrh4qZZeVQMbOsHCpmlpVDxcyycqiYWVYOFTPLyqFiZlk5VMwsK4eKmWXlUDGzrBwqZpaVQ8XMsnKomFlWAy97KulPSiVPd0j637HaP5IelvRaadyKKu0xs+YNvOxpRDwzVvKUoiLhUeDfSpN8sVQSdUfF9phZw+oue3oD8GREHK34umbWUnWVPR2zGnh03LCvSHpR0vqx+kBm1l11lT0lVTD8XWBzafCdFGE0h6Kk6e3APT3mXwesA1i8ePFUzTazhtRS9jT5NPD9iPhV6bnHejnHJH0T+MIk7XAtZbMOGHjZ05IbGbfrk4IISaI4HrOrYnvMrGF1lD1F0jnA2cC/j5v/O5J2AjuBBcCXK7bHzBo25e7PZCLiTeCKCYaPAreUHv8ncOYE011e5fXNrH38jVozy8qhYmZZOVTMLCuHipll5VAxs6wcKmaWlUPFzLJyqJhZVg4VM8vKoWJmWTlUzCwrh4qZZeVQMbOsHCpmlpVDxcyycqiYWVYOFTPLyqFiZlk5VMwsq6q1lP9M0kuS3pc0Msl010h6RdJeSXeUhp8r6fk0/HFJc6q0x8yaV7Wnsgv4U+C5XhNIOgn4OnAtsBy4UdLyNPo+YH1EnAe8Bayt2B4za1ilUImIlyPilSkmuxjYGxH7I+KXwGPAdanWz+XAxjRdP7WYzazlKpXo6NOZwM9Kj18HLgFOB96OiPdKw08o4zGmXPaUoqLhMBYeWwD8oulGDMiwLtuwLtdvz3TGSrWUI2KyioRZlcueShqNiJ7HcLpqWJcLhnfZhnm5ZjpvpVrKfTpAUZ1wzFlp2JvAaZJOTr2VseFm1mF1nFJ+ATg/nemZA6wGNkVEAM8AN6TppqrFbGYdUPWU8iclvQ78AfBDSZvT8I9KegIg9UJuBTYDLwPfjYiX0lPcDtwmaS/FMZYH+3zpDVXa3WLDulwwvMvm5RpHRYfBzCwPf6PWzLJyqJhZVp0IlaqXA7SVpPmStkjak/7O6zHdcUk70m1T3e3s11Tvv6S56XKMvenyjHPqb+XM9LFsN0s6XPqcbmmindMh6SFJh3p950uFr6VlflHSRX09cUS0/gZcQPFlnGeBkR7TnATsA5YCc4CfAMubbvsUy/VV4I50/w7gvh7T/U/Tbe1jWaZ8/4G/Ar6R7q8GHm+63RmX7WbgH5tu6zSX62PARcCuHuM/ATwJCLgUeL6f5+1ETyUqXA4w+NZVch3F5QnQ/csU+nn/y8u7EbgiXa7Rdl1ct6YUEc8BRyaZ5DrgkShso/he2aKpnrcTodKniS4H6Pm1/5Y4IyIOpvs/B87oMd0pkkYlbZPU1uDp5/3/YJoovmrwDsVXCdqu33XrU2k3YaOksycY3zUz2qbquPanL225HCC3yZar/CAiQlKv8/tLIuKApKXAVkk7I2Jf7rZaJT8AHo2IY5L+gqJHdnnDbWpEa0IlBnc5QKMmWy5Jb0haFBEHU7fyUI/nOJD+7pf0LHAhxT5+m/Tz/o9N87qkk4GPUFyu0XZTLltElJfjAYrjZV03o21qmHZ/JrwcoOE2TWUTxeUJ0OMyBUnzJM1N9xcAlwG7a2th//p5/8vLewOwNdIRwZabctnGHWtYRfHt8a7bBNyUzgJdCrxT2l3vrekj0H0epf4kxf7cMeANYHMa/lHgiXFHq1+l+C9+V9Pt7mO5TgeeBvYATwHz0/AR4IF0/w+BnRRnHHYCa5tu9yTLc8L7D9wDrEr3TwG+B+wFfgwsbbrNGZft74CX0uf0DLCs6Tb3sUyPAgeBX6Xtay3wOeBzabwofmBtX1r3JjzzOv7mr+mbWVbDtPtjZi3gUDGzrBwqZpaVQ8XMsnKomFlWDhUzy8qhYmZZ/R/X9q+Wg9JdSgAAAABJRU5ErkJggg==\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAARUAAAD8CAYAAABZ0jAcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy86wFpkAAAACXBIWXMAAAsTAAALEwEAmpwYAAAP2ElEQVR4nO3df6hk9X3G8eepcf/IVmrVVTfGTdOwxGygWveiay2N21TRJam1tOBSNITAYslCI6UgCDZ/toY0YEkMm1YaoTG0xE3ErD9DyzYNpt4r7g+7GlerZLOLu5qwNrUktf30jzlXx7szd87c85nz675fMNyZc8535nvmnHnu98zMmY8jQgCQ5Rea7gCAfiFUAKQiVACkIlQApCJUAKQiVACkSgkV2/fYPm774Jj5tn2X7cO299u+dGjetbafK+bdltEfAM3JGqn8naRrl5l/naSNxWWHpLslyfZpkr5YzN8kabvtTUl9AtCAlFCJiL2SfrzMItdLujcGnpB0pu31ki6TdDgiXoyIn0v6erEsgI56V02Pc4GkHw7dPlJMGzX98lF3YHuHBqMcrV27dvNFF100m56ivIUZ3e/mGd0vSltYWHg1ItatpG1doeIR02KZ6adOjNglaZckzc3Nxfz8fF7vMNqorVOHSWHFmSUzZ/vllbatK1SOSLpw6PZ7JR2VtGbMdNStqQBZidL/itCEuj5SfkDSzcWnQFsknYyIY5KelLTR9vttr5F0Y7EsZs1LLl3Xt/XpsJSRiu37JF0l6RzbRyT9uaTTJSkivixpj6Rtkg5LekPSJ4t5b9reKekRSadJuicinsnoE4asxhcZo5nGpIRKRGyfMD8kfXrMvD0ahA4yrcYgmWT4OSFgZoZv1AJIVdcbtagLI5RyFp8nRizpCJU+IEhWjkOidIRKVxEk+QiYFIRKlxAk9SFgVoxQ6QLCpFm8/zIVQqXNCJN2IVxK4SPltiJQ2ottsyxGKm3CztodjFrGIlTagDDpLt7QPQWHP00jUPqDbSmJkUpz2AH7icMiRioAcjFSqRsjlNVhFY9YGKnUiUBZfVbhNidU6rIKdy4UVtm25/Bn1lbZDoUxVtHhEKEyK4QJRlkF4ZJV9nTZ0qW2/8z208XloO3/tX1WMe8l2weKef2ou0GgYJIe7yOVRypDpUuv1qAUx5O2H4iIf19cJiI+J+lzxfIfl3RrRAxXNNwaEa9W7Usr9HhnQTKrlyOWjJHKtKVLt0u6L+Fx24XSEFiJHu43GaEyrqTpKWy/W4NC7t8YmhySHrW9UJQ27Z6e7RRoQI/2oYw3aqepsPJxSf+65NDnyog4avtcSY/ZfrYo+P7OBxmqpbxhw4aqfc7To50BDevJ4VDGSGVcSdNRbtSSQ5+IOFr8PS5ptwaHU6eIiF0RMRcRc+vWrahudD4CBdl6sE9lhEqp0qW2f0nSRyR9a2jaWttnLF6XdI2kgwl9AtCQyoc/40qX2r6lmP/lYtEbJD0aEf811Pw8SbttL/blaxHxcNU+1aIH/1HQUh0/DPKgImm3zM3Nxfx8g19pIVBQhwZfmrYXImJuJW0592daBArq0tF9ja/pl9XRDYyO6+DX+hmplEGgoGkd2gcJlUk6tDHRcx3ZFwmV5XRkI2IV6cA+SagASEWojNOB/whYpVq+bxIqo7R8owFt3kcJlaVavLGAd2jpvkqoAEhFqAxrafIDY7VwnyVUAKQiVBa1MPGBUlq27xIqUus2CjC1Fu3DhEqLNgZQSUv2ZUIFQKrV+9MHLUl1IFULfiqBkQqAVKszVBiloO8a3MfrqqV8le2TQ/WU7yjbNh2BgtWioX29llrKhX+JiI+tsC2AjmiilnJW2+kxSsFq08A+X2ct5Sts77P9kO0PT9lWtnfYnrc9f+LEiYRuA5iFjFApU0v5KUnvi4iLJf21pG9O0XYwsY1lTwGcopZayhHxekT8tLi+R9Lpts8p0zYNhz5YrWre92uppWz7fBe1TW1fVjzua2XapiBQsNrV+Bqoq5byH0j6Y9tvSvpvSTfGoN7qyLZV+wSgOf2vpcwoBXhbyZc7tZQBtEa/Q4VRCvBONbwm+h0qAGpHqABI1d9Q4dAHGG3Gr43+hgqARhAqAFL1M1Q49AGWN8PXSD9DBUBjCBUAqfoXKhz6AOXM6LXSv1AB0Kj+1P1hhAJMbwZ1ghipAEhFqABI1Y9Q4dAHqCbxNdSPUAHQGoQKgFR1lT39I9v7i8v3bF88NO8l2weKcqglfyMSQFvVVfb0PyR9JCJ+Yvs6SbskXT40f2tEvFq1LwCaV0vZ04j4XkT8pLj5hAb1fXLwJi2QI+m1VGfZ00WfkvTQ0O2Q9KjtBds7xjWi7CnQDRnfqC1dutT2Vg1C5TeHJl8ZEUdtnyvpMdvPRsTeU+4wYpcGh02am5vrXl0RYJWopeypJNn+NUl/I+n6iHhtcXpEHC3+Hpe0W4PDKQAdVVfZ0w2S7pd0U0T8YGj6WttnLF6XdI2kg6UfmfdTgFwJr6m6yp7eIelsSV8qSiq/WVQ/O0/S7mLauyR9LSIertonAM3pdtlTRipAvqDsKYAWIVQApOpuqHDoA8xGxddWd0MFQCsRKgBSESoAUhEqAFIRKgBSESoAUhEqAFJ1M1QWmu4A0G+btXnzStt2M1QAtBahAiAVoQIgFaECIBWhAiAVoQIgFaECIFVdZU9t+65i/n7bl5ZtC6BbKofKUNnT6yRtkrTd9qYli10naWNx2SHp7inaAuiQWsqeFrfvjYEnJJ1pe33JtgA6JKNC4aiyp5eXWOaCkm0lDcqeajDK0YYNG6SXq3UawHgLXljxyTAZI5UyZU/HLVO6ZGpE7IqIuYiYW7du3ZRdBFCXjJFKmbKn45ZZU6ItgA6ppexpcfvm4lOgLZJORsSxkm0BdEhdZU/3SNom6bCkNyR9crm2VfsEoDndLnsKYCYoewqgNQgVAKkIFQCpCBUAqQgVAKkIFQCpCBUAqQgVAKkIFQCpCBUAqQgVAKkIFQCpCBUAqQgVAKkIFQCpCBUAqQgVAKkIFQCpKoWK7bNsP2b7+eLvL49Y5kLb/2T7kO1nbP/J0LzP2v6R7aeLy7Yq/QHQvKojldskfSciNkr6TnF7qTcl/WlEfEjSFkmfXlLa9AsRcUlx2VOxPwAaVjVUrpf01eL6VyX93tIFIuJYRDxVXP9PSYc0qEwIoIeqhsp5Rf0eFX/PXW5h278i6dclfX9o8k7b+23fM+rwaajtDtvztudPnDhRsdsAZmViqNh+3PbBEZepCqnb/kVJ35D0mYh4vZh8t6QPSLpE0jFJnx/XnrKnQDdMLCYWEb8zbp7tV2yvj4hjttdLOj5mudM1CJS/j4j7h+77laFlviLpwWk6D6B9qh7+PCDpE8X1T0j61tIFbFvS30o6FBF/tWTe+qGbN0g6WLE/ABpWNVT+QtLVtp+XdHVxW7bfY3vxk5wrJd0k6bdHfHR8p+0DtvdL2irp1or9AdCwSrWUI+I1SR8dMf2oBrWTFRHfleQx7W+q8vgA2odv1AJIRagASEWoAEhFqABIRagASEWoAEhFqABIRagASEWoAEhFqABIRagASEWoAEhFqABIRagASEWoAEhFqABIRagASEWoAEg187KnxXIvFb9F+7Tt+WnbA+iOOsqeLtpalDadW2F7AB0w87KnM24PoGXqKnsakh61vWB7xwraU/YU6IiJJTpsPy7p/BGzbp/ica6MiKO2z5X0mO1nI2LvFO0VEbsk7ZKkubm5mKYtgPrUUva0qAOkiDhue7ekyyTtlVSqPYDuqKPs6VrbZyxel3SN3i5vOrE9gG6po+zpeZK+a3ufpH+T9O2IeHi59gC6q46ypy9Kunia9gC6i2/UAkhFqABIRagASEWoAEhFqABIRagASEWoAEhFqABIRagASEWoAEhFqABIRagASEWoAEhFqABIRagASEWoAEhFqABIRagASDXzsqe2P1iUO128vG77M8W8z9r+0dC8bVX6A6B5My97GhHPFeVOL5G0WdIbknYPLfKFxfkRsWdpewDdUnfZ049KeiEiXq74uABaqq6yp4tulHTfkmk7be+3fc+owycA3TIxVGw/bvvgiMv10zyQ7TWSflfSPw5NvlvSByRdIumYpM8v055aykAH1FL2tHCdpKci4pWh+37ruu2vSHpwmX5QSxnogJmXPR2yXUsOfYogWnSD3i6HCqCj6ih7KtvvLubfv6T9nbYP2N4vaaukWyv2B0DDZl72tLj9hqSzRyx3U5XHB9A+fKMWQCpCBUAqQgVAKkIFQCpCBUAqQgVAKkIFQCpCBUAqQgVAKkIFQCpCBUAqQgVAKkIFQCpCBUAqQgVAKkIFQCpCBUAqQgVAKkIFQKqqtZT/0PYztv/P9twyy11r+znbh23fNjR9Yi1mAN1SdaRyUNLvS9o7bgHbp0n6ogZ1fzZJ2m57UzF7Yi1mAN1SKVQi4lBEPDdhscskHY6IFyPi55K+rkENZmn6WswAWq5SiY6SLpD0w6HbRyRdXlx/Ry1m22NrMdveIWlHcfNntvtYeOwcSa823YkZ6eu69XW9PrjShhNDxfbjks4fMev2iFiuIuFbdzFi2tRlS4fLntqej4ix7+F0VV/XS+rvuvV5vVbatlIt5ZKOSLpw6PZ7JR0trk9TixlAB9TxkfKTkjbafr/tNZJu1KAGszRdLWYAHVD1I+UbbB+RdIWkb9t+pJj+Vi3liHhT0k5Jj0g6JOkfIuKZ4i5G1mIuYVeVfrdYX9dL6u+6sV5LOGLqtzcAYCy+UQsgFaECIFUnQqXq6QBtVfY0Bdsv2T5g++kqH/XN2qTn3wN3FfP32760iX6uRIl1u8r2yWIbPW37jib6OS3b99g+Pu57XyvaZhHR+oukD2nwZZx/ljQ3ZpnTJL0g6VclrZG0T9Kmpvs+Yb3ulHRbcf02SX85ZrmXJJ3TdH8nrMvE51/SNkkPafDdpS2Svt90vxPX7SpJDzbd1xWs229JulTSwTHzp95mnRipRPXTAdqqT6cplHn+r5d0bww8IenM4vtJbdfFfauUiNgr6cfLLDL1NutEqJQ06nSACxrqS1nvOE1B0rjTFELSo7YXitMV2qjM89/FbSSV7/cVtvfZfsj2h+vp2sxNvc3qOPenlLacDpBtufWa4m6ujIijxblRj9l+tvgP0yZlnv9WbqMSyvT7KUnvi4if2t4m6ZuSNs66YzWYepu1JlRitqcDNGa59bJd6jSFiDha/D1ue7cGw/G2hUqZ57+V26iEif2OiNeHru+x/SXb50RE1082nHqb9enwZ7nTAdpq4mkKttfaPmPxuqRrNPgdm7Yp8/w/IOnm4hOFLZJOLh7+tdzEdbN9vm0X1y/T4LX1Wu09zTf9Nmv63eeS71DfoEFi/kzSK5IeKaa/R9KeJe9U/0CDd+pvb7rfJdbrbA1+nOr54u9ZS9dLg08c9hWXZ9q8XqOef0m3SLqluG4NfrDrBUkHNOaTvDZeSqzbzmL77JP0hKTfaLrPJdfrPknHJP1P8Rr7VNVtxtf0AaTq0+EPgBYgVACkIlQApCJUAKQiVACkIlQApCJUAKT6fzYsyUmXyjVrAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] @@ -787,9 +787,9 @@ "metadata": {}, "outputs": [], "source": [ - "fuel_or = openmc.ZCylinder(r=0.39)\n", - "clad_ir = openmc.ZCylinder(r=0.40)\n", - "clad_or = openmc.ZCylinder(r=0.46)" + "fuel_outer_radius = openmc.ZCylinder(r=0.39)\n", + "clad_inner_radius = openmc.ZCylinder(r=0.40)\n", + "clad_outer_radius = openmc.ZCylinder(r=0.46)" ] }, { @@ -805,9 +805,9 @@ "metadata": {}, "outputs": [], "source": [ - "fuel_region = -fuel_or\n", - "gap_region = +fuel_or & -clad_ir\n", - "clad_region = +clad_ir & -clad_or" + "fuel_region = -fuel_outer_radius\n", + "gap_region = +fuel_outer_radius & -clad_inner_radius\n", + "clad_region = +clad_inner_radius & -clad_outer_radius" ] }, { @@ -826,9 +826,9 @@ "name": "stderr", "output_type": "stream", "text": [ - "/home/sam/openmc/openmc/openmc/mixin.py:71: IDWarning: Another Cell instance already exists with id=1.\n", + "/home/master/anaconda3/envs/OMC/lib/python3.8/site-packages/openmc/mixin.py:68: IDWarning: Another Cell instance already exists with id=1.\n", " warn(msg, IDWarning)\n", - "/home/sam/openmc/openmc/openmc/mixin.py:71: IDWarning: Another Cell instance already exists with id=2.\n", + "/home/master/anaconda3/envs/OMC/lib/python3.8/site-packages/openmc/mixin.py:68: IDWarning: Another Cell instance already exists with id=2.\n", " warn(msg, IDWarning)\n" ] } @@ -879,7 +879,7 @@ "metadata": {}, "outputs": [], "source": [ - "water_region = +left & -right & +bottom & -top & +clad_or\n", + "water_region = +left & -right & +bottom & -top & +clad_outer_radius\n", "\n", "moderator = openmc.Cell(4, 'moderator')\n", "moderator.fill = water\n", @@ -928,7 +928,7 @@ "metadata": {}, "outputs": [], "source": [ - "water_region = box & +clad_or" + "water_region = box & +clad_outer_radius" ] }, { @@ -965,14 +965,14 @@ } ], "source": [ - "root = openmc.Universe(cells=(fuel, gap, clad, moderator))\n", + "root_universe = openmc.Universe(cells=(fuel, gap, clad, moderator))\n", "\n", - "geom = openmc.Geometry()\n", - "geom.root_universe = root\n", + "geometry = openmc.Geometry()\n", + "geometry.root_universe = root_universe\n", "\n", "# or...\n", - "geom = openmc.Geometry(root)\n", - "geom.export_to_xml()\n", + "geometry = openmc.Geometry(root_universe)\n", + "geometry.export_to_xml()\n", "!cat geometry.xml" ] }, @@ -991,8 +991,9 @@ "metadata": {}, "outputs": [], "source": [ + "# Create a point source\n", "point = openmc.stats.Point((0, 0, 0))\n", - "src = openmc.Source(space=point)" + "source = openmc.Source(space=point)" ] }, { @@ -1008,11 +1009,16 @@ "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.source = src\n", - "settings.batches = 100\n", - "settings.inactive = 10\n", - "settings.particles = 1000" + "settings.batches = batches\n", + "settings.inactive = inactive\n", + "settings.particles = particles" ] }, { @@ -1030,11 +1036,6 @@ " 1000\r\n", " 100\r\n", " 10\r\n", - " \r\n", - " \r\n", - " 0 0 0\r\n", - " \r\n", - " \r\n", "\r\n" ] } @@ -1065,8 +1066,8 @@ "source": [ "cell_filter = openmc.CellFilter(fuel)\n", "\n", - "t = openmc.Tally(1)\n", - "t.filters = [cell_filter]" + "tally = openmc.Tally(1)\n", + "tally.filters = [cell_filter]" ] }, { @@ -1082,8 +1083,8 @@ "metadata": {}, "outputs": [], "source": [ - "t.nuclides = ['U235']\n", - "t.scores = ['total', 'fission', 'absorption', '(n,gamma)']" + "tally.nuclides = ['U235']\n", + "tally.scores = ['total', 'fission', 'absorption', '(n,gamma)']" ] }, { @@ -1117,7 +1118,7 @@ } ], "source": [ - "tallies = openmc.Tallies([t])\n", + "tallies = openmc.Tallies([tally])\n", "tallies.export_to_xml()\n", "!cat tallies.xml" ] @@ -1168,165 +1169,164 @@ "\n", " | The OpenMC Monte Carlo Code\n", " Copyright | 2011-2020 MIT and OpenMC contributors\n", - " License | http://openmc.readthedocs.io/en/latest/license.html\n", - " Version | 0.12.0-dev\n", - " Git SHA1 | dd74b2f43f1d2060486de2823ee30058b8971dbd\n", - " Date/Time | 2020-03-03 13:58:53\n", - " MPI Processes | 1\n", - " OpenMP Threads | 8\n", + " License | https://docs.openmc.org/en/latest/license.html\n", + " Version | 0.12.0\n", + " Git SHA1 | 3d90a9f857ec72eae897e054d4225180f1fa4d93\n", + " Date/Time | 2020-08-15 06:49:11\n", + " OpenMP Threads | 4\n", "\n", " Reading settings XML file...\n", " Reading cross sections XML file...\n", " Reading materials XML file...\n", " Reading geometry XML file...\n", - " Reading U235 from /home/sam/openmc/libs/nndc_hdf5/U235.h5\n", - " Reading U238 from /home/sam/openmc/libs/nndc_hdf5/U238.h5\n", - " Reading O16 from /home/sam/openmc/libs/nndc_hdf5/O16.h5\n", - " Reading Zr90 from /home/sam/openmc/libs/nndc_hdf5/Zr90.h5\n", - " Reading Zr91 from /home/sam/openmc/libs/nndc_hdf5/Zr91.h5\n", - " Reading Zr92 from /home/sam/openmc/libs/nndc_hdf5/Zr92.h5\n", - " Reading Zr94 from /home/sam/openmc/libs/nndc_hdf5/Zr94.h5\n", - " Reading Zr96 from /home/sam/openmc/libs/nndc_hdf5/Zr96.h5\n", - " Reading H1 from /home/sam/openmc/libs/nndc_hdf5/H1.h5\n", - " Reading O17 from /home/sam/openmc/libs/nndc_hdf5/O17.h5\n", - " Reading c_H_in_H2O from /home/sam/openmc/libs/nndc_hdf5/c_H_in_H2O.h5\n", - " Maximum neutron transport energy: 20000000.000000 eV for U235\n", + " Reading U235 from /home/master/data/nuclear/endfb71_hdf5/U235.h5\n", + " Reading U238 from /home/master/data/nuclear/endfb71_hdf5/U238.h5\n", + " Reading O16 from /home/master/data/nuclear/endfb71_hdf5/O16.h5\n", + " Reading Zr90 from /home/master/data/nuclear/endfb71_hdf5/Zr90.h5\n", + " Reading Zr91 from /home/master/data/nuclear/endfb71_hdf5/Zr91.h5\n", + " Reading Zr92 from /home/master/data/nuclear/endfb71_hdf5/Zr92.h5\n", + " Reading Zr94 from /home/master/data/nuclear/endfb71_hdf5/Zr94.h5\n", + " Reading Zr96 from /home/master/data/nuclear/endfb71_hdf5/Zr96.h5\n", + " Reading H1 from /home/master/data/nuclear/endfb71_hdf5/H1.h5\n", + " Reading O17 from /home/master/data/nuclear/endfb71_hdf5/O17.h5\n", + " Reading c_H_in_H2O from /home/master/data/nuclear/endfb71_hdf5/c_H_in_H2O.h5\n", " Minimum neutron data temperature: 294.000000 K\n", " Maximum neutron data temperature: 294.000000 K\n", " Reading tallies XML file...\n", " Preparing distributed cell instances...\n", " Writing summary.h5 file...\n", + " Maximum neutron transport energy: 20000000.000000 eV for U235\n", " Initializing source particles...\n", "\n", " ====================> K EIGENVALUE SIMULATION <====================\n", "\n", " Bat./Gen. k Average k\n", " ========= ======== ====================\n", - " 1/1 1.32572\n", - " 2/1 1.46138\n", - " 3/1 1.46068\n", - " 4/1 1.39592\n", - " 5/1 1.37519\n", - " 6/1 1.38777\n", - " 7/1 1.50242\n", - " 8/1 1.42042\n", - " 9/1 1.47458\n", - " 10/1 1.49148\n", - " 11/1 1.39339\n", - " 12/1 1.40637 1.39988 +/- 0.00649\n", - " 13/1 1.42972 1.40983 +/- 0.01063\n", - " 14/1 1.46319 1.42317 +/- 0.01531\n", - " 15/1 1.41538 1.42161 +/- 0.01196\n", - " 16/1 1.38163 1.41494 +/- 0.01182\n", - " 17/1 1.41257 1.41461 +/- 0.01000\n", - " 18/1 1.43455 1.41710 +/- 0.00901\n", - " 19/1 1.33136 1.40757 +/- 0.01241\n", - " 20/1 1.41560 1.40837 +/- 0.01113\n", - " 21/1 1.38911 1.40662 +/- 0.01021\n", - " 22/1 1.28621 1.39659 +/- 0.01370\n", - " 23/1 1.45693 1.40123 +/- 0.01343\n", - " 24/1 1.46839 1.40603 +/- 0.01333\n", - " 25/1 1.46738 1.41012 +/- 0.01306\n", - " 26/1 1.43977 1.41197 +/- 0.01236\n", - " 27/1 1.44066 1.41366 +/- 0.01173\n", - " 28/1 1.39358 1.41254 +/- 0.01112\n", - " 29/1 1.39142 1.41143 +/- 0.01057\n", - " 30/1 1.38525 1.41012 +/- 0.01012\n", - " 31/1 1.38025 1.40870 +/- 0.00973\n", - " 32/1 1.45348 1.41074 +/- 0.00949\n", - " 33/1 1.35893 1.40848 +/- 0.00935\n", - " 34/1 1.32332 1.40493 +/- 0.00963\n", - " 35/1 1.46285 1.40725 +/- 0.00952\n", - " 36/1 1.33760 1.40457 +/- 0.00953\n", - " 37/1 1.41117 1.40482 +/- 0.00917\n", - " 38/1 1.45574 1.40664 +/- 0.00903\n", - " 39/1 1.43472 1.40760 +/- 0.00876\n", - " 40/1 1.30110 1.40405 +/- 0.00918\n", - " 41/1 1.41765 1.40449 +/- 0.00889\n", - " 42/1 1.45300 1.40601 +/- 0.00874\n", - " 43/1 1.40491 1.40597 +/- 0.00847\n", - " 44/1 1.42053 1.40640 +/- 0.00823\n", - " 45/1 1.38805 1.40588 +/- 0.00801\n", - " 46/1 1.34293 1.40413 +/- 0.00798\n", - " 47/1 1.35441 1.40279 +/- 0.00787\n", - " 48/1 1.29370 1.39991 +/- 0.00818\n", - " 49/1 1.48467 1.40209 +/- 0.00826\n", - " 50/1 1.41759 1.40248 +/- 0.00806\n", - " 51/1 1.37151 1.40172 +/- 0.00790\n", - " 52/1 1.42403 1.40225 +/- 0.00773\n", - " 53/1 1.38826 1.40193 +/- 0.00755\n", - " 54/1 1.48944 1.40392 +/- 0.00764\n", - " 55/1 1.41452 1.40415 +/- 0.00747\n", - " 56/1 1.47337 1.40566 +/- 0.00746\n", - " 57/1 1.35700 1.40462 +/- 0.00738\n", - " 58/1 1.40305 1.40459 +/- 0.00722\n", - " 59/1 1.41608 1.40482 +/- 0.00708\n", - " 60/1 1.47254 1.40618 +/- 0.00706\n", - " 61/1 1.36847 1.40544 +/- 0.00696\n", - " 62/1 1.34103 1.40420 +/- 0.00694\n", - " 63/1 1.39510 1.40403 +/- 0.00681\n", - " 64/1 1.40228 1.40399 +/- 0.00668\n", - " 65/1 1.29401 1.40200 +/- 0.00686\n", - " 66/1 1.42693 1.40244 +/- 0.00675\n", - " 67/1 1.36447 1.40177 +/- 0.00666\n", - " 68/1 1.37498 1.40131 +/- 0.00656\n", - " 69/1 1.36958 1.40077 +/- 0.00647\n", - " 70/1 1.38585 1.40053 +/- 0.00637\n", - " 71/1 1.42133 1.40087 +/- 0.00627\n", - " 72/1 1.44900 1.40164 +/- 0.00622\n", - " 73/1 1.37696 1.40125 +/- 0.00613\n", - " 74/1 1.48851 1.40261 +/- 0.00619\n", - " 75/1 1.38933 1.40241 +/- 0.00610\n", - " 76/1 1.41780 1.40264 +/- 0.00601\n", - " 77/1 1.41054 1.40276 +/- 0.00592\n", - " 78/1 1.38194 1.40246 +/- 0.00584\n", - " 79/1 1.38446 1.40219 +/- 0.00576\n", - " 80/1 1.37504 1.40181 +/- 0.00569\n", - " 81/1 1.40550 1.40186 +/- 0.00561\n", - " 82/1 1.49785 1.40319 +/- 0.00569\n", - " 83/1 1.35613 1.40255 +/- 0.00565\n", - " 84/1 1.41786 1.40275 +/- 0.00557\n", - " 85/1 1.38444 1.40251 +/- 0.00550\n", - " 86/1 1.40459 1.40254 +/- 0.00543\n", - " 87/1 1.39923 1.40249 +/- 0.00536\n", - " 88/1 1.44540 1.40304 +/- 0.00532\n", - " 89/1 1.45962 1.40376 +/- 0.00530\n", - " 90/1 1.37057 1.40335 +/- 0.00525\n", - " 91/1 1.38115 1.40307 +/- 0.00519\n", - " 92/1 1.35758 1.40252 +/- 0.00516\n", - " 93/1 1.34508 1.40182 +/- 0.00514\n", - " 94/1 1.31471 1.40079 +/- 0.00519\n", - " 95/1 1.41434 1.40095 +/- 0.00513\n", - " 96/1 1.33895 1.40023 +/- 0.00512\n", - " 97/1 1.44716 1.40077 +/- 0.00509\n", - " 98/1 1.38455 1.40058 +/- 0.00503\n", - " 99/1 1.52127 1.40194 +/- 0.00516\n", - " 100/1 1.35488 1.40141 +/- 0.00513\n", + " 1/1 1.42066\n", + " 2/1 1.39831\n", + " 3/1 1.46207\n", + " 4/1 1.44888\n", + " 5/1 1.42595\n", + " 6/1 1.35549\n", + " 7/1 1.36717\n", + " 8/1 1.45095\n", + " 9/1 1.36061\n", + " 10/1 1.36554\n", + " 11/1 1.36973\n", + " 12/1 1.44276 1.40625 +/- 0.03652\n", + " 13/1 1.35512 1.38920 +/- 0.02711\n", + " 14/1 1.54216 1.42744 +/- 0.04277\n", + " 15/1 1.39353 1.42066 +/- 0.03382\n", + " 16/1 1.38650 1.41497 +/- 0.02820\n", + " 17/1 1.38760 1.41106 +/- 0.02415\n", + " 18/1 1.38413 1.40769 +/- 0.02118\n", + " 19/1 1.39088 1.40582 +/- 0.01877\n", + " 20/1 1.47468 1.41271 +/- 0.01815\n", + " 21/1 1.45695 1.41673 +/- 0.01690\n", + " 22/1 1.40308 1.41559 +/- 0.01547\n", + " 23/1 1.40821 1.41503 +/- 0.01424\n", + " 24/1 1.32301 1.40845 +/- 0.01473\n", + " 25/1 1.36702 1.40569 +/- 0.01399\n", + " 26/1 1.30968 1.39969 +/- 0.01440\n", + " 27/1 1.38099 1.39859 +/- 0.01357\n", + " 28/1 1.42103 1.39984 +/- 0.01285\n", + " 29/1 1.39741 1.39971 +/- 0.01216\n", + " 30/1 1.36548 1.39800 +/- 0.01166\n", + " 31/1 1.41573 1.39884 +/- 0.01112\n", + " 32/1 1.39788 1.39880 +/- 0.01061\n", + " 33/1 1.35942 1.39709 +/- 0.01028\n", + " 34/1 1.40483 1.39741 +/- 0.00985\n", + " 35/1 1.39418 1.39728 +/- 0.00944\n", + " 36/1 1.41492 1.39796 +/- 0.00910\n", + " 37/1 1.49392 1.40151 +/- 0.00945\n", + " 38/1 1.45114 1.40329 +/- 0.00928\n", + " 39/1 1.42619 1.40408 +/- 0.00899\n", + " 40/1 1.35249 1.40236 +/- 0.00885\n", + " 41/1 1.35401 1.40080 +/- 0.00870\n", + " 42/1 1.40220 1.40084 +/- 0.00842\n", + " 43/1 1.36437 1.39974 +/- 0.00824\n", + " 44/1 1.33642 1.39787 +/- 0.00821\n", + " 45/1 1.36953 1.39706 +/- 0.00801\n", + " 46/1 1.30034 1.39438 +/- 0.00824\n", + " 47/1 1.44097 1.39564 +/- 0.00811\n", + " 48/1 1.37981 1.39522 +/- 0.00790\n", + " 49/1 1.34870 1.39403 +/- 0.00779\n", + " 50/1 1.41247 1.39449 +/- 0.00761\n", + " 51/1 1.33382 1.39301 +/- 0.00756\n", + " 52/1 1.37043 1.39247 +/- 0.00740\n", + " 53/1 1.38754 1.39236 +/- 0.00723\n", + " 54/1 1.40160 1.39257 +/- 0.00707\n", + " 55/1 1.37511 1.39218 +/- 0.00692\n", + " 56/1 1.38589 1.39204 +/- 0.00677\n", + " 57/1 1.40630 1.39234 +/- 0.00663\n", + " 58/1 1.29944 1.39041 +/- 0.00677\n", + " 59/1 1.40019 1.39061 +/- 0.00663\n", + " 60/1 1.42384 1.39127 +/- 0.00653\n", + " 61/1 1.36502 1.39076 +/- 0.00643\n", + " 62/1 1.37042 1.39037 +/- 0.00631\n", + " 63/1 1.42295 1.39098 +/- 0.00622\n", + " 64/1 1.40042 1.39116 +/- 0.00611\n", + " 65/1 1.36382 1.39066 +/- 0.00602\n", + " 66/1 1.31659 1.38934 +/- 0.00606\n", + " 67/1 1.36101 1.38884 +/- 0.00597\n", + " 68/1 1.46359 1.39013 +/- 0.00601\n", + " 69/1 1.41012 1.39047 +/- 0.00591\n", + " 70/1 1.27411 1.38853 +/- 0.00613\n", + " 71/1 1.45399 1.38960 +/- 0.00612\n", + " 72/1 1.40455 1.38984 +/- 0.00603\n", + " 73/1 1.33020 1.38890 +/- 0.00601\n", + " 74/1 1.44599 1.38979 +/- 0.00598\n", + " 75/1 1.34985 1.38917 +/- 0.00592\n", + " 76/1 1.36183 1.38876 +/- 0.00584\n", + " 77/1 1.41080 1.38909 +/- 0.00576\n", + " 78/1 1.43991 1.38984 +/- 0.00573\n", + " 79/1 1.35613 1.38935 +/- 0.00566\n", + " 80/1 1.31659 1.38831 +/- 0.00568\n", + " 81/1 1.51344 1.39007 +/- 0.00587\n", + " 82/1 1.38404 1.38999 +/- 0.00579\n", + " 83/1 1.39613 1.39007 +/- 0.00571\n", + " 84/1 1.43037 1.39061 +/- 0.00566\n", + " 85/1 1.47316 1.39172 +/- 0.00569\n", + " 86/1 1.39220 1.39172 +/- 0.00561\n", + " 87/1 1.44400 1.39240 +/- 0.00558\n", + " 88/1 1.42419 1.39281 +/- 0.00552\n", + " 89/1 1.30930 1.39175 +/- 0.00556\n", + " 90/1 1.46976 1.39273 +/- 0.00557\n", + " 91/1 1.38334 1.39261 +/- 0.00550\n", + " 92/1 1.35260 1.39212 +/- 0.00546\n", + " 93/1 1.38505 1.39204 +/- 0.00539\n", + " 94/1 1.38290 1.39193 +/- 0.00533\n", + " 95/1 1.42597 1.39233 +/- 0.00528\n", + " 96/1 1.41624 1.39261 +/- 0.00523\n", + " 97/1 1.42053 1.39293 +/- 0.00518\n", + " 98/1 1.36268 1.39258 +/- 0.00513\n", + " 99/1 1.39175 1.39258 +/- 0.00507\n", + " 100/1 1.38148 1.39245 +/- 0.00502\n", " Creating state point statepoint.100.h5...\n", "\n", " =======================> TIMING STATISTICS <=======================\n", "\n", - " Total time for initialization = 9.7663e-01 seconds\n", - " Reading cross sections = 8.7037e-01 seconds\n", - " Total time in simulation = 1.8046e+00 seconds\n", - " Time in transport only = 1.7523e+00 seconds\n", - " Time in inactive batches = 2.0849e-01 seconds\n", - " Time in active batches = 1.5961e+00 seconds\n", - " Time synchronizing fission bank = 6.7422e-03 seconds\n", - " Sampling source sites = 4.9715e-03 seconds\n", - " SEND/RECV source sites = 1.1323e-03 seconds\n", - " Time accumulating tallies = 5.3810e-04 seconds\n", - " Total time for finalization = 5.1810e-04 seconds\n", - " Total time elapsed = 2.7828e+00 seconds\n", - " Calculation Rate (inactive) = 47962.8 particles/second\n", - " Calculation Rate (active) = 56386.7 particles/second\n", + " Total time for initialization = 6.9980e-01 seconds\n", + " Reading cross sections = 6.8788e-01 seconds\n", + " Total time in simulation = 1.9251e+00 seconds\n", + " Time in transport only = 1.9072e+00 seconds\n", + " Time in inactive batches = 1.5794e-01 seconds\n", + " Time in active batches = 1.7672e+00 seconds\n", + " Time synchronizing fission bank = 4.2579e-03 seconds\n", + " Sampling source sites = 3.4948e-03 seconds\n", + " SEND/RECV source sites = 5.9671e-04 seconds\n", + " Time accumulating tallies = 7.9162e-05 seconds\n", + " Total time for finalization = 5.7305e-05 seconds\n", + " Total time elapsed = 2.6291e+00 seconds\n", + " Calculation Rate (inactive) = 63313.6 particles/second\n", + " Calculation Rate (active) = 50928.0 particles/second\n", "\n", " ============================> RESULTS <============================\n", "\n", - " k-effective (Collision) = 1.39737 +/- 0.00470\n", - " k-effective (Track-length) = 1.40141 +/- 0.00513\n", - " k-effective (Absorption) = 1.39596 +/- 0.00308\n", - " Combined k-effective = 1.39719 +/- 0.00286\n", + " k-effective (Collision) = 1.39516 +/- 0.00457\n", + " k-effective (Track-length) = 1.39245 +/- 0.00502\n", + " k-effective (Absorption) = 1.40443 +/- 0.00333\n", + " Combined k-effective = 1.40145 +/- 0.00319\n", " Leakage Fraction = 0.00000 +/- 0.00000\n", "\n" ] @@ -1356,10 +1356,10 @@ "\r\n", " Cell 1\r\n", " U235\r\n", - " Total Reaction Rate 0.731003 +/- 0.00253759\r\n", - " Fission Rate 0.547587 +/- 0.00210114\r\n", - " Absorption Rate 0.657406 +/- 0.0024539\r\n", - " (n,gamma) 0.109821 +/- 0.000368054\r\n" + " Total Reaction Rate 0.726151 +/- 0.00251702\r\n", + " Fission Rate 0.543836 +/- 0.00205084\r\n", + " Absorption Rate 0.652874 +/- 0.002424\r\n", + " (n,gamma) 0.10904 +/- 0.000385793\r\n" ] } ], @@ -1382,12 +1382,12 @@ "metadata": {}, "outputs": [], "source": [ - "p = openmc.Plot()\n", - "p.filename = 'pinplot'\n", - "p.width = (pitch, pitch)\n", - "p.pixels = (200, 200)\n", - "p.color_by = 'material'\n", - "p.colors = {uo2: 'yellow', water: 'blue'}" + "plot = openmc.Plot()\n", + "plot.filename = 'pinplot'\n", + "plot.width = (pitch, pitch)\n", + "plot.pixels = (200, 200)\n", + "plot.color_by = 'material'\n", + "plot.colors = {uo2: 'yellow', water: 'blue'}" ] }, { @@ -1420,7 +1420,7 @@ } ], "source": [ - "plots = openmc.Plots([p])\n", + "plots = openmc.Plots([plot])\n", "plots.export_to_xml()\n", "!cat plots.xml" ] @@ -1467,12 +1467,11 @@ "\n", " | The OpenMC Monte Carlo Code\n", " Copyright | 2011-2020 MIT and OpenMC contributors\n", - " License | http://openmc.readthedocs.io/en/latest/license.html\n", - " Version | 0.12.0-dev\n", - " Git SHA1 | dd74b2f43f1d2060486de2823ee30058b8971dbd\n", - " Date/Time | 2020-03-03 13:58:56\n", - " MPI Processes | 1\n", - " OpenMP Threads | 8\n", + " License | https://docs.openmc.org/en/latest/license.html\n", + " Version | 0.12.0\n", + " Git SHA1 | 3d90a9f857ec72eae897e054d4225180f1fa4d93\n", + " Date/Time | 2020-08-15 06:49:14\n", + " OpenMP Threads | 4\n", "\n", " Reading settings XML file...\n", " Reading cross sections XML file...\n", @@ -1534,7 +1533,7 @@ "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAMgAAADIAgMAAADQNkYNAAAABGdBTUEAALGPC/xhBQAAACBjSFJNAAB6JgAAgIQAAPoAAACA6AAAdTAAAOpgAAA6mAAAF3CculE8AAAADFBMVEUAAP9yEhL//////wAZPRNOAAAAAWJLR0QCZgt8ZAAAAAd0SU1FB+QDAw06OBFn074AAAI7SURBVGje7ZmxkcIwEEUhcAnuhxIIEIGvAwioxiU4QASU4GougA4gsM6SfNxh766sfx7N3Ix+zJv9X7KxtLtaZWVlZWWl1Fr12sQQpXLaRRMRTKFe2kbk+Na8POpNkbbmWlMqtkwxRsJllIotMxQ56F7neWX8LmrdGtNpPWc//Z4cLsbpfp6xN85XNRA904Sd2V/stXlJ16EFcL4O7Q/SnUPOrK/9xfzSvQ44K0dFhjLCmq0nRYYyvDPrq2rfka4RnVlfoyJ9GdGZ9XUdI89aWGYbpTITNUKYgvLlnXFhSsqXd7bjEcKXc8YhvYEjhdzY/Db9lUKebP4+/b6lkK7m8nNRhDAlE8WFoREuig/DpKej+DAbGmGiuDAUUsgItWR8ejZ/yaZ3+SmET+/zL4JIC8Ys2VpI7/JPkSKEbCnkyiNPCilDyI5A+AVzSzZFVAhRBCKssVvlBRB5W8iNwZCThDymiLyT5F5iSCsh3RSRN5/cfgwxoiikkpFmgqgwopZAPmTkcxnkKCO3BZDQs088/RhykpFHRjLyP5AE70uadz/Z/1iK/2TgYwF9klJ8K4GPeJrTBXDsSXMeS3O2RA69wNEaOMAD1wTgMgJceYCLVZobH3IVBS68wLUauLwDLQKgEYG0O4CmCtC6ARpEQBsKaXYBLTWgcQe0B4EmJNLqBBqqQNsWaA4jLWig0Q2005GmPTAaAAYQyJgDGKYgIxtgMISMn4AhFzJKAwZ2yFgQGT4iI86srKysrD/rC4LWcCSWwIp+AAAAJXRFWHRkYXRlOmNyZWF0ZQAyMDIwLTAzLTAzVDEzOjU4OjU2KzAwOjAw73cO+QAAACV0RVh0ZGF0ZTptb2RpZnkAMjAyMC0wMy0wM1QxMzo1ODo1NiswMDowMJ4qtkUAAAAASUVORK5CYII=\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAMgAAADIAgMAAADQNkYNAAAABGdBTUEAALGPC/xhBQAAACBjSFJNAAB6JgAAgIQAAPoAAACA6AAAdTAAAOpgAAA6mAAAF3CculE8AAAADFBMVEUAAP9yEhL//////wAZPRNOAAAAAWJLR0QCZgt8ZAAAAAd0SU1FB+QIDwcxDg24lJMAAAI7SURBVGje7ZmxkcIwEEUhcAnuhxIIEIGvAwioxiU4QASU4GougA4gsM6SfNxh766sfx7N3Ix+zJv9X7KxtLtaZWVlZWWl1Fr12sQQpXLaRRMRTKFe2kbk+Na8POpNkbbmWlMqtkwxRsJllIotMxQ56F7neWX8LmrdGtNpPWc//Z4cLsbpfp6xN85XNRA904Sd2V/stXlJ16EFcL4O7Q/SnUPOrK/9xfzSvQ44K0dFhjLCmq0nRYYyvDPrq2rfka4RnVlfoyJ9GdGZ9XUdI89aWGYbpTITNUKYgvLlnXFhSsqXd7bjEcKXc8YhvYEjhdzY/Db9lUKebP4+/b6lkK7m8nNRhDAlE8WFoREuig/DpKej+DAbGmGiuDAUUsgItWR8ejZ/yaZ3+SmET+/zL4JIC8Ys2VpI7/JPkSKEbCnkyiNPCilDyI5A+AVzSzZFVAhRBCKssVvlBRB5W8iNwZCThDymiLyT5F5iSCsh3RSRN5/cfgwxoiikkpFmgqgwopZAPmTkcxnkKCO3BZDQs088/RhykpFHRjLyP5AE70uadz/Z/1iK/2TgYwF9klJ8K4GPeJrTBXDsSXMeS3O2RA69wNEaOMAD1wTgMgJceYCLVZobH3IVBS68wLUauLwDLQKgEYG0O4CmCtC6ARpEQBsKaXYBLTWgcQe0B4EmJNLqBBqqQNsWaA4jLWig0Q2005GmPTAaAAYQyJgDGKYgIxtgMISMn4AhFzJKAwZ2yFgQGT4iI86srKysrD/rC4LWcCSWwIp+AAAAJXRFWHRkYXRlOmNyZWF0ZQAyMDIwLTA4LTE1VDA2OjQ5OjE0KzAxOjAwlLKpXwAAACV0RVh0ZGF0ZTptb2RpZnkAMjAyMC0wOC0xNVQwNjo0OToxNCswMTowMOXvEeMAAAAASUVORK5CYII=\n", "text/plain": [ "" ] @@ -1563,7 +1562,7 @@ "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAMgAAADIAgMAAADQNkYNAAAABGdBTUEAALGPC/xhBQAAACBjSFJNAAB6JgAAgIQAAPoAAACA6AAAdTAAAOpgAAA6mAAAF3CculE8AAAADFBMVEUAAP9yEhL//////wAZPRNOAAAAAWJLR0QCZgt8ZAAAAAd0SU1FB+QDAw06OBFn074AAAI7SURBVGje7ZmxkcIwEEUhcAnuhxIIEIGvAwioxiU4QASU4GougA4gsM6SfNxh766sfx7N3Ix+zJv9X7KxtLtaZWVlZWWl1Fr12sQQpXLaRRMRTKFe2kbk+Na8POpNkbbmWlMqtkwxRsJllIotMxQ56F7neWX8LmrdGtNpPWc//Z4cLsbpfp6xN85XNRA904Sd2V/stXlJ16EFcL4O7Q/SnUPOrK/9xfzSvQ44K0dFhjLCmq0nRYYyvDPrq2rfka4RnVlfoyJ9GdGZ9XUdI89aWGYbpTITNUKYgvLlnXFhSsqXd7bjEcKXc8YhvYEjhdzY/Db9lUKebP4+/b6lkK7m8nNRhDAlE8WFoREuig/DpKej+DAbGmGiuDAUUsgItWR8ejZ/yaZ3+SmET+/zL4JIC8Ys2VpI7/JPkSKEbCnkyiNPCilDyI5A+AVzSzZFVAhRBCKssVvlBRB5W8iNwZCThDymiLyT5F5iSCsh3RSRN5/cfgwxoiikkpFmgqgwopZAPmTkcxnkKCO3BZDQs088/RhykpFHRjLyP5AE70uadz/Z/1iK/2TgYwF9klJ8K4GPeJrTBXDsSXMeS3O2RA69wNEaOMAD1wTgMgJceYCLVZobH3IVBS68wLUauLwDLQKgEYG0O4CmCtC6ARpEQBsKaXYBLTWgcQe0B4EmJNLqBBqqQNsWaA4jLWig0Q2005GmPTAaAAYQyJgDGKYgIxtgMISMn4AhFzJKAwZ2yFgQGT4iI86srKysrD/rC4LWcCSWwIp+AAAAJXRFWHRkYXRlOmNyZWF0ZQAyMDIwLTAzLTAzVDEzOjU4OjU2KzAwOjAw73cO+QAAACV0RVh0ZGF0ZTptb2RpZnkAMjAyMC0wMy0wM1QxMzo1ODo1NiswMDowMJ4qtkUAAAAASUVORK5CYII=\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAMgAAADIAgMAAADQNkYNAAAABGdBTUEAALGPC/xhBQAAACBjSFJNAAB6JgAAgIQAAPoAAACA6AAAdTAAAOpgAAA6mAAAF3CculE8AAAADFBMVEUAAP9yEhL//////wAZPRNOAAAAAWJLR0QCZgt8ZAAAAAd0SU1FB+QIDwcxDg24lJMAAAI7SURBVGje7ZmxkcIwEEUhcAnuhxIIEIGvAwioxiU4QASU4GougA4gsM6SfNxh766sfx7N3Ix+zJv9X7KxtLtaZWVlZWWl1Fr12sQQpXLaRRMRTKFe2kbk+Na8POpNkbbmWlMqtkwxRsJllIotMxQ56F7neWX8LmrdGtNpPWc//Z4cLsbpfp6xN85XNRA904Sd2V/stXlJ16EFcL4O7Q/SnUPOrK/9xfzSvQ44K0dFhjLCmq0nRYYyvDPrq2rfka4RnVlfoyJ9GdGZ9XUdI89aWGYbpTITNUKYgvLlnXFhSsqXd7bjEcKXc8YhvYEjhdzY/Db9lUKebP4+/b6lkK7m8nNRhDAlE8WFoREuig/DpKej+DAbGmGiuDAUUsgItWR8ejZ/yaZ3+SmET+/zL4JIC8Ys2VpI7/JPkSKEbCnkyiNPCilDyI5A+AVzSzZFVAhRBCKssVvlBRB5W8iNwZCThDymiLyT5F5iSCsh3RSRN5/cfgwxoiikkpFmgqgwopZAPmTkcxnkKCO3BZDQs088/RhykpFHRjLyP5AE70uadz/Z/1iK/2TgYwF9klJ8K4GPeJrTBXDsSXMeS3O2RA69wNEaOMAD1wTgMgJceYCLVZobH3IVBS68wLUauLwDLQKgEYG0O4CmCtC6ARpEQBsKaXYBLTWgcQe0B4EmJNLqBBqqQNsWaA4jLWig0Q2005GmPTAaAAYQyJgDGKYgIxtgMISMn4AhFzJKAwZ2yFgQGT4iI86srKysrD/rC4LWcCSWwIp+AAAAJXRFWHRkYXRlOmNyZWF0ZQAyMDIwLTA4LTE1VDA2OjQ5OjE0KzAxOjAwlLKpXwAAACV0RVh0ZGF0ZTptb2RpZnkAMjAyMC0wOC0xNVQwNjo0OToxNCswMTowMOXvEeMAAAAASUVORK5CYII=\n", "text/plain": [ "" ] @@ -1574,7 +1573,7 @@ } ], "source": [ - "p.to_ipython_image()" + "plot.to_ipython_image()" ] } ], @@ -1595,7 +1594,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.6.9" + "version": "3.8.5" } }, "nbformat": 4,