Updated to Python 3.6

This commit is contained in:
Roger Labbe committed 2020-10-13 14:50:11 -07:00
1 parent 91be73db4e
commit 9e9121bdfe
1 file changed
+54 -47
+54 -47
View File
@@ -18,17 +18,16 @@
},
{
"cell_type": "code",
"execution_count": 1,
"execution_count": 5,
"metadata": {},
"outputs": [],
"source": [
"from __future__ import division, print_function\n",
"%matplotlib inline"
]
},
{
"cell_type": "code",
"execution_count": 2,
"execution_count": 6,
"metadata": {
"scrolled": true
},
@@ -40,7 +39,7 @@
" <style>\n",
" .output_wrapper, .output {\n",
" height:auto !important;\n",
" max-height:100000px; \n",
" max-height:100000px;\n",
" }\n",
" .output_scroll {\n",
" box-shadow:none !important;\n",
@@ -53,7 +52,7 @@
"<IPython.core.display.HTML object>"
]
},
"execution_count": 2,
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
@@ -124,17 +123,19 @@
},
{
"cell_type": "code",
"execution_count": 3,
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAA1wAAAGACAYAAAC5hmO5AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nOzdeVxU570/8M8w7CAqCIIbuIIrWdwSYzRG3FeMu5xTe5ve/pqksU2apomJWqPNbWNftzXLfdm09RzcYtziTtSYzSWxWlFRURHiggiMCLIMzHJ+fyAnDIwICswz8Hm/Xn3d3PPMeXxmPnOOfuc85zkGTdM0EBERERERUb3zcPUAiIiIiIiImioWXERERERERA2EBRcREREREVEDYcFFRERERETUQFhwERERERERNRAWXERERERERA3E09UDcDW73Q673e6wzWAwwGAwuGhEREREREQkIk3TUPWpWh4eHvDwuPd1LBZcdjuKiopcPQwiIiIiInJDAQEBNRZcnFJIRERERETUQFhwERERERERNRAWXERERERERA2k2d/D5WxxjPvNwyQiIiIioubH2foP91tsjwWXkw/ofiuNEBERERERAfcvuFhVEBERERERNZBmf4VLVBcuXIDVaoWnpyd69Ojh6uHQXcxFPMxETMxFTMxFPMxETMxFTO6aCwsuQZWUlMBiscDLy8vVQ6FKmIt4mImYmIuYmIt4mImYmIuY3DUXTikUlMFg0P9H4mAu4mEmYmIuYmIu4mEmYmIuYnLXXAyapmmuHoQr2e123Llzx2FbixYtuGgGERERERE5eJDagVUFERERERG5BXe8VsR7uIiIiIiISFiFpVa8l5SK/eduwmLT4GU0YGTPtnh1dDQCfcQvZzilkFMKiYiIiIiEVFhqRfyHh3ApuxD2SlWLhwHoFhaILb8c0qhF14PUDuKXhM1UTk4ObDYbjEYjQkNDXT0cuou5iIeZiIm5iIm5iIeZiIm5iOO9pNRqxRYA2DXgUnYhViSlYtGk3q4ZXC2x4BJUZmamvuwlD3RxMBfxMBMxMRcxMRfxMBMxMRdx7D93s1qxVcGuAfvO3RS+4OK8OSIiIiIiEo6mabDYar77yWrThF9Ig1e4BBUZGQm73c57yQTDXMTDTMTEXMTEXMTDTMTEXMRgMBjgZaz5mVueRvGfy8WCS1CtWrVy9RDICeYiHmYiJuYiJuYiHmYiJuYijpE920I9kuF0WqGHAYjr2bbRx1RXLNuJiIiIiEhIr46ORrewQHhUuYhVsUrhK6OjXTOwOuCy8FwWnoiIiIhIWIWlVqxISsW+czdhtWnwNBoQ17MtXnHBc7gepHZgwSVowWW1WqFpGgwGAzw9OfNTFMxFPMxETMxFTMxFPMxETMxFTFarVb+3zlW58DlcTUhKSoq+HGlsbKyrh0N3MRfxMBMxMRcxMRfxMBMxMRcxuWsunDdHRERERETUQHiFS1BBQUGwWq28jC0Y5iIeZiIm5iIm5iIeZiIm5iImd82F93AJeg8XERERERGJ5UFqB1YVREREREREDYQFFxERERERUQNhwUVERERERNRA3OuOs2bk8uXL+k2BXbp0cfVw6C7mIh5mIibmIibmIh5mIibmIiZ3zYUFl6Du3LmjP2eAxMFcxMNMxMRcxMRcxMNMxMRcxOSuuXBKIRERERERUQPhsvCCLgtvs9n0/zYajS4cCVXGXMTDTMTEXMTEXMTDTMTEXMQkQi4PUjtwSqGgeHCLibmIh5mIibmIibmIh5mIibmIyV1z4ZRCIiIiIiKiBsKCi4iIiIiIqIG4vOD64osv8NOf/hQxMTEICAhA+/btMXnyZBw/frzaa0+cOIGRI0ciMDAQrVq1Qnx8PC5fvuy035UrVyImJgY+Pj7o3LkzlixZAovF0tBvp97k5eXBZDIhLy/P1UOhSpiLeJiJmJiLmJiLeJiJmJiLmEwmk1vm4vJ7uD766COYTCa8/PLL6NWrF3JycrBixQoMHjwYSUlJGDFiBADg/PnzGD58OB555BFs3LgRZrMZb7/9NoYOHYqTJ08iNDRU73PZsmV466238Prrr2PUqFE4duwYFi5ciOvXr2PVqlWueqt1cuXKFX3Zy9atW7t6OHQXcxEPMxETcxETcxEPMxETcxFHVlYW1q5dC1VVMXv2bIwcOdLtcnF5wfXBBx8gLCzMYduYMWPQrVs3LF++XC+43n77bfj4+GDnzp0ICgoCADz++OPo3r073nvvPfzP//wPgPLK95133sHzzz+P5cuXAwCGDx8Oi8WChQsXYsGCBejVq1cjvkMiIiIiIqqtkpISbN++HaqqIikpSV+dsGXLlhg5cqSLR1d3Li+4qhZbABAYGIhevXrh6tWrAACr1YqdO3dCkiS92AKAyMhIPPPMM9i6datecO3duxdmsxnz58936HP+/Pl48803sW3bNrcouNq3bw+73e7y5enJEXMRDzMRE3MRE3MRDzMRE3NpfJqm4dChQ1BVFRs3bkR+fr7e9sQTT0CSJIwcORJBQUFul4vLCy5n8vPzceLECf3qVlpaGkpKStCvX79qr+3Xrx/27dsHs9kMX19fnDlzBgDQt29fh9dFRESgTZs2ervo2rRp4+ohkBPMRTzMREzMRUzMRTzMREzMpfGkp6dDVVUkJiYiLS1N396pUydIkoSEhAT06NHDhSN8eEIWXC+88AKKiorw5ptvAiifJggAwcHB1V4bHBwMTdOQl5eHiIgImEwm+Pj4ICAgwOlrK/qqSUpKCiIjIx2uppWWluL8+fMAgNatW6NTp04O+1y8eBHFxcUAgNjYWIe23NxcXL9+HUD5l6fynFObzaYXgS1atECXLl0c9k1PT0dBQQEAoHfv3vD0/DGy27dv44cffgAAtGvXzuE+NgA4deoUNE2Dn59ftS/q1atXcevWLQBAdHQ0fH199bbCwkL9Cx8WFoaIiAiHfc+ePavPa656tfDGjRvIzs4GAHTt2hWBgYF6m9lsRmpqKoDyLDp27Oiw74ULF1BSUgKDwVCtuM7JyUFmZiaA8iubrVq10tusVitSUlIAAEFBQejcubPDvpcvX9YfUNenTx+HZzjk5eXhypUrAMp/zap6gk1OTgYA+Pv7o3v37g5tV65c0W/arFigpUJBQQHS09MBAG3btkV4eLjDvikpKbBarfDx8UFMTIxDW2ZmJnJycgAA3bp1c/guFxcX4+LFiwCAkJAQdOjQwWHf1NRUmM1mGI1G9OnTx6EtOzsbN27cAABERUWhZcuWeltZWRnOnTsHoPxyfVRUlMO+aWlpKCwsBFD+Y0blX5ZMJhOuXbsGAOjQoQNCQkL0NrvdjtOnTwMov3LdtWtXh34zMjL0X7B69uwJb29vvS0/Px8ZGRkAyn8wqXo1/MyZM7DZbPD19UV0dLRD27Vr1/RjvXv37vD399fbioqKcOnSJQBAaGgo2rVr57Dv+fPnUVpaCk9PT/Tu3duhLSsrCzdv3gQAdO7cmecI8BzBc0Q5niPK8RxRjueIH/EcUa7qOaKwsBCffvopFEXBN998o78uICAA06dPhyzLePrppwEAp0+fRnJysjDniIrjtC6EK7jeeustrF27FitXrsTjjz/u0GYwGO65X+W22r7uXqxWKzRNc9imaZq+yqHVanW6z71WQbTb7Xqb3W6v1l7bfquOqXK/lZ+8XblfTdPg5eVVrc1msz1Uv/d6r5X7rfpeK3+GzvqteK/OMqptv6Jkc7/3arFYYLVanV4Srymb2vRrsVicttXUb8W+9+vXmdp+hs72rzymqu73XsvKymC3250+BLE+PkNnn1Ft+xXle8hzRPV+RcmG54jq/fIcwXNE1X55jmha5wiz2YzvvvsO7777LrZt2waz2Qyg/N/mAwcOxLhx4yBJkkOxVvnzFekcUVdCFVxLlizBO++8g2XLluHFF1/Ut1f8Gubs6tStW7dgMBj0XypCQkJgNptRXFzsUI1WvLZqEeeMp6dntYPVYDDoJ5zKvw5V3sfZCQkAPDw89DZnB0Zt+606psr9OvuyeHl5QdM0p/0ajcaH6rfy/71Xv1Xfa+XP0Fm/Fe/V2Ymytv2Kks393mvF+7zfZ/gg/dpsNqdtNfVbsW9N/T7sZ3i/91rV/d6rt7c3bDZbg3yGdru9zseNiN9DniOq9ytKNjxHVO+X5wieI6r2y3NE0zhHnDlzBn/961/x2Wef6VfegPIrUrIsY9y4cXqx2qJFi3v2K9I5oq4M2oOUaQ1gyZIlWLx4MRYvXoxFixY5tFmtVgQFBUGWZXz00UcObWPGjMHly5dx4cIFAMC6deswd+5cHD16FIMGDdJfl5WVhYiICCxbtgxvvPGGvt1ut+uXiSu0aNHC5TfjJScn65fbq04tINdhLuJhJmJiLmJiLuJhJmJiLg8nOzsb69evh6qqOHHihL49JCQEs2fPhizLePzxx+tcvIiQy4PUDkIs8bF06VIsXrwYCxcurFZsAeW/AEycOBFbtmxxeINXrlzBwYMHER8fr28bM2YMfH19sXr1aoc+Vq9eDYPBgClTpjTY+yAiIiIiao5KS0uxadMmTJo0Ce3bt8eCBQtw4sQJeHl5YcqUKdi6dSsyMzOxcuVK9O/f/4GuFLkrl08pXLFiBd5++22MGTMG48ePx9GjRx3aBw8eDKD8CtiAAQMwYcIEvP766/qDj9u0aYNXXnlFf31wcDAWLlyIt956C8HBwfqDjxcvXoyf/exnbrEkPFB+c6XVanV6CZpch7mIh5mIibmIibmIh5mIibnUjqZp+O6776CqKjZs2KAvAgIAAwYMgCzLmDlzZr2t+uiuubh8SuHw4cPx1Vdf3bO98vCOHz+O3/3udzhy5Ag8PT0xYsQIvPfee9VWLAGAv/3tb/jggw+QkZGB8PBw/TlcVedpijqlkIiIiIhIRFeuXEFiYiJUVdVv6wHKV2pMSEiAJEno2bOLine truncated
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAiwAAAEGCAYAAABGqM4kAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/d3fzzAAAACXBIWXMAAAsTAAALEwEAmpwYAAA1k0lEQVR4nO3deXxU5dn/8c+VsAYEWQRkdwFUXAHBugYFBERwAYUiW6b1V5/2ae0jPqit2lq1LlSl1drSngkBApFVKEIRRHApiAJuICiChkV2EAIJkOT+/TETnxATCJgzZ5J8368Xr2TuM3Pmm8swXJ5zn3Obcw4RERGReJYQdAARERGRE1HDIiIiInFPDYuIiIjEPTUsIiIiEvfUsIiIiEjcqxJ0gB+iYcOGrnXr1r7s++DBg9SqVcuXfZcXqkGE6qAagGoAqgGoBgX8rMOKFSt2OefOKDperhuW1q1b88EHH/iy78WLF5OcnOzLvssL1SBCdVANQDUA1QBUgwJ+1sHMvi5uXKeEREREJO6pYREREZG4p4ZFRERE4p4aFhEREYl7alhEREQk7qlhERERkVLZtPsQv8pYxVPLs/lVxio27T4Us/cu15c1i4iISGxs2n2Iwd57ZO6JNClr92xlVeY+0kNdaNEgyff31xEWEREROaHRC9Z916wUyNxziNEL1sXk/dWwiIiIyAlt359T7PiOEsbLmhoWEREROaHGdWoUO96ohPGypoZFRERETmhk93a0rH/sXJWW9ZMY2b1dTN5fDYuIiIicUIsGSaSHutDv0qacVz+Bfpc2jdmEW9BVQiIiIlJKLRokMWbgZdHFDy+L6XvrCIuIiIjEPTUsIiIiEvfUsIiIiEjcU8MiIiIicU8Ni4iIiMQ9NSwiIiIS99SwiIiISNxTwyIiIiJxTw2LiIiIlFpWVhY5ObFZ8LAw3xoWM2thZm+a2WdmttrMfhUdr29mC8zsi+jXeoVe86CZrTezdWZ2o1/ZREREpPSccyxdupSf/OQnnHnmmcyfPz/mGfy8NX8ucJ9zbqWZnQasMLMFwHDgDefcU2b2APAAMMrMLgAGAu2BpsBCM2vrnMvzMaOIiIiUYMeOHYwfP55wOMxnn31GrVq1uPPOO7ngggtinsW3IyzOuW+ccyuj3x8APgOaAf2AtOjT0oBbot/3AzKcc4edcxuB9UBnv/KJiIjI9+Xm5vLaa69x22230axZM+6//35OP/10/vnPf/LNN9/geR5t2rSJeS5zzvn/JmatgbeAC4FM59zphbbtdc7VM7MXgWXOuYnRcQ+Y55ybVmRfdwN3AzRu3LhjRkaGL5mzsrKoXbu2L/suL1SDCNVBNQDVAFQDqNg12LJlC/PmzWP+/Pns2rWLevXq0aNHD3r16kWrVq2Oea6fdejatesK51ynouO+r9ZsZrWB6cC9zrn9ZlbiU4sZ+1435ZwbC4wF6NSpk0tOTi6jpMeKrETpz77LC9UgQnVQDUA1ANUAKl4NDh06xPTp0/E8jyVLlpCQkECvXr0IhUL06dOHqlWrFvu6IOrga8NiZlWJNCvpzrkZ0eHtZnamc+4bMzsT2BEd3wy0KPTy5sBWP/OJiIhUNs45VqxYged5TJo0if3793POOefwxBNPMGzYMJo1axZ0xGL51rBY5FCKB3zmnHuu0KbZwDDgqejXWYXGJ5nZc0Qm3bYBlvuVT0REpDLZvXs3EydOJBwO8/HHH1OzZk369+9PKBTi2muv5ThnQOKCn0dYrgKGAJ+Y2YfRsYeINCpTzCwEZAIDAJxzq81sCrCGyBVGP9cVQiIiIqcuPz+fhQsX4nker776KkeOHKFTp068/PLLDBo0iLp16wYdsdR8a1icc+9Q/LwUgBtKeM0TwBN+ZRIREakMvv76a1JTU0lNTSUzM5P69etzzz33kJKSwsUXXxx0vFPi+6RbERER8V9OTg6vvvoqnufxxhtvANC9e3eeffZZ+vXrR/Xq1QNO+MOoYRERESnHPvroIzzPY+LEiezdu5dWrVrx6KOPMnz48O9djlyeqWEREREpZ/bt28ekSZMIh8OsWLGCatWqcdtttxEKhbj++utJSKh4SwWqYRERESkH8vPzWbJkCZ7nMX36dHJycrjkkkv485//zODBg6lfv37QEX2lhkVERCSObdmyhXHjxhEOh9mwYQN169YlJSWFlJQUOnToEPeXI5cVNSwiIiJx5siRI/zrX//C8zzmz59Pfn4+Xbt25bHHHuO2226jZs2aQUeMOTUsIiIicWLNmjV4nseECRPYuXMnzZo148EHH2TEiBGcc845QccLlBoWERGRAB04cICMjAzC4TDLli2jatWq9O3bl1AoRI8ePUhMTAw6YlxQwyIiIhJjzjneffddPM9jypQpHDp0iAsuuIA//elPDBkyhDPOOCPoiHFHDYuIiEiMbNu2jfHjxxMOh1m3bh21a9dm8ODBpKSk0KVLl0ozgfZUqGERERHxUW5uLnPnzsXzPF577TXy8vK4+uqreeCBBxgwYAC1atUKOmK5oIZFRETEB59//jnhcJi0tDS2bdtG48aNue+++0hJSaFdu3ZBxyt31LCIiIiUkYMHDzJ16lTC4TBvv/02iYmJ3HTTTYRCIXr16kXVqlWDjlhuqWERERH5AZxzLF++HM/zyMjI4MCBA7Rt25annnqKoUOHcuaZZwYdsUJQwyIiInIKdu7cycSJE/E8j9WrV5OUlMQdd9xBSkoKV199tSbQljE1LCIiIqWUl5fHvHnz8DyP2bNnc/ToUbp06cLYsWO58847qVOnTtARKyw1LCIiIiewYcMGUlNT+fvf/87OnTtp2LAhv/jFLwiFQrRv3z7oeJWCGhYREZFiZGdnM2PGDMLhMIsWLSIhIYHLL7+cl19+mZtvvplq1aoFHbFSUcMiIiJSyMqVK/E8j0mTJrFv3z7OOuss/vCHPzB8+HDWr19PcnJy0BErJTUsIiJS6e3Zs4dJkybheR4ffvghNWrU4PbbbyclJYXk5GQSEhIAWL9+fcBJKy81LCIiUinl5+ezaNEiPM9j5syZHD58mA4dOvDSSy8xaNAg6tWrF3REKUQNi4iIVCqZmZmMGzeO1NRUvvrqK+rVq8dPf/pTQqEQl156adDxpARqWEREpMI7fPgws2bNwvM8FixYgHOObt268cc//pFbbrmFGjVqBB1RTkANi4iIVFiffPIJnucxceJEdu/eTYsWLXj44YcZMWIErVu3DjqenAQ1LCIiUqF8++23TJ48mXA4zPvvv0+1atW45ZZbSElJoVu3biQmJgYdUU6BGhYRESn3nHO89dZbeJ7HtGnTyM7O5qKLLuKFF15g8ODBNGzYMOiI8gP51rCYWRjoA+xwzl0YHbsU+BtQA8gF/ss5tzy67UEgBOQBv3TOzfcrm4iIVAxbt24lLS2NcDjM+vXrqVOnDkOHDiUUCtGpU6cyW89n0+5DjF6wjnWZ2czctoqR3dvRokFSmexbSsfPIyzjgBeB8YXGngF+75ybZ2a9o4+TzewCYCDQHmgKLDSzts65PB/ziYhIOXT06FHmzJmD53nMmzeP/Px8rrvuOh555BFuv/12kpLKtpHYtPsQg733yNxzCIC1e7ayKnMf6aEualpiKMGvHTvn3gL2FB0GClaGqgtsjX7fD8hwzh12zm0E1gOd/comIiLlz9q1a7n//vtp3rw5t912GytXrmTUqFF88cUXLF68mCFDhpR5swIwesG675qVApl7IkdcJHbMOeffzs1aA3MKnRI6H5gPGJFm6Urn3Ndm9iKwzDk3Mfo8D5jnnJtWzD7vBu4GaNy4cceMjAxfsmdlZVG7dm1f9l1eqAYRqoNqAKoBBFOD7OxsFi1axLx581i9ejWJiYlceeWV9OrVi86dO8dkAu1Ty7NZuyf/e+Pn1U/ggc41fX//eOTn70LXrl1XOOc6FR2P9aTbe4BfO+emm9kdgAd0I9LAFFVsJ+WcGwuMBejUqZPza02HxYsXV/r1IlSDCNVBNQDVAGJXA+ccS5cuxfM8XnnlFQ4ePMh5553Hs88+y5AhQ2jcuLHvGQqbuW0Va/ds/d54u5ZNSE6+LKZZ4kUQfx9i3bAMA34V/X4q8M/o95uBFoWe15z/O10kIiKVwPbt25kwYQKe57F27Vpq1arFnXfeSSgU4kc/+lGZTaA9WSO7t2NV5r5jTgu1rJ/EyO7tAslTWcW6YdkKXAcsBq4HvoiOzwYmmdlzRCbdtgGWxzibiIjEWG5uLv/+97/xPI85c+aQm5vLlVdeied5DBgwgNNOOy3oiLRokER6qEv0KqFttGvZRFcJBcDPy5onA8lAQzPbDDwK/BQYY2ZVgByic1Gcc6vNbAqwhsjlzj/XFUIiIhXX+vXrCYfDpKWlsXXrVho1asS9995LSkoK559/ftDxvqdFgyTGDLwseiqkcp4GCppvDYtzblAJmzqW8PwngCf8yiMiIsE6dOgQ06ZNIxwOs2TJEhISEujduzcvvvgiffr0oWrVqkFHlDimO92KiIhvnHN88MEHeJ7H5MmT2b9/P+eeey5PPvkkw4YNo2nTpkFHlHJCDYuIiJS5Xbt2kZ6ejud5fPLJJ9SsWZP+/fsTCoW49tprA5tAK+WXGhYRESkTeXl5LFy4EM/zmDVrFkeOHOHyyy/nb3/7GwMHDqRu3bpBR5RyTA2LiEicK1jHZvv+HBrXqRF3V6h89dVXpKamkpqayqZNm6hfvz733HMPoVCIiy66KOh4UkGoYRERiWNF17EB4mIdm5ycHGbOnEk4HOaNN94AoEePHowePZp+/fpRvXr1wLJJxaSGRUQkjh1vHZsxA2N/ee2HH36I53mkp6ezd+9eWrVqxe9+9zuGDx9Oy5YtY55HKg81LCIicWz7/pxix3eUMO6HTzds5b8f/wsrX5/G/i3rqV69OrfeeiuhUIjrr7+ehATf1tEV+Y4aFhGRONa4To1ixxuVMF5W8vPzWbxLine truncated
"text/plain": [
"<Figure size 900x400 with 1 Axes>"
"<Figure size 648x288 with 1 Axes>"
]
},
"metadata": {},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
@@ -252,7 +253,7 @@
},
{
"cell_type": "code",
"execution_count": 4,
"execution_count": 8,
"metadata": {},
"outputs": [
{
@@ -262,7 +263,7 @@
" [ 2.727, 327.235]])"
]
},
"execution_count": 4,
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
}
@@ -288,7 +289,7 @@
},
{
"cell_type": "code",
"execution_count": 5,
"execution_count": 9,
"metadata": {},
"outputs": [
{
@@ -298,7 +299,7 @@
" [ 2.182, 261.788]])"
]
},
"execution_count": 5,
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
}
@@ -320,7 +321,7 @@
},
{
"cell_type": "code",
"execution_count": 6,
"execution_count": 10,
"metadata": {},
"outputs": [
{
@@ -330,7 +331,7 @@
" [6.956, 6.956]])"
]
},
"execution_count": 6,
"execution_count": 10,
"metadata": {},
"output_type": "execute_result"
}
@@ -352,17 +353,19 @@
},
{
"cell_type": "code",
"execution_count": 7,
"execution_count": 11,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAxwAAAFrCAYAAABbv49YAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nOzde1yUZf74/9cwwxlBVPCEHDwECoqtVlhWaq5ZWaZlu9vW5qGjpjb8+ti32tK23U+2+wnQ8tBWa5qWW5ttlrW2fsw20+xkKCioCAGCgsj5MDCH3x985m4GhvMNMwPv5+Pho7yvua6575n3XM57rsOtsVgsFoQQQgghhBCiG3g4+wSEEEIIIYQQvZckHEIIIYQQQohuIwmHEEIIIYQQottIwiGEEEIIIYToNpJwCCGEEEIIIbqNztknoDaz2YzZbLY7ptFo0Gg0TjojIYQQQggheg+LxULTjW49PDzw8HA8ltErE47q6mpnn4YQQgghhBB9hr+/f4sJh0ypEkIIIYQQQnQbSTiEEEIIIYQQ3UYSDiGEEEIIIUS36XVrOBwtDm9tTpkQQgghhBCi/RytmW5tg6Y+kXC0tmpeCCGEEEII0TWtJRzyLVwIIYQQQgjRbXrdCIdwT6dOncJoNKLT6bjsssucfTrCjUksCTVJPAk1STwJtbhbLEnCIVxCbW0tDQ0NeHp6OvtUhJuTWBJqkngSapJ4Empxt1jqkSlVlZWVrFq1ilmzZhESEoJGo2HNmjUdaqOoqIiFCxcyaNAg/Pz8mDJlCv/7v//bPScsepz1bvByR3jRVRJLQk0ST0JNEk9CLe4WSz0ywlFSUsJf//pX4uPjuf3223n99dc7VN9gMHDDDTdQVlbGunXrCA0NZcOGDcyePZt9+/Zx/fXXd9OZi54yYcIEZ5+C6CUkloSaJJ6EmiSehFrcLZZ6JOGIiIigtLQUjUbDxYsXO5xwvPHGG6SlpXHo0CGmTJkCwPTp04mPj2fVqlUcOXKkO05bCCGEEEII0UU9MqWqq0M+H3zwAdHR0UqyAaDT6bjnnnv45ptvOHfunBqnKYQQQgghhFCZWywaT0tL49prr2123DqclJ6ezvDhw1usn56eTkREBIGBgcoxg8FARkYGAMHBwYSHh9vVOX36NDU1NQDEx8fblV28eFFJcsLDwwkODlbKTCYTaWlpAPTr14+RI0fa1c3OzqaiogKA2NhYdLqf34KysjJ++uknAIYNG0ZISIhd3WPHjmGxWPD19W22I0FeXh6XLl0CIDo6Gh8fH6WsqqqKrKwsAEJDQxk6dKhd3RMnTigLj8aNG2dXVlhYSFFREQCjRo0iICBAKaurqyMzMxOAAQMGMGLECLu6p06dora2Fo1G02zor7i4mIKCAqBxBKx///5KmdFoJD09HYDAwECioqLs6p49e5bKykoA4uLi0Gq1SllpaSm5ubkADB8+nEGDBtnVTU1NBcDPz48xY8bYleXm5lJaWgpATEwM3t7eSllFRQXZ2dkADB48mCFDhtjVTU9Px2g04u3tTUxMjF1ZQUEBxcXFAIwePRp/f3+lrKamhtOnTwMwcOBAwsLC7OpmZmZSV1eHVqslLi7OrqyoqIjCwkIAIiMjCQoKUsrq6+s5efIkAEFBQURGRtrVzcrKoqqqCoDx48fb3aempKSE/Px8AMLCwhg4cKBSZjabOX78OAABAQGMGjXKrt2cnBzKy8sBGDt2LF5eXkpZeXk5OTk5AAwdOpTQ0FC7umlpaZhMJnx8fIiOjrYry8/Pp6SkBIAxY8bg5+enlFVXV3PmzBkAQkJCGDZsmF3djIwMDAYDOp2O2NhYu7Lz589z4cIFAKKioqSPQPoI6SMaSR/RSPqIRtJH/Ez6iEY5OTl2n/22uEXCUVJSwoABA5odtx6zdjItMRqNWCwWu2MWi4WGhgal3FEda3lTZrNZKTObzc3K29tu03OybddkMjls12KxONyRwGQydandlq7Vtt2m12r7Gjpq13qtjka3mrZbXFyMyWRCq9XSv39/l3tv2rrWhoYGjEajwxtMtvbetKfdhoYGh2WttWut21a7jrT3NXRU3/acmmrrWuvr6zGbzXYdv6N223oNbWMpJCREuVZHr1F723WVOOzLfURL7Xb3e2MbT9JHNHL3PsLR+fZUH2GNp5qaGukj6B19REvtdncfYY0l2+fsyT7C0fW3psMJx4EDB5g+fXq7Hnv06FEmTpzY0adwqLUpWW1N19LpdM0eo9FolA+c7a8DtnVa2mrMw8NDKXMUGE3btf5j1bTdpudk266jDtXT0xOLxeLwfLVabZfatf1vS+02vVbb19BRu9ZrdfT+NG03Ly9P+XUkODi4x96bltptLV5aeg1tH9PStXamXdv4aW+71rqttdvV17Cta22qrWv18vLCZDJ1+TUsKChQYikkJARPT0/MZnOHPzc92Ue01K70EW23293vjW082fZN0ke4bx/h6Hx7qo+wxlNr8W09p/a0K32E8/uIltrt7j7CGku2z9mTfYSj62+NxuIopWlFYWEhe/bsaddj58+f32xk4uLFi4SEhLB69ep2b407dOhQrr32Wt59912743v27GHOnDns3buXWbNmAY3ZqHWYzKpfv34Og6Yn1NbWEhsby9y5c1mxYkWzYT3RKDU1VflHvenQsxAdIbEk1CTxJNQk8STU4uxY6uj37Q6PcAwdOpT777+/c2fXSePHj1fmhNqyHms6J82V7Nq1i+zsbFJSUli/fj3z5s0jMTGRKVOmuM3eyT0hIiICs9nstMRQ9B4SS0JNEk9CTRJPQi3uFktucZbz5s0jIyPDbvtbo9HI9u3bueqqq5otAHMld999N59++imzZs3CbDbz/vvvc80115CQkMDOnTtbnEPY1/Tv358BAwbYLfoSojMkloSaJJ6EmiSehFrcLZY6PKWqsz799FOqq6uprKxk8eLFLFiwgLvuuguAm2++WdlRYsmSJWzdupWsrCwiIiKAxp0gJk2aREVFBWvXriU0NJSNGzfy0UcfNbvxn6tNqbKVlpZGSkoK27dvx2AwADBixAiWL1/OAw884DZBI4QQQggh+q6Oft/usYQjMjJS2aqtqezsbGW7rYULF7J161a7YwAXLlxg1apVfPzxx9TU1DBx4kSef/55Zs6cadeWKyccVkVFRWzevJkNGzYoW8X5+/uzePFiVq5c2WwbQSGEEEIIIVyFyyYcPcUdEg6ruro63nnnHZKSkpQ9tzUaDXPnzkWv13Pttdf2mXUe1q2LNRpNh3c+EMKWxJJQk8STUJPEk1CLs2NJEg4XTDisAdFa+f/+7/+SlJTEp59+qhyfNGkSiYmJLFiwoMXt23oLZ++2IHoPiSWhJoknoSaJJ6EWZ8dSR79vu97P/r1ElcHImt3pTH1xPwkv7Gfqi/tZszudKkPzG6VoNBpmzpzJJ598wokTJ3jggQfw8fHh+++/57e//S1RUVG8+OKLyt0rhRBCCCGEcBcywtENqgxG5m/8ijNFVZhtXl0PDYwODWDX0msI8G59+Ku4uJhXX32VV155hQsXLgDg5+fHokWLWLlyJWPGjOnOS+hx2dnZGI1GdDqd3KtEdInEklCTxJNQk8STUIuzY0mmVLlAwrFmdzrbDufYJRtWHhq4b0okq2+LBdqebmUwGNi5cydJSUkcO3YMaBwRufXWW0lMTOS6667rM+s8hBBCCCGE80nC4QIJx9QX95NfWtti+bD+PswaN4R9Jy/QYLLgqdUwc+xgHr8xusWRD4vFwueff05SUpLdnd4vv/xyEhMTueuuu/Dy8lL9WoQQQgghhLAlCYeTEw6LxULCC/u5UFHX4mN0HhrMFkunp1tlZGSwbt06tm7dSm1tY2IzbNgwHn30UR566CEGDBigyrUIIYQQQgjRlCwadzKNRoOntvUpTkazpdl0K7MFzhRV8dLezDafIyYmhk2bNpGXl8ef/vQnhg4dSkFBAU899RRhYWE88sgjZGa23Y4QQgghhBDdTUY4ukFrazjaEhbsy8EnZrS5tsNWfX09f//730lOTubo0aPK8VtuuYXExESmT5/u8us8zp49qyx+GjlypLNPR7gxiSWhJoknoSaJJ6EWZ8eSjHC4gMdvjGZ0aAAeTb7ja2icTtWakioDU9e2vZWuLS8vL+69916+//57Pv/8c2677TY0Gg179uzhhhtu4PLLL2fr1q0YDIYuXln3qayspKKiolnwCtFREktCTRJPQk0ST0It7hZLknB0gwBvHbuWXsN9UyIJC/ZlSKAPYcG+LLw6ktBAn1br1jaYyS+r5UJFHfmltWw7nMP8jV+1mXRA43SuadOm8eGHH5KZmcmyZcvw8/MjNTWVhQsXEhkZyR//+EcuXryo0pUKIYQQQgjROplS1QNsp0d1ZrpVR7fStXXp0iVee+011q9fT0FBAQA+Pj787ne/47HHHmPs2LEdupbuYjKZlP/XarVOPBPh7iSWhJoknoSaJJ6EWpwdS7JLlQsmHLZauilgWzqzla6t+vp63nvvPZKTk/n++++V4zfddBN6vZ6ZM2e6/DoPIYQQQgjhfJJwuHjCAY1Jx0t7M/n3yQsYTRa0HlBa00BNvanFOl3dStfKYrFw8OBBkpKS+PDDD7G+/XFxcej1eu6++258fFqf9iWEEEIIIfouSTjcIOGwZZ0e1dbNAlvSdLpVR5w5c4b169fzt7/9jerqagBCQ0NZunQpjzzyCKGhoR1uUwghhBBC9G6yS5WbsU5jmjl2cLNdrdrDbIF/n7yg/L0j+ePo0aNZv349eXl5/PnPfyYsLIyioiLWrFlDeHg4999/P+np6R0/qU4oLS2lpKSE0tLSHnk+0XtJLAk1STwJNUk8CbW4WyzJCIeLaGlthwbQemgwtrLgI7SfNzfHDWFfRlGn1ndYNTQ08P7775OUlMS3336rHJ81axaJiYnMmjWr29Z5pKaLine truncated
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAhgAAAD4CAYAAACuX5+7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/d3fzzAAAACXBIWXMAAAsTAAALEwEAmpwYAAA7pElEQVR4nO3dd3hUVf7H8fchoQekKEHIhF4ElCChuaiJgEZAA1IWZLFSRCkG+CHKurqsjQXpKLKwyK5ilCLdAioKKCVSpAsETGIQJJQQQvf8/iBhQ5iQwEwyk8zn9Tx5yJ05c++5nCfhyznf873GWouIiIiIOxXydAdERESk4FGAISIiIm6nAENERETcTgGGiIiIuJ0CDBEREXE7f0934FpuvvlmW7VqVU93wyudOnWKkiVLerobgsbC22g8vIfGwnvk1lj8+OOPR6y1tzh7z6sDjKpVqxITE+PpbnillStXEhYW5uluCBoLb6Px8B4aC++RW2NhjPklq/e0RCIiIiJupwBDRERE3E4BhoiIiLidSwGGMaaLMWa7MeYPY0zoNdpFGGN2G2P2GmOGu3JNERER8X6uzmBsAx4BvsuqgTHGD5gCPAjUA7obY+q5eF0RERHxYi7tIrHW7gQwxlyrWVNgr7U2Nq1tNBAJ7HDl2iIiIuK98mKbamUgPsNxAtAsq8bGmD5AH4DAwEBWrlyZq53Lr1JSUvR34yU0Ft5F4+E9NBbewxNjkW2AYYxZAVR08tYIa+3CHFzD2fRGls+It9ZOA6YBhIaGWnfv2x0zZgx169albdu2FCqUf3Nctb/ce2gsvIvGw3toLLyHJ8Yi2wDDWtvaxWskAI4Mx0FAoovnvCHnzp1j6tSp7Nu3j9q1a/P888/z2GOPqdKciIiIm+XFf+E3ALWMMdWMMUWAbsCiPLjuVYoUKcLOnTv56KOPuOmmm3j22WdxOBy89NJL/Prrr57okoiISIHk6jbVjsaYBKAFsNQY80Xa65WMMcsArLUXgP7AF8BO4BNr7XbXun3jChcuTLdu3Vi3bh2rV68mPDycUaNGUbVqVXr27MnGjRs91TUREZECw6UAw1r7qbU2yFpb1FobaK19IO31RGtt2wztlllra1tra1hrX3e10+5gjOFPf/oT8+bNY8+ePTz33HMsWLCAxo0bExYWxsKFC7l48aKnuykiIpIv5d8sxxsQn5TKoOhNdJv2A4OiNxGflApA9erVGT9+PAkJCYwZM4b9+/fToUMH6taty+TJk0lJSfFwz0VERPIXnwkw4pNS6TFjHQs3J7I29igLNyfSY8a6y0EGwE033cSQIUPYt28fH3/8MeXLl2fAgAE4HA5eeOEFEhISPHgHIiIi+YfPBBhjlu8m7mjqFa/FHU2l+/S1V81o+Pv707VrV9auXcv3339PmzZtGDNmDNWqVePRRx9lw4YNnrgFERGRfMNnAoxDyWecvp5w7HSWMxoALVq04JNPPmHfvn0MHDiQJUuW0LRpU+6++27mz5+vPA0REREnfCbACCxdLNs2cUdTGbN8t9P3qlatyttvv01CQgJjx44lPj6eTp06Ubt2bSZOnMjJkyfd3WUREZF8y2cCjKFt6hBcrkS27eKOpDhNBE1XunRpoqKi2Lt3L3PmzKFixYoMGjQIh8PBsGHDiI+Pz+LMIiIivsNnAgxH+RJ8+HQzIkMq0aJ6OYLKFnfabvfhU9dMBE3n7+9P586dWbNmDWvXriUiIoKxY8dSrVo1unfvzvr163P7lkRERLyWzwQYcCnImNCtER/1acFHvZpfNaNRoogfqeeuzKnIKhE0o2bNmhEdHU1sbCxRUVEsW7aMZs2aXa6zoTwNERHxNT4VYGSUeUYjMqQSdQJLOW2bXSJouuDgYEaPHk1CQgITJkzgt99+o3PnztSsWZPx48eTnJycm7ckIiLiNXw2wIArZzQmdGtEcPkc5GgcTeXvS7ZdM0+jVKlSDBw4kJ9//pn58+cTFBREVFQUDoeDIUOGcODAgVy6IxEREe/g0wFGZjlNBF21JylHeRp+fn507NiRVatWsX79etq1a8eECROoUaMGXbt25YcffsiN2xAREfE4BRgZ5DQR9OyFP644zkmeRpMmTZg9ezb79+9n6NChLF++nLvuuovmzZvzySefcOHChVy5JxEREU9QgJFJdomgRf2d/5XlNE/D4XAwatQo4uPjmTRpEklJSfz5z3+mRo0avP3225w4ccLt9yQiIpLXFGBcg7NE0Htq3ZLt53KSpxEQEED//v3ZtWsXCxYsoFq1agwdOpSgoCCef/55YmNjc+u2REREcp2/pzvg7dJnNNLFJ6Wy+9DJq55rktmqPUlXLKVsijvOh083w5EpkdTPz4/IyEgiIyPZuHEj48aNY8qUKUyaNIkOHTowePBg7rrrLowx7r0xERGRXKQZjOuUm3kad955J//97385cOAAL7zwAt988w0tW7a8XGfj/PnzuXJPIiIi7uZSgGGM6WKM2W6M+cMYE3qNdgeMMVuNMZuNMTGuXNMb5HaeRuXKlXnjjTeIj4/nnXfe4cSJE3Tv3p3q1aszevRojh8/7u5bEhERcStXZzC2AY8A3+Wgbbi1NsRam2Ugkh/lZp5GyZIl6devHzt37mTx4sXUqlWLYcOGERQUxMSJE9m3b19u3ZaIiIhLXMrBsNbuBHw+PyC38zQKFSpE+/btad++PZs3b2bcuHHMnj2bBQsWEBkZyeDBg2nZsqXPj4OIiHgPY611/STGrASGWmudLn8YY/YDxwALvGetnXaNc/UB+gAEBgY2jo6Odrl/nnA49SLz95zn+FnL76mWpDM5+3suX8xwSwlDmaKGR2oVpkIJP6ftfvnlF1asWMGiRYtITk6mdu3adOnShbCwMPz9lbubl1JSUggICPB0NySNxsN7aCy8R26NRXh4+I9ZrUxkG2AYY1YAFZ28NcJauzCtzUquHWBUstYmGmMqAMuBAdbabJdVQkNDbUxMvk/ZID4plR4z1l0xo1HUv9BViaCZBZcr4XRGA2DlypWEhYWRmprKf//7X8aNG8fu3bupXLky/fv3p0+fPpQrV87t9yJXSx8L8Q4aD++hsfAeuTUWxpgsA4xsczCsta2ttQ2cfC3MaQestYlpfx4GPgWa5vSzBUFu5mmUKFGCvn37smPHDpYuXUrdunV58cUXcTgcPPfcc+zZsye3bktERCRLuT6XbowpCRSy1p5M+/5+YGRuX9fbuDtPI7NChQrRtm1b2rZty5YtWxg/fjzTp0/n3XffpX379gwePJh7771XeRoiIpInXN2m2tEYkwC0AJYaY75Ie72SMWZZWrNAYLUxZguwHlhqrf3clesWBK7W03hr/eks62k0bNiQmTNn8ssvv/DXv/6VH374gfDw8Mt1Ns6dO5cr9yQiIpLOpQDDWvuptTbIWlvUWhtorX0g7fVEa23btO9jrbUN077qW2tfd0fHCwJX6mnsOvrH5Xoa6/YlOV1GqVixIiNHjiQuLo5//etfnD17lscee4yqVavyxhtvkJSUlOv3KCIivkmVPL2EK3kaT87acM3HxxcvXpxevXqxfft2PvvsM26//XZGjBiBw+GgX79+7N69OzdvTUREfJD2M3qRG83TSD138Yrj9GWUoLLFCSxdjKFt6uAoXwJjDBEREURERLBt2zbGjx/PzJkzmTp1Ku3atWPw4MGEh4crT0NERFymGQwvltM8DWcylyXPvIxy063VmT59OnFxcbz66qts2LCBVq1aERISwvvvv8/Zs2dz8c5ERKSg0wyGl8s4q+GsnkaJIn5XzWBklr6MkrHd+tgkGlQuw8mz5wms8zCrNw5g1RcLGDduHE8++STDhw/nueee45lnnuGWW7JfqhEREclIMxj5SMYZjbrlChEZUomZjze5KjnUmcxByMHksyzfeejyLMfTH26hTWQ3fvrpJ7788kvuvPNO/va3vxEcHEyfPn3YsWNHbt2WiIgUQAow8pn0GY3hTYszoVsjmtUof8PLKBmlF/V6/uPNzNgfQK3HXmfFmhh69uzJf/7zH+rXr0/btm1Zvnw57igvLyIiBZuWSAoAdyyjwNVFvdbHFqXBn/oQUacDh9YtYcPKudx///3cfvvtREVF0b17d4oVK+b+GxIRkXxPMxgFjLPtrjldRslc1Ct9GWXz75aD1dtR5/n/8tdRk0g8fpqnnnqKcoGVGfzCXzl8+HBu3Y6IiORTmsEogDJvdwX48OlmjFm+m8PJZyhZ1I+dB1P49fjpy+/n5OFrCckXmFekJiW7j8Pvly2c3LCAcf98nXcmjOEvf/kLUVFR1K9fP1fuSURE8hcFGD7CWY2N9ICjQulipJ69yPKdh7I9T+q5ixhjKF41hOJVQzifFM8fW5fx/n8+YMaMGdx7X2teHDaU+++/X/U0RER8mAIMH3WjRb0yK1zeAWF9ubVZd1I2f84PG5YSERFBucrVqX5fV5q2fpjh7Ro6feS8iIgUXAowBPhf7sa1llGulSzqV7w0N7XoSummHbmwZw2//zCfmP++xca577Cw+cPc2+EvXCxW+orKoiIiUnApwJDLsltGebRJMP8376drznIYv8IUrhvGrXXu5Wz8VpI3LODXbz5k9nfRlKwXRukmHdgUd5zRne5g9oY4DiWfUdAhIlIAKcCQLGWXLBp/7DQJx047/awxhmLBd1As+A7OH/2V5JhFnNq2glNbV3CsSghdNnXEBDfCmEsbma6oLKqAQ0Qk31OAIdflRmpuFC5XmfL396PM3X8hZcvnnPxxMfHRr+BfLojSTSIpWT+cg8lwMPl/Saaa5RARyd8UYMgNy5y3kd0yil/xUtzUvAulm3Tg1K7VnNywgKNfTOH4d/8lIORBSt3ZDv+AckAOnp+igENExKu5FGAYY0YDDwHngH3Ak9ba407aRQATAD9gurX2LVeuK97jRmpulCxeDFM/nJL1wjibsJ3kDQtI/uETktfNo2S9eygLine truncated
"text/plain": [
"<Figure size 900x400 with 1 Axes>"
"<Figure size 648x288 with 1 Axes>"
]
},
"metadata": {},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
},
{
@@ -392,25 +395,27 @@
},
{
"cell_type": "code",
"execution_count": 8,
"execution_count": 12,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAv0AAAFrCAYAAABCETbXAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nOy9d2Ac1bn3/53ZJmlVbatYlm25Se4lNAMGHDA2lwsBk0JCQiDhzZtA8rshgdclBWyS2LIDKb83Ae5NwnVMCOEGCIHE4ILBYFwD2MayrC7LaqtitV1py+zM+8d61rO703d3drQ6n38MmpmzZ545c+Y5z3kKxXEcBwKBQCAQCAQCgZC20KnuAIFAIBAIBAKBQEguROknEAgEAoFAIBDSHKL0EwgEAoFAIBAIaQ5R+gkEAoFAIBAIhDSHKP0EAoFAIBAIBEKaY011B6JhWRYsy0b8jaIoUBSVoh4RCAQCgUAgEAjmgeM4RCfgpGkaNC1tzzel0u/xeFLdDQKBQCAQCAQCYczgdDpllX7i3kMgEAgEAoFAIKQ5ROknEAgEAoFAIBDSHKL0EwgEAoFAIBAIaY7pfPrFAnaVfJQIBAKBQCAQCITxglgMrFLSmzGh9CtFIxMIBAKBQCAQCOMZJaWfaNIEAoFAIBAIBEKaYzpL/3inrq4ODMPAarWioqIi1d1JW4icjYHI2RiInI2ByNkYiJyNgcjZGMwkZ6L0m4zR0VEEAgHYbLZUdyWtIXI2BiJnYyByNgYiZ2MgcjYGImdjMJOciXuPyeCrD5MKxMmFyNkYiJyNgcjZGIicjYHI2RiInI3BTHKmuOgavimGZVkMDw9H/C0nJ4cE8hIIBAKBQCAQCNCnLxNNmkAgEAgEAoFASHOI0k8gEAgEAmHMYjKHBQLBtJBAXgKBQCAQCGMKt4/Bk7trsa/GhUCQg81CYdW8Yjy6phLZDqLaEAhiEJ9+k9HT04NgMAiLxYLCwsJUdydtIXI2BiJnYyByNgYiZ2NQkrPbx+Cupz9AQ7cbrECDoSlgdlE2Xn3oWqL4q4CMZ2NIlpz16MvkrTAZHR0d4dRO5CVMHkTOxkDkbAxEzsZA5GwMSnJ+cndtjMIPACwHNHS78dTuWjz+mQUG9XbsQsazMZhJzuPDfE4gEAgEQpIx2cZ52rKvxhWj8POwHLC3xmVshwiEMQKx9JuM6dOng2XZcePOlCqInI2ByNkYiJyNQUzOxLc88ciNZ47jEAjKL66YIAeO40yRF93MkHnDGMwkZzIjmYz8/PxUd2FcQORsDETOxkDkbAzRcpbyLd95uAWHGnuJb7lOhHKOVt4pioLNIq/MWy3qCiGN94UBmTeMwUxyJrMRgUAgENKaZCl3xLc8OSjtnqyaV4ydh1tEXXxoCrh5XrHutgmEdIZk7yEQCARC2mGEcrdi23609Y9KHi8ryMTB9Tcm5Ld40t06rSYzDwBd2XtI1h9COkGy96QBDMOEJ3WrlTyeZEHkbAxEzsZA5BxJstxuhHK2WCyG+ZaPJ+u02t2TVx+6Fk/trsXeGheYIAerhcLN84rxiIxMyM5MJGTeMAYzyZk8ZZNRXV0dTu20ZMmSVHcnbSFyNgYiZ2Mgco4kWcpdtJwT5Vsux3iLG1CTmefxzyxAtsOKxz+zAI9/ZoHqhZXatscLZN4wBjPJmfjMEAgEAiGtMCql46p5xaAldE0l33K1qFnAGEWyvYG1ZOYRojZoV0/bBEI6kT7mgTQhNzcXDMOkfAso3SFyNgYiZ2Mgcr5EMlM6Rsv50TWVONTYK+kj/siaSs39jybV1mkjXYsSmZnHyLbHKmTeMAYzyTn1PSBEMGPGjFR3YVxA5GwMRM7GQOR8iWQqd9FyznZYdfmWqyXVOelT4VoUT2aeVLY9FiHzhjGYSc5E6ScQCARCWpEs5U5MudbjW66WVFuntcRGJOrek7l7YsTODIFgZojSTyAQCIS0IpHKnRb3lmQo36m0Tiu5Fu054wJ38bxEuf4kc/ck2TszBILZIXn6CQQCgZB2uH1M3MqdGfK6p6oPHMdh+db9cA15Jc+x0hSCLAehEpHofiWzLkG61zwgpDd69GWi9JuMpqamcMDHzJkzU92dtIXI2RiInI2ByFkevcrdpterZa3s911dbkiKx0QsYPSgVHxMinhlQ8ZziGQvSoicjSFZcibFudKA4eHhcD5XQvIgcjYGImdjIHKWR6/ilKjMOfEqb8mMG5BDzrVIjnizCo3n8WxktqTxLGcjMZOcidJPIBAIBEIU8WbOSZbyZqQ7ilxsBEWFXHukSGZWoXRlvBViIxiPoe49Bw8exJYtW3D48GF4vV6UlZXhq1/9Kn784x+Hzxnv7j3BYDD83xaLJYU9SW+InI2ByNkYiJyTg5J7S1lBJg6uvzHm72aIBUgUUq5Fu890oWNA2t9fSjZqGK/j2Wh3svEqZ6NJlpxN7d7z5z//Gffeey++8IUvYOfOncjOzkZjYyM6OjqM6sKYgLx4xkDkbAxEzsaQznJOpbVYb+YcLakuzY6UaxEHqJaN1meYzuNZDqMLsY1XORuNmeRsiNLf3t6O//2//ze++c1v4umnnw7//dOf/rQRP08gEAiEMYSRfs1yqE39Ga3UprqKbrIQ3qOSbL61chY2vV6d8mc4Vkh1ITbC+MCQN+/3v/89PB4P1q9fb8TPEQgEAmEMIKbAmMmvWS6v+7dWzhJdmDyyumJcKG9Ksrn3D0dN8QzHCqkuxEYYHxji03/TTTfhxIkT+POf/4z169fj9OnTmDBhAu666y5s374dubm54XPFfJRaW1sxffr0iPN8Ph/Onj0LACgoKMC0adMirqmvr8fIyAgAYMmSJRHHent70d7eDgCYNm0aCgoKwseCwSBOnz4NIOQbFZ1eqbm5GUNDQwCABQsWwGq9NHENDAzg3LlzAIDS0lIUFhZGXHvq1ClwHIfMzExUVFREHDt//jwuXLgAlmVRUlKCjIyMcL/cbjcaGxsBAEVFRZg8eXLEtWfOnAlHhs+fPz/iWGdnJ7q7uwEAs2bNQnZ2dviY1+tFbW0tAGDChAmYOnVqxLV1dXUYHR0FRVFYvHhxxLGenp6wa9b06dORn58fPsYwDKqrqwEAubm5MSWom5qaws944cKFEVtf/f39aG1tBQBMmTIFkyZNirj25MmTAICsrCzMmTMn4lhrayv6+/sBAHPnzoXD4QgfGxoaQnNzMwCguLgYDocDLMuCpmkUFBSguroaDMPA4XBg7ty5Ee12dHSgp6cHADB79mw4nc7wsZGREdTX1wMAJk6ciLKysohra2tr4fV6YbFYsHDhwohj3d3d6OzsBACUl5cjLy8vfMzv96OmpgYAkJeXh/Ly8ohrGxsb4Xa7AQCLFi2K8OHr6+tDW1sbAKCsrAwTJ04MH2NZFp988gkAIDs7G7NmzYpot6WlBYODgwCAefPmwW63h48NDg6ipaUFADB58mQUFRVFXHv69GkEg0FkZGSgsjJkAe3v7wfLsujv74fH4wEAzJkzB1lZWeHrPB4PGhoaAACFhYUoLS2NaPfs2bPw+XywWq1YsCDSOtrV1QWXywUgVOp8PMwRAFBZWYmMjIzwsc7OTnR1dQEASkpKTD1HuH0Mtu06g92nO8CwHOwWC25ZPCVsAd70ejX+eKgFYh8mmgLu/lQJvlgRkmky54iSkpKIa6urq+Hz+WCz2VBQOBlff7E61soNYGZhFkYCrGp/93SZIziOw9DQEFpaWvC7D/uxq84DVuTeaQCfmZeHX923IuLvbW1t6OvrAxCSv81mC8/Pwjli0qRJmDJlSsS16TJHqPHpX1se1DVHiOkR/PxM0zQ6OztNM0cAY0OPEJsjxPSI/v5+9PX1YXh4GDRNJ0yPMK1Pf3t7O0ZGRvD5z38eGzduxK9+9SscP34cjz/+OE6fPo33339fdvXKMAyi1yYcxyEQCISPi13DH4+GZdnwMZaNnZbUthvdJ2G7wsANYbscx4mmbQoGg+Fr29vbQx+Wi5OImnal7lXYbvS9CmUo1i5/r2LPRm27Znk20ffa2toanuAKCgoQCATAMIzoyyK8V6V2xfobCAREj8m1K7xXrc9crQzFrhf2KRqle/X7/WBZNmLy5eVM03S4L3plKCYjtc/GLOMwUXNEdLtdXV3h9sw6R7h9DLbuqsFLx8+DidBq2AgL8L4al6jCD4RcYw40XMBnZ1yaG8Xule+bVH/5e9Nyr4FAACzLwufzoWrXaTR0e2J99gE09Y6gsjgHXZRXlb97uswRFEWFZXi0bVRU4QdCMjrc6pZtt6urCwzDhOdntzeAZ4/34V8dPnBUFzLstRGuQukyRzy6phLv17nQ3BspP6E7We3pk7rmCLG5R/gdFPZZrl2iR8i/c2J6BC9n/rpE6hFaMUTpZ1kWXq8Xjz/+ODZs2AAAWLlyJex2Ox5++GG8/fbbWLVqlXQnrdaYAUNRVHigCi1pwmukcqLSNB0+JqbkqW03uk/CdsUCN2w2GziOE23XYrHAZrOJDmI17Qr/FWuXb0eIUIZi7fL3Kvayqm3XLM9G6V75+1SSoZ52g8Gg6DG5dvlr5dqNV4ZK9xqN0r3a7XYEg0HJ6+ORIcuysu+NUrtmGYeJmCPE2o0+T6q/euYIfv5VmiOi3VWE7Y4GWHzx6Q9Q54pV9oBLAa5P7j6rwjXm0jfB6DmCVxKOtYsr9Py9DHkZzC7KVowF4NtNpznCarVC4RGC4WJdu4TtCpUit4/Bl3d8jMae0YuLQRbwMBELxXSZI7IdVuy8dzG2/vM0jrV7wVE0HDZrRCE2NXOEmOuY3NzDcVx4x4boEcnRI6LP0dqu1ByhFUPce66++mocOXIEH330EZYtWxb+e11dHSorK7Ft2zasW7cOAEnZ2dvbG95ui96SIiQOImdjIHI2hkTLWW0grdrz5NwWhJQVZAKArjSZRtDb24tgMIjbfncSPW5xKyMLine truncated
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAaEAAAD6CAYAAAAIlvhCAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/d3fzzAAAACXBIWXMAAAsTAAALEwEAmpwYAABfqklEQVR4nO29eZxU1Z33/zm1dlVX7003Wzeb0EJUEGiRoAYXEifDjIlmQcyEJC7Jk2UwER/1l99PH51fHs2AUfM8M8mgJONMVCZuY4YnEjtGYkBAZBWBBgHpBppuurqb6uraq87zR/Utbt065y5Vt7qrqs/79fIlXXXr3nO38z3fnVBKIRAIBALBaGAZ7QEIBAKBYOwihJBAIBAIRg0hhAQCgUAwagghJBAIBIJRQwghgUAgEIwaQggJBAKBYNQwRQgRQn5ICPmIEHKQEPISIaTMjP0KBAKBoLQhueYJEUImAdgKYA6lNEgI+S2A31NK/5X3m/r6ejp16tScjquHoaEhlJeX5/04+aQUzgEojfMQ51AYiHMoDPSew+7du3sppeN439tMGo8NgIsQEgXgBnBWbeOpU6figw8+MOnQfLZs2YKlS5fm/Tj5pBTOASiN8xDnUBiIcygM9J4DIeSU2vc5m+MopWcArAPQAaALwAVK6Vu57lcgEAgEpY8Z5rgaAK8C+CqAAQAvA3iFUvobxXb3ALgHABobGxds3Lgxp+Pqwe/3w+Px5P04+aQUzgEojfMQ51AYiHMoDPSew/XXX7+bUrqQuwGlNKf/AHwZwAbZ318H8M9qv1mwYAEdCd55550ROU4+KYVzoLQ0zkOcQ2EgzqEw0HsOAD6gKvLAjOi4DgBXE0LchBAC4EYAh03Yr0AgEAhKHDN8QjsBvAJgD4APh/e5Ptf9CgQCgaD0MSU6jlL6CIBHzNiXQCAofDq9Aaxra0e3L4TGyjKsWdaCpjr3aA9LUISYFaItEAjGCJ3eAO7YsBMdfYHUZ3s7BvDCnYuEIBIYRpTtEQgEhljX1p4mgACgoy+pGQkERhFCSCAQGKLbF2J+3sP5XCBQQwghgUBgiMZKdmnIBs7nAoEaQggJBAJDrFnWgubadN9Pc60ba5a1jNKIBMWMCEwQCASGaKpz44U7F2FdWzt6fCE0iOg4QQ4IISQQCAzTVOfGMyuuHO1hCEoAIYQEAkFBo8xJWlIZH+0hCUxECCGBQFCwsHKStrsIFi8KCPNfiSACEwQCwYjQ6Q1g9ca9WLF+O1Zv3ItOb0DzN6ycpJ4gFTlJJYTQhAQCQd7JtsqCyEkqfYQmJBAI8k62VRZETlLpIzQhgUCQd7LVaNYsa8HejoE0AdbgIiOekyQKtuYPIYQEggKnFCbAbDUaVk7Sksq+ET1/UbA1vwghJBAUMKUyAbI0Gr1VFpQ5SVu2bEn9eyQEtJopUeRK5Y4QQgJBAVMqE2A+qiyMlIAWwRH5RQghgaCAKaUJ0OwqCyMloEVwRH4R0XECQQEjJkA+IyWgRcHW/CI0IYGggMnFl1JKyH0/CIYw4/LAiAloUbA1vwghJBDkCTOc5qUyAeZyLVi+nzs27MTa264wRUDrGZso2Jo/hBASCPKAmU7zYp8Ac70WPN/Pi7s6chbQpRJ9WMwIISQQ5IFSiWozg1yvhZrvJ1cBrWdspZCnVcgIISQQ5IFSimrLlVyvRT59P7yxdfT6sXrjXnR4A2jvHkQgcrF9hNCUzEVExwkEeUBEtV0k12uRz+g03tjae4bwxr6z2Ns5kCaAAH017wT6EUJIIMgDIxHWm01rhNEg12shBWfcMm8iFk+vxdUTrKZpIqyxuR3WDMGjZCxqtPlCmOMEgjyQ76g2vQ71QuhKasa1kPt+tmzZYtp1ZI2twxvA3s4B1d+NRY02XwghJBDkiXxGtel1qBdKV9J8XAuzAgaUY1u9ca+qEBqLeVr5RAghgaAI0ePsV+tKKp90iyX6Sxpne0cQL5z6AAfPDKDLF059b1bAACtB2O2woqWhHM31noK9PsWKEEICQYEjFxIepxUEFhzr8TO3lZuJ9AiqYsmTUY7zSF93xjZmhcCXSoJwsSCEkEBQwLCEBA+lmUhPVFqx5MmwxsnCrICBYk8QLiaEEBIIChg9k29duQPXzKzPEA56upJqaUuFoinxxqlEErCFIDgF+jBFCBFCqgE8B+AyABTAtyil283Yt0AwltEz+UbjCebnerqSamlLhVL5gTdOOZImWCiCU6APs/KEngGwmVJ6KYC5AA6btF+BYEyjZ/L1hWJ4Y99Z3LFhJzdXiHJ+u7K1GW6HNe0zuVmvUCo/sPJ5JlW7cNPscVg8vRa3zJuYEjJqglNQeOSsCRFCKgFcB+AbAEApjQCI5LpfgaAUyNUsxDKp8TAaot3pDeD+Vw+kJWa6HVasve2K1BhHuvID73rJtbr2jnNoaR7PvZaFIjgF+iCU8tZIOndAyDwA6wEcQlIL2g1gNaV0SLHdPQDuAYDGxsYFGzduzOm4evD7/fB4PHk/Tj4phXMASuM8jJ5DTyCOdbvC6AlefMcaXARrWp1ocFtVfpm5n9eORTEQpiizUoAQtPclEIxlbntprQUPXuUCAPxyfwg7ujKTU6+eYMV35pZpfm/mOSjPo9pJcOtMe9o+9B5L6z7oOS8zxpsLY+l9uP7663dTShfyvjfDJ2QDMB/ADyilOwkhzwB4EMD/J9+IUroeSWGFhQsX0qVLl5pwaHW2bNmCkThOPimFcwBK4zyMnsPqjXvREzyb9llPkGKbrxbPfN6YP+UrjH2/se9sxnYtzeOxdGly3898tA3AQMY2EVsFli5dgl8e3Q6gL+N74qrC0qWL0ekN4PW2dkys88MxFEWDx4mmOndWTv5ObwAPb9iJjr6LwuFs2IkX7rwqtS+910vrPsy4PFMDbK51Y+3X9PuE9Iw3F8bi+8DDDJ/QaQCnKaU7h/9+BUmhJBCMafJpFtKqx9bpDeBI9yD7+P5kgqeaqU0y5b2x7yz2nfbhdH8QR7oHsbK1OTUJG6ldp8dPw7teW4/1GqqPp6w1J/cX6UX4lUaOnDUhSuk5QkgnIaSFUtoO4EYkTXMCwZjGDH+KHh+JMqFSEiBBThHOceV2APzKAKd6/bj9uR043R9M+10gEsfKDTtx+cQK1Fc4cbjLjzMDF7dRi0DTI5B518s7FIH3RF/qGN+/LD0akHeNconeE36lkcOsPKEfAHiBEOIAcALAN03ar0BQkOgJOGBN8kbqjmmFGvMmWq3coub6pB1fLsgkzSkQiWPfaR/3t/EE5X6vFrqtRyDrCcLo6AvgtWNWfAXJ6/PopoP4yzEvwrGLgsmMcGzRimPkMEUIUUr3AeA6ngSCUkJvHoqe8i9qwkxvjo5yHx0qJitlsioABMJxHDx7AZF4bkFK0lhY8LSula3Nqb+V1+tojx9ef2ag7UCYqlaS6OgL4LFNh/DsquynpFwWECJR1hiiYoJAMIzeycNIAidLW5GOc6rXj8PnBhGOXZz83z/hxW+//Wk01bmzrv2mzPuRmFzjwt9fTtN8Oiue3ZFmUssVyd+kpKnOjbW3XYFvPr8rFRIeiMRx/6sH0oS3/Hrxgi+qnURT23v32HnsPO7Fi7s6shIG2daPK6VE2UgkArvdDkJIXo8jmtoJBECaI37HiT7V5M9c/AVKh79cAAFAly+MRzcdBJB97bdAJM5MQH3prqvTQozXtbVrCiCrxdgENBSKcoMIXtzVYahLKSv4wkKAD8/H8c6RHtVxhGMJfPP5XbruJw9JIL50z2I8s+JKXUKk2AMa3nrrLcyZMweEEDidTtx99915P6bQhAQCGDN9KR32Enr8BXpqwW1p78XqjXuxsrVZ0yTEE4hS2wHlKv64bButkkDNtUnt5cVdHWk+IzX6gzHskAURyDUALeG987gX972yH/1DESQATKt1o6WxAk01Lrx/0otoAkhQYCgGIMZIklLAE3h6AxayMasVW0DD+fPn8cgjj+AXv/gF8/tP6hdh9ca9eTUpCiEkEEDf5NETiOPhDTuZQshmIWn+DYA9iempBRdLULyx7yz2dgykhADPJMTTlprrPczJticQx+qNe9HtC3GFKQA4rBZMrnHiua0nMRiOoqnOjQduvjQ1lnKnNSM6Toly0lfT7HYe9+KODTsRS1zUDA+dG8Shc4OwWgBOeTzD6BUG2ZrVCj2ggVKKl19+Gffeey+6uroyvm9pacGPH3sCz35ShY6+AD4G8PHws5gvk6IQQgIB9E0erx2LpiUvyoklKF7c1YFFM+oA8CexWY3lusfU0RfAi7s6VFfuRhzond7AcEWCi34WKwFY8QiReALvHe9P+0w5EUlCtscXwtFuP7xDmUEE8kmfNdbGSid2nPDid/vOcuvbZSOA3A4rU2vTKwyyLdyaa0RkPjh16hQeeughvPTSS8zvH3jgATz00EOoqqoCkPTFdfSl++LyWbRW+IQEAmgnfwLJqCw1tLqadvQFQGDJOI7efbIwkpi5rq09rSQOkBRA4yudqCt3wGlTnw6Uvg25z+SamfXM38gnfeVYr5lRh/O+MLp9Ya4AMsLkGlfqGvx6Vavm/VQjW7OaGYmyuRKLxbB+/Xq43W4QQjB16tQ0AbRo0SJs27YNlFJQSvHEE0+kBBAw8iZFoQkJBNAXDVXtVHfS6+lq6g9H045TPtwpdU9HP1OT0LNy15uYyRvTtPpyvHTPYqxYvz3lz+HBm4hWtjaj7VB3mvZhIUDvYBid3gAz+u2an/4JJlnZ0Fzrzpjsc+mOytOMO/uDWLF+u6qPaDQa4h06dAhr1qzBm2++yfz+8ccfx+rVq+FyuTT3NdImRSGEBIJhtCaPW2facTbsZAYWGOlqygvbZtU7W9nanPLh5JpzojW56GkbwZqIWNW4gWQQwbbjXqx4dgc23n11xrgvBKN6h67K5BoXU9vIRRiwzGo2C8Hp/mDKlzaaodehUAg///nP8cADDzC//9znPodLine truncated
"text/plain": [
"<Figure size 900x400 with 1 Axes>"
"<Figure size 648x288 with 1 Axes>"
]
},
"metadata": {},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"[[10.888 2.922]\n",
" [ 2.922 2.259]]\n"
"[[10.433 2.934]\n",
" [ 2.934 2.317]]\n"
]
}
],
@@ -433,25 +438,27 @@
},
{
"cell_type": "code",
"execution_count": 9,
"execution_count": 13,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAWcAAAFrCAYAAAANekj2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nOy9eXgc1Z33+63q6kVqrTZaLVteZMmyzJIQ8ILDYgx2yIQEZ5gZJmGALO8bm/d5hwTPAIHEYnIJkGvPZG7eQObezAyYAMkkGGImgMAxhnjBBDBeZEuyZcu2JFuStbfU6u7qqvtHu9rV3bV2V1e3pN/nefRgurqqTp+q+tY5v/NbGFEURRAEQRBZBZvpBhAEQRCJkDgTBEFkISTOBEEQWQiJM0EQRBZC4kwQBJGFcJluQDyCIEAQhJjPGIYBwzAZahFBEETqiKKIeOc4lmXBsspj5KwU57GxsUw3gyAIIu14vV5VcSazBkEQRBZC4kwQBJGFkDgTBEFkIVlnc1Za+NOyyxAEQUwGlNbTtBwdJoU4a61oEgRBTFa0xJkUjyAIIgvJupFzJmlrawPP8+A4DrW1tZluTsagfrgE9UUE6ocIdvYDibMMv9+PUCgEp9OZ6aZkFOqHS1BfRKB+iGBnP5BZQ4YUiTjdoxGpHy5BfRGB+iGCnf3AZFuyfUEQMDo6GvNZfn4+LQgSBDGpMattpHgEQRBZCIkzQWigN7HMsoknMYWgBUGCiMMX4LG5qRU7jvUgFBbhdDBYXV+GjWvqkOfmdLcThBWQzVlGX18fwuEwHA4HSkpK0n6+bGU694MvwGPdM3twotcHQfZksAxQU5qHF765FHf/+37V7ds2XDclBXo63xNyUukHs9o29e6iFOju7o66yUznG3A698PmptYE4QUAQQRO9Ppw7398qLl9S1MrNt3eYF+DbWI63xNy7OwHsjkThIwdx3oShFdCEIE2BWGWb3/nWE/6GkdMK2jkLKO6uhqCIEx7t73p2g+iKCIUTm0BkA9Hql1MNX/g6XpPxGNnP5A4yygqKsp0E7KC6doPDMPA6dAWVYZhAA2B5hxTM1Bjut4T8djZD9P7NUgQcayuLwOroq0sA9SW5mluv6W+LH2NI6YVJM4EIWPjmjrUKAiw5I3x3Deu1dz+4Jo6+xpLTGnIlU4Gz/NReyHHTV+Lz3TvB1+Ax5amVrxzrAehsAAny+DmRaX4hy/UR/2cpe18WATnYHBLfRkenMJ+ztP9npBIpR/MahuJs4yDBw9G3WSuvPLKtJ8vW6F+uMSnn34KnudV+2IqLv4pQfdEhFT6gXJrEISF6AnvdBBmIjNM3/mJAgUFBdFE2tMZ6odLUF9EoH6IYGc/kFmDIAjCBsisQRAEMQUgcSYIgshCSJwJgiCykOlt3Y/j5MmTUWP//PnzM92cjEH9cImp1BepuP1NpX5IBTv7IWPi/Mtf/hLf/va34fV64fP5MtWMGEZHR6nCMKgf5Ez2vrCqMMBk7wersLMfMiLOXV1d2LhxIyorKzE8PJyJJhDElEetcMDWfR3Y235hyhYGmCpkxJXuS1/6EhiGwYwZM/C73/0uZuScSVe6cDgc/bfD4Uj7+bIV6odLTOa+aNzejK37OhTzT7MMcM/yuYYLA0zmfrCSVPoh613pfvWrX+G9997DM888Y/epdXE4HNG/6cxk6gcjY4tUxh+Z6gsrxkx6hQPMFAaYTPdEOrGzH2yd0/T29uKBBx7AU089haqqKsP7NTc3o7q6GgUFBdHPAoEAWlpaAADFxcWYM2dOzD7Hjx/H+Pg4ACTEwF+4cAFdXV0AgDlz5qC4uDi6LRwO48iRIwAib7V4o/+pU6cwMjICAGhoaIiJFBoaGsLp06cBAJWVlQllbA4dOgRRFJGTk4Pa2tqYbWfPnsXAwAAAoK6uDh6PJ7rN5/Ohvb0dAFBaWoqKioqYfY8ePRq1gy1evDhm27lz59Db2wsAWLBgAfLy8qLbJiYm0NraCgCYMWMGZs+eHbNvW1sb/H4/GIbBFVdcEbOtr68P3d3dACIJyOV5bnmeR3NzM4BIRNW8efNi9j158iRGRkbAMAyWLFkSc6MPDg7izJkzAIBZs2bhsssui9n34MGDGA8J+M3RcXx0LhBjR/2bxbkIjEWuTdW8GvxsV0fU1upgRHy2jMPXrijEvKoKlJeXxxy3ubkZPM/D7XZj0aJFMdu6u7vR19cHAKipqYHX641uGx8fx/HjxwEAM2fOTLivW1tbMTExAYfDgSVLlsRs6+3txblz5wAAc+fOhcPjvWQf5gWIAo9rZ3mwYeVsLK5dELNve3t7dMZ5+eWXx4y++vv7cfbsWYxPBKDFRCCUsEjY0dERNTXW19fD5XJFtw0PD6OjowMAUFFRgdLS0pjjHTlyBOFwGB6PB3V1sdn5Ojs70d/fDwBYuHAhcnNzo9vGxsZw4sQJAEBJSQkqKytj9m1paUEgEADHcWhoiB3pnz9/Hj09kZfMvHnzsl4jzGCrOG/YsAF1dXVYv369qf2kTFByRFFEKBSKblfaR9oejyAI0W2CICRsN3rc+DbJjyuf/siPK4qi4mJCOBxO6bhqv1V+3PjfKu9DpeNKv1Vphd/oceV9KC1OvXmoEyFegIMFbjvFRrO9xf9WpWszMh7A998dQOdIGPJe2rqvA7uOufCj6wsBAH/9//0ZJy+MxYwczw8HcPD8BP79rpmIleZIH/I8rzjF1Lo2en0oXZv4baIoxhzXNxHCff+RaB9+o20MR/pO4r8fqI6xD2tdc0EQwPM8dOoGwMEm5gaRtykevd8aDAYhCILiqNKKPlSaTRg9brZohBlsE+dXXnkFr7/+Og4cOGDanYfjuIR9GIaJipzSm4njONUVVZZlo9vkD+Pg4GDkpnY4wLKs7nHj2yQ/rtIN6nQ6IYqi4nEdDkdKx5X/V+248cIj78P40asgCNEXidL1Mnpc6beqLU69sP8M9p0aiC5OqV0bid8cG08QZiAyTT8zFMRvjkVGQvHCDAACgK6RMH65/zx+Mi92FCX9zvg+HBwcRCAQiN4TWveh2rWRqjXHe06wEHB1uRNfu6IQ/+f9M8qFYxH5XfGFY51Op+79vbQqB28cT+wHIGLPXFGdl/C5/LrG94PP54v2g9JvdblcCIfDuvdhMn0oCILp5yYbNcIMtiwI+nw+1NTU4Otf/zoee+yx6OcbNmzA9u3b0dnZCafTCa/XSylDs4B09INVi1Mrn96JzkG/6vaq4hwA0P3O7odW6Z4LsK4v1F5OUpJ+X4BH99CEJW02ek4z3hr0bESYcilDL1y4gJ6eHmzZsgXFxcXRv5dffhljY2MoLi7G1772NTuaQmQIKxanjBRgDfECQuHEaagcqQirnWxualUeGYvAiV4fBseCmvsn0+Y8N4dtG67DPcvnoqo4B+UFHlQV5+Ce5XPJjW4SYMvVKS8vx7vvvpvw+VNPPYX33nsPb775ZsLCTyaYNWsWVRiG9f1gRFSNVK02UoDVyem32UwRVqv6Qu/lFNTpn2QLx+a5OWy6vQGbbm9IKUKQno0IdvaDLeLs8Xhw4403Jnz+3HPPweFwKG7LBNnwgsgGrO4HI6JqRHx8AR75GqM9qcCqCGiaUMwUYbWiL4y8nNwcC38w0ZYOWFc4NhX7Jz0bEezsh+n9GiRsQ6+qtZ74SPbTlvOjitvlBVb1irTaXYTVyMupKNeFhWXZ02Yi82RUnJ977rmsyatBpJdUBVOy2aqNP+vK8qN21Gy0teq9nNYsLsu6NhOZhSqhELaRStVqI14aat4M2VCE1aznRDa0mbAWqr6dAuQuFMFoP6QiIGb2FUURy57ciZ4RdVez8gIP9j2yylR7jLTBynsilZeTElrtt1rc6dmIYKcrHc2VCFNYlYLSjHBYtaAIWNf+ZDDiOaEnqlrtBxDdFuQFuDg25rfRaHxyQeIsIzc3lyoMQ70f7ExBGS8kq+vLUvbASKb9Vt0T8b9H/m+jLwyt9u8+0QchLOJk/3jMeZ/b24FXP+mE18MhLCDplxE9GxHs7AcyaxCGsTIFpRIJIsUyWL340qgwVZtt4/ZmPL+3Q3FRkQFw74rU2q/7exSE0YwtWqv/zZBMhCCROmRzJtJGKotyeqiJFAAU5TjR9N3r4XVzmjZbNTFcf+MCPLurHS98cBphDWVzsAzuXlZtiYnDqOiaeeHp9b8Z1F6mZPpIHyTORFpI16KchN6osCiHw+6Hb46KZryIqIkhg8hUng+L0A7qjmB2VKkmZnqie9e1c+B0sLovDOmFZ6T/zSIdO5N2+OkELQgSacHKRbl4RFHUDG8GgCE/H5OZLf48arkrROiHRsuRcl3EZ4GTY0TM9MK1f/PnswgLoqrftoQ8rF2v/83Ch0WMToTw1Wf3UimrLISGozLOnDmDkydPRpO9T1fU+iHVKD85vgCPxu3NWPn0Tiz98R/RPaQ/XX/76HnVbXribgZ5Iqb4vpBG6Fv3daBz0I+ekQl0DvqxdV8H1j2zB74AbyyXiAFhBmJfeFr9nwycg8GWt9s0EzJtaYoUY6BnI4Kd/UDiLGNwcBADAwMYHBzMdFMyilo/WBUWHS9wvaMBQ8IaFiKjbMWk6iZGx0aQRqzxfaGXXW5LU6tlo9z4F976GxegwGNN1Wfp2EazBdKzEcHOfqD5CmEYKSw61UAKNYHTwxcI4fM/eVfRlGD1lF/NRKMnZs/vi3iD3FBbgpc/PJP0aD7+hecL8Lj73/dj2K9cuSOZY3/v1lq8cUR9NgJkJr0qEYHEWcaiRYtotRra/WBFCspkTRC+QBi+wCXzh2QXfeGbSzWz1SVD/kXvD3lfGBmhC2LEt7jQ48DcmV509I8leGuwDANez2tk6RxsXLso+sLTyy1iBJYBKgo9uHVxefRlanQdgZ6NCHb2A4mzDLfbnekmZAVG+yHZxT+rTBCCCBzv8WHNv7yPkQn1EaXLERFLine truncated
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAPEAAAD4CAYAAAA9xADBAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/d3fzzAAAACXBIWXMAAAsTAAALEwEAmpwYAAAtx0lEQVR4nO2dfXQb5Z3vv48kS5YsK7bll7zZBEIw5oaQgPMCoTQU3HZbX9Lb9IXC7sn2loRtt+ek2zrb9vYcWHp3t/eUtFt6tr00JLfl3gIubWjTm95ua1pcGiAhATtpSOKEQLAcO36R7cjyyJIlPfcPeRS9PDOakUbSjPR8zuEQWZrR4/F85/k9v7eHUErB4XCMi6nYA+BwOLnBRczhGBwuYg7H4HARczgGh4uYwzE4lmJ8aX19PV2xYoWqY2ZnZ1FVVZWfAeURPu7CYcQxA8rG/cYbb0xQShuYb1JKC/7fbbfdRtXy0ksvqT5GD/BxFw4jjplSZeMGcJxK6Imb0xyOweEi5nAMDhcxh2NwuIg5HIPDRczhGJyihJg4+sTjFbCnZwCjvjkgMIeVNwtodjuKPSxOBriIC0CiOJpclejqaNWFOBLHVW2rwKlL0xjxBePvP7j/KJ753EZdjJUjDRdxnhkTInhk/1EMTgrxn/UNThddHB6vgAdTxpXK4GRM5E/cv66AI+OoRbM1MSHETAjpI4Qc0uqcpcAL5+fThCKKo5js6RmQFbDImG+uAKPh5IKWM/EuAGcAuDQ8p+GZDrKbLhRbHKMKv7/RVZnnkeh3uWEUNBExIWQ5gI8C+BcAX9binKVCjY0wf14IccjRpOD7W+oc6Opozes4WGa9HpYbRoJQDdrzEEJ+AeBbAKoBdFFKOxmf2QlgJwA0NTXd1t3dreo7/H4/nE5nzmMtNBfH/fjhaRPGAlevc6OdoGu9DY0Oc9HGNSZEsOdYMGlcdZVASzXBXITAaY7gU232vI/xyRNzODISSfv5piVm/N0t6h50Rr1HlIz77rvvfoNS2s56L+eZmBDSCWCMUvoGIWSL1OcopXsB7AWA9vZ2umWL5EeZ9Pb2Qu0xeqC3txcHvrgBe3oGMOabQ6OrEg+sb8GzxwYxOpSb+ZirGXr7RiFpXInHF+p6P3nuNQCTaT8n9kXYsuV2Vecy8j2Sy7i1MKc3A7iPEPIRAJUAXISQn1JK/1qDc5cEzW5H3MOrlfmoxXkSx5UrSh8oqZ9z2tgzfbGXG0YiZxFTSr8O4OsAsDATd3EBS8PyCmcTysn1PFo6k5Q+UFifW1ZjxxKXLSk+XYi1eCnB48QFRsorrNZbnct5lIhuTIhgV3efIpErfaCwPndpOoCOtiZsuM7MNOs5mdFUxJTSXgC9Wp7T6KSKQSvzUcq7rOQ8mUTn8QoLTq/h+PtyprrSB4rU5/zBeTy1nemzSYKHotjwmTiPsMSglfn4wPoW/ObkCMLR5OjCq29P4OZ/+h2qrBZc12BHlJK0Gz6T6Pb0DCR5rQF5U13pAyWXBw8PRUnDRZxHWGLQynx89thgmoABYNwfAgDMzIVxOUGsiTd8JjGpNdW7OlrRNzidJDDWg0np51hk6wMoh9mbiziP5Go+shBvyj+eHVN1XOINn0lMamfMZrcDz3xuIzNclSqix7etwbPHBlU9wDxeAYfPTzDfk/MBlMvszUWcR3IxH1koKVqQQ7zh5UQHxGbM1wZGkqyIhqoKHH1nEjf/0++wyF6B73ziFmxc6Y6/zwpXaSEi8Rze2RDzfblrqVUkQO/wpgB5pKujFY325LTLxBnP4xWwq7sP9+99Dbu6++DxSovT4xXwmX1HshYwkHzDi6L79rZbAAC7D5yIj6HZ7UDXehu2rl2K26+rwx0razEpzOOybw4zc2EMTQXw4P6jOHrBK/t9ciJSilyhRiZTXKtIgN7hM3EeEcXwiq+OaWayZinR3ExcwwGx2t6hqQDze8wEqHFUwDs7LzkWh9WcdsOzxnDo5DD+7ZNr0egwo2tjK/b0DOC3py4jkrL8DkcpvvKLEzj81Q9IfqeUiH5zcgSH357AuuZFeLRzteysLHUOEwEe37ZG9litLSG9wkWcZxodZjzxkXTTTWqW+tufHENg/moucd/gNG5oqpKdgSMUsgIGgGvdVWk3PGsMkSiw62f92HmzBY+8Lm+6+wLJ3ymufwe9Asb9QcwE2GMKRym8/hBePDOOE4OH8cIX7pQUo5QQozTm3FtaY5fsRqLGkWZkBxgXcYHJ5JhKFDAQE/ZMgL0eVMPZ0Rn8uu8S/jAwFr9R35vwMz9LAez9SxgUYdlzuuwV8X9nu14fn53HY4dOYd/2DWnvebwChGAEZGFMrPdTvzO1G0lrUzVmg2EQAqyVmPmN7gDjIs4RNU/wbG/02VB6lY9aIlGKXT/rTxKDwypdoZSpts1iIvjOJ26Jv87UZEBKiADQ77mS9jMl12rMH0xbYohr7q6O1rTjz43OMs9jdAcYF3EOqH2CK+2mkUplhRmhiPysqIRUEQkh6VlODhMBHBUEO//PcditFixZVInBKfnfS+47WNWwma7Vsho76qsqmH6CMd+cKmEa3QHGvdM5oNb7qrSbRiL2ChMCwdwFLEWVRBqoHFEK+IJRXFlIKOnzTMPrz97kv3FxddrPTnnSyxMTaVvixDX17Bpcz1RAcqnAEqbRHWB8Js4BtU9wJd00gNhMJyZjBeajWY1NKQENTPVcuegV4uG1PT0DOOWZxAWv/APPOxMEgQk2iwnBcPI1GpoKYFJFXDmXTDI9wEWcA2qf4KybxUyQFL6xWQiC4dy7rSglNXRUDC5NB/DNQ6cxMDqjeLkxMDaLviGf5PtCKAKH1Qwh4SElJcxMyS96h4tYJam9mtUUM7BulntaG/Ht3w/AF5iHy14BE41icDrIPL6UeXNwSjIri4WgwIK4ts6BK8Fw/NrKxZW1bJBQaLiIVSBV1H5vWwNmgxFFT3BWlw/ROeOby9/aV+/M5GHd/+6kEBe7by6M3QdOSjodeZy4REn9w84G03tIX5oOoH1FbVqcUzx2YDCAX17uS7opPF4BX3vhJF57xwtGIVJZEgpru/Y3m0jabC3lneZx4hKF9Ye1WdjO/FRHVuqxZyeH4zcFAGx78lWMzZSfyZxPzADMJiC08CyISDwdXzo7hl3dyQ9VHicuMcQZ9PDbE2lhk1QvqEiqI0vqpvjY/3wF9gozF3AeaHDZcNmX+br65sI42D+M4xen0L1jE5rdDh4nLiXEGfRg/7Bk3DN1NmY5sqRuCq8/JFnEwMmMiQC3NS/CEpct6ectdQ4sWWRXdS7RIw7wOHFJoSSj6q5VDXDY5LtyKI0Hc9QRpbFc6+99el1aY4E9PQPo80ynHbPYVYkJf5DZBeXNwSkAPE5cUmTKqGqpc+CRzpvS2rCmdoVk3RQcbRicFPCVX5zAcw9tSvo7dHW04s9nhzGZ8iecFkJwWM1Mzz9ZKPXOJk6sJ282F3ECUjOou8qKO1fVp/2h5Lya4k3x0lvD8MlXCXJUMjQVwIeeeBk/3r4+3l2k2e1Akx1pIp4LRyV9GWYTwf17X4uLUKkTS2/ebC7iBFgzqMNqxg8fuDWpFY2IlAPrM/uOYHmtHU2uSnx+rRXPXbDwWVljhFAEn336GH68fT2ePTaIQa+As1Psz1IA1gTPtcioL4jRBWeYGhHqzZvNRZxAs9uBx7etwWefPhaPMQqhCL70sz6sXlaDmeA8qm0VoIjCH4zg/Bg7yX5oKhB3YP2HGfjqh1fif716ESPTAV2kOZYKopCVZG/FbGfpi69GhHrzZnMRp7Dv8LtpN8WIL4gR32hW5wtGgG/+5owWQ+MwUCRgACEFT89BicqnVPTmzeYhpgQ8XgF/OqeuFSxH/5jYW0SnMZ6hxZFIV0crWuqSze5ierP5TJzANw+dVvTETmVRpUXWgcIpLkpTWxudtswfgv6qnkpexGpCAWLcUC2BcFTz3F9O4XlvSkhLyZRCT1VPJS1i1aEACbMrUwsbLmBjUWO3oMpWgUvTydlzXn8IB/uHDVX8AJS4iJWEAhJn6gqJxZPVDNx1QxP8wXk4F7zTx96dwpUyLh00Mu9vbcQD61skPduJzfb0ktAhR0mLOFMoQGn3yWAEcNjMSfsn3fbPPdoNlFNQvDNBPHX4gqxne3DCr6uEDjlKWsSZQgFquk+mxgDD3IQ2LIcveKVWTnHGZ+cl2+HG93DWySxd0iLOlNiupvtkagzQbrVwc9rAyPk4WuoccFdZJdvh6i3tsqTjxGIoQNwYbOvapUkXWmm1ESsGuGQRr1QqNUwEuLetAc98biNaJMTY6KrUZKM4LSnpmRiQDwUoqTa6c6Ub3/r4mrS9dkeu8LrgUmORvSK+zYucFbf7wAnm8Tztskjc0FSFmUAIUwG2afzaO5P4qyf+hFCEYj5CVe+WwDEOU8J80l5OUgkdeku7LFsRK/VMRyiFP8SlWy4kOq+krDi9NRHIWcSEkGYA/xvAYgBRAHsppU/kel4tYXkSs90XiVP6/GlgDDuePo6Z4DzT81yKaZdhAF+hlL5JCKkG8AYhpIdSelqDc+eMlCexzlG2RggnA9OBMHrOXK1aY3me9ZR2mbN3mlI6Qil9c+HfMwDOAFiW63m1QsqT+NbIDPPzCgteOGVEMT3PStB0OiKErACwDsBRxns7AewEgKamJvT29qo6t9/vV30MAAwMsr3I8xLVSnz1y2Hx+1PDuP+JUXx8VWxj9RfOz2M6SFFjI/j4qgo0OtTvLimS7b0topmICSFOAAcAfIlSmrbTFaV0L4C9ANDe3k63bNmi6vy9vb2QOkZc8w56BYz7g6ivqsA19U50dbSi9fIAzk4Oq/xtOJxkAmHgyEgEnoAVofkwxmevpmy+N2vB8w9vyHpNLHdvK0ETERNCKhAT8DOU0he0OKdSWGveoakA+od86BucxuPb1vDOkxzNSK18AmKdX776Qj+e3XFHEUakwZqYEEILine truncated
"text/plain": [
"<Figure size 900x400 with 1 Axes>"
"<Figure size 648x288 with 1 Axes>"
]
},
"metadata": {},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"[[1.004 0.001]\n",
" [0.001 1.003]]\n"
"[[1.005 0.005]\n",
" [0.005 1.005]]\n"
]
}
],
@@ -493,19 +500,21 @@
},
{
"cell_type": "code",
"execution_count": 10,
"execution_count": 14,
"metadata": {
"scrolled": true
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAyMAAAFYCAYAAABAqSPuAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nOy9aYwk6Xnn93/jyvuorMo6uqq6q7v60AyHHAmkzREkUcSQEgag9oOwltcwiMVKFNaGSJnrizZ4DHtEUksTGksCl+R+sUVapBda2BZkytBwRVEHzDG58NKwSWo4Rx/T1d11ZF0ZeURmnP6Q874dV0ZmVVdWZlY9P6DR3RWRkVGZEW+8//d5nv/DPM/zQBAEQRAEQRAEccpI4z4BgiAIgiAIgiDOJyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgCIIgCIIYCyRGCIIgjoDneXAcB57njftUCIIgCGLqUcZ9AgRBEJMOFyCO48B1XXS7XTDGIn8kSYIsy5AkKbKNIAiCIIgoJEYIgiBCeJ4Hz/Pguq6IgnieFxAdfoHhui4AwLbt2OOFBYv/T9zxCIIgCOK8QGKEIAgCCIgPLi44/uiGX5j4t/v/7ofrunBdd+AxwoIlHGkhCIIgiLMC8yjxmSCIcwpPv+IigRM34fc8D5ubm2g0GjBNE+vr61CUR+s5pmnCMAwoioJUKhXYdtRzCv8dFi3+PzwtzB9l6fc7EARBEMSkQZERgiDODYPSr8L7hkVAs9mEaZoAAFmWA/u3Wi3s7OwAABYWFlAulwPH2tzchKIoSKfTKBaLfc/xKFEWID41LBxl6ZceRoKFIAiCGDckRgiCONNwwWHbtkiR8k/Ww7Uf+/v70HUdiqLg4sWLgWMVi0XUajVks9nIRN4vClRVjWxrNBoAeiImn8+LCAYAPHz4EKZpQlEUXLhwIbCNiyZZlo+VGuY4TuBzCL82KcpCqWEEQRDEqCExQhDEmaNf+hWAwEQ/DGMMOzs7cBwHjDG4rhvYf25uDnNzc2KC76dQKEBRFNi2DU3TAtv8QoWfl/+4pmmi2+0Kly4/uq5jZ2cHjDEsLi4Goiqe56HRaEBRFKiqGhFB/Hfy/x3HsFGWQfUsBEEQBHFUSIwQBHEm8NvvJqVfua6LRqMhJvGLi4tiG2MMxWIRBwcHUBQFpmkinU6L7YqiiDSvMOl0OrBveFs2m0W73U78HRRF6Rtx4dGR8LbNzU0AQD6fx/LycmD7/v4+XNeFoigolUp9BQMV4BMEQRDjgsQIQRBTyVHSr8JsbGzA8zwoioKFhYXAvtVqFXNzc0ilUrHHOM5kmjGGpaUluK4r0qH8rK2t9RU5mqYhl8vBsqxIUbw/khFXMF+v12GaJhhjKJVKkW2Hh4dQFAWVSgWZTEZs459tOIo0jGjhotAvovoV4PezOR70HgRBEMTZgcQIQRBTw7DpV57nodPpQNd1qKqKSqUS2C+fz6PRaMBxHJimiVQqJbb7/32SDHLXihMpAFAqlSJCwn/MarUK27YDYoLDBUFcxMU0TXQ6HQAIFNvz192+fRuSJKFUKmF+fj6w3TAMMMagKMqxa1n8UZYwSQX4/P0oykIQBHE2IDFCEMREM2z6lR/btvHGG28A6KVI+cUI0It+VCqVSCH5tBEWWmFWV1f7NmL0i7mwULIsS+wTJxa2trZExOXatWuBbe12W1gc53K5viKMCvAJgiAIgMQIQRAThud56Ha7ME0TjuMEVvzDE0vuUqUoCgqFgvi5qqrIZDIwDAOdTgeWZQWKu3O53On8MmOmXw0L0LMfrlarcBwnIhgYY8hkMrHF+EByxKXVamF/fx8AsLKyEji2ZVnY2toSQsVfjG9ZFm7fvg2g51q2tLQkzsX/dxxUgE8QBDG9kBghCGLshLufv/766+h2u5AkCW9729tiX9PtdvHaa68B6BVv+8UIAOF6VSgUYl2mTptutwvLsuB5HrLZbGxK1mnDJ+NhMplMxNaY43keZmZmYNt27O+QVMdiWZYo4pdlOeIMxmm1WpHjWpYlivG5gOCMqgCfR1ooykIQBDE6SIwQBDEWktKv+ISPTxxbrRZkWQ5ESTRNg6qqsCwLrVYLjuMEJsfhOohxU6/XcXBwAAC4ePFibI3HNMAYw9zcXN/tPP3Ntu2ICPRbIsdFYzhxqWGHh4eBiIs/uuWPkKVSqdhojv89qACfIAhiciAxQhDEqcAFB49++Cec4RXnYrGITCYDx3HwyiuvwHEclMtlrK6uBl4zOzsL27ZRLBanqvYjbrJ9VkilUn1NAPL5PK5evRobVeHd6S3Lin09r2Ph+/oxTRM7OzsAgJmZmUjB/fb2NjzPg6ZpiTU2wGgK8LlQiYu0EARBnHdIjBAEMTLC6Vd+4tydgF7EY2FhAUBvwvfKK68AABqNRmSVulqtjvL0+8IYO7KgyOVyYiI6CWlj44BPyOPSuxhjuHTpEoB4sca73sdFXAZZHOu6Dtd1Ywv+d3d3oeu6sHn2CyF+7cbVxhwlyhKO/vmPkVSAT1EWgiDOAyRGCII4UY7qfmUYBjY2NtDtdjE7O4sLFy6IbZIkoVwuw3EcFIvFyGRumsjlcuemcP5xifuOy+Vy39S7TCaDxcVF2LaNbDYb2Oa3gY4TKqZpwrIsWJYVed9Op4ONjQ0AwOzsbCQ9Tdd1SJIEVVVjozlxNSlxDFuA74/+8doZv9UxQRDENEJihCCIx4ILDr6K3K/+A3hUA+JfGdc0Dd1uF0B89CPcVZwgwqiq2rcXC2MM6+vrsWKDb5ckSRTH+/GLg7g0wK2tLXieB1VVceXKlcA2XdfR6XSgKAqKxWJin5mj2Bx3Oh3cu3cPQC+dkUcRw1EWKsAnCGJaIDFCEMSRGZR+FbfCvLW1hWazifn5+UBOvyzLyOVycF034K5EnC9s20atVgPQi3SclAEBb87YTwxwC2HXdSOCQ5ZlUYwfLor314vEpd01m000Gg0AvVoZP4ZhYGtrC6qqolgsRq57x3EijmH8d/ET54aWFGXhxyCbY4IgJgkSIwRBDEU4/YoTnriEG9MBvUkTn5g1Go1IgfHly5cD++/u7uLg4ACe52F1dXVqnaeI4XFdF7qui/+fthtaXORjUGrdwsICbNuOFTqDLI5N04RpmpG0Ms/zcOvWLWEB7Tdt4Mfl4iauj8yobI7jBAuJFoIgTgISIwRBxDIo/SpMt9vF7u6uEBv+QmFN05BKpeA4DtLpdN/JD8eyLHQ6HQCIRF4mAcYYXNc90mSsXq9jZ2cHnudhcXGRokBTDq9n6sfS0hIsy4Jt2xGhk2Rx3M+hi9NoNNBsNgFEDRy63a4oxs9ms31dzYCTtzmmAnyCII4LiRGCIARHTb/y4ziO6APRaDQirkWXL1+OdSWKg09iJEk6Mza4/LPl/yaCKIqCtbU1APFRimlDVdW+rmkzMzPCmCH8u3qeh1wuF5saBiRHXDqdjrgHq9VqQIx4noeNjQ3Isox0Oo3Z2dnI+x7HMQw4WgF+v/4sJFgI4vxCYoQgzjn+3h/D2O/W63Xouo5qtRpY3c9kMlAURaz6hic3R7GznZ+fF4W5ZwVZlqGqasQVieghSVLiSv5Zg9eyhFEUBSsrK31fVyqVkE6nY2tc/GIgfL+5rgvDMMS/w2Jkc3MTrVYLiqJgdXU1cG62bcOyLKiqKty7wr+L/+9++MeGcDrnoCgLpYYRxNmFxAhBnDOOmn7lhxeiAz23IL8YYYxhbW0NmqbF9pE4CmdxwlEoFFAoFMZ9GsSU0881DOi5a6VSKdi2HaknGdSLxbZtuK4L0zQjIqfVaon7fn5+HjMzM2Kb53nY39+HoijQNC2xvmsY0TJslIUK8Ani7EBihCDOAVxw8AmHP2oRfng7joNGo4FGo4G5ubnA5CKfz4uGf7xJoZ/zUmhOkx1iXDQaDezu7oIxhtnZ2YDATUoNS6VSuH79el+XLUVRoKoqPM9LjLjE1bjs7u4CQGzB/e7uLizLgqIomJ2dHRgVPOkC/G63C1mWoSgK0uk0FeATxARCYoQgzig8/apWq6HZbMJ1XSwvL4tmaf2o1+t48OABgN7kxi8wJEnC6uoq0un0uUqpIUaPP41IURS6vvrgOI5YCDiquQNjrK9Y8TcbDZNOp1EqleA4TuT1gyIurVZLmFGEU8Pq9TpqtRoURcHc3FzAApkvePDx6jipYZ7nodlsitqcmZmZvgX4cREWKsAniNOBxAhBnCHiup+3Wi3U63UAEJMJz/NgGAZ0XUelUgkUyvpXWlutVuQ9ktJETgrTNKHrOlzXpc7l5wTLsnD//n0AvWtscXFxzGc0uXBjh9OaJCfdg4qi4MKFCwML7nkNSHhbnF040BNad+/eFe8frqNpNpvwPE9EPPo1tAwLoDj8UZa4Y/C/wwX4vH6GoiwE8XiQGCGIKYYLDv8DPS6vmsNXUnd3d0UOOF+V5KiqisXFRWQymUgPhNOi2+1ic3MTQC9H/SyIEdM0cXBwAKA3uQo3wiOIYSiXy6fegyUJWZYTa6HW1tZEemgciqLE9mqxLCvwHmH29vZExOXLine truncated
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAPgAAADuCAYAAAAUXsqNAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/d3fzzAAAACXBIWXMAAAsTAAALEwEAmpwYAACZzUlEQVR4nOz9d3gc6XXni3/eqs7IORJgAEgQBBjBzEnSjDTSyppZr9ZykqX1Oq2la/v6Zz93tbp+HPba6yDbq71eB1mSbelakkeSPRrLShM1nGFOCCQBEDlnNIDOXVXv74+3q9EAkUlwOCOc55lnQKCrurq6vu857znf8z1CSsmWbdmWvTNNe6svYMu2bMs2z7YAvmVb9g62LYBv2Za9g20L4Fu2Ze9g2wL4lm3ZO9gcq/x9K8W+ZVu2+SY268RbHnzLtuwdbFsA37ItewfbFsC3bMvewbYF8C3bsnewbQF8y7bsHWxbAN+yLXsH2xbAt2zL3sG2BfAt27J3sG0BfMu27B1sWwDfsi17B9sWwLdsy97BtgXwLduyd7BtAXzLtuwdbFsA37ItewfbFsC3bMvewbYF8C3bsnewrSb4sGWbYKZpEg6HEULgcDhwOBxomoYQm9b3v2U/pCZW0UXfUnS5jyalxDAMTNPEMAyklMn/ZmZmcLlc5OTk4HQ6EUJsAf6Hxzbti97y4A/IUkEdDAbp7OwkLS2N7OxsMjIymJubw+Vy4fP5iEQiaJqGpmk4HI4twG/Zhm3Lg2+ySSkXgHt4eJje3l527dpFJBJhenqaQCCAlJL09HQqKytJS0tDCEHqdyOEWAB4TdtKn7yDbNNW7i2Ab6JJKYnH41iWhWma3Lx5E03T2Lt3bzI0t1/X1dVFNBpNeniv10tOTg45OTn4fL4twL+zbQvgbzczTZN4PA7A7OwsN2/eZPv27ZSWlgIQj8cXAHZgYABN0ygtLUVKSSgUYnp6Gr/fTzAYJC0tLQl4r9e7Bfh3lm0B/O1iqYk0KSV9fX2MjIxQX19PWlpa8nUrAXypcwaDQaanp5meniYSidwFePt1Usrkfl3TNJxOZzJLv2UPrW0B/O1gqSF5PB6nubmZtLQ0du/efRfAFgN8cHAQgLKysjW9TyAQSAI+Go2SkZFBdnY2OTk5eDye5OuklHR3d1NWVkZaWloS8FsJu4fKtrLoD7ulhuRTU1O0traye/duCgoK1nS8EALLstb82oyMDDIyMqioqMCyrCTgW1tbicViZGZmkpOTQ3Z2NuFwGMuyMAwjeY26ri/w8FuAf2falge/R5NSMjc3R09PD9XV1XR2duL3+6mvr0960qVssQcfGhrCsizKy8vv+Zosy2Jubi7p4efm5sjLy6OgoIDs7GxcLlfy2u2Q3i7LbQH+LbGtEP1hNDsUDwaD3L59G8uyyMvLY+fOnasCZDHAh4eHMQyDbdu23ffrbG5upqCggFAohN/vxzRNsrKykh7e6XQCW4B/C20rRH+YLLW2DSok9/v9HDlyhJycnLf46u42IQSZmZkUFxcDajsxMzPD9PQ0fX19SCkXAN7hcCSThfF4fAvwb2PbAvg6LTWRZlkWbW1thEKhJEDWakuVuVaJpu6b6bpObm4uubm5ABiGkQR8T08PQgiys7OT/2matgX4t6ltAXwdtphu2tzcTFlZGbt27aKpqemezr2ZALfD7uXM4XCQl5dHXl4eoADv9/uZmpqiu7sbIUTSu2dlZS0AfH9/P7quU1xcjK7rOJ1OdF3fAvxDYlsAX4Mtrm0PDQ3R19dHfX09GRkZd+2nN2IPEyAcDgf5+fnk5+cDKl/g9/uZmJigs7MTXdeTNXh7m2JHNktl6bcA/9bZFsBXseXopseOHcPhULfvfnnfBxWir9ecTicFBQXJkl8sFmN6eprR0VEmJiZwOBzEYjFycnLIyMhA07Tk/UoFvK7rOByOLcA/QNsC+Aq2Gt3UNjtkvRd7kHvwezWXy0VRURFFRUVJ2qzT6WRoaIi5uTncbneSdJORkZGs8ZumSTQaRQixBfgHZFsAX8JSQ3KAnp4eRkdHOXjwID6f767XbwScbxcwr2ZSStxuN0VFRcksvd0lNzAwQCAQwOPxJPfw6enpW4B/gLYF8EVm17allMRiMVpaWkhLS+PYsWMr8rl/mDx4qi2VwPN4PJSUlFBSUoKUknA4jN/vp6+vj0AggM/nSwLebo1dCvBbajf3blsATzG/3080GiUrK2tddNP1PnymadLV1YXT6VzQDvp2tNUy9EIIfD4fPp9vQaec3++np6dn2U45y7KIRqM0Njayc+dOvF7vlvjFBmwL4MyH5DMzMwQCAcbGxpiZmeHIkSMr0k03YoFAgObmZoqKigDo6uoiFArhcrnQNI1oNIrb7b6v77mZtt6oQwhBWloaaWlplJWVLeiU6+zsJBQKkZ6engS8rW5jAz4SiSzw8FuAX9l+6AGeWtuOx+MMDAywbds2Ghoa7vtDMzQ0RE9PD3V1dfh8viT3XErJ4OAgo6Oj3Lp1C8MwyMzMJDc3dwGV9GG1e7lPQgjS09NJT09n27ZtCzrl2tvbmZmZob29ndzc3AWdcjbg7ZB+qxd+afuhBfjiRNro6ChdXV3k5OSwa9eu+/pepmly+/ZtDMNIltfs+jHMh7FZWVlUVVVhmiazs7NMTU0lqaTZ2dnk5uaSlZWFruv39fruxVYL0ddrizvlLl++TFlZGTMzM8lOuYyMjKSHt6Md0zQX7OG3AK/shxLgi+mm9oOzd+9eJicn7+t7BYNBmpqaKCsrY9u2bcuCITXJlkokgXlm2WKiSW5ubrLu/FbZ/Qb4UpaVlUV2djaVlZULOuVSox07aWd3ym0BXtkPHcBTa9updNNt27YxMzNzXzPZw8PDdHd3s2/fPrKyslZ87WpU0lRmWSwWY2pqakHd2Q5h7ax0qm1mdv5BADz1/JqmkZWVRVZWFtu3b8eyLGZmZvD7/QwMDCzbKWdvxVIB/8OgdvNDA/DFdNPBwUH6+/uTdFNYn+jCSmZHBdFolKNHj655D71WILpcLoqLi5N153A4zNTU1IKstA14W85ps0D4Vpf2NE1LRjs7duxYU6cczAN+eHg4GfK/ExtnfigAnlrbNgyDmzdv4nA4OH78+IL9rJ2tvRcLhUI0NTVRUlLC3r171/yw3Esd3Ov1UlZWtiArPTU1RXt7O9FolHg8zsTEBIWFhckQ9n7awwSI1TrlgCTLLjs7m1AohM/ne8eq3bzjAZ4aks/MzHDz5k127txJSUnJXa+9V4CPjIzQ2dlJXV3dqiH5YrtfD1FqVtqWc7p+/TrhcJiWlhZM01zwgNsebaP2IEL0e7HVOuUikUhy8bcTmO+k1th3LMCXopuOjY1x6NChJemmsHFOuWVZRCIRhoaGOHbs2LpC8tSHZTPCXfvB3LZtGx6PB9M08fv9C3q/7RDXbgVdjz3sAF9si/MZLS0tpKen35XAtFtj7cjq7Qr4dyTAU2vbsViM5uZmMjIyOHr06IoP8Eb24HZILoTg4MGDawbI4ofiQVFVdV1f4NHi8TjT09OMjY1x584dnE5ncv9uN4qsZG83gC9l+fn5yUU/Fovh9/sX3A97AbTvx9sJ8O8ogNs3vrm5mdraWiYnJ2lra1uzuul6Q3T7Idi3bx+3b9++l0t/y8zpdFJYWEhhYSGwsFFkbm4uyRtfjlL7dge4ZVkLFmWXy7XgfkSjUaanp5MVC3tAZE5OTrJxZinAf/nLX+ZHfuRHltwKPkh7xwA8tbbt9/tpb29nbm6OhoaGNVM/1wpwy7Job28nGAxy9OhRXC7XPXvgh6XZZHGjiD1hxabUppJMPB7PQ3HN92Kmaa4Ydbnd7gUVi8ULoN0ptxjw3/rWt3jPe97zoD7GsvaOAHhqIi0cDhMKhXA6nRw5cmRd3mUtAA+HwzQ1NVFYWMiePXuS578fAH/YLJU3blNqbZKJTQ4yTZPp6Wncbvd9p9Q+iMXDsqx1MQMXL4A24FM75VpaWohEIsvmepYyIcQ24ItAMWABn5VSfmbRax4Hvgl0J371z1LK31vpvG9rgC9FN+3s7MTtdrNz5851n281kI6Pj9Pe3k5tbe1dAov3Q/RhMx/o+7GA2OqsmZmZSVbZjRs3knkIm1JrJ6nulVK7OHzeDLuX9xBC4PV68Xq9Czrlzp07R2dnJ+9+97upra3lc5/7XJJrsYIZwP9PSnlNCJEBXBVCvCilvLXodWellB9Y6zW+bQG+mG6ayvW+fPnyhs65HAgsy+LOnTvMzc0lQ/Kljr2XEtvDEqKvx+zEUkVFBV6vd0EJqqur654ptZZlPZDI5n6WKNPS0vjFX/xF/vEf/5GrV6/S2tpKenr6qsdKKYeB4cTPc0KI20AZsBjg67K3JcBTQ/JAIEBLSwvl5eWUl5ff9wciEonQ1NREXl7eiiH/2xGg98NSk2xLUWqnp6cZHh6mra0Nt9udBPxSlNqlzv12pZHaybba2tqNHLsdOARcXOLPJ4UQjcAQ8BtSypsrnettBfDFIfnAwAADAwPU1dWtJQRat9kh+d69e5PMqOVsPQA3DIPbt28Tj8eTrKvNlk3eLFspi56q3QYqfzE9PU1vby+BQCAp9JCbm5uk1KbagwjRN8s2es+FEOnAN4Bfk1LOLvrzNaBSShkQQrwfeB6oXul8bxuAL0U3dTqdHDt2bMm93r2UbyzLoqOjg5mZmTVn4dcK0EAgQFNTE9u2bcPn8yUTVtFoFMuymJiYuC8Ms4fRFu9XbaEHm1JrZ+hzc3NxuVxv+xLcek0I4USB+x+llP+8+O+pgJdSflsI8ZdCiHwp5cRy53zon6LFY4JsuumuXbuSpYvFpuv6urOjtlmWxZUrV8jNzV2X6MNaMvAjIyN0dXVRX1+f5D/bfc+hUIhbt24lpYyEEEnv/la3hK5kGwXhYqEHuw10amoqSalNS0sjHo9jGMbbasHbyD0R6oDPA7ellH+2zGuKgVELine truncated
"text/plain": [
"<Figure size 900x400 with 1 Axes>"
"<Figure size 648x288 with 1 Axes>"
]
},
"metadata": {},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
@@ -539,7 +548,7 @@
},
{
"cell_type": "code",
"execution_count": 11,
"execution_count": 15,
"metadata": {},
"outputs": [],
"source": [
@@ -559,7 +568,7 @@
},
{
"cell_type": "code",
"execution_count": 12,
"execution_count": 16,
"metadata": {},
"outputs": [],
"source": [
@@ -575,7 +584,7 @@
},
{
"cell_type": "code",
"execution_count": 13,
"execution_count": 17,
"metadata": {},
"outputs": [],
"source": [
@@ -591,7 +600,7 @@
},
{
"cell_type": "code",
"execution_count": 14,
"execution_count": 18,
"metadata": {},
"outputs": [],
"source": [
@@ -608,7 +617,7 @@
},
{
"cell_type": "code",
"execution_count": 15,
"execution_count": 19,
"metadata": {},
"outputs": [
{
@@ -617,7 +626,7 @@
"0.0315"
]
},
"execution_count": 15,
"execution_count": 19,
"metadata": {},
"output_type": "execute_result"
}
@@ -636,7 +645,7 @@
},
{
"cell_type": "code",
"execution_count": 16,
"execution_count": 20,
"metadata": {},
"outputs": [
{
@@ -649,10 +658,8 @@
],
"source": [
"import scipy\n",
"try:\n",
" print('{:.4f}'.format(scipy.stats.multivariate_normal(mu, P).pdf(x)))\n",
"except:\n",
" print('you have an old version of scipy, upgrade it now!')"
"from scipy.stats import multivariate_normal\n",
"print(f'{multivariate_normal(mu, P).pdf(x):.4f}')\n"
]
},
{
@@ -1891,7 +1898,7 @@
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.7.6"
"version": "3.7.4"
},
"widgets": {
"application/vnd.jupyter.widget-state+json": {