Skip to content

Instantly share code, notes, and snippets.

@paulromano
Created January 7, 2025 21:51
Show Gist options
  • Save paulromano/d761c96c22f5693855a492dad943146e to your computer and use it in GitHub Desktop.
Save paulromano/d761c96c22f5693855a492dad943146e to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 13,
"id": "1f50597f-d130-4a91-a975-63b991396158",
"metadata": {},
"outputs": [],
"source": [
"import openmc\n",
"import numpy as np"
]
},
{
"cell_type": "code",
"execution_count": 14,
"id": "a14b5134-a3aa-4495-9ccc-45f7fc16057f",
"metadata": {},
"outputs": [],
"source": [
"v = np.array((-124.0, 250.0, -119.0))\n",
"a1 = np.array((0.0, -190.0, 0.0))\n",
"a2 = np.array((0.0, 0.0, 399.065125))\n",
"a3 = np.array((640.0, 0.0, 0.0))\n",
"box = openmc.model.OrthogonalBox(v, a1, a2, a3)"
]
},
{
"cell_type": "code",
"execution_count": 17,
"id": "d77182cc-7088-4b4d-b26f-8f6a209057bb",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<Axes: xlabel='x [cm]', ylabel='z [cm]'>"
]
},
"execution_count": 17,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAR8AAAEHCAYAAACX0kK7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy80BEi2AAAACXBIWXMAAA9hAAAPYQGoP6dpAAAeZ0lEQVR4nO3dfVRUdf4H8DcIMzzO8DwDhYqiIoJpPuCg+RCsZLitm2seNUNzfSCoEI/HSDesTuKSbR7bfGh3xU5rltSWZlqHYNXUSZPEQJFVV0PBASqZwS1BnO/vj37c04goowxfHt6vc+453O/9zr2fLzhv79w7914nIYQAEVE7c5ZdABF1TwwfIpKC4UNEUjB8iEgKhg8RScHwISIpGD5EJAXDh4ikcJFdQGdjtVpRWVkJb29vODk5yS6HqEMRQqCurg4hISFwdr71vg3Dx06VlZUIDQ2VXQZRh3bhwgXce++9t+zD8LGTt7c3gF9+uRqNRnI1RB2LxWJBaGio8j65FYaPnZo+amk0GoYPUQtac0iCB5yJSAqGDxFJwfAhIik6VfhUVFTg8ccfh7+/P9zd3REdHY2jR48qy4UQeOGFFxAcHAx3d3fEx8fj9OnTNuv48ccfMWvWLGg0Gvj4+GDevHm4cuVKew+FqNvrNOFz+fJljB49Gq6urtizZw9OnjyJ1157Db6+vkqf7OxsrFu3Dhs3bsThw4fh6emJhIQEXL16Vekza9YsnDhxAnl5edi1axf279+PBQsWyBgSUfcmOolly5aJMWPGtLjcarUKvV4vXn31VaWttrZWqNVqsW3bNiGEECdPnhQAxNdff6302bNnj3BychIVFRWtqsNsNgsAwmw23+FIiLoue94fnWbPZ+fOnRg+fDimTZuGoKAgDB06FH/729+U5efOnYPJZEJ8fLzSptVqERMTA6PRCAAwGo3w8fHB8OHDlT7x8fFwdnbG4cOH228wRNR5Pnb997//xYYNG9CvXz98/vnnSE5OxjPPPIO3334bAGAymQAAOp3O5nU6nU5ZZjKZEBQUZLPcxcUFfn5+Sp8b1dfXw2Kx2ExEdPc6zZcMrVYrhg8fjlWrVgEAhg4dipKSEmzcuBFJSUkO225WVhZefPFFh62fqLvqNHs+wcHBiIyMtGkbOHAgysvLAQB6vR4AUFVVZdOnqqpKWabX61FdXW2zvLGxET/++KPS50YZGRkwm83KdOHChTYZD1F312nCZ/To0SgrK7Np+89//oNevXoBAMLCwqDX65Gfn68st1gsOHz4MAwGAwDAYDCgtrYWhYWFSp+CggJYrVbExMTcdLtqtVq5lIKXVBC1oXY4AN4mjhw5IlxcXMQrr7wiTp8+LbZu3So8PDzEP//5T6XP6tWrhY+Pj9ixY4f49ttvxe9+9zsRFhYmfv75Z6XPQw89JIYOHSoOHz4sDhw4IPr16ydmzJjR6jp4touoZfa8PzpN+AghxCeffCKioqKEWq0WERER4q233rJZbrVaxZ/+9Ceh0+mEWq0WcXFxoqyszKbPDz/8IGbMmCG8vLyERqMRc+fOFXV1da2ugeFD1DJ73h9OQvCJpfawWCzQarUwm838CEZ0A3veH53mmA8RdS0MHyKSguFDRFIwfIhICoYPEUnB8CEiKRg+RCQFw4eIpGD4EJEUDB8ikoLhQ0RSMHyISAqGDxFJwfAhIikYPkQkBcOHiKRg+BCRFAwfIpKC4UNEUjB8iEgKhg8RScHwISIpGD5EJAXDh4ikYPgQkRQMHyKSguFDRFIwfIhICoYPEUnRacNn9erVcHJyQlpamtJ29epVpKSkwN/fH15eXpg6dSqqqqpsXldeXo7ExER4eHggKCgIS5cuRWNjYztXT0SdMny+/vprbNq0CYMHD7ZpX7x4MT755BPk5uZi3759qKysxKOPPqosv379OhITE9HQ0IBDhw7h7bffxpYtW/DCCy+09xCISHQydXV1ol+/fiIvL0+MGzdOPPvss0IIIWpra4Wrq6vIzc1V+paWlgoAwmg0CiGE2L17t3B2dhYmk0nps2HDBqHRaER9fX2rtm82mwUAYTab225QRF2EPe+PTrfnk5KSgsTERMTHx9u0FxYW4tq1azbtERER6NmzJ4xGIwDAaDQiOjoaOp1O6ZOQkACLxYITJ060zwCICADgIrsAe7z33nv45ptv8PXXXzdbZjKZoFKp4OPjY9Ou0+lgMpmUPr8OnqblTctupr6+HvX19cq8xWK5myEQ0f/rNHs+Fy5cwLPPPoutW7fCzc2t3bablZUFrVarTKGhoe22baKurNOET2FhIaqrq3H//ffDxcUFLi4u2LdvH9atWwcXFxfodDo0NDSgtrbW5nVVVVXQ6/UAAL1e3+zsV9N8U58bZWRkwGw2K9OFCxfafnBE3VCnCZ+4uDgUFxejqKhImYYPH45Zs2YpP7u6uiI/P195TVlZGcrLy2EwGAAABoMBxcXFqK6uVvrk5eVBo9EgMjLypttVq9XQaDQ2ExHdvU5zzMfb2xtRUVE2bZ6envD391fa582bh/T0dPj5+UGj0eDpp5+GwWDAqFGjAAATJ05EZGQkZs+ejezsbJhMJqxYsQIpKSlQq9XtPiai7qzThE9rvP7663B2dsbUqVNRX1+PhIQErF+/Xlneo0cP7Nq1C8nJyTAYDPD09ERSUhJeeukliVUTdU9OQgghu4jOxGKxQKvVwmw28yMY0Q3seX90mmM+RNS1MHyISAqGDxFJwfAhIikYPkQkBcOHiKRg+BCRFAwfIpKC4UNEUjB8iEgKhg8RScHwISIpGD5EJAXDh4ikYPgQkRQMHyKSguFDRFIwfIhICoYPEUnB8CEiKRg+RCQFw4eIpGjVc7vS09PtXvGKFSvg5+dn9+uIqHto1XO7nJ2dYTAYoFKpWrXSAwcOoKysDH369LnrAjsaPreLqGX2vD9a/cTSjz76CEFBQa3q6+3t3drVUhvYvW+a7BJIoofH5cou4Y606phPTk4OtFptq1e6adMm6HS6Oy6KiLq+Vu35JCUl2bXSmTNn3lExRNR9tPpj181cuXIFVqvVpo3HQYioNew+1X7u3DkkJibC09MTWq0Wvr6+8PX1hY+PD3x9fR1RIxF1QXaHz+OPP47Lly9j8+bNyM/PR0FBAQoKCvDvf/8bBQUFjqgRAJCVlYURI0bA29sbQUFBmDJlCsrKymz6XL16FSkpKfD394eXlxemTp2Kqqoqmz7l5eVITEyEh4cHgoKCsHTpUjQ2NjqsbiK6Obs/dh0/fhyFhYUYMGCAI+pp0b59+5CSkoIRI0agsbERzz//PCZOnIiTJ0/C09MTALB48WJ8+umnyM3NhVarRWpqKh599FEcPHgQAHD9+nUkJiZCr9fj0KFDuHTpEp544gm4urpi1apV7Toeou6uVd/z+bUJEyZg+fLliI+Pd1RNrVJTU4OgoCDs27cPY8eOhdlsRmBgIN5991384Q9/AACcOnUKAwcOhNFoxKhRo7Bnzx5MnjwZlZWVytm4jRs3YtmyZaipqWnV95g64vd8eKq9e+tIp9od8j2fJn//+9+xaNEiVFRUICoqCq6urjbLBw8ebO8q74jZbAYA5VvUhYWFuHbtmk0oRkREoGfPnkr4GI1GREdH23wNICEhAcnJyThx4gSGDh3aLrUT0R2ET01NDc6ePYu5c+cqbU5OThBCwMnJCdevX2/TAm/GarUiLS0No0ePRlRUFADAZDJBpVLBx8fHpq9Op4PJZFL63Pj9o6b5pj43qq+vR319vTJvsVjaahhE3Zrd4fPkk09i6NCh2LZtG3Q6HZycnBxR1y2lpKSgpKQEBw4ccPi2srKy8OKLLzp8O0Tdjd3h891332Hnzp0IDw93RD23lZqail27dmH//v249957lXa9Xo+GhgbU1tba7P1UVVVBr9crfY4cOWKzvqazYU19bpSRkWFzYa3FYkFoaGhbDYeo27L7VPuDDz6I48ePO6KWWxJCIDU1FR999BEKCgoQFhZms3zYsGFwdXVFfn6+0lZWVoby8nIYDAYAgMFgQHFxMaqrq5U+eXl50Gg0iIyMvOl21Wo1NBqNzUREd8/uPZ/f/va3WLx4MYqLixEdHd3sgPMjjzzSZsX9WkpKCt59913s2LED3t7eyjEarVYLd3d3aLVazJs3D+np6fDz84NGo8HTTz8Ng8GAUaNGAQAmTpyIyMhIzJ49G9nZ2TCZTFixYgVSUlKgVqsdUjcR3Zzdp9qdnVveWXLkAeeWji3l5ORgzpw5AH75kuGSJUuwbds21NfXIyEhAevXr7f5SPXdd98hOTkZe/fuhaenJ5KSkrB69Wq4uLQuh3mqnTqaznqq3e7w6e4YPtTRdNbw4W1UiUgKu8PnmWeewbp165q1//Wvf0VaWlpb1ERE3YDd4fPhhx9i9OjRzdpjY2PxwQcftElRRNT12R0+P/zww03vaqjRaPD999+3SVFE1PXZHT7h4eH47LPPmrXv2bOnS94wnogcw+7v+aSnpyM1NRU1NTV48MEHAQD5+fl47bXXsHbt2rauj4i6qDu6tqu+vh6vvPIKXn75ZQBA7969sWHDBjzxxBNtXiARdU13dA/n5ORkJCcno6amBu7u7vDy8mrruoioi7urG8gHBga2VR1E1M206oDz/fffj8uXL7d6pWPGjEFFRcUdF0VEXV+r9nyKiopw/PjxVj97vaioyOYGXEREN2r1x664uDi09jIwGTcYI6LOpVXhc+7cObtX/OsbfRER3ahV4dOrVy9H10FE3QyvaiciKRg+RCQFw4eIpGD4EJEUdodPUlIS9u/f74haiKgbsTt8zGYz4uPj0a9fP6xatYrfZCaiO2J3+Hz88ceoqKhAcnIy3n//ffTu3RuTJk3CBx98gGvXrjmiRiLqgu7omE9gYCDS09Nx/PhxHD58GOHh4Zg9ezZCQkKwePFinD59uq3rJKIu5q6uar906RLy8vKQl5eHHj164OGHH0ZxcTEiIyORnZ2NxYsXt1WddAsd6dEpRK1l957PtWvX8OGHH2Ly5Mno1asXcnNzkZaWhsrKSrz99tv44osvsH37drz00kuOqJeIugi793yCg4NhtVoxY8YMHDlyBEOGDGnWZ8KECfDx8WmD8oioq7I7fF5//XVMmzYNbm5uLfbx8fG5o4tRiaj7sDt8Zs+e7Yg6iKib4TeciUgKhg8RSdFtw+fNN99E79694ebmhpiYGBw5ckR2SUTdSrcMn/fffx/p6enIzMzEN998g/vuuw8JCQmorq6WXRpRt9Etw+cvf/kL5s+fj7lz5yIyMhIbN26Eh4cHNm/eLLs0om6j24VPQ0MDCgsLER8fr7Q5OzsjPj4eRqNRYmVE3ctdXV7RGX3//fe4fv06dDqdTbtOp8OpU6ea9a+vr7d5DJDFYnF4jUTdQbfb87FXVlYWtFqtMoWGhsouiahL6HbhExAQgB49eqCqqsqmvaqqCnq9vln/jIwMmM1mZbpw4UJ7lUrUpXW78FGpVBg2bBjy8/OVNqvVivz8fBgMhmb91Wo1NBqNzUREd6/bHfMBgPT0dCQlJWH48OEYOXIk1q5di//973+YO3eu7NKIuo1uGT7Tp09HTU0NXnjhBZhMJgwZMgSfffZZs4PQROQ4TqK1D2AnAL+c7dJqtTCbzfwIRnQDe94f3e6YDxF1DAwfIpKC4UNEUjB8iEgKhg8RScHwISIpGD5EJAXDh4ikYPgQkRQMHyKSguFDRFIwfIhICoYPEUnB8CEiKRg+RCQFw4eIpGD4EJEUDB8ikoLhQ0RSMHyISAqGDxFJwfAhIikYPkQkBcOHiKRg+BCRFAwfIpKC4UNEUjB8iEgKhg8RSdEpwuf8+fOYN28ewsLC4O7ujr59+yIzMxMNDQ02/b799ls88MADcHNzQ2hoKLKzs5utKzc3FxEREXBzc0N0dDR2797dXsMgol/pFOFz6tQpWK1WbNq0CSdOnMDrr7+OjRs34vnnn1f6WCwWTJw4Eb169UJhYSFeffVVrFy5Em+99ZbS59ChQ5gxYwbmzZuHY8eOYcqUKZgyZQpKSkpkDIuoexOdVHZ2tggLC1Pm169fL3x9fUV9fb3StmzZMjFgwABl/rHHHhOJiYk264mJiRELFy5s9XbNZrMAIMxm811UT9Q12fP+6BR7PjdjNpvh5+enzBuNRowdOxYqlUppS0hIQFlZGS5fvqz0iY+Pt1lPQkICjEZj+xRNRIpOGT5nzpzBG2+8gYULFyptJpMJOp3Opl/TvMlkumWfpuU3U19fD4vFYjMR0d2TGj7PPfccnJycbjmdOnXK5jUVFRV46KGHMG3aNMyfP9/hNWZlZUGr1SpTaGiow7dJ1B24yNz4kiVLMGfOnFv26dOnj/JzZWUlJkyYgNjYWJsDyQCg1+tRVVVl09Y0r9frb9mnafnNZGRkID09XZm3WCwMIKI2IDV8AgMDERgY2Kq+FRUVmDBhAoYNG4acnBw4O9vutBkMBixfvhzXrl2Dq6srACAvLw8DBgyAr6+v0ic/Px9paWnK6/Ly8mAwGFrcrlqthlqttnNkRHRb7XAA/K5dvHhRhIeHi7i4OHHx4kVx6dIlZWpSW1srdDqdmD17tigpKRHvvfee8PDwEJs2bVL6HDx4ULi4uIg1a9aI0tJSkZmZKVxdXUVxcXGra+HZLqKW2fP+6BThk5OTIwDcdPq148ePizFjxgi1Wi3uuecesXr16mbr2r59u+jfv79QqVRi0KBB4tNPP7WrFoYPUcvseX84CSGErL2uzshisUCr1cJsNkOj0cguh6hDsef90SlPtRNR58fwISIpGD5EJAXDh4ikYPgQkRQMHyKSguFDRFIwfIhICoYPEUnB8CEiKRg+RCQFw4eIpGD4EJEUDB8ikoLhQ0RSMHyISAqGDxFJwfAhIikYPkQkBcOHiKRg+BCRFAwfIpJC6hNLO6OmJw1ZLBbJlRB1PE3vi9Y8kYvhY6e6ujoA4PPaiW6hrq4OWq32ln340EA7Wa1WVFZWwtvbG05OTrfsa7FYEBoaigsXLnTqBwxyHB1LRx6HEAJ1dXUICQmBs/Otj+pwz8dOzs7OuPfee+16jUaj6XD/SO4Ex9GxdNRx3G6PpwkPOBORFAwfIpKC4eNAarUamZmZUKvVsku5KxxHx9JVxsEDzkQkBfd8iEgKhg8RScHwISIpGD5EJAXD5y6dP38e8+bNQ1hYGNzd3dG3b19kZmaioaHBpo+Tk1Oz6auvvrJZV25uLiIiIuDm5obo6Gjs3r27vYfTzJtvvonevXvDzc0NMTExOHLkiOySbGRlZWHEiBHw9vZGUFAQpkyZgrKyMps+48ePb/a7X7RokU2f8vJyJCYmwsPDA0FBQVi6dCkaGxvbbRwrV65sVmNERISy/OrVq0hJSYG/vz+8vLwwdepUVFVVdagx2E3QXdmzZ4+YM2eO+Pzzz8XZs2fFjh07RFBQkFiyZInS59y5cwKA+OKLL8SlS5eUqaGhQelz8OBB0aNHD5GdnS1OnjwpVqxYIVxdXUVxcbGMYQkhhHjvvfeESqUSmzdvFidOnBDz588XPj4+oqqqSlpNN0pISBA5OTmipKREFBUViYcfflj07NlTXLlyRekzbtw4MX/+fJvfvdlsVpY3NjaKqKgoER8fL44dOyZ2794tAgICREZGRruNIzMzUwwaNMimxpqaGmX5okWLRGhoqMjPzxdHjx4Vo0aNErGxsR1qDPZi+DhAdna2CAsLU+abwufYsWMtvuaxxx4TiYmJNm0xMTFi4cKFjirztkaOHClSUlKU+evXr4uQkBCRlZUlrabbqa6uFgDEvn37lLZx48aJZ599tsXX7N69Wzg7OwuTyaS0bdiwQWg0GlFfX+/IchWZmZnivvvuu+my2tpa4erqKnJzc5W20tJSAUAYjUYhRMcYg734scsBzGYz/Pz8mrU/8sgjCAoKwpgxY7Bz506bZUajEfHx8TZtCQkJMBqNDq21JQ0NDSgsLLSpydnZGfHx8dJqag2z2QwAzX7/W7duRUBAAKKiopCRkYGffvpJWWY0GhEdHQ2dTqe0JSQkwGKx4MSJE+1TOIDTp08jJCQEffr0waxZs1BeXg4AKCwsxLVr12z+FhEREejZs6fyt+goY7AHLyxtY2fOnMEbb7yBNWvWKG1eXl547bXXMHr0aDg7O+PDDz/ElClT8PHHH+ORRx4BAJhMJpt/OACg0+lgMpnatf4m33//Pa5fv37Tmk6dOiWlptuxWq1IS0vD6NGjERUVpbTPnDkTvXr1QkhICL799lssW7YMZWVl+Ne//gWg5d9907L2EBMTgy1btmDAgAG4dOkSXnzxRTzwwAMoKSmByWSCSqWCj49Psxqb6usIY7AXw6cFzz33HP785z/fsk9paanNQcGKigo89NBDmDZtGubPn6+0BwQEID09XZkfMWIEKisr8eqrryrhQ3cvJSUFJSUlOHDggE37ggULlJ+jo6MRHByMuLg4nD17Fn379m3vMm9q0qRJys+DBw9GTEwMevXqhe3bt8Pd3V1iZY7D8GnBkiVLMGfOnFv26dOnj/JzZWUlJkyYgNjYWLz11lu3XX9MTAzy8vKUeb1e3+zsRVVVFfR6vX2Ft5GAgAD06NGjQ9V0K6mpqdi1axf2799/21uexMTEAPhlL7Vv377Q6/XNzuI1jVvWWH18fNC/f3+cOXMGv/nNb9DQ0IDa2lqbvZ9f/y064hhuh8d8WhAYGIiIiIhbTiqVCsAvezzjx4/HsGHDkJOTc9ubKAFAUVERgoODlXmDwYD8/HybPnl5eTAYDG07sFZSqVQYNmyYTU1WqxX5+fnSaroZIQRSU1Px0UcfoaCgAGFhYbd9TVFREQAov3+DwYDi4mJUV1crffLy8qDRaBAZGemQum/nypUrOHv2LIKDgzFs2DC4urra/C3KyspQXl6u/C064hhuS/YR787u4sWLIjw8XMTFxYmLFy/anCptsmXLFvHuu++K0tJSUVpaKl555RXh7OwsNm/erPQ5ePCgcHFxEWvWrBGlpaUiMzOzQ5xqV6vVYsuWLeLkyZNiwYIFwsfHx+aMimzJyclCq9WKvXv32vzuf/rpJyGEEGfOnBEvvfSSOHr0qDh37pzYsWOH6NOnjxg7dqyyjqbT1BMnThRFRUXis88+E4GBge16mnrJkiVi79694ty5c+LgwYMiPj5eBAQEiOrqaiHEL6fae/bsKQoKCsTRo0eFwWAQBoOhQ43BXgyfu5STkyMA3HRqsmXLFjFw4EDh4eEhNBqNGDlypM1p0ybbt28X/fv3FyqVSgwaNEh8+umn7TmUm3rjjTdEz549hUqlEiNHjhRfffWV7JJstPS7z8nJEUIIUV5eLsaOHSv8/PyEWq0W4eHhYunSpTbf8xFCiPPnz4tJkyYJd3d3ERAQIJYsWSKuXbvWbuOYPn26CA4OFiqVStxzzz1i+vTp4syZM8ryn3/+WTz11FPC19dXeHh4iN///vc2/8F1hDHYi7fUICIpeMyHiKRg+BCRFAwfIpKC4UNEUjB8iEgKhg8RScHwISIpGD7Uof36LpBDhgxx6La2bNmibCstLc2h2yKGD3USX3zxRbNr39ra9OnTcenSpQ517VpXxqvaqVPw9/eHv7+/Q7fh7u4Od3d35YJhcizu+VC7qampgV6vx6pVq5S2Q4cOQaVS3dFezebNmzFo0CCo1WoEBwcjNTVVWebk5IRNmzZh8uTJ8PDwwMCBA2E0GnHmzBmMHz8enp6eiI2NxdmzZ9tkbGQ/hg+1m8DAQGzevBkrV67E0aNHUVdXh9mzZyM1NRVxcXF2rWvDhg1ISUnBggULUFxcjJ07dyI8PNymz8svv4wnnngCRUVFiIiIwMyZM7Fw4UJkZGTg6NGjyu04SBLJF7ZSN/TUU0+J/v37i5kzZ4ro6Ghx9erVFvu2dPP9kJAQsXz58hZfB0CsWLFCmTcajQKA+Mc//qG0bdu2Tbi5uTV77e1uOE9tg3s+1O7WrFmDxsZG5ObmYuvWrVCr1Xa9vrq6GpWVlbfdWxo8eLDyc9P9jKOjo23arl69CovFYtf2qW0wfKjdnT17FpWVlbBarTh//rzdr2/tPY1dXV2Vn52cnFpss1qtdtdAd4/hQ+2qoaEBjz/+OKZPn46XX34Zf/zjH21u/dka3t7e6N27t8NPvZNj8VQ7tavly5fDbDZj3bp18PLywu7du/Hkk09i165ddq1n5cqVWLRoEYKCgjBp0iTU1dXh4MGDePrppx1UObU17vlQu9m7dy/Wrl2Ld955BxqNBs7OznjnnXfw5ZdfYsOGDXatKykpCWvXrsX69esxaNAgTJ48GadPn3ZQ5eQIvI0qdWjnz59HWFgYjh075vDLK5qMHz8eQ4YMwdq1a9tle90V93yoU4iNjUVsbKxDt7F161Z4eXnhyy+/dOh26Bfc86EOrbGxUTkjplarERoa6rBt1dXVKQ/a8/HxQUBAgMO2RQwfIpKEH7uISAqGDxFJwfAhIikYPkQkBcOHiKRg+BCRFAwfIpKC4UNEUjB8iEiK/wN8zucdfGcCqwAAAABJRU5ErkJggg==",
"text/plain": [
"<Figure size 258.065x259.74 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAR8AAAEHCAYAAACX0kK7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy80BEi2AAAACXBIWXMAAA9hAAAPYQGoP6dpAAAc1ElEQVR4nO3dfVBU1/0G8GdRdgFlF3kXA75GUUTTouIS40ukoCF1qNQ40VpMLDF2sb7VGqIRtWlJSKwZbXxJGyHpxKA01SSKaRAjNrqOkWgCBBkxWlBc0IzsalIB4fz+yI87riDuCuxh4fnM3Cn33LP3fg9mn17u3XtWJYQQICJyMBfZBRBR98TwISIpGD5EJAXDh4ikYPgQkRQMHyKSguFDRFIwfIhIip6yC3A2jY2NqKyshKenJ1QqlexyiDoVIQRu3LiBoKAguLi0fm7D8LFTZWUlgoODZZdB1KlVVFTgoYcearUPw8dOnp6eAH785Wq1WsnVEHUuFosFwcHByvukNQwfOzX9qaXVahk+RPdgyyUJXnAmIikYPkQkBcOHiKRg+BCRFAwfIpKC4UNEUjB8iEgKhg8RScHwISIpGD5EJAXDh4ikYPgQkRQMHyKSguFDRFIwfIhICoYPEUnB8CEiKRg+RCQFw4eIpGD4EJEUDB8ikoLhQ0RSMHyISAqGDxFJwfAhIikYPkQkBcOHiKRg+BCRFAwfIpLCacPnlVdegUqlwtKlS5W2W7duwWAwwMfHB71790ZCQgKqqqqsXldeXo64uDh4eHjA398fK1euxO3btx1cPRE5Zfh88cUX2LFjB0aNGmXVvmzZMnz88cfIzs5Gfn4+KisrMXPmTGV7Q0MD4uLiUFdXh+PHj+Odd95BZmYm1q5d6+ghEJFwMjdu3BAPP/ywyM3NFZMmTRJLliwRQghRU1MjXF1dRXZ2ttK3pKREABBGo1EIIUROTo5wcXERJpNJ6bNt2zah1WpFbW2tTcc3m80CgDCbze03KKIuwp73h9Od+RgMBsTFxSE6OtqqvaCgAPX19VbtoaGhCAkJgdFoBAAYjUaEh4cjICBA6RMbGwuLxYLi4uIWj1dbWwuLxWK1EFHb9ZRdgD2ysrLw5Zdf4osvvmi2zWQyQa1Ww8vLy6o9ICAAJpNJ6XNn8DRtb9rWkrS0NKxfv74dqieiOznNmU9FRQWWLFmC9957D25ubg47bkpKCsxms7JUVFQ47NhEXZnThE9BQQGqq6vx05/+FD179kTPnj2Rn5+PzZs3o2fPnggICEBdXR1qamqsXldVVYXAwEAAQGBgYLO7X03rTX3uptFooNVqrRYiajunCZ+pU6eisLAQZ86cUZYxY8Zg7ty5ys+urq7Iy8tTXlNaWory8nLo9XoAgF6vR2FhIaqrq5U+ubm50Gq1GDFihMPHRNSdOc01H09PT4wcOdKqrVevXvDx8VHaFyxYgOXLl8Pb2xtarRaLFy+GXq/H+PHjAQAxMTEYMWIE5s2bh/T0dJhMJqxZswYGgwEajcbhYyLqzpwmfGyxadMmuLi4ICEhAbW1tYiNjcXWrVuV7T169MD+/fuxaNEi6PV69OrVC4mJidiwYYPEqom6J5UQQsguwplYLBbodDqYzWZe/yG6iz3vD6e55kNEXQvDh4ikYPgQkRQMHyKSguFDRFIwfIhICoYPEUnB8CEiKRg+RCQFw4eIpGD4EJEUDB8ikoLhQ0RSMHyISAqGDxFJwfAhIikYPkQkRZeaRpWcR07+LNkltOiJSdmyS+g2eOZDRFIwfIhICoYPEUnB8CEiKRg+RCQFw4eIpGD4EJEUDB8ikoLhQ0RSMHyISAqGDxFJ4TThk5aWhrFjx8LT0xP+/v6Ij49HaWmpVZ9bt27BYDDAx8cHvXv3RkJCAqqqqqz6lJeXIy4uDh4eHvD398fKlStx+/ZtRw6FiOBE4ZOfnw+DwYATJ04gNzcX9fX1iImJwffff6/0WbZsGT7++GNkZ2cjPz8flZWVmDlzprK9oaEBcXFxqKurw/Hjx/HOO+8gMzMTa9eulTEkom5NJYQQsot4EFevXoW/vz/y8/MxceJEmM1m+Pn5YdeuXfjlL38JADh79iyGDx8Oo9GI8ePH4+DBg3jyySdRWVmJgIAAAMD27duxatUqXL16FWq1+r7HtVgs0Ol0MJvN0Gq1HTrGroxPtXdN9rw/nObM525msxkA4O3tDQAoKChAfX09oqOjlT6hoaEICQmB0WgEABiNRoSHhyvBAwCxsbGwWCwoLi5u8Ti1tbWwWCxWCxG1nVOGT2NjI5YuXYpHH30UI0eOBACYTCao1Wp4eXlZ9Q0ICIDJZFL63Bk8TdubtrUkLS0NOp1OWYKDg9t5NETdk1OGj8FgQFFREbKysjr8WCkpKTCbzcpSUVHR4cck6g6cbibD5ORk7N+/H0ePHsVDDz2ktAcGBqKurg41NTVWZz9VVVUIDAxU+pw8edJqf013w5r63E2j0UCj0bTzKIjIac58hBBITk7G3r17cfjwYQwcONBqe0REBFxdXZGXl6e0lZaWory8HHq9HgCg1+tRWFiI6upqpU9ubi60Wi1GjBjhmIEQEQAbz3zuvF1tq+3bt8Pf39/u192LwWDArl278OGHH8LT01O5RqPT6eDu7g6dTocFCxZg+fLl8Pb2hlarxeLFi6HX6zF+/HgAQExMDEaMGIF58+YhPT0dJpMJa9asgcFg4NkNkYPZFD779u3DU089BXd3d5t2umvXLty8ebNdw2fbtm0AgMmTJ1u1Z2RkYP78+QCATZs2wcXFBQkJCaitrUVsbCy2bt2q9O3Rowf279+PRYsWQa/Xo1evXkhMTMSGDRvarU4iso1Nn/NxcXGByWSyOUw8PT3x1VdfYdCgQW0usLPh53zaBz/n0zW1++d8PvvsM+XzNLY4ePAg+vXrZ3N/Iup+bPqza9KkSXbtdMKECQ9UDBF1Hw98q726uhrV1dVobGy0ah81alSbiyKirs/u8CkoKEBiYiJKSkrQdLlIpVJBCAGVSoWGhoZ2L5KIuh67w+fZZ5/F0KFD8fbbbyMgIAAqlaoj6iKiLs7u8Pn222/xwQcfYMiQIR1RDxF1E3Z/wnnq1Kn46quvOqIWIupG7D7z+fvf/47ExEQUFRVh5MiRcHV1tdo+Y8aMdiuOiLouu8PHaDTi2LFjOHjwYLNtvOBMRLay+8+uxYsX41e/+hWuXLmCxsZGq4XBQ0S2sjt8vvvuOyxbtqzZpFxERPawO3xmzpyJzz77rCNqIaJuxO5rPkOHDkVKSgo+//xzhIeHN7vg/Lvf/a7diiOirsvub6+4exIvq52pVPj222/bXFRnxqfa2wefau+a7Hl/2H3mc+HChQcujIioidNMo0pEXYvd4ZOQkIBXX321WXt6ejpmzeqcp9JE1PnYHT5Hjx7FE0880ax9+vTpOHr0aLsURURdn93hc/PmzRa/VtjV1ZXf5klENrM7fMLDw7F79+5m7VlZWfz6GSKymd13u1566SXMnDkT58+fx+OPPw4AyMvLw/vvv4/sbN6mJCLb2B0+P//5z7Fv3z78+c9/xj//+U+4u7tj1KhROHTokN1zPRNR9/VAczjHxcUhLi6uvWshom6En/MhIilsCh9vb29cu3bN5p2GhITgv//97wMXRURdn01/dtXU1ODgwYPQ6XQ27fS7777j3D5E1Cqbr/kkJiZ2ZB1E1M3YFD53fzEgEVFb8YIzEUnRbcPnzTffxIABA+Dm5obIyEicPHlSdklE3Uq3DJ/du3dj+fLlSE1NxZdffonRo0cjNjYW1dXVsksj6ja6Zfj85S9/QVJSEp555hmMGDEC27dvh4eHB3bu3Cm7NKJuo9uFT11dHQoKChAdHa20ubi4IDo6GkajsVn/2tpaWCwWq4WI2s7u8Hn88cexfv36Zu3Xr19XHjTtzK5du4aGhoZmX/0TEBAAk8nUrH9aWhp0Op2yBAcHO6pUoi7N7vA5cuQI/vrXvyI+Ph7ff/+90l5XV4f8/Px2La4zSElJgdlsVpaKigrZJRF1CQ/0Z9ehQ4dgMpkwfvx4XLx4sZ1L6li+vr7o0aMHqqqqrNqrqqoQGBjYrL9Go4FWq7VaiKjtHih8+vbti/z8fISHh2Ps2LE4cuRIO5fVcdRqNSIiIpCXl6e0NTY2Ii8vD3q9XmJlRN2L3VNqqFQqAD+eEezatQsvv/wypk2bhlWrVrV7cR1l+fLlSExMxJgxYzBu3Di88cYb+P777/HMM8/ILo2o27A7fO7+jsE1a9Zg+PDhTvXs1+zZs3H16lWsXbsWJpMJjzzyCD755BN+/zyRAz3Qlwb6+flZtSUkJCA0NBSnTp1qt8I6WnJyMpKTk2WXQdRt2R0+/fv3b7E9LCwMYWFhbS6IiLqHbvchQyLqHBg+RCQFw4eIpGD4EJEUDB8ikoLhQ0RSMHyISAqGDxFJwfAhIikYPkQkBcOHiKRg+BCRFAwfIpKC4UNEUjB8iEgKhg8RSWH3ZGJE7eGJSdmySyDJeOZDRFIwfIhICoYPEUnB8CEiKRg+RCQFw4eIpGD4EJEUDB8ikoLhQ0RSMHyISAqGDxFJ4RThc/HiRSxYsAADBw6Eu7s7Bg8ejNTUVNTV1Vn1+/rrr/HYY4/Bzc0NwcHBSE9Pb7av7OxshIaGws3NDeHh4cjJyXHUMIjoDk4RPmfPnkVjYyN27NiB4uJibNq0Cdu3b8eLL76o9LFYLIiJiUH//v1RUFCA1157DevWrcNbb72l9Dl+/DiefvppLFiwAKdPn0Z8fDzi4+NRVFQkY1hE3ZtwUunp6WLgwIHK+tatW0WfPn1EbW2t0rZq1SoxbNgwZf2pp54ScXFxVvuJjIwUCxcutPm4ZrNZABBms7kN1RN1Tfa8P5zizKclZrMZ3t7eyrrRaMTEiROhVquVttjYWJSWluL69etKn+joaKv9xMbGwmg03vM4tbW1sFgsVgsRtZ1Thk9ZWRm2bNmChQsXKm0mkwkBAQFW/ZrWTSZTq32atrckLS0NOp1OWYKDg9trGETdmtTweeGFF6BSqVpdzp49a/Way5cvY9q0aZg1axaSkpI6vMaUlBSYzWZlqaio6PBjEnUHUmcyXLFiBebPn99qn0GDBik/V1ZWYsqUKYiKirK6kAwAgYGBqKqqsmprWg8MDGy1T9P2lmg0Gmg0mvuOhYjsIzV8/Pz84OfnZ1Pfy5cvY8qUKYiIiEBGRgZcXKxP2vR6PVavXo36+nq4uroCAHJzczFs2DD06dNH6ZOXl4elS5cqr8vNzYVer2+fARGR7RxwAbzNLl26JIYMGSKmTp0qLl26JK5cuaIsTWpqakRAQICYN2+eKCoqEllZWcLDw0Ps2LFD6XPs2DHRs2dP8frrr4uSkhKRmpoqXF1dRWFhoc218G4X0b3Z8/5wivDJyMgQAFpc7vTVV1+JCRMmCI1GI/r16ydeeeWVZvvas2ePGDp0qFCr1SIsLEwcOHDArloYPkT3Zs/7QyWEELLOupyRxWKBTqeD2WyGVquVXQ5Rp2LP+8Mpb7UTkfNj+BCRFAwfIpKC4UNEUjB8iEgKhg8RScHwISIpGD5EJAXDh4ikYPgQkRQMHyKSguFDRFIwfIhICoYPEUnB8CEiKRg+RCQFw4eIpGD4EJEUDB8ikoLhQ0RSMHyISAqGDxFJwfAhIikYPkQkBcOHiKRg+BCRFAwfIpKC4UNEUjB8iEgKpwuf2tpaPPLII1CpVDhz5ozVtq+//hqPPfYY3NzcEBwcjPT09Gavz87ORmhoKNzc3BAeHo6cnBwHVU5Ed3K68PnDH/6AoKCgZu0WiwUxMTHo378/CgoK8Nprr2HdunV46623lD7Hjx/H008/jQULFuD06dOIj49HfHw8ioqKHDkEIgIA4URycnJEaGioKC4uFgDE6dOnlW1bt24Vffr0EbW1tUrbqlWrxLBhw5T1p556SsTFxVntMzIyUixcuNDmGsxmswAgzGbzgw+EqIuy5/3hNGc+VVVVSEpKwj/+8Q94eHg02240GjFx4kSo1WqlLTY2FqWlpbh+/brSJzo62up1sbGxMBqN9zxubW0tLBaL1UJEbecU4SOEwPz58/H8889jzJgxLfYxmUwICAiwamtaN5lMrfZp2t6StLQ06HQ6ZQkODm7LUIjo/0kNnxdeeAEqlarV5ezZs9iyZQtu3LiBlJQUh9eYkpICs9msLBUVFQ6vgagr6inz4CtWrMD8+fNb7TNo0CAcPnwYRqMRGo3GatuYMWMwd+5cvPPOOwgMDERVVZXV9qb1wMBA5X9b6tO0vSUajabZcYmo7aSGj5+fH/z8/O7bb/PmzXj55ZeV9crKSsTGxmL37t2IjIwEAOj1eqxevRr19fVwdXUFAOTm5mLYsGHo06eP0icvLw9Lly5V9pWbmwu9Xm9zzUIIAOC1H6IWNL0vmt4nreroq98d4cKFC83udtXU1IiAgAAxb948UVRUJLKysoSHh4fYsWOH0ufYsWOiZ8+e4vXXXxclJSUiNTVVuLq6isLCQpuPXVFRIQBw4cKllaWiouK+7yWpZz7tSafT4dNPP4XBYEBERAR8fX2xdu1aPPfcc0qfqKgo7Nq1C2vWrMGLL76Ihx9+GPv27cPIkSNtPk5QUBAqKirg6ekJlUrVppotFguCg4NRUVEBrVbbpn3JxHF0LjLHIYTAjRs3Wvws3t1UQthyfkQdwWKxQKfTwWw2O/1/7BxH5+Es43CKW+1E1PUwfIhICoaPRBqNBqmpqU5/K5/j6FycZRy85kNEUvDMh4ikYPgQkRQMHyKSguFDRFIwfBzg4sWLWLBgAQYOHAh3d3cMHjwYqampqKurs+rT0lP9J06csNpXZ5wG9s0338SAAQPg5uaGyMhInDx5UnZJirS0NIwdOxaenp7w9/dHfHw8SktLrfpMnjy52e/9+eeft+pTXl6OuLg4eHh4wN/fHytXrsTt27cdNo5169Y1qzE0NFTZfuvWLRgMBvj4+KB3795ISEho9hC17DE0Y/NDTfTADh48KObPny/+/e9/i/Pnz4sPP/xQ+Pv7ixUrVih9mp5XO3TokLhy5Yqy1NXVKX2OHTsmevToIdLT08U333wj1qxZY/ezae0tKytLqNVqsXPnTlFcXCySkpKEl5eXqKqqklbTnWJjY0VGRoYoKioSZ86cEU888YQICQkRN2/eVPpMmjRJJCUlWf3e75yJ7/bt22LkyJEiOjpanD59WuTk5AhfX1+RkpLisHGkpqaKsLAwqxqvXr2qbH/++edFcHCwyMvLE6dOnRLjx48XUVFRnWoMd2P4SJKeni4GDhyorLf0sOzd2mMa2PY2btw4YTAYlPWGhgYRFBQk0tLSpNXUmurqagFA5OfnK22TJk0SS5YsuedrcnJyhIuLizCZTErbtm3bhFartZq2tyOlpqaK0aNHt7itpqZGuLq6iuzsbKWtpKREABBGo1EI0TnGcDf+2SWJ2WyGt7d3s/YZM2bA398fEyZMwEcffWS17UGmge1IdXV1KCgosKrJxcUF0dHR0mq6H7PZDADNfvfvvfcefH19MXLkSKSkpOCHH35QthmNRoSHh1vNghkbGwuLxYLi4mLHFA7g3LlzCAoKwqBBgzB37lyUl5cDAAoKClBfX2/17xAaGoqQkBDl36GzjOFOXeapdmdSVlaGLVu24PXXX1faevfujY0bN+LRRx+Fi4sLPvjgA8THx2Pfvn2YMWMGgAebBrYjXbt2DQ0NDS3WdPbsWSk1taaxsRFLly7Fo48+ajWTwZw5c9C/f38EBQXh66+/xqpVq1BaWop//etfAGyborejRUZGIjMzE8OGDcOVK1ewfv16PPbYYygqKoLJZIJarYaXl1ezGu83hbAjx3A3hk8bvPDCC3j11Vdb7VNSUmJ1YfDy5cuYNm0aZs2ahaSkJKXd19cXy5cvV9bHjh2LyspKvPbaa0r4UNsYDAYUFRXh888/t2q/c9qV8PBw9O3bF1OnTsX58+cxePBgR5fZounTpys/jxo1CpGRkejfvz/27NkDd3d3iZU9OIZPG9g6DWyTyspKTJkyBVFRUVbfJ3YvkZGRyM3NVdYfZBrYjuTr64sePXp0qpruJTk5Gfv378fRo0fx0EMPtdq3aXbMsrIyDB48GIGBgc3u4N09Ra+jeXl5YejQoSgrK8PPfvYz1NXVoaamxurs585/h844Bl7zaQM/Pz+Ehoa2ujR9lc/ly5cxefJkREREICMjAy4u9//VnzlzBn379lXWm6aBvZO908C2J7VajYiICKuaGhsbkZeXJ62muwkhkJycjL179+Lw4cMYOHDgfV/T9E24Tb97vV6PwsJCVFdXK31yc3Oh1WoxYsSIDqn7fm7evInz58+jb9++iIiIgKurq9W/Q2lpKcrLy5V/h844Bt7tcoBLly6JIUOGiKlTp4pLly5Z3S5tkpmZKXbt2iVKSkpESUmJ+NOf/iRcXFzEzp07lT7tMQ1se8vKyhIajUZkZmaKb775Rjz33HPCy8vL6q6KTIsWLRI6nU4cOXLE6vf+ww8/CCGEKCsrExs2bBCnTp0SFy5cEB9++KEYNGiQmDhxorKPptvUMTEx4syZM+KTTz4Rfn5+Dr1NvWLFCnHkyBFx4cIFcezYMREdHS18fX1FdXW1EOLHW+0hISHi8OHD4tSpU0Kv1wu9Xt+pxnA3ho8DZGRk3HOu2yaZmZli+PDhwsPDQ2i1WjFu3DirW6dN9uzZI4YOHSrUarUICwsTBw4ccORQWrRlyxYREhIi1Gq1GDdunDhx4oTskhT3+r1nZGQIIYQoLy8XEydOFN7e3kKj0YghQ4aIlStXNvvGzYsXL4rp06cLd3d34evrK1asWCHq6+sdNo7Zs2eLvn37CrVaLfr16ydmz54tysrKlO3/+9//xG9/+1vRp08f4eHhIX7xi19Y/Z9bZxjD3TilBhFJwWs+RCQFw4eIpGD4EJEUDB8ikoLhQ0RSMHyISAqGDxFJwfChTu3OGR4feeSRDj1WZmamcqylS5d26LGI4UNO4tChQ82ea2tvs2fPxpUrVzrNc2ldHZ9qJ6fg4+MDHx+fDj2Gu7s73N3dlYeBqWPxzIcc5t1334WPjw9qa2ut2uPj4zFv3jy797dz506EhYVBo9Ggb9++SE5OVrapVCrs2LEDTz75JDw8PDB8+HAYjUaUlZVh8uTJ6NWrF6KionD+/Pk2j4seDMOHHGbWrFloaGiwmh62uroaBw4cwLPPPmvXvrZt2waDwYDnnnsOhYWF+OijjzBkyBCrPn/84x/x61//GmfOnEFoaCjmzJmDhQsXIiUlBadOnVKm2yBJpD3SSt3SokWLxPTp05X1jRs3ikGDBonGxsYW+99rYv2goCCxevXqex4HgFizZo2ybjQaBQDx9ttvK23vv/++cHNza/ba+00oT+2DZz7kUElJSfj0009x+fJlAD/eYZo/fz5UKpXN+6iurkZlZSWmTp3aar9Ro0YpPzfNVxweHm7VduvWLVgsFnuGQO2EF5zJoX7yk59g9OjRePfddxETE4Pi4mIcOHDArn3YOmexq6ur8nNTuLXU1tjYaNfxqX3wzIcc7je/+Q0yMzORkZGB6OhoBAcH2/V6T09PDBgwoMNvvVPHYviQw82ZMweXLl3C3/72N7svNDdZt24dNm7ciM2bN+PcuXP48ssvsWXLlnaulDoSw4ccTqfTISEhAb1790Z8fPwD7SMxMRFvvPEGtm7dirCwMDz55JM4d+5c+xZKHYrXfEiKy5cvY+7cudBoNA+8j4ULF2LhwoUtbhN3zQ48YMCAZm2TJ09u1kaOwzMfcqjr169j7969OHLkCAwGg82vi4qKQlRUVAdW9uNXJvfu3Rv/+c9/OvQ49CNOIE8ONWDAAFy/fh0vvfQSfv/739+3/+3bt3Hx4kUAgEajsfvitD1u3LihfJGel5cXfH19O+xYxPAhIkn4ZxcRScHwISIpGD5EJAXDh4ikYPgQkRQMHyKSguFDRFIwfIhICoYPEUnxf3uQJ2me3b8fAAAAAElFTkSuQmCC",
"text/plain": [
"<Figure size 258.065x259.74 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAR8AAAEHCAYAAACX0kK7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy80BEi2AAAACXBIWXMAAA9hAAAPYQGoP6dpAAAcxElEQVR4nO3de1BTZ/4G8CdREkBJELlJBUWtooh1FhWD1ktlRYvbYWVbpxcXW9dLG1xv4yrVFXvZ0qK7Orr10tkKjqOlZbu1rcXuIlZsNY6VigWKTKG6oBiwnZJguwLK+/ujP84YQUyE8BJ4PjNnSt7z5pzvG8zTwzk5b1RCCAEiok6mll0AEfVMDB8ikoLhQ0RSMHyISAqGDxFJwfAhIikYPkQkBcOHiKToLbsAV9PU1ISqqip4eXlBpVLJLoeoSxFCoK6uDkFBQVCr2z62Yfg4qKqqCsHBwbLLIOrSKisrMXDgwDb7MHwc5OXlBeCXF1en00muhqhrsVqtCA4OVt4nbWH4OKj5Ty2dTsfwIboLe05J8IQzEUnB8CEiKRg+RCQFw4eIpGD4EJEUDB8ikoLhQ0RSMHyISAqGDxFJwfAhIikYPkQkBcOHiKRg+BCRFAwfIpKC4UNEUjB8iEgKhg8RScHwISIpGD5EJAXDh4ikYPgQkRQMHyKSguFDRFIwfIhICoYPEUnB8CEiKRg+RCQFw4eIpGD4EJEULhs+r7/+OlQqFVasWKG03bhxA0ajEf3790ffvn2RkJCA6upqm+dVVFQgLi4Onp6e8Pf3x5o1a3Dz5s1Orp6IXDJ8vvzyS+zZswdjxoyxaV+5ciU+/vhjZGVlIS8vD1VVVZg7d66y/tatW4iLi0NDQwNOnTqFffv2ISMjAxs3buzsIRCRcDF1dXXiwQcfFDk5OWLq1Kli+fLlQgghamtrhZubm8jKylL6lpSUCADCZDIJIYTIzs4WarVamM1mpc+uXbuETqcT9fX1du3fYrEIAMJisXTcoIi6CUfeHy535GM0GhEXF4eYmBib9vz8fDQ2Ntq0h4WFISQkBCaTCQBgMpkQERGBgIAApU9sbCysViuKi4tb3V99fT2sVqvNQkTt11t2AY7IzMzEV199hS+//LLFOrPZDI1GA29vb5v2gIAAmM1mpc/twdO8vnlda1JTU/HSSy91QPVEdDuXOfKprKzE8uXLceDAAbi7u3fafpOTk2GxWJSlsrKy0/ZN1J25TPjk5+ejpqYGv/rVr9C7d2/07t0beXl52L59O3r37o2AgAA0NDSgtrbW5nnV1dUIDAwEAAQGBra4+tX8uLnPnbRaLXQ6nc1CRO3nMuEzY8YMFBYWoqCgQFnGjRuHp59+WvnZzc0Nubm5ynNKS0tRUVEBg8EAADAYDCgsLERNTY3SJycnBzqdDqNGjer0MRH1ZC5zzsfLywujR4+2aevTpw/69++vtC9cuBCrVq2Cj48PdDodli1bBoPBgIkTJwIAZs6ciVGjRmH+/PlIS0uD2WzGhg0bYDQaodVqO31MRD2Zy4SPPbZu3Qq1Wo2EhATU19cjNjYWO3fuVNb36tULhw8fxvPPPw+DwYA+ffogMTERL7/8ssSqiXomlRBCyC7ClVitVuj1elgsFp7/IbqDI+8PlznnQ0TdC8OHiKRg+BCRFAwfIpKC4UNEUjB8iEgKhg8RScHwISIpGD5EJAXDh4ikYPgQkRQMHyKSguFDRFIwfIhICoYPEUnB8CEiKRg+RCRFt5pGtafKzntcdgkk0aNTs2SXcF945ENEUjB8iEgKhg8RScHwISIpGD5EJAXDh4ikYPgQkRQMHyKSguFDRFIwfIhICoYPEUnhMuGTmpqK8ePHw8vLC/7+/oiPj0dpaalNnxs3bsBoNKJ///7o27cvEhISUF1dbdOnoqICcXFx8PT0hL+/P9asWYObN2925lCICC4UPnl5eTAajTh9+jRycnLQ2NiImTNn4qefflL6rFy5Eh9//DGysrKQl5eHqqoqzJ07V1l/69YtxMXFoaGhAadOncK+ffuQkZGBjRs3yhgSUY+mEkII2UXcj2vXrsHf3x95eXmYMmUKLBYL/Pz8cPDgQfzud78DAFy4cAEjR46EyWTCxIkTceTIEcyZMwdVVVUICAgAAOzevRtr167FtWvXoNFo7rlfq9UKvV4Pi8UCnU7n1DHai3e192xd6a52R94fLnPkcyeLxQIA8PHxAQDk5+ejsbERMTExSp+wsDCEhITAZDIBAEwmEyIiIpTgAYDY2FhYrVYUFxe3up/6+npYrVabhYjazyXDp6mpCStWrMCkSZMwevRoAIDZbIZGo4G3t7dN34CAAJjNZqXP7cHTvL55XWtSU1Oh1+uVJTg4uINHQ9QzuWT4GI1GFBUVITMz0+n7Sk5OhsViUZbKykqn75OoJ3C5mQyTkpJw+PBhnDhxAgMHDlTaAwMD0dDQgNraWpujn+rqagQGBip9zpw5Y7O95qthzX3upNVqodVqO3gUROQyRz5CCCQlJeGDDz7AsWPHEBoaarM+MjISbm5uyM3NVdpKS0tRUVEBg8EAADAYDCgsLERNTY3SJycnBzqdDqNGjeqcgRARADuPfG6/XG2v3bt3w9/f3+Hn3Y3RaMTBgwfx4YcfwsvLSzlHo9fr4eHhAb1ej4ULF2LVqlXw8fGBTqfDsmXLYDAYMHHiRADAzJkzMWrUKMyfPx9paWkwm83YsGEDjEYjj26IOpld4XPo0CE88cQT8PDwsGujBw8exPXr1zs0fHbt2gUAmDZtmk17eno6FixYAADYunUr1Go1EhISUF9fj9jYWOzcuVPp26tXLxw+fBjPP/88DAYD+vTpg8TERLz88ssdVicR2ceuz/mo1WqYzWa7w8TLywvnz5/HkCFD2l1gV8PP+VBX060/5/PZZ58pn6exx5EjR/DAAw/Y3Z+Ieh67/uyaOnWqQxudPHnyfRVDRD3HfV9qr6mpQU1NDZqammzax4wZ0+6iiKj7czh88vPzkZiYiJKSEjSfLlKpVBBCQKVS4datWx1eJBF1Pw6Hz3PPPYfhw4fj7bffRkBAAFQqlTPqIqJuzuHw+e677/D+++9j2LBhzqiHiHoIhz/hPGPGDJw/f94ZtRBRD+Lwkc8//vEPJCYmoqioCKNHj4abm5vN+scee6zDiiOi7svh8DGZTDh58iSOHDnSYh1POBORvRz+s2vZsmV45plncPXqVTQ1NdksDB4ispfD4fPDDz9g5cqVLSblIiJyhMPhM3fuXHz22WfOqIWIehCHz/kMHz4cycnJ+OKLLxAREdHihPMf//jHDiuOiLovh7+94s5JvGw2plLhu+++a3dRXRnvaqeuxlXvanf4yOfixYv3XRgRUTOXmUaViLoXh8MnISEBb7zxRov2tLQ0PP44D/+JyD4Oh8+JEyfw6KOPtmifPXs2Tpw40SFFEVH353D4XL9+vdWvFXZzc+O3eRKR3RwOn4iICLz77rst2jMzM/n1M0RkN4evdv35z3/G3LlzUV5ejkceeQQAkJubi3feeQdZWV3nkh8RdW0Oh89vfvMbHDp0CK+99hr++c9/wsPDA2PGjMHRo0cdnuuZiHqu+5rDOS4uDnFxcR1dCxH1IPycDxFJYVf4+Pj44Pvvv7d7oyEhIfjvf/9730URUfdn159dtbW1OHLkCPR6vV0b/eGHHzi3DxG1ye5zPomJic6sg4h6GLvC584vBiQiai+ecCYiKXps+Lz55psYPHgw3N3dERUVhTNnzsguiahH6ZHh8+6772LVqlVISUnBV199hYceegixsbGoqamRXRpRj9Ejw+dvf/sbFi1ahGeffRajRo3C7t274enpib1798oujajH6HHh09DQgPz8fMTExChtarUaMTExMJlMLfrX19fDarXaLETUfg6HzyOPPIKXXnqpRfuPP/6o3GjalX3//fe4detWi6/+CQgIgNlsbtE/NTUVer1eWYKDgzurVKJuzeHwOX78OP7+978jPj4eP/30k9Le0NCAvLy8Di2uK0hOTobFYlGWyspK2SURdQv39WfX0aNHYTabMXHiRFy6dKmDS3IuX19f9OrVC9XV1Tbt1dXVCAwMbNFfq9VCp9PZLETUfvcVPgMGDEBeXh4iIiIwfvx4HD9+vIPLch6NRoPIyEjk5uYqbU1NTcjNzYXBYJBYGVHP4vCUGiqVCsAvRwQHDx7Eq6++ilmzZmHt2rUdXpyzrFq1ComJiRg3bhwmTJiAbdu24aeffsKzzz4ruzSiHsPh8LnzOwY3bNiAkSNHutS9X/PmzcO1a9ewceNGmM1mjB07Fp9++im/f56oE93Xlwb6+fnZtCUkJCAsLAxnz57tsMKcLSkpCUlJSbLLIOqxHA6fQYMGtdoeHh6O8PDwdhdERD1Dj/uQIRF1DQwfIpKC4UNEUjB8iEgKhg8RScHwISIpGD5EJAXDh4ikYPgQkRQMHyKSguFDRFIwfIhICoYPEUnB8CEiKRg+RCQFw4eIpHB4MjHqeh6dmiW7BCKH8ciHiKRg+BCRFAwfIpKC4UNEUjB8iEgKhg8RScHwISIpGD5EJAXDh4ikYPgQkRQMHyKSwiXC59KlS1i4cCFCQ0Ph4eGBoUOHIiUlBQ0NDTb9vv76azz88MNwd3dHcHAw0tLSWmwrKysLYWFhcHd3R0REBLKzsztrGER0G5cInwsXLqCpqQl79uxBcXExtm7dit27d+PFF19U+litVsycORODBg1Cfn4+Nm/ejE2bNuGtt95S+pw6dQpPPvkkFi5ciHPnziE+Ph7x8fEoKiqSMSyink24qLS0NBEaGqo83rlzp+jXr5+or69X2tauXStGjBihPH7iiSdEXFyczXaioqLEkiVL7N6vxWIRAITFYmlH9UTdkyPvD5c48mmNxWKBj4+P8thkMmHKlCnQaDRKW2xsLEpLS/Hjjz8qfWJiYmy2ExsbC5PJdNf91NfXw2q12ixE1H4uGT5lZWXYsWMHlixZorSZzWYEBATY9Gt+bDab2+zTvL41qamp0Ov1yhIcHNxRwyDq0aSGz7p166BSqdpcLly4YPOcK1euYNasWXj88cexaNEip9eYnJwMi8WiLJWVlU7fJ1FPIHUmw9WrV2PBggVt9hkyZIjyc1VVFaZPn47o6GibE8kAEBgYiOrqapu25seBgYFt9mle3xqtVgutVnvPsRCRY6SGj5+fH/z8/Ozqe+XKFUyfPh2RkZFIT0+HWm170GYwGLB+/Xo0NjbCzc0NAJCTk4MRI0agX79+Sp/c3FysWLFCeV5OTg4MBkPHDIiI7NcJJ8Db7fLly2LYsGFixowZ4vLly+Lq1avK0qy2tlYEBASI+fPni6KiIpGZmSk8PT3Fnj17lD4nT54UvXv3Flu2bBElJSUiJSVFuLm5icLCQrtr4dUuortz5P3hEuGTnp4uALS63O78+fNi8uTJQqvVigceeEC8/vrrLbb13nvvieHDhwuNRiPCw8PFJ5984lAtDB+iu3Pk/aESQghZR12uyGq1Qq/Xw2KxQKfTyS6HqEtx5P3hkpfaicj1MXyISAqGDxFJwfAhIikYPkQkBcOHiKRg+BCRFAwfIpKC4UNEUjB8iEgKhg8RScHwISIpGD5EJAXDh4ikYPgQkRQMHyKSguFDRFIwfIhICoYPEUnB8CEiKRg+RCQFw4eIpGD4EJEUDB8ikoLhQ0RSMHyISAqGDxFJwfAhIikYPkQkhcuFT319PcaOHQuVSoWCggKbdV9//TUefvhhuLu7Izg4GGlpaS2en5WVhbCwMLi7uyMiIgLZ2dmdVDkR3c7lwudPf/oTgoKCWrRbrVbMnDkTgwYNQn5+PjZv3oxNmzbhrbfeUvqcOnUKTz75JBYuXIhz584hPj4e8fHxKCoq6swhEBEACBeSnZ0twsLCRHFxsQAgzp07p6zbuXOn6Nevn6ivr1fa1q5dK0aMGKE8fuKJJ0RcXJzNNqOiosSSJUvsrsFisQgAwmKx3P9AiLopR94fLnPkU11djUWLFmH//v3w9PRssd5kMmHKlCnQaDRKW2xsLEpLS/Hjjz8qfWJiYmyeFxsbC5PJdNf91tfXw2q12ixE1H4uET5CCCxYsABLly7FuHHjWu1jNpsREBBg09b82Gw2t9mneX1rUlNTodfrlSU4OLg9QyGi/yc1fNatWweVStXmcuHCBezYsQN1dXVITk7u9BqTk5NhsViUpbKystNrIOqOesvc+erVq7FgwYI2+wwZMgTHjh2DyWSCVqu1WTdu3Dg8/fTT2LdvHwIDA1FdXW2zvvlxYGCg8t/W+jSvb41Wq22xXyJqP6nh4+fnBz8/v3v22759O1599VXlcVVVFWJjY/Huu+8iKioKAGAwGLB+/Xo0NjbCzc0NAJCTk4MRI0agX79+Sp/c3FysWLFC2VZOTg4MBoPdNQshAIDnfoha0fy+aH6ftMnZZ7+d4eLFiy2udtXW1oqAgAAxf/58UVRUJDIzM4Wnp6fYs2eP0ufkyZOid+/eYsuWLaKkpESkpKQINzc3UVhYaPe+KysrBQAuXLi0sVRWVt7zvST1yKcj6fV6/Oc//4HRaERkZCR8fX2xceNGLF68WOkTHR2NgwcPYsOGDXjxxRfx4IMP4tChQxg9erTd+wkKCkJlZSW8vLygUqna7Gu1WhEcHIzKykrodLr7HptsHEfX0pXHIYRAXV1dq5/Fu5NKCHuOj+h+WK1W6PV6WCyWLvePxBEcR9fSXcbhEpfaiaj7YfgQkRQMHyfSarVISUlx+Uv1HEfX0l3GwXM+RCQFj3yISAqGDxFJwfAhIikYPkQkBcOnnS5duoSFCxciNDQUHh4eGDp0KFJSUtDQ0GDTp7U79k+fPm2zra44xeubb76JwYMHw93dHVFRUThz5ozskmykpqZi/Pjx8PLygr+/P+Lj41FaWmrTZ9q0aS1e+6VLl9r0qaioQFxcHDw9PeHv7481a9bg5s2bnTaOTZs2tagxLCxMWX/jxg0YjUb0798fffv2RUJCQoubpGWPwWF239RErTpy5IhYsGCB+Pe//y3Ky8vFhx9+KPz9/cXq1auVPs33oh09elRcvXpVWRoaGpQ+J0+eFL169RJpaWnim2++ERs2bHD4vrOOlpmZKTQajdi7d68oLi4WixYtEt7e3qK6ulpaTXeKjY0V6enpoqioSBQUFIhHH31UhISEiOvXryt9pk6dKhYtWmTz2t8+097NmzfF6NGjRUxMjDh37pzIzs4Wvr6+Ijk5udPGkZKSIsLDw21qvHbtmrJ+6dKlIjg4WOTm5oqzZ8+KiRMniujo6C41BkcxfJwgLS1NhIaGKo9buxH2Th0xxWtHmzBhgjAajcrjW7duiaCgIJGamiqtpnupqakRAEReXp7SNnXqVLF8+fK7Pic7O1uo1WphNpuVtl27dgmdTmczLa8zpaSkiIceeqjVdbW1tcLNzU1kZWUpbSUlJQKAMJlMQoiuMQZH8c8uJ7BYLPDx8WnR/thjj8Hf3x+TJ0/GRx99ZLPufqZ4daaGhgbk5+fb1KRWqxETEyOtJntYLBYAaPH6HzhwAL6+vhg9ejSSk5Px888/K+tMJhMiIiJsZrmMjY2F1WpFcXFx5xQO4Ntvv0VQUBCGDBmCp59+GhUVFQCA/Px8NDY22vwuwsLCEBISovwuusoYHNFt7mrvKsrKyrBjxw5s2bJFaevbty/++te/YtKkSVCr1Xj//fcRHx+PQ4cO4bHHHgNwf1O8OtP333+PW7dutVrThQsXpNR0L01NTVixYgUmTZpkM1PBU089hUGDBiEoKAhff/011q5di9LSUvzrX/8CYN8UvM4WFRWFjIwMjBgxAlevXsVLL72Ehx9+GEVFRTCbzdBoNPD29m5R472mCO7MMTiK4XMX69atwxtvvNFmn5KSEpuTgleuXMGsWbPw+OOPY9GiRUq7r68vVq1apTweP348qqqqsHnzZiV8qP2MRiOKiorwxRdf2LTfPq1KREQEBgwYgBkzZqC8vBxDhw7t7DJbNXv2bOXnMWPGICoqCoMGDcJ7770HDw8PiZU5D8PnLuyd4rVZVVUVpk+fjujoaJvvCrubqKgo5OTkKI/vZ4pXZ/L19UWvXr26VE1tSUpKwuHDh3HixAkMHDiwzb7Ns1+WlZVh6NChCAwMbHEV784peDubt7c3hg8fjrKyMvz6179GQ0MDamtrbY5+bv9ddMUx3AvP+dyFn58fwsLC2lyav6bnypUrmDZtGiIjI5Geng61+t4va0FBAQYMGKA8bp7i9XaOTvHakTQaDSIjI21qampqQm5urrSaWiOEQFJSEj744AMcO3YMoaGh93xO8zfdNr/+BoMBhYWFqKmpUfrk5ORAp9Nh1KhRTqn7Xq5fv47y8nIMGDAAkZGRcHNzs/ldlJaWoqKiQvlddMUx3JPsM96u7vLly2LYsGFixowZ4vLlyzaXSptlZGSIgwcPipKSElFSUiL+8pe/CLVaLfbu3av06YgpXjtaZmam0Gq1IiMjQ3zzzTdi8eLFwtvb2+aKimzPP/+80Ov14vjx4zav/c8//yyEEKKsrEy8/PLL4uzZs+LixYviww8/FEOGDBFTpkxRttF8mXrmzJmioKBAfPrpp8LPz69TL1OvXr1aHD9+XFy8eFGcPHlSxMTECF9fX1FTUyOE+OVSe0hIiDh27Jg4e/asMBgMwmAwdKkxOIrh007p6el3nce2WUZGhhg5cqTw9PQUOp1OTJgwweayabP33ntPDB8+XGg0GhEeHi4++eSTzhxKq3bs2CFCQkKERqMREyZMEKdPn5Zdko27vfbp6elCCCEqKirElClThI+Pj9BqtWLYsGFizZo1Lb5R89KlS2L27NnCw8ND+Pr6itWrV4vGxsZOG8e8efPEgAEDhEajEQ888ICYN2+eKCsrU9b/73//Ey+88ILo16+f8PT0FL/97W9t/gfXFcbgKE6pQURS8JwPEUnB8CEiKRg+RCQFw4eIpGD4EJEUDB8ikoLhQ0RSMHyoS7t9FsixY8c6dV8ZGRnKvlasWOHUfRHDh1zE0aNHW9z71tHmzZuHq1evdql717oz3tVOLqF///7o37+/U/fh4eEBDw8P5YZhci4e+VCnuXbtGgIDA/Haa68pbadOnYJGo7mvo5q9e/ciPDwcWq0WAwYMQFJSkrJOpVJhz549mDNnDjw9PTFy5EiYTCaUlZVh2rRp6NOnD6Kjo1FeXt4hYyPHMXyo0/j5+WHv3r3YtGkTzp49i7q6OsyfPx9JSUmYMWOGQ9vatWsXjEYjFi9ejMLCQnz00UcYNmyYTZ9XXnkFv//971FQUICwsDA89dRTWLJkCZKTk3H27FllOg6SRPKNrdQDvfDCC2L48OHiqaeeEhEREeLGjRt37Xu3yfeDgoLE+vXr7/o8AGLDhg3KY5PJJACIt99+W2l75513hLu7e4vn3mvCeeoYPPKhTrdlyxbcvHkTWVlZOHDgALRarUPPr6mpQVVV1T2PlsaMGaP83DyfcUREhE3bjRs3YLVaHdo/dQyGD3W68vJyVFVVoampCZcuXXL4+fbOaezm5qb8rFKp7trW1NTkcA3Ufgwf6lQNDQ145plnMG/ePLzyyiv4wx/+YDP1pz28vLwwePBgp196J+fipXbqVOvXr4fFYsH27dvRt29fZGdn47nnnsPhw4cd2s6mTZuwdOlS+Pv7Y/bs2airq8PJkyexbNkyJ1VOHY1HPtRpjh8/jm3btmH//v3Q6XRQq9XYv38/Pv/8c+zatcuhbSUmJmLbtm3YuXMnwsPDMWfOHHz77bdOqpycgdOoUpd26dIlhIaG4ty5c06/vaLZtGnTMHbsWGzbtq1T9tdT8ciHXEJ0dDSio6Oduo8DBw6gb9+++Pzzz526H/oFj3yoS7t586ZyRUyr1SI4ONhp+6qrq1O+aM/b2xu+vr5O2xcxfIhIEv7ZRURSMHyISAqGDxFJwfAhIikYPkQkBcOHiKRg+BCRFAwfIpKC4UNEUvwfz+cqhPyNO0MAAAAASUVORK5CYII=",
"text/plain": [
"<Figure size 258.065x259.74 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"center = v + (a1 + a2 + a3)/2\n",
"\n",
"(-box).plot(origin=center, width=(1000, 1000), basis='xy')\n",
"(-box).plot(origin=center, width=(1000, 1000), basis='yz')\n",
"(-box).plot(origin=center, width=(1000, 1000), basis='xz')"
]
},
{
"cell_type": "code",
"execution_count": 18,
"id": "fd707bd6-0541-4e32-9141-854dbcb67ffd",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<Axes: xlabel='x [cm]', ylabel='z [cm]'>"
]
},
"execution_count": 18,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAR8AAAEHCAYAAACX0kK7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy80BEi2AAAACXBIWXMAAA9hAAAPYQGoP6dpAAAez0lEQVR4nO3de1hUdf4H8DcoDBeZQbnMwIaIl0IUw1BxsEyTlZS2dfMxSzMsVpOgQlg3SVdMN2lJy7U1L22CT6tdqPWyZhbiqqWjBokhKiumgeAAbTrjZeUi398f/ThPIyKMMnzBeb+e5zwP53u+55zPF/LdmXMbByGEABFRO3OUXQAR2SeGDxFJwfAhIikYPkQkBcOHiKRg+BCRFAwfIpKC4UNEUnSVXUBn09DQgIqKCnh4eMDBwUF2OUQdihACFy9ehL+/Pxwdb35sw/CxUkVFBQICAmSXQdShlZWV4a677rppH4aPlTw8PAAA67PHw83NSXI1RB3LlSt1iJ20Xfl3cjMMHys1ftRyc3OCmzvDh+hGWnNKgieciUgKhg8RScHwISIpOlX4lJeX46mnnoKXlxdcXV0RGhqKvLw8ZbkQAgsWLICfnx9cXV0RFRWFkydPWmzjp59+wtSpU6FWq+Hp6Ym4uDhcunSpvYdCZPc6TficP38eI0aMgJOTEz7//HMcO3YMy5YtQ/fu3ZU+GRkZWLFiBVavXo2DBw/C3d0d0dHRuHr1qtJn6tSpKCoqQk5ODrZt24a9e/di5syZMoZEZNccOsubDOfOnYt9+/bhq6++uuFyIQT8/f2RkpKCP/zhDwAAk8kErVaLrKwsPPHEEzh+/DhCQkLwzTffYMiQIQCAHTt2YPz48Th79iz8/f1brMNsNkOj0SD7s9/yahfRda5crsOkmC0wmUxQq9U37dtpjny2bt2KIUOGYNKkSfD19cXgwYPx7rvvKstPnz4No9GIqKgopU2j0SAiIgIGgwEAYDAY4OnpqQQPAERFRcHR0REHDx5sv8EQUecJn++//x6rVq1Cv3798MUXXyA+Ph4vvvgi1q9fDwAwGo0AAK1Wa7GeVqtVlhmNRvj6+los79q1K3r06KH0uV5NTQ3MZrPFRES3r9PcZNjQ0IAhQ4ZgyZIlAIDBgwfj6NGjWL16NWJjY2223/T0dLz66qs22z6Rveo0Rz5+fn4ICQmxaOvfvz9KS0sBADqdDgBQWVlp0aeyslJZptPpUFVVZbG8vr4eP/30k9LneqmpqTCZTMpUVlbWJuMhsnedJnxGjBiB4uJii7b//Oc/CAwMBAAEBQVBp9MhNzdXWW42m3Hw4EHo9XoAgF6vx4ULF5Cfn6/02bVrFxoaGhAREXHD/apUKqjVaouJiG5fp/nYNXv2bERGRmLJkiV4/PHHcejQIaxduxZr164F8POzJElJSfjzn/+Mfv36ISgoCH/605/g7++PCRMmAPj5SOnhhx/GjBkzsHr1atTV1SExMRFPPPFEq650EVHb6TThM3ToUGzatAmpqalYtGgRgoKCsHz5ckydOlXp88c//hGXL1/GzJkzceHCBdx///3YsWMHXFxclD4bNmxAYmIixowZA0dHR0ycOBErVqyQMSQiu9Zp7vPpKHifD1Hz7sj7fIjozsLwISIpGD5EJAXDh4ikYPgQkRQMHyKSguFDRFIwfIhICoYPEUnB8CEiKRg+RCQFw4eIpGD4EJEUDB8ikoLhQ0RSMHyISAqGDxFJwfAhIikYPkQkBcOHiKRg+BCRFAwfIpKC4UNEUjB8iEgKhg8RScHwISIpGD5EJAXDh4ikYPgQkRSdNnxef/11ODg4ICkpSWm7evUqEhIS4OXlhW7dumHixImorKy0WK+0tBQxMTFwc3ODr68v5syZg/r6+naunog6Zfh88803WLNmDQYNGmTRPnv2bPzrX/9CdnY29uzZg4qKCjz22GPK8mvXriEmJga1tbXYv38/1q9fj6ysLCxYsKC9h0Bk9zpd+Fy6dAlTp07Fu+++i+7duyvtJpMJ7733Ht5880089NBDCA8PR2ZmJvbv348DBw4AAL788kscO3YM//jHPxAWFoZx48Zh8eLFWLlyJWpra2UNicgudbrwSUhIQExMDKKioiza8/PzUVdXZ9EeHByMnj17wmAwAAAMBgNCQ0Oh1WqVPtHR0TCbzSgqKmqfARARAKCr7AKs8eGHH+Lbb7/FN99802SZ0WiEs7MzPD09Ldq1Wi2MRqPS55fB07i8cdmN1NTUoKamRpk3m823MwQi+n+d5sinrKwML730EjZs2AAXF5d22296ejo0Go0yBQQEtNu+ie5knSZ88vPzUVVVhfvuuw9du3ZF165dsWfPHqxYsQJdu3aFVqtFbW0tLly4YLFeZWUldDodAECn0zW5+tU439jneqmpqTCZTMpUVlbW9oMjskOdJnzGjBmDwsJCFBQUKNOQIUMwdepU5WcnJyfk5uYq6xQXF6O0tBR6vR4AoNfrUVhYiKqqKqVPTk4O1Go1QkJCbrhflUoFtVptMRHR7es053w8PDwwcOBAizZ3d3d4eXkp7XFxcUhOTkaPHj2gVqvxwgsvQK/XY/jw4QCAsWPHIiQkBNOmTUNGRgaMRiPmz5+PhIQEqFSqdh8TkT3rNOHTGm+99RYcHR0xceJE1NTUIDo6Gu+8846yvEuXLti2bRvi4+Oh1+vh7u6O2NhYLFq0SGLVRPbJQQghZBfRmZjNZmg0GmR/9lu4uTvJLoeoQ7lyuQ6TYrbAZDK1eIqi05zzIaI7C8OHiKRg+BCRFAwfIpKC4UNEUjB8iEgKhg8RScHwISIpGD5EJAXDh4ikYPgQkRQMHyKSguFDRFIwfIhICoYPEUnB8CEiKRg+RCQFw4eIpGD4EJEUDB8ikoLhQ0RSMHyISIpWfW9XcnKy1RueP38+evToYfV6RGQfWhU+y5cvh16vh7Ozc6s2+vXXXyMxMZHhQ0TNavU3lm7atAm+vr6t6uvh4XHLBZH1xj+YLbsEkmj7nkmyS7glrTrnk5mZCY1G0+qNrlmzBlqt9paLIqI7X6uOfGJjY63a6JQpU26pGCKyH63+2HUjly5dQkNDg0VbS9/PTEQE3MKl9tOnTyMmJgbu7u7QaDTo3r07unfvDk9PT3Tv3t0WNRLRHcjq8Hnqqadw/vx5rFu3Drm5udi1axd27dqFf//739i1a5ctagQApKenY+jQofDw8ICvry8mTJiA4uJiiz5Xr15FQkICvLy80K1bN0ycOBGVlZUWfUpLSxETEwM3Nzf4+vpizpw5qK+vt1ndRHRjVn/sOnLkCPLz83HPPffYop5m7dmzBwkJCRg6dCjq6+vxyiuvYOzYsTh27Bjc3d0BALNnz8Znn32G7OxsaDQaJCYm4rHHHsO+ffsAANeuXUNMTAx0Oh3279+Pc+fO4emnn4aTkxOWLFnSruMhsndWh8/QoUNRVlbW7uGzY8cOi/msrCz4+voiPz8fI0eOhMlkwnvvvYeNGzfioYceAvDzVbr+/fvjwIEDGD58OL788kscO3YMO3fuhFarRVhYGBYvXoyXX34ZCxcubPV9TER0+6wOn7///e+YNWsWysvLMXDgQDg5OVksHzRoUJsVdzMmkwkAlBsZ8/PzUVdXh6ioKKVPcHAwevbsCYPBgOHDh8NgMCA0NNTiNoDo6GjEx8ejqKgIgwcPbpfaiegWwqe6uhqnTp3CM888o7Q5ODhACAEHBwdcu3atTQu8kYaGBiQlJWHEiBEYOHAgAMBoNMLZ2Rmenp4WfbVaLYxGo9Ln+vuPGucb+1yvpqYGNTU1yrzZbG6rYRDZNavD59lnn8XgwYPxwQcfQKvVwsHBwRZ13VRCQgKOHj2Kr7/+2ub7Sk9Px6uvvmrz/RDZG6vD54cffsDWrVvRt29fW9TTosTERGzbtg179+7FXXfdpbTrdDrU1tbiwoULFkc/lZWV0Ol0Sp9Dhw5ZbK/xalhjn+ulpqZaPFhrNpsREBDQVsMhsltWX2p/6KGHcOTIEVvUclNCCCQmJmLTpk3YtWsXgoKCLJaHh4fDyckJubm5SltxcTFKS0uh1+sBAHq9HoWFhaiqqlL65OTkQK1WIyQk5Ib7ValUUKvVFhMR3T6rj3x+85vfYPbs2SgsLERoaGiTE86PPvpomxX3SwkJCdi4cSO2bNkCDw8P5RyNRqOBq6srNBoN4uLikJycjB49ekCtVuOFF16AXq/H8OHDAQBjx45FSEgIpk2bhoyMDBiNRsyfPx8JCQlQqVQ2qZuIbsxBCCGsWcHRsfmDJVuecG7u3FJmZiamT58O4OebDFNSUvDBBx+gpqYG0dHReOeddyw+Uv3www+Ij4/H7t274e7ujtjYWLz++uvo2rV1OWw2m6HRaJD92W/h5u7U8grtgE+127eO9FT7lct1mBSzBSaTqcVPCVYf+Vz/LFd7aU1Guri4YOXKlVi5cmWzfQIDA7F9+/a2LI2IbgFfo0pEUlgdPi+++CJWrFjRpP1vf/sbkpKS2qImIrIDVofPp59+ihEjRjRpj4yMxCeffNImRRHRnc/q8Pnvf/97w7caqtVq/Pjjj21SFBHd+awOn759+zZ5yBMAPv/8c/Tu3btNiiKiO5/VV7uSk5ORmJiI6upq5enx3NxcLFu2DMuXL2/r+ojoDnVLz3bV1NTgtddew+LFiwEAvXr1wqpVq/D000+3eYFEdGe6pXc4x8fHIz4+HtXV1XB1dUW3bt3aui4iusPd1gvkfXx82qoOIrIzrTrhfN999+H8+fOt3uj999+P8vLyWy6KiO58rTryKSgowJEjR1r99ccFBQUWL+AiIrpeqz92jRkzplXPVwHNPwRKRNSoVeFz+vRpqzf8yxd9ERFdr1XhExgYaOs6iMjO8Kl2IpKC4UNEUjB8iEgKhg8RSWF1+MTGxmLv3r22qIWI7IjV4WMymRAVFYV+/fphyZIlvJOZiG6J1eGzefNmlJeXIz4+Hh999BF69eqFcePG4ZNPPkFdXZ0taiSiO9AtnfPx8fFBcnIyjhw5goMHD6Jv376YNm0a/P39MXv2bJw8ebKt6ySiO8xtPdV+7tw55OTkICcnB126dMH48eNRWFiIkJAQZGRkYPbs2W1VJ91ER/reJqLWsvrIp66uDp9++ikeeeQRBAYGIjs7G0lJSaioqMD69euxc+dOfPzxx1i0aJEt6iWiO4TVRz5+fn5oaGjAk08+iUOHDiEsLKxJn9GjR8PT07MNyiOiO5XV4fPWW29h0qRJcHFxabaPp6fnLT2MSkT2w+rwmTZtmi3qICI7wzuciUgKhg8RSWG34bNy5Ur06tULLi4uiIiIwKFDh2SXRGRX7DJ8PvroIyQnJyMtLQ3ffvst7r33XkRHR6Oqqkp2aUR2wy7D580338SMGTPwzDPPICQkBKtXr4abmxvWrVsnuzQiu2F34VNbW4v8/HxERUUpbY6OjoiKioLBYJBYGZF9ua3HKzqjH3/8EdeuXYNWq7Vo12q1OHHiRJP+NTU1Fl8DZDabbV4jkT2wuyMfa6Wnp0Oj0ShTQECA7JKI7gh2Fz7e3t7o0qULKisrLdorKyuh0+ma9E9NTYXJZFKmsrKy9iqV6I5md+Hj7OyM8PBw5ObmKm0NDQ3Izc2FXq9v0l+lUkGtVltMRHT77O6cDwAkJycjNjYWQ4YMwbBhw7B8+XJcvnwZzzzzjOzSiOyGXYbP5MmTUV1djQULFsBoNCIsLAw7duxochKaiGzHLsMHABITE5GYmCi7DCK7ZXfnfIioY2D4EJEUDB8ikoLhQ0RSMHyISAqGDxFJwfAhIikYPkQkBcOHiKRg+BCRFAwfIpKC4UNEUjB8iEgKhg8RScHwISIpGD5EJAXDh4ikYPgQkRQMHyKSguFDRFIwfIhICoYPEUnB8CEiKRg+RCQFw4eIpGD4EJEUDB8ikoLhQ0RSMHyISIpOET5nzpxBXFwcgoKC4Orqij59+iAtLQ21tbUW/b777js88MADcHFxQUBAADIyMppsKzs7G8HBwXBxcUFoaCi2b9/eXsMgol/oFOFz4sQJNDQ0YM2aNSgqKsJbb72F1atX45VXXlH6mM1mjB07FoGBgcjPz8cbb7yBhQsXYu3atUqf/fv348knn0RcXBwOHz6MCRMmYMKECTh69KiMYRHZNQchhJBdxK144403sGrVKnz//fcAgFWrVmHevHkwGo1wdnYGAMydOxebN2/GiRMnAACTJ0/G5cuXsW3bNmU7w4cPR1hYGFavXt2q/ZrNZmg0GmR/9lu4uTu18aiIOrcrl+swKWYLTCYT1Gr1Tft2iiOfGzGZTOjRo4cybzAYMHLkSCV4ACA6OhrFxcU4f/680icqKspiO9HR0TAYDO1TNBEpOmX4lJSU4O2338Zzzz2ntBmNRmi1Wot+jfNGo/GmfRqX30hNTQ3MZrPFRES3T2r4zJ07Fw4ODjedGj8yNSovL8fDDz+MSZMmYcaMGTavMT09HRqNRpkCAgJsvk8ie9BV5s5TUlIwffr0m/bp3bu38nNFRQVGjx6NyMhIixPJAKDT6VBZWWnR1jiv0+lu2qdx+Y2kpqYiOTlZmTebzQwgojYgNXx8fHzg4+PTqr7l5eUYPXo0wsPDkZmZCUdHy4M2vV6PefPmoa6uDk5OP58IzsnJwT333IPu3bsrfXJzc5GUlKSsl5OTA71e3+x+VSoVVCqVlSMjopZ0inM+5eXlGDVqFHr27ImlS5eiuroaRqPR4lzNlClT4OzsjLi4OBQVFeGjjz7CX//6V4ujlpdeegk7duzAsmXLcOLECSxcuBB5eXlITEyUMSwiuyb1yKe1cnJyUFJSgpKSEtx1110WyxrvFNBoNPjyyy+RkJCA8PBweHt7Y8GCBZg5c6bSNzIyEhs3bsT8+fPxyiuvoF+/fti8eTMGDhzYruMhok58n48svM+HqHl2cZ8PEXVuDB8ikoLhQ0RSMHyISAqGDxFJwfAhIikYPkQkBcOHiKRg+BCRFAwfIpKC4UNEUjB8iEgKhg8RScHwISIpGD5EJAXDh4ikYPgQkRQMHyKSguFDRFIwfIhICoYPEUnB8CEiKTrF93Z1JI3fNHTlSp3kSog6nsZ/F635Ri5+b5eVzp49y+9qJ2pBWVlZky/4vB7Dx0oNDQ2oqKiAh4cHHBwcbtrXbDYjICAAZWVlLX6BWkfGcXQsHXkcQghcvHgR/v7+cHS8+VkdfuyykqOjY4uJfj21Wt3h/iO5FRxHx9JRx6HRaFrVjyeciUgKhg8RScHwsSGVSoW0tDSoVCrZpdwWjqNjuVPGwRPORCQFj3yISAqGDxFJwfAhIikYPkQkBcPnNp05cwZxcXEICgqCq6sr+vTpg7S0NNTW1lr0cXBwaDIdOHDAYlvZ2dkIDg6Gi4sLQkNDsX379vYeThMrV65Er1694OLigoiICBw6dEh2SRbS09MxdOhQeHh4wNfXFxMmTEBxcbFFn1GjRjX53c+aNcuiT2lpKWJiYuDm5gZfX1/MmTMH9fX17TaOhQsXNqkxODhYWX716lUkJCTAy8sL3bp1w8SJE1FZWdmhxmA1Qbfl888/F9OnTxdffPGFOHXqlNiyZYvw9fUVKSkpSp/Tp08LAGLnzp3i3LlzylRbW6v02bdvn+jSpYvIyMgQx44dE/PnzxdOTk6isLBQxrCEEEJ8+OGHwtnZWaxbt04UFRWJGTNmCE9PT1FZWSmtputFR0eLzMxMcfToUVFQUCDGjx8vevbsKS5duqT0efDBB8WMGTMsfvcmk0lZXl9fLwYOHCiioqLE4cOHxfbt24W3t7dITU1tt3GkpaWJAQMGWNRYXV2tLJ81a5YICAgQubm5Ii8vTwwfPlxERkZ2qDFYi+FjAxkZGSIoKEiZbwyfw4cPN7vO448/LmJiYizaIiIixHPPPWerMls0bNgwkZCQoMxfu3ZN+Pv7i/T0dGk1taSqqkoAEHv27FHaHnzwQfHSSy81u8727duFo6OjMBqNStuqVauEWq0WNTU1tixXkZaWJu69994bLrtw4YJwcnIS2dnZStvx48cFAGEwGIQQHWMM1uLHLhswmUzo0aNHk/ZHH30Uvr6+uP/++7F161aLZQaDAVFRURZt0dHRMBgMNq21ObW1tcjPz7eoydHREVFRUdJqag2TyQQATX7/GzZsgLe3NwYOHIjU1FRcuXJFWWYwGBAaGgqtVqu0RUdHw2w2o6ioqH0KB3Dy5En4+/ujd+/emDp1KkpLSwEA+fn5qKurs/hbBAcHo2fPnsrfoqOMwRp8sLSNlZSU4O2338bSpUuVtm7dumHZsmUYMWIEHB0d8emnn2LChAnYvHkzHn30UQCA0Wi0+A8HALRaLYxGY7vW3+jHH3/EtWvXbljTiRMnpNTUkoaGBiQlJWHEiBEYOHCg0j5lyhQEBgbC398f3333HV5++WUUFxfjn//8J4Dmf/eNy9pDREQEsrKycM899+DcuXN49dVX8cADD+Do0aMwGo1wdnaGp6dnkxob6+sIY7AWw6cZc+fOxV/+8peb9jl+/LjFScHy8nI8/PDDmDRpEmbMmKG0e3t7Izk5WZkfOnQoKioq8MYbbyjhQ7cvISEBR48exddff23RPnPmTOXn0NBQ+Pn5YcyYMTh16hT69OnT3mXe0Lhx45SfBw0ahIiICAQGBuLjjz+Gq6urxMpsh+HTjJSUFEyfPv2mfXr37q38XFFRgdGjRyMyMhJr165tcfsRERHIyclR5nU6XZOrF5WVldDpdNYV3ka8vb3RpUuXDlXTzSQmJmLbtm3Yu3dvi688iYiIAPDzUWqfPn2g0+maXMVrHLessXp6euLuu+9GSUkJfv3rX6O2thYXLlywOPr55d+iI46hJTzn0wwfHx8EBwffdHJ2dgbw8xHPqFGjEB4ejszMzBZfogQABQUF8PPzU+b1ej1yc3Mt+uTk5ECv17ftwFrJ2dkZ4eHhFjU1NDQgNzdXWk03IoRAYmIiNm3ahF27diEoKKjFdQoKCgBA+f3r9XoUFhaiqqpK6ZOTkwO1Wo2QkBCb1N2SS5cu4dSpU/Dz80N4eDicnJws/hbFxcUoLS1V/hYdcQwtkn3Gu7M7e/as6Nu3rxgzZow4e/asxaXSRllZWWLjxo3i+PHj4vjx4+K1114Tjo6OYt26dUqfffv2ia5du4qlS5eK48ePi7S0tA5xqV2lUomsrCxx7NgxMXPmTOHp6WlxRUW2+Ph4odFoxO7duy1+91euXBFCCFFSUiIWLVok8vLyxOnTp8WWLVtE7969xciRI5VtNF6mHjt2rCgoKBA7duwQPj4+7XqZOiUlRezevVucPn1a7Nu3T0RFRQlvb29RVVUlhPj5UnvPnj3Frl27RF5entDr9UKv13eoMViL4XObMjMzBYAbTo2ysrJE//79hZubm1Cr1WLYsGEWl00bffzxx+Luu+8Wzs7OYsCAAeKzzz5rz6Hc0Ntvvy169uwpnJ2dxbBhw8SBAwdkl2Shud99ZmamEEKI0tJSMXLkSNGjRw+hUqlE3759xZw5cyzu8xFCiDNnzohx48YJV1dX4e3tLVJSUkRdXV27jWPy5MnCz89PODs7i1/96ldi8uTJoqSkRFn+v//9Tzz//POie/fuws3NTfzud7+z+B9cRxiDtfhKDSKSgud8iEgKhg8RScHwISIpGD5EJAXDh4ikYPgQkRQMHyKSguFDHdov3wIZFhZm031lZWUp+0pKSrLpvojhQ53Ezp07mzz71tYmT56Mc+fOdahn1+5kfKqdOgUvLy94eXnZdB+urq5wdXVVHhgm2+KRD7Wb6upq6HQ6LFmyRGnbv38/nJ2db+moZt26dRgwYABUKhX8/PyQmJioLHNwcMCaNWvwyCOPwM3NDf3794fBYEBJSQlGjRoFd3d3REZG4tSpU20yNrIew4fajY+PD9atW4eFCxciLy8PFy9exLRp05CYmIgxY8ZYta1Vq1YhISEBM2fORGFhIbZu3Yq+ffta9Fm8eDGefvppFBQUIDg4GFOmTMFzzz2H1NRU5OXlKa/jIEkkP9hKduj5558Xd999t5gyZYoIDQ0VV69ebbZvcy/f9/f3F/PmzWt2PQBi/vz5yrzBYBAAxHvvvae0ffDBB8LFxaXJui29cJ7aBo98qN0tXboU9fX1yM7OxoYNG6BSqaxav6qqChUVFS0eLQ0aNEj5ufF9xqGhoRZtV69ehdlstmr/1DYYPtTuTp06hYqKCjQ0NODMmTNWr9/adxo7OTkpPzs4ODTb1tDQYHUNdPsYPtSuamtr8dRTT2Hy5MlYvHgxfv/731u8+rM1PDw80KtXL5tfeifb4qV2alfz5s2DyWTCihUr0K1bN2zfvh3PPvsstm3bZtV2Fi5ciFmzZsHX1xfjxo3DxYsXsW/fPrzwwgs2qpzaGo98qN3s3r0by5cvx/vvvw+1Wg1HR0e8//77+Oqrr7Bq1SqrthUbG4vly5fjnXfewYABA/DII4/g5MmTNqqcbIGvUaUO7cyZMwgKCsLhw4dt/nhFo1GjRiEsLAzLly9vl/3ZKx75UKcQGRmJyMhIm+5jw4YN6NatG7766iub7od+xiMf6tDq6+uVK2IqlQoBAQE229fFixeVL9rz9PSEt7e3zfZFDB8ikoQfu4hICoYPEUnB8CEiKRg+RCQFw4eIpGD4EJEUDB8ikoLhQ0RSMHyISIr/AxXhTfWnzyZ5AAAAAElFTkSuQmCC",
"text/plain": [
"<Figure size 258.065x259.74 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAR8AAAEHCAYAAACX0kK7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy80BEi2AAAACXBIWXMAAA9hAAAPYQGoP6dpAAAdIElEQVR4nO3df1RUZf4H8PegzAACg/yWBH9mKGBuqDhUasmCRtuysuam62KxaO5Qma6bJIn92GXD7JeVWptgu5HGtmn5oxUxcdPxlKQlpJwwXQgcsI7O+GNlEJ7vH325xxHEGYF5GOb9Ouee5T73mXs/D+R779x75xmVEEKAiMjB3GQXQESuieFDRFIwfIhICoYPEUnB8CEiKRg+RCQFw4eIpGD4EJEUfWUX4GxaWlpQV1cHHx8fqFQq2eUQ9ShCCJw7dw5hYWFwc+v43IbhY6e6ujqEh4fLLoOoR6upqcHAgQM77MPwsZOPjw8AYEPRPfDycpdcDVHPcvFiE9JmbFf+nXSE4WOn1rdaXl7u8OrH8CFqjy2XJHjBmYikYPgQkRQMHyKSguFDRFIwfIhICoYPEUnB8CEiKRg+RCQFw4eIpGD4EJEUDB8ikoLhQ0RSMHyISAqGDxFJwfAhIikYPkQkBcOHiKRg+BCRFAwfIpKC4UNEUjB8iEgKhg8RScHwISIpGD5EJAXDh4ikYPgQkRQMHyKSguFDRFIwfIhICqcNn7/+9a9QqVRYuHCh0nbp0iXo9XoEBATA29sbqampqK+vt3pddXU1kpOT4eXlheDgYCxZsgSXL192cPVE5JTh88UXX2DdunUYPXq0Vfvjjz+Ojz/+GEVFRSgtLUVdXR2mT5+ubG9ubkZycjIsFgv279+PDRs2oKCgAMuXL3f0EIhcntOFz/nz5zF79my89dZb6N+/v9JuMpnw9ttv48UXX8Tdd9+N2NhY5OfnY//+/Thw4AAAYOfOnfjmm2/wj3/8A2PGjMG0adPw7LPP4vXXX4fFYpE1JCKX5HTho9frkZycjISEBKv2srIyNDU1WbVHRkYiIiICBoMBAGAwGBATE4OQkBClT1JSEsxmMyoqKto9XmNjI8xms9VCRJ3XV3YB9ti4cSO+/PJLfPHFF222GY1GqNVq+Pn5WbWHhITAaDQqfa4Mntbtrdvak5ubi6effroLqieiKznNmU9NTQ0ee+wxvPvuu/Dw8HDYcbOysmAymZSlpqbGYccm6s2cJnzKysrQ0NCA2267DX379kXfvn1RWlqKV199FX379kVISAgsFgvOnj1r9br6+nqEhoYCAEJDQ9vc/Wpdb+1zNY1GA19fX6uFiDrPacJnypQpOHLkCA4fPqwsY8eOxezZs5Wf3d3dUVJSorymsrIS1dXV0Ol0AACdTocjR46goaFB6VNcXAxfX1+MGjXK4WMicmVOc83Hx8cH0dHRVm39+vVDQECA0p6eno5FixbB398fvr6+eOSRR6DT6TBhwgQAQGJiIkaNGoU5c+YgLy8PRqMR2dnZ0Ov10Gg0Dh8TkStzmvCxxUsvvQQ3NzekpqaisbERSUlJeOONN5Ttffr0wdatW7FgwQLodDr069cPaWlpeOaZZyRWTeSaVEIIIbsIZ2I2m6HValG07Zfw6ucuuxyiHuXihSbMSN4Ck8l03eujTnPNh4h6F4YPEUnB8CEiKRg+RCQFw4eIpGD4EJEUDB8ikoLhQ0RSMHyISAqGDxFJwfAhIikYPkQkBcOHiKRg+BCRFAwfIpKC4UNEUjB8iEiKXjWNKjmPeyYVyS6hXdtLZ8guwWXwzIeIpGD4EJEUDB8ikoLhQ0RSMHyISAqGDxFJwfAhIikYPkQkBcOHiKRg+BCRFAwfIpLCacInNzcX48aNg4+PD4KDg5GSkoLKykqrPpcuXYJer0dAQAC8vb2RmpqK+vp6qz7V1dVITk6Gl5cXgoODsWTJEly+fNmRQyEiOFH4lJaWQq/X48CBAyguLkZTUxMSExNx4cIFpc/jjz+Ojz/+GEVFRSgtLUVdXR2mT5+ubG9ubkZycjIsFgv279+PDRs2oKCgAMuXL5cxJCKXphJCCNlF3IjTp08jODgYpaWlmDhxIkwmE4KCglBYWIhf//rXAIBjx45h5MiRMBgMmDBhAnbs2IF7770XdXV1CAkJAQCsXbsWTzzxBE6fPg21Wn3d45rNZmi1WhRt+yW8+rl36xh7M36qvXe6eKEJM5K3wGQywdfXt8O+TnPmczWTyQQA8Pf3BwCUlZWhqakJCQkJSp/IyEhERETAYDAAAAwGA2JiYpTgAYCkpCSYzWZUVFS0e5zGxkaYzWarhYg6zynDp6WlBQsXLsTtt9+O6OhoAIDRaIRarYafn59V35CQEBiNRqXPlcHTur11W3tyc3Oh1WqVJTw8vItHQ+SanDJ89Ho9ysvLsXHjxm4/VlZWFkwmk7LU1NR0+zGJXIHTzWSYmZmJrVu3Yu/evRg4cKDSHhoaCovFgrNnz1qd/dTX1yM0NFTp8/nnn1vtr/VuWGufq2k0Gmg0mi4eBRE5zZmPEAKZmZn48MMPsXv3bgwZMsRqe2xsLNzd3VFSUqK0VVZWorq6GjqdDgCg0+lw5MgRNDQ0KH2Ki4vh6+uLUaNGOWYgRATAxjOfK29X22rt2rUIDg62+3XXotfrUVhYiC1btsDHx0e5RqPVauHp6QmtVov09HQsWrQI/v7+8PX1xSOPPAKdTocJEyYAABITEzFq1CjMmTMHeXl5MBqNyM7Ohl6v59kNkYPZFD6bN2/G/fffD09PT5t2WlhYiPPnz3dp+KxZswYAMHnyZKv2/Px8zJ07FwDw0ksvwc3NDampqWhsbERSUhLeeOMNpW+fPn2wdetWLFiwADqdDv369UNaWhqeeeaZLquTiGxj03M+bm5uMBqNNoeJj48PvvrqKwwdOrTTBfY0fM6na/A5n96py5/z+fTTT5XnaWyxY8cO3HTTTTb3JyLXY9PbrkmTJtm10zvuuOOGiiEi13HDt9obGhrQ0NCAlpYWq/bRo0d3uigi6v3sDp+ysjKkpaXh6NGjaL1cpFKpIISASqVCc3NzlxdJRL2P3eHz0EMPYcSIEXj77bcREhIClUrVHXURUS9nd/h89913+OCDDzB8+PDuqIeIXITdTzhPmTIFX331VXfUQkQuxO4zn7/97W9IS0tDeXk5oqOj4e5u/azLfffd12XFEVHvZXf4GAwG7Nu3Dzt27GizjRecichWdr/teuSRR/Db3/4Wp06dQktLi9XC4CEiW9kdPj/++CMef/zxNpNyERHZw+7wmT59Oj799NPuqIWIXIjd13xGjBiBrKwsfPbZZ4iJiWlzwfnRRx/tsuKIqPe6obtd3t7eKC0tRWlpqdU2lUrF8CEim9gdPidOnOiOOojIxTjNNKpE1LvYHT6pqal4/vnn27Tn5eVhxgxOxEREtrE7fPbu3Yt77rmnTfu0adOwd+/eLimKiHo/u8Pn/Pnz7X6tsLu7O7/Nk4hsZnf4xMTEYNOmTW3aN27cyK+fISKb2X2366mnnsL06dNx/Phx3H333QCAkpISvPfeeygq6pmTghNRz2N3+PziF7/A5s2b8Ze//AX//Oc/4enpidGjR2PXrl12z/VMRK7rhuZwTk5ORnJyclfXQkQuhM/5EJEUNoWPv78/fvjhB5t3GhERgf/+9783XBQR9X42ve06e/YsduzYAa1Wa9NOf/zxR87tQ0QdsvmaT1paWnfWQUQuxqbwufqLAYmIOosXnIlICpcNn9dffx2DBw+Gh4cH4uLi8Pnnn8suiciluGT4bNq0CYsWLUJOTg6+/PJL3HrrrUhKSkJDQ4Ps0ohchkuGz4svvoiMjAw8+OCDGDVqFNauXQsvLy+sX79edmlELsPlwsdisaCsrAwJCQlKm5ubGxISEmAwGNr0b2xshNlstlqIqPPsDp+7774bTz/9dJv2M2fOKB807cl++OEHNDc3t/nqn5CQEBiNxjb9c3NzodVqlSU8PNxRpRL1anaHz549e/Daa68hJSUFFy5cUNotFkubCeV7g6ysLJhMJmWpqamRXRJRr3BDb7t27doFo9GICRMm4OTJk11cUvcKDAxEnz59UF9fb9VeX1+P0NDQNv01Gg18fX2tFiLqvBsKnwEDBqC0tBQxMTEYN24c9uzZ08VldR+1Wo3Y2FiUlJQobS0tLSgpKYFOp5NYGZFrsXtKDZVKBeCnM4LCwkI899xzmDp1Kp544okuL667LFq0CGlpaRg7dizGjx+Pl19+GRcuXMCDDz4ouzQil2F3+AghrNazs7MxcuRIp/rs18yZM3H69GksX74cRqMRY8aMwSeffMLvnydyoBv60sCgoCCrttTUVERGRuLgwYNdVlh3y8zMRGZmpuwyiFyW3eEzaNCgdtujoqIQFRXV6YKIyDW43EOGRNQzMHyISAqGDxFJwfAhIikYPkQkBcOHiKRg+BCRFAwfIpKC4UNEUjB8iEgKhg8RScHwISIpGD5EJAXDh4ikYPgQkRQMHyKSwu7JxIi6wvbSGbJLIMl45kNEUjB8iEgKhg8RScHwISIpGD5EJAXDh4ikYPgQkRQMHyKSguFDRFIwfIhICoYPEUnhFOFz8uRJpKenY8iQIfD09MSwYcOQk5MDi8Vi1e/rr7/GnXfeCQ8PD4SHhyMvL6/NvoqKihAZGQkPDw/ExMRg+/btjhoGEV3BKcLn2LFjaGlpwbp161BRUYGXXnoJa9euxZNPPqn0MZvNSExMxKBBg1BWVoaVK1dixYoVePPNN5U++/fvxwMPPID09HQcOnQIKSkpSElJQXl5uYxhEbk0lRBCyC7iRqxcuRJr1qzBd999BwBYs2YNli1bBqPRCLVaDQBYunQpNm/ejGPHjgEAZs6ciQsXLmDr1q3KfiZMmIAxY8Zg7dq1Nh3XbDZDq9WiaNsv4dXPvYtHReTcLl5owozkLTCZTPD19e2wr1Oc+bTHZDLB399fWTcYDJg4caISPACQlJSEyspKnDlzRumTkJBgtZ+kpCQYDIZrHqexsRFms9lqIaLOc8rwqaqqwurVqzF//nylzWg0IiQkxKpf67rRaOywT+v29uTm5kKr1SpLeHh4Vw2DyKVJDZ+lS5dCpVJ1uLS+ZWpVW1uLqVOnYsaMGcjIyOj2GrOysmAymZSlpqam249J5AqkzmS4ePFizJ07t8M+Q4cOVX6uq6vDXXfdhfj4eKsLyQAQGhqK+vp6q7bW9dDQ0A77tG5vj0ajgUajue5YiMg+UsMnKCgIQUFBNvWtra3FXXfdhdjYWOTn58PNzfqkTafTYdmyZWhqaoK7+08XgouLi3HLLbegf//+Sp+SkhIsXLhQeV1xcTF0Ol3XDIiIbOYU13xqa2sxefJkRERE4IUXXsDp06dhNBqtrtXMmjULarUa6enpqKiowKZNm/DKK69g0aJFSp/HHnsMn3zyCVatWoVjx45hxYoVOHjwIDIzM2UMi8ilOcUE8sXFxaiqqkJVVRUGDhxota31SQGtVoudO3dCr9cjNjYWgYGBWL58OebNm6f0jY+PR2FhIbKzs/Hkk0/i5ptvxubNmxEdHe3Q8RCREz/nIwuf8yG6Npd4zoeInBvDh4ikYPgQkRQMHyKSguFDRFIwfIhICoYPEUnB8CEiKRg+RCQFw4eIpGD4EJEUDB8ikoLhQ0RSMHyISAqGDxFJwfAhIikYPkQkBcOHiKRg+BCRFAwfIpKC4UNEUjB8iEgKhg8RScHwISIpGD5EJAXDh4ikYPgQkRQMHyKSguFDRFI4Xfg0NjZizJgxUKlUOHz4sNW2r7/+GnfeeSc8PDwQHh6OvLy8Nq8vKipCZGQkPDw8EBMTg+3btzuociK6ktOFz5/+9CeEhYW1aTebzUhMTMSgQYNQVlaGlStXYsWKFXjzzTeVPvv378cDDzyA9PR0HDp0CCkpKUhJSUF5ebkjh0BEcLLw2bFjB3bu3IkXXnihzbZ3330XFosF69evR1RUFH7zm9/g0UcfxYsvvqj0eeWVVzB16lQsWbIEI0eOxLPPPovbbrsNr732miOHQURwovCpr69HRkYG/v73v8PLy6vNdoPBgIkTJ0KtVittSUlJqKysxJkzZ5Q+CQkJVq9LSkqCwWC45nEbGxthNputFiLqPKcIHyEE5s6di4cffhhjx45tt4/RaERISIhVW+u60WjssE/r9vbk5uZCq9UqS3h4eGeGQkT/T2r4LF26FCqVqsPl2LFjWL16Nc6dO4esrCyH15iVlQWTyaQsNTU1Dq+BqDfqK/Pgixcvxty5czvsM3ToUOzevRsGgwEajcZq29ixYzF79mxs2LABoaGhqK+vt9reuh4aGqr8b3t9Wre3R6PRtDkuEXWe1PAJCgpCUFDQdfu9+uqreO6555T1uro6JCUlYdOmTYiLiwMA6HQ6LFu2DE1NTXB3dwcAFBcX45ZbbkH//v2VPiUlJVi4cKGyr+LiYuh0OptrFkIAAC5ebLL5NUSuovXfReu/kw4JJ3TixAkBQBw6dEhpO3v2rAgJCRFz5swR5eXlYuPGjcLLy0usW7dO6bNv3z7Rt29f8cILL4ijR4+KnJwc4e7uLo4cOWLzsWtqagQALly4dLDU1NRc99+S1DOfrqTVarFz507o9XrExsYiMDAQy5cvx7x585Q+8fHxKCwsRHZ2Np588kncfPPN2Lx5M6Kjo20+TlhYGGpqauDj4wOVStWpms1mM8LDw1FTUwNfX99O7UsmjqNnkTkOIQTOnTvX7rN4V1MJYcv5EXUHs9kMrVYLk8nk9P+xcxw9h7OMwylutRNR78PwISIpGD4SaTQa5OTkOP2tfI6jZ3GWcfCaDxFJwTMfIpKC4UNEUjB8iEgKhg8RScHwcYCTJ08iPT0dQ4YMgaenJ4YNG4acnBxYLBarPu19qv/AgQNW++qJ08C+/vrrGDx4MDw8PBAXF4fPP/9cdkmK3NxcjBs3Dj4+PggODkZKSgoqKyut+kyePLnN7/3hhx+26lNdXY3k5GR4eXkhODgYS5YsweXLlx02jhUrVrSpMTIyUtl+6dIl6PV6BAQEwNvbG6mpqW0+RC17DG3Y/KEmumE7duwQc+fOFf/+97/F8ePHxZYtW0RwcLBYvHix0qf182q7du0Sp06dUhaLxaL02bdvn+jTp4/Iy8sT33zzjcjOzrb7s2ldbePGjUKtVov169eLiooKkZGRIfz8/ER9fb20mq6UlJQk8vPzRXl5uTh8+LC45557REREhDh//rzSZ9KkSSIjI8Pq924ymZTtly9fFtHR0SIhIUEcOnRIbN++XQQGBoqsrCyHjSMnJ0dERUVZ1Xj69Gll+8MPPyzCw8NFSUmJOHjwoJgwYYKIj4/vUWO4GsNHkry8PDFkyBBlvb0Py17t/vvvF8nJyVZtcXFxYv78+d1V5nWNHz9e6PV6Zb25uVmEhYWJ3NxcaTV1pKGhQQAQpaWlStukSZPEY489ds3XbN++Xbi5uQmj0ai0rVmzRvj6+orGxsbuLFeRk5Mjbr311na3nT17Vri7u4uioiKl7ejRowKAMBgMQoieMYar8W2XJCaTCf7+/m3a77vvPgQHB+OOO+7ARx99ZLXtRqaB7U4WiwVlZWVWNbm5uSEhIUFaTddjMpkAoM3v/t1330VgYCCio6ORlZWFixcvKtsMBgNiYmKsZsFMSkqC2WxGRUWFYwoH8O233yIsLAxDhw7F7NmzUV1dDQAoKytDU1OT1d8hMjISERERyt+hp4zhSr3mU+3OpKqqCqtXr7aaCN/b2xurVq3C7bffDjc3N3zwwQdISUnB5s2bcd999wG4sWlgu9MPP/yA5ubmdms6duyYlJo60tLSgoULF+L222+3mslg1qxZGDRoEMLCwvD111/jiSeeQGVlJf71r38BsG2K3u4WFxeHgoIC3HLLLTh16hSefvpp3HnnnSgvL4fRaIRarYafn1+bGq83hbAjx3A1hk8nLF26FM8//3yHfY4ePWp1YbC2thZTp07FjBkzkJGRobQHBgZi0aJFyvq4ceNQV1eHlStXKuFDnaPX61FeXo7PPvvMqv3KaVdiYmIwYMAATJkyBcePH8ewYcMcXWa7pk2bpvw8evRoxMXFYdCgQXj//ffh6ekpsbIbx/DpBFungW1VV1eHu+66C/Hx8VbfJ3YtcXFxKC4uVtZvZBrY7hQYGIg+ffr0qJquJTMzE1u3bsXevXsxcODADvu2zo5ZVVWFYcOGITQ0tM0dvKun6HU0Pz8/jBgxAlVVVfj5z38Oi8WCs2fPWp39XPl36Ilj4DWfTggKCkJkZGSHS+tX+dTW1mLy5MmIjY1Ffn4+3Nyu/6s/fPgwBgwYoKy3TgN7JXunge1KarUasbGxVjW1tLSgpKREWk1XE0IgMzMTH374IXbv3o0hQ4Zc9zWt34Tb+rvX6XQ4cuQIGhoalD7FxcXw9fXFqFGjuqXu6zl//jyOHz+OAQMGIDY2Fu7u7lZ/h8rKSlRXVyt/h544Bt7tcoDvv/9eDB8+XEyZMkV8//33VrdLWxUUFIjCwkJx9OhRcfToUfHnP/9ZuLm5ifXr1yt9umIa2K62ceNGodFoREFBgfjmm2/EvHnzhJ+fn9VdFZkWLFggtFqt2LNnj9Xv/eLFi0IIIaqqqsQzzzwjDh48KE6cOCG2bNkihg4dKiZOnKjso/U2dWJiojh8+LD45JNPRFBQkENvUy9evFjs2bNHnDhxQuzbt08kJCSIwMBA0dDQIIT46VZ7RESE2L17tzh48KDQ6XRCp9P1qDFcjeHjAPn5+dec67ZVQUGBGDlypPDy8hK+vr5i/PjxVrdOW73//vtixIgRQq1Wi6ioKLFt2zZHDqVdq1evFhEREUKtVovx48eLAwcOyC5Jca3fe35+vhBCiOrqajFx4kTh7+8vNBqNGD58uFiyZInVcz5CCHHy5Ekxbdo04enpKQIDA8XixYtFU1OTw8Yxc+ZMMWDAAKFWq8VNN90kZs6cKaqqqpTt//vf/8Qf/vAH0b9/f+Hl5SV+9atfWf2fW08Yw9U4pQYRScFrPkQkBcOHiKRg+BCRFAwfIpKC4UNEUjB8iEgKhg8RScHwoR7tyhkex4wZ063HKigoUI61cOHCbj0WMXzISezatavN59q62syZM3Hq1Kke87m03o6faienEBAQgICAgG49hqenJzw9PZUPA1P34pkPOcw777yDgIAANDY2WrWnpKRgzpw5du9v/fr1iIqKgkajwYABA5CZmalsU6lUWLduHe699154eXlh5MiRMBgMqKqqwuTJk9GvXz/Ex8fj+PHjnR4X3RiGDznMjBkz0NzcbDU9bENDA7Zt24aHHnrIrn2tWbMGer0e8+bNw5EjR/DRRx9h+PDhVn2effZZ/O53v8Phw4cRGRmJWbNmYf78+cjKysLBgweV6TZIEmkfaSWXtGDBAjFt2jRlfdWqVWLo0KGipaWl3f7Xmlg/LCxMLFu27JrHASCys7OVdYPBIACIt99+W2l77733hIeHR5vXXm9CeeoaPPMhh8rIyMDOnTtRW1sL4Kc7THPnzoVKpbJ5Hw0NDairq8OUKVM67Dd69Gjl59b5imNiYqzaLl26BLPZbM8QqIvwgjM51M9+9jPceuuteOedd5CYmIiKigps27bNrn3YOmexu7u78nNruLXX1tLSYtfxqWvwzIcc7ve//z0KCgqQn5+PhIQEhIeH2/V6Hx8fDB48uNtvvVP3YviQw82aNQvff/893nrrLbsvNLdasWIFVq1ahVdffRXffvstvvzyS6xevbqLK6XuxPAhh9NqtUhNTYW3tzdSUlJuaB9paWl4+eWX8cYbbyAqKgr33nsvvv32264tlLoVr/mQFLW1tZg9ezY0Gs0N72P+/PmYP39+u9vEVbMDDx48uE3b5MmT27SR4/DMhxzqzJkz+PDDD7Fnzx7o9XqbXxcfH4/4+PhurOynr0z29vbGf/7zn249Dv2EE8iTQw0ePBhnzpzBU089hT/+8Y/X7X/58mWcPHkSAKDRaOy+OG2Pc+fOKV+k5+fnh8DAwG47FjF8iEgSvu0iIikYPkQkBcOHiKRg+BCRFAwfIpKC4UNEUjB8iEgKhg8RScHwISIp/g90Kg9Nk5wuVgAAAABJRU5ErkJggg==",
"text/plain": [
"<Figure size 258.065x259.74 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAR8AAAEHCAYAAACX0kK7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy80BEi2AAAACXBIWXMAAA9hAAAPYQGoP6dpAAAdE0lEQVR4nO3df1SUVf4H8PegzADCDCK/JMFQMxQwT6g4VGrJiobbsrLmprlYrD9qsEyPqyQrlm1saGlaqe0qeFpJY9u0TFsRE0vHY5CWmHLCdCFxwDo644+VQbnfP/bLc5xAnFGGyzDv1znPWebeO/N8LuR7n3meee6ohBACRETtzEN2AUTknhg+RCQFw4eIpGD4EJEUDB8ikoLhQ0RSMHyISAqGDxFJ0VV2Aa6msbERNTU18PPzg0qlkl0OUYcihMDFixcRFhYGD4/Wj20YPg6qqalBeHi47DKIOrTq6mr06tWr1TEMHwf5+fkBADYWPgofH0/J1RB1LFeuNCBt4g7l30lrGD4Oanqr5ePjCZ9uDB+ilthzSoInnIlICoYPEUnB8CEiKRg+RCQFw4eIpGD4EJEUDB8ikoLhQ0RSMHyISAqGDxFJwfAhIikYPkQkBcOHiKRg+BCRFAwfIpKC4UNEUjB8iEgKhg8RScHwISIpGD5EJAXDh4ikYPgQkRQMHyKSguFDRFIwfIhICoYPEUnB8CEiKRg+RCQFw4eIpHDZ8PnrX/8KlUqFOXPmKG1Xr16FwWBAjx494Ovri9TUVNTW1to8r6qqCsnJyfDx8UFwcDDmz5+Pa9eutXP1ROSS4fPVV19h3bp1GDRokE37Cy+8gE8++QSFhYUoKSlBTU0NJkyYoPRfv34dycnJsFqtOHDgADZu3Ij8/HwsXry4vadA5PZcLnwuXbqEKVOm4G9/+xu6d++utJvNZqxfvx5vvPEGHnnkEcTFxSEvLw8HDhzAwYMHAQC7du3Cd999h3/84x8YPHgwxo0bh6VLl+Ltt9+G1WqVNSUit+Ry4WMwGJCcnIzExESb9rKyMjQ0NNi0R0VFISIiAkajEQBgNBoRGxuLkJAQZUxSUhIsFguOHTvW4v7q6+thsVhsNiK6c11lF+CIzZs34+uvv8ZXX33VrM9kMkGtVsPf39+mPSQkBCaTSRlzY/A09Tf1tSQnJwcvvfRSG1RPRDdymSOf6upqPP/889i0aRO8vLzabb+ZmZkwm83KVl1d3W77JurMXCZ8ysrKUFdXh/vvvx9du3ZF165dUVJSglWrVqFr164ICQmB1WrFhQsXbJ5XW1uL0NBQAEBoaGizq19Nj5vG/JJGo4FWq7XZiOjOuUz4jB49GkePHsWRI0eUbciQIZgyZYrys6enJ4qLi5XnVFRUoKqqCnq9HgCg1+tx9OhR1NXVKWOKioqg1WoxcODAdp8TkTtzmXM+fn5+iImJsWnr1q0bevToobSnp6dj7ty5CAgIgFarxezZs6HX6zF8+HAAwJgxYzBw4EBMnToVubm5MJlMyMrKgsFggEajafc5Ebkzlwkfe6xYsQIeHh5ITU1FfX09kpKS8M477yj9Xbp0wfbt2/HMM89Ar9ejW7duSEtLw8svvyyxaiL3pBJCCNlFuBKLxQKdTofCT38Dn26esssh6lCuXG7AxORtMJvNtzw/6jLnfIioc2H4EJEUDB8ikoLhQ0RSMHyISAqGDxFJwfAhIikYPkQkBcOHiKRg+BCRFAwfIpKC4UNEUjB8iEgKhg8RScHwISIpGD5EJAXDh4ik6FTLqLqrR0cWyi6BJNpRMlF2CbeFRz5EJAXDh4ikYPgQkRQMHyKSguFDRFIwfIhICoYPEUnB8CEiKRg+RCQFw4eIpGD4EJEULhM+OTk5GDp0KPz8/BAcHIyUlBRUVFTYjLl69SoMBgN69OgBX19fpKamora21mZMVVUVkpOT4ePjg+DgYMyfPx/Xrl1rz6kQEVwofEpKSmAwGHDw4EEUFRWhoaEBY8aMweXLl5UxL7zwAj755BMUFhaipKQENTU1mDBhgtJ//fp1JCcnw2q14sCBA9i4cSPy8/OxePFiGVMicmsqIYSQXcTtOHfuHIKDg1FSUoIRI0bAbDYjKCgIBQUF+N3vfgcAOHHiBAYMGACj0Yjhw4dj586dGD9+PGpqahASEgIAWLt2LRYsWIBz585BrVbfcr8WiwU6nQ6Fn/4GPt08nTpHe/GudvfWke5qv3K5AROTt8FsNkOr1bY61mWOfH7JbDYDAAICAgAAZWVlaGhoQGJiojImKioKERERMBqNAACj0YjY2FgleAAgKSkJFosFx44da3E/9fX1sFgsNhsR3TmXDJ/GxkbMmTMHDzzwAGJiYgAAJpMJarUa/v7+NmNDQkJgMpmUMTcGT1N/U19LcnJyoNPplC08PLyNZ0PknlwyfAwGA8rLy7F582an7yszMxNms1nZqqurnb5PInfgcisZZmRkYPv27di3bx969eqltIeGhsJqteLChQs2Rz+1tbUIDQ1Vxhw6dMjm9ZquhjWN+SWNRgONRtPGsyAilznyEUIgIyMDH330Efbs2YPIyEib/ri4OHh6eqK4uFhpq6ioQFVVFfR6PQBAr9fj6NGjqKurU8YUFRVBq9Vi4MCB7TMRIgJg55HPjZer7bV27VoEBwc7/LybMRgMKCgowLZt2+Dn56eco9HpdPD29oZOp0N6ejrmzp2LgIAAaLVazJ49G3q9HsOHDwcAjBkzBgMHDsTUqVORm5sLk8mErKwsGAwGHt0QtTO7wmfr1q14/PHH4e3tbdeLFhQU4NKlS20aPmvWrAEAjBo1yqY9Ly8P06ZNAwCsWLECHh4eSE1NRX19PZKSkvDOO+8oY7t06YLt27fjmWeegV6vR7du3ZCWloaXX365zeokIvvY9TkfDw8PmEwmu8PEz88P33zzDfr06XPHBXY0/JwPdTSd+nM+n3/+ufJ5Gnvs3LkTd911l93jicj92PW2a+TIkQ696IMPPnhbxRCR+7jtS+11dXWoq6tDY2OjTfugQYPuuCgi6vwcDp+ysjKkpaXh+PHjaDpdpFKpIISASqXC9evX27xIIup8HA6fp59+Gv3798f69esREhIClUrljLqIqJNzOHx++OEHfPjhh+jXr58z6iEiN+HwJ5xHjx6Nb775xhm1EJEbcfjI5+9//zvS0tJQXl6OmJgYeHraftblsccea7PiiKjzcjh8jEYj9u/fj507dzbr4wlnIrKXw2+7Zs+ejSeffBJnz55FY2OjzcbgISJ7ORw+P//8M1544YVmi3IRETnC4fCZMGECPv/8c2fUQkRuxOFzPv3790dmZia+/PJLxMbGNjvh/Nxzz7VZcUTUed3W1S5fX1+UlJSgpKTEpk+lUjF8iMguDofPqVOnnFEHEbkZl1lGlYg6F4fDJzU1Fa+99lqz9tzcXEyc2HEWNSKijs3h8Nm3bx8effTRZu3jxo3Dvn372qQoIur8HA6fS5cutfi1wp6envw2TyKym8PhExsbiy1btjRr37x5M79+hojs5vDVrj//+c+YMGECTp48iUceeQQAUFxcjPfffx+FhVzInIjs43D4/PrXv8bWrVvx6quv4p///Ce8vb0xaNAg7N692+G1nonIfd3WGs7JyclITk5u61qIyI3wcz5EJIVd4RMQEICffvrJ7heNiIjAf/7zn9suiog6P7vedl24cAE7d+6ETqez60V//vlnru1DRK2y+5xPWlqaM+sgIjdjV/j88osBiYjuFE84E5EUbhs+b7/9Nu6++254eXkhPj4ehw4dkl0SkVtxy/DZsmUL5s6di+zsbHz99de47777kJSUhLq6OtmlEbkNtwyfN954A9OnT8dTTz2FgQMHYu3atfDx8cGGDRtkl0bkNtwufKxWK8rKypCYmKi0eXh4IDExEUajsdn4+vp6WCwWm42I7pzD4fPII4/gpZdeatZ+/vx55UbTjuynn37C9evXm331T0hICEwmU7PxOTk50Ol0yhYeHt5epRJ1ag6Hz969e/HWW28hJSUFly9fVtqtVmuzBeU7g8zMTJjNZmWrrq6WXRJRp3Bbb7t2794Nk8mE4cOH4/Tp021cknMFBgaiS5cuqK2ttWmvra1FaGhos/EajQZardZmI6I7d1vh07NnT5SUlCA2NhZDhw7F3r1727gs51Gr1YiLi0NxcbHS1tjYiOLiYuj1eomVEbkXh5fUUKlUAP53RFBQUIBXXnkFY8eOxYIFC9q8OGeZO3cu0tLSMGTIEAwbNgwrV67E5cuX8dRTT8kujchtOBw+Qgibx1lZWRgwYIBL3fs1adIknDt3DosXL4bJZMLgwYPx2Wef8fvnidrRbX1pYFBQkE1bamoqoqKiUFpa2maFOVtGRgYyMjJkl0HkthwOn969e7fYHh0djejo6DsuiIjcg9t9yJCIOgaGDxFJwfAhIikYPkQkBcOHiKRg+BCRFAwfIpKC4UNEUjB8iEgKhg8RScHwISIpGD5EJAXDh4ikYPgQkRQMHyKSguFDRFI4vJgYdTw7SibKLoHIYTzyISIpGD5EJAXDh4ikYPgQkRQMHyKSguFDRFIwfIhICoYPEUnB8CEiKRg+RCQFw4eIpHCJ8Dl9+jTS09MRGRkJb29v9O3bF9nZ2bBarTbjvv32Wzz00EPw8vJCeHg4cnNzm71WYWEhoqKi4OXlhdjYWOzYsaO9pkFEN3CJ8Dlx4gQaGxuxbt06HDt2DCtWrMDatWvx4osvKmMsFgvGjBmD3r17o6ysDMuWLcOSJUvw7rvvKmMOHDiAJ554Aunp6Th8+DBSUlKQkpKC8vJyGdMicmsqIYSQXcTtWLZsGdasWYMffvgBALBmzRosWrQIJpMJarUaALBw4UJs3boVJ06cAABMmjQJly9fxvbt25XXGT58OAYPHoy1a9fatV+LxQKdTofCT38Dn26ebTwrItd25XIDJiZvg9lshlarbXWsSxz5tMRsNiMgIEB5bDQaMWLECCV4ACApKQkVFRU4f/68MiYxMdHmdZKSkmA0Gm+6n/r6elgsFpuNiO6cS4ZPZWUlVq9ejZkzZyptJpMJISEhNuOaHptMplbHNPW3JCcnBzqdTtnCw8PbahpEbk1q+CxcuBAqlarVrektU5MzZ85g7NixmDhxIqZPn+70GjMzM2E2m5Wturra6fskcgdSVzKcN28epk2b1uqYPn36KD/X1NTg4YcfRkJCgs2JZAAIDQ1FbW2tTVvT49DQ0FbHNPW3RKPRQKPR3HIuROQYqeETFBSEoKAgu8aeOXMGDz/8MOLi4pCXlwcPD9uDNr1ej0WLFqGhoQGenv87EVxUVIR7770X3bt3V8YUFxdjzpw5yvOKioqg1+vbZkJEZDeXOOdz5swZjBo1ChEREVi+fDnOnTsHk8lkc65m8uTJUKvVSE9Px7Fjx7Blyxa8+eabmDt3rjLm+eefx2effYbXX38dJ06cwJIlS1BaWoqMjAwZ0yJyay6xgHxRUREqKytRWVmJXr162fQ1fVJAp9Nh165dMBgMiIuLQ2BgIBYvXowZM2YoYxMSElBQUICsrCy8+OKLuOeee7B161bExMS063yIyIU/5yMLP+dDdHNu8TkfInJtDB8ikoLhQ0RSMHyISAqGDxFJwfAhIikYPkQkBcOHiKRg+BCRFAwfIpKC4UNEUjB8iEgKhg8RScHwISIpGD5EJAXDh4ikYPgQkRQMHyKSguFDRFIwfIhICoYPEUnB8CEiKRg+RCQFw4eIpGD4EJEUDB8ikoLhQ0RSMHyISAqGDxFJ4XLhU19fj8GDB0OlUuHIkSM2fd9++y0eeugheHl5ITw8HLm5uc2eX1hYiKioKHh5eSE2NhY7duxop8qJ6EYuFz5/+tOfEBYW1qzdYrFgzJgx6N27N8rKyrBs2TIsWbIE7777rjLmwIEDeOKJJ5Ceno7Dhw8jJSUFKSkpKC8vb88pEBFcLHx27tyJXbt2Yfny5c36Nm3aBKvVig0bNiA6Ohq///3v8dxzz+GNN95Qxrz55psYO3Ys5s+fjwEDBmDp0qW4//778dZbb7XnNIgILhQ+tbW1mD59Ot577z34+Pg06zcajRgxYgTUarXSlpSUhIqKCpw/f14Zk5iYaPO8pKQkGI3Gm+63vr4eFovFZiOiO+cS4SOEwLRp0zBr1iwMGTKkxTEmkwkhISE2bU2PTSZTq2Oa+luSk5MDnU6nbOHh4XcyFSL6f1LDZ+HChVCpVK1uJ06cwOrVq3Hx4kVkZma2e42ZmZkwm83KVl1d3e41EHVGXWXufN68eZg2bVqrY/r06YM9e/bAaDRCo9HY9A0ZMgRTpkzBxo0bERoaitraWpv+psehoaHK/7Y0pqm/JRqNptl+iejOSQ2foKAgBAUF3XLcqlWr8MorryiPa2pqkJSUhC1btiA+Ph4AoNfrsWjRIjQ0NMDT0xMAUFRUhHvvvRfdu3dXxhQXF2POnDnKaxUVFUGv19tdsxACAHDlSoPdzyFyF03/Lpr+nbRKuKBTp04JAOLw4cNK24ULF0RISIiYOnWqKC8vF5s3bxY+Pj5i3bp1ypj9+/eLrl27iuXLl4vjx4+L7Oxs4enpKY4ePWr3vqurqwUAbty4tbJVV1ff8t+S1COftqTT6bBr1y4YDAbExcUhMDAQixcvxowZM5QxCQkJKCgoQFZWFl588UXcc8892Lp1K2JiYuzeT1hYGKqrq+Hn5weVStXqWIvFgvDwcFRXV0Or1d723GTjPDqWjjwPIQQuXrzY4mfxfkklhD3HR3Q7LBYLdDodzGZzh/uPxBGcR8fSWebhEpfaiajzYfgQkRQMHyfSaDTIzs52+Uv1nEfH0lnmwXM+RCQFj3yISAqGDxFJwfAhIikYPkQkBcPnDp0+fRrp6emIjIyEt7c3+vbti+zsbFitVpsxLd2xf/DgQZvX6ohLvL799tu4++674eXlhfj4eBw6dEh2STZycnIwdOhQ+Pn5ITg4GCkpKaioqLAZM2rUqGa/+1mzZtmMqaqqQnJyMnx8fBAcHIz58+fj2rVr7TaPJUuWNKsxKipK6b969SoMBgN69OgBX19fpKamNrtJWvYcHGb3TU3Uop07d4pp06aJf//73+LkyZNi27ZtIjg4WMybN08Z03Qv2u7du8XZs2eVzWq1KmP2798vunTpInJzc8V3330nsrKyHL7vrK1t3rxZqNVqsWHDBnHs2DExffp04e/vL2pra6XV9EtJSUkiLy9PlJeXiyNHjohHH31UREREiEuXLiljRo4cKaZPn27zuzebzUr/tWvXRExMjEhMTBSHDx8WO3bsEIGBgSIzM7Pd5pGdnS2io6Ntajx37pzSP2vWLBEeHi6Ki4tFaWmpGD58uEhISOhQc3AUw8cJcnNzRWRkpPK4pRthf+nxxx8XycnJNm3x8fFi5syZzirzloYNGyYMBoPy+Pr16yIsLEzk5ORIq+lW6urqBABRUlKitI0cOVI8//zzN33Ojh07hIeHhzCZTErbmjVrhFarFfX19c4sV5GdnS3uu+++FvsuXLggPD09RWFhodJ2/PhxAUAYjUYhRMeYg6P4tssJzGYzAgICmrU/9thjCA4OxoMPPoiPP/7Ypu92lnh1JqvVirKyMpuaPDw8kJiYKK0me5jNZgBo9vvftGkTAgMDERMTg8zMTFy5ckXpMxqNiI2NtVnlMikpCRaLBceOHWufwgF8//33CAsLQ58+fTBlyhRUVVUBAMrKytDQ0GDzt4iKikJERITyt+goc3BEp7mrvaOorKzE6tWrbRa59/X1xeuvv44HHngAHh4e+PDDD5GSkoKtW7fiscceA3B7S7w6008//YTr16+3WNOJEyek1HQrjY2NmDNnDh544AGblQomT56M3r17IywsDN9++y0WLFiAiooK/Otf/wJg3xK8zhYfH4/8/Hzce++9OHv2LF566SU89NBDKC8vh8lkglqthr+/f7Mab7VEcHvOwVEMn5tYuHAhXnvttVbHHD9+3Oak4JkzZzB27FhMnDgR06dPV9oDAwMxd+5c5fHQoUNRU1ODZcuWKeFDd85gMKC8vBxffvmlTfuNy6rExsaiZ8+eGD16NE6ePIm+ffu2d5ktGjdunPLzoEGDEB8fj969e+ODDz6At7e3xMqch+FzE/Yu8dqkpqYGDz/8MBISEmy+K+xm4uPjUVRUpDy+nSVenSkwMBBdunTpUDW1JiMjA9u3b8e+ffvQq1evVsc2rX5ZWVmJvn37IjQ0tNlVvF8uwdve/P390b9/f1RWVuJXv/oVrFYrLly4YHP0c+PfoiPO4VZ4zucmgoKCEBUV1erW9DU9Z86cwahRoxAXF4e8vDx4eNz613rkyBH07NlTedy0xOuNHF3itS2p1WrExcXZ1NTY2Iji4mJpNbVECIGMjAx89NFH2LNnDyIjI2/5nKZvum36/ev1ehw9ehR1dXXKmKKiImi1WgwcONApdd/KpUuXcPLkSfTs2RNxcXHw9PS0+VtUVFSgqqpK+Vt0xDnckuwz3q7uxx9/FP369ROjR48WP/74o82l0ib5+fmioKBAHD9+XBw/flz85S9/ER4eHmLDhg3KmLZY4rWtbd68WWg0GpGfny++++47MWPGDOHv729zRUW2Z555Ruh0OrF3716b3/2VK1eEEEJUVlaKl19+WZSWlopTp06Jbdu2iT59+ogRI0Yor9F0mXrMmDHiyJEj4rPPPhNBQUHtepl63rx5Yu/eveLUqVNi//79IjExUQQGBoq6ujohxP8utUdERIg9e/aI0tJSodfrhV6v71BzcBTD5w7l5eXddB3bJvn5+WLAgAHCx8dHaLVaMWzYMJvLpk0++OAD0b9/f6FWq0V0dLT49NNP23MqLVq9erWIiIgQarVaDBs2TBw8eFB2STZu9rvPy8sTQghRVVUlRowYIQICAoRGoxH9+vUT8+fPt/mcjxBCnD59WowbN054e3uLwMBAMW/ePNHQ0NBu85g0aZLo2bOnUKvV4q677hKTJk0SlZWVSv9///tf8eyzz4ru3bsLHx8f8dvf/tbm/+A6whwcxSU1iEgKnvMhIikYPkQkBcOHiKRg+BCRFAwfIpKC4UNEUjB8iEgKhg91aDeuAjl48GCn7is/P1/Z15w5c5y6L2L4kIvYvXt3s3vf2tqkSZNw9uzZDnXvWmfGu9rJJfTo0QM9evRw6j68vb3h7e2t3DBMzsUjH2o3586dQ2hoKF599VWl7cCBA1Cr1bd1VLNhwwZER0dDo9GgZ8+eyMjIUPpUKhXWrVuH8ePHw8fHBwMGDIDRaERlZSVGjRqFbt26ISEhASdPnmyTuZHjGD7UboKCgrBhwwYsWbIEpaWluHjxIqZOnYqMjAyMHj3aoddas2YNDAYDZsyYgaNHj+Ljjz9Gv379bMYsXboUf/jDH3DkyBFERUVh8uTJmDlzJjIzM1FaWqosx0GSSL6xldzQs88+K/r37y8mT54sYmNjxdWrV2869maL74eFhYlFixbd9HkARFZWlvLYaDQKAGL9+vVK2/vvvy+8vLyaPfdWC85T2+CRD7W75cuX49q1aygsLMSmTZug0Wgcen5dXR1qampuebQ0aNAg5eem9YxjY2Nt2q5evQqLxeLQ/qltMHyo3Z08eRI1NTVobGzE6dOnHX6+vWsae3p6Kj+rVKqbtjU2NjpcA905hg+1K6vViieffBKTJk3C0qVL8cc//tFm6U97+Pn54e6773b6pXdyLl5qp3a1aNEimM1mrFq1Cr6+vtixYweefvppbN++3aHXWbJkCWbNmoXg4GCMGzcOFy9exP79+zF79mwnVU5tjUc+1G727t2LlStX4r333oNWq4WHhwfee+89fPHFF1izZo1Dr5WWloaVK1finXfeQXR0NMaPH4/vv//eSZWTM3AZVerQTp8+jcjISBw+fNjpt1c0GTVqFAYPHoyVK1e2y/7cFY98yCUkJCQgISHBqfvYtGkTfH198cUXXzh1P/Q/PPKhDu3atWvKFTGNRoPw8HCn7evixYvKF+35+/sjMDDQafsihg8RScK3XUQkBcOHiKRg+BCRFAwfIpKC4UNEUjB8iEgKhg8RScHwISIpGD5EJMX/AUcjEmhyB+h5AAAAAElFTkSuQmCC",
"text/plain": [
"<Figure size 258.065x259.74 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"(+box).plot(origin=center, width=(1000, 1000), basis='xy')\n",
"(+box).plot(origin=center, width=(1000, 1000), basis='yz')\n",
"(+box).plot(origin=center, width=(1000, 1000), basis='xz')"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "065dd4f3-27e2-4d61-9349-bd72c6ed6c55",
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.12.4"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment