diff --git a/cw1/spce0038_coursework_sklearn_MCMQ7.ipynb b/cw1/spce0038_coursework_sklearn_MCMQ7.ipynb index 41b1f46..f5e928d 100644 --- a/cw1/spce0038_coursework_sklearn_MCMQ7.ipynb +++ b/cw1/spce0038_coursework_sklearn_MCMQ7.ipynb @@ -66,7 +66,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 1, "metadata": {}, "outputs": [], "source": [ @@ -139,7 +139,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 2, "metadata": { "deletable": false, "editable": false, @@ -161,7 +161,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 3, "metadata": { "deletable": false, "editable": false, @@ -230,7 +230,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 4, "metadata": { "deletable": false, "editable": false, @@ -251,7 +251,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 5, "metadata": { "deletable": false, "editable": false, @@ -293,7 +293,7 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 8, "metadata": { "deletable": false, "nbgrader": { @@ -307,12 +307,25 @@ "solution": true } }, - "outputs": [], + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjMAAAHHCAYAAABKudlQAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjAsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvlHJYcgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAWKlJREFUeJzt3Xl4FFX6NuCnk6YTSEjCjghZ2BIim4BA2JdgJMiIMgMyDIuso6iAwwzJjBA2DYwooIOK+ckyiLKKgxoMyCJrWAJRJBBQCARZMopJgAghdH1/5Oue9F7VW1V1P/d19QVdXak+Xd3pevOe95yjEQRBABEREZFKBcjdACIiIiJXMJghIiIiVWMwQ0RERKrGYIaIiIhUjcEMERERqRqDGSIiIlI1BjNERESkagxmiIiISNUYzBAREZGqMZghsmHOnDnQaDRyN8OvFRQUQKPRYPXq1U79vEajwZw5c9zaJl+1evVqaDQaFBQUiP4Zw/uzePFih/ta+32qqKjA3/72NzRp0gQBAQEYMmSIxFYTVWIwQ37B8EVtuAUHB6NRo0ZISkrC22+/jVu3brnlea5evYo5c+YgNzfXLccj8mUrV67EG2+8gd///vdYs2YNpk+fjry8PMyZM0dSUEXEYIb8yrx587B27Vq89957eOmllwAA06ZNQ5s2bfDdd9+Z7Pvqq6/it99+k3T8q1evYu7cuQxmiMxY+33avXs3Hn74YSxZsgSjRo1C7969kZeXh7lz5zKYIUm0cjeAyJsGDhyITp06Ge+npqZi9+7dePLJJ/G73/0OZ86cQfXq1QEAWq0WWi1/RUjZKioqoNfrodPp5G6KXdZ+n4qKihARESFPg8inMDNDfq9fv36YNWsWLl26hI8++si43Vof/86dO9GjRw9EREQgNDQUsbGx+Pvf/w4A2Lt3Lx577DEAwHPPPWfs0jLUe+zfvx9/+MMfEBkZiaCgIDRp0gTTp0+3+Gt17NixCA0NxU8//YQhQ4YgNDQU9erVw4wZM/DgwQOTffV6PZYtW4Y2bdogODgY9erVwxNPPIHjx4+b7PfRRx+hY8eOqF69OmrXro1nn30WhYWFDs+N4RycO3cOf/rTnxAeHo569eph1qxZEAQBhYWFeOqppxAWFoaGDRvizTfftDhGUVERxo8fjwYNGiA4OBjt2rXDmjVrLPYrLi7G2LFjER4ejoiICIwZMwbFxcUW+/Xp0wd9+vSx2D527FhER0fbfT229pH6XtvSunVr9O3b12K7Xq/Hww8/jN///vfGbevXr0fHjh1Rs2ZNhIWFoU2bNli2bJnd41etUVm6dCmaNWuGoKAg5OXlAQDOnj2L3//+96hduzaCg4PRqVMnbNu2zeI4p0+fRr9+/VC9enU0btwYCxYsgF6vt9jv+PHjSEpKQt26dVG9enXExMRg3LhxVtv2wQcfGNvz2GOP4dixYyaPVz3HhtexZ88enD592uR35Q9/+AMAoG/fvsbte/futXteiPhnJxGAUaNG4e9//zt27NiBiRMnWt3n9OnTePLJJ9G2bVvMmzcPQUFB+OGHH3Dw4EEAQKtWrTBv3jzMnj0bkyZNQs+ePQEA3bp1AwBs2rQJZWVleP7551GnTh0cPXoU77zzDq5cuYJNmzaZPNeDBw+QlJSELl26YPHixfj666/x5ptvolmzZnj++eeN+40fPx6rV6/GwIEDMWHCBFRUVGD//v3Izs42ZqBee+01zJo1C8OGDcOECRPw3//+F++88w569eqFkydPivrLePjw4WjVqhUWLlyIL7/8EgsWLEDt2rWxYsUK9OvXD4sWLcK6deswY8YMPPbYY+jVqxcA4LfffkOfPn3www8/4MUXX0RMTAw2bdqEsWPHori4GFOnTgUACIKAp556CgcOHMCf//xntGrVClu3bsWYMWMkvIvu4+i9tmX48OGYM2cOrl+/joYNGxq3HzhwAFevXsWzzz4LoDJQGjFiBPr3749FixYBAM6cOYODBw8az4k9q1atwt27dzFp0iQEBQWhdu3aOH36NLp3746HH34YKSkpCAkJwcaNGzFkyBBs2bIFTz/9NADg+vXr6Nu3LyoqKoz7ffDBB8aMpEFRUREef/xx1KtXDykpKYiIiEBBQQE+/fRTi/Z8/PHHuHXrFiZPngyNRoN//vOfeOaZZ3DhwgVUq1bNYv969eph7dq1eO2113D79m2kp6cDAFq0aIGXX34Zb7/9Nv7+97+jVatWAGD8l8gmgcgPrFq1SgAgHDt2zOY+4eHhwqOPPmq8n5aWJlT9FVmyZIkAQPjvf/9r8xjHjh0TAAirVq2yeKysrMxiW3p6uqDRaIRLly4Zt40ZM0YAIMybN89k30cffVTo2LGj8f7u3bsFAMLLL79scVy9Xi8IgiAUFBQIgYGBwmuvvWby+KlTpwStVmux3ZzhHEyaNMm4raKiQmjcuLGg0WiEhQsXGrf/+uuvQvXq1YUxY8YYty1dulQAIHz00UfGbeXl5UJCQoIQGhoqlJaWCoIgCJ999pkAQPjnP/9p8jw9e/a0OJ+9e/cWevfubdHWMWPGCFFRUSbbAAhpaWl296n6Og3EvNfW5OfnCwCEd955x2T7Cy+8IISGhho/A1OnThXCwsKEiooKSce/ePGiAEAICwsTioqKTB7r37+/0KZNG+Hu3bvGbXq9XujWrZvQokUL47Zp06YJAIQjR44YtxUVFQnh4eECAOHixYuCIAjC1q1bHf7OGNpTp04d4ebNm8bt//nPfwQAwueff27cZn6OBaHyvXzkkUdMtm3atEkAIOzZs8fxCSH6/9jNRPT/hYaG2h3VZMhg/Oc//7Gaknek6l++d+7cwc8//4xu3bpBEAScPHnSYv8///nPJvd79uyJCxcuGO9v2bIFGo0GaWlpFj9rSOd/+umn0Ov1GDZsGH7++WfjrWHDhmjRogX27Nkjqu0TJkww/j8wMBCdOnWCIAgYP368cXtERARiY2NN2piZmYmGDRtixIgRxm3VqlXDyy+/jNu3b+Obb74x7qfVak2yToGBgcYibW9z9r1u2bIl2rdvjw0bNhi3PXjwAJs3b8bgwYONn4GIiAjcuXMHO3fudKp9Q4cORb169Yz3b968id27d2PYsGG4deuW8X3+5ZdfkJSUhPPnz+Onn34CUHmuu3btis6dOxt/vl69ehg5cqTJcxjOwRdffIH79+/bbc/w4cNRq1Yt431DVrLqZ4HIkxjMEP1/t2/fRs2aNW0+Pnz4cHTv3h0TJkxAgwYN8Oyzz2Ljxo2iL3aXL1/G2LFjUbt2bWMdTO/evQEAJSUlJvsa6l+qqlWrFn799Vfj/R9//BGNGjVC7dq1bT7n+fPnIQgCWrRogXr16pnczpw5g6KiIlFtj4yMNLkfHh6O4OBg1K1b12J71TZeunQJLVq0QECA6VeNodvg0qVLxn8feughhIaGmuwXGxsrqn3u5sp7PXz4cBw8eNAYPOzduxdFRUUYPny4cZ8XXngBLVu2xMCBA9G4cWOMGzcOX331lej2xcTEmNz/4YcfIAgCZs2aZfE+G4Jdw3tteE/MmZ/r3r17Y+jQoZg7dy7q1q2Lp556CqtWrcK9e/csftb882EIbKp+Fog8iTUzRACuXLmCkpISNG/e3OY+1atXx759+7Bnzx58+eWX+Oqrr7Bhwwb069cPO3bsQGBgoM2fffDgAQYMGICbN29i5syZiIuLQ0hICH766SeMHTvW4iJp71hS6PV6aDQabN++3eoxzYMHW6z9rK02CoIgrZESaTQaq89hXhxt62etMf9ZV97r4cOHIzU1FZs2bcK0adOwceNGhIeH44knnjDuU79+feTm5iIrKwvbt2/H9u3bsWrVKowePdpqcbQ58/oWw+dnxowZSEpKsvoz9j7b1mg0GmzevBnZ2dn4/PPPkZWVhXHjxuHNN99Edna2yWdHrs8CkQGDGSIAa9euBQCbFwKDgIAA9O/fH/3798dbb72F119/Hf/4xz+wZ88eJCYm2rxYnjp1CufOncOaNWswevRo43ZnuxkAoFmzZsjKysLNmzdtZmeaNWsGQRAQExODli1bOv1czoqKisJ3330HvV5vkp05e/as8XHDv7t27cLt27dNLpL5+fkWx6xVq5bV7gtDlseeWrVqWR0hZe1nHb3XtsTExKBz587YsGEDXnzxRXz66acYMmQIgoKCTPbT6XQYPHgwBg8eDL1ejxdeeAErVqzArFmzJAceTZs2BVDZhWevbUDluT5//rzFdmvnGgC6du2Krl274rXXXsPHH3+MkSNHYv369SZdj+7EWbfJGexmIr+3e/duzJ8/HzExMRZ1A1XdvHnTYlv79u0BwJh6DwkJAQCLC6bhL9eqf6kKguBwKK49Q4cOhSAImDt3rsVjhud55plnEBgYiLlz51r8lSwIAn755Renn1+M5ORkXL9+3aSGpKKiAu+88w5CQ0ON3WzJycmoqKjAe++9Z9zvwYMHeOeddyyO2axZM5w9exb//e9/jdu+/fZbhyONDD9bUlJiMkHitWvXsHXrVpP9xLzX9gwfPhzZ2dlYuXIlfv75Z5MuJgAW5z0gIABt27YVfXxz9evXR58+fbBixQpcu3bN4vGq5yo5ORnZ2dk4evSoyePr1q0z+Zlff/3V4jMj5Rw4y9bvEJE9zMyQX9m+fTvOnj2LiooK3LhxA7t378bOnTsRFRWFbdu2ITg42ObPzps3D/v27cOgQYMQFRWFoqIivPvuu2jcuDF69OgBoPJiGRERgffffx81a9ZESEgIunTpgri4ODRr1gwzZszATz/9hLCwMGzZssWlmoK+ffti1KhRePvtt3H+/Hk88cQT0Ov12L9/P/r27YsXX3wRzZo1w4IFC5CamoqCggIMGTIENWvWxMWLF7F161ZMmjQJM2bMcLoNjkyaNAkrVqzA2LFjkZOTg+joaGzevBkHDx7E0qVLjTVKgwcPRvfu3ZGSkoKCggLEx8fj008/taglAoBx48bhrbfeQlJSEsaPH4+ioiK8//77eOSRR1BaWmq3Pc8++yxmzpyJp59+Gi+//DLKysrw3nvvoWXLljhx4oRxPzHvtT3Dhg3DjBkzMGPGDNSuXdsiWzJhwgTcvHkT/fr1Q+PGjXHp0iW88847aN++vdPDkJcvX44ePXqgTZs2mDhxIpo2bYobN27g8OHDuHLlCr799lsAwN/+9jesXbsWTzzxBKZOnWocmm3IohmsWbMG7777Lp5++mk0a9YMt27dQkZGBsLCwpCcnOxUG8Vo3749AgMDsWjRIpSUlCAoKAj9+vVD/fr1Pfac5ANkGEFF5HWGodmGm06nExo2bCgMGDBAWLZsmXGIcFXmQ0l37dolPPXUU0KjRo0EnU4nNGrUSBgxYoRw7tw5k5/7z3/+I8THxwtardZkWHFeXp6QmJgohIaGCnXr1hUmTpwofPvttxZDj8eMGSOEhIQ4bI8gVA5ffuONN4S4uDhBp9MJ9erVEwYOHCjk5OSY7LdlyxahR48eQkhIiBASEiLExcUJU6ZMEfLz8+2eN8Nzmg9RttVGa0Ntb9y4ITz33HNC3bp1BZ1OJ7Rp08bq0PVffvlFGDVqlBAWFiaEh4cLo0aNEk6ePGl1qPtHH30kNG3aVNDpdEL79u2FrKwsUUOzBUEQduzYIbRu3VrQ6XRCbGys8NFHHzn9XtvTvXt3AYAwYcIEi8c2b94sPP7440L9+vUFnU4nREZGCpMnTxauXbtm95iGodBvvPGG1cd//PFHYfTo0ULDhg2FatWqCQ8//LDw5JNPCps3bzbZ77vvvhN69+4tBAcHCw8//LAwf/584cMPPzQZmn3ixAlhxIgRQmRkpBAUFCTUr19fePLJJ4Xjx4+Lao/5uRc7NFsQBCEjI0No2rSpEBgYyGHaJIpGEFihRUREROrFmhkiIiJSNQYzREREpGoMZoiIiEjVZA9mDCupVr3FxcUZH7979y6mTJmCOnXqIDQ0FEOHDsWNGzdkbDEREREpiezBDAA88sgjuHbtmvF24MAB42PTp0/H559/jk2bNuGbb77B1atX8cwzz8jYWiIiIlISRcwzo9Vq0bBhQ4vtJSUl+PDDD/Hxxx+jX79+AIBVq1ahVatWyM7ORteuXb3dVCIiIlIYRQQz58+fR6NGjRAcHIyEhASkp6cjMjISOTk5uH//vsmEU3FxcYiMjMThw4dtBjP37t0zmaFSr9fj5s2bqFOnDqfKJiIiUglBEHDr1i00atTIYsHaqmQPZrp06YLVq1cjNjYW165dw9y5c9GzZ098//33uH79OnQ6nXEpeoMGDRrg+vXrNo+Znp5udYp3IiIiUp/CwkI0btzY5uOyBzMDBw40/r9t27bo0qULoqKisHHjRouVYcVKTU3FK6+8YrxfUlKCyMhIFBYWIiwszOU2ExERkeeVlpaiSZMmxqVPbJE9mDEXERGBli1b4ocffsCAAQNQXl6O4uJik+zMjRs3rNbYGAQFBVmsUAsAYWFhDGaIiIhUxlGJiCJGM1V1+/Zt/Pjjj3jooYfQsWNHVKtWDbt27TI+np+fj8uXLyMhIUHGVhIREZFSyJ6ZmTFjBgYPHoyoqChcvXoVaWlpCAwMxIgRIxAeHo7x48fjlVdeQe3atREWFoaXXnoJCQkJHMlEREREABQQzFy5cgUjRozAL7/8gnr16qFHjx7Izs5GvXr1AABLlixBQEAAhg4dinv37iEpKQnvvvuuzK0mIiIipfCLVbNLS0sRHh6OkpIS1swQERGphNjrt+JqZoiIiIikYDBDREREqsZghoiIiFSNwQwRERGpGoMZIiIiUjUGM0RERKRqDGaIiIhI1RjMEBERkaoxmCEiIiJVYzBDREREqsZghoiIiFSNwQwRERGpGoMZIiIiUjUGM0RERKRqDGaIiIhI1RjMEBERkaoxmCEiIiJVYzBDREREqsZghoiIiFSNwQwRERGpGoMZIiIiUjUGM0RERKRqDGaIiIhI1RjMEBERkaoxmCEiIiJVYzBDREREqsZghoiIiFSNwQwRERGpGoMZIiIiUjUGM0RERKRqDGaIiIhI1RjMEBERkaoxmCEiIiJVYzBDREREqsZghoiIiFSNwQwRERGpGoMZIiIiUjUGM0RERKRqDGaIiIhI1RjMEBERkaoxmCEiIiJVYzBDREREqqaoYGbhwoXQaDSYNm2acdv169cxatQoNGzYECEhIejQoQO2bNkiXyOJiIhIURQTzBw7dgwrVqxA27ZtTbaPHj0a+fn52LZtG06dOoVnnnkGw4YNw8mTJ2VqKRERESmJIoKZ27dvY+TIkcjIyECtWrVMHjt06BBeeukldO7cGU2bNsWrr76KiIgI5OTkyNRaIiIiUhJFBDNTpkzBoEGDkJiYaPFYt27dsGHDBty8eRN6vR7r16/H3bt30adPH5vHu3fvHkpLS01uRERE5Ju0cjdg/fr1OHHiBI4dO2b18Y0bN2L48OGoU6cOtFotatSoga1bt6J58+Y2j5meno65c+d6qslERESkILJmZgoLCzF16lSsW7cOwcHBVveZNWsWiouL8fXXX+P48eN45ZVXMGzYMJw6dcrmcVNTU1FSUmK8FRYWeuolEBERkcw0giAIcj35Z599hqeffhqBgYHGbQ8ePIBGo0FAQADy8/PRvHlzfP/993jkkUeM+yQmJqJ58+Z4//33RT1PaWkpwsPDUVJSgrCwMLe/DiIiInI/sddvWbuZ+vfvb5Fhee655xAXF4eZM2eirKwMABAQYJpACgwMhF6v91o7iYiISLlkDWZq1qyJ1q1bm2wLCQlBnTp10Lp1a9y/fx/NmzfH5MmTsXjxYtSpUwefffYZdu7ciS+++EKmVhMREZGSKGI0ky3VqlVDZmYm6tWrh8GDB6Nt27b497//jTVr1iA5OVnu5hEREZECyFoz4y2smSEiInJNWXkF4mdnAQDy5iWhhs7znTtir9+KzswQEREROcJghoiIiFSNwQwRERGpGoMZIiIiUjUGM0RERKRqDGaIiIhI1RjMEBERkaoxmCEiIiJVYzBDREREqsZghoiIiFSNwQwRERGpGoMZIiIiUjUGM0RERKRqDGaIiIhI1RjMEBERkaoxmCEiIiJVYzBDRETkBWXlFYhO+RLRKV+irLxC7ub4FAYzREREpGoMZoiIiEjVGMwQERGRqjGYISIiIofKK/TG/685WGByX24MZoiIiMiu9Mw8dJi/03h/UVY+4mZtR3pmnoyt+h+t3A0gIiIi5UrPzMOKfRcttusFGLenJsd7u1kmmJkhIiJSACUO3S6v0CNjv2UgU1XG/ouydzkxmCEiIiKr1h4ugF6wv49eqNxPTgxmiEgxlPiXKZE/u3SzzK37eQqDGSIiBWOAR3KKql3Drft5CoMZIiIismpUQjQCNPb3CdBU7icnBjNERERklU4bgIk9Y+zuM7FnDHRaecMJBjNEROTX2JVnX2pyPCb3irHI0ARogMm9YmQflg1wnhkiIvIxZeUViJ+dBQDIm5eEGjpe6lyVmhyP5/s0R/t5lRPnzUyKxfieTWXPyBgooxVERESkaFUDlzHdoxUTyAAMZoiIiEjlGMwQERGRqjGYISIiIlVjMENERESqxmCGiMgHcbgx+RMGM0RERKRqDGaIiIi8oLxCb/z/moMFJvfJNQxmiIiIPCw9Mw8d5u803l+UlY+4WduRnpknY6t8B6dFJCIi8qD0zDys2HfRYrtegHG7EpYEUDNFZWYWLlwIjUaDadOmmWw/fPgw+vXrh5CQEISFhaFXr1747bff5GkkERGRSOUVemTstwxkqsrYf5FdTi5STDBz7NgxrFixAm3btjXZfvjwYTzxxBN4/PHHcfToURw7dgwvvvgiAgIU03QiIiKr1h4ugF6wv49eqNyPnKeIbqbbt29j5MiRyMjIwIIFC0wemz59Ol5++WWkpKQYt8XGxnq7iURERJJdulnmlv24eKZ9ikhvTJkyBYMGDUJiYqLJ9qKiIhw5cgT169dHt27d0KBBA/Tu3RsHDhyQqaVERETiRdWuIXo/jnZynuzBzPr163HixAmkp6dbPHbhwgUAwJw5czBx4kR89dVX6NChA/r374/z58/bPOa9e/dQWlpqciMiIvK2UQnRCNDY3ydAA1wrvsvRTi6QNZgpLCzE1KlTsW7dOgQHB1s8rtdXRqWTJ0/Gc889h0cffRRLlixBbGwsVq5cafO46enpCA8PN96aNGnisddARJacnX2Wf5mSr9FpAzCxZ4zdfVo/HIb/O3jRorbGMNqJAY1jsgYzOTk5KCoqQocOHaDVaqHVavHNN9/g7bffhlarRYMGDQAA8fGmQ9ZatWqFy5cv2zxuamoqSkpKjLfCwkKPvg4ich3n4bCOAZ76pSbHY3KvGIsMTYAGGN8jGt//ZL/3gKOdHJM1mOnfvz9OnTqF3Nxc461Tp04YOXIkcnNz0bRpUzRq1Aj5+fkmP3fu3DlERUXZPG5QUBDCwsJMbkSkXIZ5OPiXqSkGeL4jNTkeJ2YNMN6fmRSLs/MHolF4dVGjnT45YvsPeJJ5NFPNmjXRunVrk20hISGoU6eOcftf//pXpKWloV27dmjfvj3WrFmDs2fPYvPmzXI0mYjcTOw8HH95PA46reXfX746yoMTrfmeqp/fMd2jodMGiB7tdPlXcfv5K9kLgB2ZNm0aUlNTMX36dLRr1w67du3Czp070axZM7mbRkRuwHk4LKlxojU1r9ItZ1ee2NFOkbXE7eevFBfM7N27F0uXLjXZlpKSgsLCQty5cweHDh1Cjx495GkcEbmdu+bh8CUM8LxH7q48saOdRnSJ9Ep71EpxwQwR+Rcp83D4CwZ4rhGbaVFCrZaY0U4Te8ZY7WL1tho6LQoWDkLBwkGK686V/+wQkV8T+5fpqIRor7RHCRjgOU9spkVJXXn2RjtN7hXD2igRGMwQkazU9JeptzDAc46UTIvSuvJsjXZiICOO/3w7EJFDchVxpibHY1x3y+kW/PUvUwZ40knNtCixK8/aaCcSh2eKiBRhRlKcyX1//8uUXQ/SSM20sCvPtzCYISJF4l+mll0PAHBi1gBRgYy/zRwsNdPCrjzf4t/fFERECmce0IkJ8OQebiwHqZkWduX5Fr5LREQ+RAnDjeNnZ3l94jxnMi3syvMdDGaIiHyEkoYbe5uzmRaOInJMDbM7M5ghIvIRShtu7G3OZlo4ikj9lDWFHxEROU2Jw42d4crioanJ8Xi+T3O0n1dZMzQzKRbjezZlgOLj+O4SkSKJTWdX3U+pKXBv4XDjSsy0+B++w0REKmKvfoHDjX2Xvw21l4rBDBGRj1DqcGM1FJAqmT8OtZeKwQwRyUbMRU5JF0IltcUWDjf2LUoYaq8GDGaIFEINF0qxmBKXF4cb+wZ/HmovFYMZInIrpsSVgUWw6ufvQ+2l4NBsInIbQ0rcnCElDsDt2YGqf5Wuy76M5/s054Xbi6oOoyb7aui0KFg4SPT+vjLU3hv4G0+kEGofYixHSjw9Mw89Fu0x3l/y9XlmgchncKi9eAxmiMgtXE2Jmwc5joIeFkYqG+s4XMeh9uIxmCEit3A1JV41w2K4bysgkZIF8qXCajXpMH8nA0oXKXWovRKxZoaI3MLVlLitDMv9B5Z/4UvJAo3oEimqXeScxVlnrW6vWic1NbGlN5vkUww1Zhn7TbOQAZrKQIYj1CoxnCPyMH/JDDiTEhfTFbH60CWLbSyMVIbyCr3V96cqDh12HYfaO8ZghoiMXJkfxpmU+CdHLjs8rrUMDAsjlUFshkzM++wuah5dZRjtVLBwkMXimmodau+tP+bUcTaIyOPcMT9ManI8xnWPsvrYuO5RFn9JXv7VucyJPxdGrsv2XmDgiNjMl7PvM5FYDGZIlfyl68Zb3DkyaEZSnOjtkbWcy5z4U2GkeU3Kkq/Py9QSS2IzX86+z0Riqf83nYhcIueU6WKKc21lYGytQQRYzwKpUXpmHlYetF2TYqv41pyngv+hHRs73CdAI+59JnIFgxkiPyfnlOliMidju1nvtgIqR8lYa7ut7JCaiAkyVx+6JGtxrZj3z1cyZKRs/IQR+TmljAyytcqzs4GJIRuh1mJQta/LowEwoTuHDpN3cJ4ZIj+nlJFBO6f3Qv+39hnvn5g1ABE1dH5bE6WUINNZAoCVhy4iMFD588xIXTOJlIeZGVIEFvTKRykjg8y7Ivy9a0IpQaYrDAXkYmt7SJlcmbLBW/z724KIFDsyyPClqfYFOJ2llCDTHRxNrGfOXRfPvHlJFvO1kDTumLLBGxjMEPkwsRkvW90AhroVb9Q9vLv3B5P7hi/NZbvOefy5lUinDbBb/AxUFkerIYPlqPanKrVcPP2BmhZzVf5vARHJwtkp053Nnnxy9IrFtsoC10KbP2Or+8JXujUcFT/XCw1WZMrfWWq6ePo6OadscAaDGSKySulTptubg2XlwUs+E9DYo8asha1uI7VdPH2d2kbTKfebisjPVP2SXpd92S++tKW+RsN5ETsHizuZr/mjlEJItWUtbAVgart4+jq1jaZjZRSRAqRn5plcnJd8fR7Ldp3HxJ6+O0+H+WsWw3BeOsfUFnXhcwdbCxcuysrHGzvyFfMeZey/iL88HmfMplUdbqy0wmlDAAbAeO7UdvH0dWobTcfMDJHM/LFOwNZrFkMvANkXbrq/UU5Q0nukxqxF1W4jtV08fZ3aRtMxmCGSkT/WCYh5zUoi5twr5T1SUtaiX2w9h/tUDcDUdvH0dUqdssEWZbSCyE+puU7A1rBvWxd1Q42JmNfsDo4ujGJ9cuSyw32U8h55O2thr8i6QXiwqGMYAjC1XTz9ga3FXL05ZYNY/FQQycjX6gTSM/PQY9Eeq48ZCj+3nLAcgu0J1uZocaZo9/Kv6niPvJ21cLSi949Ft0Udp2oApqaLp79ITY7HiVkDjPednbLB0xjMEMno4Yjqbt1PTmLqYPQCkHftltues3N0Lavbx3WPsjpHizNDmSNrqaOWw5tZCzFdhccKfnWq20gtF09/UvVzpdQpGxTVooULF0Kj0WDatGkWjwmCgIEDB0Kj0eCzzz7zetuInGVvanaxPSFu6jHxGE/UwTh6zQEa4IPRnaw+Vi80GIu2n7H6mNSi3RFdIh3uI2cthxxZCzFdhQKALtF17O5jKwBTw8WTlEUxQ7OPHTuGFStWoG3btlYfX7p0KTQapX+lE5kyH35sPpz3SvFvoo4jdj+5fHLkstvrYJ7rHmW3G8NeJmJRVr7D45sPZbZFzIVUrlqOmUmxGN+zqdefW2yXWouGoWjbJAwZ+00zdgEaKGZIu1pwZW/7FBHu3r59GyNHjkRGRgZq1bJMG+fm5uLNN9/EypUrZWgdkXPEDLn2leGoYutKzNnrhpiRFOfR+gl3FO3KXcshV9ZCyufW291Galjh2Z946/1QRDAzZcoUDBo0CImJiRaPlZWV4Y9//COWL1+Ohg0bytA6UiKlf2GJHXI9/LFIjw1HtTXZmyc0DBM3csXcgZl9jf9/vpfpSJbyCj1Sk+NN9pme2MLkQujq++5K0a4/13JIHUbtrW4jLlKpLN58P2QPZtavX48TJ04gPT3d6uPTp09Ht27d8NRTT4k+5r1791BaWmpyI9+hhi8ssUOuNxy77BPDUf/QqbFTP1f1da0wC/46zN+J9Mw8k31Gdo00uS9m2LQ9rmS8/LmWQ8wwam+v6O2Pk08qmbffD1l/EwsLCzF16lSsW7cOwcGWf9lt27YNu3fvxtKlSyUdNz09HeHh4cZbkyZN3NRikptavrCkDLn2heGo7rho2XpPl+06Z/NnnO3eAvxrAjZPZDJTk+MxJsF2cbSjFb/dyR8nn1QyOd4PWYOZnJwcFBUVoUOHDtBqtdBqtfjmm2/w9ttvQ6vVYufOnfjxxx8RERFhfBwAhg4dij59+tg8bmpqKkpKSoy3wsJCL70i8iQ1fWFJrYVx1J0ihWEyO291MZmTMlmdmPdqXbbt39+Hw50fsq6GjJc7eCqTmZ6Zh7XZrmXG3EXNk0/6IjneD1l/k/v3749Tp04hNzfXeOvUqRNGjhyJ3Nxc/OMf/8B3331n8jgALFmyBKtWrbJ53KCgIISFhZncSP3U9IXlzNTs9rpT1GTblO7G/0/oEW38//TEFhb7bjrueAI9W+95emYe3txpO2tji5oyXq4G5p7KZLqytpYn+Nrkk2onx/sh69DsmjVronXr1ibbQkJCUKdOHeN2a0W/kZGRiImx319L6mKeBrc23FRNX1iGmgLDysDW+GpmoOpr+r8DBcb/j+waiSVfnzfZ19kh54aLqTPGdVNHIOPMquJVic1kihmeLvW4hv1q6EQf1iW+MirQV8jxfvjeNympjtg0uNq+sHyhFsbTGkuc2Xhd9mXcvlvh0kX+w4PK6Iq0xx2ZD09lMsWureVqcbYUXKRSWeR4PxQXzOzdu9duwa8gCBgyZIjX2kOeJSUNrsYvLCVNzW6v8FOui7vUUVBLvj6PNnOyXLrICwBGrzzi/AE8zF2zKXsqkyl2f1eKs6XiIpXKIsf7wXeWZCO1oFetX1hKmZrdVsbLPDMGVAaFy76WXo8CSAuMnDkX7ijTyL5wU7HZGXetKu6pTKbY/cWuaeUuzIQqi7ffD2V965NfcSYNzi8s15hnvDxRIPq75QedapuUUVDuoIRCcWukZEoWZ521+ZinMplijguIW9PK3ZSUCSXvvh8MZkg2zqbB/fULyzDkOjrlS5SVV9jcT0zGIWP/RVG1J84MdbcVoNqbLwYAdk7vJel5XKWEQnFrpGRKVh68ZDPg9FQmU8xxDfvJQSmZUKrkrfeD7zLJxpU0OL+wbBNTeKkXgOkbTrqtQNTZ+WI6Ldhl/L+330Oxnz9v1xmJzXwY2As4PZXJtHdcIjnwCkCyUWNBrzuIzbA4S2zh5eVfxQ2LFpPBcGW+GDmI/VxZqycyMCy34G5iMx8GjgJOT2UyzY8LmK61ReRNDGZINmot6FWysvIKfHJU3IzXkbXEDYsWk8Fwdr4Yezz5V76Yz5Wj4dGeXELDkPkQy1HA6alMpvlx+LtKcuEnj2TFgl55BGiAJcMfdSkzVrVr40bJXTe2rtLYblF2H3cU61h7XOznSsrwaE8toZGaHI+/DLCcNdkapcytVEOnRcHCQShYOAg1dLLOyUp+hsEMyc5fC3rlNLFnDEKDtU5nxsy7X/ad/9nhc0rNtMxIirMZ6HZtWtvhEG1rj5+YNUDU50rK8GhPLqEhpivMF7tiiaRiMEOKwIJe93BUg2OemXAmM+bs7LQju0pfvd7WApwtG9SUfCxAfDeIpyaSk0pMe9kVS56mhoybMltF5KPiZ2fh+Kv9ZXnumUmxVte8Sk2Ox/N9mqP9vJ129wNcm512av+WWHtYXD3P8Vf7G780rS3A6eluFU9NJOeKAI1lEfW47lE+mcE0XDyJxGI4T+Qn7GW8xGbGXJmd1tE8M1KIHQnnyeNXfR5vdPNUzVDNTIrFuQUDMXtwazs/QeQ/GMwQkWiudKdYm2fGWWJGwjkqIHb1+Abe6uZhVyyRbfxtIJ/m6Tld3KlqX7Sjfun42VmSX487Rty40p3iKKPzwf4frW63dV5Sk+Mxrrv1gGVc9yjMSIqT3sgqbNUTGXDEnXhqqLkgdWMwQz5FacFLWXkF4mdnyd0MAO6Z5E3q7LRSfFSlnmZd9mVRwZetgMXVQMbA2sRwAEfcESkNgxkiP+GOSd6kzk4rRdXEzZKvz1td4VsO1rpz2M1DpCzM95HqVc1+yDVSyB3cMYJD7CKTf3k8zumLsSEb8cG+iw7nenGFIfi6/8AzayARke+Q9G1WWOi+Aj4icr+XPjnpcB93TPJWOTttS5eOIdbqQ5e88jzkmrx5SayHIdlI+uRFRUWhdu3aaNeuHdq3b2+8lZeX4+2338aaNWs81U4iRauaHZLrSz09Mw85l4pF7euOSd5+KhG/HlObh8Nw6qdSp55HSQtUEpEySfrGvXjxIk6ePInc3FycPHkSGzduxNWrVwEAYWFhHmkgqYMSLub+TOpkdu6Y5C2ylrhjNIkIdjqQcSd+Lol8l+TMTFRUFIYMGWLcdvjwYYwZMwbz5s1zd9uI/I756CexF2Apk9m5a5K3EV0isSgr3+F+hcXuX4TSYHHWWZvb3TWiiYiUz+Vy/ISEBCxbtgyLFy92R3uIyAlSuo3cNcmbt0bz2BoKnp6Zh5UHrdfTrDx4yWagQ0S+R1Jmpry8HDqdzmJ7ixYtcPr0abc1ivyDK6OQqo78UcJ8MnIT222UEFPHI3OjWFs3yF3GdouyCFrEdKuxcJjIf0j60yo0NBTt27fHc889h2XLlmHfvn344Ycf8M477yAxMdFTbSQiWA67rnpfzGR2GgBrxnf2QMtM1w1yRkJMHZsrd1vrLhLTrcbCYf/E2Yb9k6R3evfu3fj222/x7bffYt26dUhNTcXdu5X94U888QRmz56NNm3aoE2bNoiLY381eZ6SZtj1pPTMPItMRIf5OzGxZ+V0+obJ7Fbss52teK57lMe6hlw5rgZAxphOqNDrra7cbS3z5o7RWETked5aAV1SMNOjRw/06NHDeF+v1yM/Px+5ubnIzc3F0aNHkZGRgaKiIjx48MDtjSXyBVUzKmsOFhgv2rakZ+ZZDVIMk8oBlfO+pCbHo+x+BdYetj4flFILYgUAbedmmSwM6WiGXXeMxiIi3+FSDi4gIACtWrVCq1atMGLECOP2GzduuNwwIjWpmj0oK6+wm97usWiP8f+LsvLxxo58Y4bFnJjakKoz+j7fu7nNYEbJ9AJM6mKqDvO3ZlRCNF7LPGO3K8mTdTxq46m/jr31VzeRIx7pUGzQoIEnDks+zJ++FM0vsFUzLFMTTWfV/eTIZVG1IWsPF6Do1l1Jc82omZhutYk9YzA1saVfdEMS+TtWR5FP6bRgl1efz501Oxn7L+L5Ps1Ntl3+VVxtyJYTV5B37ZaofT09weGIzo3xydErJtsCNECnqFo4WvCr257HkMnK2H/RJOAL0MCY6VL7SDeO2iMSh8u+EimEXqjMxFQldpbdMyICGVuLULq7iPoFs4BsZlIszs4fiGb1Q932HAapyfEmI6mmJ7bA2fkDjYEOR7YQ+QcGM+SznP1Ltqy8AtEpXyI65Uuv/zVsnokZ0SXS4ZBrAKJWr96Sc8XhPp7I1BiKecUGZuYcBSRVC4VHdo302mR+RKQc/K0nvyZn4GKN+QXfUBtizyONaoo6ttzDmUd0iZT1+YnIdzGYIbIjfnaWS0GOlCn1AzTWL/ipyfGY3CvG5qRyzzzaWNTxvTWcubou0Op2Zkxcwy4zItv47UJeIVcGxHxOF1t1I556bilT6usFGCeNM5eaHI8TswaYbHslsSX+8ngchj8mLuMhdj85BGiAcd2jbD6utAwaESkLgxkV8tUvdne/rslrc/BoleBgUVY+4mZtR3pmnsvHFkPsStYa2L+QG7y965zJ/cU7zyFu1nZM+PcxUe3ZcOyy451kYCgQdnZSv6pZCmYsiPwTf/PJZ+VcKrbYZj5rridtOeG44BYA/tilCWYkxdlcARqwvUK0XgCyL9wU9Txy18zYYigQrij3XtZMCn+aA4lIrZiZIb+Usf+iR7uc0jPzRM/70rSu/SHLYmYBFoNLABCRr2JmhlRNSoFtVYZZcz0xwkZK8BGgqZyav0JvO7AS210l5nnciRO6EZFSMDNDbiFHoa2trhexxHa7xM/OklTHk7LllOjgY2TXJg5H+bije+jRyAiXj0FEpFQMZshl6Zl56DDfu4W27uh68VS3y77zP4ved2r/lg73EdvOztG1bD6Wc6nYa8XPNXRamwtEknv56mAAIqkYzJBL0jPzsGLfRZuLJ3rq4ulq14u7ul2c7eYS68P9FzD8McezAAdogA9Gd7K7j+E98XSbiYi8jcEMOU1MdsRThbaudr1M7Bnj8iRuUueRccabO8+j9ZwstH44zO5+Ul6Pp9tMRORtDGZ8lDfSz2KyI4ZCW2tcqbNxtovIMGuuO4Zlu6Mw18BRtuS7K6VWAxpnXo+72uyK46/2l7sJRORDFBXMLFy4EBqNBtOmTQMA3Lx5Ey+99BJiY2NRvXp1REZG4uWXX0ZJSYm8DSUA4rMj1vazVmdT9b4joxKiRS3AWHUfw+RszgYy5ksbuGveFrEZntM/lZrcd/X1uBOn2iciOSkmmDl27BhWrFiBtm3bGrddvXoVV69exeLFi/H9999j9erV+OqrrzB+/HgZW0oGYrMj5vstzjprs87GwFE2SacNcNj18miTcJMlAAyTs7mLuwqIxWZ4zHdx9+shIlIrRXwT3r59GyNHjkRGRgZq1frfiIzWrVtjy5YtGDx4MJo1a4Z+/frhtddew+eff46KClbuy01MdsRaoa2YLISjLqf0zDx8d6XU7j5vj3jU7Rf7qt13Qzs2FvX6HfH2zLxi2iQHQ3aHI6GISCpFBDNTpkzBoEGDkJiY6HDfkpIShIWFQau1ncq+d+8eSktLTW7kfjptACb2jLG7j7XCVDFZiE3HbS8FIHZYtrOFx2JrjHTaAIztZn9NJUePA96fmVdMm4iI1ET2YGb9+vU4ceIE0tPTHe77888/Y/78+Zg0aZLd/dLT0xEeHm68NWnSxF3NJTOpyfFWF0l0tdD2SvFvNh8T2y2z9eRPTj23N/1W/gDzvzwjal93JFRaPxzm9IKOrhJbVyPnSudEpE6yBjOFhYWYOnUq1q1bh+DgYLv7lpaWYtCgQYiPj8ecOXPs7puamoqSkhLjrbCw0I2tJnPmF0dDYerUxJZOj6hqHFHd5mNiu2Wyvr8u6TmlElO4685h0M85WFlbzEU/72opbt/933uhtGDB1gSMnBuHiOyRddhBTk4OioqK0KFDB+O2Bw8eYN++ffjXv/6Fe/fuITAwELdu3cITTzyBmjVrYuvWrahWrZrd4wYFBSEoKMjTzScbHK2CHKBx3NX0h06NbT72UJj9wNfg/H/vePQiOOnfx0UNTbdnwzFxgfbkXjGYmtjS7vINnxy57PA4egHo8vou4/1FWfl4Y0e+8f6agwUY37Opy7VGzqw0bZiA0ZxegN3XzVWtiUjWzEz//v1x6tQp5ObmGm+dOnXCyJEjkZubi8DAQJSWluLxxx+HTqfDtm3bHGZwSPnE1GzYu5hKWRzSUWbEWpeG2EzF0YJfRbfDln/t+dHu438Z0ALnFogbfn3xl9uintM8vqoacHljKQpr3LUyOBH5J1kzMzVr1kTr1q1NtoWEhKBOnTpo3bq1MZApKyvDRx99ZFLMW69ePQQGBsrRbHLRjKQ4VAsMQMZ+0+HZYjI2gP1Ax5y946Vn5plcQBdl5eOfWfkWF3s5ic2SpGfmYdNx99QIGZY9AOC1OWzcOQEhEfkf2QuA7Tlx4gSOHDmCU6dOoXnz5njooYeMN9bBqFtqcrzJHDAzk2JxYGZfrz2/rTWl3H099cYwaMNrcXfbPbUUhTXeHp5ORL5FccHM3r17sXTpUgBAnz59IAiC1Vt0dLSs7ZSTGkd7WGtz1YzDmO7RJvuvy77s9tdlOJ4rXRpS2+TOYdDWVqP2ZPeMvaUo3E3K8HS1fOaJyHsUF8yQfbZGexhqHAyTusXPzpKriVb1WLTH+P9FWflo+ep2pG//X13G4qyzJvss+fq8qNoNKZkPQ4GsK10aL39yUvS+CTF1PD4M2tPdM97KmIhdngKQr66HiJSLwYyK2OoaMdQ4KPnL3doFd+3h/3UVrjx4yanXFdcwVHQbLv9aeWF25QJ9/FKxqP00ANaM7+z08xg4Go3l6WDDUxP6mc85I2YCxqrU8JknIu9hMKMSYroTvFnjIIa72rJi30XjXDVl5RUmWaez18WN4AGAyFo1UFZegX8fdt/cL7ZM6lU587GrQ8NXH7pk9zx6cvZga0tReFJqcjwm94qRlG1T2meeiOTBYEYlxHQn6AVxc414izfaIqWLRcqQbjHGdY+yeuEdldAEqcnxSM/Mszs/ihiO6lakdM8AtttsjbWlKDzNvDDcEW/W9RCRcjGYUQmx3QmGrhQlUFJbAGlDusWYkRRn9cI7tX9Ltxbm2nvvpXbPmLe5c3Qti2USXF2KwlVS3yeOhCIiWeeZIfHEdidE1vLuooWAZXdSeYUeNXTytEWKcd2jsPqQZa2OFLYuvO4szHX03huCDvN5e2yp2ubV4zqjvEKP9vMqi8pnJsW6ZQZgb/L2Qp1EpDzq+cbyc2K6EwI01rtSPDmU1Xx0FQB0mL8T6Zl5bu/WsUZKF4v5OTDPUsxMisWYBPe02V3ZArF1K+bdM9MTW4h+DvMh8koKZMR85r1Z16M0apymgcgTlPOtRXaJ6U6wVePgqaGsjkZXvb3rnFufz5qRXcWviG6thsf8Qj5zYCvRx/NGYe7YblGig4uq+43s6vlA0hsczdMjR12PUjiapoHIn/jnt4BK2RrtUbXGwdbomapDWQ1z0TizmrWBmJqQqusiOfoL21phaoCmcrs9U/u3tH/gKq6V3hW9rxj2CpzFFuZ+9XJPu497ep4apZuRFOfwM++P1DxNA5EnMJhRGWvLAJydX7kQYXmF3uHCiu4ayip2dJWBo6UKZiTFYc9f+hjv92pRF9+lJVlczJ0NvgD311bYK3AWW5gbEixtfTHz+Vn8gb3PvD9S4zQNRJ7GYEaFbNU4eHP4ttSaEEddAYuzzqLvm3uN9/ed/xlt52a5PE+LgSdqKxwVODszb4o5XpAqKbmux9vE/p5zyDr5E//9RvBBnhy+bd415e4sh60ZgF2dp8XAE7UVu87ccBhsuFKYC/yvmJrIQOzvOYeskz9hMONDvDl8W+zoKrl5srbi+OViPDpvp8P9XCnMZQ0EmRP7e84h6+RPGMz4EFeGb0slpiZkVJUL97ps789MPD2xheTaCqndWq5OJSO2G8ldNRAcvqt+Yn/P/XnIOvkfBjM+RKcNkDyU1ZWRTanJ8VZHGwVogLaNw7C2SgCz5Ovzko7tDiO7RkrqWlqcddYt3VpSCnMHLNknaj9XaiCqXvgWZeVbzAtE6uLKNA1Evoqfdh9jayivp7pbzJ9vZlIsxnWLwXdXSt02A6497sz4OBoJ5glSzpGzNRDWapFI3cRM00DkT/xjbKef8+YU9SO6RHr0L/9lZhPxuTPjo/SLvJQaCLFdSeUVev4Fr1KpyfF4vk9zVS9FQeQu/NT7kLLyCsTPzrLY7s2hrJ8cueyxoGBx1lmsPVwoen/zwEfNpNZAbDp+RdR+SlplnaTjkHWiSvzkk1s5u1K2mBmApXYDrcsu9JliV6k1EFeKfxO1n9JWNicicgaDGR9iaySOuyaeE8PZYd8zkuJMZgk2jESqWpMjNeMjtWhWCUPJzTlbA9E4orqo/ZS+sjkRkRgMZnxEemaezZE4Kw9ewrzPv/dKO0Z0iXQ6KDCfj8UdKXMpRbOORoIBlgGP1NcqpetrxoCWTk/b/4dOjUXt542VzYmIPI3BjA8Qs1bLqoOXrHa5dFqwy61tETM83JukFM3OSIrDmAT7F/eq2aOZSbEms/vmzUtyuGaSlJqfcS4Mr42oocPkXo7XhmKNBRH5An6T+QAxa7UIAEavPOKV9jiz0rOtrrCqAZjUhI+1otmqxzOfQG5x1lmTuXGsUVPBpb3hu0REvkS538QkmtiulOwLNxVbELv6kGXmKD0zz2SYt9RBUiO7NjEJNsyPZz6BnLX1odTO2orTVe8rkT+uDE5EruE3hQ+Q0pWy9nCBIusk9ILpcOJlu85J6pKxZmr/lsb/p2fmYcU+y644V4MXw4VXycyzSeZtljrzMxGR0jCY8QGjEqKx4MszojIXSl5Jt+pw4nXZrgUywP+WFRBTU+TP1BCQERHZw2DGB+i0AegaUweHL/7icF8lr6R74b+3jf93Z3ePmJoiX8CghIj8FYMZFbJ20VozvjNavrrd7s8ZCmKldCvYmlVYqgCN4wAl51Kxy89jjZKzUURE5DoWAPsInTbA6grWVU3sGYM3d5yVZdVkOYdrKykbZW2mYyIicg2DGR9ib0i0IdBZse+izQyJ2JmCy8orEJ3yJaJTvhSd5ZmRFIfO0bVE7euK46/2t9g2KiFaMQHEjKQ4k9FE0xNbGP+fNy8JefOS5GgWEZGqMZjxExN6NnVYBGtteLQ16xzMxWLLD0W3He/kATptACb2dDyBnK31oTzRHoORXZU3skyJOFybiOxhMOMnNh2/4rBmRexaRku+Pi/5+W/frcDNsvuSf85dwYSYCeTMsyaA6Yy/zuKFmIjIsxjMKIwzXThiiF1F2bxY1lHX01s7TNcaspXZSdnynajnN+eOYMJAzARy5jP6KnmGXyIiqsQ/E/2E2FWUq+5XXqHH6kPWF680+Pio6XwwHebvtFrse+JysajnNyc1mHA0PNl8AjkiIlI//tnpJ/7QqbGoLpv07WeRnpkHwLn5WfQCrK7eXfyb9C4mJaih07Iol4hI4ZiZ8SH2incNRbDWpvSvSgCM+/x2X5nrOBEREVXFzIyPMF9E0RpbRbDWZOy/KLprioiISE4MZnyAYRFFR11CZeUVovYDKruLBHhmaLKzbA2ddjRZIJGv4kg5okoMZlRO7CKKYuaPMfdT8W+yztxrbkZSnMXophOzBtidLNAdeMEgIlI2BjMqJ7ZIN2XLKckBTVTtGh4PFKyxlw3i0GkiIjLHK4HKiV1Ecd/5nyWtyWRYlNKaqlPwe4KSskHexiwQEZF0/LZUKWdWs5YyzHpizxjotAGoKLfM5ozsGunULMCOBGgqn3dqYkurw7uVwNE8NkRE5H2KyswsXLgQGo0G06ZNM267e/cupkyZgjp16iA0NBRDhw7FjRs35GukArmzSFcDYHKvGKQmx7vvoCLMTIrF2fkDvf689uTNS2J2hIhIBRTzTX3s2DGsWLECbdu2Ndk+ffp0fPnll9i0aRPCw8Px4osv4plnnsHBgwdlaqnyjO0W5dZMxl8ed65OJkBjP/tj7/Ex3aP9sv7FG5keZpOIyNcp4upx+/ZtjBw5EhkZGahVq5Zxe0lJCT788EO89dZb6NevHzp27IhVq1bh0KFDyM7OlrHFyjIjKU70/DGOCBC32KQ1jmpd3F0Lc/zV/qih07LOhIjIzykimJkyZQoGDRqExMREk+05OTm4f/++yfa4uDhERkbi8OHDNo937949lJaWmtx8nfkiiq4QW1RszlZQFaCp7LqSY2SUI1VHeK05WODUEHYiIpKX7H/Grl+/HidOnMCxY8csHrt+/Tp0Oh0iIiJMtjdo0ADXr1+3ecz09HTMnTvX3U1VPHd100TVrmH8v9SLe2pyPJ7v0xzt51WOnJqZFIvxPZtCpw1w6yrg7rA466zJQpqLsvLxxo58GVtERETOkDUzU1hYiKlTp2LdunUIDg5223FTU1NRUlJivBUWFjr+IQJgOiTb1hIJy3ads3sM85WplVoLs/LgJYsaHqkLaxIRkfxkvcrk5OSgqKgIHTp0gFarhVarxTfffIO3334bWq0WDRo0QHl5OYqLi01+7saNG2jYsKHN4wYFBSEsLMzkRuJGPRmGZNtbImHtYc8Gh1WzQeuyL3v0ucS2g4iIlEvWYKZ///44deoUcnNzjbdOnTph5MiRxv9Xq1YNu3btMv5Mfn4+Ll++jISEBBlbrnzjukfBPHYRHGQdxnWPQmpyvOglEjzBPBvkiflsxPrkiGcCKRYsExG5l6zfpDVr1kTr1q1NtoWEhKBOnTrG7ePHj8crr7yC2rVrIywsDC+99BISEhLQtWtXOZrsceYFqYZ6E6msDdV21INiKNAVu0SCucVZZzF7cGvHO9pgyAbZO/4LfZs7fXypLv/qXCE0ERF5lzKLGapYsmQJnnzySQwdOhS9evVCw4YN8emnn8rdLI8wz0osyspH3KztSM/M82o7nB3NtPrQJae7ZsRkg1YfugRtwP8+sp7OajSvF+rR4xMRkXsoLse9d+9ek/vBwcFYvnw5li9fLk+DvMRWVkIvwLjdW7PjVh3NJIVeqMzqjOgSKflnPzly2WE2SC+43vVTQ6fFuQUDETdru8MJ/mytTUVERMqi+MyMPxCTlcjYf9FrBamjEqKdnoDP2ayO2C4dd3T96LQBmNgzxu4+hkJoIiJSPn5bK4CYGhVD1sMbxFzsbXE2qxNZS9zPid3PkdTkeLsT/ClpjSgiIrKPwYwCiM1mOJv1sOZQSl+7j9u62NvjStfMiC6RDp8rQAOnurBsMZ81WYmLXRIRkWMMZhRAbDbD2ayHNZuOX3G4j9QlEsZ2i3K6a0aurh+1TPBHRES28ZtbAcTUqEjNeozrHmX3mGLnb7F2cbd1XClrL+XNS0LevCSTbWK6fjhHCxERmWMwowCeyErMSIrDgZn2u5Kc9exjja1uX5x11uVjs+uHiIikYjCjEFILUsWs9uzuLpO8eUmY3CsGHx+13kW18uAlt8yJw64fIiKSglcJBRGblbA1uZ47MiP2KG0IOREREaDASfP8naOshL3J9awtYeBOrk5sZ6h3qcr8PhERkVTMzKiIJxeAnDWolcOCWm9ObEdERCQWgxkVcXYBSEfEjpTy9sR2REREYjCYURGxk+aNToiSNGxZ7PwwckxsR0RE5AiDGRVxZXK9UQlNXJ4fRuwQ8ogaOs4FQ0REXsNgRkVcmVxvav+WkmbztSU1OR7jukdZfV6uaURERHJgMKMiUifXq5oVqaHTum2+FvNMDie2IyIiObEPQGUMAUPG/osmxcABmspARo6AghPbERGRnBjMyKysvALxs7MAwGKtIltSk+PxfJ/maD+vcuK8mUmxGN+zqeSAIm9eEmtaiIhI9fjntEpxyn8iIqJKvAISERGRqjGYISIiIlVjMENERESqxmCGiIiIVI1DWcguaytd+xJff31ERP6AmRkiIiJSNWZmfJh51qGsvMJtx82bl2ScH8fdmC0hIiIpGMwoDC/kRERE0rCbyY+tOViA8gq93M0gIiJyCYMZP7I466zJ/UVZ+YibtR3pmXkytYiIiMh1DGb8RHpmHlYevGSxXS8AK/ZdZEBDRESqxWDGD5RX6JGx/6LdfTL2X2SXExERqRKDGT+w9nAB9IL9ffRC5X5ERERqw2DGD1y6WebW/YiIiJSEwYwfiKpdw637ERERKQnnmfEDoxKi8VrmGbtdTQGayv3E4nw4RESkFMzMyKxq0a2n5n3RaQMwsWeM3X0m9oyBTsuPAxERqY9GEAQHpaHqV1paivDwcJSUlCAsLEzu5hilZ+YhY/9Fk4xJgKYysEhNjlf98xEREblC7PWbwYxM0jPzsGKf7eHSk3uJDzDKyiuM6yTlzUtCDZ3t3sPisnK0n7cTADAzKRbjezZlRoaIiBRJ7PWbVzEZyDnvS9XAZUz3aAYyRESkerySyYDzvhAREbkPgxkZcN4XIiIi92EwIwPO+0JEROQ+DGZkMCohGgEa+/tInfeFiIjIX8kezLz33nto27YtwsLCEBYWhoSEBGzfvt34+PXr1zFq1Cg0bNgQISEh6NChA7Zs2SJji13HeV+IiIjcR/arZePGjbFw4ULk5OTg+PHj6NevH5566imcPn0aADB69Gjk5+dj27ZtOHXqFJ555hkMGzYMJ0+elLnlrklNjsfkXjEWGZoAjbRh2URERP5O9mBm8ODBSE5ORosWLdCyZUu89tprCA0NRXZ2NgDg0KFDeOmll9C5c2c0bdoUr776KiIiIpCTkyNzy12XmhyPE7MGGO/PTIrF2fkDGcgQERFJIHswU9WDBw+wfv163LlzBwkJCQCAbt26YcOGDbh58yb0ej3Wr1+Pu3fvok+fPjaPc+/ePZSWlprclIrzvhAREblGEQtNnjp1CgkJCbh79y5CQ0OxdetWxMdXZic2btyI4cOHo06dOtBqtahRowa2bt2K5s2b2zxeeno65s6d663mExERkYwUkQaIjY1Fbm4ujhw5gueffx5jxoxBXl4eAGDWrFkoLi7G119/jePHj+OVV17BsGHDcOrUKZvHS01NRUlJifFWWFjorZdCREREXqaIzIxOpzNmWjp27Ihjx45h2bJl+Nvf/oZ//etf+P777/HII48AANq1a4f9+/dj+fLleP/9960eLygoCEFBQV5rPxEREclHEZkZc3q9Hvfu3UNZWeUMuAEBps0MDAyEXu/+dYuIiIhIfWTPzKSmpmLgwIGIjIzErVu38PHHH2Pv3r3IyspCXFwcmjdvjsmTJ2Px4sWoU6cOPvvsM+zcuRNffPGF3E0nIiIiBZA9mCkqKsLo0aNx7do1hIeHo23btsjKysKAAZVDljMzM5GSkoLBgwfj9u3baN68OdasWYPk5GSZW05ERERKIHsw8+GHH9p9vEWLFqqf8ZeIiIg8R5E1M0RERERiMZghIiIiVZO9m4lcV17xv5Fdaw4WYHzPpjZnEq6h06Jg4SBvNY2IiMjjmJlRufTMPHSYv9N4f1FWPuJmbUd6Zp6MrSIiIvIeZmZULD0zDyv2XbTYrhdg3M5FK4mIyNcxM6NS5RV6ZOy3DGSqyth/0aQLioiIyBcxmFGptYcLoBfs76MXKvcjIiLyZQxmVOrSzTK37kdERKRWDGZUKqp2DbfuR0REpFYsAJaZs0OlRyVE47XMM3a7mgI0lfsRERH5MmZmVEqnDcDEnjF295nYM8bmfDNERES+gpkZFTMMu87Yf9EkQxOgqQxkOCybiIj8gUYQBAdjYtSvtLQU4eHhKCkpQVhYmNzNcbvisnK0n1c5cd7MpFi7MwATERGphdjrN694PqBq4DKmezQDGSIi8iu86hEREZGqMZghIiIiVWMwQ0RERKrGYIaIiIhUjcEMERERqRqDGSIiIlI1BjNERESkagxmiIiISNUYzBAREZGqMZghIiIiVWMwQ0RERKrGYIaIiIhUjcEMERERqRqDGSIiIlI1BjNERESkagxmiIiISNW0cjeAXFdDp0XBwkFyN4OIiEgWzMwQERGRqjGYISIiIlVjMENERESqxmCGiIiIVI3BDBEREakagxkiIiJSNQYzREREpGoMZoiIiEjVGMwQERGRqjGYISIiIlWTPZh577330LZtW4SFhSEsLAwJCQnYvn27yT6HDx9Gv379EBISgrCwMPTq1Qu//fabTC0mIiIiJZE9mGncuDEWLlyInJwcHD9+HP369cNTTz2F06dPA6gMZJ544gk8/vjjOHr0KI4dO4YXX3wRAQGyN52IiIgUQCMIgiB3I8zVrl0bb7zxBsaPH4+uXbtiwIABmD9/vtPHKy0tRXh4OEpKShAWFubGlhIREZGniL1+Kyq98eDBA6xfvx537txBQkICioqKcOTIEdSvXx/dunVDgwYN0Lt3bxw4cEDuphIREZFCaOVuAACcOnUKCQkJuHv3LkJDQ7F161bEx8cjOzsbADBnzhwsXrwY7du3x7///W/0798f33//PVq0aGH1ePfu3cO9e/eM90tKSgBURnhERESkDobrtsNOJEEB7t27J5w/f144fvy4kJKSItStW1c4ffq0cPDgQQGAkJqaarJ/mzZthJSUFJvHS0tLEwDwxhtvvPHGG28+cCssLLQbRyiyZiYxMRHNmjVDSkoKmjZtirVr1+JPf/qT8fHhw4dDq9Vi3bp1Vn/ePDNTXFyMqKgoXL58GeHh4R5vv1KVlpaiSZMmKCws9OvaIZ6HSjwPlXge/ofnohLPQyUlnAdBEHDr1i00atTI7sAfRXQzmdPr9bh37x6io6PRqFEj5Ofnmzx+7tw5DBw40ObPBwUFISgoyGJ7eHi4X38wDQzD4P0dz0MlnodKPA//w3NRieehktznQUwSQvZgJjU1FQMHDkRkZCRu3bqFjz/+GHv37kVWVhY0Gg3++te/Ii0tDe3atUP79u2xZs0anD17Fps3b5a76URERKQAsgczRUVFGD16NK5du4bw8HC0bdsWWVlZGDBgAABg2rRpuHv3LqZPn46bN2+iXbt22LlzJ5o1ayZzy4mIiEgJZA9mPvzwQ4f7pKSkICUlxennCAoKQlpamtWuJ3/C81CJ56ESz0Mlnof/4bmoxPNQSU3nQZEFwERERERiKWrSPCIiIiKpGMwQERGRqjGYISIiIlVjMENERESq5jPBzPLlyxEdHY3g4GB06dIFR48etbv/pk2bEBcXh+DgYLRp0waZmZleaqlnSTkPGRkZ6NmzJ2rVqoVatWohMTHR4XlTC6mfB4P169dDo9FgyJAhnm2gl0g9D8XFxZgyZQoeeughBAUFoWXLlj7xuyH1PCxduhSxsbGoXr06mjRpgunTp+Pu3bteaq1n7Nu3D4MHD0ajRo2g0Wjw2WefOfyZvXv3okOHDggKCkLz5s2xevVqj7fT06Seh08//RQDBgxAvXr1EBYWhoSEBGRlZXmnsR7kzOfB4ODBg9BqtWjfvr3H2ieVTwQzGzZswCuvvIK0tDScOHEC7dq1Q1JSEoqKiqzuf+jQIYwYMQLjx4/HyZMnMWTIEAwZMgTff/+9l1vuXlLPw969ezFixAjs2bMHhw8fRpMmTfD444/jp59+8nLL3UvqeTAoKCjAjBkz0LNnTy+11LOknofy8nIMGDAABQUF2Lx5M/Lz85GRkYGHH37Yyy13L6nn4eOPP0ZKSgrS0tJw5swZfPjhh9iwYQP+/ve/e7nl7nXnzh20a9cOy5cvF7X/xYsXMWjQIPTt2xe5ubmYNm0aJkyYoPoLudTzsG/fPgwYMACZmZnIyclB3759MXjwYJw8edLDLfUsqefBoLi4GKNHj0b//v091DInubpIpBJ07txZmDJlivH+gwcPhEaNGgnp6elW9x82bJgwaNAgk21dunQRJk+e7NF2eprU82CuoqJCqFmzprBmzRpPNdErnDkPFRUVQrdu3YT/+7//E8aMGSM89dRTXmipZ0k9D++9957QtGlToby83FtN9Aqp52HKlClCv379TLa98sorQvfu3T3aTm8CIGzdutXuPn/729+ERx55xGTb8OHDhaSkJA+2zLvEnAdr4uPjhblz57q/QTKRch6GDx8uvPrqq0JaWprQrl07j7ZLCtVnZsrLy5GTk4PExETjtoCAACQmJuLw4cNWf+bw4cMm+wNAUlKSzf3VwJnzYK6srAz3799H7dq1PdVMj3P2PMybNw/169fH+PHjvdFMj3PmPGzbtg0JCQmYMmUKGjRogNatW+P111/HgwcPvNVst3PmPHTr1g05OTnGrqgLFy4gMzMTycnJXmmzUvji96Q76PV63Lp1S9Xfk85atWoVLly4gLS0NLmbYkH2GYBd9fPPP+PBgwdo0KCByfYGDRrg7NmzVn/m+vXrVve/fv26x9rpac6cB3MzZ85Eo0aNLL7A1MSZ83DgwAF8+OGHyM3N9UILvcOZ83DhwgXs3r0bI0eORGZmJn744Qe88MILuH//viK/vMRw5jz88Y9/xM8//4wePXpAEARUVFTgz3/+s+q7maSy9T1ZWlqK3377DdWrV5epZfJavHgxbt++jWHDhsndFK86f/48UlJSsH//fmi1ygsdVJ+ZIfdYuHAh1q9fj61btyI4OFju5njNrVu3MGrUKGRkZKBu3bpyN0dWer0e9evXxwcffICOHTti+PDh+Mc//oH3339f7qZ51d69e/H666/j3XffxYkTJ/Dpp5/iyy+/xPz58+VuGsns448/xty5c7Fx40bUr19f7uZ4zYMHD/DHP/4Rc+fORcuWLeVujlXKC68kqlu3LgIDA3Hjxg2T7Tdu3EDDhg2t/kzDhg0l7a8GzpwHg8WLF2PhwoX4+uuv0bZtW0820+Oknocff/wRBQUFGDx4sHGbXq8HAGi1WuTn56tyUVNnPg8PPfQQqlWrhsDAQOO2Vq1a4fr16ygvL4dOp/Nomz3BmfMwa9YsjBo1ChMmTAAAtGnTBnfu3MGkSZPwj3/8AwEB/vE3oK3vybCwML/Myqxfvx4TJkzApk2bVJ29dsatW7dw/PhxnDx5Ei+++CKAyu9JQRCg1WqxY8cO9OvXT9Y2qv63UqfToWPHjti1a5dxm16vx65du5CQkGD1ZxISEkz2B4CdO3fa3F8NnDkPAPDPf/4T8+fPx1dffYVOnTp5o6keJfU8xMXF4dSpU8jNzTXefve73xlHcDRp0sSbzXcbZz4P3bt3xw8//GAM5gDg3LlzeOihh1QZyADOnYeysjKLgMUQ4Al+tJSdL35POuuTTz7Bc889h08++QSDBg2SuzleFxYWZvE9+ec//xmxsbHIzc1Fly5d5G6ib4xmWr9+vRAUFCSsXr1ayMvLEyZNmiREREQI169fFwRBEEaNGiWkpKQY9z948KCg1WqFxYsXC2fOnBHS0tKEatWqCadOnZLrJbiF1POwcOFCQafTCZs3bxauXbtmvN26dUuul+AWUs+DOV8ZzST1PFy+fFmoWbOm8OKLLwr5+fnCF198IdSvX19YsGCBXC/BLaSeh7S0NKFmzZrCJ598Ily4cEHYsWOH0KxZM2HYsGFyvQS3uHXrlnDy5Enh5MmTAgDhrbfeEk6ePClcunRJEARBSElJEUaNGmXc/8KFC0KNGjWEv/71r8KZM2eE5cuXC4GBgcJXX30l10twC6nnYd26dYJWqxWWL19u8j1ZXFws10twC6nnwZzSRjP5RDAjCILwzjvvCJGRkYJOpxM6d+4sZGdnGx/r3bu3MGbMGJP9N27cKLRs2VLQ6XTCI488Inz55ZdebrFnSDkPUVFRAgCLW1pamvcb7mZSPw9V+UowIwjSz8OhQ4eELl26CEFBQULTpk2F1157TaioqPByq91Pynm4f/++MGfOHKFZs2ZCcHCw0KRJE+GFF14Qfv31V+833I327Nlj9ffd8NrHjBkj9O7d2+Jn2rdvL+h0OqFp06bCqlWrvN5ud5N6Hnr37m13f7Vy5vNQldKCGY0g+FHelIiIiHyO6mtmiIiIyL8xmCEiIiJVYzBDREREqsZghoiIiFSNwQwRERGpGoMZIiIiUjUGM0RERKRqDGaIiIhI1RjMEBERkaoxmCEiIiJVYzBDRKr0+uuvQ6PRWNyWLl0qd9OIyMu4NhMRqdKtW7dw584d4/3Zs2djx44dOHDgABo3bixjy4jI27RyN4CIyBk1a9ZEzZo1AQCzZs3Cjh07sHfvXgYyRH6I3UxEpGqzZ8/G2rVrsXfvXkRHR8vdHCKSAYMZIlKttLQ0/Pvf/2YgQ+TnGMwQkSqlpaVhzZo1DGSIiDUzRKQ+CxYswHvvvYdt27YhODgY169fBwDUqlULQUFBMreOiLyNo5mISFUEQUBERARKS0stHjt69Cgee+wxGVpFRHJiMENERESqxpoZIiIiUjUGM0RERKRqDGaIiIhI1RjMEBERkaoxmCEiIiJVYzBDREREqsZghoiIiFSNwQwRERGpGoMZIiIiUjUGM0RERKRqDGaIiIhI1RjMEBERkar9Pwusemg7q/oVAAAAAElFTkSuQmCC", + "text/plain": [ + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], "source": [ "# Plot data\n", "def plot_dist_mod():\n", " # YOUR CODE HERE\n", - " raise NotImplementedError()\n", + " plt.scatter(z_sample, mu_sample)\n", + " plt.errorbar(z_sample, mu_sample, yerr=dmu, fmt=\"o\")\n", + " #raise NotImplementedError()\n", " plt.xlabel('$z$')\n", " plt.ylabel('$\\mu$')\n", " plt.title('Distance modulus vs redshift')\n", @@ -403,6 +416,13 @@ "YOUR ANSWER HERE" ] }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] + }, { "cell_type": "markdown", "metadata": { @@ -4235,9 +4255,9 @@ ], "metadata": { "kernelspec": { - "display_name": "MLBDenv", + "display_name": "Python 3 (ipykernel)", "language": "python", - "name": "mlbdenv" + "name": "python3" }, "language_info": { "codemirror_mode": { @@ -4249,10 +4269,9 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.11.7" + "version": "3.11.11" } }, "nbformat": 4, "nbformat_minor": 4 } -