You cannot select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.

150 lines
34 KiB
Plaintext

{
"cells": [
{
"cell_type": "code",
"execution_count": 12,
"id": "0f872066",
"metadata": {},
"outputs": [],
"source": [
"# data for the OLS\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"\n",
"theta0 = 2.6486\n",
"theta1 = 4.98\n",
"m = 50\n",
"disp = 4\n",
"x = np.linspace(0,5,m)\n",
"y = theta0+theta1*x+disp*np.random.rand(m)"
]
},
{
"cell_type": "code",
"execution_count": 13,
"id": "c4906c75",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAh8AAAGdCAYAAACyzRGfAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAIElJREFUeJzt3QmMVtX5P/AzoFCqAgVZA1qUzSXalrrFJYjolDYKVVNrTNRm0qEtmiratCSNS2KC1aTV5icUEiIxKWo1RauJGkWWmkKrNMQlkYCFAGHRmrA2oGH459x/ZzqDA8478773vcvnk9xc3oWZm+s475dznvOchiNHjhwJAAAp6ZXWNwIAED4AgNQZ+QAAUiV8AACpEj4AgFQJHwBAqoQPACBVwgcAkKoTQsa0tLSE7du3h1NOOSU0NDTU+3IAgC6IPUv37dsXRo4cGXr16pWv8BGDx+jRo+t9GQBAN2zdujWMGjUqX+Ejjni0Xnz//v3rfTkAQBfs3bs3GTxo/RzPVfhonWqJwUP4AIB86UrJhIJTACBVwgcAkCrhAwBIlfABAKRK+AAAUiV8AACpEj4AgFQJHwBAqoQPACBVwgcAkCrhAwAybtu2bWH58uXJuQiEDwDIsEWLFoXTTz89TJkyJTnHx3nXcOTIkSMhY7viDRgwIOzZs8fGcgCU2rZt25LA0dLS0vZc7969w+bNmzvdtj6+f8OGDWHcuHFfuq19PT+/jXwAQEZt2LChQ/CIDh8+HDZu3JjrEZKKwsf8+fPDeeed17bd/SWXXBJeeeWVttcPHjwYZs2aFQYPHhxOPvnkcMMNN4Rdu3bV4roBoPDGjRsXevXq+FEdRz7Gjh37hRGP5ubmtqASzzNnzsxsjUhF4SMO4Tz88MNh7dq14Z133knS1fTp08MHH3yQvH733XeHl156KTz33HNh5cqVYfv27eH666+v1bUDQKGNGjUqLFy4MAkcUTwvWLDgC1MqlYyQFKLmY9CgQeHRRx8NN954YxgyZEhYsmRJ8ufoww8/DGeddVZYvXp1uPjii7v09dR8AEBHcQQjBok44nGsWo9KakNyW/MRE9UzzzwTDhw4kEy/xNGQzz//PEydOrXtPRMnTgynnXZaEj6O5dChQ8kFtz8AgP+JAWLy5MnHDBJdHSHJihMq/QvvvfdeEjZifUes61i6dGk4++yzw7p160KfPn3CwIEDO7x/2LBhYefOncf8enPnzg0PPvhg964eAEg0NTWFxsbG446Q5DZ8TJgwIQkacVjl+eefD7fddltS39Fdc+bMCbNnz257HEc+Ro8e3e2vBwBlNWrUqEyHjm6Hjzi60VplO2nSpPD222+Hxx9/PNx0003hs88+C7t37+4w+hFXuwwfPvyYX69v377JAQCUQ4/7fMTilli3EYPIiSeeGJYtW9b22vr168OWLVuSaRoAgIpHPuIUybRp05Ii0n379iUrW1asWBFee+21pMI1zjfFKZS4AiZWut55551J8OjqShcAKJttdexKmovw8fHHH4dbb7017NixIwkbseFYDB5XX3118vrvfve7pBlKbC4WR0Ni4cu8efNqde0AkGuLFi1qaw4WPz/jipX4D/mis7cLANTBtgz05qgme7sAQMZtyFlX0mqysRwAZHjfliISPgCgDkbVqStpnO5Zvnx5XTedU/MBABnetyUvBa6V1HwIHwBQAttqXOCq4BQAyGyBq5oPACiBcRkqcBU+AKAERtWpwLUzaj4AoES21ajAtZKaj4p3tQUA8mvUqFF176Bq2gUASJXwAQCkSvgAAFIlfAAAqRI+AIBUCR8AQKqEDwAgVcIHAJAq4QMAathNdPny5cmZ/xE+AKAGFi1alGxhP2XKlOQcH/P/2dsFAKosjnTEwNF+C/u4kdvmzZvr3tq8VirZ28XIBwBU2YYNGzoEj+jw4cPJhm4IHwBQdePGjQu9enX8930c+Yg7ySJ8AEDVxamVhQsXJoEjiucFCxYUdsqlUmo+AKCGtR9xqiWOeBQ9eOytoObjhNSuCgBKJgaOooeO7lBwCkBp6LuRDcIHAKWg70Z2qPkAoPDK2Hcjbfp8AEA7+m5ki2kXAApP341sET4AKDx9N7JFzQcApVGmvhtp0+cDADqh70Y2mHYBAFIlfAAAqRI+AIBUCR8AQKqEDwAgVcIHAIVg07j8ED4AyD2bxuWLJmMA5JpN47LBxnIAlIZN4/LHtAsAuWbTuPwRPgDItXptGqfAtfvUfABQCGluGhcLXJubm0NLS0vo1atXEn6amppCme3duzcMGDAg7NmzJ/Tv3/+47xU+AKACClw7p+AUAGpEgWvPqfkAoGaKWBehwLXnhA8AaqKojb/qVeBaJGo+AKi6vNdFxOuP0ytxlONY15tmgWseqPkAoK7yXBfR1RGbGDgmT54seNR62mXu3LnhggsuCKecckoYOnRomDFjRli/fn2H98T/EA0NDR2On/zkJ925NgByKq91EXE0o3UJbRTPM2fOLFTNSu7Cx8qVK8OsWbPCmjVrwuuvvx4+//zzcM0114QDBw50eN+Pf/zjsGPHjrbjkUceqfZ1A5Bhea2LyPOITZ6cUMmbX3311Q6PFy9enIyArF27NlxxxRVtz3/1q18Nw4cPr95VApA7selWY2NjruoiWkdsjq5VyfqITalWu8QuZtGgQYM6PP/HP/4xnHrqqeHcc88Nc+bMCf/5z3+O+TUOHTqUFKm0PwAohrzVReR1xKY0q11iKrzuuuvC7t27w1tvvdX2fPyPFgt0Ro4cGd59993wy1/+Mlx44YXhz3/+c6df54EHHggPPvjgF57vSntWAKgFK1ky2l79pz/9aXjllVeS4HG8RPjmm2+Gq666Khl2O/PMMzsd+YhH+4sfPXq08AEABQ0fFdV8tLrjjjvCyy+/HFatWvWlQ1EXXXRRcj5W+Ojbt29yAFBOXempQYlrPuIgSQweS5cuTUY0xowZ86V/Z926dcl5xIgR3b9KAAqpqF1QCdWbdvnZz34WlixZEl588cUwYcKEtufjMEu/fv3CRx99lLz+3e9+NwwePDip+bj77ruTJBuX6VZ72AaA/Mp7F1RS6nA6f/785IvGyuU4ktF6PPvss8nrffr0CW+88UbS+2PixInhnnvuCTfccEN46aWXKvk2AJSAnhrlVVHNx5cNksRC0a6OcABQbnpqlJddbQGoCz01ysuutgDUlZ4axVDzpbYAUM0REAWm5WLaBQBIlfABAKRK+ADgCzUYy5cvT85ZkcVrovuEDwAy3XE0i9dEz1jtAkBmO45m8ZpIucMpAMWVxY6jWbwmek74AKBDx9H24ijD2LFj63aHsnhN9JzwAUBmO45m8ZroOTUfACUSayjiVEYcUTjWB3gWO45m8ZroSIdTAL4grhJpbm5OaijiVEYcUWhqaspFx9EsXhPdZ+QDoASsGqHWrHYBoAOrRsgSBacAJWDVCFkifACUgFUjZImaD4ASsWqEWrHaBYBOWTVCFph2AQBSJXwAUDFb3NMTwgcAFbHFPT2l4BSALtOsjGPRZAyAmtCsjGow7QJAl2lWRjUIHwB0mWZlVIOaDwAqplkZR9NkDICa0qyMnjDtAgCkSvgAqCPNuigj4QOgIM26BBnyQvgAqIMYFJqbm0NLS0vyOJ5nzpyZPN8duo6SJ8IHQM6bdVU7yECtCR8AdZjiqGazLl1HyRvhA6BC1ZjiqGazLl1HyRtNxgDquLFatZp1xQAUp1ri1E1rkGlqaur214NKaTIGUCPHm+LoTnioVrOuGDQaGxurEmSg1k6o+XcAKJDWKY6jRz66U6tRbbqOkhdqPgAqYGM16Dk1HwDdYGM16EjNB0CNmeKA7jPtAgCkSvgAAFIlfAAAqRI+AIBUCR8AJdhLBrJE+AAowV4ykCX6fACUaC8ZyEKfDyMfABmeKjneXjKQV8IHQIanSlr3kmkvK3vJQHcJHwBVFkc6mpub20Ys4jlud9+dERB7yVBEdrUFqLLjTZV0p06jqakpNDY2Jn8/jnio9aBUIx9z584NF1xwQTjllFPC0KFDw4wZM8L69es7vOfgwYNh1qxZYfDgweHkk08ON9xwQ9i1a1e1rxsgs2oxVRIDx+TJkwUPyhc+Vq5cmQSLNWvWhNdffz18/vnn4ZprrgkHDhxoe8/dd98dXnrppfDcc88l79++fXu4/vrra3HtAJlkqgRquNT2k08+SUZAYsi44oorkuU1Q4YMCUuWLAk33nhj8p4PP/wwnHXWWWH16tXh4osvrupSHYAsizUepkooi70VfH73qOYjfoNo0KBByXnt2rXJaMjUqVPb3jNx4sRw2mmnHTN8HDp0KDnaXzxAUUZA1GdAFVe7xGKqu+66K1x66aXh3HPPTZ7buXNn6NOnTxg4cGCH9w4bNix57Vh1JDEptR6jR4/u7iUBAEUOH7H24/333w/PPPNMjy5gzpw5yQhK67F169YefT0AINu6Ne1yxx13hJdffjmsWrWqw5Di8OHDw2effRZ2797dYfQjrnaJr3Wmb9++yQEAlENFIx+xNjUGj6VLl4Y333wzjBkzpsPrkyZNCieeeGJYtmxZ23NxKe6WLVvCJZdcUr2rBqghO8hChkY+4lRLXMny4osvJr0+Wus4Yq1Gv379knNshjN79uykCDVWu955551J8OjKSheAY4WB2Lgr9s+odQFnbIPe2p009upYuHBh8nsNqNNS24aGhk6ff/LJJ8Ptt9/e1mTsnnvuCU8//XSyiiV25Zs3b94xp12OZqkt0J0wUI2AYgdZ6L5KPr971OejFoQPoNIwUK3RirgDbdwIrrPnY3dRoDqf3zaWAzKrK9vJV3MTNzvIQjqEDyCzuhIGuhJQukpbdEiH8AFkVlfCQLVHK+J0TZzWiVMt8azYFKpPzQeQ+z1SYs1HnGqJIx6tAUVogHQpOAVKxyZuUJKN5QCywiZukB9qPgCAVAkfAECqhA8AIFXCBwCQKuEDAEiV8AEApEr4AABSJXwAAKkSPgCAVAkfAECqhA8AQPgAAIrLyAcAkCrhAwBIlfABAKRK+AAAUiV8AACpEj4AgFQJHwBAqoQPACBVwgfQLdu2bQvLly9PzgCVED6Aii1atCicfvrpYcqUKck5Ps4LoQnqT/gAKv7wbm5uDi0tLcnjeJ45c2YuRkDyHJqgSIQPoCIbNmxoCx6tDh8+HDZu3JjpO5nn0ARFI3wAFRk3blzo1avjr47evXuHsWPHZvpO5jU0QREJH0BFRo0aFRYuXJgEjiieFyxYkDyfZXkNTVBEwgdQsaamprB58+ZktUs8x8dZLwDNa2iCImo4cuTIkZAhe/fuDQMGDAh79uwJ/fv3r/flADUUCz5b6zDiqEQMBz0JMl0RQ06caokjHoIH1OfzW/gAavpBH2st4pTH0R/08bW44qR9HUYcjYgjKUIBFDt8mHYB6rKsVQEolJfwAdRlWasCUCgv4QOouq6MaigAhfI6od4XAGSjBqOaWkc1jq7nOHpZaywubWxsVAAKJWPkA0oizdbilYxqxOcmT56syBRKxGoXKIF6rSyxrBXKY28Fq11Mu0DJazBqGT7i17ZsFjiaaRcoAStLgCwRPiAHetqC3MoSIEvUfECJWpCrwQBqRXt1KAgtyIG80F4dCkILcqCI1HxAhikUBYpI+IAMq0ehaE+LWwG+jIJTyIG0CkWrWdwKlMveCpqMCR9AQnEr0BMKTqGEejpdorgVyGzNx6pVq8K1114bRo4cGRoaGsILL7zQ4fXbb789eb798Z3vfKea1wzUYNM4xa1AZsPHgQMHwvnnnx+eeOKJY74nho0dO3a0HU8//XRPrxM4hjjS0VqnEcXzzJkzKx4B0QUVSEvFG8tNmzYtOY6nb9++Yfjw4T25LqAOm8bF4tLGxsZUiluB8qrJrrYrVqwIQ4cODV/72teSYeCHHnooDB48uNP3Hjp0KDnaF6wAlU+XtA8gcUluDA/dYSdaIHd9PuKUy1NPPRWWLVsWfvOb34SVK1cmIyXxX2KdmTt3brI0p/UYPXp0tS8JCs10CZA3PVpqG4tJly5dGmbMmHHM9/zrX/8KZ555ZnjjjTfCVVdd1aWRjxhAurJOGPgfm8YBeVlqW5Npl/bOOOOMcOqppyZzyJ2Fj1gfEg/I2wd9rLWIUx5ZqYswXQLkRa80fkl/+umnYcSIEbX+VpCpZa3alANUKXzs378/rFu3LjmiTZs2JX/esmVL8tovfvGLsGbNmrB58+ak7mP69OlJ4VusoIeyLGutRt8NgKKqOHy888474Zvf/GZyRLNnz07+fN999yUV9u+++2647rrrwvjx45Nle5MmTQp//etfTa1QCF3pAlqtvhsARVVxzcfkyZPD8WpUX3vttZ5eE+R6WWs1+24AFFHNaz6gbMtatSkHOD7hAyoUpxNjTVPcxC2ej95yXt8NgBr2+aj3OmHIMn03gDLZm6U+H1BW+m4AdM60CwCQKuED6kwzMqBshA+oI83IgDJScAp1HPGI3U+P7hkSV9DoBwIUueDUyAdkuFsqQBEJH1AnmpEBZSV8QJ1oRgaUlZoPqDPNyIAi0GQMckQzMqBsTLtQGPplAOSD8EEh6JcBkB9qPsg9/TIA6k+fD0pFvwyAfDHtQu7plwGQL8IHuadfBkC+qPmgMPTLAKgffT4oJf0yAPLBtAsAkCrhAwBIlfABAKRK+AAAUiV8AACpEj4AgFQJHwBAqoQPACBVwgcAkCrhAwBIlfABAKRK+AAAUiV8AACpEj4AgFQJHwBAqoQPACBVwgcAkCrhAwBIlfABAKRK+ICjbNu2LSxfvjw5A1B9wge0s2jRonD66aeHKVOmJOf4GIDqajhy5MiRkCF79+4NAwYMCHv27An9+/ev9+VQInGkIwaOlpaWtud69+4dNm/eHEaNGlXXawPIuko+v418wH9t2LChQ/CIDh8+HDZu3OgeAVSR8AH/NW7cuNCrV8f/JeLIx9ixY90jgCoSPuC/4tTKwoULk8ARxfOCBQtMuQBUmZoP6KT2I061xBEPtR4A1a/5OKGLXxNKIwYOoQOgdky7AACpEj4AgFQJHwBAqoQPACDb4WPVqlXh2muvDSNHjgwNDQ3hhRde6PB6bJh63333hREjRoR+/fqFqVOnJs2bAAC6FT4OHDgQzj///PDEE090+vojjzwSfv/734c//OEP4e9//3s46aSTQmNjYzh48KA7DgBUvtR22rRpydGZOOrx2GOPhV//+tdh+vTpyXNPPfVUGDZsWDJC8sMf/tAtB4CSq2rNx6ZNm8LOnTuTqZZWseHIRRddFFavXt3p3zl06FDSmKT9AbVsILZ8+fLkDEABwkcMHlEc6WgvPm597Whz585NAkrrMXr06GpeErRZtGhRsmvtlClTknN8DEAJV7vMmTMnacXaemzdurXel0QBRyvi321ubm7btTaeZ86caQQEIO/hY/jw4cl5165dHZ6Pj1tfO1rfvn2THvDtD6j2aEVccdUaPFodPnw42cMFgByHjzFjxiQhY9myZW3PxRqOuOrlkksuqea3ogSqOVoxbty40KtXxx/3uGtt3DwOgIyHj/3794d169YlR2uRafzzli1bkr4fd911V3jooYfCX/7yl/Dee++FW2+9NekJMmPGjFpcPwVWzdGKuFHcwoULk8ARxfOCBQtsIAdQBw1H4vrYCqxYsSJceeWVX3j+tttuC4sXL06W295///3JL/rdu3eHyy67LMybNy+MHz++6lvyUmxxhCNOtbQPIDE0bN68uduhIX7NGF7iiIedawGqp5LP74rDR60JH7QXazziVEsc8WgdrWhqanKTADJG+KBQjFYAFCt8VNzhFNIWp0dMkQAUR937fAAA5SJ8UFfanQOUj/BB3Wh3DlBOVrtQmGW0AOSj4NTIB3Wh3TlAeQkf1IV25wDlJXxQF9qdA5SXmg86rceI0yJxdKLW9RcaiAEUg5oPcrMCJYabyZMnKzIFKBEjH3RrBUqaoyMAZJ+RD7rV0KurK1D05wCgJxSclkRXAkNXVqDE4NLc3NwWUuI57jrbWaDRvRSAzggfJdDVwNCVFShGRwDoKeGjBCpp6NXU1JTUeMTpmXiOj2s5OgJA+QgfJVBpQ6/jrUCp5ugIAOUkfJRAtRt6VWN0BIDystS2RNJs6BULWuNUSxzxaA07R4cUAMq51Fb4oGZ0LwUoj70VhI8TUrsqSieOrmhABsDR1HwAAKkSPgCAVAkfAECqhA8AIFXCBwCQKuGjRmyqBgCdEz5qwJbzAHBsmozVYMQjblnffm+T2OEztiHX8wKAoqqkyZiRjyqzqRoAHJ/wUWU2VQOA4xM+Mr6DLAAUjZqPGrGpGgBlstfGcvVnUzUA6JxpFwAgVcIHAJAq4QMASJXwAQCkSvgAAFIlfAAAqRI+CsIuugDkhfBRAHbRBSBPdDjNObvoApAFdrUt0ZSKXXQByBvTLjmfUrGLLgB5I3xkWBzpaG5uDi0tLcnjeJ45c2aHERC76AKQNyfU+wI4tuNNqcTQ0aqpqSk0NjYmz48dO7bDawCQNaUKH3HEIH6gx6mKPHxAt06ptA8gvXv3TgLG0eyiC0BelGbaJY/LUU2pAFBEpVhqm/flqPH6TakAkGWVfH6XYtqlq7UTWWVKBYAiKcW0i+WoAFDg8PHAAw+EhoaGDsfEiRNDPamdAIDsqMm0yznnnBPeeOON/32TE+o/u2M5KgBkQ01SQQwbw4cPD1mjdgIAClrzEQs8R44cGc4444xwyy23hC1bthzzvYcOHUoqZNsfAEBxVT18XHTRRWHx4sXh1VdfDfPnzw+bNm0Kl19+edi3b1+n7587d26yNKf1GD16dCiT420aBwBFVPM+H7t37056bPz2t79N6i46G/mIR6s48hEDSDX7fGRVbHTWundL7GS6cOHCTu8RABSpz0fNl9oOHDgwjB8/Pump0Zm+ffsmF9n+KIOubBoHAEVU8/Cxf//+8NFHH4URI0bU+lsVpvEZABRZ1cPHvffeG1auXJm0Lv/b3/4Wvv/97yetzG+++eZqf6tc0/gMgLKqeviI0wYxaEyYMCH84Ac/CIMHDw5r1qwJQ4YMqfa3yjWNzwAoq1JsLJdlNo0DoAhsLJcjGp8BUDal2FgOAMgO4QMASJXwAQCkSvgAAFIlfAAAqRI+urHRm83gAKD7hI+jNnqLm+BNmTIlOcfH3XkPAHBsmoy1G82IYaL9fiuxLXxsEx97cXT1PQBQRnuztKttkTZ6sxkcAPSc8FHBRm82gwOAnhM+KtjozWZwANBzaj66sdGbzeAAoCMby9V4ozebwQFA95l2AQBSJXwAAKkSPgCAVAkfAECqhA8AIFXCBwCQKuEDAEiV8AEApEr4AABSJXwAAMIHAFBcJ4SMOXLkSNsGNQBAPrR+brd+jucqfOzbty85jx49ut6XAgB043N8wIABx31Pw5GuRJQUtbS0hO3bt4dTTjklNDQ0VD2VxVCzdevW0L9//6p+bdzntPl5dp+Lxs90vu9zjBMxeIwcOTL06tUrXyMf8YK/bEv7noo3W/ioPfc5He6z+1w0fqbze5+/bMSjldUuAECqhA8AIFWlCh99+/YN999/f3LGfc47P8/uc9H4mS7Pfc5cwSkAUGylGvkAAOpP+AAAUiV8AACpEj4AgFSVJnw88cQT4etf/3r4yle+Ei666KLwj3/8o96XVDirVq0K1157bdLdLnanfeGFF+p9SYU0d+7ccMEFFyRdgIcOHRpmzJgR1q9fX+/LKpz58+eH8847r60R0yWXXBJeeeWVel9W4T388MPJ74+77rqr3pdSKA888EByX9sfEydOrNv1lCJ8PPvss2H27NnJ0qJ//vOf4fzzzw+NjY3h448/rvelFcqBAweSexuDHrWzcuXKMGvWrLBmzZrw+uuvh88//zxcc801yf2nemKn5fhBuHbt2vDOO++EKVOmhOnTp4cPPvjAba6Rt99+OyxYsCAJfVTfOeecE3bs2NF2vPXWWzX4Ll1TiqW2caQj/kvx//7v/9r2j4l97e+8887wq1/9qt6XV0gxVS9dujT5Vzm19cknnyQjIDGUXHHFFW53DQ0aNCg8+uijoampyX2usv3794dvfetbYd68eeGhhx4K3/jGN8Jjjz3mPldx5COORq9bty5kQeFHPj777LPkXy5Tp07tsH9MfLx69eq6XhtUw549e9o+GKmNw4cPh2eeeSYZXYrTL1RfHM373ve+1+F3NdW1YcOGZFr8jDPOCLfcckvYsmVLqJfMbSxXbf/+97+TXxzDhg3r8Hx8/OGHH9btuqAa4ihenBu/9NJLw7nnnuumVtl7772XhI2DBw+Gk08+ORnNO/vss93nKovBLk6Jx2kXajcDsHjx4jBhwoRkyuXBBx8Ml19+eXj//feT+rG0FT58QNH/tRh/edRz7rbI4i/qOEwdR5eef/75cNtttyXTWwJI9cRt3X/+858n9UtxQQC1MW3atLY/x5qaGEZOP/308Kc//aku04iFDx+nnnpq6N27d9i1a1eH5+Pj4cOH1+26oKfuuOOO8PLLLyerjGJxJNXXp0+fMHbs2OTPkyZNSv5l/vjjjydFkVRHnBaPxf+x3qNVHK2OP9exTu/QoUPJ73Cqa+DAgWH8+PFh48aNoR56leGXR/ylsWzZsg5D1fGxuVvyKNaIx+ARpwDefPPNMGbMmHpfUmnE3x3xw5Dqueqqq5LprTjC1Hp8+9vfTmoS4p8Fj9oV+H700UdhxIgRoR4KP/IRxWW2cbg0/kBfeOGFSQV1LBz70Y9+VO9LK9wPc/sUvWnTpuSXRyyEPO200+p6bUWbalmyZEl48cUXk7nanTt3Js8PGDAg9OvXr96XVxhz5sxJhqrjz+6+ffuSe75ixYrw2muv1fvSCiX+DB9dr3TSSSeFwYMHq2OqonvvvTfpwxSnWrZv3560nojB7uabbw71UIrwcdNNNyXLEe+7777kF3VcwvXqq69+oQiVnom9EK688soOoS+KwS8WOlG95lfR5MmTOzz/5JNPhttvv91trpI4FXDrrbcmxXkx2MV58hg8rr76aveY3Nm2bVsSND799NMwZMiQcNlllyW9guKf66EUfT4AgOwofM0HAJAtwgcAkCrhAwBIlfABAKRK+AAAUiV8AACpEj4AgFQJHwBAqoQPACBVwgcAkCrhAwBIlfABAIQ0/T9mBNNmCydnCgAAAABJRU5ErkJggg==",
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.plot(x,y, '.k')\n",
"plt.show()"
]
},
{
"cell_type": "code",
"execution_count": 18,
"id": "896d7633",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"2.5\n",
"17.115760540070784\n"
]
}
],
"source": [
"xSum = 0\n",
"ySum = 0\n",
"\n",
"for i in range(m):\n",
" #print(x[i])\n",
" xSum = xSum + x[i]\n",
" ySum = ySum + y[i]\n",
" pass\n",
"xMean = xSum/m\n",
"yMean = ySum/m\n",
"print(xMean)\n",
"print(yMean)"
]
},
{
"cell_type": "code",
"execution_count": 20,
"id": "3670b14b",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"theta 0: 4.012173816371986\n",
"theta 1: 5.24143468947952\n"
]
}
],
"source": [
"num = 0\n",
"den = 0\n",
"for i in range(m):\n",
" num = num + x[i]*(yMean-y[i])\n",
" den = den + x[i]*(xMean-x[i])\n",
" pass\n",
"theta1H = num/den\n",
"theta0H = yMean-theta1H*xMean\n",
"print(\"theta 0: \", theta0H)\n",
"print(\"theta 1: \", theta1H)"
]
},
{
"cell_type": "code",
"execution_count": 21,
"id": "729e6b3b",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAh8AAAGdCAYAAACyzRGfAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAO55JREFUeJzt3Qd4VHXWx/FfQAhVkA4GqdIs4CIiFgxFsCyW5bWigEYDLrIioCu7rsLKim0VCygoqLgqrgUQVyz0VURZEAQLDyARkKKglNDFvM+518kkIWWSzMydmfv9PM/sMDeTmeskmzlz/ud/TlJWVlaWAAAAoqRMtJ4IAACA4AMAAEQdmQ8AABBVBB8AACCqCD4AAEBUEXwAAICoIvgAAABRRfABAACi6hjFmF9//VWbN29W1apVlZSU5PXpAACAEFjP0j179qhBgwYqU6ZMfAUfFng0bNjQ69MAAAAlsHHjRqWkpMRX8GEZj8DJH3vssV6fDgAACMHu3bud5EHgfTyugo/AUosFHgQfAADEl1BKJig4BQAAUUXwAQAAoorgAwAARFXM1XyEup3nl19+0ZEjR7w+lYRVtmxZHXPMMWx3BgCEXdwFH4cOHdKWLVu0b98+r08l4VWqVEn169dX+fLlvT4VAEACiavgwxqQrV+/3vlUbk1M7E2RRmSRySxZkPfjjz86r/eJJ55YZMMYAAASMviwN0QLQGwfsX0qR+RUrFhR5cqV03fffee87hUqVODlBgCERVx+nOVTOK8zACB+xWXwAQAA4hfBBwAAiCqCjwQxf/58p/h2586dIX9P48aNNXbs2IieFwCg9DZt2qR58+Y514mA4CNK+vfv7wQHAwcOPOprgwYNcr5m9wEAIKdJkyapUaNG6tq1q3Ntt+MdwUcU2S6dqVOnav/+/dnHDhw4oFdeeUUnnHBCNE8FABAHNm3apPT0dGenp7HrAQMGFJgBiZcMSeIEH3v3upesrOCxQ4fcYwcP5n/f336YjsOH3WMHDoR23xL43e9+5wQgb731VvYx+7cFHqeddlr2sYMHD+pPf/qT6tSp42xxPeecc7RkyZJcj/Xuu++qRYsWzpbYLl26KCMj46jn++ijj3Tuuec697Hntcfca/8tAIC4sGbNmuzAI8C6e69duzauMyTFCj6efvppnXrqqdnj7jt16qRZs2bl+hRvSwg1a9ZUlSpV1Lt3b23btk1RUaWKe9m+PXjs4YfdY7femvu+deq4xzdsCB4bN849lpaW+76NG7vHv/46eOyFF0p8mjfeeKOef/757NuTJ0/WDTfckOs+d955p9588029+OKLWrZsmZo3b66ePXvqp59+cr6+ceNG/eEPf1CvXr20fPly3XTTTbrrrrtyPca6det0wQUXOD+DL774Qq+99poTjNya97UAAMSsE/Np8miNNu19oUQZEvt6ZqbiKvhISUnRAw88oKVLl+p///ufE11deuml+vLLL52v33777Zo5c6Zef/11LViwQJs3b3beJBF03XXXOUGANe+yy8cff+wcC7DMhAV5Dz/8sC688EK1adNGzz77rJO9CESx9vVmzZrpn//8p1q2bKk+ffocVS8yZswY5/iQIUOcX96zzjpLTzzxhKZMmeIEiQCA2JeSkqKJEyc6AYex6wkTJjjHS5QhefRR6aqrcmfzY73DqX3Szukf//iH80a4ePFi54WwN0erX7CgxNgn/NatWztfP/PMMxVRgUguZ+fTO+6QhgyRjsnzn/nDD+51xYrBY4MGSTffbD/Z3PcNLGfkvG8pCkNr166tiy++WC+88ILTxtz+XatWrVwZi8OHD+vss8/OPmadRs844wx9/Vv2xa47duyY63EtC5XTihUrnIzHyy+/nH3Mni/Qot5+LgCA2JeWluZkvy2QsIxH3sAjZ4YkZwCSX4ZE6enSiy9Kq1ZJp56quGuvbhGVZTjsk7q98Vk2xN40u3fvnn2fVq1aOfUMn3zySYHBh9U32CVg9+7dJTuhypWPPmYD0fIbipbffcuVcy+h3rcUbOklsPwxzpZ7IiAzM9NJuVmdR14UtwJAfElJSck36MibIbG/+/b+nJ0hqVtXeucd6fe/d+947LHS558f/aE8yor97CtXrnSCDUvdW13HtGnTnKUBqz2wQW/Vq1fPdf+6detq69atBT6eLQ+MGjVKfmK1GDYvxbbXWjSbky2n2OtoyzFWMGQsqLOCU1tCMZa1ePvtt3N9n2WX8ha3fvXVV0dHvQAAf2RIateWzjtP+uQTaebMYADiceBRot0uVmNggcann36qW265Rf369XPe5EpqxIgR2rVrV/bFiikTnUWktnRir1tgHS+gcuXKzut6xx136L333nPuc/PNN2vfvn3OL5axXiG2vmf3Wb16tbPUZcs4Of35z3/WokWLnAyL/bzs/jNmzKDgFAASWEpKilJTU90sSXKydMYZUrVq0i+/KJYUO/yxT+WBT9Pt27d3PpE//vjjuuqqq5xP89ZhM2f2w3a71KtXr8DHS05Odi5+Y7uFCmJFvbZud/3112vPnj06/fTT9f777+u4447LXjax3TBW4Pvkk0869SD333+/s5wTYLuSrOj3r3/9q7Pd1uo9LKtiPycAQILat8/qIqSqVd3bDz5ou0Gk3zLpsSIpy96VSsGKS+3N0AIQK6Z89dVXne2dxj6VW91HYTUfeVnNR7Vq1ZwsSN43aFvqsWLJJk2aMOI9Cni9ASCOfPmlu5PFCklts0FSUlSfvrD371JlPmyJxLZ/WrBhn8gt3W8zRexTuT2hLQsMHTpUNWrUcJ548ODBTn1IxHe6AAAQpzZt2uQsjduOlcKKSotkGza++UbasUOyWsv69RWrihV8/PDDD+rbt6+2bNniBBuW2rfA4/zzz3e+/thjjzlbfSzzYTtYrPBl/PjxkTp3AADi2qRJk7Kbg9n7p+1YCdT3hcQWLwIZDmu58OqrUufOtttDsazUyy7hxrJL7GDZBQAim/Fo1KjRUb05bFxGSBmQRYvceo4ZM6RCaitjcdklcWa7AACQoHNbjmLfN3iw9Nln0t13K94QfAAAEMNzW/Jl3/evf0k2G8xapscZgg8AAGJ4bks2G+T6xhvB2zYmY/Jkt2tpMZd75s2bd/TQuSjyvs0ZAAA+lRbC3BbHBx9IF13kTlk/7TRrh+1NgWuYEHwAABDDc1scNrDVdrFYD4/jj1dJWKYjEHgYu7ZZMBb8lGqLbwkQfAAAEItmz3aDDqvvsHksH36Y/7DUMBS4Rjv4oOYjSvr37+8MkrNLuXLlnIF71h9l8uTJR/0yFMZmuOQd3gcASDC33CJZD63HHw8eK0XgUeoC1zAj+IjyNFtr0GZ7uGfNmqUuXbrotttu0+9//3v9EmNDfwAAHjr1VDfjYbNavCpwjSBfBx/Rrvi1AXo2ZO/44493Rt7/5S9/cSbNWiASmEr76KOP6pRTTnGm2zZs2FB//OMflZmZ6XzNWtnfcMMNTgOXQBZl5MiRztdeeuklZwBd1apVnee49tprnY60AIA4kJUl7dwZvD1woLR8ufTXv4b1aay41D4A23ufXXtRbOrr4MMqfq2znA3Gs2u77QV7/rZt2+qtt95ybltK7IknntCXX36pF198UXPnztWdd97pfO2ss87S2LFjnc5xlkGxy/Dhw52vHT58WPfdd59WrFih6dOnO79UttQDAIhxO3ZIl14qXXih/TF3j1nL9FNOicjTWaYjNTXVk4yHrwtOY6ni19jk3y+++ML595AhQ7KPN27cWKNHj9bAgQOdGTnly5d3WtdaxsOyGzndeOON2f9u2rSpE8B06NDByZpUsa1ZAIDYtGePtHChtH+/tGSJfdJUovNl5qNULW0jwMbrWEBhZs+erW7dujlLM7aEcv3112vHjh3aV8S639KlS9WrVy9n4rB933nnnecc37BhQ1T+GwAAJdS4sTsQ7tNPfRF4+Db4iKWKX/P111+rSZMmzlKJFZ/atOA333zTCSjGjRvn3OfQoUMFfv/evXudrI0tx7z88stasmSJpk2bVuT3AQA8sHGj1KOHtGJF8JgtubRr55sfhy+Dj1iq+LWajpUrV6p3795OsGEZmX/+858688wz1aJFC23evDnX/W3pxbI0OX3zzTdOduSBBx7Queee6yzjUGwKADFqxAi3Z8fNN7uFpj7ky5qPYrW0DaODBw9q69atTvCwbds2vffeexozZoyT7ejbt69WrVrlFI4++eSTzhLKxx9/rGeeeSbXY1gdiNVxzJkzxylUrVSpkrPUYkGJfZ/Vh9jjWPEpACAGjR3r1nnYQLjfltx9JyvG7Nq1y8JA5zqv/fv3Z3311VfOdbzp16+f899ll2OOOSardu3aWd27d8+aPHly1pEjR7Lv9+ijj2bVr18/q2LFilk9e/bMmjJlivM9P//8c/Z9Bg4cmFWzZk3n+L333usce+WVV7IaN26clZycnNWpU6est99+2/n6559/XuJzjufXGwBixtdfZ2VNnJiV6HYV8v6dV5L9j2LI7t27nR0d1svCahhyOnDggNavX+/UR1SoUMGzc/QLXm8AKKWMDOmkk+wPqq2zS79tBkhEhb1/5+XbZRcAAKKyk+Wqq2zrodSiBS/4bwg+AAAIp2XLpNatpYoVnb5Sa666Sie2bq2U+vV5nf282wUAgIiYMEHq2FEaOjTYSfuCC9SoSRPPOmnHIoIPAADCpUkT6ZdftGnjxnw7aUdrllisI/gAAKA0fhv+6bDmYZ9+qjVDh8ZUJ+1YE5fBR4xt0ElYvM4AUAgbAmcNw2wA3M8/B4+fcYZObNEipjppx5q4Cj7KlSvnXBc15wThEXidA687ACAH2z77+uvudto334zZTtqxKK76fBgbI79z507VqVPH6e4ZGMiG8LFfCQs8rEV79erVVZ8KbQDI3//+J333ndS7d75fthqPaHbSjpc+H3EXfNjpWotyC0AQWRZ41KtXjwAPAIxlg4cMkXr1ci/wT5Mxy3TYJ3HLfNgcFESGLbUE0oUAkCicvhtr1jjTzYudiXjiCenZZ6Xp06X166XKlSN1mgkv7oKPAHtj5M0RABAq67MR2P5qxaBWk2FDRkN2++3Sxx+72Q8Cj1KJu2UXAABKkvGwhl85t7/aB9iMjIyCMyC2g+Wll6TBg/07fbYYEnrZBQCA4rKlloL6buQbfNhOlg4dpHXrpCpVpBtv5EX361ZbAABKwmo8itV3wyan33ST1KyZ1LYtL3qYEXwAABJeSH03Nm+Wtm4N3r7zTunzz6X27T0448RGzQcAwDcK7Lsxf750xRVSu3bS++9LebIkKBo1HwAA5MMCjnxrPOrWdft4bN8u7dgh1a7N6xdBFJwCAPzJikqttsO0bi3Nni2ddlrwGCKGvBIAwH/+9S+paVPbBhM81qkTgUeUEHwAAPzFttxOmmTDwqTHH/f6bHyJZRcAgL9YMallPqyB2B13eH02vkTwAQBIbNbI+8kn3S6l1q3UHH+8dNddXp+ZbxF8AAASe2jcf/4j3XabTcyUevSQWrb08jRBzQcAIFGGxtnslq5duzrXdjvbxRdLffpIjz0mtWjh5WniNzQZAwAk3tC4MmWUYXNbmjQJLr0wHC5mmoyx2wUAkHhD4379VWv/+tfgAQKPmELwAQBIvKFxkpqfdZZn54TCEXwAAOJaSq1amnj//bmHxo0dq5Rbb434cs+8efOcaxQPu10AAPHr22+lyy9X2sGD6vnVV1q7efPRQ+MiwApa09PTneUey7rYxNy0tLSIPmcioeAUABC/bAhc27bSoUPSBx+4U2m9KHAtW1YZGRkRD3piGQWnAIDEdfhw8N81a0ozZkgrVkQl8CiwwPXIEa1duzYqz58IqPkAAMRPXcRnn0lt2kjz5gWPtW8v1a8vTwtcy5Z1lnsQGoIPAED0G3+V1HPPSZZh+Nvf3N4dHrClFavxyFXgOmGCr5dciouaDwBA/NRF7N0r3XuvdPfdUvXqinqr9jz3saWWaBS4xgNqPgAAngpbXcSHH0o5m4VVriw98khEA49QMzYWcKSmphJ4RHrZZcyYMerQoYOqVq2qOnXq6LLLLtPq1atz3cd+EElJSbkuAwcOLMm5AQDiVFjqImwb7QUXSPff7w6HiwLLZgS20Bq7HjBgAL08vAw+FixYoEGDBmnx4sX68MMPdfjwYfXo0UN7LQ2Ww80336wtW7ZkXx566KFwnzcAIIaFpS6iaVPpzjsl+wDbtauigZ0scVDz8eOPPzoZEAtKOnfunJ35aNeuncaOHRvxNSMAQGwrdl3EW2/ZG4lUo4YnA+Ho4REHNR/2BKZG4JfkNy+//LJq1aqlk08+WSNGjNC+ffsKfIyDBw86J5zzAgBIDMWqi/j736XevaWbbgruZInyQDh2ssR45sPWwS655BLt3LlTH330UfZxS7NZgU6DBg30xRdf6M9//rPOOOMMvWXRbD5GjhypUaNGHXWczAcA+MzSpdI550jDhkn2vvDbko0X2MkS2cxHiYOPW265RbNmzXICj8Ii2rlz56pbt25O2q1Zs2b5Zj7skvPkGzZsSPABAInO3n6++05q3Dh4bMuWqDYMgzfBR4kGy91666165513tHDhwiJTaR07dnSuCwo+kpOTnQsAwEdsif2GG6T587Xpvfe0JjOz0J4aSCzFqvmwJIkFHtOmTXMyGk2aNCnye5YvX+5c1yeSBQAEVKjgZD0m7dqlRh07hrcLKmJesZZd/vjHP+qVV17RjBkz1LJly+zjlmapWLGi1q1b53z9oosuUs2aNZ2aj9tvv92JZG1HTCjY7QIACcp6Z1gB6W9FpJs++kiNzjuP6bAJImK7XZ5++mnnQa1y2TIZgctrr73mfL18+fKaPXu20/ujVatWGjZsmHr37q2ZM2eW7r8IABDfvv9e6tbN2odmH1pz+DDTYX2qWDUfRSVJrFA01AwHAMBH/v1vp75DX30lXXutVKlSdhfUvPNfmA6b+JhqCwCIvNtucy///a8TeBh6avgXU20BAOG3Zo301FPSY49JeWa85EVPjcQQ8a22AAAUaP9+6dxzpW3bbD1eGj680BfLMiBssfUXll0AAOFVsaKNQZfOO0+65hpeXRyF4AMAUHrLlrlLLQH9+1uLa+n443l1cRSCDwDAUTUY8+bNc65DYrO7zjxTuuoqm5nhHrNeHkXUekT0nBDTCD4AANmsw6h1Gi1Wx1Ebo2EFhjajJcesLk/PCTGN3S4AAIdlFezNPW/fjYyMjKMLQjdvlho0CN7OyJAaNcruXurJOSExO5wCABLXmjVriu44euSI9Le/SU2bSkuWBI9b1iPMgUfI54S4Q/ABAHAEOo7mdFTHUfv611+7yytRGJ0R0jkh7hB8AACK7jgaGK9h2Y1nn5XeeEP6+98j/srRBTUxUfMBAD5iNRS2lGEZhYJqJnJ1HK1ZUxo2TKpaVXrwwaifb77nRK1H3Nd8EHwAgE/YLpH09HSnhsKWMizLkZaWVvg3vf++dMEF7nLL6tUSyx0oAMEHACB8u0buuUc6+2ypZ09eVRSI3S4AgJLtGtm1y11mycwMHrPaDgIPhBGD5QDABwK7RvJmPo7aNXLJJdLChW4Q8txz0T9R+AK7XQDAB0LeNXLffW5dx803e3Oi8AUKTgHAR47aNbJ1qx2UTj89eKfDh6Vy5bw8TSR4zQfLLgDgIxZwZGc7VqyQevSwNIj779q13eMEHogwll0AwK9secX6eNSq5dZ4AFFC5gMA/GT7djfYMJUrS+++K9WtK1WsGPZmZUBByHwAgF+89prUrJn0+uu5B8IVM/BgxD1Ki+ADAPziiy+sKlB68cUSP4RlPAJdUo1dDxgwwDkOhIplFwBIZDYQLjDqfuRIqWFD6aabItKsjOUXhIrMBwAkatAxbpx07bXBibS2i2XgQOmYkn/uZMQ9woHgAwASkbVNv/12aepUaebMsD0sI+4RDjQZA4BEZZkPaxh2223BpZcwYcQ98mKqLQD4zZEj0sMPS336uHUdQJQx1RYA/MaWWEaMcGs8LBABYhg1HwDgIVu+mDdvXum3qv7pT1KDBlJ6utsuHYhhBB8A4JFSNes6eFBatCh4u3lzbVqwQPNSUui5gZhH8AEAHihVsy5rkd6pk9S9u/TVV8FApmXLkgUyQJQRfACABwpr1lWkGjWkOnXc2Sxbt9J1FHGH4AMAPKjVKHazrszMYCGpfd+UKdLy5VLXrqULZAAPEHwAQDGFY7BasZp1LV0qnXaa9MADwWOW+Tj+eOefdB1FvKHJGAAUg2U6LODImWmwwCEjI6NEs01Catb1wgvSDTdITZpIX36Z7xRaC4CsZsQyHoFAJi0tjZ8tYrLPB4PlAMDDwWr2PUV+X79+0p49bgOxfAIPY4FGz549iw5kgBjAsgsAFENUljjmzpUuvtjdTmusNfrgwW6haSEs4EhNTSXwQMwj+ACAWBqstm+fdM010rvvSo88ws8GCYmaDwCItcFq//mPNH26NHasu50WiAMMlgOAePLWW+4wuA4dvD4ToMQYLAcA8eL556XevaWrr5Z27fL6bICooOYDALx0+eXuFtorrpAqVeJnAV9gqy0ARFNWlvTZZ1LHju7t6tWlL76QqlTh5wDfIPMBANFi7dFteeXMM6VZs4LHCTzgMwQfABAttj23dm3pmGOkjIyozpIBYgnBBwBEknVDPXAgeNt6d3z6qXTLLVGdJQPEEvp8AECkbN0q9e0r1a8vvfhiTMySASKFrbYAEAM2LVqkebNna9O//y19+23YZ8kA8YplFwCIAGep5Ior1DUrS40OHdKkefNid5YMEGUEHwAQLpbduPRSbVqxQunp6dkZC7u2cfclKRaN+CwZwAP0+QCAcPXvuO466ZNPtCYzs8ClkpIEDWlpaerZs2fkZskAsZz5GDNmjDp06KCqVauqTp06uuyyy7R69epc9zlw4IAGDRqkmjVrqkqVKurdu7e2bdsW7vMGgNhiY+8nTpTOP18n3n9/2JdKLOBITU0l8ID/go8FCxY4gcXixYv14Ycf6vDhw+rRo4f27t2bfZ/bb79dM2fO1Ouvv+7cf/PmzfrDH/4QiXMHAG9ZZ1KbQBtw8snSBx8opWNHlkqASG21/fHHH50MiAUZnTt31q5du1S7dm298sor+r//+z/nPt98841at26tTz75RGdaV78wbtUBAM8sXiylpkrJydLy5e58ljysxoOlEvjF7mK8f5eq5sOewNSoUcO5Xrp0qZMN6d69e/Z9WrVqpRNOOKHA4OPgwYPOJefJA0DMO/10qX17+wMoVa1a4FIJ9RlAGHe7WDHVkCFDdPbZZ+tkSzU6/XS2qnz58qpug5JyqFu3rvO1gupILFIKXBo2bFjSUwKAyFq50i0sNdYi/d13pbfflmrV4pUHohF8WO3HqlWrNHXqVJXGiBEjnAxK4LJx48ZSPR4ARMQDD0jt2knjxgWPVavmFpoCKJYSLbvceuuteuedd7Rw4cJcKcV69erp0KFD2rlzZ67sh+12sa/lJzk52bkAQEyrVMmd02LZDwDRy3xYbaoFHtOmTdPcuXPVJE+BVfv27VWuXDnNmTMn+5htxd2wYYM6depUujMFgCjJniCbsyX64MHS7NnShAn8HIBoZj5sqcV2ssyYMcPp9RGo47BajYoVKzrX1gxn6NChThGqVbsOHjzYCTxC2ekCAAUFAzbjxFqNR7qA09qiB7qT2qeziRMmKC093V1e6daNHxAQ7a22SQWsbT7//PPq379/dpOxYcOG6dVXX3V2sVhXvvHjxxe47JIXW20BFBgMlCnj9M+wDzmRCFDynSBbpowyvvuOXStAGN+/S9XnIxIIPgAUd5x8qAFKUWyppWvXrvket+6iAMLz/s1gOQAxK5Rx8haglHqI2549NkBFJx45wgRZIAoIPgDErFDGyYcSoBRp2DBp8mSlDB+uic88wwRZIMIIPgDErFDGyYcSoBRp9GipY0enh0fazTc7yzq21GLXJVm+AVA4aj4AxLyiZqRYzYcttVjGIxCgFBo0/PijFXJIV14ZPGblbzQMA0qMglMAvhPyELdt26TTTpN++EGaP18655xoniaQsKI2WA4AYkXIQ9zq1nX7dSxbJuWZQwUgOgg+ACS+DRukOnWkChXc208/LVmdiLVMBxB1FJwCSGw2dfbUU6U77wweq1KFwAPwEMEHgMRWrpy0a5e0ZIm1YPb6bAAQfABISL/8Evz3hRdK77wjLVwYXHYB4CkyHwASh22XnThRattW2rkzePzii90MCICYQPABIHHs2yeNGSN99ZX0zDNenw2AArDbBUDiqFxZmjpV+u9/paFDvT4bAAUg+AAQv2ymy4MPSu3aubUdxtqk2wVAzCL4ABC/xo+X/vIXqXZt6ZtvpBo1vD4jACGg5gNA/LrpJumss6QHHpCOO87rswEQIjIfAOLHoUPSm29K11zj3rats1bfkWeqLYDYRvABIH56d5x3nrR4sXT4sNS3r3ucwAOIO3xcABAfjjlGuugid3mFJRYgriVlZVlXnvgcyQsgwe3d67ZEr1nTvX3kiPTDD1L9+l6fGYBSvH+T+QAQm1aulNq3l/r0cbfUmrJlCTyABEDNB4DYZLUc330n7dkjbdwoNWrk9RkBCBOCDwCxwzIcgQLSk06Spk2TTj9dqlXL6zMDEEYsuwCIDTZ19pRTpG+/DR674AICDyABEXwA8J7Vvd9zjzsQ7u67vT4bABFG8AHAe0lJ0pQp0qBB0sSJXp8NgAgj+ABQIps2bdK8efOc6xJ5+23pueeCt084QXrqKalKFX4iQIIj+ABQbJMmTVKjRo3UtWtX59puF8uCBdKll7qZjlWr4itoAlBqBB8AisXetNPT0/Xrb7037HrAgAHFezPv3NkNPgYPllq0iJ+gCUBYEHwAKJY1a9ZkBx4BR44c0dq1awsvKLVlFpvPEqjxsAFxjzwilS8fP0ETgLAg+ABQLCeeeKLK5BnmVrZsWTVv3rzgb7LlFct0jByZ85tiP2gCEBEEHwCKJSUlRRMnTnQCDmPXEyZMcI4XyKbR2v09LCYtUdAEICIYLAegRGy5wrIG9uZ9VOBhGYbt26U6dYLH1q2TmjXL93EsK2HBQaEBTBhYjYcttVjGIxA0paWlRfQ5Ab/YXYzBcgQfAMLrxx+lvn2lDRukJUukSpUKDQYCdRiWlbCMSqSDgUKDJgAlRvABwDs28r5tW2nnTm2aMkVratXKN6thQYDtOMlZh2HZiIyMDIICIMGDD2o+AJSe7WYJsKWW11/XpBEj1Ojqqwvc1koBKOBfBB8ASicjw+3b8d//Zh/a1Lix0keNKnRbKwWggH8RfAAonQcflD76yN1O+1uwEUpWo0S7ZgAkhGO8PgEA0RORnSUPPyzt2SONHi39tpU1kNXIW8+Rd1urFZf27NmTAlDAZ8h8AD4RttbiNovloYeCt613x7/+JTVuXKKshh1LTU0l4wH4CFttAR8I286SzZsly17s3y/NmCFdckmRz8u2VsAfdhdjtwvLLoAPFFaDUazgo0ED6Y9/lL78UjrzzCLvbo9NDQeAvAg+AB8ItQYjX59+KrVsKVWv7t4eM8ZtlZ6nVTkAhIq/HkAcsOWLefPmlXgCa4l3lkycKJ19tpSeHuzlUa4cgQeAUiH4AHxSKGo7S6zGw4IYuw6pjXm7dlJSkhtwHD5coucFgLwoOAVimCctyH/6SapRI3jb6jvatHGDEAAoAO3VgQQR1Rbkhw5Jw4ZJrVq5u1oCTjqJwANAWLHsAsSwqLYgt5qOuXPdqbRvvx3+xweA3xB8ADEsqi3Ik5OlqVO16bnnNK9lyxIXtwJAUQg+gBhXokLRUGRmSjfcIL3wQvahSR99pEbp6aXvggoAhaDgFPCrJ5+U/vQnyToRfvedNmVmRr+4FUDCoOAU8KFi9wKxTqXXXSfNnOk0EItqcSsAXyv2ssvChQvVq1cvNWjQQElJSZo+fXqur/fv3985nvNywQUXhPOcAZSkF8j27dJ992WPvXe6lL70ktS5c/SLWwH4WrGDj71796pt27YaN25cgfexYGPLli3Zl1dffbW05wmgAJbpSE9Pz85a2PWAAQNyZ0CsQdhZZ0n33CM9+qj3xa0AfK3Ys10uvPBC51KY5ORk1atXrzTnBSCcQ+OsQ+nw4W7g0b17gY9lxaw9e/ZkEi2A+BssN3/+fNWpU0fHHXeckwYePXq0atasme99Dx486FxyFqwACMPQuEqVpIwMqXFj9+DNN0vXXy9VrFjo4zGJFkDcbbW1JZcpU6Zozpw5evDBB7VgwQInU2KfxPIzZswYVatWLfvSsGHDcJ8SkNDyXS4ZOlQpVmt15ZVu51Jj7dGLCDwAIOa32lox6bRp03TZZZcVeJ9vv/1WzZo10+zZs9WtW7eQMh8WgOzatUvH2hZAACGxGg9barEC0RT7v3XbtlKzZu5uFpZBAcTQVtuILLvk1LRpU9WqVcv5o5hf8GH1IXYB4u2N3motbMkjVgoyU6pVU0pqavDA/PnunJby5b08LQCIfodT+yO9Y8cO1a9fP9JPBcTUiPti990oKctyTJ4sNWokLVsWPH7qqQQeABIj+MjMzNTy5cudi1m/fr3z7w0bNjhfu+OOO7R48WKnK6LVfVx66aVOGtgq6AFfbGstRoASNra08vPP0jPPRPZ5AMCLmg/bydKlS5ejjvfr109PP/20U//x+eefa+fOnU4jsh49eui+++5T3bp1w75mBESbZTIsoMjveOpvSx4WiES9TflPP0lTprjt0vM0CgOAuK/5sD+whcUr77//fnEfEoj/ba05uoCG1HejNOyxrV+HDYYbOdI9VqOGNGRI6R8bAKKAj0hAMYTSBTTibcoXLpTuuEP6+9+lFSv4+QGIOwQfQJhH3Ee8Tbkt79jyitV3WFEpAPipz0ckUPOBhOy7UZrAw+ayPP64dMstUuXK4TxFAEjMPh+AX4WtTfnVV0tvvSWtXi09+2w4Tg0APMWyCxDrBg1yC0rZrg4gQZD5AGKtW+q+fe5AuDZt3DvY1l67XbWq16cKAGFB5gPw0FHNyB54QOrQQerRQ9qxI3hHAg8ACYTgA4ilbql3361N+/e7vTws2wEACYjgA/BIgc3IrHGY9e9o356fDYCERPABeKTAZmRW41G7Nj8XAAmL4APwiNOMrEcPua3IItCMDABiFLtdAA+lzZihnnfdpbXduql527YEHgB8geADiKZ335VmzpTGj5eSkqTy5ZXy6KMi1wHATwg+kLj9MmLN999Ll18uHToknXee27kUAHyImg8kZr+MSZMUc44/Xrr/fnconAUhAOBTDJZDQmQ8LODIuW3Vijdt4qznGZBXXnGzHBZ4AEAC212MwXJkPpC4/TLWrpWnRo+W+vSRrrvOTsjbcwGAGELwgcTtl9G8uTx15ZWSRf+pqd6eBwDEGIIPJEa/jIkTnYDD034ZWVnSV18Fb7doIa1fL917r51UdM8FAGIYNR9IqNoPW2qxjEfUA4/MTOnaa6U5c6Rly6SWLaP7/AAQRzUfbLVFwrCAw7MC04oVpX373NoOm8tC8AEABSL4AErql18kqzWxiy2rvPSS9MMPUtu2vKYAUAhqPoCS2LBB6tJFeuSR4LH69Qk8ACAEBB9ASVhtx0cfSQ8+aAudvIYAUAwsuwAl0b+/9N130vXXu9tpAQAhI/MBhOLrr6V+/dy5LMaGwo0cKTVrxusHAMVE5gMoigUcPXtKGzdKJ5wg3XcfrxkAlAKZD6Ao5ctLTz4pnX++NGgQrxcAlBKZDyA/n30mlSsnnXaae/vSS6VLLnGXWwAApULmA8hrxgzp7LPd2Sx79gSPE3gAQFgQfAB5nXuu27PDsh55puUCAEqPZRfAfPut1LSp+1rUqOEuu9StS7YDACKAzAf8zSbR3nWXO4F29uzg8Xr1CDwAIEIIPuBvVsfx88/uQLgFC7w+GwDwBZZd4E8WbNgwOPPYY+5ulosu8vqsAMAXyHzAX/bulW66SUpPDx6rVInAAwCiiOAD/vL559LkydLzz0srV3p9NgDgSyy7wF/OOUd65BGpXTvplFO8PhsA8CUyH0hsP/0kDRggbd8ePDZ0qNS1q5dnBQC+RuYDie2aa6QPPpB27JDeeMPrswEAEHwg4T34oLR5szRiRMjfsmnTJq1Zs0YnnniiUlJSInp6AOBHLLsgsXz/vTR/fvC21XasWCG1bx/St0+aNEmNGjVS165dnWu7DQAIr6SsLGvxGDt2796tatWqadeuXTr22GO9Ph3Eky+/lM47z+3hsXy51KhRsb7dMh4WcPyaY55L2bJllZGRQQYEAML4/k3mA4nDWqQ3by41aSL98kuxv92WWnIGHubIkSNau3ZtGE8SAEDBKeKbLbM0aOC2SS9XTpo+XTruOCk5udgPZTUeZcqUOSrz0dwCGgBA2JD5QPx6+WWpVStpwoTcA+FKEHgYKy6dOHGiE3AYu54wYQJLLgAQZgQfiF9bt0qZmdLbb7vTacMgLS3NqfGYN2+ec223AQDhxbIL4osFGbbEYm6/3c10XH118FgYWAaELbYAEDlkPhA/QcfYsdIFF7i7WUyZMlKfPsHptACAuEDwgfiwaZP0t7+53Upff93rswEAlALLLogPDRu6haU7d0pXXeX12QAASoHgA7HJ+nSMHu0GGq1bu8euvdbrswIAeLHssnDhQvXq1UsNGjRQUlKSpltfhRysYeo999yj+vXrq2LFiurevbvTvAkoFpvFMmqUW0x6+DAvHgD4OfjYu3ev2rZtq3HjxuX79YceekhPPPGEnnnmGX366aeqXLmyevbsqQMHDoTjfOEXNvbemnv95S9u8zAAQMIo1WwXy3xMmzZNl112mXPbHsoyIsOGDdPw4cOdY9bjvW7dunrhhRd0tX2KLQKzXXxq/37pv/+VevTIvfRyDCuDABAPPJvtsn79em3dutVZagmwE+nYsaM++eSTfL/n4MGDzgnnvMBnrIi0Y0fpooukxYuDxyMQeNjwOGsgZtcAAG+ENfiwwMNYpiMnux34Wl5jxoxxApTApaHtaoC/VKsmtWkj1azpZkAiZNKkSc7U2q5duzrXdhsA4MM+HyNGjHBSNIHLxo0bvT4lRMOuXZb2cv9t3UltG+2KFVKXLhHJVtj3pqenZw+Ns+sBAwaQAQGAeA8+6lmra0nbtm3LddxuB76WV3JysrM2lPOCBPfZZ9Jpp7k7WnJmP/L5HQlXtsJ2XOWcVmuOHDmitWvXlujxAAAxEnw0adLECTLmzJmTfcxqOGzXS6dOncL5VIhnFpyuXy/NmCHt2ROVbMWJJ56oMtaOPQebWtvcdtQAAGI7+MjMzNTy5cudS6DI1P69YcMGZ/fLkCFDNHr0aL399ttauXKl+vbt6+yACeyIgU/l3FTVq5c0ZYq0bJlUtWpUshU2KG7ixIlOwGHsesKECQyQA4B42Go7f/58dclnXb5fv37Odlp7uHvvvdf5Q79z506dc845Gj9+vFq0aBHS47PVNgG9/750333SrFmFBht5WYbDllpyBiAWNNio+5JOnbXHtODFMh5MrgWA8CnO+3ep+nxEAsFHgrGi0latpIwMt2HYP/5RrG+3Gg9barGMRyBbkZaWFrHTBQCUDMEHYsvHH0uvvio9/LBUsWKxv51sBQDEPoIPeGvqVKlWLSlHszkAQGLbXYxlF3pXI7z+9S/p+uulOnWklSvdawAAYqnJGBJM795S27bSgAFSjRpF3p125wDgPwQfKB2rV87R18Wp6fj0U+nvfy9yNgvtzgHAnwg+UHK2BfaKK9zaDisoDUhOLvJbaXcOAP5F8IFS/PaUkU46SSpXzp3VUgy0OwcA/6LgFMVz5Ii0b1+wWdjf/iZdeaUbhJSg3XneBmK0OweAxEfmA6H7/nupWzfp2muD7dKtrqOYgYeh3TkA+BeZD+Rbj2HLIpadyNWCfPt2afFid5ll9Wq3c2kpWKfSnj170u4cAHyGzAdC34FiW2hfeskdCFfKwCPAgpvU1FTmrACAjzDbBYUPcpOUsWiRUjp1Ci07AgDwpd3F6HBK5sNHimrole8OFElrhw7NdYz+HACA0iD48IlQAobADpScyiYlqfkTT5SoPwfdSwEA+SH48IFQA4aUbds0MS3N2fJqnBH2zz6rlA4dit2fg+wIAKAg7HbxgcIChux6jc8/lzp1UlqZMuo5a5bWlivn9NzIW88RSn+OgoId29lCfQgAgMyHD+S7nJK3oVe7dlKPHtLFFyvl9NML3IESSn8OupcCAArDbhefsGUQyz5YxiMQMKSdfLJ0+ukWQbh3ss6lNhguKanIx7PshmVO8suO5LtrpmxZZWRkkPkAgATFbhfk29DL3vxtt4tdp23b5iyzaPTo4J0qVQop8CiqPwfdSwEAhaHmw0csKMgOFho1clukb9niXocYdISK7qUAgIIQfPjJ3r1S5cruv/v0kZo1k848MzrBDgAAv6Hg1A/275duucVdZrF/B0Qw8AAAoCAEH36wa5f01lvSypXSBx94fTYAAJ9j2cUP6tWTXn7Zre04/3yvzwYA4HNkPhLRzp3StddKixYFj3XvTuABAIgJZD4S0ahR0quvSkuWSF9/LR3DjxkAEDvIfESIp0PVLPiw5RVbaiHwAADEGIKPCIj6ULWtW6Vnnw3ePvZYt7D0jDMi+7wAAJQA7dXDLOqtxX/6SWrTRrKOpTNnSr//ffifAwCAItBe3UNRH6pWo4Z09dXSqae6TcMAAIhxVCKGWSgj50tt3Tqpdm13ecU8+KC7jbZChfA9BwAAEULNR5hFfKja9OnSaadJAwe6AYdJTibwAADEDTIfERDRoWp16kj79kkbN7rXgVktAADECYKPeBiqdvCgm90wZ50lzZ3rXrONFgAQh1h2iWW2rPLUU1KLFu522oDOnQk8AABxi+Ajlh06JE2cKG3Y4F4DAJAAWHaJZbbU8tpr0ocfSoMHe302AACEBcFHLDlyRBo92vbruoPhTOvW7gUAgARB8BFLXnhBGjlSqlpV6tZNqlvX6zMCACDsCD5iSb9+0rRpbsdSAg8AQIKi4NRLtoV2woRgszDbOmvzWa67Lr6m6AIAUAwEH16x9uupqW6n0rFjg8eTkmJ/ii4AAKVA8OGVMmWk66+XatWSWrYs8cNYpiM9PT17loxdDxgwgAwIACBmEXxE0+7d0vffB2/fcov09dfSRReVeEkl6lN0AQAoJYKPaPniC+l3v5P+7/+kw4eDSyyW+SjFkkpgim5OYZ+iCwBAGBF8RMuxx0rbt0ubN7tD4cK0pBLxKboAAIQZW20j6ZdfgjNYGjeW/vMfqU0b6bjjQvr2wpZUcgYXEZ2iCwBAmPkq8xHV7ahz5rgD4Wy5JeDss0MOPIq7pGIBR2pqKoEHACDm+Sb4iPp21CeekNavl0aNKvFDsKQCAEhESVlZgQ5XsWH37t2qVq2adu3apWOtTiIMLNNhAUfOJQzLIGRkZEQuU2D1HQ895LZLr1Sp1OfPkgoAIJYV5/3bFzUfodZOlMobb7iZjjvucG/bLhYLPsLAzpE6DgBAovBF8BGoncib+QjbdtTPPpOuuMJtHNa5s9SxY3geFwCABBT2mo+RI0cqKSkp16VVq1byUsRrJ844Q7rxRunOO91eHgAAILqZj5NOOkmzZ88OPklgu6mHwrod1cpkpk6VLrtMqljRPfbccyWaywIAgN9EJCqwYKNevXqKNWGrnbj1Vmn8eLc9ul0bAg8AALzbamsFng0aNFDTpk3Vp08fbdiwocD7Hjx40KmQzXmJeZbxKF9eatrUzYIAAADvttrOmjVLmZmZatmypbZs2aJRo0bp+++/16pVq1S1atV8a0TsPnmFc6ttqR054g6EO+GE4DFrk96gQakf2rbRWrBmRbHsaAEA+GGrbcT7fOzcudPpsfHoo486dRf5ZT7skvPkGzZsGDvBxw8/SNdea+kcafnyYnUoLYo1OgvMbrHdOFYUm99rBABAIgUfEe9wWr16dbVo0aLAEe/JycnOSea8xBQrKM3IcJuGff552B42lKFxAAAkoogHH7YEs27dOtWvX19xI2dDMlsqsgZiS5dKXbtGpfEZAACJLOzBx/Dhw7VgwQKndfmiRYt0+eWXO301rrnmGsUFe/M/80xp5szgsXbtpDD3KinO0DgAABJJ2IMPWzawQMMKTq+88krVrFlTixcvVu3atRUXnn1WWrJEGjbMLTSNEIbGAQD8yheD5YrFil9vv10aMUJq2DDiT8fQOABAIoip3S4xH3zYDpbXXpPuv59GYQAAlBBTbUNlO1jOOUfau1dq00a6/vqSvuYAACBE3g9d8ZKNvb/7bmnRIumii7w+GwAAfMF/wcfChbbVRAps/bVJtDaXhdksAAAkRp+PmNvJ0qWL1LdvsJeHbXcl8AAAIGr8FXxYfUeFCrbPVTp0yOuzAQDAl/y17NK6tbRqldSkSakGvTEMDgCAkvNX5sMUEnjYoDcbgte1a1fn2m6X5D4AAKBg9PnIkc2wYCLnvBVrd25t4gMZkFDuAwCAH+2Opam28SKUQW8MgwMAoPQIPoox6I1hcAAAlB7BRzEGvTEMDgCA0qPmowSD3hgGBwBAbsx2KQULOIoqHg3lPgAAIH8suwAAgKgi+AAAAFFF8AEAAKKK4AMAAEQVwQcAAIgqgg8AABBVBB8AACCqCD4AAEBUEXwAAICoIvgAAAAEHwAAIHEdoxiTlZWVPaAGAADEh8D7duB9PK6Cjz179jjXDRs29PpUAABACd7Hq1WrVuh9krJCCVGi6Ndff9XmzZtVtWpVJSUlhT0qs6Bm48aNOvbYY8P62OB1jjZ+n3mdEw2/0/H9Ols4YYFHgwYNVKZMmfjKfNgJR3pcvb3YBB+Rx+scHbzOvM6Jht/p+H2di8p4BLDbBQAARBXBBwAAiCpfBR/Jycm69957nWvwOsc7fp95nRMNv9P+eZ1jruAUAAAkNl9lPgAAgPcIPgAAQFQRfAAAgKgi+AAAAFHlm+Bj3Lhxaty4sSpUqKCOHTvqs88+8/qUEs7ChQvVq1cvp7uddaedPn2616eUkMaMGaMOHTo4XYDr1Kmjyy67TKtXr/b6tBLO008/rVNPPTW7EVOnTp00a9Ysr08r4T3wwAPO348hQ4Z4fSoJZeTIkc7rmvPSqlUrz87HF8HHa6+9pqFDhzpbi5YtW6a2bduqZ8+e+uGHH7w+tYSyd+9e57W1QA+Rs2DBAg0aNEiLFy/Whx9+qMOHD6tHjx7O64/wsU7L9ka4dOlS/e9//1PXrl116aWX6ssvv+RljpAlS5ZowoQJTtCH8DvppJO0ZcuW7MtHH30UgWcJjS+22lqmwz4pPvXUU9nzY6yv/eDBg3XXXXd5fXoJyaLqadOmOZ/KEVk//vijkwGxoKRz58683BFUo0YNPfzww0pLS+N1DrPMzEz97ne/0/jx4zV69Gi1a9dOY8eO5XUOY+bDstHLly9XLEj4zMehQ4ecTy7du3fPNT/Gbn/yySeenhsQDrt27cp+Y0RkHDlyRFOnTnWyS7b8gvCzbN7FF1+c6281wmvNmjXOsnjTpk3Vp08fbdiwQV6JucFy4bZ9+3bnD0fdunVzHbfb33zzjWfnBYSDZfFsbfzss8/WySefzIsaZitXrnSCjQMHDqhKlSpONq9Nmza8zmFmgZ0tiduyCyK3AvDCCy+oZcuWzpLLqFGjdO6552rVqlVO/Vi0JXzwAST6p0X74+Hl2m0isz/Ulqa27NIbb7yhfv36OctbBCDhY2Pdb7vtNqd+yTYEIDIuvPDC7H9bTY0FI40aNdK///1vT5YREz74qFWrlsqWLatt27blOm6369Wr59l5AaV166236p133nF2GVlxJMKvfPnyat68ufPv9u3bO5/MH3/8cacoEuFhy+JW/G/1HgGWrbbfa6vTO3jwoPM3HOFVvXp1tWjRQmvXrpUXyvjhj4f90ZgzZ06uVLXdZu0W8chqxC3wsCWAuXPnqkmTJl6fkm/Y3w57M0T4dOvWzVnesgxT4HL66ac7NQn2bwKPyBX4rlu3TvXr15cXEj7zYWybraVL7Rf6jDPOcCqorXDshhtu8PrUEu6XOWcUvX79euePhxVCnnDCCZ6eW6IttbzyyiuaMWOGs1a7detW53i1atVUsWJFr08vYYwYMcJJVdvv7p49e5zXfP78+Xr//fe9PrWEYr/DeeuVKleurJo1a1LHFEbDhw93+jDZUsvmzZud1hMW2F1zzTXygi+Cj6uuusrZjnjPPfc4f6htC9d77713VBEqSsd6IXTp0iVX0Gcs8LNCJ4Sv+ZVJTU3Ndfz5559X//79eZnDxJYC+vbt6xTnWWBn6+QWeJx//vm8xog7mzZtcgKNHTt2qHbt2jrnnHOcXkH2by/4os8HAACIHQlf8wEAAGILwQcAAIgqgg8AABBVBB8AACCqCD4AAEBUEXwAAICoIvgAAABRRfABAACiiuADAABEFcEHAACIKoIPAAAQVQQfAABA0fT/R7EKCJi01cMAAAAASUVORK5CYII=",
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"yH = theta0H+theta1H*x\n",
"plt.plot(x, yH, ':r', label=\"Model\")\n",
"plt.plot(x, y, '.k', label=\"Data\")\n",
"plt.legend()\n",
"plt.show()"
]
}
],
"metadata": {
"kernelspec": {
"display_name": ".venv (3.14.3)",
"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.14.3"
}
},
"nbformat": 4,
"nbformat_minor": 5
}