Skip to content

Instantly share code, notes, and snippets.

@fomightez
Created March 15, 2024 21:11
Show Gist options
  • Save fomightez/e92d4e1abf2657676cda84560d5d8f2b to your computer and use it in GitHub Desktop.
Save fomightez/e92d4e1abf2657676cda84560d5d8f2b to your computer and use it in GitHub Desktop.
Adding a trend line after the fact in a different cell of a notebook.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "markdown",
"id": "94dcaf35-a686-48ff-a624-2c7d3048217e",
"metadata": {},
"source": [
"# SO Plot with trend line\n",
"based on https://stackoverflow.com/q/78168545/8508004"
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "ad86baa1-e5b8-420b-80e2-564496ba42a2",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAh8AAAGdCAYAAACyzRGfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAwOElEQVR4nO3df3RV9Znv8c9JlBCZ5AxkFpxEY025TjVEBYaq/FiOo6AwFLF3Zjq2MItx7sIOxlFkpiN2ipaiRpxVh2u1IK52nFkM2v5RpNTKXRQsFicKSuOURf2FsWWQYAt6DsZL9Obs+0fcKUnOr33O/vHde79fa+WPJCfZX6PL73Oe7/M834RlWZYAAAB8UhX0AgAAQLwQfAAAAF8RfAAAAF8RfAAAAF8RfAAAAF8RfAAAAF8RfAAAAF8RfAAAAF+dEfQChstms3rnnXdUV1enRCIR9HIAAEAJLMvSyZMn1dTUpKqqwrkN44KPd955R83NzUEvAwAAlOHw4cM655xzCr7GuOCjrq5O0sDi6+vrA14NAAAoRSaTUXNz8+A+XohxwYd91FJfX0/wAQBAyJRSMkHBKQAA8BXBBwAA8BXBBwAA8BXBBwAA8BXBBwAA8BXBBwAA8BXBBwAA8BXBBwAA8JVxQ8YAAIiL/qylvd0n9O7JUxpfN1qXtoxTdVX07zUj+AAAIADbDxzV6m0HdTR9avBrjcnRuntBq+a2NQa4Mu9x7AIAgM+2HziqZZv2Dwk8JKknfUrLNu3X9gNHA1qZPwg+AADwUX/W0uptB2Xl+J79tdXbDqo/m+sV0UDwAQCAj/Z2nxiR8TidJelo+pT2dp/wb1E+I/gAAMBH757MH3iU87owIvgAAMBH4+tGu/q6MCL4AADAR5e2jFNjcrTyNdQmNND1cmnLOD+X5SuCDwAAfFRdldDdC1olaUQAYn9+94LWSM/7IPgAAMBnc9satX7xVKWSQ49WUsnRWr94auTnfDBkDACAAMxta9Sc1hQTTgEAgH+qqxKaPrEh6GX4juADAIAc4nrvih8IPgAAGCbO9674gYJTAABOE/d7V/xA8AEAwCe4d8UfBB8AAHyCe1f8QfABAMAnuHfFHxScAgDwCbfuXaFTpjCCDwAAPmHfu9KTPpWz7iOhgSmkhe5doVOmOI5dAAD4hNN7V/qzljoPHdfWriPqPHRcP/4vOmVKQeYDAIDT2PeuDM9epIZlL3JlOKoSytspk9BAp8yc1lTsj2AcZz6ee+45LViwQE1NTUokEnrqqaeGfN+yLN11111qbGxUbW2tZs+erTfeeMOt9QIA4Lm5bY3ac8dVemLp5frfN0zWE0sv1547rhoSeOTKcBTqwKVT5nccBx+9vb265JJL9Mgjj+T8/gMPPKCHHnpIGzZs0IsvvqgxY8bo2muv1alTVAYDAMLDvndl4eSzNX1iw5CjlnyzQEpBp0wZxy7z5s3TvHnzcn7PsiytW7dOX/va17Rw4UJJ0r//+79rwoQJeuqpp3TDDTdUtloAAAJWbBZIMaV21ESZqwWn3d3d6unp0ezZswe/lkwmddlll6mzs9PNRwEAEIhyMxcJDXS9FOqUiQtXC057enokSRMmTBjy9QkTJgx+b7i+vj719fUNfp7JZNxcEgAArionc5GrUybOAm+17ejoUDKZHPxobm4OekkAgBgZ3i5b7N4WexZIoRBieHyRSo7W+sVTmfPxCVczH6lUSpJ07NgxNTb+7g987NgxTZ48OefP3HnnnVqxYsXg55lMhgAEAOCLcgaC2bNAlm3ar4SGttbaMcfDX5yisWNqmHCah6uZj5aWFqVSKe3cuXPwa5lMRi+++KKmT5+e82dqampUX18/5AMAAK/la5ctZSCYPQsklRx6BGNnOP704qacnTIY4Djz8cEHH+jNN98c/Ly7u1tdXV0aN26czj33XC1fvlz33HOPzj//fLW0tGjVqlVqamrS9ddf7+a6AQAoW6F22VIHgs1ta9Sc1hR3uJTBcfDx0ksv6U/+5E8GP7ePTJYsWaLHH39c//iP/6je3l7ddNNNev/99zVr1ixt375do0fTWgQAMEOxdtnTB4JNn9iQ93X2LBA44zj4uPLKK2VZ+YtxEomEvvGNb+gb3/hGRQsDAMArPzmYuwNzOAaCeSPwbhcAAPzUn7W0petISa9lIJg3CD4AALGyt/uETvR+XPR1DWNGMRDMIwQfAIBYKfUoZeHkJopHPULwAQCIlVKPUua0pjxeSXwRfAAAYqXUCaXv9X7k25rihuADABAr9oTSQrKW1L658KAxlI/gAwAQKk7vYsllblujHvnSlBF3sAy3etvBsn4/CnP1bhcAALxUzl0s+YwdU6NCcUWpg8bgHJkPAEDZ3MhClKqSu1hyKbXrhUFj7iPzAQAoi5tZiGLcuItluFK7Xhg05j4yHwAAx9zOQhTj5C6WUhXrekloIJhi0Jj7CD4AAI4Uy0JI7hdqenFEcnrXy/AAxP787gWtkRo05ucxWSEcuwAAHHHrRlgnvDoimdvWqPWLp444Pkp5dHwUJD+PyYoh+AAAOBJEoaZ9RNKTPpUz45LQQMBQzhHJ3LZGzWlNaW/3Cb178pTG1w38nihlPOxjsuF/O/uYbP3iqb4GIBy7AAAcCaJQ0+sjkuqqhKZPbNDCyWdr+sSGSAUeQRyTFUPwAQBwJKhCTfuIJJUcGtSkkqN9f+deClPqK7wo1q0Uxy4AAEfsLMSyTfuVkIa8o/a6UDMsRyQm1VeYOM+EzAcAwLEgsxDFjkiCzjj43YZcjInzTMh8AADKYmIWIuiMgxfD0CrlZbFuuch8AACGcJI5MKlQ04SMg4n1FSbOMyHzAQAYFHTmoFymZBxMrK+QzJtnQvABAJBk3iwIJ4IYfJaLifUVNpOOyQg+AADGZA7KZUrGwcT6itPZx2RBo+YDAGBkrYITpmQcTKyvMBHBBwDAmMxBuUy6oTZsw9CCwLELACDwzEF/1qqoFiHIwWe5mFRfYSKCDwBAoLUKbnXYmNbRYUp9hYkSlmUFM2w+j0wmo2QyqXQ6rfr6+qCXAwCxYXe7SLkzB14cGeTrsKnkmZVmUVAeJ/s3NR8AAEn+1yp4dduqSYPPkBvHLgCAQX7WKpgymwP+I/gAAAzhV61C2DtsUD6OXQAAgQi6wwbBIfgAAATCpNkc8BfBBwAgEEwDjS+CDwBAYJgGGk8UnAIAAsU00Pgh+AAABI5poPHCsQsAAPAVwQcAAPAVxy4AAHiM+2aGIvgAAHjKr43X1A3erVt7o4TgAwDgilyb/46DPb5svKZu8Plu7e1Jn9KyTftj206csCzL2XWBHnNyJS8AwAy5Nv/fP+tMvf/hxyNea+ci3Np4823wbj/Hqf6spVlrd+W9PC+hgXkme+64yogMTaWc7N8UnAIAKmJv/sM32VyBh6TBIGH1toPqz1b2/rc/a2n1toMjAg+3n1MOJ7f2xg3BBwCgbP1ZS1//Ye7NvxC3Nl6TN3hu7c2P4AMAULaHd72hnkz5m2elG6/JGzy39uZH8AEAKMv2A0f1Lz95o6LfUenGa/IGz629+RF8AAAcs2styuXWxmvyBs+tvfkRfAAAHCtWa1GImxuv6Rs8t/bmxpwPAIBjTmoohrfcplyev2Fv8MNbfd1+Trm4tXckgg8A8JGpUzidKrWG4vdqztB917dp7JgaT/+ZTd/gubV3KIIPAPCJqVM4y2HXWvSkTxVss+3t+39q3/xzrV88VQsnn+3pmtjgw4OaDwDwQb5BXPaY7e0Hjga0svKcXmtRSFCDvvqzljoPHdfWriPqPHQ8kCFjyI/MBwB4rNgUzoQGNuc5rSljjglKYddafHXLL3SiN/c0U2nooC8/MhNRyjBFFZkPAPCYyVM4KzW3rVGrPjeppNf6MegrahmmqHI9+Ojv79eqVavU0tKi2tpaTZw4UWvWrJFh99cBgG9MnsLphlS9GYO+TL7nBUO5fuyydu1arV+/Xv/2b/+mSZMm6aWXXtKNN96oZDKpW2+91e3HAYDxTJ7C6YZixaf27a1eD/pykmGiMDVYrmc+/vM//1MLFy7U/Pnzdd555+nP//zPdc0112jv3r1uPwoAQqGSKZxhKJw0ZdBX1DNMUeJ68DFjxgzt3LlTr7/+uiTplVde0Z49ezRv3rycr+/r61MmkxnyAQBRUu7mvP3AUc1au0tffOwF3fZkl7742AuatXaXkXULJkzyjHqGKUpcP3ZZuXKlMpmMLrjgAlVXV6u/v1/33nuvFi1alPP1HR0dWr16tdvLAACjOJnC2Z+19PCuN3Je2mYXTpo4mntOa0p1NWeq863fShqYuXH5pxt86+Ax5fgHxSUslytBn3zySX3lK1/RP//zP2vSpEnq6urS8uXL9eCDD2rJkiUjXt/X16e+vr7BzzOZjJqbm5VOp1VfX+/m0gAgcMUmnG4/cFRf/+HBgtfU25vonjuu8nxjL3UiayXtrW5OfbW7XSQNCUDs32Zi0BYVmUxGyWSypP3b9eCjublZK1euVHt7++DX7rnnHm3atEmvvvpq0Z93sngACAMnG/iyTfsLTgw93RNLL/e0cLLUgCLfukvZ8L2YycGcj2A42b9dP3b58MMPVVU1tJSkurpa2WzW7UcBgPFK3QgLtYnm42XhZL6AYvixTyUD1Ep9hlOm3/MCDwpOFyxYoHvvvVdPP/203n77bW3ZskUPPvigPv/5z7v9KAAwmpOBV+VcUT++brQn3TBO5mWUO0DN65kc9j0vCyefrekT/as7QWlcz3x861vf0qpVq3TzzTfr3XffVVNTk7785S/rrrvucvtRAOAbp3UJTjMCTrIYds3He719mrV2l+vHC04CinLbW5nJEW+uBx91dXVat26d1q1b5/avBoBAlFND4HRzddr+ed0ljWrf/POyjywKBVNOAopy21uZyRFvXCwHAAWUW5fgdHMt9Yr6VH2N7vrcJK15uvyL6ooFU04CinLbW5nJEW9cLAcAeVRSl+B0cy00iMx2++w/1PMrr9bYMaPKvqiulDoUJxNZyx2gVsnUV4QfwQcA5FHJbbTlbK75poQ2Jkdrw+Kpum32+Y7qQ4a/rtRgSpKjgKKc6aamjGRHMDh2AYA8KqlLsDfXZZv2K6HcA69yba6ltImWe2ThJJhyMpG11HUP5/QZiA6CDwDIo9K6hHI3V7tNNJ9y6yycBlNOA4pi686FmRzxRPABAHm4cVeIF5truVmVcoKpcgIKp/x4BsxCzQcA5OFWXYIXA6/KqbOgyBOmcP1ul0pxtwuAIOWaf7HjYI+xd4U4HX7GxWvwSqAXy1WK4ANAUArNv4hSXQIXr8ELBB8A4FAlN7OGkZvX2ANSwLfaAkDYVHIza1i5VeRJEINyEHwAiD0uOSsPxzcoF90uAGKPS86cK2VMO5APwQeA2OOSM2cqufMGkAg+AID5Fw5VcucNIBF8AACXnDnEMRUqRfABACpvYqhJ+rOWOg8d19auI+o8dNzTIw+OqVApul0A4BNhveTM764TN+68QbyR+QCA03hxD4uXgug64ZgKlSL4AICQCrLrJOzHVAgWxy4AIicsUzcrXWfQw9HCekyF4BF8AIiUsEzddGOdJnSduDWmHfHCsQuAyAjL1E231knXCcKK4ANAJIRl6qab62Q4GsKK4ANAJIRl6qab66TrBGFF8AEgEkyof3Dz+aW+jq4ThBEFpwAiISz1D16sk64ThA3BBwCjlNt+Gpapm16tk64ThAnBBwBjVNJ+atc/LNu0XwlpyMZuUv1DWNYJeImaDwBGcKP9NCz1D2FZJ+CVhGVZwfadDZPJZJRMJpVOp1VfXx/0cgD4oD9radbaXXm7QOyjiD13XFVSRiAuE04BkzjZvzl2ARA4t8eEh6X+ISzrBNxG8AHAdU7f0YelTRaAOwg+ALjqx/91VF/bekAnej8a/FqxotGwtMkCcAcFpwBc0/Hjg7p58/4hgYc0cGRSqGiUMeHl689a6jx0XFu7jqjz0PHAx8cDpSDzAcAVP/6vd/Toc915v29p4M6SOa2pEUcwtJ+WJyw3+ALDkfkAULH+rKWvbT1Q9HWF7iyh/dSZsNzgmw8Zm3gj8wGgYnu7T+hE78clvbZQ0ShjwktT7GbchPJnmUxAxgZkPgBUzEkXSrGiUbv9dOHkszV9YoORm2fQwnKDby5hz9jAHQQfACpWahdKw5hRFI26IKytycUyNtJAxoYjmOgj+ABQMbtbpZg1C9vIZLggrK3JYc7YwF0EHwAqZnerFAorvnxFi/70Ys7z3RDW1uSwZmzgPoIPAK6wu1WGZ0DGjTlT3/7SFN35p60BrSx67GBP0ogAxOTW5LBmbOA+ul0AuIZuFf/Ywd7wrpGUwV0jdsamJ30qZ92HfYGgaRkbuI9bbQEgxMJ2M67d7SLlHibHTJfwcrJ/E3wAAHzFnI9ocrJ/c+wCACEQtgxHIRzPgeADAAwXxUyBPUwO8US3CwAYjImgiCKCDwAwFBNBEVUEHwBgKCaCIqoIPgDAUEwERVQRfACAoZgIiqgi+AAAQ4X1DhegGE+CjyNHjmjx4sVqaGhQbW2tLrroIr300ktePAoAIiusd7gAxbgefLz33nuaOXOmzjzzTD3zzDM6ePCgvvnNb2rs2LFuPwpABPRnLXUeOq6tXUfUeeg4nRvD2He4pIZd2JdKjmYUOULL9fHqK1eu1PPPP6+f/exnZf0849WB+Iji8CyvRGnCKaIp0LtdWltbde211+q///u/tXv3bp199tm6+eabtXTp0pJ+nuADiAd7eNbw/wHF+YIxAgyEWaB3u7z11ltav369VqxYoa9+9avat2+fbr31Vo0aNUpLliwZ8fq+vj719fUNWTyAaCs2PCuhgeFZc1pTsdl8yQIhTlyv+chms5o6daruu+8+TZkyRTfddJOWLl2qDRs25Hx9R0eHksnk4Edzc7PbSwJgGIZnDcUIdcSN68FHY2OjWltbh3ztwgsv1K9//eucr7/zzjuVTqcHPw4fPuz2kgAYhuFZv8MIdcSR68cuM2fO1GuvvTbka6+//ro+9alP5Xx9TU2Nampq3F4GAIM5GZ4V9ToIJ1kgboFFVLgefNx+++2aMWOG7rvvPn3hC1/Q3r17tXHjRm3cuNHtRwEIKXt4Vk/6VM53/AkNtJK+19unWWt3RboOgiwQ4sj1Y5fPfvaz2rJli5544gm1tbVpzZo1WrdunRYtWuT2owCEVCnDs667pFHtm38e+ToIRqgjjlxvta0UrbYIo6gfDXglX4fHqvmtWvP0wbzHEXZmZM8dV4X+79yftTRr7a6iWaAo/LMi2gJttQXihhbJ8s1ta9Sc1tSIwC1OdRB2FmjZpv1KSEMCEEaoI6q4WA6oAC2SlauuSmj6xAYtnHy2pk9sUHVVInZ1EIxQR9yQ+QDKFMVBWaYcH8WxDiJfFigs/+0AThB8AGWK2tGAScdHpXbDRO0qeTsLBEQdxy5AmaJ0NGDa8ZGXV8lziy4QPDIfQJmicjRg6vGRXQcxPBuTqiAbY1J2B4gzgg+gTFE5GjD5+MjNOoh8t+ja2R0KOwH/cOwClMnLowE/mX58lKsbxinuTwHMQvABVCAKLZJROT4qhFt0AbNw7AJUKOwtklE5PirE9OwOEDcEH0ABpc69CHOLZBwmbMYhuwOECcEHkIcfnRGmDPXyorPEJHHI7gBhwsVyQA75OiPssMCNeg4T2z5LDYZMCZqcsP+dSrmzO2Gp0QFM5WT/JvgAhrFvGfXyRlWvgxsvg4NKg6YgAxcTAz4gKrjVFqiA13MvvB7q5eUGW+msjKA3/7AXBwNRQastMIzXnRFetn16OSa90lkZpoxwd2NuCIDKEHwAw3jdGeFVcOP1IK1KgiaGfAE4HcEHMIzdGZHv/XBCA0cF5XZGeBXceD1Iq5KgiSFfAE5H8AEM4/XYdK+CG6+PiyoJmhjyBeB0BB9ADl6OTfcquPH6uKiSoIkhXwBOR7cLkIeXnRGVDvXK1a7q9SCtSiahMuQLwOmY8wEEqJyZF4XaVSV5Pkir3HZZhnwB0caQMSCiShlOJsnYsfBBz/kA4B2CDyCCnExelWTsIK1cgYtk7noBlIYJp0AEOZ28auotu8NvACYbAsQP3S5ASESxXdWUqacA/EXwAYRE1NpVmXoKxBfBBxASXk9e9RtTT4H4IvgAQsLryat+i+IxEoDSEHwAIeLl5FW/Re0YCUDp6HYBQsbLyat+YuopEF8EH0AIDW9XDaNyx7WXO+AMgDkIPgAExukdN8wEAaKBCacIBd7tRlsp/35LGS1PAAIEhwmniBTe7UZfsWOkYjNBEhqYCTKnNUVQCoQA3S4wGhMwITETBIgagg8YiwmYsDETBIgWgg8Yy4t3u/1ZS52Hjmtr1xF1HjpO4BISzAQBooWaDxjL7Xe71I6EFzNBgGgh8wFjuflul9oRc5WSjYraaHkg7sh8wFhuvdulU8JcTrJRTmeCADAXwQeMZM99mNeW0neff9vRBMzhnNSOhH1qaJjkm9thZ6Nyze2Iymh5IO4IPmCcXO+GEwnp9HF4Tt7t0ilhnkqyUdVVCV3aMm4wANnbfYIABAgZgg8YJd+7YbsM4H/NPE+zW1OONhs6JczSn7X0+PPdZWejKBwGwo+CUxij0LthaeDd8I8P9Dh+l2vXjuT7iYQGNq+odUqY2Fa8/cBRzVq7S2ue/mVJrx+ejaJwGIgGMh8whle1GeXenhpmJmYH8mW1Cjk9G0XhMBAdZD5gDC9rM+xOiVRy6NFKKjk6cheSmZgdKJbVGi5XNooR60B0kPmAMbyuzYhDp4Sp2YFigcPp8mWjKBwGooPgA8bwY4plsdtTw87UtmInAUG+TiYKh4Ho4NgFxmCKZeVMzQ6UGhCsmn+h9txxVc5jsLgWDgNRRPABo7hdm2Fix4eXTM0OlBo4/PXMlrzBJcEpEB0cu8A4btVmmNjx4TVTL2Bzq+OIEetANCQsyzLqrWAmk1EymVQ6nVZ9fX3Qy0FI5WvrtLc2Ezpc7BHybhe/2v/sknIGIOUManOLWwGhV387AOVzsn8TfCBy+rOWZq3dlbfw0n73v+eOq1zfsErdFL3OyuT6/VWJ302Kdft5ThA4ANFE8IFY6zx0XF987IWir3ti6eWudnyUGlD4lZWxN/kdB3v03effHvF9k7JAAMLPyf7tecHp/fffr0QioeXLl3v9KEBSMB0fpQ72KjaHQxqYw+FGYax9AdszB3pyft/t5wFAqTwNPvbt26dHH31UF198sZePAYbwu+PDSUDh95ROpoICMJFnwccHH3ygRYsW6bHHHtPYsWO9egwwgt/zIJxs8H5nZUyd+wEg3jwLPtrb2zV//nzNnj274Ov6+vqUyWSGfACV8HsehJMN3u+sjElzP+I2cwVAfp7M+XjyySe1f/9+7du3r+hrOzo6tHr1ai+WgRjzcx6Ekw3e7zkcpsz9iOPMFQD5ud7tcvjwYU2bNk07duwYrPW48sorNXnyZK1bt27E6/v6+tTX1zf4eSaTUXNzM90ucIUfbZ12a2+xDd5u7c03h8Or7hO/n5fv+SbPXAFQuUBbbZ966il9/vOfV3V19eDX+vv7lUgkVFVVpb6+viHfG45WW4SR0w3e70xAUJmHIGeuAPBXoMHHyZMn9atf/WrI12688UZdcMEFuuOOO9TW1lbw5wk+EFZON3i/h20FMdwrqJkrAPznZP92veajrq5uRIAxZswYNTQ0FA08gDBzeidNdVXC1w3X7+dJdNsAyI2L5QAXBbHBm8ykbhsA5vAl+PjpT3/qx2OAUOnPWnrh0HF1vvVbSQNBy+Wfbij7KMTEO1NM6bYBYBYyH0AAth84qpU/+IXe//Djwa89/Oyb+v2zztT9//Mix0Wgpray2jNXlm3ar4RyF+O6OXMFQDh4frcLgKG2Hziqv920f0jgYXv/w4/1t6fdBVPq7yvlXpmg2DNXUsmhRyup5GjabIGY4lZbwEf9WUsz79+pnkxfwdc1lth+GqZWVhOPhQC4J9BuFwD57e0+UTTwkH53F0yx4lUn98oEXQhLMS4AG8cugI+ctJSW8lpaWQGEEcEH4CMnLaWlvJZWVgBhRPAB+OjSlnFK1dcUfV1jie2nditrvsqJhIPfBQB+IfgAfFRdldDXr5tU9HWltp/araySRgQgtLICMBXBB+CzuW2N2rB4qn7/rDNHfG/sWWdqg8P2U1pZAYQNrbZAQOIw4RRAfAR6q22lCD4AAAgf5nwAEUJGA0DUEHyEGJtS9Jl6ZwsAVILgI6TYlKLPvrNl+LmofWcLxaQAwopulxAy/SIxVK4/a2n1toM5r6G3v7Z620H1Z40q2QKAkhB8hAybUjw4ubMFAMKG4CNk2JSC15+11HnouLZ2HVHnoeOeBHrc2QIgyqj5CBk2pWD5VWvDnS0AoozMR8gEtSn58W7fdH7W2nBnC4AoI/MRMvam1JM+lbPuI6GBsdpubkp01hSvtUlooNZmTmvKlXZn+86WZZv2KyENeS53tgAIOzIfIeP3RWJ01gwIotaGO1sARBWZjxCyN6Xh2YiUy9kIv9/tB6WUYW1B1drMbWvUnNYUw+QARArBR0j5sSk5ebc/fWKDa8/1U6lHSkEWgFZXJUL79wWAXAg+QszrTSnqnTVOJogGUWsDAFFFzQfycuvdvomdMk6HtfldawMAUUbmA3m58W7f1E6Zco6U/Kq1AYCoI/iIKDduvK203dOti9G8uL233CMlCkABoHIEHxHkZrah3Hf7bnXKeJU5qeRIiQJQAKgMwUfEeHENeznv9t3olPHySnkKSAEgOBScRkixbIOl8m+8td/tL5x8tqZPbCh6zFBpp4zXt/dSQAoAwSH4iJBi2QZpINvw8K43PF9LJcca/VlLjz/f7flEUSaIAkAwOHaJkFKzDf/ykzf0mVSdp5truccauWo8Cql0xggFpADgPzIfEeJkumYlRxalKOdYI989MoW4MVHU6ZESAKAyBB8RYmcbSuH2JWi5ODnWKFTjkQtXygNAeHHsEiF2tuFvN+0v6fV+jEUv9VijlHoVGwWhABBuBB8RM7etUbfP/kP9y09eL/paLy5By6WUuRhOAiEmigJAuBF8RNAtV/0PPbH3V+rJ9OV9TVVCeq/3Ix9XVVipgdCq+Rfqr2e2kPEAgBCj5iOCqqsS+vp1k0YUep4ua0ntm/dr+4Gjvq2rELteJd+a7RoPAg8ACD+Cj4ia29aoR740RcX2aa+7Xgo5/bbbvd0ntGo+Q78AIA44domwsWNqVCiuKGXEuVfy3dly0xUt+uErR7k1FgAijOAjwiodce6VQne2bHyuW498aYrGjqlRT+aUTnzQp3FjRilZO0r9Wcv1zIcXN+YCAAoj+IiwSkace6WU227XPP1LrZrfqge2v+r6bban8+rG3HIQBAGIk4RlWcEc+OeRyWSUTCaVTqdVX18f9HJCrT9radbaXUVHnO+54yrfNrrOQ8f1xcdeKOtn7RW6ce9KvuyLm89wshZTgiAAKJeT/ZuC0wgz8ebWSo543LjNVvL+xlwn8o2U70mf0rJN5nQjAYCbCD4i5vQOks5DxzWnNWXUza2VHvG4cZttsWmqbjyjFCYFQQDgJ2o+IqRQ+n7PHVcZUVNQ7LbbUlWSQTGlENdJEOR3NxIAeInMR0QUS9/vONhjxM2tpRwFlaKSDIophbimBEEA4DeCjwhwO30//OjG7bR/odtuv/2lKSVNOq3kNttSp6l6fWOuKUEQAPiNY5cIcDN971fnRaHbbquqElq2ab8S0pCAyq0iWTv7ku8ZlqR5bQNr8/J4qtgRlN2N5HUQBAB+I/MRAW6l7/3uvLBvux1+FFQoM+JWkWy+ZyQ+iTO++/zb+uJjL2jW2l2edZyY2I0EAH5gzkcElDo744mll+fNfNgzQfJlUIKYCeLH4C37GT852KPvPP/2iO/7MfeDOR8AosDJ/s2xi4+82kzdSN+b2HlhZ0a8fsalLeO04vtdOb9vT11dve2g5rSmPAm8Ch1BAUAUEXz4xMt3t8VqGKTi6fs4d16YEHj5EWgBgCmo+fCBH7UUldZJxLnzIs6BFwAEwfXMR0dHh37wgx/o1VdfVW1trWbMmKG1a9fqM5/5jNuPCoVSLlJzK6VfSfo+zp0XcQ68ACAIrmc+du/erfb2dr3wwgvasWOHPv74Y11zzTXq7e11+1Gh4Oco70pqSuLceWHK3A8AiAvXMx/bt28f8vnjjz+u8ePH6+WXX9YVV1zh9uOM51dK342aEvvoZvjvSUW888KNmhkAQOk8LzhNp9OSpHHj4vmu0Y+Ufr7r4e2aEidtonHtvIhr4AUAQfB0zkc2m9V1112n999/X3v27Mn5mr6+PvX19Q1+nslk1NzcHJk5H/b8jGK1FOXOzzBxPkeY+TFbBACiyMmcD0+7Xdrb23XgwAE9+eSTeV/T0dGhZDI5+NHc3OzlkjyT7z4Ur2spTLkePiryTV0FALjHs2OXW265RT/60Y/03HPP6Zxzzsn7ujvvvFMrVqwY/NzOfIRJsXoLL1P6tIkCAMLG9eDDsiz93d/9nbZs2aKf/vSnamlpKfj6mpoa1dTUuL0M35Rab+FVLQVtogCAsHE9+Ghvb9fmzZu1detW1dXVqaenR5KUTCZVW1vr9uMC5XSGhxdTLOM8nwMAEE6u13ysX79e6XRaV155pRobGwc/vve977n9qMCZUG8R5/kcAIBw8uTYJS5MqbegTRQAECZcLFcBk+ot4jqfAwAQPgQfFTCt3sL0m1GZoQEAkAg+KsJY7tK5Mf4dABANng4Zi4NKr7KPA7sdeXhxrt2OvP3A0YBWBgAIApkPF7hdbxGl4wmn7cgAgOgj+HCJW/UWJh5PVBIMOWlHNrleBQDgntgFHyZnFdy8ndbNNVUSDJnSjgwAMEesgg8Tswo2E48n3AiGTGpHBgCYITYFp6YXPZowLfV0xYIhaSAYsm/vzcduR84XLiU0EAAy/h0A4iMWwYdbG6mXTDuecCsYYvw7AGC4WAQfpmUVcjHteMLNYIh2ZADA6WJR82FaViEX06aluh0MMf4dAGCLRebDtKxCLqYdT3hRq2G3Iy+cfLamT2wg8ACAmIpF8BGWokeTjidMC4YAANGRsCwruCrLHDKZjJLJpNLptOrr6137vXa3i5T7DhaTag9MmkVicnsyAMAcTvbv2AQfEhtpuUwKhgAAZnKyf8ei4NR2etFjT/r/6kTvRxr3ezVK1o5Sf9ZiQ83DrdHxAABIMQs+pIGNNP1/P9ID/+c1MiAAAAQgFgWnpzN90ikAAFEXq+AjDJNOAQCIulgFH2GYdAoAQNTFquYjDJNOvUC3CgDAJLEKPryedGriJk97MQDANLEKPry8P8XETd4urh3+z2oX15o0WA0AEB+xqvnwamS4iR00FNcCAEwVq+BDcv/+FFM3eYprAQCmitWxi83N692dbPJ+TgmNa3EtAMB8sQw+JPdGhpu6yXtdXAsAQLlid+ziNlM3ebu4Nl8uJ6GBgthyimsBAKgEwUeFTN3kvSquBQCgUgQfFTJ5k3e7uBYAADckLMsyqtcyk8komUwqnU6rvr4+6OWUzMQ5HzYTh58BAKLFyf5N8OEiNnkAQFw52b9j2+3iBbc6aAAAiDJqPgAAgK8IPgAAgK8IPgAAgK8IPgAAgK8IPgAAgK8IPgAAgK8IPgAAgK9iM+eDAWAAAJghFsGHyaPPAQCIm8gfu2w/cFTLNu0fEnhIUk/6lJZt2q/tB44GtDIAAOIp0sFHf9bS6m0HlevyGvtrq7cdVH/WqOttAACItEgHH3u7T4zIeJzOknQ0fUp7u0/4tygAAGIu0sHHuyfzBx7lvA4AAFQu0sHH+LrRrr4OAABULtLBx6Ut49SYHK18DbUJDXS9XNoyzs9lAQAQa5EOPqqrErp7QaskjQhA7M/vXtDKvA8AAHwU6eBDkua2NWr94qlKJYceraSSo7V+8VTmfAAA4LNYDBmb29aoOa0pJpwCAGCAWAQf0sARzPSJDUEvAwCA2Iv8sQsAADALwQcAAPCVZ8HHI488ovPOO0+jR4/WZZddpr1793r1KAAAECKeBB/f+973tGLFCt19993av3+/LrnkEl177bV69913vXgcAAAIEU+CjwcffFBLly7VjTfeqNbWVm3YsEFnnXWWvvvd73rxOAAAECKuBx8fffSRXn75Zc2ePft3D6mq0uzZs9XZ2Tni9X19fcpkMkM+AABAdLkefPz2t79Vf3+/JkyYMOTrEyZMUE9Pz4jXd3R0KJlMDn40Nze7vSQAAGCQwLtd7rzzTqXT6cGPw4cPB70kAADgIdeHjP3BH/yBqqurdezYsSFfP3bsmFKp1IjX19TUqKamxu1lAAAAQ7kefIwaNUp/9Ed/pJ07d+r666+XJGWzWe3cuVO33HJL0Z+3LEuSqP0AACBE7H3b3scL8WS8+ooVK7RkyRJNmzZNl156qdatW6fe3l7deOONRX/25MmTkkTtBwAAIXTy5Eklk8mCr/Ek+PjLv/xL/eY3v9Fdd92lnp4eTZ48Wdu3bx9RhJpLU1OTDh8+rLq6OiUSXPxWjkwmo+bmZh0+fFj19fVBLyc2+LsHh799MPi7B8fEv71lWTp58qSampqKvjZhlZIfQahkMhklk0ml02lj/qOMA/7uweFvHwz+7sEJ+98+8G4XAAAQLwQfAADAVwQfEVRTU6O7776bFmaf8XcPDn/7YPB3D07Y//bUfAAAAF+R+QAAAL4i+AAAAL4i+AAAAL4i+AAAAL4i+IiQjo4Offazn1VdXZ3Gjx+v66+/Xq+99lrQy4qd+++/X4lEQsuXLw96KZF35MgRLV68WA0NDaqtrdVFF12kl156KehlRV5/f79WrVqllpYW1dbWauLEiVqzZk1Jd3rAmeeee04LFixQU1OTEomEnnrqqSHftyxLd911lxobG1VbW6vZs2frjTfeCGaxDhB8RMju3bvV3t6uF154QTt27NDHH3+sa665Rr29vUEvLTb27dunRx99VBdffHHQS4m89957TzNnztSZZ56pZ555RgcPHtQ3v/lNjR07NuilRd7atWu1fv16Pfzww/rlL3+ptWvX6oEHHtC3vvWtoJcWOb29vbrkkkv0yCOP5Pz+Aw88oIceekgbNmzQiy++qDFjxujaa6/VqVOnfF6pM7TaRthvfvMbjR8/Xrt379YVV1wR9HIi74MPPtDUqVP17W9/W/fcc48mT56sdevWBb2syFq5cqWef/55/exnPwt6KbHzuc99ThMmTNB3vvOdwa/92Z/9mWpra7Vp06YAVxZtiURCW7ZsGbwx3rIsNTU16e///u/1D//wD5KkdDqtCRMm6PHHH9cNN9wQ4GoLI/MRYel0WpI0bty4gFcSD+3t7Zo/f75mz54d9FJi4Yc//KGmTZumv/iLv9D48eM1ZcoUPfbYY0EvKxZmzJihnTt36vXXX5ckvfLKK9qzZ4/mzZsX8Mripbu7Wz09PUP+n5NMJnXZZZeps7MzwJUV58mttgheNpvV8uXLNXPmTLW1tQW9nMh78skntX//fu3bty/opcTGW2+9pfXr12vFihX66le/qn379unWW2/VqFGjtGTJkqCXF2krV65UJpPRBRdcoOrqavX39+vee+/VokWLgl5arPT09EjSiBvjJ0yYMPg9UxF8RFR7e7sOHDigPXv2BL2UyDt8+LBuu+027dixQ6NHjw56ObGRzWY1bdo03XfffZKkKVOm6MCBA9qwYQPBh8e+//3v6z/+4z+0efNmTZo0SV1dXVq+fLmampr426MkHLtE0C233KIf/ehHevbZZ3XOOecEvZzIe/nll/Xuu+9q6tSpOuOMM3TGGWdo9+7deuihh3TGGWeov78/6CVGUmNjo1pbW4d87cILL9Svf/3rgFYUH1/5yle0cuVK3XDDDbrooov0V3/1V7r99tvV0dER9NJiJZVKSZKOHTs25OvHjh0b/J6pCD4ixLIs3XLLLdqyZYt27dqllpaWoJcUC1dffbV+8YtfqKura/Bj2rRpWrRokbq6ulRdXR30EiNp5syZI1rJX3/9dX3qU58KaEXx8eGHH6qqauj2UV1drWw2G9CK4qmlpUWpVEo7d+4c/Fomk9GLL76o6dOnB7iy4jh2iZD29nZt3rxZW7duVV1d3eCZXzKZVG1tbcCri666uroRdTVjxoxRQ0MD9TYeuv322zVjxgzdd999+sIXvqC9e/dq48aN2rhxY9BLi7wFCxbo3nvv1bnnnqtJkybp5z//uR588EH9zd/8TdBLi5wPPvhAb7755uDn3d3d6urq0rhx43Tuuedq+fLluueee3T++eerpaVFq1atUlNT02BHjLEsRIaknB//+q//GvTSYueP//iPrdtuuy3oZUTetm3brLa2Nqumpsa64IILrI0bNwa9pFjIZDLWbbfdZp177rnW6NGjrU9/+tPWP/3TP1l9fX1BLy1ynn322Zz/X1+yZIllWZaVzWatVatWWRMmTLBqamqsq6++2nrttdeCXXQJmPMBAAB8Rc0HAADwFcEHAADwFcEHAADwFcEHAADwFcEHAADwFcEHAADwFcEHAADwFcEHAADwFcEHAADwFcEHAADwFcEHAADwFcEHAADw1f8HBKlHSzQf0w0AAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"#based on https://stackoverflow.com/a/75562058/8508004 & https://stackoverflow.com/a/26447505/8508004\n",
"import numpy as np\n",
"x_data, y_data = np.repeat(np.linspace(0, 9, 100)[None,:], 2, axis=0) + np.random.rand(2, 100)*2\n",
"import matplotlib.pyplot as plt\n",
"fig, ax = plt.subplots()\n",
"ax.scatter(x_data, y_data)\n",
"\n",
"# next make the fit BUT DON'T ADD IT YET\n",
"#z = np.polynomial.polynomial.polyfit(x_data, y_data, 1) comment https://stackoverflow.com/questions/26447191/how-to-add-trendline-to-a-scatter-plot#comment100202969_26447505 seems wrong\n",
"z = np.polyfit(x_data, y_data, 1)\n",
"p = np.poly1d(z)"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "5b57b5fd-5af6-4c14-940e-607e5a7f8bb1",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAh8AAAGdCAYAAACyzRGfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAABB5klEQVR4nO3de3xU1bn/8e8kQBJpMnIpuShIpCqGq0hFRK1FVKxQ1B4toD3UnpetFqtIL6ItKvWC2JZSreKlF3t+NKBVEaEVq6goGi6CkeZQAREUkYCCJhAl0Jn5/REmZDJ77vs2ez7v14s/Zs3O7JXYV/czaz3Ps3yhUCgkAAAAm+Q5PQEAAJBbCD4AAICtCD4AAICtCD4AAICtCD4AAICtCD4AAICtCD4AAICtCD4AAICtOjg9gfaCwaA++ugjFRcXy+fzOT0dAACQhFAopH379qmiokJ5efHXNlwXfHz00Ufq2bOn09MAAABp2L59u4499ti417gu+CguLpbUMvmSkhKHZwMAAJLR2Nionj17tj7H43Fd8BHeaikpKSH4AAAgyySTMkHCKQAAsBXBBwAAsBXBBwAAsBXBBwAAsBXBBwAAsBXBBwAAsBXBBwAAsBXBBwAAsJXrmowBAJArAsGQVm/dq937DqhHcaFOq+yq/Dzvn2tG8AEAgAOW1u3UjMUbtLPhQOtYub9Qt42t0uj+5Q7OzHpsuwAAYLOldTt17bx1EYGHJNU3HNC189Zpad1Oh2ZmD4IPAABsFAiGNGPxBoUM3guPzVi8QYGg0RXeQPABAICNVm/dG7Xi0VZI0s6GA1q9da99k7IZwQcAADbavS924JHOddmI4AMAABv1KC409bpsRPABAICNTqvsqnJ/oWIV1PrUUvVyWmVXO6dlK4IPAABslJ/n021jqyQpKgAJv75tbJWn+30QfAAAYLPR/cs198ohKvNHbq2U+Qs198ohnu/zQZMxAAAcMLp/uc6rKqPDKQAAsE9+nk/D+3Rzehq2I/gAAMBArp67YgeCDwAA2snlc1fsQMIpAABtePrclX//Rqr2SQsrHJ0GKx8AAByW6NwVn1rOXTmvqiy7tmD+0yQ98aUjr79wNoBi5QMAgMM8ee7KB09FBh6SNO59Z+ZyGCsfAAAc5qlzV4IBaclJ0v4tR8aOGy+NmO/cnA4j+AAA4DCzzl1xvFLmk1XSP0+PHBu9Vuo6xL45xEHwAQDAYeFzV+obDhjmffjU0oU03rkrjlfKLB8n7Xj2yOujB0gX1ko+92RauGcmAAA4LNVzVwLBkGq27NGi2h2q2bJH/1jvYKXM/vdaKlnaBh5nPyt9Y72rAg+JlQ8AACKEz11pv3pR1m71wmiFI88nZyplaqdJG2a1GfBJlzdJHYrMvY9JUg6FXn31VY0dO1YVFRXy+Xx65plnIt4PhUK69dZbVV5erqKiIo0aNUqbN282a74AAFhudP9yrbhppOZffbp+N36w5l99ulbcNDIi8DBa4QgaRR6HWVIp07y3ZbWjbeBx6v3SxKBrAw8pjeCjqalJgwYN0gMPPGD4/r333qv77rtPDz30kFatWqXOnTvrggsu0IEDWZAZDADAYeFzV8YNPkbD+3SL2GqJ1QskGaZVyrz7qPRUu3NhLt0tnXSdOZ9voZS3XS688EJdeOGFhu+FQiHNmTNHv/jFLzRu3DhJ0v/+7/+qtLRUzzzzjMaPH5/ZbAEAcFiiXiCJJFtRE1OgWXqyixT44sjYSTdIp87J7HNtZGoGytatW1VfX69Ro0a1jvn9fg0bNkw1NTVm3goAAEeku3LhU0vVS7xKmYQ+el56vDAy8BizMasCD8nkhNP6+npJUmlpacR4aWlp63vtNTc3q7m5ufV1Y2OjmVMCAMBU6axcGFXKpCQUkp4fJu1dc2Ss7Hzp60slXxa1eT/M8dqbmTNnyu/3t/7r2bOn01MCAOSQ9uWygXhZozrSCyTeI799fFHmL9TcK4ek1+fj0/XS/LzIwGPUa9LI57My8JBMXvkoKyuTJO3atUvl5Uf+wLt27dLgwYMNf+bmm2/W1KlTW183NjYSgAAAbJFOQ7BwL5Br562TT5GlteFQ4PcTTlGXzgWZdzhdeZX03mNHXhcdI43bJuVld6cMU1c+KisrVVZWpmXLlrWONTY2atWqVRo+fLjhzxQUFKikpCTiHwAAVotVLptMQ7BwL5Ayf+QWTHiF4xsDKwwrZZL2+Y6WEtq2gccZ1dIlH2Z94CGlsfKxf/9+vfvuu62vt27dqtraWnXt2lW9evXSlClTdOedd+qEE05QZWWlpk+froqKCl188cVmzhsAgLTFK5dNtiHY6P7lOq+qzPwzXP5vpvT2LZFjlzVKHYsz+1wXSTn4ePPNN/X1r3+99XV4y2TSpEl67LHH9LOf/UxNTU36/ve/r88++0xnnnmmli5dqsLCDEuLAAAwSaJy2bYNwYb36RbzunAvEFMc2if9rd3q/6C7pH63GF+fxVIOPs455xyFQrGTcXw+n375y1/ql7/8ZUYTAwDAKi9uMK7AbM+0hmCJbFsgvTEhcuziD6WjjrHn/jbL/o0jAABSEAiGtLB2R1LXZtwQLOFkDkqPF0SOVf63NPwv1t7XYQQfAICcsnrrXu1tOpTwum6dO2XWECyRbfOlNyZGjl1YK3UZZN09XYLgAwCQU5LdShk3uML802elloZh8w2KTScEJJ/j7bdskRu/JQAAhyW7lXJeVZn5N/9kZXTgcer90sRQzgQeEisfAIAcE+5QWt9wIObJtHk+6dOmg+be+Jle0ufbI8cub5I6HGXufbJA7oRZAADoSIfSeIIhaXJ1/EZjSdu7rqVhWNvA4yvfb1ntyMHAQyL4AABkmVTPYjEyun+5Hph4StQZLO3NWLwhrc9vVe2Tlp4aOXbxh9JpD6f/mR7AtgsAIGukcxZLLF06FyheXJFsozFDn38oPWNwTtnEDAIZD2HlAwCQNjNWIZKVyVksRpKtekm50Vi1LzrwOP0vBB5tsPIBAEiLmasQiZhxFkt7yVa9JN1oLHhIWtApepygIworHwCAlJm9CpFIKmexJCtc9RIrVPGpJZhKqtHYa9+KDjy6DyfwiIHgAwCQkkSrEJIJiZrtWLFF0rbqpX0AEn5929iqxCsp1T5p+9ORY5c3See/kfRc7GLnNlk8bLsAAFJi1omwqTB9i+Sw0f3LNffKIVHbR2XJbB9tflhac030uEtXO+zcJkuE4AMAkBLLEjXjSNQYzKeWgCGds1hG9y/XeVVlWr11r3bvO6AexS2fE3fFo9rgvTEbpZITU76/HcLbZO3/duFtsrlXDrE1AGHbBQCQEqtWIeIxbYskzucP79NN4wYfo+F9usX+nHDDsPYmhlwbeDixTZYIwQcAICWmJmqmILxFUuaPDGrK/IX2fHM3ahh25hMxt1nckl9hRbJupth2AQCkJLwKce28dfJJEd+ozViFiCetLZJMHfxMerJL9Hic3A435Vc4sU2WCCsfAICUObkKkWiLxNQVh0W9owOPr1yTMPCwsww5ESe2yRJh5QMAkBZHViESMG3FIRSS5ht8P58QkHyxv7db0QwtU1Ym66aLlQ8AQIRUVg6STtS0gWkrDmunRgceBd1bVjviBB6SO/MrrE7WTQcrHwCAVm7KVUiFaSsORpUsl34sFXZPah5uzK+QMuxnYgGCDwCAJPf1gkhFxo3PPlwkvXpx9HiKDcPcmF8R5qZtMoIPAIArcxVSkdGKg9Fqx/krpe7DUp6HG/Mr2gpvkzmNnA8AgCtzFVKR1orDF/WxG4alEXhI7syvcCOCDwCAa3MVkpVy47Pnh0sL220hnfWUKeeyON4MLQuw7QIAcDxXIRAMZZSLkHTjs9Ahqbog+gNMPgzOTfkVbkTwAQBwNFfBrAqbhBUd/7lfenxm5A/1+4U06I6MfwcjbsmvcCNfKBRy1dm/jY2N8vv9amhoUElJidPTAYCcEa52kYxXDqzYMohVYZPJPQ1XURYYZBmMPyjldUxr3oiWyvObnA8AgCT7cxWsOm01ovFZp1ejA4/uZ7RssxB4OIZtFwBAKztzFTLuzZGIUSXLJfVSUWnqnwVTEXwAACLYlatgWYVNw7+lv1dFj5ucVIr0EXwAABxhSYXN452lwOeRY2k2DIN1CD4AAI4wtcLmYIP05NHR46x2uBIJpwAAR5jWDXTF5dGBx2mPEni4GCsfAADHZHTaaigozc+PHp8QlHw083Izgg8AgKPSqrDZeJ+09obIsT7/Iw37g7WThSkIPgAAjkupwsaohPby/VKHzuZOCpYh5wMAkB12LY8OPIoqWnI7CDyyCisfAAD3M1rtGPuuVNzH/rkgYwQfAAD3atouLeoVPZ5llSyZntrrNQQfAABLpf3g/Xt/qeH/IsfOeU6qGG3ufSxm1qm9XkLwAQAwhdHD/4UN9ak/eAMHpMeLosfjrHa49QEf69Te+oYDunbeOksO7MsGvlAo5Kq1q1SO5AUAuIPRw//oozrqs88PRV0bXoswfPCunSJt/F3k2KC7pH63xL230QM+7n1sEAiGdOasl2Ienhfu4LrippGuWKHJVCrPb6pdAAAZCT/82z9kjQIPSa1BwozFGxQIHn4VCrUklbYPPMYfiht4BIIhzVi8wbA9u+F9bJTKqb25huADAJC2QDCk2581fvjHE/Hg3TZfmt/ucVR2fss2S1787AA3P+AtO7XXA8j5AACk7fcvbVZ9Y/oPz+GrukcPfusTqSC5hmNufsBbcmqvR7DyAQBIy9K6nfrti5vT+tnhnddr28AxkYMlJ7esdiQZeEjufsCHT+2Nlc3hU0tSbFKn9noMKx8AgJSFcy3SERV0SNKFb0ldBqf8WeEHfH3DAcOtn3BSpxMP+PCpvdfOWyefFDG/lE7t9SBWPgAAKUuUa2GkR4c9xoHHxFBagYd05AEvKWqFwQ0P+PCpvWX+yJWXMn9hzpbZSqx8AADSkEoOxdFHdVTtVy6IfmPA7dKA2zKeS/gB377Ut8wFfT6kNE/t9TiCDwCwkVu7cKYq2RwKf4HPOPCYEJR85v3ebn/Ap3Rqbw4g+AAAm7i1C2c6EuVaSNJ9ve7VN49+NXKwyxDpwrWWzIkHfPYg5wMAbBCrEVe4zfbSup0OzSw9bXMtjGwbOCYq8Ah8q9GywKO9QDCkmi17tKh2h2q27HGkyRhiY+UDACyWqAunTy1dOM+rKnPNNkEywrkWtyz8l/Y2tXQz/VGP+fpx2V+jru29fonmDzuo4X2sn5eXVpi8ipUPALCYm7twZmp0/3JNH9NPUstqR/vA49yNc9V7/RJJ9jT68toKk1eZHnwEAgFNnz5dlZWVKioqUp8+fXTHHXfIZefXAYBt3NyF0wzHd9hiWELbe/0SbWnu2fra6kZfbj7nBZFM33aZNWuW5s6dq7/85S/q16+f3nzzTV111VXy+/26/vrrzb4dALiem7twZuyZ4zTo8w8ihu786Hv6wyeXtr62q9FXKitMJKY6y/Tg44033tC4ceN00UUXSZJ69+6t+fPna/Xq1WbfCgCyQiZdOF1bmvufJumJL0UNV65f4lgnT6+vMHmJ6cHHGWecoUceeUSbNm3SiSeeqLffflsrVqzQ7NmzDa9vbm5Wc3Nz6+vGxkazpwQAjkq3zbZrEydXXS1t+UPk2JDZWvqf8Sp737lGX55eYfIY04OPadOmqbGxUX379lV+fr4CgYDuuusuXXHFFYbXz5w5UzNmzDB7GgDgKql04QwEQ/r9S5sND20LJ0460po7FJLmG6QKTghIvjydFwypuKCjat77RFJLz43Tj+9m20qNm895QSRfyORM0AULFuinP/2pfvWrX6lfv36qra3VlClTNHv2bE2aNCnqeqOVj549e6qhoUElJSVmTg0AHJdoG2Vp3U7d/uyGuMfUhx+iK24aafmDPTzfwu1/0Skf/TjyzWMvkc5+WlJmqzRmbi2Fq10k4xWmXD5PxWqNjY3y+/1JPb9NDz569uypadOmafLkya1jd955p+bNm6d33nkn4c+nMnkAyAbJPlzDD85k/095/tWnW5o4GQ4oao4bFf3mf30qdTq69TqjeSfzwLdia8m121Uel8rz2/Rtl88//1x5eZHLcvn5+QoGg2bfCgBcL9kHYbwy0VisTJxcWrdTc59aoJoTpkaMNwc7qm/dQs0d+IVG9z86owZqsYKWTLeW3H7OCywIPsaOHau77rpLvXr1Ur9+/fTWW29p9uzZ+t73vmf2rQDA1VJ5uKZ1RH1xoSXVMIFgSF+rPV6jT4icz3kbH9Dm5uMiAop0y1ut7vrKOS/uZnrwcf/992v69On64Q9/qN27d6uiokI/+MEPdOutt5p9KwCwTaoP+VQfrqmsYoRzPj5tataZs14yd3vh0D7l/61ERe3ySsNdSsPzDwcU6Za30pMjt5kefBQXF2vOnDmaM2eO2R8NAI5IJ4cg1YdrquWf3xxUrsnVb6W9ZWEYTL3za6n2ZxHXXbPtZi1tHGH4GeGfTUb76+jJkds4WA4A4kg3LyHVh2syR9RLUllJgW4d0093/D39LYv2wVSeAnpv4Lio69qudhgJBy3plLfSkyO3cbAcAMSQyVkhqT5c2x5RH2sz58ZRJ+r1aeeqS+dOaR9U1/7gtVElq6ICj8CIJzX8/RdjzsOnlpWf8NZTrHnHa6AWDlqSuQe8h+ADAGLI5DTadB6u4UZkZf7IwKXcX6iHrhyiG0adkFJ+SPvr2gdT2waO0R963xFxzZnvPy/1vDSlgCLWvMv8hTFXhtINWuANbLsAQAyZ5CWk21I9mTLRdLcswsFU38KtWnrijyLeu3fnf+vBjy+XdEirt+5NqSNrsvNuL9V7wDsIPgAghkzzEtJ9uCYqE003z2L3vgN6os/PdFrnDRHjA+sWqDH4pYjrwvNPJaBIp7yVnhy5ieADAGIw46wQKx6uaa2qfFGvcRuOlTofGVry2Zm67oNpUZ/fNpiyo18GPTlyDzkfABCDWXkJ4YfruMHHaHgfcw5aSynPYu0UaWHkKsvZ7zwaFXiQ5Am7sPIBAG20739xXlWZa/MSEq6qBA5IjxdF/MwXHctVtfbRqM8iyRN2IvgAgMPiNRNbcdNIV+YlxNyyePcP0uqrI8fOfUVFpV/T3JOjf083BFPIHaafapspTrUF4IRMTmZ1lVBImm+woz4hKPmOBEtWnAmD3OboqbYAkG2sPuTMNrtekZZ9PXJs2B+lPtEHe5qV5EkQg3QQfADIeZ445GzhMdIXH0WOffsLKd+69uTpnHkDSFS7AEB2H3L28etStS8y8Og7VZoYsjzwaNumPSx85s3Sup2W3RvZj5UPADkvaw85qzbY3rikXioqtfS2ntmmgmNY+QCQ87LukLPPdxgHHhNDlgceUmZn3gASwQcAZNchZ8+dIj1zbOTY15a0BB42yeptKrgCwQcAKL2TWW31ny9aVjs+rY0cnxiSjrlIgWBINVv2aFHtDtVs2aNA0LpgJGu3qeAa5HwAwGGuPeTshbOlj1+LHBswQxpwqyT7q07MOPMGuY0mYwDgZka5HeMPSXkt3x2dao4Wvq9kfLCdK1aLYKtUnt9suwCAG22YFTup9HDgkajqRGqpOrFiC8b121RwNbZdAHhOtnTdjDlPo6Bj3AdS554RQ043R3PtNhVcj+ADgKdkS9dNo3mO7rFZD5XdGH1xjEoWN1SdmNWmHbmF4AOAZ8TKfwh33XTLdoDRPLcNHBN94fkrpe7DYn4OVSfIVuR8APAEJ/MfUtF+nt07fGoYeATGB+MGHlIWNkcDDiP4AOAJ2dJ1s+083+l/qd6s+k7E+7/fdbl6r1+S1Dyzqjka0AbBBwBPcEP+Q7L3z1NA2waOUWHewYj3eq9frF/v+u/W65JB1QmyETkfADwhW/IfTml8QO8NnBkxtvFAL12w6cGIsVTmSdUJsg3BBwBXSbdMNiu6blb71KvdUL+6J9QUPKr1dbrzpOoE2YTgA4BrZFImG85/uHbeOvlk3HXTsfyHrfOkmu9EDVeuX+KueQI2IecDgCuEy0/bJ42Gy2SX1u1M+BmuzH+o9kUHHt/cqqUDP3LXPAEbcbYLAMcFgiGdOeulmNUq4a2IFTeNTGpFwBUdThs2SH/vFz3epmGYK+YJmCSV5zfbLgAcZ3abcMfzH4zao5/+mHT8pIghx+cJOITgA4DpUv1Gny1lsgkd2if9zeAbX4z26ECuIvgAYKp/rN+pXyyq096mIz0sEiWNZkuZbFz/GCx99nbkWO8rpDPmOTIdwM0IPgCYZuY/NujhV7dGje9McLZKVpTJxhIKSfMNcvfHH5LyrP+/WPJGkI0IPgCY4h/rPzIMPMJCajlb5byqsqiHo6vLZONZf5tU98vIMV8HacIhW26fLSf4Au1RagsgY4FgSL9YVJfwunhnq7iyTDaeal904HHJTlsDj0xLk50UCIZUs2WPFtXuUM2WPY4f+Ad7sfIBIGOrt+7V3qbkHrrxkkazok34zhekl8+PHrcxqTTRCb4+xV5lcgNWbEDwASBjqVShJEoadXX5qVEJ7bmvSKVfs3UaZpcm2ym8YtM+cKpPkBcEb2HbBUDGkq1C6da5kzuTRhM58Ilx4DExZHvgIWVvaXKiFRupZcWGLRjvI/gAkLFwtUoid4zr78ptgLiWnSs9/eXIsUF3O9q7I1tLk1NZsYG3EXwAyFi4WiVeWPGDsyv1jYFZtJwe/E/LaseulyLHJ4akfjc7M6fDwsFerL+3Ty05FG5bZcrWFRuYj+ADgCnC1SrtV0C6du6oByeeopu/UeXQzNKw6UFpQcfIsZNudE2n0nCwJykqAHFzaXK2rtjAfCScAjBNVlSrJLKgQAoejBz79gEpv8CZ+cQQDvbaV42UubhqJKubycFUnGoLAJL0WZ30jwGRY0cPkL6x3pn5JCnbOpyGq10k42ZyVLtkr1Se3wQfALDs3Ojcjv/aK3Xq4sx8PI4+H96UyvObbRcAuevAbunp0sixYy+Wzl7oyHTiybYVjng8sT2HjBB8AMhNb/1M+vevIsfGbJJKTnBmPnF4caXA1c3kYDmqXQDklv80tZTQtg08CktbKllcGnhk8xkugBGCDwC5o+a70hNfihw79yXp0npHppMIHUHhVWy7APC+UEiab/Bda0JQ8rk3zyCbz3AB4mHlA4C3vftodODR69st2ywuDjwkOoLCu1j5AOBdRofBXdYgdcyOMn46gsKrWPkA4D2frIwOPDoUt6x2ZEngIWXvGS5AIpYEHzt27NCVV16pbt26qaioSAMGDNCbb75pxa0AIFK1T/rn8Mixi/4tXd7ozHwykK1nuACJmB58fPrppxoxYoQ6duyo5557Ths2bNBvfvMbdelCp0AA0QLBkGq27NGi2h2q2bIn/cqNL+qNt1kmhiR/38wm6aDwGS5l7Q7sK/MX0oocWcv09urTpk3T66+/rtdeey2tn6e9OpA7TGue9fxwac/KyLGzn5GOHWfORF3ASx1O4U2Onu1SVVWlCy64QB9++KGWL1+uY445Rj/84Q919dVXJ/XzBB9Abgg3z2r/f0ApHTAWOCg9bnDa7MTs7HtBgIFs5ujZLu+9957mzp2rqVOn6pZbbtGaNWt0/fXXq1OnTpo0aVLU9c3NzWpubo6YPABvS9Q8y6eW5lnnVZXFfvjW3ixtuCdyrN/PpUF3mjxbe3ixhToQi+krH506ddLQoUP1xhtvtI5df/31WrNmjWpqaqKuv/322zVjxoyocVY+AO+q2bJHEx5dmfC6+Vefbtw8yyi3Y/xBKa+jCbOznymrQIDDUln5MD3htLy8XFVVVRFjJ598sj744APD62+++WY1NDS0/tu+fbvZUwLgMmk3z3rnt9GBR/czWrZZsjTwoIU6cpHp2y4jRozQxo0bI8Y2bdqk4447zvD6goICFRQY7NkC8KxUmmeF8yCGr+oefcE3t0pf6m3u5GxGC3XkItNXPm688UatXLlSd999t959911VV1frkUce0eTJk82+FYAslWzzrE+bmvXT3800DjwmhrI+8JBooY7cZHrw8dWvflULFy7U/Pnz1b9/f91xxx2aM2eOrrjiCrNvBSBLJdM865uDyvWNumM0u/TnEe//z9ZbVbl+iWeOkqeFOnKR6QmnmaLUFtmIEsn0xKrwuPP8rjp345Co63uvXyKpJUAp8xdqxU0js/7vHAiGdOasl1TfcMAw78NLvyu8zdFSWyDXUCKZvtH9y3VeVVlE4Hb6WyfIt/HTiOteaBimq9+f3vraS3kQ4VWga+etk0+KCEBooQ6vIvgAMhCrRLK+4YCunbeOEskk5Of5WgKIUFCanx/1fu/1ixW9OdPCK3kQ4Rbq7YPYMoJYeBTBB5AmUxpluYxj20ervi9teTRiaFtzuc7Z+GiMH2jhpTwIo1Ugtu/gVQQfQJq8ViLp2PaRQcOwwKWfasJv18qn+HkQXjtKvnUVCPA406tdgFzhpRLJ8PZR+2AqvH1kSWXJruUxT6HNLzzasqPkTTtFF0DaWPkA0uSVEklHto+Mgo6x70rFfVpfWpEHQXIw4A4EH0Cawo2yEpVIun1rwNbto6bt0qJe0eMxTqE1Mw+C5GDAPdh2AdKUTKOsbCiRtG37qNoXHXicszRm4BEWzoMYN/gYDe/TLe2tFs5PAdyD4APIQHhroMwfubVS5i/Mmm/Slm8f/acpZm6HKi5I7zNTlMrqDgDrse0CZCjbSyQt3T56fpi0Z3Xk2HHjpRHz05lq2ryUHAx4AcEHEEeyfS+yuUTSkg6boZA032BhdfxBKa9jBrNNj1eSgwGvIPgAYrCjMsItZ8KYWlmybb70xsTo8QS5HVbySnIw4BUcLAcYiFUZEQ4LzMjncGPZZ7LBUMzrjHI7Lv5QOuoYG2YfX/i/qWS8upMtOTqAW6Xy/Cb4ANoJnzIaK0HRjFNGrQ5urFxRMQqazuq+Q/+v4gfRFxusdji52uPGgA/wCk61BTJgdd8Lq5t6WfmANQqatg0cE33hhW9JXQbbOrdkZHtyMOAVlNoC7VhdGWFl2aeVbdLbB03+/H2GgUdgfDBm4GF7C3cDZvQNAZAZgg+gHasrI6wKbqxupNU2aPpz79v0dr8JEe9f/8FP1Hv9EsOgiSZfANpi2wVox+rKCKuCG6u3i3bvO6A8BfTewHFR7/Vev1jhjBWjoMlrJwADyAwrH0A7VrdNDwc3sX7ap5Y8iFSDG6u3iwbveyQq8PjzJ2PVe/0Stf1LGQVNNPkC0BbBB2DAyrbpVgU3lm4XVft03PZfRgyd9K+nNeOjIxUu8YImmnwBaIttFyAGKysjMm3qZVSuasl20UfPS6+MjhjafKCnzt80N6VOqDT5AtAWwQcQh5Vt09MNbuKVq5raJt2oYdi497Vla0eV7UotaLKkhTuArEWTMSCLJNOcTFJmvTT2b5WePT56vE3DsHQbhTnd5wOAdehwCnhQKp1XJaW3XfRML+nz7ZFj574slZ6T2eTbMApc0p4vANegwyngQamWq6a0XfSfJumJL0WPW3AYXPutLFZDgNxDtQuQJSwrV111dXTgMeS3tpxC65aupwDsxcoHkCVML1cNhaT5Bt8/JgQkn/XfS6w+4waAe7HyAWQJU5uTbflTdOBx7CUtqx02BB6StWfcAHA3Vj6ALGFauapRCe1/fSZ18psz0STR9RTIXax8AFkko86re9ZEBx75hS2rHTYHHhJdT4FcxsoHkGXSak5mtNpx0f9J/irrJpoAXU+B3EXwAWShpDuvHtgtPV0aPW5DJUsi6W4jpdvgDIB7EHwAXmW02nHWU1LPS+2fSwypnnFDTxDAG+hwiqzAt90UBA9JCzpFj7tgtSOWZP77JtNangAEcA4dTuEpfNtNwfrbpLpfRo6VjpTOXebMfJKUaBuJniCAt1DtAlejA2YKqn3Rgce3m10feCSDniCAtxB8wLUSfduVWr7tBoLu3U6wxaYHovM7ug5t2WbJN9h+yUL0BAG8hW0XuFaqB6klw3O5I0ZJpWPflYr72D8XC9ETBPAWgg+4ltnfdj2VO7Lzn9LLF0SPuzipNBP0BAG8hW0XuJaZ33Y9lTtS7YsOPM56OmsDj0AwpJote7Sododqtuwx3EYL9wSRFHW2TUqt5QG4AisfcC2zvu16plLii3ppocEKTZYGHVJqq1Gp9gQB4F4EH3ClcG7Ghf3L9KfXt2V0kJoVuSO2M8rtKDtfGvm8/XMxSay+HeHVKKO+HWm1lgfgOgQfcB2jb8M+n9S2HV4q33azulIiFJTm50ePTwhIvuzdNc1kNSo/z6fTKru2BiCrt+4lAAGyDMEHXCXWt+FwGsD/jOitUVVlKT1ssrZSYtlIadfL0eNZvM0itQQej72+Ne3VKE8lDgM5iuADrhHv27DU8m34H3X1uuWi1BILs7JSwmib5dJdUmGPpD/CjWXFRoFDPO1Xo9LZqgHgPgQfcA2rcjPSPT3VEZselN6cHD2e4mqHG1cHYgUO8bRdjfJM4jAASm3hHlbmZoQrJcr8kVsrZf5C93xbrvZFBx4jl6UVeLitrDjRqlZ7PrUES21Xo2ixDngHKx9wDatzM1xbKbF7hfTiWdHjaeR2uHV1IFHg0Fas1aisThwGEIHgA65hR25GotNTbWeU23HSFOnU36b1cW4tK04lIIhVyZS1icMAohB8wDWyKjcjUwc/lZ40CKIyrGRx6+pAsgHB9ItO1ndHVBr+N87KxGEAhsj5gKuYnZuRTOtu21X7LAk8JPeuDoQDh1hhYzjHI1bgIdFiHfASVj7gOmblZrix4sNwm+XyJqnDUaZ8vFtXB8xa1aLFOuANvlAo5IKvgkc0NjbK7/eroaFBJSUlTk8HWSpWWWf40WZ7hcuzX5H2b4karhn2ienJr+HfXZJhAJJOozazmBUQurGHCZDrUnl+E3zAcwLBkM6c9VLMxMvwt/8VN400/YFl+FBcEL27WfOVv2vqsgLLVmWMHvJ5viOdYs2+XyoIHABvIvhATqvZskcTHl2Z8Lr5V59uasVH+wf+/3R/RtMr/hB93cCPbFmVCT/kX9hQrz+9vi3qfcdWgQB4UirPb8sTTu+55x75fD5NmTLF6lsBkpyp+Gjf2GvbwDHRgUe/nyswPhi3D4fU0ofDjMTY8AFsz9XVG75v9v0AIFmWBh9r1qzRww8/rIEDB1p5GyCC3RUfbRt79Svcom0Dx0RdM/z9FxUYcIftXTrpCgrAjSwLPvbv368rrrhCjz76qLp06WLVbYAoyZZ1mlXxEX7Abxs4Rn8/8YaI974IFqj3+iWtD3i7V2Xc2vcDQG6zLPiYPHmyLrroIo0aNSrudc3NzWpsbIz4B2TC7n4QHzc2Gq529F6/WCfXPdX6OpxgmQyzVmXc1PfDlT1XADjCkj4fCxYs0Lp167RmzZqE186cOVMzZsywYhrIYbb1g6j26ZsGw73XL4kaC1d22NmHwy19P1zZcwWAY0yvdtm+fbuGDh2qF154oTXX45xzztHgwYM1Z86cqOubm5vV3Nzc+rqxsVE9e/ak2gWmsLSs06Bh2Fc3/D99/J/Ibcb2pb2x+nBYVX1i9/1i3d81PVcAWMLRUttnnnlGl1xyifLz81vHAoGAfD6f8vLy1NzcHPFee5TawvVeu0za/mTUcOXh1Y5kHvB2rwQ4tfLgZM8VAPZyNPjYt2+f3n///Yixq666Sn379tVNN92k/v37x/15gg+4mlF79LOeknpemvID3u5mW04093Kq5woA+6Xy/DY956O4uDgqwOjcubO6deuWMPAAXOvDRdKrF0ePtzkMLtUzafLzfLY+cO2+n0S1DQBjHCwHJGK02nHsxdLZC6OGnXjAu5mbqm0AuIctwccrr7xix20Acx34RHr6y9HjE83ZqQwEQ1q5ZY9q3vtEUkvQcvrx3dLeCnHjmSluqbYB4C6sfABGjFY7JNMCj6V1OzXt6X/ps88PtY79/uV3dfRRHXXPpQNSTgJ1aylruOfKtfPWySfjZFwze64AyA6Wn+0CZJVQyDjw+HazqYHHNfPWRQQeYZ99fkjXzFunpXU7U/q8tufKhNU3HNC1KX6WFcI9V8r8kVsrZf5CymyBHMWptkDY0qHS3rXR4yYFHVLL1siIe5apvrE57nXlSZafZlMpqxu3hQCYx9FqFyArGa12jNkklZxg6m1Wb92bMPCQjhz2lih5NZWD45xOhCUZF0AY2y7IbRvvMw48JoZMDzyk1EpKk7mWUlYA2YiVD+Quo6Bj2B+lPt+z7JaplJQmcy2lrACyEcEHcs9n/5L+MTB63MTcjlhOq+yqspKCpHI+kik/pZQVQDZi2wW5pdoXHXiUjrQl8JBa8h5u/2a/hNclW34aLmWVjpSuhlHKCsCtCD6QGwLNsXM7zl1m61RG9y/XQ1cO0dFHdYx6r8tRHfVQiuWnlLICyDaU2sL7nh8u7TE43Mym1Y5YcqHDKYDc4eiptpki+ICpjFY7LmuUOhbbPxcA8DD6fACbH5bWXBM97vBqRzpY0QDgNQQfWYyHUgyGDcM2SiUn2j+XDLn1zBYAyATBR5bioWRg7zpp6anR41m42iEdObOl/ezDZ7aQTAogW1HtkoXcfpCYI6p90YHHmU9kbeARCIY0Y/EGw94d4bEZizcoEMzO3w9AbiP4yDI8lNo5+FnsEtpel9k+HbOkcmYLAGQbgo8sw0OpjUW9pSe7RI595RrLVzsCwZBqtuzRotodqtmyx5JAjzNbAHgZOR9ZhoeSpFBImm8QN08ISD5r42m7cm04swWAl7HykWWceijZ8W0/KWunRgceBd1bVjtsCDzsyrUJn9kSq3bJp+TPfwEAt2HlI8s4cZCYayprjHI7Lv1YKuxu+a0T5dr41JJrc15VmSnlzuEzW66dt04+KeK+nNkCINux8pFl7D5IzBWVNR8uip1UakPgITmTa8OZLQC8ipWPLBR+KLVfjSgzeTXC7m/7hoyCjvNXSt2HmXaLZJq1OZVrM7p/uc6rKqOZHABPIfjIUnY8lFL5tj+8TzfT7itJ2veutPiE6HGTK1mS3VJyMgE0P89n/t8XABxE8JHFrH4oOVZZs/Sr0t43I8eGzJb63mjubVLoIOpErg0AeBU5H4jJrG/7SVfKBA62bLO0DzwmhkwPPFJt1mZ3rg0AeBkrH4jJjG/7SVfK1N4ibZgZ+cMnXi8N/V1mv0QM6Wwp2ZVrAwBeR/DhUWaceJtpuWfS2xpGSaXjD0p5HU37XdpLd0uJBFAAyBzBhweZ2Zcj3W/7yVTKvPLiIxq9/vbIN7ufIZ3/uiW/S1uZbCmRAAoAmfGFQiFXnUDW2Ngov9+vhoYGlZSUOD2drBNrtSH8vTzd/hCprj7UbNmjCY+ujPn+toFjogcvqZeKSltfWvW7SC2/z5mzXkq4pbTippGsagBAElJ5fpNw6iGJVhtCSv/E2/C3/XGDj9HwPt0SPpBjbWt8peAD48BjYigi8LD69F4SSAHAOQQfHpIoiVJqSaL8/UubLZ+L0XbFtoFj9OJJP4wY+1e/f0b17ggEQ3rs9a2WdxSlgygAOIOcDw9JNonyty9u1kllxZY+XNtWyny5wx6trpoUdc3w91/UivEjI8aMcjziybTHCAmkAGA/gg8PSaW7ptVt0cPbGqPXV0S99/DuS3VP/fc098rIbY1YOR7xmNFRlARSALAXwYeHhFcbklk1sKwtelgwYBh49F6/WOX+Is29MrJaJV6OhxE6igJA9iL48JDwasM189Yldb3pbdHDnh8u7YmudFlU9aHmDzPe1kgmXyWMhFAAyG4EHx4zun+5bhx1on774qaE11pxCJphw7BvfSIVdNO4OD+WSiBER1EAyG5Uu3jQdSO/orKSgrjX5PmkT5sOmnfTd35rHHhMDEkFibd2kg2Epl90slbcNJLAAwCyGMGHB+Xn+XT7N/tF9a9oKxiSJlev09K6nZnfsNonrZsaOTZqeVQJbTzhfJVYc/appbPpd0dUstUCAFmO4MOjRvcv1wMTT1Gi53Qmjbq0d23s1Y4eZyf88ban3a7eulfTL6LpFwDkAnI+PKxL5wLFiyuMTm5NmlHQMeB2acBtSf14rDNbvn92pZ59eyenxgKAhxF8eFi6J7fGdahR+ps/ejyFLZZ4p90+8upWPTDxFHXpXKD6xgPau79ZXTt3kr+okwLBkOkrH1acmAsAiI/gw8MyObnVkNFqR/EJ0tjElTVhyZx2e8ff/63pF1Xp3qXvmH6abVtWnZibDoIgALmEU209zLSTW0Mhab5BetC3m6X8TinNKdFpt/GYcZptmJUn5qYzF7cEQQCQLk61hSSTTm6t+a5x4DExlHLgIWXW2MyM02wl60/MTUU4CGrfYK2+4YCunWdSNRIAuAzBh8e0rSCp2bJH51WVpX9ya7VP2vqXyLFx21LK72gv08ZmZpxmm6ibqhn3SIabgiAAsBM5Hx4Sb/l+xU0jk88p2PywtOaa6PEMgo6wtqfdZvJpmaygWJKIm4ZUgiAOvgPgJQQfHhGvguTaeeuSz2EwSiod8bh03OWmzDO8FXTtvHXySRHzbf86nkxWUExPxE2TW4IgALAb2y4eYMryfcM7hoFHzbBPFOh5mSnzDBvdvzzmVtCDE09JqtNpJqfZJttN1eoTc90SBAGA3Vj58ICMl+8Ngo61TX31rS2/ltavtKTyYnT/cp1XVWa4FZSX54u5MiJl3uk0mdWXC/u3zM3KktdEW1DhaiSrgyAAsBsrHx6Q9vJ94KBh4NF7/eKWwOMwqyov8vN8Gt6nm8YNPkbD+3RrfcjHWxkxqwQ21j18h/8cf3p9myY8ulJnznrJsooTU6qRACAL0efDA5LtnTH/6tOPrHwY5XZI6r1+ieF40j1BTGRH463wPV7cUK8/vr4t6n07+n7Q5wOAF6Ty/GbbxUZWPUxTXr43CDzWnLpZl/15Y8x7OFF5EV4Zsfoep1V21dQnag3fD3ddnbF4g86rKrMk8Iq3BQUAXkTwYRMrv90mymGQDi/fr/mBtOXR6A+YGNJHtTuSupcXKy/cUPJqR6AFAG5BzocN7OhimTBPYn1FdOBx7sutvTtyufKCklcAsJfpKx8zZ87U008/rXfeeUdFRUU644wzNGvWLJ100klm3yorJHOQmllL+obL90dvV/7zFdEXt2sYlsuVF7kceAGAE0xf+Vi+fLkmT56slStX6oUXXtChQ4d0/vnnq6mpyexbZQU7W3m3zykZvqq78p8/JfKiU35l2Kk0lysv3NL3AwByhekrH0uXLo14/dhjj6lHjx5au3atzj77bLNv53p2Lem3zSkp9B3QOwP+K/qiBO3Rw1s37XNTyjxeeZF0zowHAy8AcILlCacNDQ2SpK5dc/Nbox1L+m1bq999zP2a2O35iPc/7DJex144P6nPytXKi1wNvADACZYGH8FgUFOmTNGIESPUv39/w2uam5vV3Nzc+rqxsdHKKdnO6lyKIzklIW0bODbq/ePXP6tS/1FacUEo6QAiVysvcjXwAgC7WVrtMnnyZNXV1WnBggUxr5k5c6b8fn/rv549e1o5Jcu0P8o+fI6K1bkUq7fu1Zl5i6MCj+cazlDv9UsUVJ4tx8N7RayuqwAA81jW4fS6667TokWL9Oqrr6qysjLmdUYrHz179syqDqfJ9PCwrM+HQcOwAXWPa1+wc8TY78YP1rjBx6R/HwAA4nC0w2koFNKPfvQjLVy4UK+88krcwEOSCgoKVFBQYPY0bJPsUfamL+nv3yY9G/m3PRDspL51TxteTpkoAMAtTA8+Jk+erOrqai1atEjFxcWqr6+XJPn9fhUVFZl9O0el2sPDtFyKN66Utv01Yui8jQ9qc3OvqEu93J8DAJCdTM/5mDt3rhoaGnTOOeeovLy89d/jjz9u9q0cZ2cPD0nSocaWbZa2gUe307R04Ed6t7lXzvXnAABkJ0u2XXKFrW25N/xKqv1Z5NjodVLXUzRaokwUAJA1OFguA7a05Q4GpAUG/5naNQyjTBQAkC0IPjJg+XkoHz4rvToucuysp6Welxhe7vb+HO3bvxMcAUBuIvjIgKVtuZ/6stT8SeTY+ENSXnb+J7Os1BgAkHUsbTKWCxIeZZ/qg/XgZ9KK8ZGBx6C7WrZZsjjwuHbeuqjk3HA58tK6nQ7NDADghOx8mrmMafkWu16WaiZJn29vHVp96ns69YTeyjd5znZJtRwZAOB9BB8mySjfInBAevvn0juzJUnbDx2jG96fonWfnyyt36By/3uObk9kkquRSjmym/NVAADmybngw3VJj5/WtjQNa/g/SdJf94zWXTv/R58HjzRka98t1U6Z5mrYWo4MAMgKORV8uCrpMRiQ3vm1tH66FDykUGGpfvLBj/RU/eCoS53anki2dXw8tpQjAwCySs4knLoq6XH/VmnZOVLtNCl4SDr2Er1Ztdww8AgzvVtqAolyNaSWYCh8em8s4XLkWOGSTy0BIO3fASB35ETwYdaDNGOhkPTaZdKzx0sfr5A6FEun/1k66yl9dOBLSX2EXdsTZrWOD5cjS6L9OwBAUo4EH7afwWJk1yvS/Dxp+5Mtr798lvSNt6Xjvyv5fK7bnjAzV8P0cmQAQFbLiZwPx5MeFx4jffFR5Ni5L0t5RwpoLe+WmiKzgyHavwMAwnJi5cOxVYXGzS2n0LYNPPr++HDDsMjOHW7bnrAiVyNcjjxu8DEa3qcbgQcA5KicCD4cSXp89VJpyYmRY5fukob8OuaPuGl7wm3BEADAO3yhUMjiLMvUNDY2yu/3q6GhQSUlJaZ9brjaRTI+g8W0h/vBz6Qnu0SO9ThHGvVy0h/hpl4kripPBgC4VirP75wJPiQbHqTLRra0SG/rG+ulowdk/tkOclMwBABwp1Se3zmRcBrWNumxvuEL7W06qK5fKpC/qJMCwVD6D9RAs/R4u3wRXwdpwqHMJ+0CGbWOBwCgnZwKPqSWB2nDFwd17/MbzVkBWfcT6Z3fRI6d+YTU6zITZgsAgPfkXPBhRsvwVtUGKyXjD0l5OfdnBQAgaTlR7RJmWqfTD/4WHXj0OOdwCS2BBwAA8eTUk9KU492NVjsu3S0VftmcSQIA4HE5FXxk1On0szrpHwZVKxNdVSxkiGoVAICb5FTwkXan0/kdpFAgcmz0m1LXUyOG3PiQp08HAMBtcir4SPn8lIOfSk8adD01WO1w40Pe1ORaAABMklMJpym1DF8+LjrwOP2xmIHHtfPWReWThB/yS+t2mvMLpMC05FoAAEyWU8GHlMT5Kf1KW5JKdzwb+YMTgtLxk6I+z60P+VSSawEAsFNObbuExTzefeNsaf5PIi/+yg+k0x6K+VmmVNBYIKPkWgAALJSTwYdk0DLcqIT28iapw1FxP8etD/m0k2sBALBYzm27RKlfFh14dK5sye1IEHhI7n3Ih5NrY9Xa+NSSENuaXAsAgE1yO/h49gTppVGRY998Txr3XtIf4daHfErJtQAA2Cg3g4/mPS2rHfvfjRyfGJK+VJnSR7n5IZ8wuZYyWwCAA3yhUMhVtZaNjY3y+/1qaGhQSUmJ+TfYtkB6Y0Lk2MhlUtnIjD7WjX0+wtzY/AwA4C2pPL9zK/j45wjpkzeOvO5QLF3eaNrH85AHAOSqVJ7fuVPtcqgxMvAY9arU4yxTbxFVQQMAAKLkTvDRsUT66lyp6QNp0J2SLzfTXQAAcFruBB+SdMI1Ts8AAICcx9d/AABgK4IPAABgK4IPAABgK4IPAABgK4IPAABgq5ypdqEBGAAA7pATwYebW58DAJBrPL/tsrRup66dty4i8JCk+oYDunbeOi2t2+nQzAAAyE2eDj4CwZBmLN4go8NrwmMzFm9QIOiq420AAPA0Twcfq7fujVrxaCskaWfDAa3eute+SQEAkOM8HXzs3hc78EjnOgAAkDlPBx89igtNvQ4AAGTO08HHaZVdVe4vVKyCWp9aql5Oq+xq57QAAMhpng4+8vN8um1slSRFBSDh17eNraLfBwAANvJ08CFJo/uXa+6VQ1Tmj9xaKfMXau6VQ+jzAQCAzXKiydjo/uU6r6qMDqcAALhATgQfUssWzPA+3ZyeBgAAOc/z2y4AAMBdCD4AAICtLAs+HnjgAfXu3VuFhYUaNmyYVq9ebdWtAABAFrEk+Hj88cc1depU3XbbbVq3bp0GDRqkCy64QLt377bidgAAIItYEnzMnj1bV199ta666ipVVVXpoYce0lFHHaU//elPVtwOAABkEdODj4MHD2rt2rUaNWrUkZvk5WnUqFGqqamJur65uVmNjY0R/wAAgHeZHnx88sknCgQCKi0tjRgvLS1VfX191PUzZ86U3+9v/dezZ0+zpwQAAFzE8WqXm2++WQ0NDa3/tm/f7vSUAACAhUxvMta9e3fl5+dr165dEeO7du1SWVlZ1PUFBQUqKCgwexoAAMClTA8+OnXqpFNPPVXLli3TxRdfLEkKBoNatmyZrrvuuoQ/HwqFJIncDwAAskj4uR1+jsdjSXv1qVOnatKkSRo6dKhOO+00zZkzR01NTbrqqqsS/uy+ffskidwPAACy0L59++T3++NeY0nw8e1vf1sff/yxbr31VtXX12vw4MFaunRpVBKqkYqKCm3fvl3FxcXy+Tj4LR2NjY3q2bOntm/frpKSEqenkzP4uzuHv70z+Ls7x41/+1AopH379qmioiLhtb5QMusjyCqNjY3y+/1qaGhwzf8ocwF/d+fwt3cGf3fnZPvf3vFqFwAAkFsIPgAAgK0IPjyooKBAt912GyXMNuPv7hz+9s7g7+6cbP/bk/MBAABsxcoHAACwFcEHAACwFcEHAACwFcEHAACwFcGHh8ycOVNf/epXVVxcrB49eujiiy/Wxo0bnZ5Wzrnnnnvk8/k0ZcoUp6fieTt27NCVV16pbt26qaioSAMGDNCbb77p9LQ8LxAIaPr06aqsrFRRUZH69OmjO+64I6kzPZCaV199VWPHjlVFRYV8Pp+eeeaZiPdDoZBuvfVWlZeXq6ioSKNGjdLmzZudmWwKCD48ZPny5Zo8ebJWrlypF154QYcOHdL555+vpqYmp6eWM9asWaOHH35YAwcOdHoqnvfpp59qxIgR6tixo5577jlt2LBBv/nNb9SlSxenp+Z5s2bN0ty5c/X73/9e//73vzVr1izde++9uv/++52emuc0NTVp0KBBeuCBBwzfv/fee3XffffpoYce0qpVq9S5c2ddcMEFOnDggM0zTQ2lth728ccfq0ePHlq+fLnOPvtsp6fjefv379eQIUP04IMP6s4779TgwYM1Z84cp6flWdOmTdPrr7+u1157zemp5JwxY8aotLRUf/zjH1vHvvWtb6moqEjz5s1zcGbe5vP5tHDhwtYT40OhkCoqKvTjH/9YP/nJTyRJDQ0NKi0t1WOPPabx48c7ONv4WPnwsIaGBklS165dHZ5Jbpg8ebIuuugijRo1yump5IRnn31WQ4cO1WWXXaYePXrolFNO0aOPPur0tHLCGWecoWXLlmnTpk2SpLffflsrVqzQhRde6PDMcsvWrVtVX18f8f85fr9fw4YNU01NjYMzS8ySU23hvGAwqClTpmjEiBHq37+/09PxvAULFmjdunVas2aN01PJGe+9957mzp2rqVOn6pZbbtGaNWt0/fXXq1OnTpo0aZLT0/O0adOmqbGxUX379lV+fr4CgYDuuusuXXHFFU5PLafU19dLUtSJ8aWlpa3vuRXBh0dNnjxZdXV1WrFihdNT8bzt27frhhtu0AsvvKDCwkKnp5MzgsGghg4dqrvvvluSdMopp6iurk4PPfQQwYfFnnjiCf31r39VdXW1+vXrp9raWk2ZMkUVFRX87ZEUtl086LrrrtOSJUv08ssv69hjj3V6Op63du1a7d69W0OGDFGHDh3UoUMHLV++XPfdd586dOigQCDg9BQ9qby8XFVVVRFjJ598sj744AOHZpQ7fvrTn2ratGkaP368BgwYoO985zu68cYbNXPmTKenllPKysokSbt27YoY37VrV+t7bkXw4SGhUEjXXXedFi5cqJdeekmVlZVOTyknnHvuufrXv/6l2tra1n9Dhw7VFVdcodraWuXn5zs9RU8aMWJEVCn5pk2bdNxxxzk0o9zx+eefKy8v8vGRn5+vYDDo0IxyU2VlpcrKyrRs2bLWscbGRq1atUrDhw93cGaJse3iIZMnT1Z1dbUWLVqk4uLi1j0/v9+voqIih2fnXcXFxVF5NZ07d1a3bt3It7HQjTfeqDPOOEN33323Lr/8cq1evVqPPPKIHnnkEaen5nljx47VXXfdpV69eqlfv3566623NHv2bH3ve99zemqes3//fr377rutr7du3ara2lp17dpVvXr10pQpU3TnnXfqhBNOUGVlpaZPn66KiorWihjXCsEzJBn++/Of/+z01HLO1772tdANN9zg9DQ8b/HixaH+/fuHCgoKQn379g098sgjTk8pJzQ2NoZuuOGGUK9evUKFhYWh448/PvTzn/881Nzc7PTUPOfll182/P/1SZMmhUKhUCgYDIamT58eKi0tDRUUFITOPffc0MaNG52ddBLo8wEAAGxFzgcAALAVwQcAALAVwQcAALAVwQcAALAVwQcAALAVwQcAALAVwQcAALAVwQcAALAVwQcAALAVwQcAALAVwQcAALAVwQcAALDV/wc04PKToZF/2AAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"ax.plot(x_data,p(x_data), '-', color= \"orange\")\n",
"fig"
]
},
{
"cell_type": "markdown",
"id": "83da75ff-1807-40fd-b161-f17f07b152cd",
"metadata": {},
"source": [
"\n",
"---- \n",
"\n",
"Enjoy!"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.10.6"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment