update notebooks

notebooks/plots.ipynb

@@ -102,7 +102,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 57,
+   "execution_count": 2,
    "metadata": {},
    "outputs": [
     {
@@ -144,7 +144,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 55,
+   "execution_count": 3,
    "metadata": {},
    "outputs": [
     {
@@ -187,7 +187,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 61,
+   "execution_count": 4,
    "metadata": {},
    "outputs": [
     {
@@ -230,12 +230,12 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 2,
+   "execution_count": 5,
    "metadata": {},
    "outputs": [
     {
      "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAABW0AAAPdCAYAAADxjUr8AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8fJSN1AAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOzdd3gUVd/H4c/uplcICUloSegl9Bp670h5UZoCAoIiIAJSFIQAgoBUJYhKEUFFRVTKAwQhdAQEFAi9l9AhBRKy2Z33jyEry24aqSS/22su2bMzZ86cbDaz3z1zRqMoioIQQgghhBBCCCGEEEKIHEGb3Q0QQgghhBBCCCGEEEII8R8JbYUQQgghhBBCCCGEECIHkdBWCCGEEEIIIYQQQgghchAJbYUQQgghhBBCCCGEECIHkdBWCCGEEEIIIYQQQgghchAJbYUQQgghhBBCCCGEECIHkdBWCCGEEEIIIYQQQgghchAJbYUQQgghhBBCCCGEECIHkdBWCCGEEEIIIYQQQgghchAJbYUQL4Xo6GhGjx5Ny5Yt8fLyQqPRMGnSpFRvf/ToUdq1a0exYsVwdHTEw8ODoKAgVq5cmartw8LC0Gg0Vpf9+/e/4FGZ8/f3p3379hlSF8CNGzeYNGkSR48ezbA6c6vGjRvTuHHj7G6GEEIIIZ6zfPlyNBoNly5dSvO2iedvv/zyS4a1J7HOsLCwDKszu0yaNAmNRpPdzRApyMjXXOPGjQkMDEx/o556/PgxkyZNyhW/D0LkRDbZ3QAhhEiNe/fu8dVXX1G5cmU6derEN998k6btHz58SNGiRenRoweFCxfm0aNHrFq1ijfeeINLly4xfvz4VNUzbdo0mjRpYlaWkSc+GenGjRsEBwfj7+9PlSpVsrs5QgghhBBp1q5dO/bt24evr292N0UI8ZzHjx8THBwMIAMghMgEEtoKIV4Kfn5+PHjwAI1Gw927d9Mc2lobSdm+fXsuXrzIV199lerQtlSpUtSpUydN+xZpoygKcXFxODo6ZndThBBCCJHNvLy88PLyyu5miJeAnEMKIXIbmR5BCPFSSJyKIKN5enpiY5P5319duHCB7t27U6hQIezt7fH29qZZs2ZWpy7YtGkT1apVw9HRkbJly7J06VKLdY4fP07Hjh3Jnz8/Dg4OVKlShW+//db0fFhYGDVr1gTgzTffNPVf4pQSaWnPs/r27YuLiwvnzp2jbdu2uLi4ULRoUUaOHMmTJ0/M1r1//z6DBw+mcOHC2NnZUbx4cT766COL9TQaDUOGDOHLL7+kXLly2Nvb8+2335ouh9y2bRtvvfUWBQoUwM3Njd69e/Po0SNu3rzJa6+9Rr58+fD19WXUqFHo9XqzuoODg6lduzYeHh64ublRrVo1lixZgqIoyR6nEEIIIXIGa9MjJF7iffDgQRo0aICTkxPFixfn008/xWg0WtQRFxfHiBEj8PHxwdHRkUaNGnHkyBGzdQ4dOkT37t3x9/fH0dERf39/evToweXLl1NsY2q3TTyW7du388477+Dp6UmBAgXo0qULN27csKj3+++/JygoCBcXF1xcXKhSpQpLliwxW2fr1q00a9YMNzc3nJycqFevHn/++adFXRs2bKBKlSrY29sTEBDAZ599luJxvYiMOH9L7zkkwNmzZ+nZsycFCxbE3t6ecuXKsXDhwlQdw8OHD+nfvz8eHh64uLjQrl07Lly4YHV6tt27d9OsWTNcXV1xcnKibt26bNiwwaLOlM7dE506dYrWrVvj5OSEp6cnb7/9NtHR0alq9507dxg4cCBFixbF3t4eLy8v6tWrx9atWy3WTc3vzpUrV3j99dfN+nD27Nmm9S5dumT6QiU4ONj0eaNv375pbo8QwjoZaSuEyFOMRiNGo5EHDx7w888/s3nzZr744otUb//uu+/SvXt3nJycCAoKYsKECdSvXz/F7dq2bYvBYGDmzJkUK1aMu3fvsnfvXh4+fGi23j///MPIkSMZO3Ys3t7efPPNN/Tv35+SJUvSsGFDAE6fPk3dunUpWLAgCxYsoECBAqxcuZK+ffty69YtRo8eTbVq1Vi2bBlvvvkm48ePp127dgAUKVIkTe2xRq/X88orr9C/f39GjhzJzp07mTJlCu7u7nz88ceA+uGoSZMmnD9/nuDgYCpVqsSuXbuYPn06R48etTiZ/e2339i1axcff/wxPj4+FCxYkIMHDwIwYMAAunTpwo8//siRI0f48MMPSUhI4PTp03Tp0oWBAweydetWZsyYQaFChRgxYoSp3kuXLjFo0CCKFSsGwP79+xk6dCjXr183tVUIIYQQL5+bN2/Sq1cvRo4cycSJE1m7di3jxo2jUKFC9O7d22zdDz/8kGrVqvHNN98QGRnJpEmTaNy4MUeOHKF48eKAes5QpkwZunfvjoeHBxERESxatIiaNWsSHh6Op6dnkm1J67YDBgygXbt2fP/991y9epUPPviA119/nW3btpnW+fjjj5kyZQpdunRh5MiRuLu7c/z4cbMgeOXKlfTu3ZuOHTvy7bffYmtry+LFi2nVqhWbN2+mWbNmAPz555907NiRoKAgfvzxR9M54K1bt9L9c0jKi56/ZcQ5ZHh4OHXr1qVYsWLMnj0bHx8fNm/ezLBhw7h79y4TJ05Mst1Go5EOHTpw6NAhJk2aRLVq1di3bx+tW7e2WHfHjh20aNGCSpUqsWTJEuzt7QkJCaFDhw788MMPdOvWDUjduTvArVu3aNSoEba2toSEhODt7c2qVasYMmRIqvr8jTfe4PDhw3zyySeULl2ahw8fcvjwYe7du2e2Xmp+d+7cuUPdunWJj49nypQp+Pv7s379ekaNGsX58+cJCQnB19eXTZs20bp1a/r378+AAQMATEFuatsjhEiGIoQQL5k7d+4ogDJx4sQ0bzto0CAFUADFzs5OCQkJSdV2hw8fVt577z1l7dq1ys6dO5WlS5cq5cqVU3Q6nbJp06Zkt717964CKPPmzUt2PT8/P8XBwUG5fPmyqSw2Nlbx8PBQBg0aZCrr3r27Ym9vr1y5csVs+zZt2ihOTk7Kw4cPFUVRlIMHDyqAsmzZshdqjzV9+vRRAOWnn34yK2/btq1SpkwZ0+Mvv/zS6nozZsxQAGXLli2mMkBxd3dX7t+/b7busmXLFEAZOnSoWXmnTp0UQJkzZ45ZeZUqVZRq1aol2XaDwaDo9Xpl8uTJSoECBRSj0Wh6rlGjRkqjRo2SP3ghhBBCZLnE84GLFy+ayho1aqQAyl9//WW2bvny5ZVWrVqZHm/fvl0BlGrVqpn93b906ZJia2urDBgwIMn9JiQkKDExMYqzs7Myf/58izq3b9+e5m0Tj2Xw4MFm68+cOVMBlIiICEVRFOXChQuKTqdTevXqleQ+Hj16pHh4eCgdOnQwKzcYDErlypWVWrVqmcpq166tFCpUSImNjTWVRUVFKR4eHkpGRwLpPX/LiHPIVq1aKUWKFFEiIyPNyocMGaI4ODhYrP+sDRs2KICyaNEis/Lp06dbfP6oU6eOUrBgQSU6OtpUlpCQoAQGBipFihQxveZSe+4+ZswYRaPRKEePHjVbr0WLFim+5hRFUVxcXJThw4cnu05qf3fGjh1rdb133nlH0Wg0yunTpxVFSf5zWWraI4RInkyPIITIUz788EMOHjzIhg0b6NevH0OGDEnV5WFVq1Zl3rx5dOrUiQYNGvDmm2+yd+9efH19Td+OJ8XDw4MSJUowa9Ys5syZw5EjR6xeugdQpUoV06hQAAcHB0qXLm02qmLbtm00a9aMokWLmm3bt29fHj9+zL59+zKsPdZoNBo6dOhgVlapUiWLNjo7O9O1a1eLNgIWl+01bdqU/PnzW91f+/btzR6XK1cOwDR6+Nny5y9D3LZtG82bN8fd3R2dToetrS0ff/wx9+7d4/bt2ykcqRBCCCFyKh8fH2rVqmVW9vz5SKKePXuaTbPl5+dH3bp12b59u6ksJiaGMWPGULJkSWxsbLCxscHFxYVHjx5x8uTJZNuS1m1feeUVi3YDpraHhoZiMBh49913k9zn3r17uX//Pn369CEhIcG0GI1GWrduzcGDB3n06BGPHj3i4MGDdOnSBQcHB9P2rq6uFudz1iiKYlZ/QkJCitvAi5+/pfccMi4ujj///JPOnTvj5ORk1u62bdsSFxfH/v37k2z3jh07AHjttdfMynv06GH2+NGjR/z111907doVFxcXU7lOp+ONN97g2rVrnD592nRMqTl33759OxUqVKBy5cpm6/Xs2TPJ9j6rVq1aLF++nKlTp7J//36LaScSpeZ3Z9u2bZQvX95ivb59+6Ioitmo8PS2RwiRNAlthRB5SrFixahRowZt27Zl0aJFDBw4kHHjxnHnzp0015UvXz7at2/Pv//+S2xsbJLraTQa/vzzT1q1asXMmTOpVq0aXl5eDBs2zGKOqgIFClhsb29vb1b/vXv3rN5BuVChQqbnk5OW9ljj5ORkdtKf2Ma4uDizNvr4+FjMQ1ywYEFsbGws2pjcHaE9PDzMHtvZ2SVZ/mwbDhw4QMuWLQH4+uuv2bNnDwcPHuSjjz4CSPZnJoQQQoicLTXnTIl8fHyslj17PtKzZ0+++OILBgwYwObNmzlw4AAHDx7Ey8srxXOGtG77fNvt7e2B/85NEs9LE6e1siZxaoOuXbtia2trtsyYMQNFUbh//z4PHjzAaDQm2QcpSZx24dklNV70/C2955D37t0jISGBzz//3KLdbdu2BeDu3btJtvvevXvY2NhYtNPb29vs8YMHD1AUJVXn5Kk9d0889uel5ucEsHr1avr06cM333xDUFAQHh4e9O7dm5s3b5qtlxWfN9LSHiFE0mROWyFEnlarVi2+/PJLLly48EJ3Jlae3tAqpZuk+fn5mW4ccebMGX766ScmTZpEfHw8X375ZZr2WaBAASIiIizKE29gkdyca5nRnqTa+Ndff6Eoilnf3L5Line truncated
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAABW0AAAPdCAYAAADxjUr8AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/GU6VOAAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOzdd3gUVd/H4c/uplcICUloSegl9Bp670h5UZoCAoIiIAJSFIQAgoBUJYhKEUFFRVTKAwQhdAQEFAi9l9AhBRKy2Z33jyEry24aqSS/22su2bMzZ86cbDaz3z1zRqMoioIQQgghhBBCCCGEEEKIHEGb3Q0QQgghhBBCCCGEEEII8R8JbYUQQgghhBBCCCGEECIHkdBWCCGEEEIIIYQQQgghchAJbYUQQgghhBBCCCGEECIHkdBWCCGEEEIIIYQQQgghchAJbYUQQgghhBBCCCGEECIHkdBWCCGEEEIIIYQQQgghchAJbYUQQgghhBBCCCGEECIHkdBWCCGEEEIIIYQQQgghchAJbYUQL4Xo6GhGjx5Ny5Yt8fLyQqPRMGnSpFRvf/ToUdq1a0exYsVwdHTEw8ODoKAgVq5cmartw8LC0Gg0Vpf9+/e/4FGZ8/f3p3379hlSF8CNGzeYNGkSR48ezbA6c6vGjRvTuHHj7G6GEEIIIZ6zfPlyNBoNly5dSvO2iedvv/zyS4a1J7HOsLCwDKszu0yaNAmNRpPdzRApyMjXXOPGjQkMDEx/o556/PgxkyZNyhW/D0LkRDbZ3QAhhEiNe/fu8dVXX1G5cmU6derEN998k6btHz58SNGiRenRoweFCxfm0aNHrFq1ijfeeINLly4xfvz4VNUzbdo0mjRpYlaWkSc+GenGjRsEBwfj7+9PlSpVsrs5QgghhBBp1q5dO/bt24evr292N0UI8ZzHjx8THBwMIAMghMgEEtoKIV4Kfn5+PHjwAI1Gw927d9Mc2lobSdm+fXsuXrzIV199lerQtlSpUtSpUydN+xZpoygKcXFxODo6ZndThBBCCJHNvLy88PLyyu5miJeAnEMKIXIbmR5BCPFSSJyKIKN5enpiY5P5319duHCB7t27U6hQIezt7fH29qZZs2ZWpy7YtGkT1apVw9HRkbJly7J06VKLdY4fP07Hjh3Jnz8/Dg4OVKlShW+//db0fFhYGDVr1gTgzTffNPVf4pQSaWnPs/r27YuLiwvnzp2jbdu2uLi4ULRoUUaOHMmTJ0/M1r1//z6DBw+mcOHC2NnZUbx4cT766COL9TQaDUOGDOHLL7+kXLly2Nvb8+2335ouh9y2bRtvvfUWBQoUwM3Njd69e/Po0SNu3rzJa6+9Rr58+fD19WXUqFHo9XqzuoODg6lduzYeHh64ublRrVo1lixZgqIoyR6nEEIIIXIGa9MjJF7iffDgQRo0aICTkxPFixfn008/xWg0WtQRFxfHiBEj8PHxwdHRkUaNGnHkyBGzdQ4dOkT37t3x9/fH0dERf39/evToweXLl1NsY2q3TTyW7du388477+Dp6UmBAgXo0qULN27csKj3+++/JygoCBcXF1xcXKhSpQpLliwxW2fr1q00a9YMNzc3nJycqFevHn/++adFXRs2bKBKlSrY29sTEBDAZ599luJxvYiMOH9L7zkkwNmzZ+nZsycFCxbE3t6ecuXKsXDhwlQdw8OHD+nfvz8eHh64uLjQrl07Lly4YHV6tt27d9OsWTNcXV1xcnKibt26bNiwwaLOlM7dE506dYrWrVvj5OSEp6cnb7/9NtHR0alq9507dxg4cCBFixbF3t4eLy8v6tWrx9atWy3WTc3vzpUrV3j99dfN+nD27Nmm9S5dumT6QiU4ONj0eaNv375pbo8QwjoZaSuEyFOMRiNGo5EHDx7w888/s3nzZr744otUb//uu+/SvXt3nJycCAoKYsKECdSvXz/F7dq2bYvBYGDmzJkUK1aMu3fvsnfvXh4+fGi23j///MPIkSMZO3Ys3t7efPPNN/Tv35+SJUvSsGFDAE6fPk3dunUpWLAgCxYsoECBAqxcuZK+ffty69YtRo8eTbVq1Vi2bBlvvvkm48ePp127dgAUKVIkTe2xRq/X88orr9C/f39GjhzJzp07mTJlCu7u7nz88ceA+uGoSZMmnD9/nuDgYCpVqsSuXbuYPn06R48etTiZ/e2339i1axcff/wxPj4+FCxYkIMHDwIwYMAAunTpwo8//siRI0f48MMPSUhI4PTp03Tp0oWBAweydetWZsyYQaFChRgxYoSp3kuXLjFo0CCKFSsGwP79+xk6dCjXr183tVUIIYQQL5+bN2/Sq1cvRo4cycSJE1m7di3jxo2jUKFC9O7d22zdDz/8kGrVqvHNN98QGRnJpEmTaNy4MUeOHKF48eKAes5QpkwZunfvjoeHBxERESxatIiaNWsSHh6Op6dnkm1J67YDBgygXbt2fP/991y9epUPPviA119/nW3btpnW+fjjj5kyZQpdunRh5MiRuLu7c/z4cbMgeOXKlfTu3ZuOHTvy7bffYmtry+LFi2nVqhWbN2+mWbNmAPz555907NiRoKAgfvzxR9M54K1bt9L9c0jKi56/ZcQ5ZHh4OHXr1qVYsWLMnj0bHx8fNm/ezLBhw7h79y4TJ05Mst1Go5EOHTpw6NAhJk2aRLVq1di3bx+tW7e2WHfHjh20aNGCSpUqsWTJEuzt7QkJCaFDhw788MMPdOvWDUjduTvArVu3aNSoEba2toSEhODt7c2qVasYMmRIqvr8jTfe4PDhw3zyySeULl2ahw8fcvjwYe7du2e2Xmp+d+7cuUPdunWJj49nypQp+Pv7s379ekaNGsX58+cJCQnB19eXTZs20bp1a/r378+AAQMATEFuatsjhEiGIoQQL5k7d+4ogDJx4sQ0bzto0CAFUADFzs5OCQkJSdV2hw8fVt577z1l7dq1ys6dO5WlS5cq5cqVU3Q6nbJp06Zkt717964CKPPmzUt2PT8/P8XBwUG5fPmyqSw2Nlbx8PBQBg0aZCrr3r27Ym9vr1y5csVs+zZt2ihOTk7Kw4cPFUVRlIMHDyqAsmzZshdqjzV9+vRRAOWnn34yK2/btq1SpkwZ0+Mvv/zS6nozZsxQAGXLli2mMkBxd3dX7t+/b7busmXLFEAZOnSoWXmnTp0UQJkzZ45ZeZUqVZRq1aol2XaDwaDo9Xpl8uTJSoECBRSj0Wh6rlGjRkqjRo2SP3ghhBBCZLnE84GLFy+ayho1aqQAyl9//WW2bvny5ZVWrVqZHm/fvl0BlGrVqpn93b906ZJia2urDBgwIMn9JiQkKDExMYqzs7Myf/58izq3b9+e5m0Tj2Xw4MFm68+cOVMBlIiICEVRFOXChQuKTqdTevXqleQ+Hj16pHh4eCgdOnQwKzcYDErlypWVWrVqmcpq166tFCpUSImNjTWVRUVFKR4eHkpGRwLpPX/LiHPIVq1aKUWKFFEiIyPNyocMGaI4ODhYrP+sDRs2KICyaNEis/Lp06dbfP6oU6eOUrBgQSU6OtpUlpCQoAQGBipFihQxveZSe+4+ZswYRaPRKEePHjVbr0WLFim+5hRFUVxcXJThw4cnu05qf3fGjh1rdb133nlH0Wg0yunTpxVFSf5zWWraI4RInkyPIITIUz788EMOHjzIhg0b6NevH0OGDEnV5WFVq1Zl3rx5dOrUiQYNGvDmm2+yd+9efH19Td+OJ8XDw4MSJUowa9Ys5syZw5EjR6xeugdQpUoV06hQAAcHB0qXLm02qmLbtm00a9aMokWLmm3bt29fHj9+zL59+zKsPdZoNBo6dOhgVlapUiWLNjo7O9O1a1eLNgIWl+01bdqU/PnzW91f+/btzR6XK1cOwDR6+Nny5y9D3LZtG82bN8fd3R2dToetrS0ff/wx9+7d4/bt2ykcqRBCCCFyKh8fH2rVqmVW9vz5SKKePXuaTbPl5+dH3bp12b59u6ksJiaGMWPGULJkSWxsbLCxscHFxYVHjx5x8uTJZNuS1m1feeUVi3YDpraHhoZiMBh49913k9zn3r17uX//Pn369CEhIcG0GI1GWrduzcGDB3n06BGPHj3i4MGDdOnSBQcHB9P2rq6uFudz1iiKYlZ/QkJCitvAi5+/pfccMi4ujj///JPOnTvj5ORk1u62bdsSFxfH/v37k2z3jh07AHjttdfMynv06GH2+NGjR/z111907doVFxcXU7lOp+ONN97g2rVrnD592nRMqTl33759OxUqVKBy5cpm6/Xs2TPJ9j6rVq1aLF++nKlTp7J//36LaScSpeZ3Z9u2bZQvX95ivb59+6Ioitmo8PS2RwiRNAlthRB5SrFixahRowZt27Zl0aJFDBw4kHHjxnHnzp0015UvXz7at2/Pv//+S2xsbJLraTQa/vzzT1q1asXMmTOpVq0aXl5eDBs2zGKOqgIFClhsb29vb1b/vXv3rN5BuVChQqbnk5OW9ljj5ORkdtKf2Ma4uDizNvr4+FjMQ1ywYEFsbGws2pjcHaE9PDzMHtvZ2SVZ/mwbDhw4QMuWLQH4+uuv2bNnDwcPHuSjjz4CSPZnJoQQQoicLTXnTIl8fHyslj17PtKzZ0+++OILBgwYwObNmzlw4AAHDx7Ey8srxXOGtG77fNvt7e2B/85NEs9LE6e1siZxaoOuXbtia2trtsyYMQNFUbh//z4PHjzAaDQm2QcpSZx24dklNV70/C2955D37t0jISGBzz//3KLdbdu2BeDu3btJtvvevXvY2NhYtNPb29vs8YMHD1AUJVXn5Kk9d0889uel5ucEsHr1avr06cM333xDUFAQHh4e9O7dm5s3b5qtlxWfN9LSHiFE0mROWyFEnlarVi2+/PJLLly48EJ3Jlae3tAqpZuk+fn5mW4ccebMGX766ScmTZpEfHw8X375ZZr2WaBAASIiIizKE29gkdyca5nRnqTa+Ndff6Eoilnf3L5Line truncated
       "text/plain": [
        "<Figure size 1400x1000 with 6 Axes>"
       ]
@@ -270,6 +270,50 @@
     "plt.tight_layout()\n",
     "plt.show()"
    ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 8,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnYAAAHWCAYAAAD6oMSKAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/GU6VOAAAACXBIWXMAAA9hAAAPYQGoP6dpAABzXUlEQVR4nO3deVxU5eIG8OfMwrAPAsqwiaiksrhvaJpW4lJuLVqZmUtlm3lp05s3t7pq91YuN/m1qKRWWplmXTKxcruaprlj5pbKLigM6zDL+f0xMDDOgDM4MDA+3z7zwTnnPee8ox14eN/zvq8giqIIIiIiImr2JM6uABERERE5BoMdERERkYtgsCMiIiJyEQx2RERERC6CwY6IiIjIRTDYEREREbkIBjsiIiIiF8FgR0REROQiGOyIiIiIXASDHRE1muTkZAiCYHrJZDKEhYVh8uTJyMjIuOXz79y503Tu5ORkq2XuvvtuCIKANm3a3PL1AKBNmzZmn6nmq7i42CHXcLR58+ZBEATk5eXdtOygQYMwaNCghq8UETmEzNkVIKLbz5o1a9CxY0eUlZVh9+7dWLRoEXbt2oUTJ07Ay8vrls/v4+ODVatW4cknnzTbfvHiRezcuRO+vr63fI2a+vfvj3//+98W2z09PR16HSKim2GwI6JGFxsbi549ewIABg8eDL1ej4ULF2LLli2YMGECysrKkJeXh/Dw8Hqdf/z48fjkk09w9uxZREVFmbavXr0aoaGhiIuLQ1pamtkxZ8+eRfv27SEIgt3X8/PzQ9++fetVVyIiR2JXLBE5XVUounTpEgAgJycHERERGDhwIJKSkmzqMqxpyJAhCA8Px+rVq03bDAYDPv30U0yaNAkSieW3vqeeegoRERF47bXXcOTIkVv4NA3n888/R3x8PLy9veHt7Y2uXbti1apVpv2pqakYPXo0wsLC4O7ujvbt2+OZZ56p9e/vypUreOCBB+Dr6wulUonHH38cV69evWk9Kioq8NZbb6Fjx45QKBRo2bIlJk+ebNOxRNSwGOyIyOnOnTsHAGjZsiUAIDw8HJs3b0ZoaCheffVVBAcHY8SIEVi3bh2Kiopuej6JRIInn3wSa9euhV6vBwBs374d6enpmDx5stVjVq5ciccffxybNm1C9+7d0alTJyxYsABnz5696fVEUYROpzN7GQwGWz++Td58801MmDABISEhSE5OxubNmzFp0iRTGAaA8+fPIz4+HklJSdi+fTvefPNNHDhwAHfeeSe0Wq3FOceOHYv27dvj66+/xrx587BlyxYMHTrUatkqBoMBo0ePxuLFi/HYY4/hv//9LxYvXozU1FQMGjQIZWVlDv3cRGQnkYiokaxZs0YEIP7666+iVqsVi4qKxO+//15s2bKl6OPjI2ZnZ1scU1JSIm7cuFEcO3as6O7uLnp4eIgPP/yw+M0334jl5eVmZX/55RcRgPjVV1+JFy5cEAVBEL///ntRFEXx4YcfFgcNGiSKoijed999YkRERK31PHjwoJiYmCiGhYWJAMQePXqI7777rpienm5RNiIiQgRg8XrjjTdu4W/K3IULF0SpVCpOmDDB5mMMBoOo1WrFS5cuiQDEb7/91rRv7ty5IgDxb3/7m9kxn332mQhAXL9+vWnbXXfdJd51112m91988YUIQNy0aZPZsb/99psIQFy5cqWdn46IHIktdkTU6Pr27Qu5XA4fHx/cf//9UKlU+OGHHxAUFGRR1tPTE+PGjcM333yDnJwcfPjhhygpKcG4ceMQFBSE//3vf1avERkZiUGDBmH16tXIz8/Ht99+iylTpthUv169euHdd9/F5cuXsXv3bvTt2xdLlixB69at8fLLL1uUv/POO/Hbb7+ZvZ577rk6r6HX621u4UtNTYVer8fzzz9f5zlzc3Mxffp0hIeHQyaTQS6XIyIiAgBw+vRpi/ITJkwwez9u3DjIZDL88ssvtV7j+++/h5+fH0aOHGlW/65du0KlUmHnzp111pGIGhYHTxBRo1u7di06deoEmUyGoKAgBAcH23RcSUkJCgoKUFhYCL1eDx8fH7i5udVafurUqZg8eTLee+89eHh44KGHHrKrnlqtFoWFhSgoKEBZWRnc3Nzg4+NjUU6pVJoGg9iqXbt2Zt2oc+fOxbx586yWrXp2LSwsrNbzGQwGJCQkIDMzE//4xz8QFxcHLy8vGAwG9O3b12oXqUqlMnsvk8kQEBCA/Pz8Wq+Tk5ODgoKCWv/e7X0ekogci8GOiBpdp06dbA5CeXl52LRpEzZs2IDdu3ejRYsWePDBB/H2229j4MCBdY5ifeCBB/D8889j8eLFeOqpp+Dh4XHT6+l0Ovz000/YuHEjNm/ejKKiItx9991YtmwZHnzwQYdNlfLdd99Bo9GY3oeEhNRaturZw/T09FpHCp88eRLHjh1DcnIyJk2aZNpe9fyiNdnZ2QgNDTW91+l0yM/PR0BAQK3HBAYGIiAgANu2bbO631rwJaLGw2BHRE1ORUUF1q9fj40bN+Lnn3+GQqHAqFGj8O2332Lo0KGQy+U2ncfDwwNvvvkmdu/ejWeffbbOsvv378enn36KTZs2IS8vD3379sX8+fMxfvx4q13EtyouLs7msgkJCZBKpUhKSkJ8fLzVMlUBV6FQmG3/8MMPaz3vZ599hh49epjef/nll9DpdHVOSHz//fdjw4YN0Ov16NOnj82fgYgaB4MdETU5mZmZmD59OoYOHYq1a9di9OjR9Z7sNzExEYmJiTctN3v2bOTl5WHmzJl47LHHEBkZWa/rNYQ2bdrg73//OxYuXIiysjI8+uijUCqVSEtLQ15eHubPn4+OHTuiXbt2mDVrFkRRhL+/P7777jukpqbWet5vvvkGMpkMQ4YMwalTp/CPf/wDXbp0wbhx42o95pFHHsFnn32GESNG4KWXXkLv3r0hl8uRnp6OX375BaNHj8bYsWMb4q+BiGzAYEdETU5wcDCys7Ph7+/faNf84osvbH7WzxkWLFiAqKgorFixAhMmTIBMJkNUVBRmzJgBAJDL5fjuu+/w0ksv4ZlnnoFMJsO9996LHTt2oHXr1lbP+c0332DevHlISkqCIAgYOXIkli5dWudzi1KpFFu3bsWyZcuwbt06LFq0yLQ03F133WVXSyQROZ4giqLo7EoQERER0a3jdCdERERELoLBjoiIiMhFMNgRERERuQgGOyIiIiIXwWBHRERE5CIY7IiIiIhcBOexs8JgMCAzMxM+Pj51LldERERE1NBEUURRURFCQkIgkdTdJsdgZ0VmZmat6zESEREROcOVK1cQFhZWZxkGOyuqFrG+cuWKwxb8tpVWq8X27duRkJBg83qYRNQ08X4mch3OvJ/VajXCw8NN+aQuDHZWVHW/+vr6OiXYeXp6wtfXlz8IiJo53s9ErqMp3M+2PB7GwRNERERELoLBjoiIiMhFMNgRERERuQgGOyIiIiIXwWBHRERE5CIY7IiIiIhcBIMdERERkYtgsCMiIiJyEQx2RERERC6CwY6IiIjIRTDYERE1AL1Bj0M5h3Cs4hgO5RyC3qB3dpWI6DbAtWKJiBxsx6UdWHxwMXJKcwAAX/30FYI8gzCr9yzcG3Gvk2tHRK6MLXZERA6049IOJO5MNIW6KrmluUjcmYgdl3Y4qWZEdDtgsCMichC9QY/FBxdDhGixr2rbkoNL2C1LRA2GXbFERDehM+hQqClEgaYA18uvG79qrqOg3PxrelG6RUtdTSJEZJdmY87/5qBnUE8Eewcj2Mv4cpe5N+InIiJXxWBHRLcVg2hAUUVRdUCrI6hV7VdXqB1ah+8vfI/vL3xvts3f3R8hXiGmsBfiHQKVlwohXiEI8Q6Br5svBEFwaD2IyPUw2BFRsyWKIoq1xRZBrObXG/9coCmAQTTYfS0BApQKJfwUfsaXux9aKFpUf1X4Ia8sD8uPLL/puQaHD4bOoENWSRYyizNRqivFtfJruFZ+DSfzT1o9xlPmaWzd8w62CIDBXsFo6dESUonU7s9FRK6FwY6ImgRRFFGmK6u99cxKeCsoL4BO1NXrej5yH7Nw5qe4Iajd8NXXzfemwUlv0GPjmY3ILc21+pydAAFBnkF4f9D7pnOJogh1hdoU8rJKspBVnIXMkkxkFWchqyQL+eX5KNWV4nzheZwvPG/12jJBhiCvIIvWvqogqPJSsbuX6DbAYEdEDUKj1+B6+XUUagptDmoavaZe1/KQeVgPZQo/tHC3/Kp0U0IulTv4EwNSiRSzes9C4s6/QRBFiDW6TgVRBATg9d6vmwVEQTC2BCoVSnT072j1vOW6cmSXZCOzJNP4tTIAVn3NKcmBTtQhozgDGcUZQC2P+Vnr7q16xo/dvUSugcGOiG5Ka9AaA9qNz6XV8Xxama6sXtdyk7ihhXuL6iB2Q2Cr6gqt2u+n8GtSLVH3lpTivZw8LA7wQ46s+ltskF6P1/MLcG9Jqd3ndJe5o42yDdoo21jdrzfocbXsqin01Wzts7e711prX1UQZHcvUdPHYEd0m9Eb9FBXqG/ailZzW5G2qF7Xkgmy6m7OWoLajS1sHjKP5ttqZNAD217HvaWlGFxait/dFbgqlaKlXo/u5RpIAeD7mYDcA5ApAEEKSGSVL4nxq2mb1PgyK3PjNuNXqUQKlZcKKi8VurbqalGtm3X3ZpZk4lr5NZTqSnGu4BzOFZyz+vFu7O6tau2rGQAVUkVD/g0TOYdBD+HSXoRe2w/hki/QdqDx/muCGOyImrGqEZ62tKJVDRwo1BRaff7rZiSCBEo35U27OWuGN2+5d/MNafVxaR+gzgQASAH0KrfStVyaD3z2kIMvLNQZ/gSJDEpBAqVEho5WyyhQLglFtgTIhAFZggFZgh5Z0CNTrEAWtMgRtbZ190oUCJF5I1jmjRC5D4LdfBEsVyLETYlghR98ZZ4QpPIaIVZqCqjV9ZFUf54bt90YfG0Jw7fT/4PkeGlbgW2Line truncated
+      "text/plain": [
+       "<Figure size 640x480 with 1 Axes>"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    }
+   ],
+   "source": [
+    "# Create subplots\n",
+    "fig, ax = plt.subplots()\n",
+    "\n",
+    "# Loop through each measurement type\n",
+    "model = \"P>M>F\"\n",
+    "data = cable_data\n",
+    "imgclass = \"cable\"\n",
+    "\n",
+    "# Plot both bottle and cable data\n",
+    "ax.plot([1,3,5], data[\"CAML\"].get(\"2 ways only detect if faulty or not\", []), marker='o', label=f'CAML', linestyle='-')\n",
+    "ax.plot([1,3,5], data[\"P>M>F\"].get(\"2 ways only detect if faulty or not\", []), marker='o', label=f'P>M>F', linestyle='-')\n",
+    "ax.plot([1,3,5], data[\"ResNet50\"].get(\"2 ways only detect if faulty or not\", []), marker='o', label=f'ResNet50', linestyle='-')\n",
+    "\n",
+    "\n",
+    "ax.set_title(f'{model} - {imgclass}')\n",
+    "ax.set_xlabel(\"Shots per class\")\n",
+    "ax.set_ylabel(\"Accuracy\")\n",
+    "ax.legend()\n",
+    "ax.grid(True)\n",
+    "\n",
+    "ax.set(xlim=(0.75, 5.25), xticks=[1,3,5])\n",
+    "\n",
+    "# Adjust layout\n",
+    "plt.tight_layout()\n",
+    "plt.savefig(f\"comparison-2way-{imgclass}.png\")"
+   ]
   }
  ],
  "metadata": {