added shading to amorphous regions

This commit is contained in:
swim committed 2026-11-08 18:16:12 +09:00
1 parent 360d190439
commit 53d30f714e
1 file changed
+44 -31
@@ -160,7 +160,7 @@
{ {
"data": { "data": {
"application/vnd.jupyter.widget-view+json": { "application/vnd.jupyter.widget-view+json": {
"model_id": "0a6acad0be424dec83a5234a9362898b", "model_id": "6398ba34148e4fda86ce67d9eef2de5c",
"version_major": 2, "version_major": 2,
"version_minor": 0 "version_minor": 0
}, },
@@ -181,7 +181,7 @@
{ {
"data": { "data": {
"application/vnd.jupyter.widget-view+json": { "application/vnd.jupyter.widget-view+json": {
"model_id": "497da14575304610bc40120a1e2043ff", "model_id": "97335076ddcc44d4a97df4f9813c5588",
"version_major": 2, "version_major": 2,
"version_minor": 0 "version_minor": 0
}, },
@@ -202,7 +202,7 @@
{ {
"data": { "data": {
"application/vnd.jupyter.widget-view+json": { "application/vnd.jupyter.widget-view+json": {
"model_id": "2ed5546a726842ffb5ed46d0e850da3b", "model_id": "d4511d64692b4e149b429378572bbd28",
"version_major": 2, "version_major": 2,
"version_minor": 0 "version_minor": 0
}, },
@@ -223,7 +223,7 @@
{ {
"data": { "data": {
"application/vnd.jupyter.widget-view+json": { "application/vnd.jupyter.widget-view+json": {
"model_id": "30832de7b6c04805971d9fcaba83108b", "model_id": "3d399fea54864be38895c0f5ef00add9",
"version_major": 2, "version_major": 2,
"version_minor": 0 "version_minor": 0
}, },
@@ -284,7 +284,7 @@
{ {
"data": { "data": {
"application/vnd.jupyter.widget-view+json": { "application/vnd.jupyter.widget-view+json": {
"model_id": "ec7460a8807f4512a2016bdf2ba06e96", "model_id": "7e240cf27c6e403ea776427cb5711009",
"version_major": 2, "version_major": 2,
"version_minor": 0 "version_minor": 0
}, },
@@ -298,7 +298,7 @@
{ {
"data": { "data": {
"application/vnd.jupyter.widget-view+json": { "application/vnd.jupyter.widget-view+json": {
"model_id": "c872d0c0ba644ce6bdf6fb325c5ca2a1", "model_id": "4322bffbb7be4e8680874fce071332f1",
"version_major": 2, "version_major": 2,
"version_minor": 0 "version_minor": 0
}, },
@@ -312,7 +312,7 @@
{ {
"data": { "data": {
"application/vnd.jupyter.widget-view+json": { "application/vnd.jupyter.widget-view+json": {
"model_id": "251920645cdf4d1684abe653d54bab04", "model_id": "bc6040d669a248f48c8f6bac761a9840",
"version_major": 2, "version_major": 2,
"version_minor": 0 "version_minor": 0
}, },
@@ -326,7 +326,7 @@
{ {
"data": { "data": {
"application/vnd.jupyter.widget-view+json": { "application/vnd.jupyter.widget-view+json": {
"model_id": "d6a6853af3f148068e67d1f8be0fd1f6", "model_id": "1bec199574bc403bac57c8e8c72d6900",
"version_major": 2, "version_major": 2,
"version_minor": 0 "version_minor": 0
}, },
@@ -880,15 +880,15 @@
}, },
{ {
"cell_type": "code", "cell_type": "code",
"execution_count": 20, "execution_count": 40,
"id": "76404cd6", "id": "76404cd6",
"metadata": {}, "metadata": {},
"outputs": [ "outputs": [
{ {
"data": { "data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAA94AAAKyCAYAAAAjLAa+AAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAPYQAAD2EBqD+naQAA9ntJREFUeJzs3Xt8jvUfx/HXvbHNsDmMOYy12sqpnBMikVMOSYUk1AodCSVU6CT90klOYw7JKQklqaURoYNDJaqJxV3mmM1xY/f1++ObMRs2du+6t72fj8f92H1f93Xf92eXuT/X53t9Dw7LsixERERERERExC287A5AREREREREJD9T4S0iIiIiIiLiRiq8RURERERERNxIhbeIiIiIiIiIG6nwFhEREREREXEjFd4iIiIiIiIibqTCW0RERERERMSNVHiLiIiIiIiIuJEKbxERERERERE3UuEtkgu+/vprHnzwQapUqULRokWpWLEid9xxBxs2bMjS63v37k2xYsUu+Pz+/fvx8fGhW7duF9wnKSkJf39/OnbsmO34s+Oqq66id+/eaY/j4+NxOBzMmDEjbduMGTNwOBzEx8fn2OeOHDkSh8Nxyf0sy2LevHk0adKEsmXL4ufnR0hICK1bt2bq1Knp9nU4HIwcOTLHYhQREc+3cuVKHA5Hprf169df8vXK2ZemnC0FUSG7AxApCCZOnMjBgwfp378/1apVY//+/YwdO5abbrqJL774gubNm1/R+5cpU4aOHTuyePFi/v33X0qWLJlhn3nz5nHixAkiIyOv6LMuZdGiRQQEBLj1M67E0KFDGTNmDA8//DBPP/00xYsX56+//uLrr79myZIlPPTQQ2n7rlu3jpCQEBujFRERu7z66qvceuut6bbVqFHjit9XOTvrlLMlP1HhLZILxo8fT9myZdNta9OmDeHh4bz66qtXXHgDREZGsnDhQmbPns3jjz+e4flp06YRHBxMu3btrvizLqZ27dpuff8rceLECd5++2169uxJVFRUuud69+6Ny+VKt+2mm27KzfBERMSDREREuC0PKGdfmnK25Dfqai6SC84vugGKFStGtWrV2L1792W957fffktQUBDt27fn2LFjtG7dmpCQEKZPn55h323btvHdd9/Rs2dPChW6cHvbme51c+bMYciQIZQvX55ixYrRoUMH9u7dy5EjR+jTpw9BQUEEBQXxwAMPcPTo0XTvcX63tez46quvaNGiBQEBAfj7+9O4cWNWrFiRYb/PPvuMWrVq4evrS1hYGG+88UaW3v/YsWMkJydTvnz5TJ/38kr/lahuayIikhOUs5WzRVR4i9gkMTGRjRs3Ur169Wy/9sMPP6RFixZ06dKFJUuWULRoUby8vOjduzcbN27kp59+Srf/mcT+4IMPZun9hw0bxr59+5gxYwZjx45l5cqV3Hvvvdx1110EBgYyd+5cnnnmGWbNmsWwYcOyHX9mPvjgA1q1akVAQAAzZ87kww8/pFSpUrRu3TpdIl+xYgV33HEHxYsXZ968efzvf//jww8/zPTk5XxBQUGEh4czYcIE3nzzTX777Tcsy8qR+EVEJH957LHHKFSoEAEBAbRu3Zo1a9Zc1vsoZytniwBgiYgt7rvvPqtQoULWjz/+eMl9e/XqZRUtWtSyLMt67bXXLG9vb2vMmDEZ9tuxY4flcDisJ598Mm3bqVOnrHLlylmNGze+5OfExsZagNWhQ4d02wcMGGAB6d7XsiyrU6dOVqlSpdJtCw0NtXr16pX2eOfOnRZgTZ8+PW3b9OnTLcDauXOnZVmWdezYMatUqVIZPjc1NdWqWbOmdeONN6Zta9CggVWhQgXrxIkTaduSkpKsUqVKWVn5Svv++++typUrW4AFWMWLF7fat29vvf/++5bL5Uq3L2CNGDHiku8pIiL5x8aNG63+/ftbixYtsr755htr2rRpVtWqVS1vb29r+fLll3y9crZytkhmdMVbxAbPP/88s2fP5q233qJu3bpZeo1lWfTt25cRI0YwZ84cnnnmmQz7hIWFceuttzJ79mxSUlIA+Pzzz0lISMhyyzlA+/bt0z2uWrUqQIaxZlWrVuXQoUMZuq5l19q1azl06BC9evXi9OnTaTeXy0WbNm344YcfOHbsGMeOHeOHH36gc+fO+Pn5pb2+ePHidOjQIUufVb9+fbZv387y5csZNmwYDRs2ZMWKFfTs2ZOOHTuqNV1EpICrXbs2b7/9Np06daJJkyY88MADrF27lvLly2eaezOjnK2cLXI+Fd4iuWzUqFG8/PLLvPLKK5lOqHIhKSkpzJ8/n+rVq9O2bdsL7hcZGcnBgwf55JNPANNlrVixYnTp0iXLn1WqVKl0j318fC66/eTJk1l+78zs3bsXgLvvvpvChQunu40ZMwbLsjh06BD//vsvLpeLcuXKZXiPzLZdSOHChWndujWvvPIKX3zxBbt376ZZs2YsXbqUzz///Ip+FxERyX9KlChB+/bt+fnnnzlx4sQl91fOVs4WOZ8Kb5FcNGrUKEaOHMnIkSOzPc7K19eX2NhYdu/ezW233ca///6b6X6dO3emZMmSTJs2jf3797N06VK6du160TVF7RYUFATAuHHj+OGHHzK9BQcHU7JkSRwOBwkJCRneI7NtWVW6dGkGDBgAwJYtWy77fUREJP86c3U1K+tPK2crZ4ucT4W3SC556aWXGDlyJM899xwjRoy4rPeoXbs2q1atwul00qxZM/bt25dhHz8/P7p3786XX37JmDFjOHXqVLa6rNmhcePGlChRgq1bt1KvXr1Mbz4+PhQtWpQbb7yRjz/+OF2L/ZEjR/j0008v+TmnTp3i4MGDmT63bds2ACpUqJAzv5SIiOQb//77L0uXLqVWrVrpuk1fjHK2crbIubSOt0guGDt2LC+88AJt2rShXbt2rF+/Pt3z2Vl7smrVqqxevZrbbruNpk2b8tVXXxESEpJun8jISMaPH8+bb75JlSpVaNSoUY78Hu5SrFgxxo0bR69evTh06BB33303ZcuWZf/+/fz000/s37+fiRMnAqYBo02bNrRs2ZJBgwaRmprKmDFjKFq0KIcOHbro5yQmJnLVVVdxzz33cNttt1GpUiWOHj3KypUreeedd6hatSqdO3fOjV9ZREQ8VPfu3alcuTL16tUjKCiIuLg4xo4dy969e5kxY0a23ks5Wzlb5AwV3iK54EzL7vLly1m+fHmG57M7OcjVV1+dlsibNGnCihUruPrqq9Oer127NrVr12bTpk0e33J+Ro8ePahcuTKvv/46ffv25ciRI5QtW5ZatWqlW2O0ZcuWLF68mOeee46uXbtSrlw5Hn30UU6cOMGoUaMu+hkBAQGMGjWKFStWMGzYMPbu3YvD4SAsLIwBAwYwZMgQ/P393fybioiIJ7vhhhuYP38+kyZN4ujRo5QqVYqbb76ZWbNmUb9+/Wy/n3K2crYIgMPSdIAiIiIiIiIibqMx3iIiIiIiIiJupMJbRERERERExI1UeIuIiIiIiIi4ke2F94QJEwgLC8PPz4+6deuyevXqi+6fnJzM8OHDCQ0NxdfXl2uuuYZp06blUrQiIiIFl3K2iIjI5bF1VvP58+czYMAAJkyYQOPGjZk8eTJt27Zl69atVK5cOdPXdOnShb179xIdHU14eDj79u3j9OnTuRy5iIhIwaKcLSIicvlsndW8QYMG1KlTJ22tPzDrHXbq1InRo0dn2H/58uV069aNHTt2UKpUqdwMVUREpEBTzhYREbl8tl3xTklJYcOGDTz77LPptrdq1Yq1a9dm+ppPPvmEevXq8frrrzNr1iyKFi1Kx44deemllyhSpEimr0lOTiY5OTntscvl4tChQ5QuXRqHw5Fzv5CIiEgOsSyLI0eOUKFCBby8bB8VppwtIiKSiezka9sK7wMHDpCamkpwcHC67cHBwSQkJGT6mh07drBmzRr8/PxYtGgRBw4c4NFHH+XQoUMXHDM2evRoRo0alePxi4iIuNvu3bsJCQmxOwzlbBERkYvISr62dYw3kKEF27KsC7Zqu1wuHA4Hs2fPJjAwEIA333yTu+++m/Hjx2fagj506FAGDhyY9jgxMZHKlSuze/duAgICcvA3ERERyRlJSUlUqlSJ4sWL2x1KOsrZIiIiZ2UnX9tWeAcFBeHt7Z2hpXzfvn0ZWtTPKF++PBUrVkxL4GDGl1mWhdPpJCIiIsNrfH198fX1zbA9ICBASVxERDyap3SvVs4WERG5sKzka9sGjvn4+FC3bl1iYmLSbY+JiaFRo0aZvqZx48b8888/HD16NG3bH3/8gZeXl0d0xRMREcmPlLNFRESujK0ztgwcOJCpU6cybdo0tm3bxlNPPcWuXbvo168fYLqc9ezZM23/7t27U7p0aR544AG2bt3KN998w9NPP82DDz54wYlaRERE5MopZ4uIiFw+W8d4d+3alYMHD/Liiy+yZ88eatSowbJlywgNDQVgz5497Nq1K23/YsWKERMTwxNPPEG9evUoXbo0Xbp04eWXX7brVxARESkQlLNFREQun63reNshKSmJwMBAEhMTNV5MREQ8knKVoeMgIiKeLDt5yv7FQUVERERERETyMRXeIiIiIiIiIm6kwltERERERETEjVR4i4iIiIiIiLiRCm8RERERERERN1LhLSIiIiIiIuJGKrxFRERERERE3EiFt4iIiIiIiIgbqfAWERERERERcSMV3iIiIiIiIiJupMJbRERERERExI1UeIuIiIiIiIi4kQpvERERERERETdS4S0iIiIiIiLiRiq8RURERERERNxIhbeIiIiIiIiIG6nwFhEREREREXEjFd4iIiIiIiIibqTCW0RERERERMSNVHiLiIiIiIiIuJEKbxERERERERE3UuEtIiIiIiIi4kYqvEVERERERETcSIW3iIiIiIiIiBup8BYRERERERFxIxXeIiIiIiIiIm6kwltERERERETEjVR4i4iLine truncated "image/png": "iVBORw0KGgoAAAANSUhEUgAACBUAAAW+CAYAAAAyN4pZAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAuIwAALiMBeKU/dgABAABJREFUeJzs3XlcVFUfx/EvIAoK7vu+m+WCSAUYWmnmmmmSe1rS8lBqPT6ZbVqZlu25VmpZmRtpLqhZpikKmoCkuW+57yuoqMA8f9wghxlghhlWP+/X677kHs6550eZec/8zu+4mEwmkwAAAAAAAAAAAAAAANJxzesAAAAAAAAAAAAAAABA/kRSAQAAAAAAAAAAAAAAsIqkAgAAAAAAAAAAAAAAYBVJBQAAAAAAAAAAAAAAwCqSCgAAAAAAAAAAAAAAgFUkFQAAAAAAAAAAAAAAAKtIKgAAAAAAAAAAAAAAAFaRVAAAAAAAAAAAAAAAAKwiqQAAAAAAAAAAAAAAAFhFUgEAAAAAAAAAAAAAALCKpAIAAAAAAAAAAAAAAGAVSQUAAAAAAAAAAAAAAMAqkgoAAAAAAAAAAAAAAIBVJBUAAAAAAAAAAAAAAACrSCoAAAAAAAAAAAAAAABWkVQAAAAAAAAAAAAAAACsIqkAAAAAAAAAAAAAAABYRVIBAAAAAAAAAAAAAACwiqQCAAAAAAAAAAAAAABgFUkFAAAAAAAAAAAAAADAKpIKAAAAAAAAAAAAAACAVSQVAAAAAAAAAAAAAAAAq0gqAAAAAAAAAAAAAAAAVpFUAAAAAAAAAAAAAAAArCKpAAAAAAAAAAAAAAAAWEVSAQAAAAAAAAAAAAAAsIqkAgAAAAAAAAAAAAAAYBVJBQAAAAAAAAAAAAAAwKoieR0AAAC5JSkpSbGxsdqxY4fOnDmjmzdvqmTJkqpfv75atmypChUq5HWIKGQOHDigF154QevXr1e1atU0evRo9e7dO6/DgpOlpKRo+/bt2rlzp44fP64rV66oSJEiKlWqlOrWravmzZurUqVK2X5+bGyshg0bpri4ODVo0EAffPCB2rVr58SfAAAAAADgLKdPn9awYcP0888/q0yZMho2bJiGDRuW12EhB+zZs0fbt2/XkSNHlJCQIFdXV5UsWVK1a9dWkyZNVLNmzWw/mzUlAPmOCQCAQm7Xrl2mp59+2lS6dGmTJKuXq6urqU2bNqaffvrJlJKSkmuxjR492iyOUqVKOX2Ozz//3OLn/fDDD50+j8lkMrVp08Zirq1bt+bIXPldQkKCqV69ehb/PFasWJHhmG+++cai/4ULF7KcKzd+H+U38fHxFv+svvnmm1yNYcOGDaZBgwaZypUrl+GfLalX06ZNTePHjzedPXvWrjmOHz9uKlu2rNmzihYtavrzzz9z6KcCAAAACq5ff/01y7+b23rl5HtVdt/97LF48WKLOZ5//nmnzpFq4MCBFnMtWbIkR+bK71JSUkyBgYEW/zy++OKLDMesWbPGov+WLVuynCs3fh/lR+nfwUePHp2r82/dutX0wgsvmKpUqZLlnyP169c3vfHGG6YjR47YNUd21pQAIKdx/AEAoNC6efOmXnnlFTVp0kTTpk3TxYsXM+ybkpKitWvXqnv37nrwwQd16NCh3As0h/Xu3VtFipgXJ5o9e7bT5zlx4oQiIiLM2po2baqmTZs6fa6CYP369dq/f79F+7fffpsH0cCZ9u/frw4dOqhVq1aaOXOmzp07l+WYbdu26ZVXXlHdunX1wQcfKDk52aa5lixZovPnz5u13bhxQ3PmzMlW7AAAAEBhdvjw4bwOId/o2LGjypcvb9YWFhampKQkp85z48YNLVq0yKytXLly6tChg1PnKSj27NmjyMhIi/aZM2fmfjBwqtOnT6tv375q3ry5Jk2apBMnTmQ5Zt++fXr33XdVv359jRgxQomJiTbNxZoSgPyIpAIAQKGUkJCghx56SB988IHdL8y///677r77bsXExORQdLmrYsWKat++vVnbli1btHv3bqfOM3/+fKWkpJi19e/f36lzFAYuLi55HQIcsGjRIvn4+GjlypXZGn/58mW98sorevDBB3XhwoVsx8HvIwAAAMASSQX/cnd31+OPP27Wdvr0af32229OnWfFihW6dOmSWdvjjz8ud3d3p85T0PEOV7BFRUWpadOmmjNnjkwmk93jr1+/rg8//FB+fn4O/TnF7yMAealI1l0AAChYUlJS1L17d61du9bie2XLllXbtm1Vo0YNubu769ixY1q7dq2OHDli1u/MmTPq2LGj/vjjD9WuXTuXIs85AwYM0PLly83aZs+erbfffttpc8ybN8/s3sXFRX369HHa8wuaoKAg1a9fX/v27TNrHzRoUN4EBIf99NNPCg4OtlploHTp0mrTpo2qVaum8uXL69q1azp27Jg2b96svXv3WvRft26d2rZtq4iICJUoUSLDObt166bXXnvNrFpB0aJF1bdvX+f8UAAAAEAhkv7DOhcXF7m6Zm9fXfqKfwXRgAEDNGXKFLO22bNn6+GHH3baHOnXAiSpX79+Tnt+QdOwYUO1atVKGzZsMGt/8skn8ygiOGrjxo1q166drl69avG94sWLq02bNqpZs6YqVqyoGzdu6OTJk4qNjdW2bdss+m/fvl1BQUHatGmTKleunOGcrCkByI8K/t+MAABI54MPPtCqVavM2kqUKKFPP/1UTzzxhIoVK2b2PZPJpJUrV+r555/XgQMH0trPnDmjgQMHWk1OKGi6deumkiVL6vLly2ltc+bMcVpSweHDh7Vx40aztjZt2qhGjRpOeX5BVLx4cf38888aMmSI1q9fr2rVqmnUqFEWVSNQMOzdu1cDBgywSCho1qyZ3nvvPbVv3z7DRcetW7dq3Lhxmj9/vtmOhi1btig0NDTT8oWVK1fWr7/+qmHDhmnLli1q2LChPvjgAzVp0sQ5PxgAAABQiKRPKujbt69mzZqVR9HkPX9/fzVo0MAs0fmnn37Sl19+KQ8PD4eff+3aNS1dutSsrU6dOmrVqpXDzy6oXFxctHDhQg0dOlQrV65U6dKl9eKLL+rpp5/O69CQDefPn1ePHj0sEgrq1KmjcePG6dFHH83wv6X9+/fr448/1rRp08yqqB4+fFi9e/fWmjVrMqw8wJoSgPyI4w8AAIXK2bNnNXbsWLO2UqVKaePGjXr66actEgok44WvQ4cO+vPPP3X33XebfW/dunVavHhxjsacGzw9PdWjRw+ztr179yo6Otopz583b55F+bfbeWdCqnr16mn58uW6fPmydu7ceVtXbijoXnrpJV25csWsbciQIYqNjVWnTp0y3cXUrFkzzZ07Vz/++KM8PT3Nvvfdd98pKioq07l9fX0VERGhhIQExcbGql27dtn/QQAAAIBCLH1Swe2c6J4q/bt5fHy8wsPDnfLs8PBwJSQkmLVRVc04hnLu3Lm6cOGCDh48qGHDhuV1SMimUaNG6cSJE2ZtPXr00I4dO9S7d+9Mk3Pq1aunKVOmaPXq1SpfvrzZ99auXWu1ykf68awpAchPSCoAABQqM2fOtHihnTx5sk27er28vLR48WJ5eXmZtU+dOtWpMeaVAQMGWLTNnj3bKc+eO3eu2X2xYsXUs2dPpzwbyGt79+7VsmXLzNp69uypCRMmyM3Nzebn9OjRQ19++aVFe/pypAAAAACy5+jRo2b31atXz6NI8o/+/ftbtOXUWkBG8wEF0eXLl/XNN9+YtQUEBGju3Ll2VfoICgrSjz/+aHEUy+TJk50SJwDkFpIKAACFypIlS8zua9asaVeWfJUqVSzOJ/vtt9907do1Z4SXp+6//36LBZV58+YpJSXFoefu27dPsbGxZm2dO3dW6dKlHXoukF+k/3PFzc1N77//fraeNWDAAPn6+pq1pT+uBQAAAID9Tp8+bfHuTqUCY7dzYGCgWdvy5ct16dIlh54bHx+v5cuXm7X5+vrqjjvucOi5QH7x66+/Whx78N5778nd3d3uZ7Vp00bdunUza4uKirKoiAgA+RlJBQCAQsNkMmnLli1mbT179szwfLKMdOzY0ew+KSlJO3bscDi+vObq6mqRYHH8+HGtXbvWoeeyMwGFXfo/Vxo2bKh69epl+3ndu3c3uz958qTFQgUAAAAA+6Q/+kAiqSBV+nf069eva+HChQ49c/HixUpMTMx0HqAgS78WULx4cbVu3Trbz0t/LGlycrIOHTqU7ecBQG4jqQAAUGhcuHDB4uiD+vXr2/2cWrVqWbSdOXMm23HlpCVLlqhIkSJycXFJu0qWLGnx4pMqJ45ASH8GXOnSpdW5c2eHnpmRmTNnmv2s6Y+qSEpK0g8//KBu3bqpZs2aKlasmMqXL68HHnhAU6ZM0fXr1zN8dlJSkhYvXqxBgwapcePGKlOmjIoWLarq1aura9eu+vrrry0WTLLy4osvmsVbu3bt7PzYue7s2bOaNGmSevTooXr16qlkyZIqVqyYqlWrpoCAAI0cOVK///67U+eMi4vTyy+/rMDAQFWoUEHu7u4qXry46tevr969ezvtzM/sOHv2rNm9tT8j7BEaGqrNmzebXUWLFs2w/6OPPmr2++j+++93aH4AAACgMLqdkgr++OMPeXl5mb0nFC1aVD///LPV/r169bJ453D2WoCbm1uOnfn++++/m/2sLi4uFu9pS5YsUZ8+fVS3bl15enqqbNmyCggI0Pjx4xUfH5/p81etWqXQ0FA1bdpU5cuXl7u7uypXrqyHHnpIEyZMsLuqw2effWYRb0EQHx+vb775Rn369NEdd9yh0qVLq2jRoqpSpYpatmypYcOGadmyZQ5XvLzV3r17NXr0aD3wwAOqXLmyihUrJk9PT9WqVUuPPvqoZs+ereTkZKfNZ4/0v8dq1qzp0L/L7t27W6wFZLZOVFDXlAAUXkXyOgAAAJzF3d1dY8aMMWvLzodv1j54LlIk//0vMzIyUr1Line truncated
"text/plain": [ "text/plain": [
"<Figure size 1000x700 with 4 Axes>" "<Figure size 2100x1500 with 4 Axes>"
] ]
}, },
"metadata": {}, "metadata": {},
@@ -905,7 +905,7 @@
"ys = np.array(ys)\n", "ys = np.array(ys)\n",
"Y = ys.sum(axis=0)\n", "Y = ys.sum(axis=0)\n",
"\n", "\n",
"fig, axs = plt.subplots(2, 2, figsize=(10, 7))\n", "fig, axs = plt.subplots(2, 2, figsize=(7, 5), dpi=300)\n",
"\n", "\n",
"keys = ['2V', '5V', '8V', '30V']\n", "keys = ['2V', '5V', '8V', '30V']\n",
"\n", "\n",
@@ -925,6 +925,9 @@
" ax.set_xlabel(\"Depth (Å)\")\n", " ax.set_xlabel(\"Depth (Å)\")\n",
" ax.set_ylabel(\"Phase (rad)\")\n", " ax.set_ylabel(\"Phase (rad)\")\n",
" ax.set_ylim(-0.05, 0.6)\n", " ax.set_ylim(-0.05, 0.6)\n",
" shade=\"#ddd\"\n",
" ax.axvspan(0, popt[2]*dz, color=shade)\n",
" ax.axvspan(popt[3]*dz, xlim[-1], color=shade)\n",
"\n", "\n",
"fig.tight_layout()\n", "fig.tight_layout()\n",
"plt.show()\n" "plt.show()\n"
@@ -932,13 +935,23 @@
}, },
{ {
"cell_type": "code", "cell_type": "code",
"execution_count": 32, "execution_count": 21,
"id": "669830e5",
"metadata": {},
"outputs": [],
"source": [
"plt.rcParams['font.family'] = \"Inter Variable ss02\""
]
},
{
"cell_type": "code",
"execution_count": 22,
"id": "fcd12e61", "id": "fcd12e61",
"metadata": {}, "metadata": {},
"outputs": [ "outputs": [
{ {
"data": { "data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAABaYAAAVeCAYAAACKJuzmAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAuIwAALiMBeKU/dgABAABJREFUeJzs3XVcVNn7B/DPUAoCEgaoiIq5drdgrN2xNnbXWht2rq27dq2N3a6JjYliIFgooIBggAgokvf3Bz/4OtwBpphLfN6v17x275l7znkuDDg8c+5zZIIgCCAiIiIiIiIiIiIi0hE9qQMgIiIiIiIiIiIiotyFiWkiIiIiIiIiIiIi0ikmpomIiIiIiIiIiIhIp5iYJiIiIiIiIiIiIiKdYmKaiIiIiIiIiIiIiHSKiWkiIiIiIiIiIiIi0ikmpomIiIiIiIiIiIhIp5iYJiIiIiIiIiIiIiKdYmKaiIiIiIiIiIiIiHSKiWkiIiIiIiIiIiIi0ikmpomIiIiIiIiIiIhIp5iYJiIiIiIiIiIiIiKdYmKaiIiIiIiIiIiIiHSKiWkiIiIiIiIiIiIi0ikmpomIiIiIiIiIiIhIp5iYJiIiIiIiIiIiIiKdYmKaiIiIiIiIiIiIiHSKiWkiIiIiIiIiIiIi0ikmpomIiIiIiIiIiIhIp5iYJiIiIiIiIiIiIiKdYmKaiIiIiIiIiIiIiHSKiWkiIiIiIiIiIiIi0ikmpomIiIiIiIiIiIhIp5iYJiIiIiIiIiIiIiKdYmKaiIiIiIiIiIiIiHSKiWkiIiIiIiIiIiIi0ikDqQMgIt0JDw/HtWvXUo7t7OyQJ08eCSMiIiIiIiIiIiKpxMTEICAgIOXY0dERFhYWOpmbiWmiXOTatWvo3Lmz1GEQEREREREREVEWdPz4cXTq1Eknc7GUBxERERERERERERHpFBPTRERERERERERERKRTLOVBlIvY2dnJHR8/fhylS5eWKBoiIiIiIiIiIpLSq1ev5Mq+ps4dZSYmpolykdQbHZYuXRoVK1aUKBoiIiIiIiIiIspKUueOMhNLeRARERERERERERGRTjExTUREREREREREREQ6xcQ0EREREREREREREekUE9NEREREREREREREpFNMTBMRERERERERERGRTjExTUREREREREREREQ6xcQ0EREREREREREREekUE9NEREREREREREREpFNMTBMRERERERERERGRTjExTUREREREREREREQ6xcQ0EREREREREREREekUE9NEREREREREREREpFNMTBMRERERERERERGRTjExTUREREREREREREQ6xcQ0EREREREREREREekUE9NEREREREREREREpFNMTBMRERERERERERGRTjExTUREREREREREREQ6xcQ0EREREREREREREekUE9NEREREREREREREpFMGUgdAuhcfH4/AwEAEBAQgICAAX758wbdv3xAfHw9zc3NYWFigYMGCqF69OqytrXUen4+PD7y8vPDu3TtERUUhT548sLKyQvny5VGlShXkzZtX5zG9ffsWnp6eePv2LSIjI2FoaAgLCwuULVsW1apVg6mpqc5jIiIiIiIiIiIiyq6YmM7hEhIS8PjxY9y8eRP379+Hp6cnnj17hpiYGKX6lyhRAi1atMDgwYNRv379TIvz7du3WLt2LQ4cOIC3b9+meZ6JiQlatmyJ0aNH4+eff860eADg48eP2LhxI1xcXPDixYs0zzMyMoKjoyNGjBiBzp07Q19fP1PjIiIiIiIiIiIiyu5kgiAIUgdBmcPT0xMNGzZEVFSUVsarUqUK1q1bh0aNGmllPACIjY3F/PnzsXTpUsTGxqrUt1mzZti8eTMcHBy0Fg8ACIKA1atXY+bMmYiMjFSpb40aNbBlyxbUqFFDqzFpi7e3NypVqpRy7OXlhYoVK0oYERERERERERERSUXKXBFrTOdg375901pSGkhKdDdp0gRjxoxBfHy8xuN9+vQJzZo1w4IFC1ROSgPA5cuXUatWLbi6umocS7Lo6Gh07doVv/76q8pJaQB48OABGjRogN27d2stJiIiIiIiIiIiopyGiWlSiSAIWL9+PXr37o24uDi1xwkPD8fPP/+MmzdvahRPeHg4OnTogEuXLmk0DpC0ertLly44fvy4RuPExMRgwIABTE4TERERERERERGlgYnpXEQmk6FixYoYNWoU1q1bh4sXL8LX1xehoaGIjY1FVFQU3r17h6tXr2LZsmWoWbNmmmMdPnwYkyZNUjsWZ2dnPHr0SOFzrVq1wuHDhxEUFITY2FiEh4fj5s2bmDRpEvLlyyc6PyYmBj169MCbN2/UjgcAJk+ejPPnzyt8rl69eti1axf8/f0RGxuLyMhIeHh4YPbs2bCyshKdLwgChg4divv372sUExERERERERERUU7EGtM52J07d9CyZUu0bdsWXbp0QYsWLWBtba3SGBcvXsSIESPg6+srek5PTw+3bt1C3bp1VRpzy5YtGD58uKjdxMQEO3bsQI8ePdLs+/btW3Tt2hUeHh6i5xwdHXH16lWVYknm6uqKli1bitr19fWxatUqjBs3Ls2+YWFh6NOnj8KkdpkyZeDl5QUjIyO14tI21pgmIiIiIiIiIqJkrDFNmaJq1ar4+PEj9u/fj549e6qclAaAFi1a4NatW6hQoYLoucTERCxcuFCl8b58+YJp06aJ2vX09HD06NF0k9IAULx4cVy6dEnhD8i1a9dw4MABleIBgISEBIwfP17hcxs2bEg3KQ0AVlZWOHnyJJo2bSp6zsfHB6tWrVI5JiIiIiIiIiIiopyMiekczNjYGHny5NF4nMKFC8PFxQX6+vqi5y5cuICIiAilx1q3bh0+ffokap86dSpatWql1Bj58+fHvn37YGBgIHpu/vz5UPUmgAMHDuD58+ei9p49e2LYsGFKjWFkZIS9e/cif/78oueWL1+Ob9++qRQTERERERERERFRTsbENCmlevXqaN26tag9JiZG6Q0MExISsHbtWlF7gQIFMGPGDJXiqVy5MoYMGSJq9/b2hqurq0pj/f3336I2Q0NDLFu2TKVxbGxs8Mcff4jaP336hD179qg0FhERERERERERUU7GxDQpTVFiGgDevXunVP+LFy8iODhY1D5y5EiYmpqqHM/kyZMVtu/evVvpMV68eIF79+6J2nv27Ak7OzuVYxo1ahRMTEw0iomIiIiIiIiIiCinY2KalJZWovbDhw9K9T98+LDC9v79+6sVT5kyZVCvXj1R+4kTJxAXFydJTPnz50fHjh1F7Tdv3lSYlCciIiIiIiIiIsqNmJgmpRkZGSlsz5cvn1L9L168KGorU6YMypYtq3ZM7dq1E7VFRkbizp07SvVXVPYjX758Cjcy1CQmQRAUXj8REREREREREVFuxMQ0KS2tFb/FihXLsG9AQAD8/f1F7U5OThrFlFYC+dq1axn2jYuLU5jAbtCgAQwNDSWJiYiIiIiIiIiIKDdgYpqUdvnyZYXtDRs2zLDvgwcPFLbXrFlTo5iqV68OfX19UfvDhw8z7Pv06VPExMRoPaaiRYvCxsZGrZiIiIiIiIiIiIhyAyamSSkBAQE4duyYqP3nn39G4cKFM+zv6empsL1ixYoaxWViYgJ7e3tR+5MnTySLCQB++uknUZu3tzcSExM1HpuIiIiIiIiIiCi7Y2KaMhQdHY0BAwbg27dvcu0ymQwzZsxQagxfX1+F7aVLl9Y4PkVj+Pv7Z5gE1nVMMTExePfuncZjExERERERERERZXdMTFO6nj59imbNmuHKlSui5yZPnowmTZooNc6bN29EbQYGBkqtts5IkSJFRG1xcXEZJoEVxQQkleLIjJgAKKyzTURERERERERElNsYSB0AZR2xsbGIiIjAmzdv8ODBAxw9ehSurq5ISEgQnTts2DAsXbpU6bE/fPggaitUqBBkMplGMQNQWM85ec70NmZUFBMArSTL04uJiIiIiIiIiIgot2NiOhd6/vw5KlSooFZfa2trLFu2DIMGDVKpX1hYmKgtf/78asWg7DiK5szoeWNjYxgZGUkWkyo+fPiAjx8/qtTn1atXWpufiIiIiIiIiIhIXUxMk1Jq1aqF3r17Y/jw4TA1NVW5f2RkpKhNnXEUMTMzU3rOrB6TKtavX4+5c+dqbTwiIiIiIiIiIiJdYY1pypChoSGsrKxgZGSU4YaCaYmNjRW1aWNlMpAUnyIxMTHZLiYiIiIiIiIiIqLcgIlpylBcXBwuXLiAcePGoXjx4pg/fz7i4+NVHiM1AwPtLNhPKwmsaM6sHhMREREREREREVFuwFIeuVDBggUxf/58UXtUVBTCw8Ph6+uLe/fuITw8XHTOly9fMGvWLPz33384deoUChUqpNScenp6ok0U1V19nVpa4+jppf+5i6LnpY5JFaNHj0aPHj1U6vPq1St07txZazFkdYIgIDExEYIgSB0KEREREVEKmUwGPT09rWwGT0RElF0xMZ0LWVtbY8aMGemeIwgCrly5grVr1+LYsWOi593d3dG8eXNcu3YNVlZWGc5pZGSE6OhouTZVV12nJa1xMirLoeh5qWNSRaFChZT+YCC3SExMRFRUFCIjI/H161fRhyFERERERFmJvr4+8uXLBzMzM5iammp1IQsREVFWx8Q0KSSTydCsWTM0a9YMp0+fxoABAxAaGip3jpeXF0aMGIFDhw5lOF7evHlFienUx+r69u2bwnZjY+MMY0pN6phIPd+/f0doaCiioqK0tuqdiIiIiCizJSQkICIiAhEREdDT04OpqSmsra0V/q1CRETq8feXOgLllCghdQS6x8Q0Zahdu3ZwdXVF06ZN8eXLF7nnDh8+jLNnz6JNmzbpjmFpaYnPnz/LtUVERGglvsjIyDTnzCgmZcfSVUykuujoaLx9+5YJaSIiIiLK1hITExEREYGoqCgUL16ci1qIiLSkZEmpI1BObqxCyvuLine truncated "image/png": "iVBORw0KGgoAAAANSUhEUgAABacAAAVbCAYAAAA1KRZrAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAuIwAALiMBeKU/dgABAABJREFUeJzs3XdcU9f7B/BPWArKEMS9BfdAxTrq3hO31r3raOuoq65at63WWieOVusWd9UKiuJWXDhw40DcKHvJSH5/8IWfNzeEm0ES4PN+vfJq77n3nPPcEBCenDxHplAoFCAiIiIiIiIiIiIiMiAzYwdARERERERERERERLkPk9NEREREREREREREZHBMThMRERERERERERGRwTE5TUREREREREREREQGx+Q0ERERERERERERERkck9NEREREREREREREZHBMThMRERERERERERGRwTE5TUREREREREREREQGx+Q0ERERERERERERERkck9NEREREREREREREZHBMThMRERERERERERGRwTE5TUREREREREREREQGx+Q0ERERERERERERERkck9NEREREREREREREZHBMThMRERERERERERGRwTE5TUREREREREREREQGx+Q0ERERERERERERERkck9NEREREREREREREZHBMThMRERERERERERGRwTE5TUREREREREREREQGx+Q0ERERERERERERERkck9NEREREREREREREZHBMThMRERERERERERGRwTE5TUREREREREREREQGx+Q0ERERERERERERERkck9NEREREREREREREZHBMThMRERERERERERGRwVkYOwAiMqyIiAicPXs2/bhkyZLIkyePESMiIiIiIiIiIiJj+fz5M0JCQtKPmzZtCgcHB4PMzeQ0US5z9uxZdO3a1dhhEBERERERERGRCTp06BC6dOlikLlY1oOIiIiIiIiIiIiIDI7JaSIiIiIiIiIiIiIyOJb1IMplSpYsKTg+dOgQXFxcjBQNEREREREREREZU1BQkKAErHLuKCsxOU2Uyyhvfuji4oKqVasaKRoiIiIiIiIiIjIlyrmjrMSyHkRERERERERERERkcExOExEREREREREREZHBMTlNRERERERERERERAbH5DQRERERERERERERGRyT00RERERERERERERkcExOExEREREREREREZHBMTlNRERERERERERERAbH5DQRERERERERERERGRyT00RERERERERERERkcExOExEREREREREREZHBMTlNRERERERERERERAbH5DQRERERERERERERGRyT00RERERERERERERkcExOExEREREREREREZHBMTlNRERERERERERERAbH5DQRERERERERERERGRyT00RERERERERERERkcExOExEREREREREREZHBWRg7ADIdr169wr179/D8+XNERUUhKSkJtra2KFKkCKpUqYIqVarAzMww72ekpKTg7t27CAwMxKdPnxATEwN7e3s4Ozujdu3acHV1NUgcX4qPj8fVq1fx7NkzfPr0CUlJSXB0dETRokVRr149FC5c2OAxERERERERERERZVdMTudyDx8+xN9//42DBw8iKChI7bUFCxaEh4cHxo0bh5o1a2ZJPNeuXcOaNWtw4MABREdHZ3hdoUKFMHDgQIwZMwbly5fPkliA1CT5kSNHsG7dOvj5+SEpKSnDaytXroyRI0di6NChcHBwyLKYiIiIiIiIiIiIcgKZQqFQGDsIMrznz5/jp59+wt69e6HNS6BPnz5YsWIFihQpopd4QkNDMX78eOzatUujfhYWFpg0aRLmzJkDa2trvcSS5vbt2xg+fDhu3LihUT9HR0csX74cgwcP1ms8+nLv3j1Uq1Yt/TgwMBBVq1Y1YkRERERERERERGQsxswVseZ0LrR161bUqFEDXl5eWiWmAWDPnj1wc3PDpUuXdI7n1q1bqFGjhsaJaQBITk7Gr7/+ivr16+Pt27c6x5Jm+/btqFu3rsaJaQAICwvDkCFD0L9/f7UrrYmIiIiIiIiIiHIzJqdzmUWLFmHw4MGIiYnReaz379+jTZs2uHz5stZj+Pv7o2nTpnj37p1Osdy5cwcNGzZESEiITuMAwNq1azFo0CCdE8s7d+6Eh4cHE9REREREREREREQqMDmdi/z111+YOXOmynOWlpbo0qULNmzYgBs3buDVq1d4//49bt26hV27dqFbt26wsrIS9YuNjUW3bt0QGhqqcTxv3rxBt27dEBUVJTqXJ08ejBo1CseOHcPz588RFhaGJ0+eYPfu3ejRowdkMpmoz4sXL9CtWzfEx8drHEsaX19fjBs3TuWK8sKFC+Pnn3/GhQsX8OrVK4SGhuLu3btYu3YtvvrqK5XjeXt7Y/z48VrHQ0RERERERERElFOx5nQu8ezZM1SvXh1xcXGic7169cKSJUtQrlw5tWPcu3cPgwcPVlnqYvDgwdiyZYtGMbVs2RKnT58WtX/99dfYtWsXSpYsmWHf69evo0uXLnjz5o3o3Pfff49Vq1ZpFAuQWo6jYsWK+Pjxo+jcqFGjsHz5ctjY2GTYf/v27Rg4cKDKc/v27UOPHj00jikrsOY0ERERERERERGlYc1pynJz5sxRmZhesGABvLy8Mk1MA0DVqlVx7tw5uLm5ic5t27YNQUFBkuPZs2ePysR027Zt4efnpzYxDQDu7u4YMWKEynPr1q3DnTt3JMeSZubMmSoT03PnzoWnp6faxDQADBgwAMWLF1d57scff1T5/BMREREREREREeVWTE7nApGRkfDy8hK1jxs3LsMyHxmxsbHB3r17YW5uLmiXy+XYvn27pDHkcjmmT58uai9Xrhz27dsHS0vLTMf4/Pkz/vrrL5XnUlJSVI6vTlBQEDZs2CBq7927N37++WdJY9y9e1flSm4AePnyJdasWaNRTERERERERERERDkZk9O5wKVLl5CYmChos7Kykpx0Vebi4oLOnTuL2lWthFbl8OHDeP78uah9/fr1yJ8/v6QxNm7ciNevX2d4/vjx43j8+LGksQBg1apVkMvlgjYHBwesXr1a8hhz585VWas6zerVq5GSkiJ5PCIiIiIiIiIiopyMyelc4MmTJ6K2WrVqwcnJSesxW7ZsKWp78OCBpL5r164VtTVp0gStWrWS1P/z589YsmSJ2msUCgU8PT0ljRcfH4/NmzeL2seNGwdnZ2dJY9y9excHDhxQe83Lly9x7NgxSeMRERERERERERHldBbGDoCynru7OxYvXixoc3Fx0WnM0qVLi9oiIiIy7RceHo4zZ86I2sePHy95buVV061bt0aDBg0wb948wXUHDx7E8uXLMx3P19cX0dHRgjYLCwt89913kmNSXjW9cOFC+Pj44Ny5c6KYPDw8JI9LRERERERERESUUzE5nQs0bNgQDRs21OuYFhbil46UkhzHjh1DcnKyoM3e3h4dO3aUNK+qVdO//PILChYsKEpOv3jxAnfv3kX16tXVjnn48GFRW6tWrVCoUCFJMSmvmnZ0dMQPP/yAggULipLTx44dg1wuh5kZP7RARERERERERES5GzNkpBVVpUKKFy+eaT9Vq6ZbtWqFPHnySJpXedV0q1at0LBhQ1SoUAHlypUTXX/27FmtYpKaLAfEq6Z//PFH2NraokOHDqJrQ0NDcf/+fcljExERERERERER5VRMTpNWjhw5ImqrX79+pv1u3LghamvWrJmkOTNaNa1uHFXzfSkiIgLPnj3TOiZVq6bHjRsHAChRooTKhHlmMREREREREREREeUGTE6Txs6ePQtfX19R+zfffKO2X2JiIu7duydqr127tqR5Va2a/vrrr9WOk1kiOCAgQLDqGQDy5s2LKlWqSIopo1XTusRERERERERERESUGzA5TRp5/PgxevbsKWqvW7cuWrZsqbbvy5cvkZSUJGqXkgjObNV0RuOoWhX9padPn4raKlasKKkmtLpV07rERERERERERERElBswOU2S7dy5E40aNcLHjx8F7dbW1ti0aRNkMpna/iEhIaK2AgUKwMHBIdO5M1s1DUBlCY3Y2FhERkZqFJOqcVTJbNV0RmN9eR9ERERERERERES5lYWxAyDTlJycjPfv3+P169fw9fXFvn37EBAQILrOysoKe/bsQY0aNTIdU1VStnDhwpn2U7Vqes6cOZLHev36Nezt7fUak/Kq6QIFCuCHH34QXVeoUCFJcxIREREREREREeU2TE4TAODVq1coWbKkRn3Kly+PHTt2oF69epKuV7WCuWDBgpn2U1413bJlSzRq1Eh0Xd68eZEvXz7ExsZmOq+uMalaNW1nZye6ztnZWdKcREREREREREREuQ2T06Sx2rVrY8SIERg+fDisrKwk90tISBC1WVtbq+0jpdb0l/LmzStKTsfHx+s1JlWrppVrTX8Zj7LExETI5XJJda0z8+HDB4SGhmrUJygoSOd5iYiIiIiIiIiIdMXkNGmsVatWcHV1hVwu16ifqiRxnjx51PbZtGmTpFXTaVQlg1UloHWJad68eZJWTWcUT1pMNjY2aueRYu3atZg7d67O4xARERERERERERkaN0Qkjf32229o3bo1ihUrhuXLlyMxMVHrsdRtoqjpqumMxvsykaxrTIGBgdi/f3/6sbpV0+rG0jQmIiIiIiIiIiKinIYrpwkAYGtri0mTJqUfy+VyxMbG4vXr1wgMDERwcLCoT3h4OCZNmoRNmzbBx8cn05rVqlYRf/78OcPrN23ahFevXqUft2jRQu2qaUD1KumMVi9rE5NyremJEydmuGo6o3gyi4mIiIiIiIiIiCg3YHKaAAD29vZYtmxZhuefPn2KzZs3Y9WqVYiKihKce/DgARo1aoQzZ86gbNmyGY6hqpZzRslbbVZNZzSeuhrSmsSkatX0+PHjNY7H0tIS5ubmavtJNXbsWPTq1UujPkFBQejatate5iciIiIiIiILine truncated
"text/plain": [ "text/plain": [
"<Figure size 1500x1500 with 1 Axes>" "<Figure size 1500x1500 with 1 Axes>"
] ]
@@ -985,13 +998,13 @@
"handles = []\n", "handles = []\n",
"\n", "\n",
"eb = plt.errorbar(voltages, damage, damage_err,\n", "eb = plt.errorbar(voltages, damage, damage_err,\n",
" fmt=\"o\", capsize=5, capthick=1.5, elinewidth=1.5, label=\"Current study\", markersize=MARKERSIZE,\n", " fmt=\"*\", capsize=5, capthick=1.5, elinewidth=1.5, label=\"Current study\", markersize=15,\n",
" zorder=ZTOP, color='blue'\n", " zorder=ZTOP, color='red'\n",
" )\n", " )\n",
"handles.append(eb)\n", "handles.append(eb)\n",
"\n", "\n",
"for i, (v, d, mk, lab) in enumerate(literature):\n", "for i, (v, d, mk, lab) in enumerate(literature):\n",
" ln, = plt.plot(v, d, mk, label=lab, zorder=ZTOP - 1 - i, color=plt.colormaps['inferno'](i / 6), markersize=MARKERSIZE)\n", " ln, = plt.plot(v, d, mk, label=lab, zorder=ZTOP - 1 - i, color=plt.colormaps['gray'](i / 6), markersize=MARKERSIZE)\n",
" handles.append(ln)\n", " handles.append(ln)\n",
"\n", "\n",
"plt.xlabel(\"FIB milling voltage (kV)\", fontsize=14)\n", "plt.xlabel(\"FIB milling voltage (kV)\", fontsize=14)\n",
@@ -1010,23 +1023,23 @@
}, },
{ {
"cell_type": "code", "cell_type": "code",
"execution_count": 24, "execution_count": 23,
"id": "e3824127", "id": "e3824127",
"metadata": {}, "metadata": {},
"outputs": [ "outputs": [
{ {
"data": { "data": {
"text/plain": [ "text/plain": [
"[<matplotlib.lines.Line2D at 0x7fb0064d1fd0>]" "[<matplotlib.lines.Line2D at 0x7fc3e1359af0>]"
] ]
}, },
"execution_count": 24, "execution_count": 23,
"metadata": {}, "metadata": {},
"output_type": "execute_result" "output_type": "execute_result"
}, },
{ {
"data": { "data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcYAAAGoCAYAAADcnEUDAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAPYQAAD2EBqD+naQABAABJREFUeJzs/Xe4Zcdd5wt/qmqtHU8+3X06J2WpJUuWJRtnOYJJHsCADVySYV7eGeZe+73jeTyBYJjx4DswzMNzB2ZgDHPxYDDGA7ZxxkGWo2wrSy11VOfuk9NOa62q948Kq9Y+RwZzsexndOp5du8+e6+9VoVf/cL3F0oYYwxbbattta221bbaVgNAfqs7sNW22lbbalttq307tS3BuNW22lbbalttq0VtSzButa221bbaVttqUdsSjFttq221rbbVtlrUtgTjVttqW22rbbWtFrUtwbjVttpW22pbbatFbUswbrWtttW22lbbalHbEoxbbattta221bZa1LYE41bbalttq221rRa1LcG41bbaVttqW22rRe3bVjDee++9vOY1r2FycpJ2u82dd97Jn/zJn3yru7XVttpW22pb7X/xlnyrO7BZ+/SnP82rX/1qarUaP/qjP8r4+Djve9/7+LEf+zFOnz7Nv/yX//Jb3cWtttW22lbbav+LNvHtVkQ8z3Ouv/56zp07xxe+8AVuu+02AFZXV/mO7/gOHn/8cR599FGuueaab3FPt9pW22pbbav9r9i+7aDUT37yk5w4cYI3vOENQSgCjI6O8m/+zb8hz3P+8A//8FvYw6221bbaVttq/yu3bzso9dOf/jQAr3rVqzZ85z/7zGc+83e+n9aaCxcuMDo6ihDiH6SPW22rbbWtttWenmaMYXV1ld27dyPl02PLfdsJxmPHjgFsCpVOTk6ybdu2cM1mrd/v0+/3w9/nz5/nxhtv/Ifv6Fbbalttq221p62dPXuWvXv3Pi3P+rYTjMvLywCMj49v+v3Y2Bjnzp17yt+//e1v51d/9Vc3fH7r9/9LUhrIzCAKSHoFIjfIXCNyDYDAgPe4GvsSuQZjQBuE1hgBCGFfSmCUBCEw0n3mvt/guBXuflqH+4p+bv/uDxBZHu5rpMC0m5hmik4kRSMBKcpn+w4K0EqgFRglKOru+f6RGmTm3w0yNwhjEIUdpygMQhtErpH9HAqNGGSIQWZnw4+tlmDqqe2Dkhgp7O+MDvOEASQY37/w7rvrnqsNIisQ/czOq5/5eoJupSAlRV1hEomR2HkFhLZ9tp+BEaBTgUlEuVbY8dpxgszcXGv73DBtxiCzAjEoINeIwQChDUYpSOyzi3YNk7jxCzuHIisQfg21G5O2NzVKluP3dKA1orDvspdBXtjP8hykxNRrkCh0I6VoJnbdE2nXOaIb498l6ESU9IlBFm4dDQg/D7ocsxxYOhd5gegN3DrYfiAlJnG0lSaYVNnxKlGun1sj4R5qYgKTAiOEWxtTWScxyC1NG4Mo7B4ySkGaYJRdbyMtLSEjmsbSslZuHyT2GUYJjCrH59dVxutgjKXzrAhjp9B27Flm6SdJQElMqtCNxK2vXS9R2L1g505bmhEClLR7LZFhzGWf3TIVGjEo7Ht3gOgPHC+wY9QjDXS7DlLY+/jnFW7MzhDKm47+/VoLR9OZ3auqr5F+H/m5LuyEiEzbedcGssx+HniKxDTrmJqy81pTgUcJKO9XLgMGwyMn/idzi0dp1CeZmrgKo3PmFp8gyztcu+Mu9o0csXOc5RQ6496F97M2mGOyfYCR8d2sda+wuHiC1uhOjtzxs6Skbr3sy0iBTuweM4lEJ7ZHwvVB5Abd7/LlT/97RkdHebrat51g/H/b3vrWt/LmN785/L2yssK+fftor6SkSRo2lNAKgUEIg/A7zVARLAgQaSnigmCSQ4IQAiOJhZe9xn3mLhTGEXNWIFVuN5JJEOS2D24zalnDiBpGSkRihZJnksIJNStksExJCnQnEszCPzIakHLMS7ku55ZpCalRRQZoRJIgTGo3v5IR03SMOy0FoxcyVjCZchMKqoLcKwWO6QiVI40TjO5lVIKRNYwUSKnsRvZCCb9m2o5XlwzUKDc+z6gixUYIz9R1UHo8ExeFQGjbMSEEKDD1FJMm6JrEtGuWgbkZFIVBicIxM22VAkHw0ptEOmHuFQr3LHe9xDGqLEeI3DGrBIREk6BMijECjRXGZV/LebPCsSoYvTBEuPnw37vbSJkjRIFlNtoKdWz/TKKgUbfKTi2xaxspd6LQFYUiTIaj9yBATXlNoG9V2JfW0M+sUim9ViPRpFZoOKYtnPIZlKy8fIad03Jeyw3pFRP/XIPAIIUGYRAit3tbFghpid6kCSiFThUiTdy+cntTanu9LoV52O9CoKJ9H/a2nwNtkMKOV2iF1CqsEUJQyDpG1p3QV259LP8JSqMxpD2JyZ1SqKySJEw5VoG23Y1+K7VVBESeQ6btXBQKtLR7WCUYJdG1ul1nJTGpI1ztaCi6n1/Ly4uPMrd4lPHR/dz6rJ9BpjWENmS9Ne69/3c5Pns3O0evp5mOIPIep9a+xtpgjkNjz+HaXS+jGGthUsWJC5/i9OlPcvn0Fzl06BWW72qNzLWlN7fWOvH7WZTRLwZy4wyXp9EV9m0nGL2l6C3H4baysvKU1iRAvV6nXq9v+FwWJjDKWNvCa7sm+qwiHL3AE6VwSqoMpNIiDT8rupw6/TesrJ6j11sky7qkaZtWaxv7dt7JzMh1gEQ4ARQ2kin7aPy764sQBG0rMAbXb4l/tnHdEpigekVC3DMvXWpu5QbBMjABSNcv6cfq7yMw0jIjY5yQFIJef4XLC48wv/QE6905BoM10qTJ+PgB9u97EROtPeEeWhRcWXmC2dXjLPcu0M1WEULQbm5n18xt7N77XIS0pqEVjDowD69lg7Qy1wneoICIaB5x6xuNE2Mq626iNTZKWEsqYphesFsmaJ/rySTLupy48ClWOhfo9pfI8i5p2qLV2s7eXc9lx8QN9jZSWCtTCHIz4NjCF7i8fox+sU49abNj+iYO77sLqdrOQi3XN/TbM0m/hm54Iqy7KEnXKyxaWyvVKRb+PjE64RW9IFSj+RBDLh0j4HNf/S16/SU2a3tm7uD6q74PIYVFVLA028mWOTH3ZeY7p+kX66SqQbu5nb0zd7Jjx82V/eMFgVewPH0h/Zr6zhCsnEAbWoMX6H68hnK/ynLMXigGC9312V8jfB8CCiAwGK4sPsrZi1+i050lL/o06uNMjB/k4MwLaKeTbhCCY4uf48TSFzedJykSXvqiX43G4vafU6bQANoKA+PXuRy/VyTsu7a/KZziYxx1+r2sopdf77DUdkxCliiXV7ZmFx8D4OD+l6BUarshIK2PsG/38zl26kOcW32IayZfgJaCc51HUSLlqvHnRWsh2L/vRZw79wUunr+XQwdeZve05zlBgTcIEUlD38FojZ/O9m0nGL1v8dixY9x+++2V7xYXF5mbm+P5z3/+N3xf2SuQaVHCgN7iKhyUqi3kEjRkb3l5rVJJdGoJS6cyaIwl5OT+MSVfHvS6XLz4VcbG9rFt240kSZNBvs783FEeOvqnzM88hxsPfS8mS5yWXARNNVgXKORABmZtpBNoualsDs+AApcMQs5rtnLIqjLIgbbwoINTg4bsITUPJSZyiIECOMgvUijOnv0ST577LM3GFFMTV1Ortel055mbf4zZuUc5ct0PMzN9BCEFnWyZB87/JUrWmG4dYPvoteQMmF05xuMnP8j8whM869o3QGLxMw/7xoLRJAZTCLupvTXh4Ey7Bg5izLWFVbWBwsKhlnk6c8P91qTKvpJyrr1CYQ0dVVEoMNBbW+bC3P2Mje5l+/QNJGmTQbbO3PxRHn7s3eyeeQ43HP4+O39KUmQ9vnzhz1gdzDLdOsCu+g2sZLOcufQFFpdPcfstP4dK66Xy4mnBK21DIL2IyNUYB+nGik8/RwxyC+Pmhf3OwYkkykJqzoIwypu/TthIK9hMZBkZYXmWShrs2/v8kpwcaY2N7kXXVan0C5jvPMn9F94HwPbWVTRrE2QiY7U/y8LyCXbsuLkUwk4oBqYZBLnBGg6lIioiyy6gEbmlaYxx0LWxdJ9Ia6HWEoyDywNdu9vabaPC36GZ8nnHj32Is+c/R602yvbpG1FJnbX1S1y49FUuzz7EHTf8LGNywimX9ia7R26k0Zyy8+1cMELICnTphaK19uzYTS4C/W1onncV1iVD4dY3z+33Utk1riWYRs1ZjI6+vaIGbi+L0qvh195AP18HoN6asha7V7YlNEamAVjonEHvfind3hx9vc62xgESmWBcv0SuqYmUifGDzM0/Rm9tllZjWwkjO2sQrPFikthiFKANKss3mYBvbvu2E4wveclLePvb387HPvYxfvRHf7Ty3cc+9rFwzTfaRKERqiRET4SB2RqncTlcPmiKSal9UYF1qtqm8OpUgBah2ZrkhS/810iZBIvLKEGmB9z3xf+bC5e/wr5dz2NcTlQgHcBpvCYQUOndEaUVEffffuXG5zRO4ZibsNLSWzmB+fjxxlaUjCym2AcU+WJKYSEiDRZGx/bx7FveyMT4wQqUurR8ivvveydHT7yf7VPXI4RAJXVu2P2d7Bm9CYXz+yhJvusVfOXEHzO39ASzcw+zY/vNbv1iDdmEOQrWkjQIIt9uUFC8pajDfOK17Aj+9ZBZaS269fLKhr/GWAHkfVqN0Wle+MJ/jRKqoqnnRZ+vfu33whqPqWkQgtPzX2J1MMuh8Tu4bvtL7f3ThGPz93Dy8md58txnOXzw5UEoxuhBWNu4BWEpnDAxwRKOmae3HPF04ulNWqXL/t/fz9GL16Ni36ej0SRpcPDqV4b7eavd6UoB+uz3Fnng4l9ST0a5Y9cP0VSjVhA365Ao8mQzayAaq59TYxABVyzpXMT7RFuapihCTADaYIQEoZxiJ4Nvs3R1uIFGPu2Line truncated "image/png": "iVBORw0KGgoAAAANSUhEUgAAAbEAAAGhCAYAAADvBooqAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAPYQAAD2EBqD+naQABAABJREFUeJzsvXm8Xddd2Pv9rb33Ge+gK+lqsizZsuU5jhM7kyEDCSGEJHbIWAI00CaQltLHo7wHgbYQHpQCfeWFQl8THgRaSPsAJxAgAwFC5sGxHSee7ciSJdmapTudae+91vtjjfvcmxSKyPuoH63P5+jonrPP3mut32/95kGMMYaL4+K4OC6Oi+PiuACH+v97AhfHxXFxXBwXx8XxPzouMrGL4+K4OC6Oi+OCHReZ2MVxcVwcF8fFccGOi0zs4rg4Lo6L4+K4YMdFJnZxXBwXx8VxcVyw4yITuzgujovj4rg4LthxkYldHBfHxXFxXBwX7LjIxC6Oi+PiuDgujgt2XGRiF8fFcXFcHBfHBTu+4Uzswx/+MDfeeCP9fp9bb72Vu+666xs9hYvj4rg4Lo6L43+S8Q1lYvfccw+vfvWredOb3sSnP/1piqLglltu4dixY9/IaVwcF8fFcXFcHP+TjG8oE/vFX/xFXvOa1/ATP/ET3HTTTWzbtg2AX//1X/9GTuPiuDgujovj4vifZMg3sgDw1q1b+c3f/E1uv/12HnvsMW677Ta+93u/lz/+4z/mc5/73Nf9rdaaJ598ktnZWUTkGzTji+PiuDgujovjfA1jDCsrK+zatQulzo8O9Q1jYktLS2zatIkHHniAa6+9lre97W3ccsstLC4u8ta3vpUTJ040rh+Px4zH4/D30aNHue66674RU704Lo6L4+K4OP4ex+HDh9m9e/d5uVd+Xu7yNxgrKysA9Ho9Tpw4wZ//+Z/zzne+k0996lPhu3T8wi/8Au94xzvWfX7T7T9JQQdVGqSGfFQjlUFVGqk0AIIBz5qNfUmlwRjQBtEaI4CIfWWCyRSIYJT7zH2/jsOLu5/W4b4yruzf4wlSVuG+Rgmm38V0C3SuqDs5KInP9hMU0JmgMzCZULfd8/0jNajSvxtUZRBjkNquU2qDaINUGjWuoNbIpEQmpd0Nv7ZWjmkXdg6ZwiixvzM67BMGUGD8/MK7n657rjZIWSPj0u6r3/l2ju4VoBR1O8PkCqOw+wqItnO2n4ER0IVgcomwwq7XrhNU6fZa2+eGbTMGVdbIpIZKI5MJog0myyC3z677LUzu1i92D6WsEQ9D7dak7U1NpuL6PR5ojdT2XY1KqGr7WVWBUph2C/IM3Smou7mFe64snBO8Mf5dgc4l4icGVTs4GhC/DzquWU0snktVI6OJg4OdB0phcodbRY4pMrveTCL8HIzEPdSkCKYEI+JgYxpwkkllcdoYpLZnyGQZFDkms/A2yuISKsFpLC7rzJ2D3D7DZILJ4vo8XFUKB2Msnpd1WDu1tmsvS4s/eQ6ZwhQZupM7+Fp4SW3Pgt07bXFGBDJlz1quwprjnB2Yao1Mavs+nCDjiaMFdo16poPut0GJvY9/Xu3W7JSLquvw38NaHE6X9qxmY43y58jvdW03REpt910bKEv7eaApCtNtY1qZ3ddWFmiUQLxfBEOTlnnao7VdY41dZ1naPS4r+9NWYXG6XVDPtTG5QreUXbPGzt04uqON3ZvcnjGTK3RuHyZuDlIZVG2oqhFf+Ot/y+zsLOdrfMOY2MzMDACDwYB3v/vdvO1tb6PdbrO2tha+S8fb3/52fvRHfzT8vby8zKWXXkp/uaDIi4D8ojMEg4hB/KkwNJgAAlJEdhSYiJpiWhAOfcpo7DXuM3ehGId4ZY3KKosQJkeo7BzcwdGqhZEWRikktwzEEzRxDMgyBCwBUYIeJExU/COTBWWO0GRuypUlMKI0WV0CGslzxBT2oGYqIXCOyBaRiXmGYJmIiQdGaDJdz8AdgZCsQhnHxNzLZDlGtTBKUCqzh84zEDzMtF2vjsTOZG59nqgkQoiIJ8A6CCie4EotiLYTExHIwLQLTJGjWwrTb9mD53ZQakMmtSM82jJwIXiHTa4c4/XM3z3LXa9wRKWsEKkcYclBFJqczBQYI2gs44xzjftmGVmTiXnGhbj98N+72yhVIVJjCYO2DBg7P5Nn0GlbwaSVW9gmgpjUusH8w2Y4fA/MzsRrAn5ntX1pDePSCoDKSyAKTWEJvCOw4gTFIBBV8Rl2T+O+xgPphQj/XINgUKJBDCKVPduqRpRFelPkkGXoIkOK3J0rdzaVttfryHjDeRchS859ONt+D7RBiV2v6AylswAjRKhVG6PajkFnDj6W/gQBzxiKkcJUToDLrEAjJq5V0Ha6yW+VtkxbqgpKbfeizkAre4azHJMpdKtt4ZwpTOEQVzscSu4X4K2SdeaJ8Opoj6JAqcridDWyv60zQGEyRV0VGMnQJqMW5ciDhZFojaq0xTcHa5378ywx6sLhg/JC0nl0CX3DmNj8/DybN2/m7rvv5r3vfS/33nsvAA899BD79u1bd3273abdbq/7XNUmELVUisFLkSb5rMHIPHOSyEjy5mFvjERyDi/iu8FipWSWSRhAHLMISG/iHI1/d3MRIUgx4RC7eSv8s42blmCCSJMwXE9odJSIIjJjiY0Ays1LSbIef4gt4TDGMbQGE6fBwDwjatxD4bSZZO+0nbfUTq0QA0YcE9PhoHvpFZTlj45JBmFBkn3EwTdZJ8Y04G4SGJtMrIaSEDfPhC3Bss8NaOLPvJfQp/AgrFuJXa8k93WwCHAGRLv9SuAb5u0Jmoehe74EuEtEXS9caG2lZycE+PukWr8nVoEBJvshU+6HdI+nrw800IAosZYKpGE8CGctoZeNkRBtLwwF/FKGBr54Rg8RN7QGz3z9eg3xvKq4Zs/Agubr5uyvET+HoF17fE7OkodveHdILh5p3FfG2PVq938tU2tx588JPmgAbYm239Op+9kz69ZZG6sRacd8A4w8I1NBKA0CJpYpGAyiovXIC0bGXffE0U9x5KnPU1ZD5uZ2s3/fdzDT3W7hmysOnrqTJ858kUpP2Nrdy7WbX0yrNedgIVCLFRq9kJPgtAQaaRBJOJefoF+v1hsgy99tfMOYmIjw4he/mO/5nu/hx3/8x5mbmwPgjjvu4CUvecnf+D5qVKOKep00IbUzJ2prdgiSp0Ne46W1TKELiwS6UEFCiWYX949JiJtHvGQYd0i0l7ZrwZS5kz7rIAEGqZ0MNVGBsBrlmE9lGogciIXEZwdCIzhzV3OuaqKticyZFIPk6c1K3pyWqyliB+DMXgnzl3SpCQPzJgnxmpVy0re4U+2RuqqhFmty8M/KLGv2ps+UiZncYGqx9/NSujPpWRg4M1ulrWlRG6itSdASOncw3G9NkdlXHvc6rEOBUVmD+XuTbIA5kYkHHFKCMY5wZsqZftzQjmNVGim1JdjiNEQvaHhc8ALWFD4FGQUwxpk1UyFlXCGTyu5tVdvvnEmNPLNmJSeZ270mcBdL2JypJ0sIvjSfbZj6zOGecp976dmbYY2xBNdrOXbfHMP0jM4TuMB0DUYnD8Ktz58Xr+VXFqcxDp+0xSVyBeI0TmcyDnjtbmunnoW/wzDJ89bheiq4uPPp8dszGHCaUm0Fw2RP/P09A7NalF27qSLjXDc87aqtW4LawbeyZj1UZmHcyjGdltPEHH57oQrcWZYoUHjYu78fP/Qxnjj8Ca655jV0+4s8cejj3PPl3+LZz/5f6LQ7HD31RQ6c+RzXL76UTjbL54++l1PDQ7xk3w9b3Ku046cS8dedG3t24nlQtbHugRR5PH2r6g024e82vmFMDODHf/zH+cM//EPG4zH33nsv73nPe3j44Ye54447/sb3kFojWXLq3OYEwugOlrcjBwksj1INDdNGU4qzhzcSH09ETOOEY4msslKlUY46JlJhGF5SccCO3giJ0nk6/+QxOGkHzzCcVGj8JZ5Q+PWm2omSqImkPovEdxAJe5QmoxQeqWrDnOi2R8S/yzqlN2gbaNCRwFn4pZKnCdd7no0ySCoyBGHCa2A67CeeuCQm0EBQVSqgJEzJX2MsEQ4+GGUSM9j0itzsGpK8v4/TRojStEFZnPE/C9pU6q9oMrHI2NyB94DwjCDAOGHanohZe6MjrJKq8hZfPNxSX11iHk4FPhMAYd+9+U+0nmIIcX5W+JpaT1xUvBa/X9KUlAJsEyGn1pagJ5qnEQWSOcFIBV9cNPe7CSY+2CDABGEC8L5zQ9CcNxwbuBkCHLHz9T6wpu/RnlvxOGPcw1IzfQRP9EXqFMae6RPga/Is+O/XrXndIuyaDWC05vCRT7P/qttY3HEjIFxz/Ru46wu/xpPHvshlO7+ZA09+kv2XvIQdvevRdUkuLSoz4fTgEFtaV9pbJr6/sBkePzENPKJWEb+cECMp/p/H8Q1lYrfccgt/+qd/yk/8xE/wrne9ixtvvJEPf/jDf6solai1gLed2y+Si9IDBkFrCwjnzQQeqVxQBeLfiYhmICsN1DVHHv1rnjr8BcrJKrNzu9l/3Xcy01mMj24QUS8SaquZGYPCBEneuPBSSW3Xfh3pgfGHvdbBHOeZjjc5yMRJ6elGZHnURLwjW6nAuOuO/b/OJQkysM8yVcWRB/6C44e/SF1N2Lz1Kq66+jbaefRdnj77KI8d/AiD4Wl67c1ctfWFLPYud6YQdwCV8wUoFRhNw/eWwNS7pvDEX+LfJIQtaJp1HQi/s1Y2hQjHoFRl761bULsgk7rttTMwYn+vJvb+T93/MY4dcjCe3c3VV7yK2c42pNbocswDj3+Q48sPoci4ZOY6rlr4JitEQfA/WkEqQSLHyBqBKQ2Au7/9emptGawhErSqjuY1h2vWxOS0kSzxOzlGZvdRUbcEnYPJbdCQn9ZTD/w1h7/0QdLRas/xvBf9JNSGyXiFxx78AKfPPEKmCnbP3cjLine truncated
"text/plain": [ "text/plain": [
"<Figure size 640x480 with 1 Axes>" "<Figure size 640x480 with 1 Axes>"
] ]
@@ -1059,13 +1072,13 @@
}, },
{ {
"cell_type": "code", "cell_type": "code",
"execution_count": 25, "execution_count": 24,
"id": "2e8a5fad", "id": "2e8a5fad",
"metadata": {}, "metadata": {},
"outputs": [ "outputs": [
{ {
"data": { "data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAooAAATDCAYAAADx8xQSAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAuIwAALiMBeKU/dgAAs19JREFUeJzs3Xd4lFX+/vF70hNCQoISelmagIAEpCqwgBUQWJciTVelqlj4WhcURVzXsigWWEBXASkuuoAVQhBEEBcQpBmKGIK0hB7SJmV+f/gjS5gTkufJJJPyfl0Xl+bMc875TEjCnaec43C5XC4BAAAAl/HxdgEAAAAonQiKAAAAMCIoAgAAwIigCAAAACOCIgAAAIwIigAAADAiKAIAAMCIoAgAAAAjgiIAAACMCIoAAAAwIigCAADAiKAIAAAAI4IiAAAAjAiKAAAAMCIoAgAAwIigCAAAACOCIgAAAIwIigAAADAiKAIAAMCIoAgAAAAjgiIAAACMCIoAAAAwIigCAADAiKAIAAAAI4IiAAAAjAiKAAAAMCIoAgAAwIigCAAAACOCIgAAAIwIigAAADDy83YBQFnhcrm0Z88e7d69WydOnFBKSoqCg4N19dVXq1mzZmrZsqX8/PiWAgCUH/yrBhRgz549mjFjhv7zn/8oMTEx3+PCw8PVt29fPfTQQ2rfvn0JVggAQPFwuFwul7eLAEqj8+fP68knn9Ts2bOVk5Njqe+gQYM0Y8YMRUVFFVN1AAAUP4IiYHDw4EH16dNHP//8s+0xateurRUrVqhNmzYerAwAgJJDUAQuc/jwYd1www1KSEgo8liRkZFau3atWrZs6YHKAAAoWTz1DFzC6XRqwIABxpDocDg0ePBgffnll0pMTFRmZqZOnTql2NhY3X///fL393frc/r0afXr10/nzp0rifIBAPAozigCl/jrX/+ql156ya396quv1scff6zu3bvn23fnzp3q37+/Dh486PbayJEj9eGHH3qyVAAAih1BEfj/Dh48qObNmysjIyNPe6VKlbRhwwa1bt26wDEOHz6s9u3b6/jx43naHQ6HNm7cqI4dO3q0ZgAAihOXnoH/7+WXX3YLiZL02muvFSokSlKdOnX0r3/9y63d5XJp6tSpRa4RAICSxBlFQL/fS1izZk23oNi8eXPt3LlTPj7Wfqe69dZbtXLlSrf2uLg4NW3atEi1AgBQUjijCEj697//bTyb+Nhjj1kOiZL0f//3f8b2BQsWWB4LAABv4YwiIOmmm27S6tWr87QFBgYqMTFRYWFhlsfLyclRnTp1dPTo0TztTZs2VVxcXJFqBQCgpHBGERVeenq6vvvuO7f2rl272gqJkuTj46PbbrvNrX3v3r06fPiwrTEBAChpBEVUeJs3b1Z6erpb+5WWwimMP/7xj8b2devWFWlcAABKCkERFd6PP/5obG/btm2Rxm3Xrp2xfdu2bUUaFwCAkkJQRIW3Y8cOY3uLFi2KNG6jRo0UEBDg1r5z584ijQsAQEkhKKLCM+2kEhwcrFq1ahVpXF9fX9WvX79Q8wEAUBoRFFHhHTp0yK2tRo0acjgcRR67Zs2abm0JCQnKyckp8tgAABQ3giIqvMTERLe26tWre2Rs0ziZmZk6e/asR8YHAKA4+Xm7AMCbnE6nUlJS3NrDw8M9Mn5+45w+fVqRkZFFGjsxMVFJSUmW+pw/f15btmxRWFiYqlSpojp16igwMLBIdQAAik9GRkaeZdW6deumKlWqlNj8BEVUaMnJycb20NBQj4xfuXJlS/Na8e677+r5558v8jgAgLJj2bJl6tevX4nNx6VnVGhOp9PYbnpa2Q5/f39ju2m7QAAAShuCIiq0zMxMY7ufn2dOtucXFPObFwCA0oRLz6jQfHzMvyt56qnk/MbJb14rxo8fr4EDB1rqs2fPHg0aNCj348bBtyvIt/D3Y0b5WN/SMMTf3nt1yfo29CF+9uaKDLD+hLvdr5CMbOt9UrOsfy4kKT3bXr/zmVmW+5zKuWBrrhyH9U+In8v8C1hBkn3OWe/jcn/YrTB8Zf2qRDNHc1tzXR1Ucv+Uh9n4XpEkPxvdnDa/yU6lW+/4W4b7verFJchh7es3NfuMtqV+mvtxnTp1PF3SFREUUaHld4k5K8v6P5RWxvHEpe1q1aqpWrVqRRojyDdcIb5VC318mG8Vy3OE+vla7iPZC4qhNkNptUDr/4rZzGBKtxEULzjsTZZqs58rx/oZ7zSHvYeish3Wv9f8Xfa+fzJ9rP89p7vsBQg/BVnuE+a42tZckflcuSgOETaDop1vTTu/VElSTrb1oHgmq+j3jRdWiKNoP/9L+gFELj2jQgsKMv8wT0tL88j4qampxvbg4GCPjA8AQHEiKKJCCwsLM96PeP78eY+Mn9/TzRERER4ZHwCA4kRQRIVnCm2nTp3yyNgnT540thd1DUUAAEoCQREVnmlP5+PHj3tkbNM4kZGRXHoGAJQJBEVUePXr13drS0xMzPf+Qivi4+Pd2ho0aFDkcQEAKAkERVR4TZs2dWtzuVzav39/kcY9d+6ccYu9Jk2aFGlcAABKCkERFV6bNm2M7du3by/SuNu2bbM0HwAApQ1BERVe+/btje2bNm0q0rjff/+9pfkAAChtCIqo8Bo0aGC8bzAmJqZI45r6h4SEqFOnTkUaFwCAkkJQBCTddtttbm2//PJLvpePC5KYmKhvv/3Wrb1Xr14e2ZUFAICSQFAEJA0dOtTYPmvWLFvjzZkzR9nZ7vtP5TcPAAClEUERkNS5c2c1b97crf3DDz80LnFzJefOndObb77p1l6tWjX169fPbokAAJQ4giIgyeFw6PHHH3drz8jI0OjRo+VyuQo91sSJE43L4kyYMCHfvaUBACiNCIrA/zdixAi1bt3arT0mJkYPPfRQocLiq6++qvfee8+tvVatWnrkkUc8USYAACWGoAj8f76+vpo1a5b8/PzcXnvnnXd022235bsI99GjRzV8+HA98cQTxtfffvttVapUyaP1AgBQ3Nz/RQQqsI4dO+rll1/W//3f/7m9tnLlSjVt2lSdOnXS9ddfr4iICJ0/f17bt2/Xt99+q6ysLOOYjzzyiPr371/MlQMA4HkEReAyEydOVGJiol555RW311wulzZu3KiNGzcWaqzhw4fr9ddf93SJAACUCC49AwZ///vf9e677yowMNBWf19fXz333HOaN2+efHz4NgMAlE38CwbkY9y4cdq1a5fuvPNOS2GvV69e2rx5s6ZMmSKHw1GMFQIAULy49AxcQaNGjbR06VLFx8frP//5j9auXas9e/bo+PHjSktLU1BQkK666io1a9ZMN954owYMGKBmzZp5u2wAADyCoAgUQv369fXoo4/q0Ucf9XYpAACUGC49AwAAwIigCAAAACOCIgAAAIwIigAAADAiKAIAAMCIoAgAAAAjgiIAAACMCIoAAAAwIigCAADAiKAIAAAAI4IiAAAAjAiKAAAAMCIoAgAAwIigCAAAACOCIgAAAIwIigAAADAiKAIAAMCIoAgAAAAjgiIAAACM/LxdAADvcTrS5eNILfTxidn+luc4me2w3EeSInyCLffJzrE1lRLTrHf0c9h7X1cF+VruE+xnby5fe93k72P97znIGWFrrrNZGZb7ZCjL1lwROVUt9/F3BNiay45zOem2+kVk2/unvGm49XNFPja/pg6nuCz3OZZq7+8522V9rlBHoK25Lrisf/2muzItHe902fs8eApnFAEAAGBEUAQAAIARQREAAABGBEUAAAAYERQBAABgRFAEAACAEUERAAAARgRFAAAAGBEUAQAAYERQBAAAgBFBEQAAAEYERQAAABgRFAEAAGBEUAQAAIARQREAAABGBEUAAAAYERQBAABgRFAEAACAEUERAAAARgRFAAAAGBEUAQAAYERQBAAAgBFBEQAAAEYERQAAABgRFAEAAGBEUAQAAIARQREAAABGBEUAAAAYERQBAABgRFAEAACAEUERAAAARgRFAAAAGBEUAQAAYERQBAAAgBFBEQAAAEYERQAAABgRFAEAAGBEUAQAAIARQREAAABGBEUAAAAYERQBAABgRFAEAACAEUERAAAARn7eLgCA96Q5UpTjKPyPgTRHSjFWk9dxG338swJtzRWSE2q5T7CCbM2VmRNiuU/jcH9bczWqbKubLmRZP4eQlG7vvIN/isNyn+POVFtz5chluU8NVbM1V2Vf639nwX72PodVA+3187H+qVdmjq2plJZl/XN/ISvT1lwXlG65z3nHOVtzZfikWe6TpQxLxztd5y3P4UmcUQQAAIARQREAAABGBEUAAAAYERQBAABgRFAEAACAEUERAAAARgRFAAAAGBEUAQAAYERQBAAAgBFBEQAAAEYERQAAABgRFAEAAGBEUAQAAIARQREAAABGBEUAAAAYERQBAABgRFAEAACAEUERAAAARgRFAAAAGBEUAQAAYERQBAAAgBFBEQAAAEYERQAAABgRFAEAAGBEUAQAAIARQREAAABGBEUAAAAYERQBAABgRFAEAACAEUERAAAARgRFAAAAGBEUAQAAYERQBAAAgBFBEQAAAEYERQAAABgRFAEAAGBEUAQAAIARQREAAABGBEUAAAAYERQBAABgRFAEAACAEUERAAAARgRFAAAAGPl5uwAA3uN0pUou32KdI1uZtvo5bPweG+6obmsuO0Idgbb6NQ73t9yndZUsW3NdHei01c+ZY/1zf8DH3ufjnNP6XJWyAmzNFexyWe4TGWj97+v3fta/r+pUsl6fJDWqZO/vOdPlsNxn1zl7nw9njvX3liF7X/dJPsct90l1nbE1V3qW9X45yrF0fFZ2muU5PIkzigAAADAiKAIAAMCIoAgAAAAjgiIAAACMCIoAAAAwIigCAADAiKAIAAAAI4IiAAAAjAiKAAAAMCIoAgAAwIigCAAAACOCIgAAAIwIigAAADAiKAIAAMCIoAgLine truncated "image/png": "iVBORw0KGgoAAAANSUhEUgAAAooAAATCCAYAAAA6r8e3AAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAuIwAALiMBeKU/dgAAvPdJREFUeJzs3Xd4FUX//vH7pBBCCYQqTXpVQEABla4UpSNNARXRrwWRB9ujWBDsXVGagIiABaSKAjaKgoJUKVJCL4be0ss5vz/8JQ/JTpKzm3JS3q/r4pLM7sx8TijezO7Oujwej0cAAABAKn6+LgAAAAC5E0ERAAAARgRFAAAAGBEUAQAAYERQBAAAgBFBEQAAAEYERQAAABgRFAEAAGBEUAQAAIARQREAAABGBEUAAAAYERQBAABgRFAEAACAEUERAAAARgRFAAAAGBEUAQAAYERQBAAAgBFBEQAAAEYERQAAABgRFAEAAGBEUAQAAIARQREAAABGBEUAAAAYERQBAABgRFAEAACAEUERAAAARgRFAAAAGAX4ugAgr0hMTNT27du1Y8cOnT17VhERESpRooTKli2rpk2bqnbt2r4uEQCALEVQBDLw559/asKECVqwYIEuX76c5nnlypXTkCFD9PDDD6tmzZo5WCEAANnD5fF4PL4uAsiNTp8+rZEjR+rLL7+01S8gIEBPPPGExowZo+Dg4GyqDgCA7EdQBAy2bt2q2267TeHh4Y7HaNSokZYvX64KFSpkYWUAAOQcgiKQyvr169WpUyddunQp02NVq1ZNa9asUZUqVbKgMgAAchZPPQNXOHHihHr37m0MiUFBQXrwwQf13Xff6eDBgzp37pz27dunr776SnfccYdcLpelz6FDh9S7d29FR0fnRPkAAGQpVhSBK9xyyy365ZdfLO0333yzvvzyy3RXBjdu3KiePXvqxIkTlmOPPvqoPvrooyytFQCA7MaKIvD/ff3118aQ2LlzZ61cuTLDy8fXX3+97r//fuOxSZMm6a+//sqSOgEAyCmsKAKS3G63atWqpYMHD6Zor1GjhrZt26ZixYplOEZsbKxq1qyp48ePG4/ffvvt+u6777KkXgAAcgIrioCkxYsXW0KiJE2ZMsWrkChJU6dOTTMkStKyZcu0d+9exzUCAJDTCIqApIkTJ1ra2rRpo1tvvdWr/rGxsXrjjTfSPcfj8Wjy5MmO6gMAwBcIiijwzp8/r1WrVlnaR44c6fUYqVcTO3bsqBdffNFy3sKFCx3VCACAL3CPIgq82bNna8iQISnaSpQooZMnTyooKCjD/qZ7E9euXasyZcqobt26lvP/+usvNWzYMPOFAwCQzVhRRIFnWk289dZbvQqJknU18dZbb9VNN92kOnXqqEaNGpbzV69e7bhWAAByEkERBd6mTZssbe3atfOqr+nexJdeeindcUzzAQCQGxEUUaDFxcVp586dlvamTZt61d+0mnjzzTenOw5BEQCQVxAUUaAdOXJE8fHxlvYGDRpk2Dej1cS0xjlw4IC9IgEA8BGCIgq0o0ePWtpCQ0NVsmTJDPtmtJooyXiPYmRkpC5evGi/WAAAchhBEQWaaYPs8uXLZ9jPtJo4ZswYr8dKb2NuAAByC4IiCjTTyl6ZMmUy7Jd6NfGWW25Rq1atLOcVLlxYRYsW9WpeAABymwBfFwD4UkxMjKUtODg43T7e3Jt4pcKFCysyMjJFW3R0tPdFpuPUqVM6ffq0rT6XLl3Sxo0bFRISopIlS6pKlSpebwUEAMh5sbGxKW6Vatu2rVe3SGUFgiIKNFNgyyg0TZs2zavVxCSFCxe2tJkCqhMTJ07U2LFjs2QsAEDesGjRIvXs2TNH5uLSM5CKy+VK85jd1cS0xuOFSACAvICgiALNtNoXGxub5vnTpk3TsWPHkr/u0KFDuquJknn10DQvAAC5DZeeUaCZ7kdM67Kwk9XEtMbL6D5Ibz3yyCPq16+frT67du1S//79k7+uHXy7CvuX8Lp/eb8QW/NJUpFAZ/8m9cj+ymuRAGdzlSqU9kpyWtyOZpJiE+33iUpwtgodk+is36X4BNt9zrojHM3ldtn/hgR4Ah3NddnP/oNklz2nHM3lr0K2+9R3ZbyHq0nZwjn3v/MQB39WJCnAQbc4h3/IzsbY73gsNjLjk7JIYZe9379Riee1JWpB8tdVqlTJ6pLSRFBEgRYSYg09Z86cMZ5rWk1s3bp1uuPHxsYqIsL6P8/ixYvbrNSsXLlyKleuXKbGKOxfQkX8S3t9foh/SdtzFAvwt91HchYUizkMpeWC7P9fzGEGU4yDoBjhcjZZlMN+Hrd1I/qMRLucPRSV6LIfSgM99kOYJMX72f91jvE4CxABsn/lIMRV1tFcpQKdBWcnQh0GRSd/NJ38o0qS3In2g+L5hMvOJnOgiMvZ798kOfkAIpeeUaBVrFjR0hYeHm5pc7qaePLkSWN7pUqVvCsQAAAfIiiiQKtcubKl7dy5c7p06VKKNieriZJ08OBBS1twcLBKlSrloFoAAHIWQREF2tVXXy1/f+tl0b///jv5505XE6V/7wdMrVq1arZqBADAVwiKKNCCg4NVr149S/uWLVuSf556NbF9+/ZerSZK0tatWy1tTZo0sV8oAAA+QFBEgdesWTNL26pVqyRlbjVRklauXOnVfAAA5EYERRR4bdq0sbT99NNPio+PN64mms432b9/v/bt22dp93Y1EgAAXyMoosDr1q2b5e0pZ8+e1eLFizO1mjhnzhxLW4UKFXT99dc7qhMAgJxGUESBV758ed14442W9n79+jleTYyNjdWUKVMs7T169Ej3FYEAAOQmBEVA0kMPPZThOWPGjPF6vEmTJunEiROO5gEAILcgKAKSBgwYoAoVKqR5vF27dmrbtq1XYx04cEAvvPCCpb1Nmza67rrrnJYIAECOIygCkgoVKqSxY8emefz555/3apyzZ8+qa9euxtf2vfbaa47rAwDAFwiKwP83bNgw3XDDDcZjr776qvFS8pW2bNmili1bavfu3ZZjQ4YM0c0335wldQIAkFMIisD/5+fnpy+//FKhoaGWYytXrlStWrX0yCOPaPny5Tp8+LDOnz+v/fv3a968eerfv7+uv/56hYWFWfrWrVtXH330UU58BAAAslSArwsAcpOaNWtq3rx56tatm2JiYlIci46O1qRJkzRp0iSvxytXrpwWL16sEiVKZHWpAABkO1YUgVRuueUW/fTTT8aVRTtq1KihtWvXqm7dullUGQAAOYugCBjcfPPN+uuvv9SjRw/bfV0ulx588EFt3rxZtWrVyobqAADIGQRFIA2VK1fW4sWL9euvv+rOO+9U0aJF0z2/VKlSeuSRR7R9+3ZNnjyZy80AgDyPexSBDLRq1UqtWrVSQkKCtm7dql27duncuXOKjIxUSEiIypQpo8aNG6t+/fq8dQUAkK8QFAEvBQQE6Prrr+ddzQCAAoNLzwAAADAiKAIAAMCIoAgAAAAjgiIAAACMCIoAAAAwIigCAADAiKAIAAAAI4IiAAAAjAiKAAAAMCIoAgAAwIigCAAAACOCIgAAAIwIigAAADAiKAIAAMCIoAgAAAAjgiIAAACMCIoAAAAwIigCAADAiKAIAAAAI4IiAAAAjAiKAAAAMArwdQEAfCvOFSM/V5TX559KDLQ9x5lEl+0+khTqF2y7T6Lb0VQ6FW2/Y4DL2ecqU9jfdp/gAGdz+TvrpkA/+7/OheNCHc11ISHWdp9YJTiaK9Rd2nafQFchR3M5cdEd46hfaKKz/53XLWF/vcjP4e+po5Ee233+iXL265zosT9XMVeQo7kiPPZ//8Z44m2dH+dx9n3ICqwoAgAAwIigCAAAACOCIgAAAIwIigAAADAiKAIAAMCIoAgAAAAjgiIAAACMCIoAAAAwIigCAADAiKAIAAAAI4IiAAAAjAiKAAAAMCIoAgAAwIigCAAAACOCIgAAAIwIigAAADAiKAIAAMCIoAgAAAAjgiIAAACMCIoAAAAwIigCAADAiKAIAAAAI4IiAAAAjAiKAAAAMCIoAgAAwIigCAAAACOCIgAAAIwIigAAADAiKAIAAMCIoAgAAAAjgiIAAACMCIoAAAAwIigCAADAiKAIAAAAI4IiAAAAjAiKAAAAMCIoAgAAwIigCAAAACOCIgAAAIwIigAAADAiKAIAAMCIoAgAAAAjgiIAAACMAnxdAADfinZFyu3y/q+CaFdkNlaTUriDPoEJQY7mKuIuZrtPsAo7miveXcR2n9olAh3NVau4o26KSLC/jnA6xtnaQ2Cky3af8LgoR3O55bHdp4LKOZqruL/9X7PgAGffw9JBzvr52f/WK97taCpFJ9j/3kckxDuaK0Ixtvtccl10NFesX7TtPgmKtXV+nOeS7TmyCiuKAAAAMCIoAgAAwIigCAAAACOCIgAAAIwIigAAADAiKAIAAMCIoAgAAAAjgiIAAACMCIoAAAAwIigCAADAiKAIAAAAI4IiAAAAjAiKAAAAMCIoAgAAwIigCAAAACOCIgAAAIwIigAAADAiKAIAAMCIoAgAAAAjgiIAAACMCIoAAAAwIigCAADAiKAIAAAAI4IiAAAAjAiKAAAAMCIoAgAAwIigCAAAACOCIgAAAIwIigAAADAiKAIAAMCIoAgAAAAjgiIAAACMCIoAAAAwIigCAADAiKAIAAAAI4IiAAAAjAiKAAAAMCIoAgAAwIigCAAAACOCIgAAAIwIigAAADAiKAIAAMCIoAgAAACjAF8XAMC34jxRksc/W+dIVLyjfi4H/5Yt4brK0VxOFHMFOepXu0Sg7T6NSyY4mqtsUJyjfnFu+9/7MD9n34+LcfbnKppQyNFcwR6P7T6lguz/ev3bz/6fqypF7dcnSbWKOvt1jve4bPfZcdHZ9yPObf+zxcrZ7/vTfuG2+0R5zjuaKybBfj+33LbOT0iMtj1HVmFFEQAAAEYERQAAABgRFAEAAGBLine truncated
"text/plain": [ "text/plain": [
"<Figure size 600x1200 with 1 Axes>" "<Figure size 600x1200 with 1 Axes>"
] ]
@@ -1120,13 +1133,13 @@
}, },
{ {
"cell_type": "code", "cell_type": "code",
"execution_count": 26, "execution_count": 25,
"id": "b5a59943", "id": "b5a59943",
"metadata": {}, "metadata": {},
"outputs": [ "outputs": [
{ {
"data": { "data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAiwAAAGiCAYAAADEJZ3cAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAPYQAAD2EBqD+naQAARzxJREFUeJzt3Xtc1FX+x/HXMCimAVFkZoOYgpey1J9mN81balm5bbXlpTY3DLP73TUzMVuttu2qpiSaZWZW2u7aZa0EjcqyzLKLChbmaKWUAmmCMPP74wiKgs4gcObyfj4e8+DLd8bxjTjMh/M953McXq/Xi4iIiEgAi7AdQERERORwVLCIiIhIwFPBIiIiIgFPBYuIiIgEPBUsIiIiEvBUsIiIiEjAU8EiIiIiAU8Fi4iIiAS8SNsBaovH42HLli1ER0fjcDhsxxEREREfeL1eioqKaN68ORER1Y+jhEzBsmXLFhISEmzHEBERkRrYtGkTLper2vtDpmCJjo4GzBccExNjOY2IiIj4orCwkISEhIr38eqETMFSfhkoJiZGBYuIiEiQOdx0Dk26FRERkYCngkVEREQCngoWERERCXgqWERERCTgqWARERGRgKeCRURERAKeChYREREJeCpYREREJOCpYBEREZGAp4JFREREAp4KFhGxwu12k5mZidvtth1FRIKAChYRqXcZGRkkJibSp08fEhMTycjIsB1JRAKcw+v1em2HqA2FhYXExsZSUFCgzQ9FApjb7SYxMRGPx1Nxzul0kpeXd8it5UUkNPn6/q0RFhGpVzk5OZWKFYCysjJyc3MtJRKRYKCCRUTqVXJyMhERlX/0OJ1OkpKSLCUSkWCggkVE6pXL5SI9PR2n0wmYYmXGjBm6HCQih6Q5LCJihdvtJjc3l6SkJBUrImHM1/fvyHrMJCJSweVyqVAREZ/pkpCIiIgEPBUsIiIiEvBUsIhIyFD3XJHQpTksInJYbrebnJwckpOT62feiccDP/0EmzfDr79Cfj6UlEBKyr7HPPYYfPGFOV9WxvduNx+uXEkxsA5okJ5OyvXXm8d+9hns2AHHHw8tWsAxx4DDUfdfh4jUGq0SEpFDysjIIDU1FY/HQ0REBOnp6aTsXzgciR07TPFQbtQoeO89+PFHU4jsr0kT+P33fZ8PHAhvv13l03qAhhER5G3caAqsv/wFXntt3wNiYqBlS0hMhFat4MEHzTkRqXdaJSQiR8ztdlcUKwAej4eRI0cyYMAA/0datm+Hjz6C7GxYuRK++gr27IHffts32rFlC5R3vHU64cQTzajIcceZW1mZOQ8wYgT06wcNG7IuN5cpTz5JQ6Ah5gdbmcdDbm6uyZmQAB06wC+/wLZtUFho/v6vvoKGDeGf/9yX8+9/h/Xr4Zxz4NxzoUsX8xgRsUojLCJSrczMTPr06VPl+V69evn2JE88AbNnw9dfw4E/bhwOcLuheXPz+aefwh9/mNGPk06CSN9+p/Jrf6Jdu2Djxn23oiK4++599//f/5lLTeUaNYIzzoAePcyozjnn6HKSSC3SCIuIHLHyNvoHFgLVttHfsAEWL4bUVDjqKHPul19gzRpz3KYNdO8OZ58NHTvCqadC48b7/ny3bjXKWd49d+TIkZSVlR26e27jxtC+vblV5amnYMUK+PBDc8vPhw8+MLe5cyEvb99jd+82BY2I1DmNsIjIIWVkZBxUCFSaw7JlC8ybBy++aC6xALz5phmNADOysn69ubxywgl1mrXWu+d6vSZ7dja8/74Z+Zk0ydxXWmpGgbp1g6uvhkGD9hVpIuIzX9+/VbCIyGEdVAjs3AkLF5oi5f33zaoeMJdwevSAsWOhb1+7oevaZ5+ZS0XloqPh8svhr3+FXr102UjERypYRKTufPopnHnmvs/PPReuucasxjn2WHu56tt338FLL5lLRRs37jvfvj1MmQJVzP8RkcpUsIhI7fB6zeqetWsr90EZOhTatjWXQ1q3tpcvEHg85t/oxRfh5ZfNRN7PPjMrjMD8G2rERaRKKlhEwlCtNnjzemHpUhg3Dj7+2ExW3bQpvEZQaqKw0MzhGTJk37lRo6C42PxbnnyyvWwiAcjX92+15hcJERkZGSQmJtKnTx8SExPJyMio+ZNlZ5vLGeefb4qVqCgzolJcXHuBQ1VMTOVipbAQ5swxS7vbtIEbbjCFn4j4RQWLSAiorsGb33vqrF0LF15oJs5mZZmGabfeapbyPvecaeQm/omJgcxM6N/frCyaMQOSksy/a36+7XQiQUMFi0gIyMnJqdQrBaCsrIzc8q6xvoqKMoVKZKTppZKba/qSNGtWe2HD0Zlnwv/+B8uXQ8+eZtuBZ54xc4Def992OpGgoIJFJASUN3jb3yEbvJXzeMybaLmTTzarXtauNSMBCQl1kDb41Nou0D16mNGW996D0083WxOcckrthBQJcSpYREJAeadX5959dg7Z6bXcF1+YNvM9e5p5KuUuu0yrfvZTq3ODwKwW6tsXPv/cFIv7X2abNs3srSQiB9EqIZEQ4lOn19JSePhhmDDBHB99tBlNGTq0fsMGAb/2KDpSS5bAgAGmgJk1Cy64oHafXyRAaZWQSBhyuVz06tWr+jfTnBxzWWLcOFOsXH45rFunYqUatTY3yBdxcWYV0U8/mYnPN95oOgqLCKCCRSR8PP88dOpkNvaLiTFNzl59dd9OyXKQGs8NqokzzjCX6W65xXz+7LPm+7X/5TqRMKaCRSSc7NoFvXub3ZOvvlrdVw+jRnODjkTjxvD00/Duu+BymVVa3bublVoiYU5zWERCWVkZ7H2zxeuFf//b7Cocod9V/FHru0D7YscOuOkmsxP2a6+Zy3ciIUit+UXC3auvwoMPwrJlaqcfrLxe+OQTOOusfedKS02fHJEQoUm3IuGqtBRGj4Yrr4Svv4bHH7edSGrK4ahcrLjdZifoN96wFknEFhUsIqGkqAguvhgefdR8fs89kJZmNZLUon/+08xr+fOfYfJkMwIjEiY0rigSKn75BQYOhFWrzOTN2bPNKIuEjsceM6MuTz0F990Hmzeb4/J5SiIhrEYjLCtXrmTgwIHExcXRpEkTunXrxrx582ocYs+ePXTq1AmHw0G7du1q/DwiYWvDBtO1dtUqOP54sx+QipWAVqN2/w0awJNPmiLF4YCpU+Gqq2D37jrLKRIo/C5YsrKy6N69Ox988AFXXHEFo0aNIj8/n2HDhjFp0qQahZg4cWLdNGISCRdHH23ewFq1go8+Mj09JGAdcbv/W2+F+fPNbtqvv252gt6+vW7CigQIv1YJlZaW0q5dO9xuNx9//DGdO3cGoKioiLPPPpt169bx7bffkpyc7HOAVatWceaZZ/L4449z66230rZtW9auXev3F6JVQhL2vv8emjSBE06wnUQOoVbb/WdlwaWXQmKi2ZcoNrZWs4rUhzpZJbR06VI2bNjA0KFDK4oVgOjoaMaNG0dpaSmzZ8/2+flKSkoYPnw4Z511FjfffLM/UUTkrbfglVf2fd6qlYqVIFCr7f579YIPPoC331axIiHPr0m3WVlZAPTv3/+g+8rPLVu2zOfnS0tLIycnhy+//BKHOm6K+O6dd8xKkdJSaNbM7LgsQaG83f+BIyw1bvd/2mmVP3/5ZfP/QVsuSIjxa4QlJycHoMpLPnFxccTHx1c85nBWrlzJo48+yoQJE2jTpo0/MQAoLi6msLCw0k0kLCxZYi4DlJSYj+ecYzuR+KFO2/2/+KLZyLJ3b9iy5cifTySA+FWwFBQUABBbzdBjTExMxWMOpbi4mOHDh9O5c2fuuusufyJUmDx5MrGxsRW3hISEGj2PSFB57z3405+guNgUKy+/bFaOSFBJSUkhLy+PzMxM8vLySElJqZ0n7tEDWrSA9euhTx+z87NIiLDSOG7cuHHk5OQwa9asit8y/DVmzBgKCgoqbps2barllCIBZulSuOQSs4R10CAzf6VhQ9uppIZcLhe9evWq3b2JWrY0E3FbtIB160zR8ssvtff8Ihb5VbCUj6xUN4pSPtP3UFatWsXjjz/O2LFjOe3Aa69+iIqKIiYmptJNJGStX2+KlN27TSfbV19VsSJVO/lkyMyEhARYu9Y0Eywqsp1K5Ij5VbCUz12pap7K9u3byc/PP+yS5q+++oqysjLS0tJwOByVbgDr1q3D4XBwzDHH+BNNJLQlJcGIEdCvn9m5V8WKVMPtdpO5cSM/vfSSaSK4ahVccQXs2WM7msgR8WuVUM+ePZk8eTJLlixh8ODBle5bsmRJxWMOpU2bNtVer83IyCA2NpYrrriCxo0b+xNNJLRFRJgOpyUlKlakWhkZGaSmpuLxeIiIiGDh3//On558Erp21Q7PEvT8bhzXtm1bNm/ezIoVK+jUqRNQuXHcN998U7HqJz8/n/z8fOLj44mPjz98GIdDjeNEyu3aZfaOGT0aoqJsp5EAV11Duk0ffsiJZ55pMZnIodVJ47jIyEhmzpyJx+OhR48epKamcvfdd9OxY0e++eYb0tLSKi1RnjJlCu3bt2fKlCk1/0pEwlFpKQwZAuPHw9VX204jQaC6hnTr/vhj34k//jANB0WCkN+rhHr37k12djbdu3dnwYIFTJs2jeOOO465c+cyduzYusgoEl68XrNXzH/+A40awe23204kQaC8Id3+KjWk27XLzIG6+GL4738tJBQ5Mn5dEgpkuiQkIWP6dBg1ymxm+NprcNllthNJkMjIyGDkyJGUlZVVNKSrmDPo9cINN0B6OkRHw6efQrt2dgOL4Pv7twoWkUCSnW16Z+zZAw8/bOaviPjB7XaTm5tLUlLSwT1e9uyB8883GyW2aWOKFu1BJJapYBEJNm63Wc3xyy9w5ZUwf74ZZRGpTVu3mv9nmzbBRReZS48RVnqIigB1NOlWROrQ99+bxnCnnQazZqlYkbrRtCksWmTmR735ppnYLRIEVLCIBIrzzoOVK+GNN6BJE9tLine truncated "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiMAAAGdCAYAAADAAnMpAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAPYQAAD2EBqD+naQAARJhJREFUeJzt3Xlc1NX+x/HXMCEuCGqlZoOYgZmV2ZUybbmGbVaaqb8o08rMtc28lZrdsm5ZlvealqYUmS1mi+ntappadEvLFrfMXJi8mKOhpggEhgrz++MIuIACDpxZ3s/HYx7znWGAt474/XC+53yOw+v1ehERERGxJMx2ABEREQltKkZERETEKhUjIiIiYpWKEREREbFKxYiIiIhYpWJERERErFIxIiIiIlapGBERERGrTrIdoDwKCwvZtm0bdevWxeFw2I4jIiIi5eD1esnJyaFJkyaEhZU9/hEQxci2bduIiYmxHUNEREQqYcuWLbhcrjI/HhDFSN26dQHzh4mKirKcRkRERMojOzubmJiY4vN4WQKiGCm6NBMVFaViREREJMAcb4qFJrCKiIiIVSpGRERExCoVIyIiImKVihERERGxSsWIiIiIWKViRERERKxSMSIiIiJWqRgRERERq1SMiIiIiFUqRkRERMQqFSMi4nMej4fU1FQ8Ho/tKCISAFSMiIhPpaSkEBsbS2JiIrGxsaSkpNiOJCJ+zuH1er22QxxPdnY20dHRZGVlaaM8ET/m8XiIjY2lsLCw+Dmn00l6evoxtw8XkeBU3vO3RkZExGfS0tIOK0QACgoKcLvdlhKJSCBQMSIiPhMfH09Y2OH/rTidTuLi4iwlEpFAoGJERHzG5XKRnJyM0+kETCEydepUXaIRkWPSnBER8TmPx4Pb7SYuLk6FiEgIK+/5+6RqzCQiIcLlcqkIEZFy02UaERERsUrFiIiIiFilYkREAoK6uooEL80ZEQlxHo+HtLQ04uPjq2eeR2Eh/PYbbN0Ku3bB77/Dvn3Qr1/Ja8aNg5UrzfMFBWzyeFj6/ffkAxuA8ORk+vXvb177ww+wZw+ceio0bQr16oHDUfV/DhHxGa2mEQlhKSkpDBgwgMLCQsLCwkhOTqbfoUXBidizxxQGRQYPhsWL4ddfTZFxqDp14I8/Sh5fdx3Mn1/qly0EaoSFkb55syme/u//4MMPS14QFQXNmkFsLDRvDk89ZZ4TkWqn1TQickwej6e4EAEoLCxk4MCBXHPNNRUfIcnMhK+/hiVL4Pvv4ccfYf9+2L27ZJRi2zYo6sTqdMJpp5nRjJNPNreCAvM8wN13w1VXQY0abHC7efnFF6kB1MD8p1VQWIjb7TY5Y2Lg3HNh+3bYuROys833//FHqFEDXnihJOeIEbBxI3ToAJdcAm3bmteIiFUaGREJUampqSQmJpb6fMeOHcv3RcaPh2nT4Kef4Mj/ShwO8HigSRPz+LvvYO9eM2px+ulwUvl+F6rQfjd5ebB5c8ktJwceeqjk43/5i7n8U6RmTbjwQrjsMjMa06GDLvGI+JBGRkTkmIpatx95ki+zdfsvv8DcuTBgANSqZZ7bvh3WrDHHLVrApZdC+/Zw/vlwzjlQu3bJ5190UaVyFnV1HThwIAUFBcfu6lq7Npx9trmVZsIEWLYMli41t99/h6++Mre334b09JLX/vmnKVZEpMppZEQkhKWkpBx1kj9szsi2bTBjBrz1lrnsATBvnhlFADMisnGjueTRqFGVZvV5V1ev12RfsgQ++8yM2IwZYz524IAZvbnoIujdG7p2LSnARKTcynv+VjEiEuKOOsnn5sJHH5kC5LPPzOoXMJdVLrsMRo2CTp3shq5qP/xgLt8UqVsXevSA22+Hjh11KUeknFSMiEjlfPcdtGtX8viSS6BPH7NqpUEDe7mq27p18M475vLN5s0lz599Nrz8MpQy30ZEDqdiRESOz+s1q2DWrz+8z0evXnDWWeYSxZln2svnDwoLzd/RW2/Bu++aSbE//GBW4oD5O9RIiUipVIyIBBmfNifzeuHzz+Hvf4dvvjETP7dsCa2Rj8rIzjZzZm69teS5wYMhP9/8XZ5xhr1sIn6ovOdvtYMXCQApKSnExsaSmJhIbGwsKSkplf9iS5aYSwxXXmkKkYgIMxKSn++7wMEqKurwQiQ7G6ZPN8ubW7SAQYNMUSciFaJiRMTPldWcrMJ7tKxfD507m0moX3xhmn3df79Zzvrqq6YJmVRMVBSkpsLVV5sVOFOnQlyc+Xv9/Xfb6UQChooRET+XlpZ2WC8QgIKCAtxF3UzLKyLCFCEnnWR6hbjdpu9G48a+CxuK2rWDTz+FL7+Ev/7VtLp/6SUz5+azz2ynEwkIKkZE/FxRc7JDHbM5WZHCQnOCLHLGGWZ1yPr15jf4mJgqSBt4fLYb8GWXmVGSxYuhdWvTDr9VK9+EFAlyKkZE/FxRB1LnwX1bjtmBtMjKlaa1+V//auaFFOneXatjDuHTuThgVtV06gTLl5tC8NBLX5Mnm716ROQoWk0jEiDK1YH0wAF47jl48klzHBlpRkF69aresAGgQnvenKiFC+Gaa0xx8vrrcO21vv36In5Kq2lEgozL5aJjx45lnyjT0sylgr//3RQiPXrAhg0qRMrgs7k45VG/vllt89tvZhLxkCGm062IACpGRILDG29AmzZmE7ioKNOg64MPSnbMlaNUei5OZVx4obl0dt995vErr5j369BLaCIhTMWISLDIy4MrrjC76Pbura6gx1GpuTgnonZtmDgRFi0Cl8usZrr0UrOiSSTEac6ISKAqKICDJ1K8Xvj3v83usmH6HaMifL4bcHns2QP33GN2RP7wQ3NJTSQIqR28SDD74AN46in473/Vwj1Qeb3w7bdw8cUlzx04YPrAiAQJTWAVCUYHDsDw4XDzzfDTT/Cvf9lOJJXlcBxeiHg8ZkfgOXOsRRKxRcWISKDIyYEbboDnnzePH34YRo+2Gkl86IUXzDySm26CZ581IyciIULjgSKBYPt2uO46WLHCTIScNs2MjkjwGDfOjJZMmACPPgpbt5rjonlBIkFMIyMi/u6XX0w31RUr4NRTzf4yKkT8WqVazIeHw4svmgLE4YBJkyApCf78s8pyiviLShUjCxYsoHXr1tSpU4cOHTqwfPnycn1ednY2DoeDevXqVebbioSmyEhzcmreHL7+2vSsEL91wi3m778fZs40uyrPmmV2BM7MrJqwIn6iwsXIypUr6datG7169WLp0qWEh4eTkJBARkbGcT937NixXHDBBTTQ7H+R8mvUyLQT//prsz29+C2Px8OAAQOKO7sWFhYycODAim/Cd/PNZifg6GjIytJybQl6Ff4XPnbsWLp3786IESNo06YNDRs2BGDSpEnH/Lxt27Yxe/ZshgwZclTXQxE5wiefwHvvlTxu3twUJeLXfNpivmNH+OormD/fFCUiQazCVcHixYtJSkoCwO12s3btWsaMGcOiRYuO+XlPPPEE9957LydpDb3IsS1YYFZU9Opl+ohIwPB5i/nzzju8pf+778K2bSeQUMQ/VagYycrKYteuXbRo0QKAcePGMWzYMFq1asWmTZvK/Lx169Yxf/587rrrrnJ9n/z8fLKzsw+7iYSEhQuhWzfYt8/cd+hgO5FUQJW2mH/rLVOgXnGFChIJOhUqRnJycgCoXbs2O3bsYOHChfTp04fIyMjij5VmxIgRPPDAA9SsWbNc3+fZZ58lOjq6+BYTE1ORmCKBafFiuPFGyM83hci775oVFhJQ+vXrR3p6OqmpqaSnp9OvXz/ffOHLLoOmTWHjRkhMNDsAiwSJChUjkZGRAOTl5TFhwgQGDRpEREQEubm5xR870pIlS/jqq68YNGhQub/PyJEjycrKKr5t2bKlIjFFAs/nn0OXLmYZZ9euZr5IjRq2U0kluVwuOnbs6Nu9bpo1M8u6mzaFDRtMQbJ9u+++vohFFZrAER0dTYMGDVixYgUzZsxg9erVAKxfv57mzZuX+jmPPPII999/P3Xr1i3394mIiCAiIqIi0UQC18aNpgD580/TYfWDD1SISOnOOANSU83k1vXrTSO8L76ACvz/KuKPKjQy4nA4SExMpHfv3iQlJRVvejNr1iw6dep01Ot//fVXvvnmG8aMGUPNmjWpWbMmffv25ZdffqFmzZp8/fXXvvlTiASyuDi4+2646iqzg6sKESmDx+MhdfNmfnvnHdMAb8UK6NkT9u+3HU3khFR4Nc3w4cMBM8l09erVDB06lA0bNjBkyBAAunTpwj333APAaaedxrp16/jxxx9ZtWoVq1atKv78VatWccEFF/jqzyESuMLCTOfNuXNBI4JShkObqbk6duTf/fubrQESErTTrwS8Cv8LTkhIYO7cuYwYMYKpU6fSunVrFixYUHxtNCMjg1q1agEQHh5Oy5YtD/v8li1bcuaZZx71vEhIycsze5EMH15SgGhERMpQWjO1HmPHsmXpUk5r185yOpET5/B6/X9ryOzsbKKjo8nKyiq+NCQSsA4cgB494OOPzRD7Bx/YTiR+LjU1lcTExFKf79ixo3mwd6+ZT3LdddUbTuQYynv+VitUkerk9Zq9Rz7+GGrWhKFDbSeSAHDcZmp5eWbO0Q03wH/+YyGhyIlRMSJSnaZOhVdeMRvfvfMOXHKJ7UQSAI7bTK1WLTjnHFPs3nabWWkjEkB0mUakuixZYnpD7N8Pzz1n5ouIVIDH48HtdhMXF3d0D5P9++HKK+HLL6FFC/juO+1pI9aV9/ytYkSkOng8ZtXD9u1mR9aZM83oiIgv7dhh/p1t2QLXX28uB2pjUrFIc0ZE/MmmTaap2XnnweuvqxCRqtGwIcyebeYjzZsHTzxhO5FIuagYEakOl18O338Pc+ZAnTq200gwa9sWXnvNHE+eDL//bjePSDmoU45IVcrLM42pAOLj7WaR0HHbbbBzp9nv6JRTbKcROS6NjIhUleXL4cwzzW6Line truncated
"text/plain": [ "text/plain": [
"<Figure size 640x480 with 1 Axes>" "<Figure size 640x480 with 1 Axes>"
] ]
@@ -1180,13 +1193,13 @@
}, },
{ {
"cell_type": "code", "cell_type": "code",
"execution_count": 27, "execution_count": 26,
"id": "242dd227", "id": "242dd227",
"metadata": {}, "metadata": {},
"outputs": [ "outputs": [
{ {
"data": { "data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAooAAAS+CAYAAABGTwoNAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAuIwAALiMBeKU/dgABAABJREFUeJzs3Xd8VFXeBvBnkkmZ9IQSklAl9Ca9hiC9Coj0KqMEEVkja1lFBEHXdXVBFAklSJMOAoqBAEIMTTqEEoqQBEhCGqROymTm/YM3kTD3ps7cac/38+HzMuecOfeXd3fh4dx7z5FptVotiIiIiIieY2PsAoiIiIjINDEoEhEREZEgBkUiIiIiEsSgSERERESCGBSJiIiISBCDIhEREREJYlAkIiIiIkEMikREREQkiEGRiIiIiAQxKBIRERGRIAZFIiIiIhLEoEhEREREghgUiYiIiEgQgyIRERERCWJQJCIiIiJBDIpEREREJIhBkYiIiIgEMSgSERERkSAGRSIiIiISxKBIRERERIIYFImIiIhIEIMiEREREQliUCQiIiIiQQyKRERERCSIQZGIiIiIBDEoEhEREZEgBkUiIiIiEsSgSERERESCGBSJiIiISBCDIhEREREJkhu7ACJzodVqcf36dVy7dg2PHj1CdnY2FAoFatSogWbNmqFVq1aQy/k/KSIishz8W42oDNevX8eyZcvw888/IykpSXScu7s7hg0bhrfffhudOnWSsEIiIiLDkGm1Wq2xiyAyRRkZGfjggw+watUqaDSaCn13zJgxWLZsGby9vQ1UHRERkeExKBIJuHv3LoYOHYobN25Ueo7atWtj3759aNu2rR4rIyIikg6DItFz7t+/jx49eiAuLq7Kc3l5eeHYsWNo1aqVHiojIiKSFt96JnpGfn4+Ro4cKRgSZTIZxo4di99++w1JSUkoKChAamoqjhw5gtdffx12dnY630lLS8Pw4cORnp4uRflERER6xRVFomd8/PHH+OKLL3Taa9Soge3bt6NXr16i342KisKIESNw9+5dnb4pU6Zg/fr1+iyViIjI4BgUif7f3bt30bx5c+Tl5ZVod3Z2xokTJ9CmTZsy57h//z46deqExMTEEu0ymQwnT55Ely5d9FozERGRIfHWM9H/+/LLL3VCIgB8/fXX5QqJAFCnTh38+OOPOu1arRaLFi2qco1ERERS4ooiEZ4+S+jr66sTFJs3b46oqCjY2FTs31QDBw7EwYMHddqjo6PRpEmTKtVKREQkFa4oEgHYsWOH4Griu+++W+GQCAD//Oc/Bds3bdpU4bmIiIiMhSuKRAD69euHw4cPl2hzcHBAUlIS3NzcKjyfRqNBnTp1EB8fX6K9SZMmiI6OrlKtREREUuGKIlm93NxcHD9+XKe9Z8+elQqJAGBjY4NBgwbptN+8eRP379+v1JxERERSY1Akq3f27Fnk5ubqtJe2FU55vPTSS4LtERERVZqXiIhIKgyKZPUuXLgg2N6+ffsqzduhQwfB9osXL1ZpXiIiIqkwKJLVu3LlimB7ixYtqjSvv78/7O3tddqjoqKqNC8REZFUGBTJ6gmdpKJQKODn51eleW1tbVG/fv1yXY+IiMgUMSiS1YuNjdVp8/HxgUwmq/Lcvr6+Om1xcXHQaDRVnpuIiMjQGBTJ6iUlJem01apVSy9zC81TUFCAJ0+e6GV+IiIiQ5IbuwAiY8rPz0d2drZOu7u7u17mF5snLS0NXl5eVZo7KSkJycnJFfpORkYGzp07Bzc3N3h4eKBOnTpwcHCoUh1ERGQ4eXl5JbZVCwwMhIeHh2TXZ1Akq5aZmSnY7uLiopf5XV1dK3Tdivjhhx+wcOHCKs9DRETmY8+ePRg+fLhk1+OtZ7Jq+fn5gu1CbytXhp2dnWC70HGBREREpoZBkaxaQUGBYLtcrp/FdrGgKHZdIiIiU8Jbz2TVbGyE/62kr7eSxeYRu25FzJo1C6NHj67Qd65fv44xY8YUf96zZw/8/f1LjFHnq/H9q+uQr9INsx1HtcZLM7tXrmAiIqqwO3fuYMSIEcWf69SpI+n1GRTJqondYlar1XqZX2wefdzarlmzJmrWrFmlOfz9/QU3Fu/eryuiDkTrtGdeza/yRuRERFR5Ur+AyFvPZNUcHR0F21UqlV7mz8nJEWxXKBR6md9QmvdpLNieci8NyXdTJa6GiIiMhUGRrJqbm5vg84gZGRl6mV/s7WZPT0+9zG8ozV5qBLH9xq8duSVtMUREZDQMimT1hEJbaqp+Vs1SUlIE26u6h6KhuVZ3Rt0XhY8wjDpwQ+JqiIjIWBgUyeoJnemcmJiol7mF5vHy8jL5W88A0Lyv8O3nmPMPkP6o6vtAEhGR6WNQJKtXv359nbakpCTR5wsrIiYmRqetQYMGVZ5XCq0GNBPtizqo+6ILERFZHgZFsnpNmjTRadNqtbh9+3aV5k1PTxc8Yq9xY+GVOlPj7V8dtRrXEOy78tt1iashIiJjYFAkq9e2bVvB9kuXLlVp3osXL1boeqao1UDhVcW7Z+KQlap7RjYREVkWBkWyep06dRJsP336dJXmPXXqVIWuZ4raDBYOilqNFlfDb0pcDRERSY1BkaxegwYNBJ8bPHToUJXmFfq+k5MTunbtWqV5pVSrSU1UbyD8hvalX69JXA0REUmNQZEIwKBBg3Ta/vrrL9Hbx2VJSkrCH3/8odPet29fvZzKIhWZTIY2g4RXFe+cikFGEt9+JiKyZAyKRAAmTJgg2B4SElKp+VavXo3CwsJyX8eUtRkifGSfVqPFxX1cVSQismQMikQAunXrhubNm+u0r1+/XnCLm9Kkp6fj22+/1WmvWbMmhg8fXtkSjca3uTe8/asL9m37724sXrwYCQkJEldFRERSYFAkwtNbrO+9955Oe15eHmbMmAGtVlvuuebOnSu4Lc6cOXNEz5Y2ZTKZDK2GNBXsU+S54H8Ll6Ju3bqYOXMmcnNzJa6OiIgMiUGR6P9NnjwZbdq00Wk/dOgQ3n777XKFxf/+978IDQ3Vaffz88M777yjjzIlp1KpsGTnV6L9zZxbQ61WY+XKlRg4cCBUKpWE1RERkSExKBL9P1tbW4SEhEAul+v0LV++HIMGDRLdhDs+Ph6TJk3C+++/L9j//fffw9nZWa/1SiU4OBiHTh7Ew9w4wf6mzq2Lfx8REYHg4GCpSiMiIgPT/RuRyIp16dIFX375Jf75z3/q9B08eBBNmjRB165d0bFjR3h6eiIjIwOXLl3CH3/8AbVaLTjnO++8gxEjRhi4csOIj48vXiGNzo6Cn2NdnTEedl7wc6iHh3mxAIDQ0FB8+umn8PHxkbRWIiLSPwZFoufMnTsXSUlJ+Oor3dutWq0WJ0+exMmTJ8s116RJk/DNN9/ou0TJhIaGFgfgWzlX0ctrIGxltjrjWrq0Kw6KarUaoaGhmDdvnqS1EhGR/vHWM5GA//znP/jhhx/g4OBQqe/b2tri008/xYYNG2BjY77/Mzt+/Hjx71WaHNxT3RIc19ipOexlf///KjIy0uC1ERGR4Znv32BEBvbmm2/i6tWrGDVqVIXCXt++fXH27FksWLAAMpnMgBUaXnp6eonPV7MuCI6zs7FHE+eWxZ8zMjIMWhcREUmDt56JSuHv74+dO3ciJiYGP//8M44dO4br168jMTERKpUKjo6OqF69Opo1a4aAgACMHDkSzZoJn2Rijtzd3Ut8vqe6gyx1JlzkrjpjW7q0RVTWeQCAm5ubJPUREZFhMSgSlUP9+vURHBxsdW/09ujRA+Hh4cWftdDgevYldHIP0Bnr41AH1exqILUgGQEBuv1ERGR+eOuZiEQplUqd7YKuZomff93KpQPkcjmUSqWhSyMiIgkwKBKRKF9fX53Q90Sdige5sYLjW7i8COXU17k1DhGRhWBQJKJSLVmyBIGBgSXarv7/s4jPc7BxxISAqVKURUREEmBQJKJSKRQKhIWFISgoqPg29K2ca1AV5giOP7PlcoXOxiYiItPFoEhEZVIoFAgJCUFsbCwWLVqE3v16I80tUXBsQvQjxFx4IHGFRERkCHzrmYjKzdfXt/jElZTYNHzZ63sILR6e2HAWDdrXkbg6IiLSN64oElGlVK/nhSaB/oJ9V367joykLIkrIiIifWNQJKJK6z6lg2B7YYEGJzaelbgaIiLSNwZFIqq0poH+8KrtIdh3ctM55KsKpC2IiIj0ikGRiCrNxtYGPaZ1FOzLeazCuV2XJa6IiIj0iUGRiKqk89h2cHR1EOz7Y+2f0Gi4VQ4RkbliUCSiKnF0dUDncW0F+5LvpuL6kVsSV0RERPrCoEhEVRYwrTNsbGWCfRGrT0lcDRER6QuDIhFVmaefO1oPbi7Yd/dMHGLO35e4IiIi0gcGRSLSi15vdBXtO7z8uISVEBGRvjAoEpFe1Gnti4ad6wn23fj9Nh5cTZC4IiIiqioGRSLSmz5v9RDtO/IDVxWJiMwNgyIR6U3jgBdQp7WvYF9U2A08upMscUVERFQVDIpEpDcymUx0VVGrBX7/4YTEFRERUVUwKBKRXrXo1wS1mtQU7LuwNwopMWkSV0RERJXFoEhEemVjI0OfWcKrippCLcK/jZC4IiIiqiwGRSLSuxeHNkf1+l6CfRf2RCHxNp9VJCIyBwyKRKR3NrY26Pd2gGCfVgscXHJM2oKIiKhSGBSJyCDajWiFmg2rC/Zd+e0G91UkIjIDDIpEZBA2tjYYEBwo2n/gf8ekK4aIiCqFQZGIDKb14ObwbeYt2Hfj99u4xzOgiYhMGoMiERmMjY0MA+f2Eu3/9d+HodVqpSuIiIgqhEGRiAyqeZ/GqNtG+LSWmHP3cfXgTYkrIiKi8mJQJCKDkslkGPReb9H+X/9zGIUFhRJWRERE5cWgSEQG17jHC2gc8IJgX8q9NJzafF7iioiIqDwYFIlIEsM+6geZTLgv/Ns/oMrIlbYgIiIqE4MiEUnCt5k3OoxLine truncated "image/png": "iVBORw0KGgoAAAANSUhEUgAAAooAAAS2CAYAAACqHIhgAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAuIwAALiMBeKU/dgABAABJREFUeJzs3XdYVMe/BvB36UsVRRTsvWvsGEHsPUaNJbZoshrshvQYo7Gk/EzRmKigYtfErokRNSaKosbeO6KgIoKA9Law9w8vRNxzlgV2z7b38zw+152Zc/ZLfrn6ZuacGZlKpVKBiIiIiOglVoYugIiIiIiME4MiEREREQliUCQiIiIiQQyKRERERCSIQZGIiIiIBDEoEhEREZEgBkUiIiIiEsSgSERERESCGBSJiIiISBCDIhEREREJYlAkIiIiIkEMikREREQkiEGRiIiIiAQxKBIRERGRIAZFIiIiIhLEoEhEREREghgUiYiIiEgQgyIRERERCWJQJCIiIiJBDIpEREREJIhBkYiIiIgEMSgSERERkSAGRSIiIiISxKBIRERERIIYFImIiIhIEIMiEREREQliUCQiIiIiQTaGLoDIVOTl5eHKlSu4evUqEhISkJaWBjc3N1SsWBGtWrVCvXr1DF0iERGRTjEoEhXjzJkzWLp0KXbu3InU1FTRcZ6enhgzZgwmTZqEOnXqSFghERGRfshUKpXK0EUQGaP4+HjMmDEDv/76a4mus7GxwQcffIA5c+ZALpfrqToiIiL9Y1AkEnDx4kX06dMHsbGxpb5H8+bNsX//fnh5eemwMiIiIukwKBK95NSpU+jZsydSUlLKfK+aNWvi6NGjqFatmg4qIyIikhbfeiZ6QUxMDAYNGiQYEu3t7REQEIA///wT9+7dQ2JiIu7cuYPffvsNb7zxBmQymdo19+/fx6BBg5CZmSlF+URERDrFGUWiF3Tr1g3//POPWnvHjh3x66+/apwZPHv2LF5//XXExMSo9U2dOhU///yzTmslIiLSN84oEv2/LVu2CIbEXr164fDhw8UuH7dp0wbjx48X7Fu+fDkuX76skzqJiIikwhlFIgD5+fmoW7cu7t27V6S9du3auHTpEpydnYu9R3Z2NurUqYNHjx4J9vft2xd//vmnTuolIiKSAmcUiQDs2bNHLSQCQHBwsFYhEQBWrlwpGhIBIDQ0FLdv3y51jURERFJjUCQCsGzZMrW2Tp06oXv37lpdn52djW+//VbjGJVKhaCgoFLVR0REZAgMimTxkpKScOTIEbX2GTNmaH2Pl2cTe/TogdmzZ6uN27VrV6lqJCIiMgQ+o0gWb+PGjRgzZkyRNjc3Nzx58gT29vbFXi/0bOLx48fh4eGBBg0aqI2/fPkymjVrVvbCiYiI9IwzimTxhGYTu3fvrlVIBNRnE7t3745XX30V9evXR+3atdXGh4WFlbpWIiIiKTEoksU7d+6cWlvnzp21ulbo2cQvv/xS432Evo+IiMgYMSiSRcvJycG1a9fU2lu1aqXV9UKziR07dtR4HwZFIiIyFQyKZNGio6ORm5ur1t64ceNiry1uNlHsPpGRkSUrkoiIyEAYFMmiPXjwQK3N3d0d5cqVK/ba4mYTAQg+o5ieno7k5OSSF0tERCQxBkWyaEIbZFeqVKnY64RmE+fMmaP1vTRtzE1ERGQsGBTJognN7Hl4eBR73cuzid26dYOvr6/aOAcHBzg5OWn1vURERMbGxtAFEBlSVlaWWptcLtd4jTbPJr7IwcEB6enpRdoyMzO1L1KDuLg4xMfHl+ialJQUnD17Fq6urihXrhyqVaum9VZAREQkvezs7CKPSvn7+2v1iJQuMCiSRRMKbMWFplWrVmk1m1jAwcFBrU0ooJbGsmXLMHfuXJ3ci4iITMPu3bvx+uuvS/JdXHomeolMJhPtK+lsotj9eCASERGZAgZFsmhCs33Z2dmi41etWoWHDx8Wfu7atavG2URAePZQ6HuJiIiMDZeeyaIJPY8otixcmtlEsfsV9xyktiZPnoyhQ4eW6Jrr169j2LBhhZ93796NunXrFh3zzx3s/eaQ4PVvBw9FxdrFv/BDRES6ERERgYEDBxZ+rlatmmTfzaBIFs3V1VWt7enTp4JjhWYT/fz8NN4/OzsbaWlpau0uLi4lrFSYp6cnPD09y3SPunXrokmTJkXaalWtjdOLryA/T32JXBlthSavNVFrJyIiaUj5AiKXnsmieXt7q7XFxsaqtZV2NvHJkyeC7VWqVNGuQANxdJOjVpvqgn3X/74jcTVERGQoDIpk0apWrarWlpiYiJSUlCJtpZlNBIB79+6ptcnlcpQvX74U1UqrUbd6gu3RFx4i9Wm6YB8REZkXBkWyaNWrV4e1tbVa+40bNwp/X9rZROD584Avq1mzZolqNJQm3eoLtqtUwLW/bklcDRERGQKDIlk0uVyOhg0bqrVfuHCh8PcvzyZ26dJFq9lEALh48aJaW8uWLUteqAFUrF0BHrWEZz4vh94QbCciIvPCoEgWr3Xr1mptR44cAVC22UQAOHz4sFbfZ4xkMhma9VIP0QBw58Q9ZCTr5nQZIiIyXgyKZPE6deqk1nbo0CHk5uYKziYKjRdy9+5d3Lmj/uKHtrORxqBF38aC7fnKfFz767bE1RARkdQYFMni9e/fX+30lISEBOzZs6dMs4mbNm1Sa/Py8kKbNm1KVachVG3mBfcqboJ9l/dz+ZmIyNwxKJLFq1SpEjp06KDWPnTo0FLPJmZnZyM4OFitfcCAARqPCDQ2MpkMzfs0Euy7ffQuslLFT7EhIiLTx6BIBGDixInFjpkzZ47W91u+fDliYmJK9T3GpplIUFTm5OHaIb79TERkzhgUiQAMHz4cXl5eov2dO3eGv7+/VveKjIzEF198odbeqVMnvPLKK6Ut0WBqtKwK10rCJ8mc33NV4mqIiEhKDIpEAOzs7DB37lzR/lmzZml1n4SEBPTr10/w2L6vv/661PUZkpWVTPSlltvH7nLzbSIiM8agSPT/FAoF2rZtK9j31VdfCS4lv+jChQvw8fHBzZs31frGjBmDjh076qROQ2j1elPB9vw8Ff437Uc8fvxY4oqIiEgKDIpE/8/Kygq//vor3N3d1foOHz6MunXrYvLkydi/fz+ioqKQlJSEu3fvYtu2bRg2bBjatGmDiIgItWsbNGiAn3/+WYofQW886rsj31Ep2Hf3SDSqV6+OiRMnIisrS+LKiIhIn2wMXQCRMalTpw62bduG/v37q4WezMxMLF++HMuXL9f6fp6entizZw/c3IS3mDEFmZmZ6Nu3L3JiZOhQrotav7d9NTjBBcHBwbh58yZCQ0Mhl8sNUCkREekaZxSJXtKtWzccOnRIcGaxJGrXro3jx4+jQYMGOqrMMAIDAxEWFoYb6VdExzRyag4ACAsLQ2BgoFSlERGRnjEoEgno2LEjLl++jAEDBpT4WplMhoCAAJw/fx5169bVQ3XSiYmJQUhICADgmTIBsdmPBMc1dnoFwPP9IUNCQvjMIhGRmWBQJBJRtWpV7NmzB8eOHcOIESPg5OSkcXz58uUxefJkXLlyBUFBQSa93FwgJCQESuV/zybeSL8sOK6cbXlUs68BAFAqlYXhkoiITBufUSQqhq+vL3x9faFUKnHx4kVcv34diYmJSE9Ph6urKzw8PNCiRQs0atTIpE5d0UZ4eHiRzzfTL6OTe09Yy6zVxjZxboUH2fcBAMeOHZOiPCIi0jMGRSIt2djYoE2bNiZ1VnNZJScnF/mcmZ+ByIxbqOekvq9ifcfGOJy4D9mqLKSkpEhVIhER6RGXnolIlNDy+dW084Jjbaxs0cDp+X6Lrq6ueq2LiIikwaBIRKJ8fX3V2u5n3UWaUnjGsKlzSwCAn5+fXusiIiJpMCgSkSiFQgEbm6JPqKiQj2vpFwXHV7avCi95VSgUCgmqIyIifWNQJCJR3t7egqHvWtoF0WuGtBoJLy8vfZZFREQSYVAkIo0WLVoEf3//Im3PlIl4kHVPcLzjUzdkJGdKURoREekZgyIRaSSXyxEaGoqAgIAiy9AXU08Ljs/NUuLs9ktSlUdERHrEoEhExZLL5QgKCkJUVBTmz5+Pnj17wqOZG3KtcwTHH994Fvn5KomrJCIiXeM+ikSkNW9vb8yaNavw84FFR3Dwp6Nq457eS8Sd45Fo4FdHyvKIiEjHOKNIRKXmM6IVrKyFT6M5vu6MxNUQEZGuMSgSUam5VXZF014NBfuu/30b8fcSJK6IiIh0iUGRiMqk42jhIw1VKuBoyCmJqyEiIl1iUCSiMqnToSa8GlYS7Duz/SLSkzIkroiIiHSFQZGIykQmk8F/vI9gX26WEic3nZO4IiIi0hUGRSIqs5YDmsLV01mwL3zdGSizlRJXREREusCgSERlZmNnjY5j2wr2pcan4fyeqxJXREREusCgSEQ68eqoNrCT2wr2HQ4+wQ24iYhMEIMiEemEYzk52g59RbAv7u5TXNl/Q9qCiIiozBgUiUhnOinaQ2YlvAH3oZ+PQaXirCIRkSlhUCQinfGoUR6v9G8i2Bdz4wluHL4jcUVERFQWDIpEpFPdpviK9h36JZyzikREJoRBkYh0yquBJ5r0aCDYF3X+Ie6evC9tQUREVGoMikSkc92nis8qHvzpKGcViYhMBIMiEelc9RZVUN+vtmDf3VNRuHPinsQVERFRaTAoEpFedJ/qJ9q3//vDnFUkIjIBDIpEpBd12tdA3Q41BfuiLjziG9BERCaAQZGI9Kb3B51F+0K/P8LTWoiIjByDIhHpTa021dGwc13BvpjrsTythYjIyDEoEpFeaZpVPPDjEeQp86UrhoiISoRBkYj0qlozbzTr3VCw70nEU5zZflHagoiISGsMikSkd70CO0MmfAQ0Dvx4BNnpOZLWQ0RE2mFQJCK982rgiZavNxPsS4lLw5GVJyWuiIiItMGgSESS6PNBF1jbWQv2HVlLine truncated
"text/plain": [ "text/plain": [
"<Figure size 600x1200 with 1 Axes>" "<Figure size 600x1200 with 1 Axes>"
] ]
@@ -1220,13 +1233,13 @@
}, },
{ {
"cell_type": "code", "cell_type": "code",
"execution_count": 28, "execution_count": 27,
"id": "d8b29d01", "id": "d8b29d01",
"metadata": {}, "metadata": {},
"outputs": [ "outputs": [
{ {
"data": { "data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAA40AAATDCAYAAADSodCrAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAuIwAALiMBeKU/dgABAABJREFUeJzs3Xd4VFX+x/HPpPeE0EOXJl2qgEAUkSIiogsoCrpGBVl1RVZdXSws6rqu/kBXVxBjwy4i6CoSQIkUUVCqGIpI6CS0JKRPZn5/uETC3IGbZHJnMnm/nofnIeecOfc7IcZ8cs491+Z0Op0CAAAAAMBAgLcLAAAAAAD4LkIjAAAAAMAtQiMAAAAAwC1CIwAAAADALUIjAAAAAMAtQiMAAAAAwC1CIwAAAADALUIjAAAAAMAtQiMAAAAAwC1CIwAAAADALUIjAAAAAMAtQiMAAAAAwC1CIwAAAADALUIjAAAAAMAtQiMAAAAAwC1CIwAAAADALUIjAAAAAMAtQiMAAAAAwC1CIwAAAADALUIjAAAAAMAtQiMAAAAAwC1CIwAAAADALUIjAAAAAMAtQiMAAAAAwC1CIwAAAADALUIjAAAAAMAtQiMAAAAAwC1CIwAAAADALUIjAAAAAMAtQiMAAAAAwK0gbxcAVBdOp1Pbtm3TTz/9pCNHjig3N1fh4eGqW7eu2rVrp06dOikoiP+kAAAA4F/4CRc4j23btumFF17QJ598ooyMDLfjYmNjNWLECN19993q1auXhRUCAAAAVcfmdDqd3i4C8EXZ2dl68MEH9corr8jhcJTrtWPGjNELL7yg+vXrV1F1AAAAgDUIjYCB3bt366qrrtLPP/9c4TkaN26sTz/9VF27dvVgZQAAAIC1CI3AWfbt26d+/fpp7969lZ4rPj5eK1asUKdOnTxQGQAAAGA9Tk8FzlBUVKRRo0YZBkabzaaxY8fqiy++UEZGhoqLi3Xs2DEtX75ct912m4KDg11ec/z4cY0cOVJZWVlWlA8AAAB4HCuNwBn+9re/6amnnnJpr1u3rj788ENdeumlbl+7ZcsWXXPNNdq9e7dL34QJE/Tmm296slQAAADAEoRG4H92796t9u3bq7CwsEx7ZGSkVq9erS5dupx3jn379qlXr146fPhwmXabzaY1a9aod+/eHq0ZAAAAqGpsTwX+5+mnn3YJjJL07LPPmgqMktSkSRO9/vrrLu1Op1MzZsyodI0AAACA1VhpBPTbvYcJCQkuobF9+/basmWLAgLK9/uVoUOHasmSJS7taWlpatu2baVqBQAAAKzESiMg6aOPPjJcZbzvvvvKHRgl6S9/+Yth+9tvv13uuQAAAABvYqURkHTFFVdo2bJlZdpCQ0OVkZGhmJiYcs/ncDjUpEkTHTx4sEx727ZtlZaWVqlaAQAAACux0ogar6CgQKtWrXJpHzBgQIUCoyQFBARo2LBhLu3bt2/Xvn37KjQnAAAA4A2ERtR469atU0FBgUv7uR6vYcZll11m2J6amlqpeQEAAAArERpR4/3444+G7d27d6/UvD169DBs37BhQ6XmBQAAAKxEaESNt3nzZsP2Dh06VGreVq1aKSQkxKV9y5YtlZoXAAAAsBKhETXe7t27XdrCw8PVqFGjSs0bGBio5s2bm7oeAAAA4KsIjajx0tPTXdoaNmwom81W6bkTEhJc2vbu3SuHw1HpuQEAAAArEBpR42VkZLi0NWjQwCNzG81TXFyskydPemR+AAAAoKoFebsAwJuKioqUm5vr0h4bG+uR+d3Nc/z4ccXHx1dq7oyMDGVmZpbrNdnZ2Vq/fr1iYmIUFxenJk2aKDQ0tFJ1AMCZCgsLyzxaKDExUXFxcd4r6Bz4PgrAF/ni91FCI2q0nJwcw/aoqCiPzB8dHV2u65bHf/7zH02fPr3S8wBAVVq4cKFGjhzp7TIM8X0UQHXgC99H2Z6KGq2oqMiw3ejU04oIDg42bC8sLPTI/AAAAEBVIzSiRisuLjZsDwryzCK8u9Do7roAAACAr2F7Kmq0gADj35t46nRTd/O4u255TJ48WaNHjy7Xa7Zt26YxY8aUftw6/EqFBbq/f7N+QIzpuSOCzb0np5zm5wwyN2d8iPmTbs3+yxaWmJ5SeXZz76mgxNy47GK76Wsfc5wyNc5hM/+GgpzGv+w4W05Aluk5c5yuB04ZCZT5Vf52tvamxtUN8/z/6mJMfs0FleMQ5iKTX5zHCsx/f9pf6HrPdmWF2c799ZFXckIb8haUftykSROP1+Apnvg+unDhQrVq1UqSlJ15SrPHzTN83TWPDVGbfhdUvFgANcauXbt0zTXXlH7sC99HCY2o0dxtQ7Xbzf/Qfi7u5vHE9td69eqpXr16lZojLDBWEYG13fbHBMaZnisqKNDUuPKExiiTQbReqPmfzE3mNhWUIzSespmbNM/kOKfD/Ep0vs3cARwlNvNf08FOc1+fxQHmP+8FTnPhJUhhpueMsdU1NS7ezYp/ZdQyGRpNfglLMv+LCkeJ+dB4wl75+6fPFmEr3/cvXz4kxhPfR1u1aqUOHTqUfvzVBd/p+P6TLuMCMkLKjAMAs3zh+yjbU1GjhYUZ/4Can5/vkfnz8vIM28PDwz0yPwDAt7To1dSw/df1+wzbAaA6IDSiRouJiTG8fzE7O9sj87s7JbVWrVoemR8A4Fta9DDeRnZg6yEV5hkfvgYAvo7QiBrPKMAdO3bMI3MfPXrUsL2yz2gEAPimFj2NVxodJU6lb9hvcTUA4BmERtR4jRo1cmk7fPiwR+Y2mic+Pp7tqQDgp+q3qqOIWsbf439dxxZVANUToRE1XvPmzV3aMjIy3N6PWB579uxxaWvRokWl5wUA+CabzaYW3Y23qP66bq/F1QCAZxAaUeO1bdvWpc3pdGrnzp2VmjcrK0uZmZku7W3atKnUvAAA3+Zui+qv6/epuMAzp3MDgJUIjajxunbtati+cePGSs27YcOGcl0PAOAfWl7czLDdXmjXrz+w2gig+iE0osbr1auXYfvatWsrNe+3335brusBAPxD404NFR5j/Einnat+tbgaAKg8QiNqvBYtWhjeZ7h06dJKzWv0+oiICPXp06dS8wIAfFtAYIBa9W1u2Ldj1W5riwEADyA0ApKGDRvm0vbLL7+43WJ6PhkZGfrmm29c2gcNGqSQkJAKzQkAqD7a9LvAsP3A1kPKPVH5g9YAwEqERkDSuHHjDNtnz55dofnmzp2rkpIS09cBAPgXd6HR6ZR2rdljbTEAUEmERkBS37591b59e5f2N9980/CxGeeSlZWl559/3qW9Xr16GjlyZEVLBABUI7Wb1VKtRrGGfe//e74OHTpkcUUAUHGERkC/PVfr/vvvd2kvLCzUHXfcIafTaXquqVOnGj5q45577lFYmPHBCAAA/1JQUKAM20HDvowtx9W0aVNNmjRJBQUFFlcGAOVHaAT+Z/z48erSpYtL+9KlS3X33XebCo7/+te/lJyc7NLeqFEj3XvvvZ4oEwDg4/Lz8zVs2DAt2fC5YX9ccLwiFa05c+Zo6NChys/Pt7hCACgfQiPwP4GBgZo9e7aCgoJc+l566SUNGzZMO3fuNHztwYMHddNNN+mBBx4w7H/xxRcVGRnp0XoBAL5pypQpSk1N1d4C94/XaBbWSpKUmpqqKVOmWFUaAFSI60/HQA3Wu3dvPf300/rLX/7i0rdkyRK1bdtWffr0Uc+ePVWrVi1lZ2dr48aN+uabb2S32w3nvPfee3XNNddUceUAAF9w8ODB0h0nBY48HSk8qPqhCS7jLghvo82n1kmSkpOT9dhjj6lhw4aW1goAZhEagbNMnTpVGRkZeuaZZ1z6nE6n1qxZozVr1pia66abbtJzzz3n6RIBAD4qOTm5zC8R9xTsMgyNTcNbKMgWLLuzWHa7XcnJyZo2bZqVpQKAaWxPBQz885//1H/+8x+FhoZW6PWBgYF67LHH9NZbbykggP/MAKCmWLVqVZmPd+dtNxwXZAtWs7CWpR+vXLmySusCgMrgp1nAjTvvvFNbt27VddddV67gN2jQIK1bt06PP/64bDZbFVYIAPA1WVlZZT4+VHRAuSWnDMe2jGhb+vfs7OwqrQsAKoPtqcA5tGrVSvPnz9eePXv0ySefaMWKFdq2bZsOHz6s/Px8hYWFqU6dOmrXrp369++vUaNGqV27dt4uGwDgJbGxZz+b0alf83eoY1Q3l7EtwttIsklyKiYmxoryAKBCCI2ACc2bN9eUKVM44Q4AcE79+vVTSkpKmbZf8rYbhsbIwCg1DGmkQ0X71b9/f6tKBIByY3sqAACAhyQlJbk8uim94BfZncWG4y+IaKugoCAlJSVZUR4AVAihEQAAwEMSEhJcAqDdWay9+cbPbGwZ3lZJSUk8bgOATyM0AgAAeNDMmTOVmJhYpu2XfONTVOuE1NejUx+3oCoAqDhCIwAAgAeFh4dr8eLFmjhxYulW1V/dhEZJ2plqvAoJAL6C0AgAAOBh4eHhmj17ttLT0zVjxgz1vbyPTgVnGY7d9Pk2i6sDgPIhNAIAAFSRhIQETZs2TUuWLNHoP19jOCZ9wwGdOGAcKAHAFxAaAQAALNDlyvZu+zYvZrURgO8iNAIAAFig7gW1ldCuvmEfW1QB+DJCIwAAgEW6DDdebUzfcEAnDrJFFYBvIjQCAABY5JxbVL9gtRGAbyI0AgAAWKTuBbXV8ELjLao/LtxqcTUAYA6hEQAAwEIXXWW82rh/6yEd2ZVpcTUAcH6ERgAAAAt1vbqj274fF26xsBIAMIfQCAAAYKHaTWupeY8mhn0/Ltwqh8NpcUUAcG6ERgAAAIt1H9nJsP34/pPa88M+i6sBgHMjNAIAAFisy1XtFRBk/GPYj5+wRRWAbyE0AgAAWCyyVoQuvLSVYd+Gz7aquKDY4ooAwD1CIwAAgBd0v8Z4i2pBTqG2LEmzuBoAcI/QCAAA4AUdBrVVeEyYYd/3H26Line truncated "image/png": "iVBORw0KGgoAAAANSUhEUgAAA4wAAATCCAYAAAD2P2gwAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAuIwAALiMBeKU/dgABAABJREFUeJzs3Xd4VFX+x/HPTOoQSAid0JsUBaSDEEBAmqCAYsW1RBfsi2V3RRQRXV0b2BDUgNgFQVEExAKhSJMmUqSHEkInAVInM78/+JElzJ3kJpncSXm/nodH5pxzz3wTY8wn59xzbW632y0AAAAAAC5h93cBAAAAAIDiicAIAAAAADBEYAQAAAAAGCIwAgAAAAAMERgBAAAAAIYIjAAAAAAAQwRGAAAAAIAhAiMAAAAAwBCBEQAAAABgiMAIAAAAADBEYAQAAAAAGCIwAgAAAAAMERgBAAAAAIYIjAAAAAAAQwRGAAAAAIAhAiMAAAAAwBCBEQAAAABgiMAIAAAAADBEYAQAAAAAGCIwAgAAAAAMERgBAAAAAIYIjAAAAAAAQwRGAAAAAIAhAiMAAAAAwBCBEQAAAABgiMAIAAAAADBEYAQAAAAAGAr0dwFASZGVlaXNmzfrzz//1IkTJ3T27FlFRESoatWqatu2rZo0aeLvEgEAAACfIjACeVi7dq3effddzZkzR2fOnPE6rlq1arrjjjt0//33q1GjRhZWCAAAABQNm9vtdvu7CKA4OnbsmB599FF98cUX+bouMDBQjz/+uMaNGyeHw1FE1QEAAABFj8AIGNi4caMGDBigxMTEAs/RqlUrLVy4UDVr1vRhZQAAAIB1CIzAJVavXq2+ffsqOTm50HPVr19fS5cuVZ06dXxQGQAAAGAtTkkFLpKQkKChQ4cahsWQkBCNHDlSP/zwg/bu3auTJ09q586d+vLLL3XDDTfIZrN5XLNv3z4NHTpUqampVpQPAAAA+BQrjMBFevfurV9//dWjvWvXrvriiy9yXSn8/fffdf311yshIcGj76GHHtLbb7/t01oBAACAosYKI/D/vvrqK8Ow2K9fPy1evDjPbaXt27fXvffea9j33nvv6Y8//vBJnQAAAIBVWGEEJLlcLjVu3Fh79+7N0d6wYUNt2rRJ5cuXz3OO9PR0NWrUSIcOHTLsHzhwoH744Qef1AsAAABYgRVGQNLcuXM9wqIkTZ061VRYlKQPPvjAa1iUpAULFmjHjh0FrhEAAACwGoERkDR58mSPtu7du6tPnz6mrk9PT9fLL7+c6xi3260pU6YUqD4AAADAHwiMKPNOnTqlJUuWeLQ/+uijpue4dHXxmmuu0bPPPusx7ptvvilQjQAAAIA/cA8jyrxPP/1Ud9xxR462iIgIHTlyRCEhIXleb3Tv4ooVK1SlShU1bdrUY/wff/yhli1bFr5wAAAAoIixwogyz2h1sU+fPqbCouS5utinTx9dddVVuuyyy9SwYUOP8XFxcQWuFQAAALASgRFl3rp16zzaevbsaepao3sXn3vuuVznMXo/AAAAoDgiMKJMy8jI0JYtWzza27Zta+p6o9XFrl275joPgREAAAAlBYERZdr+/fuVmZnp0d6iRYs8r81rddHbPHv27MlfkQAAAICfEBhRph04cMCjLTIyUhUrVszz2rxWFyUZ3sN47tw5JSUl5b9YAAAAwGIERpRpFwe+C6pXr57ndUari+PGjTM9l9H7AgAAAMUNgRFlmtFKX5UqVfK87tLVxd69e6tbt24e40JDQxUWFmbqfQEAAIDiJtDfBQD+lJaW5tHmcDhyvcbMvYsXCw0N1blz53K0paammi8yF0ePHtWxY8fydU1ycrJ+//13hYeHq2LFiqpTp47pR4gAgBnp6ek5tvz36NHD1FZ/f+D7KIDiqDh9HyUwokwzCm55/U//ww8/NLW6eEFoaKhHm1FQLYjJkydr/PjxPpkLAIrKt99+q+uvv97fZRji+yiAksCf30fZkgpcwmazee3L7+qit/ncbneBagMAAACsRGBEmWa0+peenu51/IcffqiDBw9mv+7Vq1euq4uS8Wqi0fsCAAAAxQ1bUlGmGd2v6G27aEFWF73Nl9d9kmY98MADGj58eL6u2bp1q2666abs100cAxUaEOF1fHV7uOm5ywWZ+x2UW+ZXWMsFmpuzUrD3leFLuUyOS88yPaVSnOY+prQsc+OSM52m3/uE66ypcS6b+Q8o0B1katwZu/kDnM64j5oaF6Bg03M2t+X9zFRJqhrq+//dhZv8mgs0/6WpDJNfnCfSzH4VSwfTz+U9KJ9Cbbl/faRkndKGlDnZr+vUqePzGnzFF99Hv/32WzVu3FiS9OPEJdo0f5vHNVUaVNI9799cuGIBlBm7du3SkCFDsl/78/sogRFlWni4Zxg6fvy44Vij1cXo6Ohc509PT9fZs54/zFeoUCGflRqrVq2aqlWrVqg5QgMiVC6gstf+8ICKpucqHxhgalx+AmN5kyG0Woj5n8pNZjal5SMwnrWZmzTF5Di3K9P0e6fazB22kWUzH0KD3OZCW6bd/Oc9zW0uuATK/Ap8uK2qqXGVgswF4PyINBkYTX4JSzL/SwpXlvnAeMp5xnwBJpWzmQ/1Ut73hvuTL76PNm7cWJdffrkkKXVglg79fMJz0CGpQa2GKlfRN78wBFC2+PP7KFtSUaZFRUV5tCUmJnq0FXR18ciRI4bttWrVMlcgAKBEadihrte+fesOeO0DgOKKwIgyrXbt2h5tJ0+eVHJyco62gqwuStLevXs92hwOhypVqlSAagEAxV1k7QhF1DDeRbJn7X6LqwGAwiMwokyrW7euAgI8t1Fu2/a/+08Kuroonb/P5VL169fPV40AgJLDZrOpQXvjVca9BEYAJRCBEWWaw+FQs2bNPNo3bNiQ/fdLVxevvvpqU6uLkrRx40aPtjZt2uS/UABAidGgo3FgPPBHgjLTzN+fDADFAYERZV67du082pYsWSKpcKuLkrR48WJT7wcAKD0atjc+zTAr06UDfyRYXA0AFA6BEWVe9+7dPdp+/vlnZWZmGq4uGo03snv3bu3cudOj3ezqJACgZKrRtJpCKxifaLh7VbzF1QBA4RAYUeYNGjRINlvO4/FPnDihuXPnFmp18bPPPvNoq1mzptq3b1+gOgEAJYM9wK767YxXGXcs32NxNQBQOARGlHnVq1dXly5dPNqHDx9e4NXF9PR0TZ061aP9uuuu8winAIDSp/FV9Q3b4zccVPq5DGuLAYBCIDACkkaNGpXnmHHjxpme77333lNCgud9KmbeBwBQ8jXt1tCwPSvTpT1r2JYKoOQgMAKSbr75ZtWsWdNrf8+ePdWjRw9Tc+3Zs0fPPPOMR3v37t115ZVXFrREAEAJUqNZdZWvXM6wj22pAEoSAiMgKTg4WOPHj/faP3bsWFPznDhxQtdee63Onj3r0fef//ynwPUBAEoWu92mJl2NVxl3rthrcTUAUHAERuD/xcTEqEOHDoZ9L774ouEW04tt2LBBnTt31vbt2z367rjjDnXt2tUndQIASoYm3RoYth/eflTPj5mgw4cPW1wRAOQfgRH4f3a7XV988YUiIyM9+hYvXqzGjRvrgQce0MKFCxUfH69Tp05p9+7dmjVrlm666Sa1b99eu3bt8ri2adOmevvtt634EAAAxUjddlFe+754a5bq1q2rUaNGKS0tzcKqACB/Av1dAFCcNGrUSLNmzdKgQYM8/geempqq9957T++9957p+apVq6a5c+cqIiLC16UCAIqx1NRU3RpzsxpktlKloCoe/fVCG2n7uc2aOnWqtm/frgULFsjhcPihUgDIHSuMwCV69+6tn3/+2XClMT8aNmyoFStWqGnTpj6qDABQUowePVpxcXHan7rbsL9eaKPsv8fFxWn06NFWlQYA+UJgBAx07dpVf/zxh6677rp8X2uz2TRy5EitX79ejRs3LoLqAADFWUJCgmJjYyVJ8WnGJ6KWDwxXlaDq2a9jY2O5pxFAsURgBLyoXbu25s6dq2XLlunWW29VWFhYruMrVaqkBx54QJs3b9aUKVPYhgoAZVRsbKycTqck6WDaPrncWYbjGpX73w4Up9OZHTIBoDjhHkYgD926dVO3bt3kdDq1ceNGbd26VSdPntS5c+cUHh6uKlWqqHXr1mrevLlsNpu/ywUA+Nny5cuz/57uTtOh9P2qE+p5YmojR1OtTlqa/XrZsmWW1AcA+UFgBEwKDAxU+/bt1b59e3+XAgAoxpKSknK83p3yl2FgrBFSW2EB5XUu6/yze5OTky2pDwDygy2pAAAAPnTpLQl7Uv/yOraB47Lsv4eHhxdZTQBQUARGAAAAH+rWrVuO16edJ3Ui46jh2EaO/93HGB0dXaR1AUBBEBgBAAB8KCYmRoGBOe/62ZO6w3Bs3dCGCrQFKTAwUDExMVaUBwD5QmAEAADwoaioKI/wtzt1u+HYIHuw6oY2UExMjGrWrGlFeQCQLwRGAAAAH5s4caJ69OiR/fpw+kGlZJ0zHNu1YU9NmjTJosoAIH8IjAAAAD7mcDi0YMECjRw5UoGBgXLLrb2pOw3HNnBcpuDgEIsrBABzCIwAAABFwOFwaMqUKYqPj9eECRNU4bJQw3Fnj53Tvt/3W1wdAJhDYAQAAChCUVFRGjt2rD5ZNF2BwQGGYzbN32ZxVQBgDoERAADAAqHlQ9S0R2PDvj8WbJPL5ba4IgDIG4ERAADAIq2vbWHYnnzkjPatO2BxNQCQNwIjAACARS7vfZnXbal//LDV4moAIG8ERgAAAIuEVghR0+6NDPs2sS0VQDFEYAQAALBQbttS966Jt7gaAMgdgREAAMBCl/dpqgAv21LXfbvZ4moAIHcERgAAAAuFVghRi6ubGPZt+mGrMtOcFlcEAN4RGAEAACzWdmhLw/a0M+natmSnxdUAgHcERgAAAIs179lEoRVCDPvWf8O2VADFB4ERAADAYkGLine truncated
"text/plain": [ "text/plain": [
"<Figure size 780x1200 with 2 Axes>" "<Figure size 780x1200 with 2 Axes>"
] ]
@@ -1316,7 +1329,7 @@
}, },
{ {
"cell_type": "code", "cell_type": "code",
"execution_count": 29, "execution_count": 28,
"id": "23a55aa1", "id": "23a55aa1",
"metadata": {}, "metadata": {},
"outputs": [ "outputs": [
@@ -1333,7 +1346,7 @@
}, },
{ {
"data": { "data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAA40AAATDCAYAAADSodCrAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAuIwAALiMBeKU/dgABAABJREFUeJzs3Xd4VFX+x/HPpPdGCYQuTbpU6SCigNhYRRQF3Y2KsjZ01VURRdR1XV3RtYCIDVfFBURdRSlKBBTpTQhFIHQSWnqbzPz+8EeWMHfgJpncmUzer+fhecg5Z879JoaYz5xzz7U5nU6nAAAAAAAwEODtAgAAAAAAvovQCAAAAABwi9AIAAAAAHCL0AgAAAAAcIvQCAAAAABwi9AIAAAAAHCL0AgAAAAAcIvQCAAAAABwi9AIAAAAAHCL0AgAAAAAcIvQCAAAAABwi9AIAAAAAHCL0AgAAAAAcIvQCAAAAABwi9AIAAAAAHCL0AgAAAAAcIvQCAAAAABwi9AIAAAAAHCL0AgAAAAAcIvQCAAAAABwi9AIAAAAAHCL0AgAAAAAcIvQCAAAAABwi9AIAAAAAHCL0AgAAAAAcIvQCAAAAABwi9AIAAAAAHCL0AgAAAAAcIvQCAAAAABwi9AIAAAAAHAryNsFANWF0+nU1q1b9euvv+ro0aPKzc1VeHi46tSpozZt2qhDhw4KCuKfFAAAAPwLv+EC57F161a99tpr+vzzz5Wenu52XGxsrK666irde++96tGjh4UVAgAAAFXH5nQ6nd4uAvBFWVlZevTRR/X222/L4XCU67U33HCDXnvtNSUmJlZRdQAAAIA1CI2Agd27d+vKK6/Utm3bKjxHw4YN9eWXX6pz584erAwAAACwFqEROMv+/fvVt29f7du3r9JzJSQkaOnSperQoYMHKgMAAACsx+mpwBmKioo0YsQIw8Bos9k0atQoffPNN0pPT1dxcbGOHz+uJUuW6Pbbb1dwcLDLa06cOKFrrrlGmZmZVpQPAAAAeBwrjcAZnnjiCT3//PMu7XXq1NFnn32mgQMHun3t5s2bde2112r37t0ufWPHjtUHH3zgyVIBAAAASxAagf+3e/dutW3bVoWFhWXaIyMjtWLFCnXq1Om8c+zfv189evTQkSNHyrTbbDb99NNP6tmzp0drBgAAAKoa21OB//fCCy+4BEZJeumll0wFRklq1KiR3nvvPZd2p9OpKVOmVLpGAAAAwGqsNAL6/d7DpKQkl9DYtm1bbd68WQEB5Xt/ZejQofruu+9c2lNTU9W6detK1QoAAABYiZVGQNJ//vMfw1XGBx98sNyBUZL+8pe/GLZ/9NFH5Z4LAAAA8CZWGgFJl112mRYvXlymLTQ0VOnp6YqJiSn3fA6HQ40aNdKhQ4fKtLdu3VqpqamVqhUAAACwEiuNqPEKCgq0fPlyl/b+/ftXKDBKUkBAgIYNG+bSvn37du3fv79CcwIAAADeQGhEjbd69WoVFBS4tJ/r8RpmXHLJJYbtKSkplZoXAAAAsBKhETXeunXrDNu7du1aqXm7detm2L5+/fpKzQsAAABYidCIGm/Tpk2G7e3atavUvC1atFBISIhL++bNmys1LwAAAGAlQiNqvN27d7u0hYeHq0GDBpWaNzAwUE2bNjV1PQAAAMBXERpR46Wlpbm01a9fXzabrdJzJyUlubTt27dPDoej0nMDAAAAViA0osZLT093aatXr55H5jaap7i4WKdOnfLI/AAAAEBVC/J2AYA3FRUVKTc316U9NjbWI/O7m+fEiRNKSEio1Nzp6enKyMgo12uysrK0Zs0axcTEKC4uTo0aNVJoaGil6gCAMxUWFpZ5tNCAAQMUFxfnvYLOgZ+jAHyRL/4cJTSiRsvOzjZsj4qK8sj80dHR5bpuebz55puaPHlypecBgKo0f/58XXPNNd4uwxA/RwFUB77wc5TtqajRioqKDNuNTj2tiODgYMP2wsJCj8wPAAAAVDVCI2q04uJiw/agIM8swrsLje6uCwAAAPgatqeiRgsIMH7fxFOnm7qbx911y2P8+PEaOXJkuV6zdetW3XDDDaUftwy/QmGB7u/fTAyIMT13RLC5z8kpp/k5g8zNmRBi/qRbs/9lC0tMT6k8u7nPqaDE3LisYrvpax935Jga57CZ/4SCnMZvdpwtOyDT9JzZTtcDp4wEyvwqfxtbW1Pj6oR5/n91MSa/54LKcQhzkclvzuMF5n8+HSh0vWe7ssJs5/7+yCs5qfV580o/btSokcdr8BRP/Bx944031KRJE0lS3okC/ffRHwxf1+fPXdTgosSKFwugxkhLS9Of//zn0o994ecooRE1mrttqHa7+V/az8XdPJ7Y/lq3bl3VrVu3UnOEBcYqIrCW2/6YwDjTc0UFBZoaV57QGGUyiNYNNf+bucncpoJyhMYcm7lJ80yOczrMr0Tn28wdwFFiM/89Hew09/1ZHGD+617gNBdeghRmes4YWx1T4xLcrPhXRrzJ0GjyW1iS+TcqHCXmQ+NJe+Xvnz5bhK18P798+ZAYT/wcbdKkiVq2bFn68ZZGe5RzNM9lXOjJqDLjAMAsX/g5yvZU1GhhYca/oObn53tk/rw8118cJCk8PNwj8wMAfEu9dsZvZhz5tXyntAKALyE0okaLiYkxvH8xKyvLI/O7OyU1Pj7eI/MDAHxL/fbGofHYzpOyF3pmFwsAWI3QiBrPKMAdP37cI3MfO3bMsL2yz2gEAPimxHa1DdsddofSt5+wuBoA8AxCI2q8Bg0auLQdOXLEI3MbzZOQkMD2VADwUwlNYhUSaXwf7ZHNbFEFUD0RGlHjNW3a1KUtPT3d7f2I5bF3716XtmbNmlV6XgCAb7IF2FTPzWrjwQ1HLa4GADyD0Igar3Xr1i5tTqdTO3furNS8mZmZyshwfVe5VatWlZoXAODb6nc0PpH16NZj3NcIoFoiNKLG69y5s2H7hg0bKjXv+vXry3U9AIB/cPc8xpJih45sNb7XHQB8GaERNV6PHj0M21euXFmpeX/++edyXQ8A4B9qNY9TaLTx8ywPrmeLKoDqh9CIGq9Zs2aG9xkuWrSoUvMavT4iIkK9evWq1LwAAN8WEBigpE7GW1QPERoBVEOERkDSsGHDXNp+++03t1tMzyc9PV0//vijS/vgwYMVEmL87jMAwH8kudmimrHzpApziyyuBgAqh9AISBo9erRh+7Rp0yo034wZM1RSUmL6OgAA/+Luvkanw6nDm3j0BoDqhdAISOrdu7fatm3r0v7BBx8YPjbjXDIzM/Xqq6+6tNetW1fXXHNNRUsEAFQjcY2iFVHL+Jm8Cz5YpPT0dIsrAoCKIzQCkmw2mx5++GGX9sLCQt15551yOp2m53rooYcMH7Vx3333KSwsrFJ1AgCqh8LCQh2zHTbsy9yeo0suuUSTJk1SYWGhxZUBQPkRGoH/N2bMGHXq1MmlfdGiRbr33ntNBcd//OMfmjlzpkt7gwYN9MADD3iiTACAjysoKNAdd9yhlNTvDfvjg2srSjGaPXu2br/9dhUUFFhcIQCUD6ER+H+BgYGaNm2agoKCXPreeOMNDRs2TDt37jR87aFDh3TLLbfokUceMex//fXXFRkZ6dF6AQC+6fnnn9eqVau0N3+X2zHNwltKklatWqXnn3/eqtIAoEJcfzsGarCePXvqhRde0F/+8heXvu+++06tW7dWr1691L17d8XHxysrK0sbNmzQjz/+KLvdbjjnAw88oGuvvbaKKwcA+IKjR49q7ty5kqR8R66OFB5QvdCGLuOahbfU+uxfJElz587VPffco7p1jR/TAQDeRmgEzvLQQw8pPT1dL774okuf0+nUTz/9pJ9++snUXLfccotefvllT5cIAPBRc+bMKfMm4p78nYahsWFYUwXZgmV3Fstut2vOnDkaP368laUCgGlsTwUM/P3vf9ebb76p0NDQCr0+MDBQTz31lD788EMFBPDPDABqinXr1pX5eE++8W0NQbZgNQprVvrx2rVrq7QuAKgMfpsF3Lj77ru1ZcsWXXfddeUKfoMHD9bq1av19NNPy2azVWGFAABfk52dXebjI0WHlFeSYzj29H2NkpSTYzwGAHwB21OBc2jRooXmzJmjvXv36vPPP9fSpUu1detWHTlyRPn5+QoLC1Pt2rXVpk0b9evXTyNGjFCbNm28XTYAwEuio6PPanFqT/4utYu6yGXsBeGt9L2+liRFRUVVfXEAUEGERsCEpk2basKECZowYYK3SwEA+LAuXbpo+fLlZdr25O8wDI0xQXGqE1xPGcVH1LVrV4sqBIDyY3sqAACAh1x//fUuj25KK/hNJc4Sw/EtI9oqKChI119/vRXlAUCFEBoBAAA8JDExUdddd12ZtkJHgQ4U7DUc3yLiQl133XU8bgOATyM0AgAAeNDjjz+uHj16lGnblbfNcGztkETdc9t9VpQFABVGaAQAAPCgsLAwzZgxQ6NGjSrdqrorP9Xt+EOrM6wqDQAqhNAIAADgYWFhYXrmmWf0/fff6/7771fnXp2UGXjCcOyeFQcsrg4AyofTUwEAAKpIYmKixo8fL0na8Nk2/fLORpcx6anHlXssT5G1I6wuDwBMYaURAADAAs36NHTbt2c5q40AfBehEQAAwAKxDaKV0DTWsG/X0jSLqwEA8wiNAAAAFmnWr5Fh+9Gtx5V9JMfiagDAHEIjAACARVoMbOy2b9fSfRZWAgDmERoBAAAsEtcoRrVbxBv2ERoB+CpCIwAAgIVaXNLEsP3E7lM6sTfT4moA4PwIjQAAABZqPsD4vkZJ2vUDB+IA8D2ERgAAAAtF1Y1UvfZ1DPt+W5omp9NpcUUAcG6ERgAAAIu1uMT4QJysw7k6uvW4xdUAwLkRGgEAACx2Qb9GsgXYDPu2L9xtcTUAcG6ERgAAAIuFx4WpUff6hn2/pexTcb7d4ooAwD1CIwAAgBe0vryZYXtxnl17Vuy3uBoAcI/Line truncated "image/png": "iVBORw0KGgoAAAANSUhEUgAAA4wAAATCCAYAAAD2P2gwAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAuIwAALiMBeKU/dgABAABJREFUeJzs3Xd8VFX+//H3TDpplARI6L0oIL2GLkURBUEsuLob/YJiWdQtIoJgXddVbAhixO4iRbGBWJCm0kF67xACBFJIncz8/vBHljB3kpt2J+X1fDx4LHPOuWc+ycaQd86559pcLpdLAAAAAABcwe7tAgAAAAAAZROBEQAAAABgiMAIAAAAADBEYAQAAAAAGCIwAgAAAAAMERgBAAAAAIYIjAAAAAAAQwRGAAAAAIAhAiMAAAAAwBCBEQAAAABgiMAIAAAAADBEYAQAAAAAGCIwAgAAAAAMERgBAAAAAIYIjAAAAAAAQwRGAAAAAIAhAiMAAAAAwBCBEQAAAABgiMAIAAAAADBEYAQAAAAAGCIwAgAAAAAMERgBAAAAAIYIjAAAAAAAQwRGAAAAAIAhAiMAAAAAwBCBEQAAAABgiMAIAAAAADDk6+0CgPIiJydH27Zt0/bt23Xu3DmlpqYqPDxckZGR6tChg5o1a+btEgEAAIASRWAECrB+/Xq9+eabWrRokVJSUjyOq1mzpu68807dd999atKkiYUVAgAAAKXD5nK5XN4uAiiLzpw5o4cffliffvppoa7z9fXVo48+qqlTpyooKKiUqgMAAABKH4ERMLBlyxYNHTpU8fHxRZ6jbdu2Wrp0qaKiokqwMgAAAMA6BEbgCmvXrtWgQYOUnJxc7LkaNmyolStXql69eiVQGQAAAGAtTkkFLnPy5EmNGDHCMCwGBARo3Lhx+uabb3To0CElJiZq3759+u9//6ubb75ZNpvN7ZrDhw9rxIgRSk9Pt6J8AAAAoESxwghcZsCAAfrpp5/c2nv27KlPP/0035XCDRs26MYbb9TJkyfd+h544AG9/vrrJVorAAAAUNpYYQT+v3nz5hmGxcGDB2v58uUFbivt1KmT7rnnHsO+t956S7///nuJ1AkAAABYhRVGQJLT6VTTpk116NChPO2NGzfW1q1bFRISUuAcmZmZatKkiU6cOGHYf9111+mbb74pkXoBAAAAK7DCCEhavHixW1iUpNmzZ5sKi5I0Z84cj2FRkpYsWaK9e/cWuUYAAADAagRGQNLMmTPd2nr37q2BAweauj4zM1MvvPBCvmNcLpdmzZpVpPoAAAAAbyAwotI7f/68fv75Z7f2hx9+2PQcV64uXnvttZoyZYrbuM8//7xINQIAAADewD2MqPQ++ugj3XnnnXnawsPDdfr0aQUEBBR4vdG9i2vWrFFERIRatGjhNv73339XmzZtil84AAAAUMpYYUSlZ7S6OHDgQFNhUXJfXRw4cKB69Oih5s2bq3Hjxm7jV6xYUeRaAQAAACsRGFHpbdy40a2tb9++pq41unfxqaeeynceo/cDAAAAyiICIyq1rKws7dixw629Q4cOpq43Wl3s2bNnvvMQGAEAAFBeEBhRqR09elTZ2dlu7a1bty7w2oJWFz3Nc/DgwcIVCQAAAHgJgRGV2rFjx9zaqlWrpqpVqxZ4bUGri5IM72G8ePGikpKSCl8sAAAAYDECIyq1ywPfJbVq1SrwOqPVxalTp5qey+h9AQAAgLKGwIhKzWilLyIiosDrrlxdHDBggHr16uU2LjAwUMHBwabeFwAAAChrfL1dAOBNGRkZbm1BQUH5XmPm3sXLBQYG6uLFi3na0tPTzReZj4SEBJ05c6ZQ1yQnJ2vDhg0KCwtT1apVVa9ePdOPEAEAMzIzM/Ns+e/Tp4+prf7ewPdRAGVRWfo+SmBEpWYU3Ar6R/+dd94xtbp4SWBgoFubUVAtipkzZ2ratGklMhcAlJYvvvhCN954o7fLMMT3UQDlgTe/j7IlFbiCzWbz2FfY1UVP87lcriLVBgAAAFiJwIhKzWj1LzMz0+P4d955R8ePH8993b9//3xXFyXj1USj9wUAAADKGrakolIzul/R03bRoqwuepqvoPskzbr//vs1evToQl2zc+dO3XLLLbmvmwVdp0CfcI/ja9nDTM9dxc/c76BcMr/CWsXX3JzV/T2vDF/JaXJcZo7pKZXmMPcxZeSYG5ec7TD93uecqabGOW3mPyBfl5+pcSl28wc4pbgSTI3zkb/pOVvZCn5mqiRFBpb8P3dhJr/mfM1/aSrL5BfnuQyzX8XS8cyLBQ8qpEBb/l8faTnntTltUe7revXqlXgNJaUkvo+++eabatCggSRpwwfbdHDVcbdrwuuFavCU/H/BCACXHDlyRBMmTMh97c3vowRGVGphYe5h6OzZs4ZjjVYXY2Ji8p0/MzNTqanuP8yHhoYWslJjNWvWVM2aNYs1R6BPuKr41PDYH+ZT1fRcIb4+psYVJjCGmAyhNQPM/1RuMrMpoxCBMdVmbtI0k+NczmzT751uM3fYRo7NfAj1c5kLbdl285/3DJe54OIr8yvwYbZIU+Oq+5kLwIVRzWRgNPklLMn8LymcOeYD43lHivkCTKpiMx/qpYLvDfemkvg+2qBBAzVr1kyS5Ortq+S1We6DEqT6UQ0UEFK4zx0ASN79PsqWVFRq0dHRbm3x8fFubUVdXTx9+rRhe506dcwVCAAoV2pf7eGXGC7p9E7jX0gCQFlGYESlVrduXbe2xMREJScn52kryuqiJB06dMitLSgoSNWrVy9CtQCAsi60drCqVDdeJY/fQWAEUP4QGFGp1a9fXz4+7tsod+3alfv3oq4uSn/c53Klhg0bFqpGAED5YbPZPK4yxm8v3PMeAaAsIDCiUgsKClLLli3d2jdv3pz79ytXF/v162dqdVGStmzZ4tbWvn37whcKACg3al8VYdiesCdROdmFuDkaAMoAAiMqvY4dO7q1/fzzz5KKt7ooScuXLzf1fgCAiqP2VcYrjDlZOUrYk2hxNQBQPARGVHq9e/d2a/vhhx+UnZ1tuLpoNN7IgQMHtG/fPrd2s6uTAIDyqUaTqvKrYnwQ/YnNxoehAUBZRWBEpTds2DDZbHmPxz937pwWL15crNXFjz/+2K0tKipKnTp1KlKdAIDywe5jV1Qb40d1EBgBlDcERlR6tWrVUvfu3d3aR48eXeTVxczMTM2ePdutffjw4W7hFABQ8dS5ppZhe8Luc8pON/9MVADwNgIjIGn8+PEFjpk6darp+d566y2dPHmySO8DACj/6lxjvMLodDgVv4PTUgGUHwRGQNKYMWMUFRXlsb9v377q06ePqbkOHjyoJ5980q29d+/euuaaa4paIgCgHKneqKoCwwMM+9iWCqA8ITACkvz9/TVt2jSP/ZMnTzY1z7lz53T99dcrNTXVre+5554rcn0AgPLFZrcp2sMqI4ERQHlCYAT+v9jYWHXu3Nmw79lnnzXcYnq5zZs3q1u3btq9e7db35133qmePXuWSJ0AgPKhTjvj+xjP7E/UzBlvKSEhweKKAKDwCIzA/2e32/Xpp5+qWrVqbn3Lly9X06ZNdf/992vp0qU6cuSIzp8/rwMHDmj+/Pm65ZZb1KlTJ+3fv9/t2hYtWuj111+34kMAAJQhEa2rGrbbZNOS979Xv379NGXKFGVmZlpbGAAUgvFDgoBKqkmTJpo/f76GDRumjIyMPH3p6el666239NZbb5mer2bNmlq8eLHCw8NLulQAQBmWkZGhR6dMVFtHT4X5VnXrbxjUTHvTdmjevHk6dOiQ5syZo8DAQOsLBYACsMIIXGHAgAH64YcfDFcaC6Nx48Zas2aNWrRoUUKVAQDKi+eee07r1q/T4XT3nSeS1CioqaQ/HrO0bt067nMHUGYRGAEDPXv21O+//67hw4cX+lqbzaZx48Zp06ZNatq0aSlUBwAoy06fPq2FCxdKkg6l7zUcE+wTqpr+tXNfL1y4kHsaAZRJBEbAg7p162rx4sVatWqVbrvtNgUHB+c7vnr16rr//vu1bds2zZo1i22oAFBJLViwQA6HQ5J0NOOQclwOw3GNgprn/t3hcGjBggWW1AcAhcE9jEABevXqpV69esnhcGjLli3auXOnEhMTdfHiRYWFhSkiIkLt2rVTq1atZLPZvF0uAMDLNm3alPv3bFeWjmccUYOgJm7jGgU109qkFbmvN27caEl9AFAYBEbAJF9fX3Xq1EmdOnXydikAgDIsJSUlz+uD6XsNA2OUfx0F2qsow5kmSYbP8AUAb2NLKgAAQAkKDQ3N8/pQ+j7DcTabXQ0vC5IhISGlWhcAFAWBEQAAoAR16NAhz+sLjnO6kH3OcGyToJa5f+/YsWOp1gUARUFgBAAAKEGjRo2Sr2/eu348rTI2CmomH/nK19dXo0aNsqI8ACgUAiMAAEAJqlWrlm6++eY8bQfS9xiO9bcHqH5QI918882qWbOmFeUBQKEQGAEAAErYpEmT1KVLl9zXxzMOKyMnzXBs9/p99MQTT1hVGgAUCoERAACghAUGBmrOnDkaM2aMfH195ZRTB9L3Go5t4N9Efr5+FlcIAOYQGAEAAEpBYGCgpk+frp9++kkPP/yw/BoYj8tMzlL8jrPWFgcAJhEYAQAASlGtWrV0//3364W5z8g3wMdwzKHVxy2uCgDMITACAABYwC/QV/U6Rxn2HV5zXC6Xy+KKAKBgBEYAAACLNOxR17A99UyaTu8yflYjAHgTgREAAMAiDbpFy+5r/OPX/uVHLK4GAApGYAQAALBIQIi/6nasbdh3cMVROXOcFlcEAPkjMAIAAFioab/6hu3pFzJ1ckuCxdUAQP4IjAAAABZq2L2Ox9NS9//MtlQAZQuBEQAAwEJ+QX5q0K2OYd+h1ceVk5VjcUUA4BmBEQAAwGKetqVmXczW0fWnLK4GADwjMAIAAFisXqco+Yf4Gfbt/4ltqQDKDgIjAACAxXzLine truncated
"text/plain": [ "text/plain": [
"<Figure size 780x1200 with 2 Axes>" "<Figure size 780x1200 with 2 Axes>"
] ]