From 477fa46b33f69c87018290f3dfb9e9f53221aa06 Mon Sep 17 00:00:00 2001
From: Robert Lanzafame <R.C.Lanzafame@tudelft.nl>
Date: Tue, 17 Sep 2024 09:15:55 +0200
Subject: [PATCH] tutorial

---
 content/tutorials/Week_1_3/Analysis.ipynb     | 328 ++++++++++++++++++
 .../tutorials}/Week_1_3/data.csv              |   0
 2 files changed, 328 insertions(+)
 create mode 100644 content/tutorials/Week_1_3/Analysis.ipynb
 rename {src/students => content/tutorials}/Week_1_3/data.csv (100%)

diff --git a/content/tutorials/Week_1_3/Analysis.ipynb b/content/tutorials/Week_1_3/Analysis.ipynb
new file mode 100644
index 00000000..9e09ea32
--- /dev/null
+++ b/content/tutorials/Week_1_3/Analysis.ipynb
@@ -0,0 +1,328 @@
+{
+ "cells": [
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "# Week 1.3: Programming Tutorial\n",
+    "\n",
+    "<h1 style=\"position: absolute; display: flex; flex-grow: 0; flex-shrink: 0; flex-direction: row-reverse; top: 60px;right: 30px; margin: 0; border: 0\">\n",
+    "    <style>\n",
+    "        .markdown {width:100%; position: relative}\n",
+    "        article { position: relative }\n",
+    "    </style>\n",
+    "    <img src=\"https://gitlab.tudelft.nl/mude/public/-/raw/main/tu-logo/TU_P1_full-color.png\" style=\"width:100px\" />\n",
+    "    <img src=\"https://gitlab.tudelft.nl/mude/public/-/raw/main/mude-logo/MUDE_Logo-small.png\" style=\"width:100px\" />\n",
+    "</h1>\n",
+    "<h2 style=\"height: 10px\">\n",
+    "</h2>\n",
+    "\n",
+    "*[CEGM1000 MUDE](http://mude.citg.tudelft.nl/): Week 1.3. September 16, 2024.*\n",
+    "\n",
+    "_This notebook was prepared by Berend Bouvy and used in an in-class demonstration on Monday (the first of several programming tutorials)._"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 74,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "def add(a, b):\n",
+    "    result = a+b\n",
+    "    return result\n",
+    "\n",
+    "def gen_xhat(A, y):\n",
+    "    x_hat = np.linalg.inv(A.T @ A) @ A.T @ y\n",
+    "    return x_hat"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 6,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/plain": [
+       "3"
+      ]
+     },
+     "execution_count": 6,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
+   "source": [
+    "a = 1\n",
+    "b = 2\n",
+    "result = add(a, b)\n",
+    "\n",
+    "result"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 25,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "import numpy as np\n",
+    "import matplotlib.pyplot as plt"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 11,
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "(100, 2)\n"
+     ]
+    }
+   ],
+   "source": [
+    "# Replace 'file.csv' with the path to your CSV file\n",
+    "data = np.genfromtxt('data.csv', delimiter=',', skip_header=1)  # 'skip_header=1' skips the first row (header)\n",
+    "\n"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 37,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "t = data[:,0]\n",
+    "y = data[:,1]\n",
+    "n_rows = data.shape[0]\n",
+    "n_cols = data.shape[1]"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 61,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/plain": [
+       "Text(0, 0.5, 'y')"
+      ]
+     },
+     "execution_count": 61,
+     "metadata": {},
+     "output_type": "execute_result"
+    },
+    {
+     "data": {
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjsAAAHHCAYAAABZbpmkAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA/F0lEQVR4nO3dfXxU5Z338e8kQEKRDAQkkyBIFulKREVEIMK2VUNBMTdsqQVuUNS+dG8aWh76IHSLlIJS6G6xCIXK9lZaBLTtKgKVLgIvXWh4kIiVYnlqEBYzQXmYIJhIM3P/kXummTCTOTM5M+fMyef9euWPnJw5uTJq59vr+l2/yxUIBAICAABwqAyrBwAAAJBMhB0AAOBohB0AAOBohB0AAOBohB0AAOBohB0AAOBohB0AAOBohB0AAOBohB0AAOBohB0AAOBohB0AtvHHP/5RP/zhD3XhwgWrhwLAQQg7AGzjj3/8o+bNm0fYAWAqwg4AAHA0wg4AW/jhD3+o7373u5KkwsJCuVwuuVwunThxIuL9U6dO1TXXXKPLly9f9bMJEybI4/Govr5ekvT2229rxIgR6tq1q9q3b6/CwkI9+uijzY5n8uTJ6tq1q65cuXLVz7785S/rH//xH+P8CwFYhbADwBa+8pWvaMKECZKkJUuW6Ne//rV+/etf69prr414/7hx43Tp0iVt3rw57Prly5e1ceNGffWrX1VmZqbOnDmjL3/5yzpx4oRmzZqlZ599VhMnTtTu3bubHc+DDz6os2fP6g9/+EPYda/Xq+3bt2vSpEkt+GsBpFIbqwcAAJJ0yy23aMCAAVq3bp3GjBmjXr16NXv/sGHD1L17d7300kt64IEHQtc3b96sS5cuady4cZIa6oDOnz+v//qv/9LAgQND9y1YsKDZ599999267rrrtGbNGt1///2h6+vWrZPf7yfsAGmEmR0AacnlcumBBx7Q73//e33yySeh6y+99JK6d++uYcOGSZI6deokSdq0aVPEJaloMjIyNHHiRL322mu6ePFi6PqLL76oO++8U4WFheb8IQCSjrADIG2NGzdOn376qV577TVJ0ieffKLf//73euCBB+RyuSRJX/ziFzV27FjNmzdPXbt21ejRo/X888+rrq4u5vMfeughffrpp3rllVckSYcPH9b+/fv14IMPJu+PAmA6wg6AtDVkyBD16tVLL7/8siRp48aN+vTTT0NLWFLDDNBvf/tblZeXa+rUqTp9+rQeffRR3X777WEzQpEUFRXp9ttv15o1ayRJa9asUbt27fS1r30teX8UANMRdgDYRnA2Jh5f+9rXtGXLFtXU1Oill15Sr169NGTIkKvuGzJkiJ566im9/fbbevHFF/XnP/9Z69evj/n8hx56SNu3b1dVVZXWrl2rUaNGqXPnznGPE4B1CDsAbKNDhw6SFFdTwXHjxqmurk6rV6/Wli1brpp1OX/+vAKBQNi1/v37S5KhpawJEybI5XJp2rRp+utf/0phMpCGXIGm/ysAABbZt2+fBg0apPvuu0/jx49X27ZtVVpaGgpB0fTp00fV1dW6ePGi9u/frwEDBoR+9swzz+jnP/+5/vmf/1m9e/fWxYsXtWrVKlVVVenAgQOGCo1LS0u1adMmderUSV6vV1lZWS3+WwGkDjM7AGzjjjvu0Pz58/Xuu+/q4Ycf1oQJE/TRRx/FfN24ceN08eJF3XDDDWFBR2ooUB44cKDWr1+vb33rW1q8eLH69Omj7du3G95R9dBDD0lqWDIj6ADph5kdAIhhw4YNGjNmjN566y390z/9k9XDARAnwg4AxHD//ffr/fff17FjxxIqogZgLTooA0AU69ev15/+9Cdt3rxZP/vZzwg6QJpiZgcAonC5XLrmmms0btw4rVy5Um3a8P8PgXTEf7kAEAX/XxBwBnZjAQAARyPsAAAAR2MZS5Lf79eHH36ojh07UoAIAECaCAQCunjxogoKCpSREX3+hrAj6cMPP1SPHj2sHgYAAEjAqVOndN1110X9OWFHUseOHSU1vFk5OTkWjwYAABhRU1OjHj16hD7HoyHs6O8nLefk5BB2AABIM7FKUChQBgAAjkbYAQAAjkbYAQAAjkbYAQAAjkbYAQAAjkbYAQAAjkbYAQAAjkbYAQAAjkbYAQAAjkYHZQAAkBT1/oD2Vp7TmYu16tYxW4MKc5WZkfoDtwk7AADAdFsOVmnexkOq8tWGruW7szW3tEgj++WndCwsYwEAAFNtOVilKWsqwoKOJHl9tZqypkJbDlaldDyEHQAAYJp6f0DzNh5SIMLPgtfmbTyken+kO5KDsAMAAEyzt/LcVTM6jQUkVflqtbfyXMrGRNgBAACmOXMxetBJ5D4zEHYAAIBpunXMNvU+MxB2AACAaQYV5irfna1oG8xdatiVNagwN2VjIuwAAADTZGa4NLe0SJKuCjzB7+eWFqW03w5hBwAAmGpkv3ytmDRAHnf4UpXHna0VkwakvM8OTQUBAIDpRvbL1/AiDx2UAQCAc2VmuFTcu4vVw2AZCwAAOBthBwAAOBphBwAAOBphBwAAOBphBwAAOBphBwAAOBphBwAAOBphBwAAOBphBwAAOBphBwAAOBphBwAAOJqlYeett95SaWmpCgoK5HK59Oqrr4b9PBAI6Mknn1R+fr7at2+vkpISHT16NOyec+fOaeLEicrJyVGnTp309a9/XZ988kkK/woAAGBnloadS5cu6dZbb9Xy5csj/nzx4sVaunSpVq5cqT179qhDhw4aMWKEamtrQ/dMnDhRf/7zn7V161Zt2rRJb731lh5//PFU/QkAAMDmXIFAIGD1ICTJ5XLplVde0ZgxYyQ1zOoUFBTo29/+tr7zne9Iknw+n/Ly8vTCCy9o/Pjxev/991VUVKR9+/Zp4MCBkqQtW7bovvvu0//8z/+ooKDA0O+uqamR2+2Wz+dTTk5OUv4+AABgLqOf37at2amsrJTX61VJSUnomtvt1uDBg1VeXi5JKi8vV6dOnUJBR5JKSkqUkZGhPXv2RH12XV2dampqwr4AAIAz2TbseL1eSVJeXl7Y9by8vNDPvF6vunXrFvbzNm3aKDc3N3RPJAsXLpTb7Q599ejRw+TRAwDgbPX+gMqPn9WGA6dVfvys6v22WCiKqI3VA7DC7NmzNXPmzND3NTU1BB4AAAzacrBK8zYeUpXv7zW0+e5szRnVV507ZOnMxVp165itQYW5ysxwWTjSBrYNOx6PR5JUXV2t/Pz80PXq6mr1798/dM+ZM2fCXve3v/1N586dC70+kqysLGVlZZk/aAAAHG7LwSpNWVOhpvM4Vb5afWPtO2HX8t3ZmltapJH98mUl2y5jFRYWyuPxaNu2baFrNTU12rNnj4qLiyVJxcXFunDhgvbv3x+6Z/v27fL7/Ro8eHDKxwwAgJPV+wOat/HQVUEnGq+vVlPWVGjLwaqkjisWS2d2PvnkEx07diz0fWVlpQ4cOKDc3Fz17NlT06dP14IFC9SnTx8VFhZqzpw5KigoCO3Y6tu3r0aOHKnHHntMK1eu1JUrVzR16lSNHz/e8E4sAABgzN7Kc2FLV7EEJLkkzdt4SMOLPJYtaVkadt5++23dddddoe+DdTSTJ0/WCy+8oO9973u6dOmSHn/8cV24cEHDhg3Tli1blJ2dHXrNiy++qKlTp+qee+5RRkaGxo4dq6VLl6b8bwEAwOnOXDQedIICalji2lt5TsW9u5g/KANs02fHSvTZAQA4Ub0/oL2V50wrGC4/flYTVu1O6LU/G99fo/t3T/h3R2L089u2BcoAACBx0XZMGSkYjhaSBhXmKt+dLa+v1nDdTlC3jtmxb0oSwg4AAA4TbcdUsGB4xaQBUQNPrJA0t7RIU9ZUyCUZCjwuSR53Q2Cyim13YwEAgPg1t2MqeG3exkMRmwAGQ1LTIuTGu6pG9svXikkD5HHHnqkJLpjNLS2ytN8OMzsAADhIrB1TwYLhJVuPaOgNXUNLVEZC0qzfvaeO2W01vMij4UWesKWu85c+0/zN4TNCHpv02SHsAADgIEZ3TC3bcUzLdhwLLVG527eLua38wqdXNPE/9kSt/RnRz2NqQbRZWMYCAMBB4i0EDi5RvXEo+pmS0V7TtFlgZoZLxb27aHT/7iru3cUWQUci7AAA4CjBHVNGY0ZwieqVA6cN/45YtT92Q9gBAMBBMjNcmltaJElxBZ5zl64ot0O7uF4TbBZod4QdAAAcJp4dU42N6d9w1FI8i0+JdFVONQqUAQBwoJH98kM7pnYd+0jLdhyP+ZrhRR4NKsy9qs9Oc6xsFmgUYQcAAIcKFgwPKszV7ypOR+183LjxX2aGS8OLPNp9/KzK1lbowqdXIj7bDs0CjWIZCwAAh2uujidS47/MDJeG9umqH4+9WS6Dr7Ezwg4AAK1AtDoejzs76vERibzGjjj1XJx6DgBoPRI5Cd3s09PNwqnnAADgKsE6nmS/xk5YxgIAAI5G2AEAAI5G2AEAAI5G2AEAAI5G2AEAAI5G2AEAAI5G2AEAAI5G2AEAAI5G2AEAAI5G2AEAAI7GcREAAKRQoudMJet8Kruee2Umwg4AACmy5WCV5m08pCpfbehavjtbc0uLmj1BPNHXJWs86YZlLAAAUmDLwSpNWVMRFiwkyeur1ZQ1FdpysMrU1yVrPOmIsAMAQJLV+wOat/GQAhF+Frw2b+Mh1fvD70j0dckaT7oi7AAAkGR7K89dNYPSWEBSla9WeyvPteh19f6Ayo+f1YYDp1V+/GzUsJLoeNIVNTsAACTZmYvRg0Vz98XzunjqbxIdT7piZgcAgCTr1jE7ofuMvu7Ex5fjqr9JdDzpirADAECSDSrMVb47W9E2dLvUMAszqDA37td5crK0bu/JuOpvEh1PuiLsAACQZJkZLs0tLZKkqwJG8Pu5pUVX9bcx8roJg3rKWxNf/U2i40lXhB0AAFJgZL98rZg0QB53+NKQx52tFZMGRO1rE+t1vbp2MPT7m9bfJDqedESBMgAAKTKyX76GF3ni7ljc3OvKj5819Lsj1d8kOp50Q9gBACCFMjNcKu7dxbTXBetvvL7aiHU7LjXM1kSrv0l0POmEZSwAANJYa6u/SQRhBwCANNea6m8SwTIWAAAWMuvU8dZSf5MIwg4AABYx+9Tx1lB/kwiWsQAAsEBrOnXcaoQdAABSrLWdOm41wg4AAClm9NTxJVuPNHt6OYyhZgcAAJMYLTY2epr4sh3HtGzHsYh1PGYVNrcGhB0AAEwQT7FxvKeJB+t4gtvIzS5sdjqWsQAAaKLeH1D58bPacOC0oWWkeIuNY5063lTjOp7f/4nC5njZOuzU19drzpw5KiwsVPv27dW7d2/Nnz9fgcDf/6ULBAJ68sknlZ+fr/bt26ukpERHjx61cNQAgHS25WCVhi3argmrdmva+gOasGq3hi3aHjVEJFJs3FzX42iCdTw/2HCQwuY42TrsLFq0SCtWrNCyZcv0/vvva9GiRVq8eLGeffbZ0D2LFy/W0qVLtXLlSu3Zs0cdOnTQiBEjVFtrbD0UAICgRLaDGy023lt5Lux6tK7HsZy79Fncv6u1s3XNzh//+EeNHj1ao0aNkiT16tVL69at0969eyU1zOo888wz+sEPfqDRo0dLkn71q18pLy9Pr776qsaPH2/Z2AEA6SXWDI1LDbMmw4s8YYXARouNI93XuOvxrmMfadmO44kN3sDvas1sPbNz5513atu2bTpy5Igk6d1339XOnTt17733SpIqKyvl9XpVUlISeo3b7dbgwYNVXl4e9bl1dXWqqakJ+wIAtG6JztAYLTaOdl+w6/GM4f/YbB2PS1Juh7Yt+l2tla3DzqxZszR+/HjdeOONatu2rW677TZNnz5dEydOlCR5vV5JUl5eXtjr8vLyQj+LZOHChXK73aGvHj16JO+PAACkhURnaGIVG7vUsFNqUGFus881cnr5gtH9TPldrY2tw87LL7+sF198UWvXrlVFRYVWr16tf/u3f9Pq1atb9NzZs2fL5/OFvk6dOmXSiAEAdhdtp1WiMzRGQsrc0iJDPXBinV5+3y0Fpv2u1sTWNTvf/e53Q7M7knTzzTfrgw8+0MKFCzV58mR5PB5JUnV1tfLz/95XoLq6Wv3794/63KysLGVlZSV17AAA+2muP83wIo/y3dny+moj1u241BA6Is2aBENK02d7Euh9E+v0cjN/V2th67Bz+fJlZWSETz5lZmbK7/dLkgoLC+XxeLRt27ZQuKmpqdGePXs0ZcqUVA8XAGBjwZ1WTYNM44Z9c0uLNGVNhVxS2H1GZk1ihZR4xDq93Mzf1RrYOuyUlpbqqaeeUs+ePXXTTTfpnXfe0U9/+lM9+uijkiSXy6Xp06drwYIF6tOnjwoLCzVnzhwVFBRozJgx1g4eAGAbRnda7Xzi7hbNmsQKKWZK5e9Kd7YOO88++6zmzJmjb3zjGzpz5owKCgr0L//yL3ryySdD93zve9/TpUuX9Pjjj+vChQsaNmyYtmzZouxsKtEBAA3i2WnFrInzuAKN2xG3UjU1NXK73fL5fMrJybF6OAAAk204cFrT1h+Ied/PxvfX6P7dkz8gmMLo57etd2MBAGCGlvbCQXoj7AAAHM+sXjhIT4QdAIDjmdkLB+mHsAMAaBViNeyLtdMqWjNC2J+td2MBAGCmRHdaNdeMkCZ+9sduLLEbCwAQXbRmhMF4ZGRWCMnBbiwAAFooVjNCqaEZIUta9kbYAQAginiaEcK+CDsAAERx5mL0oJPIfbAGYQcAgChoRugMhB0AAKKgGaEzEHYAAIiCZoTOQNgBAJjGiY33WtqMENajqSAAwBRObryXaDPCen8g7tfAfDQVFE0FAaClaLx3NSeHP7ugqSAAICVS1XjPyBKZXZbRguGvaY8er69WU9ZUaMvBKkvG1VqxjAUAaJF4Gu8V9+6S0O8wMktil5mUWOHPpYbwN7zIw5JWijCzAwBokWQ33jMySxLtnipfrf7PmgrN3/jnlM300HXZfpjZAQC0SDIb7xmZJfnha3+W5Ip4T9Avd53QL3edSMlMD12X7YeZHQBAiySz8Z6RWRJvTZ28NcaCQypqZui6bD+EHQBAiySz8Z7Zsx+pOKmcrsv2Q9gBALRYshrvJWP2I9k1M3Rdth9qdgAApki08V5zgrMkXl9txJocl6S8nCxJLlXXRL4nmmTWzATDX9PdYR767FiCsAMAME1mhivh7eXRnje3tEhT1lTIJYWFmWCE+uH/ukmSIt7TnGTXzCQj/CExLGMBAGzNyBJZtHsiSWXNTDD8je7fXcW9uxB0LMJxEeK4CABIB0bOmQres/WQV/9314mos0Gt8fgKJzL6+c0yFgAgLRhZIgveU9y7iwYV5lIzA0mEHQCAQ1EzgyDCDgAgLRlZ1jJaMG3kWUhfhB0AQNox89BPuxwgiuRhNxYAIK0YORjUimfBvgg7AIC0EetgUMn4URBmPgv2RtgBAKQNIweDGj0Kwsxnwd4IOwCAtGH0iAcj95n5LNgbYQcAkDaMHvFg5D4znwV7I+wAANJG8GDQaJvC4zkKwsxnwd4IOwCAtBE8GFTSVSEl+P3c0iJDPXLMfBbsjbADAEgrRg4GteJZsC8OAhUHgQJAOjKz6zEdlNMTB4ECABzN6FEQqX4W7IdlLAAA4GiEHQAA4GiEHQAA4GiEHQAA4GiEHQAA4GjsxgIAIAFsV08ftp/ZOX36tCZNmqQuXbqoffv2uvnmm/X222+Hfh4IBPTkk08qPz9f7du3V0lJiY4ePWrhiAEALVXvD6j8+FltOHBa5cfPqt5vr5ZwWw5Wadii7ZqwaremrT+gCat2a9ii7dpysMrqoSECW8/snD9/XkOHDtVdd92l119/Xddee62OHj2qzp07h+5ZvHixli5dqtWrV6uwsFBz5szRiBEjdOjQIWVnc3gbAKSbLQerNG/jIVX5/n7aeL47W3NLi2zR0XjLwSpNWVOhpvHL66vVlDUVdF62IVt3UJ41a5Z27dql//7v/47480AgoIKCAn3729/Wd77zHUmSz+dTXl6eXnjhBY0fP97Q76GDMgDYQ7QgEVwcsjpI1PsDGrZoe1gQa8ylhqMmdj5xN0taKWD089vWy1ivvfaaBg4cqAceeEDdunXTbbfdplWrVoV+XllZKa/Xq5KSktA1t9utwYMHq7y8POpz6+rqVFNTE/YFALBWvT+geRsPXRV0JIWuzdt4yNIlrb2V56IGHalhnFW+Wu2tPJe6QSEmW4edv/71r1qxYoX69OmjP/zhD5oyZYq+9a1vafXq1ZIkr9crScrLywt7XV5eXuhnkSxcuFButzv01aNHj+T9EQAAQ9IhSJy5GH18idyH1LB12PH7/RowYICefvpp3XbbbXr88cf12GOPaeXKlS167uzZs+Xz+UJfp06dMmnEAIBEpUOQ6NbRWC2o0fuQGrYOO/n5+SoqKgq71rdvX508eVKS5PF4JEnV1dVh91RXV4d+FklWVpZycnLCvgAAqRFtp1U6BIlBhbnKd2crWjWOSw3F1IMKc1M5LMRg691YQ4cO1eHDh8OuHTlyRNdff70kqbCwUB6PR9u2bVP//v0lNRQr7dmzR1OmTEn1cAEAMTS302p4kUf57mx5fbUR63aCxb9WBonMDJfmlhZpypoKuaSwcQYD0NzSIoqTbcbWMzszZszQ7t279fTTT+vYsWNau3atnnvuOZWVlUmSXC6Xpk+frgULFui1117Te++9p4ceekgFBQUaM2aMtYMHAIQJ7rRqWpcT3LK99ZBXc0sbZvObRgU7BYmR/fK1YtIAedzhM0wed7blu8UQma23nkvSpk2bNHv2bB09elSFhYWaOXOmHnvssdDPA4GA5s6dq+eee04XLlzQsGHD9POf/1yf//znDf8Otp4DQHLFs2V76yGvrfvsBNFB2XpGP79tH3ZSgbADAMlVfvysJqzaHfO+qXfdoKE3dNXt13fW/g/OEyTQLKOf37au2QEApL96f0C7jn1s6N5lO45p2Y5joZmc0f27J3l0aA1sXbMDAEhvwTOklu04FtfrgnU8nDUFMxB2AABJEa0g2Qi7dEyGMxB2AACma+7oB6Ps0DEZzkDYAQCYLtbRD/Hg6AW0FGEHAGA6MwMKRy+gpdiNBQC4Skt7yBgNKLkd2un8pc9s2zEZzkDYAQCEae5IB6NN/YJnSMU6+mHOqCKVreXoBSQXy1gAgJBYRzoY3QoePENKav7oh/tu4egFJB8dlEUHZQDOEWn5SZKhJal4jnQwOttidJaIoxeQCDooA0ArEylYdPpcW0nShctXQteihY0XdlU2u4Oq8Vbw4t5dDI1pZL98DS/yxAwymRkuw88E4kXYAQAHCC4/NZ2qbxxygoJLUsFlokghqTnx7rQiyMBqhB0ASHPxNvALqGFJat7GQ/L7pbK1V4ek5rAVHOmGsAMAaS6RBn7BJakfbDhoOOiwFRzpit1YAJDmWtLA79ylzwzdx1ZwpDNmdgAgTUTbsZSKZSVPnH12ADsh7ABAGmhuC/fwIk+zDfwicUnq3KGtzl26uoC5qTmj+urhoYXM6CBtsYwFADYXq9Hf1kPeqA38Ignes2B0P+W7s6O+xqWGQEXQQboj7ACAjTW30yp4bd7GQxpe5InYibjT59qGeu0EBbsT33dLgaEux8GgU+8PqPz4WW04cFrlx8+q3t/qe9IiTbCMBQA2FmunVeNGf9Ea+AWfE6mp38h+Dcc1NF0ia1qjE2kZzZOTpQmDeqpX1w50PYatEXYAwMaM7rQK3hetgV9zTf1idTmO1rDQW1OnJW8cDX0f72GhQKrEvYw1efJkvfXWW8kYCwCgCaM7rVq6IysYkkb3767i3l3Clq6MNiyM97BQIFXiDjs+n08lJSXq06ePnn76aZ0+fToZ4wIASBpUmGuoiDieRn/x1N7E07CwcQ0R9Tywk7jDzquvvqrTp09rypQpeumll9SrVy/de++9+u1vf6srV2JvYQQAGJeZ4YqriDiWLQerNGzRdk1YtVvT1h/QhFW7NWzR9qizMfE2LGxcQwTYRUK7sa699lrNnDlT7777rvbs2aMbbrhBDz74oAoKCjRjxgwdPXo09kMAAIYEi4ib7rQK7qoyWiMTawt7pMCT6PJYS7o6A2ZrUYFyVVWVtm7dqq1btyozM1P33Xef3nvvPRUVFWnx4sWaMWOGWeMEgFYtVhFxLLG2sAcPBh1e5Al7ZnAZLZ6GhRKHhcJe4p7ZuXLlin73u9/p/vvv1/XXX6/f/OY3mj59uj788EOtXr1ab7zxhl5++WX96Ec/SsZ4AaDVilZEbEQ8W9ib/s54GxbGW0MEJFvcMzv5+fny+/2aMGGC9u7dq/79+191z1133aVOnTqZMDwAgBni3cLeWLRePE1xWCjsKu6ws2TJEj3wwAPKzo4+RdmpUydVVla2aGAAAPO0dAt702W0Ex9f1rq9J+Wtid6IELALVyAQaPX7A2tqauR2u+Xz+ZSTk2P1cAA4TLTTylM9hmGLtketvXGpIazsfOLuuOqArP670LoZ/fymgzIAJFFzp5WncgYkWHszZU2FXFJY4El0+Slat2bAbjgIFACSJJGt3slk1hZ2IN0wswMASZDoVu9ka+kWdiAdEXYAIAni2eqd6qUglp/Q2rCMBQBJ0JKt3gDMRdgBgCRI1WnlAGIj7ABAEiTjtHIAiSHsAEASmH1aOYDEEXYAIEnY6g3YA7uxACCJ2OoNWI+wAwBJZmSrN0cvAMlD2AEAi9nlSAnAqajZAQAL2e1ICcCJCDsAYJFYR0pIDUdK1Psj3QHAKMIOAFjE6JESL+yqJPAALUDYAQCLGD0qYv7m9zVs0XaWtIAEEXYAwCLxHBVBDQ+QuLQKOz/+8Y/lcrk0ffr00LXa2lqVlZWpS5cuuuaaazR27FhVV1dbN0gAjlXvD6j8+FltOHBa5cfPtnhpKdaREo1RwwMkLm22nu/bt0+/+MUvdMstt4RdnzFjhjZv3qzf/OY3crvdmjp1qr7yla9o165dFo0UgBNF2x4+Z1Rfde6QlVB/nOCRElPWVMglRSxUbixYw7O38lyobw/9eYDY0iLsfPLJJ5o4caJWrVqlBQsWhK77fD798pe/1Nq1a3X33XdLkp5//nn17dtXu3fv1pAhQ6waMgAHCW4PbxpGqny1+sbad8KuGemP0zSgLP/fAzR/86Fmi5UbC9b60J8HMCYtlrHKyso0atQolZSUhF3fv3+/rly5Enb9xhtvVM+ePVVeXh71eXV1daqpqQn7AoBImtseHkms2potB6s0bNF2TVi1W9PWH9CEVbs1f/MhzRnVV3NG9TX0O7p1zKY/DxAH24ed9evXq6KiQgsXLrzqZ16vV+3atVOnTp3Crufl5cnr9UZ95sKFC+V2u0NfPXr0MHvYAByg3h/QC7sqDc+4SM3X1jQXUMrWvqN8d3azNTwuNczc3H59Z/rzAHGwddg5deqUpk2bphdffFHZ2cZ3LcQye/Zs+Xy+0NepU6dMezYAZwjOwMzf/H7cr21cWxNkpIHg/M3va86oIkm6KvAEv59bWqT9H5w31J+n8e8HWjNbh539+/frzJkzGjBggNq0aaM2bdrozTff1NKlS9WmTRvl5eXps88+04ULF8JeV11dLY/HE/W5WVlZysnJCfsCgKBoMzDxatxHx2gDwc4d2mnFpAHyuMP/D57Hna0VkwZoZL98w/15jN4HOJ2tC5Tvuecevffee2HXHnnkEd1444164okn1KNHD7Vt21bbtm3T2LFjJUmHDx/WyZMnVVxcbMWQAaS5eGt0mtO4j048AWV0/+4aXuSJusvKaH+eePr4AE5m67DTsWNH9evXL+xahw4d1KVLl9D1r3/965o5c6Zyc3OVk5Ojb37zmyouLmYnFoCExJqBMapT+7byBwKq9weUmeGKO6BkZrhC28ubCvbn8fpqI4YylxpmggYV5iY4esBZbL2MZcSSJUt0//33a+zYsfrCF74gj8ej//zP/7R6WADSlFlLPxc+vaKJ/7EndMxDrAaCweJjIwEl2J8n+Lqmz5EaanvotwM0cAUCgVZfrl9TUyO32y2fz0f9DtDKlR8/qwmrdse8b86ovsp3t4/ZHycYN1ZMGiBJmrKmQlJ4A8HG98TTH4c+O2jtjH5+E3ZE2AHSmdkdhOv9AQ1btD3mEtHOJ+5WZoZL9f6Adh8/q7K1Fbrw6ZWIz2z8mq2HvKYGFDooozUz+vlt65odAGhOMmY2mjvCIdISUWaGSxkZrqhBRwrfCj6yX36zxceJjDdabQ+ABmlfswOgdUpmB+GR/fJjbv9uLN6t4MGAMrp/dxX37sJMDJBkzOwASDuxGvS51NBBeHiRJ+EgEc8MDFvBAXsj7ABIO0Yb9DU+Hbwxo3UuRpeI2AoO2BthB0DaaUkHYTvU+QBILWp2ANhWvT+g8uNnteHAaZUfPxs62DLRZSM71fkASB1mdgDYUnMzMMOLPHEvG9mtzgdA6jCzA8AS0WZtpNgzMFsPeaN2EJYawst9/RpCR/C58dT5tAQ7rQD7YWYHQMrFmrUxMgOz84m7tWLSgKuek+GS/AHpl7tO6Je7ToSeW/c3v6GxcVI44DyEHQApFZy1aRpmgrM200v6GJ6BabxstPWQV/931wn5mzy48XONYHs44DwsYwFImVh1M5L0/K4Thp7VuEHfoMJcvX7QG/G+4HPX7T0pT445B3ECSC+EHQApY6RuprljFxprPANj5LnemjpNGNRTUuSTwgOSxt/RQ5v+9OFVNUQA0hvLWABSxmg9TKf2beX79IrhnVZGn3v2Up2ml3xe6/aelLfm769xf66tJGnJG0dD1zg9HHAOZnYApIzRephHhhZKijwDI13doM/oc39V/oGWvHFEUkAzSvroZ+P7a0bJ5+W7fEUXLofPKJnReweAPRB2gFaqua3fyRI8ViFW3czUu2+Iq0FfrOc2VV1Tp2feOKq2GS6t33ey2RqieRsPsaQFpDlXIBBo9f8V19TUyO12y+fzKScnx+rhAEmXjCMT4vndU9ZUSIp8rELjMGP0DKvmnhuNS1LnDm117lLsGqF1jw0xdEYWgNQy+vnNzA7QyiTzyAQj4jlWIZ4GfdGeG01AMhR0JHrvAOmOAmWgFUnFkQlGJOtYhcbPff1glX5V/oEp46X3DpDeCDtAKxLPkQnJXrYJztrEw8iyVuPnGgk7uR3a6fylzwzv/AKQfgg7QCtidDnGjss28dYZBYuWYx0WOmdUkcrWVoR67TT+uXT1zi8A6YeaHaAVMbocY7dlm0TqjDIzXFEPC20cZO67xXgNEYD0xMwO0IoYne2w07JNS+qMgkXLTWeEPE1mhJJVQwTAHgg7QCsSnO2YsiZ9lm1aWmdkNMgkUkMEID0QdoBWxuhsh12YUWdEkAFaN8IO4ECxdi2l07JNutYZAbAPwg7gMEZ3LaXLbEc61hkBsBd2YwEOYnV35Ghacg6X0V1VdpyVAmAPzOwADmGX7shNmXEOV7rVGQGwF8IO4BB26o4cFJxpahrAgjNN8Rz6mU51RgDshbADOITduiPHM9O09ZDXUXVGAOyFmh3AIey2a8noTNOy7cdsWWcEwDkIO4BDBHctRVvUcalhtiRVu5aMziA9v6sy6uyP1DD7E09BMwA0RdgBHMJuu5aMziBd+PRK1J81rjMCgEQRdgAHCe5assOhlkZmmjq1b2voWXY8hR1A+qBAGXAYu+xaMnIO1yNDe2nJG0djPovuyABagpkdwIGCu5ZG9++u4t5dLNueHWumaerdfWxVZwTAmZjZAVIgVg8Zq56VCrFmmtLtFHYA6ccVCARa/TaHmpoaud1u+Xw+5eTkWD0cOIwZHYST8Sw7cerfBSC5jH5+E3ZE2EHyROsgHJyniFU03HgW58THl/XMG0cSflaymDXTlG4zVgCsZ/Tzm2UsIElaelZVpNmOSNL93KsguiMDSBYKlIEkieesqqainV6eyLOSxa4nrANAU4QdIEkSPauquRkhs35nS8WatZLofAzAPgg7QJIkelZVrBkhM35nS7Vk1goAUo2wAyRJomdVJTI7Y9dzr3Yd+4jZHQCWI+wASZLoWVXxzs7Y+dyrZTuOa9ii7dTvALCUrcPOwoULdccdd6hjx47q1q2bxowZo8OHD4fdU1tbq7KyMnXp0kXXXHONxo4dq+rqaotGDIRL5KyqWDNCTdnx3KvGKFgGYDVb99kZOXKkxo8frzvuuEN/+9vf9P3vf18HDx7UoUOH1KFDB0nSlClTtHnzZr3wwgtyu92aOnWqMjIytGvXLsO/hz47SLZ4e8gEdzpJV3cVDkiaUdJHvbp2sLQfTbQxRuJSQyjb+cTd9M4BYBpHNhX86KOP1K1bN7355pv6whe+IJ/Pp2uvvVZr167VV7/6VUnSX/7yF/Xt21fl5eUaMmSIoecSdmBH6dBV2GgvoKB1jw2hlw4A0ziyqaDP55Mk5eY2FGHu379fV65cUUlJSeieG2+8UT179owr7AB2ZJfTy5sTHOOSrUe0bMexmPenams8ADSWNmHH7/dr+vTpGjp0qPr16ydJ8nq9ateunTp16hR2b15enrxeb9Rn1dXVqa6uLvR9TU1NUsbcWtDm3xzR3kcrZ0Kajun26ztr/wfnrxrj0Bu6Ggo7qdoaDwCNpU3YKSsr08GDB7Vz584WP2vhwoWaN2+eCaNCOiy1pAM7vo+RxpThkhrvJA+OcXiRR/nubHl9tRHrd4I1O6naGg8Ajdl6N1bQ1KlTtWnTJu3YsUPXXXdd6LrH49Fnn32mCxcuhN1fXV0tj8cT9XmzZ8+Wz+cLfZ06dSpZQ3c0jgswhx3fx2hjatoyJzjGrYe8CW2zB4BUsHXYCQQCmjp1ql555RVt375dhYWFYT+//fbb1bZtW23bti107fDhwzp58qSKi4ujPjcrK0s5OTlhX4gPxwWYw47vYzzHVTQe4/AiT9zb7AEgFWy9jFVWVqa1a9dqw4YN6tixY6gOx+12q3379nK73fr617+umTNnKjc3Vzk5OfrmN7+p4uJiipOTLJ7jAth9E50d38d4j6toPMZ0KKoG0PrYOuysWLFCkvSlL30p7Przzz+vhx9+WJK0ZMkSZWRkaOzYsaqrq9OIESP085//PMUjbX0SPeQS4cx4H80uEE/0n1nwdVYXVQNAU7YOO0ZaAGVnZ2v58uVavnx5CkaEoEQPuUS4lr6PyShsTvSfGf+sAdiVrWt2YF+JHnKJcC15H5NV2BzvcRX8swZgd4QdJCTRQy6Tod4fUPnxs9pw4LTKj59Nq6LoRN/HZBY2NzempthpBSAdEHaQsEQOuTTbloNVGrZouyas2q1p6w9owqrdtjhlO54Alsj7GE9hcyKijalpnmGnFYB0kFZnYyULZ2O1TKQCWUlJ35ETXMZp+i9w8LdY9SGcaB1NPIXGGw6c1rT1B2KO5Wfj+2t0/+5x/w3RxhStgzIAWMGRZ2PBnpruvklFN+BYyzgu/b33Syo/jKMFsGAdTXMBLJ5dTKkqEI80JnZaAUg3LGPBVKnqBpzsZZxEpLJBIAXiAGAcYQemSeWHvR37/KQygNmpQBwA7I6wA9Ok8sPejn1+Uh3A7FAgDgDpgJodmCaVH/bBZRw7nbJtRQAz83gGszsxA4BdEHZgmlR+2AeXcaasqZBLCgs8Vi3jWBXAzDieIRVF5QBgFZaxYJpUF83abRknXetoUlVUDgBWoc+O6LNjpuAHpxR5tiUZIcRuyy/pNEtS7w9o2KLtUWutgrNRO5+423YhDQDoswNLBGdbmn7Ye5L4YW+3U7aN1tHYIaTFU1Rup/cYAOJB2IHpzCyaTVexAliqZ3+iBSs7buEHALMRdpAUdpttSVQyZl9a0mU50d8XLVjZcQs/AJiNsANEkYzZl1QfcxErWC3/37fZbgs/AJiN3ViwtXhODzdTsnYopbLxopGO1vM3v685o9JvBxkAxIOZHdhWtJmVOaP6qnOHrKTVAyVz9iWVNTJGg1XnDu1SXlQOAKlE2IEtRVt+qfLV6htr3wm7ZnZhb7w7lOKp60lljUw8wWp0/+6tvqgcgHMRdmA7zc2sRGJ2Ya/RkPD6wSrtrTyndXtPyltjrK4nlV2W4w1WTikqB4CmqNmB7cSaWWnK7BPVjYaEX5V/oCVvHAkLOlLzdT2p7LKc6o7WAGBXhB3YTiL1KmYW9sYKCUbGIkUPX6k65iJdj68AALOxjJUkyeyOa4fOu8nUknoVMwp7mztk1KhYnYdT1XjRio7WAGA3hJ0kSGZ33HQ6dylRsepammNW87toISFezYWvVNXI0NEaQGvHQaAy9yDQaLuIEj0Is/EszomPL+uZN46Y9mw7i3agaDTJOrAy+P6/frBKvyr/IO7Xr3tsCEW/AJAkHARqAbP7s0SaxYkkGZ13g6xaMotnZiWZ9SeNZ1/iCTt0HgYA+yDsmMjME6SjzRCZ8WyjrF4yi7T8cv7SZ5q/2Zz6k3iCXDxLaxT/AoC9EHZMZFZ33Hj7zBh9djwf7qk+rDKaSHUtI/q1vP4k3iAXT9Eyxb8AYC+EHROZ1R033j4zRp4dz4d7qg+rjFdLC3sTDXJRdzblZGnCoJ7q1bUDxb8AYEOEHROZ1R03ke3TzT073g93M5fj7CI4q+X1far5m99POMixswkA0g9hx0TNLXXEU8cR7/bp5p6dyCxNPMclSLL9h73RQm/p70FuydYjGnpD14h/G8cqAEB6oYOyyczojhtvB9/mnh3PLE1QPMclTFi1W8MWbY94NIIdBGe14l0WXLbjmO3/NgCAMczsJEFLlzpizRAFJM0o6WOoRiSRoul4m/qlumjZqJYUegfZ9W8DABjHzE6SBJc6RvfvruLeXeJe5mluhmjlpAGaVvJ5Q89OpGi6uTOVIjH7IE6ztKTQO8iufxsAwDhmdmzMjGLYRIum4z0uwY5Fy2ackyXZ828DABhH2LG5lhbDtqRounHYMnpcglkBwwxmnZMVZKe/DQBgHMtYrUBLiqaDYeteg/UqzQWMen9A5cfPasOB0yo/fjbpy0JGCr1zO7RV2V29DT3P7PAEAEgNZnYcIlZ35JYuibW0h5AVR08YmdV6+p9v1vAij/6z4nSL+yMBAOyJU89l7qnnVkhVkIh2EnmsU9fNPgk+Xkben0T/NgCAdYx+fhN2ZK+wE+8p46kOEvEGq3p/QMMWbY9a5BycNdn5xN1JbUxo5H21+uBTAEB8CDtxsEvYSZcgEU8gKz9+VhNW7Y75zHWPDbHFTqd4wyYAwDpGP7+p2bGJRA6ntOoMq3h2iJl1EnyqcBQEADgPu7FsINb5VVLkpnbpECTMOgkeAIBEEXZsIJHzq6T0CBKxtn+71LBUx04nAECyEHZsINEZmnQIEs0dPRGpqWGqe/EAAJyPmh0bSHSGpiXdkVMp2tETngjbv9kNBQAwG7uxZP1urOCuqlhN7aLtqkqXkNDcTiere/EAANJPq9t6vnz5cv3kJz+R1+vVrbfeqmeffVaDBg0y9Fqrw47U8qZ26bxl2i69eAAA6cXo57cjanZeeuklzZw5U3PnzlVFRYVuvfVWjRgxQmfOnLF6aIa15Pwq6e9bpkf3767i3l3SKhQkWqANAIARjpjZGTx4sO644w4tW7ZMkuT3+9WjRw9985vf1KxZs2K+3g4zO0HpPEMTr+DfavRE9Z+N76/R/bunYGQAgHTQapoKfvbZZ9q/f79mz54dupaRkaGSkhKVl5dHfE1dXZ3q6upC39fU1CR9nEa1lqZ2keqMYqEXDwAgEWm/jPXxxx+rvr5eeXl5Ydfz8vLk9XojvmbhwoVyu92hrx49eqRiqPj/gvVJRoOOHbbQAwDSV9qHnUTMnj1bPp8v9HXq1Cmrh9RqNNctOhI7baEHAKSntF/G6tq1qzIzM1VdXR12vbq6Wh6PJ+JrsrKylJWVlYrhoYlYxchNNe3FAwBAvNJ+Zqddu3a6/fbbtW3bttA1v9+vbdu2qbi42MKRIRKj3aIfKr5e6x4bop1P3E3QAQC0SNrP7EjSzJkzNXnyZA0cOFCDBg3SM888o0uXLumRRx6xemhowmiR8b398ltFoTYAIPkcEXbGjRunjz76SE8++aS8Xq/69++vLVu2XFW0DOsFz/OK1S2aYmQAgFkc0WenpezUZ6c1aGm3aAAApFbWQTldcKJ3g5Z2izaC9xoAEOSIZax0YMfDOq3s1jyyX76GF3mS8vvt+F4DAKzDMpaSv4xlxYnesYKMUwMBp6cDQOvR6k49b4lkhh0rTvSOFWScGgg4PR0AWhdqdmwi1Sd6RzuKweur1ZQ1Ffr9nz6M2sE4eG3exkNpWePC6ekAgEgIO0lmtIme0fua09xRDMFrP9hw0LGBIJXvNQAgfRB2ksxoEz0zTvQ2MrNx7tIVQ89Kx0CQyvcaAJA+CDtJFmyiF61CxMwTvc0MKOkYCFL5XgMA0gdhJ8kyM1yaW1okSVd9CJt9orfRgJLboZ0jA0Eq32sAQPog7KRAKproScZnNhaM7hf6vunPpfQOBKl6rwEA6YOt50rdcRGpaOJn9CgGp/bZCbKyYSIAIDXosxMHq87GStYHstEgQyAAAKQzwk4crAg7yZ5ZIcgAAJyOsBOHVIcdp3YwBgAgleigbFNGGv9Z2cGY08IBAE7DqecpFs+RBsW9u6RuYHLu4aAAgNaNmZ0Us+uRBrHO1NpysCql4wEAwCyEnRSz45EGdl9aAwCgJQg7KWbHIw04LRwA4GSEnRSz45EGdl1aAwDADIQdC9jtSAM7Lq0BAGAWdmNZZGS/fA0v8tii8V9wac3rq41Yt+NSQxBLx8NBAQAg7FgoM8OV8u3l0cYxt7RIU9ZUyKXIZ2ql8+GgAIDWjWUsSLLf0hoAAGZhZgchdlpaAwDALIQdhLHL0hoAAGZhGQsAADgaYQcAADgaYQcAADgaYQcAADgaYQcAADgaYQcAADgaYQcAADgaYQcAADgaYQcAADgaHZQlBQINR1/W1NRYPBIAAGBU8HM7+DkeDWFH0sWLFyVJPXr0sHgkAAAgXhcvXpTb7Y76c1cgVhxqBfx+vz788EN17NhRLpd5h17W1NSoR48eOnXqlHJyckx7LsLxPqcO73Vq8D6nBu9zaiTzfQ4EArp48aIKCgqUkRG9MoeZHUkZGRm67rrrkvb8nJwc/kNKAd7n1OG9Tg3e59TgfU6NZL3Pzc3oBFGgDAAAHI2wAwAAHI2wk0RZWVmaO3eusrKyrB6Ko/E+pw7vdWrwPqcG73Nq2OF9pkAZAAA4GjM7AADA0Qg7AADA0Qg7AADA0Qg7AADA0Qg7SbR8+XL16tVL2dnZGjx4sPbu3Wv1kBxl4cKFuuOOO9SxY0d169ZNY8aM0eHDh60eluP9+Mc/lsvl0vTp060eiuOcPn1akyZNUpcuXdS+fXvdfPPNevvtt60elqPU19drzpw5KiwsVPv27dW7d2/Nnz8/5tlKiO2tt95SaWmpCgoK5HK59Oqrr4b9PBAI6Mknn1R+fr7at2+vkpISHT16NCVjI+wkyUsvvaSZM2dq7ty5qqio0K233qoRI0bozJkzVg/NMd58802VlZVp9+7d2rp1q65cuaIvf/nLunTpktVDc6x9+/bpF7/4hW655Rarh+I458+f19ChQ9W2bVu9/vrrOnTokP793/9dnTt3tnpojrJo0SKtWLFCy5Yt0/vvv69FixZp8eLFevbZZ60eWtq7dOmSbr31Vi1fvjzizxcvXqylS5dq5cqV2rNnjzp06KARI0aotrY2+YMLICkGDRoUKCsrC31fX18fKCgoCCxcuNDCUTnbmTNnApICb775ptVDcaSLFy8G+vTpE9i6dWvgi1/8YmDatGlWD8lRnnjiicCwYcOsHobjjRo1KvDoo4+GXfvKV74SmDhxokUjciZJgVdeeSX0vd/vD3g8nsBPfvKT0LULFy4EsrKyAuvWrUv6eJjZSYLPPvtM+/fvV0lJSehaRkaGSkpKVF5ebuHInM3n80mScnNzLR6JM5WVlWnUqFFh/17DPK+99poGDhyoBx54QN26ddNtt92mVatWWT0sx7nzzju1bds2HTlyRJL07rvvaufOnbr33nstHpmzVVZWyuv1hv3vh9vt1uDBg1PyuchBoEnw8ccfq76+Xnl5eWHX8/Ly9Je//MWiUTmb3+/X9OnTNXToUPXr18/q4TjO+vXrVVFRoX379lk9FMf661//qhUrVmjmzJn6/ve/r3379ulb3/qW2rVrp8mTJ1s9PMeYNWuWampqdOONNyozM1P19fV66qmnNHHiRKuH5mher1eSIn4uBn+WTIQdOEJZWZkOHjyonTt3Wj0Uxzl16pSmTZumrVu3Kjs72+rhOJbf79fAgQP19NNPS5Juu+02HTx4UCtXriTsmOjll1/Wiy++qLVr1+qmm27SgQMHNH36dBUUFPA+OxjLWEnQtWtXZWZmqrq6Oux6dXW1PB6PRaNyrqlTp2rTpk3asWOHrrvuOquH4zj79+/XmTNnNGDAALVp00Zt2rTRm2++qaVLl6pNmzaqr6+3eoiOkJ+fr6KiorBrffv21cmTJy0akTN997vf1axZszR+/HjdfPPNevDBBzVjxgwtXLjQ6qE5WvCzz6rPRcJOErRr10633367tm3bFrrm9/u1bds2FRcXWzgyZwkEApo6dapeeeUVbd++XYWFhVYPyZHuuecevffeezpw4EDoa+DAgZo4caIOHDigzMxMq4foCEOHDr2qdcKRI0d0/fXXWzQiZ7p8+bIyMsI/+jIzM+X3+y0aUetQWFgoj8cT9rlYU1OjPXv2pORzkWWsJJk5c6YmT56sgQMHatCgQXrmmWd06dIlPfLII1YPzTHKysq0du1abdiwQR07dgyt+7rdbrVv397i0TlHx44dr6qD6tChg7p06UJ9lIlmzJihO++8U08//bS+9rWvae/evXruuef03HPPWT00RyktLdVTTz2lnj176qabbtI777yjn/70p3r00UetHlra++STT3Ts2LHQ95WVlTpw4IByc3PVs2dPTZ8+XQsWLFCfPn1UWFioOXPmqKCgQGPGjEn+4JK+36sVe/bZZwM9e/YMtGvXLjBo0KDA7t27rR6So0iK+PX8889bPTTHY+t5cmzcuDHQr1+/QFZWVuDGG28MPPfcc1YPyXFqamoC06ZNC/Ts2TOQnZ0d+Id/+IfAv/7rvwbq6uqsHlra27FjR8T/TZ48eXIgEGjYfj5nzpxAXl5eICsrK3DPPfcEDh8+nJKxuQIB2kYCAADnomYHAAA4GmEHAAA4GmEHAAA4GmEHAAA4GmEHAAA4GmEHAAA4GmEHAAA4GmEHAAA4GmEHgGN96Utf0vTp060eBgCLEXYAAICjcVwEAEd6+OGHtXr16rBrlZWV6tWrlzUDAmAZwg4AR/L5fLr33nvVr18//ehHP5IkXXvttcrMzLR4ZABSrY3VAwCAZHC73WrXrp0+97nPyePxWD0cABaiZgcAADgaYQcAADgaYQeAY7Vr10719fVWDwOAxQg7AByrV69e2rNnj06cOKGPP/5Yfr/f6iEBsABhB4Bjfec731FmZqaKiop07bXX6uTJk1YPCYAF2HoOAAAcjZkdAADgaIQdAADgaIQdAADgaIQdAADgaIQdAADgaIQdAADgaIQdAADgaIQdAADgaIQdAADgaIQdAADgaIQdAADgaIQdAADgaP8PCye4QBeHXxwAAAAASUVORK5CYII=",
+      "text/plain": [
+       "<Figure size 640x480 with 1 Axes>"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    }
+   ],
+   "source": [
+    "plt.plot(t, y,'o')\n",
+    "plt.title('t vs y')\n",
+    "plt.xlabel('t')\n",
+    "plt.ylabel('y')"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 39,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "one_vector = np.ones(n_rows)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 42,
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "[ 1.          1.1010101   1.2020202   1.3030303   1.4040404   1.50505051\n",
+      "  1.60606061  1.70707071  1.80808081  1.90909091  2.01010101  2.11111111\n",
+      "  2.21212121  2.31313131  2.41414141  2.51515152  2.61616162  2.71717172\n",
+      "  2.81818182  2.91919192  3.02020202  3.12121212  3.22222222  3.32323232\n",
+      "  3.42424242  3.52525253  3.62626263  3.72727273  3.82828283  3.92929293\n",
+      "  4.03030303  4.13131313  4.23232323  4.33333333  4.43434343  4.53535354\n",
+      "  4.63636364  4.73737374  4.83838384  4.93939394  5.04040404  5.14141414\n",
+      "  5.24242424  5.34343434  5.44444444  5.54545455  5.64646465  5.74747475\n",
+      "  5.84848485  5.94949495  6.05050505  6.15151515  6.25252525  6.35353535\n",
+      "  6.45454545  6.55555556  6.65656566  6.75757576  6.85858586  6.95959596\n",
+      "  7.06060606  7.16161616  7.26262626  7.36363636  7.46464646  7.56565657\n",
+      "  7.66666667  7.76767677  7.86868687  7.96969697  8.07070707  8.17171717\n",
+      "  8.27272727  8.37373737  8.47474747  8.57575758  8.67676768  8.77777778\n",
+      "  8.87878788  8.97979798  9.08080808  9.18181818  9.28282828  9.38383838\n",
+      "  9.48484848  9.58585859  9.68686869  9.78787879  9.88888889  9.98989899\n",
+      " 10.09090909 10.19191919 10.29292929 10.39393939 10.49494949 10.5959596\n",
+      " 10.6969697  10.7979798  10.8989899  11.        ]\n"
+     ]
+    }
+   ],
+   "source": [
+    "print(one_vector+ t)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 79,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "A = np.column_stack((one_vector, t))"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 75,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "x_hat = gen_xhat(A, y)\n",
+    "y_hat = A @ x_hat"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 76,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/plain": [
+       "Text(0, 0.5, 'y')"
+      ]
+     },
+     "execution_count": 76,
+     "metadata": {},
+     "output_type": "execute_result"
+    },
+    {
+     "data": {
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjsAAAHHCAYAAABZbpmkAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAABZlklEQVR4nO3dd3xUdb7/8ddMKi1DZ0KTiKhEUJoUwYZBUERRLLCoqKx7l4WlWVERUQTBXQuKsPaC2BelCIrAgmAoguyCIE0QBBKkJbS0mXP/OE5MyEwyk0zP+/l48Pj9zuTMmW9y7+6877d8PhbDMAxEREREopQ11AMQERERCSSFHREREYlqCjsiIiIS1RR2REREJKop7IiIiEhUU9gRERGRqKawIyIiIlFNYUdERESimsKOiIiIRDWFHREREYlqCjsiEja+++47nnjiCY4dOxbqoYhIFFHYEZGw8d133zF+/HiFHRHxK4UdERERiWoKOyISFp544gkeeOABAFJSUrBYLFgsFnbv3u32/mHDhlG9enVOnTpV4mcDBgzAbrfjcDgA+P777+nZsyd169alSpUqpKSkcM8995Q6nkGDBlG3bl3y8/NL/Ozqq6/mvPPO8/E3FJFQUdgRkbBw0003MWDAAACef/553nvvPd577z3q1avn9v7bbruNkydPMn/+/GKvnzp1irlz53LzzTcTExPDwYMHufrqq9m9ezcPP/wwL730EgMHDmTVqlWljueOO+7g8OHDfPXVV8Vez8jIYMmSJdx+++0V+G1FJJhiQz0AERGACy+8kHbt2vHBBx/Qt29fmjVrVur93bp1o1GjRnz00Ufccsstha/Pnz+fkydPcttttwHmPqCjR4/y9ddf06FDh8L7JkyYUOrzu3fvTuPGjZk5cybXXXdd4esffPABTqdTYUckgmhmR0QiksVi4ZZbbuHLL7/kxIkTha9/9NFHNGrUiG7dugFQs2ZNAObNm+d2ScoTq9XKwIEDmTNnDsePHy98/f333+eSSy4hJSXFP7+IiAScwo6IRKzbbruN06dPM2fOHABOnDjBl19+yS233ILFYgHg8ssvp1+/fowfP566detyww038NZbb5Gbm1vm8++8805Onz7N7NmzAdi6dSvr1q3jjjvuCNwvJSJ+p7AjIhGrc+fONGvWjI8//hiAuXPncvr06cIlLDBngD799FPS09MZNmwY+/bt45577qF9+/bFZoTcSU1NpX379sycOROAmTNnEh8fz6233hq4X0pE/E5hR0TChms2xhe33norCxcuJDs7m48++ohmzZrRuXPnEvd17tyZp59+mu+//57333+fH3/8kQ8//LDM5995550sWbKEAwcOMGvWLHr37k2tWrV8HqeIhI7CjoiEjWrVqgH4VFTwtttuIzc3l3feeYeFCxeWmHU5evQohmEUe61NmzYAXi1lDRgwAIvFwogRI/j555+1MVkkAlmMM/9bQEQkRNauXUvHjh259tpr6d+/P3FxcfTp06cwBHnSokULMjMzOX78OOvWraNdu3aFP3vhhRd45ZVXuPHGG2nevDnHjx/ntdde48CBA2zYsMGrjcZ9+vRh3rx51KxZk4yMDBISEir8u4pI8GhmR0TCxsUXX8xTTz3Ff//7X+666y4GDBjAb7/9Vub7brvtNo4fP84555xTLOiAuUG5Q4cOfPjhhwwfPpwpU6bQokULlixZ4vWJqjvvvBMwl8wUdEQij2Z2RETK8MUXX9C3b1+WL1/OpZdeGurhiIiPFHZERMpw3XXXsWXLFnbs2FGuTdQiElqqoCwi4sGHH37I//73P+bPn8+LL76ooCMSoTSzIyLigcVioXr16tx2223MmDGD2Fj934cikUj/yRUR8UD/t6BIdNBpLBEREYlqCjsiIiIS1bSMBTidTvbv30+NGjW0AVFERCRCGIbB8ePHadiwIVar5/kbhR1g//79NGnSJNTDEBERkXLYu3cvjRs39vhzhR2gRo0agPnHSkpKCvFoRERExBvZ2dk0adKk8HvcE4Ud/ui0nJSUpLAjIiISYcragqINyiIiIhLVFHZEREQkqinsiIiISFRT2BEREZGoprAjIiIiUU1hR0RERKKawo6IiIhENYUdERERiWoKOyIiIhLVVEFZREREAsLhNFiz6wgHj+dQv0YiHVNqE2MNfsNthR0RERHxu4WbDjB+7mYOZOUUvpZsS2Rcn1R6tUoO6li0jCUiIiJ+tXDTAYbMXF8s6ABkZOUwZOZ6Fm46ENTxKOyIiIiI3zicBuPnbsZw8zPXa+PnbsbhdHdHYCjsiIiIiN+s2XWkxIxOUQZwICuHNbuOBG1MCjsiIiLiNwePew465bnPHxR2RERExG/q10j0633+oLAjIiIiftMxpTbJtkQ8HTC3YJ7K6phSO2hjUtgRERERv4mxWhjXJxWgROBxXY/rkxrUejsKOyIiIuJXvVolM/32dthtxZeq7LZEpt/eLuh1dlRUUERERPyuV6tkeqTaVUFZREREoleM1UKX5nVCPQwtY4mIiEh0U9gRERGRqKawIyIiIlFNYUdERESimsKOiIiIRDWFHREREYlqCjsiIiIS1RR2REREJKop7IiIiEhUU9gRERGRqBbSsLN8+XL69OlDw4YNsVgsfP7558V+bhgGjz/+OMnJyVSpUoW0tDS2b99e7J4jR44wcOBAkpKSqFmzJoMHD+bEiRNB/C1EREQknIU07Jw8eZKLLrqIadOmuf35lClTmDp1KjNmzGD16tVUq1aNnj17kpOTU3jPwIED+fHHH1m0aBHz5s1j+fLl/OUvfwnWryAiIiJhzmIYhhHqQQBYLBZmz55N3759AXNWp2HDhtx3333cf//9AGRlZdGgQQPefvtt+vfvz5YtW0hNTWXt2rV06NABgIULF3Lttdfy66+/0rBhQ68+Ozs7G5vNRlZWFklJSQH5/URERMS/vP3+Dts9O7t27SIjI4O0tLTC12w2G506dSI9PR2A9PR0atasWRh0ANLS0rBaraxevTroYxYREZHwExvqAXiSkZEBQIMGDYq93qBBg8KfZWRkUL9+/WI/j42NpXbt2oX3uJObm0tubm7hdXZ2tr+GLSIiUik4nAZrdh3h4PEc6tdIpGNKbWKsllAPy62wDTuBNGnSJMaPHx/qYYiIiESkhZsOMH7uZg5k/bGHNtmWyNjeLalVLSHsAlDYhh273Q5AZmYmycnJha9nZmbSpk2bwnsOHjxY7H0FBQUcOXKk8P3ujBkzhtGjRxdeZ2dn06RJEz+OXkREJDot3HSAITPXc+aG3wNZOfxt1g/FXku2JTKuTyq9WiUTSmG7ZyclJQW73c7ixYsLX8vOzmb16tV06dIFgC5dunDs2DHWrVtXeM+SJUtwOp106tTJ47MTEhJISkoq9k9ERERK53AajJ+7uUTQ8SQjK4chM9ezcNOBgI6rLCGd2Tlx4gQ7duwovN61axcbNmygdu3aNG3alJEjRzJhwgRatGhBSkoKY8eOpWHDhoUntlq2bEmvXr249957mTFjBvn5+QwbNoz+/ft7fRJLREREvLNm15FiS1dlMQALMH7uZnqk2kO2pBXSsPP9999z5ZVXFl67lpYGDRrE22+/zYMPPsjJkyf5y1/+wrFjx+jWrRsLFy4kMTGx8D3vv/8+w4YN46qrrsJqtdKvXz+mTp0a9N9FREQk2h087n3QcTEwl7jW7DpCl+Z1/D8oL4RNnZ1QUp0dERGJRv4+MZW+8zADXltVrve+2L8NN7RpVO7Pdsfb7++w3aAsIiIi5efpxJQ3G4Y9haSOKbVJtiWSkZXj9b4dl/o1Esu+KUAUdkRERKKMpxNTrg3D029v5zHwlBWSxvVJZcjM9VjAq8BjAew2MzCFStiexhIRERHflXZiyvXa+LmbcThL3uEKSWduQi56qqpXq2Sm394Ou63smRrXgtm4PqkhrbejmR0REZEoUtaJKdeG4ecXbaPrOXULl6i8CUkPf7aRGolx9Ei10yPVXmyp6+jJPJ6aX3xGyB4mdXYUdkRERKKItyemXl66g5eX7ihcorJViS/zWPmx0/kMfH21x70/PVvZw7KFhJaxREREooivG4FdS1TfbPbcU9LTe84sFhhjtdCleR1uaNOILs3rhEXQAYUdERGRqOI6MeVtzHAtUc3esM/rzyhr70+4UdgRERGJIjFWC+P6pAL4FHiOnMyndrV4n97jKhYY7hR2REREoowvJ6aK6tvGbLXky+JTeaoqB5s2KIuIiEShXq2SC09MrdzxGy8v3Vnme3qk2umYUrtEnZ3ShLJYoLcUdkRERKKUa8Nwx5TafLZ+n8fKx0UL/8VYLfRItbNq52GGzlrPsdP5bp8dDsUCvaVlLBERkShX2j4ed4X/YqwWuraoyzP9WmPx8j3hTGFHRESkEvC0j8duS/TYPqI87wlH6nqOup6LiEjlUZ5O6P7unu4v6nouIiIiJbj28QT6PeFEy1giIiIS1RR2REREJKop7IiIiEhUU9gRERGRqKawIyIiIoEV4oPfCjsiIiISGI58+PY5ePd6cDpDNgwdPRcRERH/+3UdzB0OmZvM65/mQer1IRmKwo6IiIj4T+5xWPI0rJ4BGFClNvScCC37hGxICjsiIiLiH9u+gnmjIftX8/rC28ygU61uSIelsCMiIhJE5W29EKiWDX557vFMWPgQ/DjbvK7ZFK57Hs5Jq/D4/EFhR0REJEgWbjrA+LmbOZCVU/hasi2RcX1SS22qWd73BWo8hQwDfngPvn4McrLAYoXOf4MrH4H4auUel7+pEShqBCoiIoG3cNMBhsxcz5lfuq45FE9dxMv7vkCNp9ChHTBvJOz+1ry2XwjXvwQN2/g8lvLy9vtbR89FREQCzOE0GD93c4lgARS+Nn7uZhzO4neU932BGg8ABXmw/FmYfokZdGKrwNUT4N6lQQ06vlDYERERCbA1u44UWyo6kwEcyMphza4jFXqfw2mQvvMwX2zYR/rOwx5DUHnHw9618OrlsGQCOHKheXcYugou+TvEhO/OmPAdmYiISJQ4eNxzsCjtPl/e58v+G5/Hk3scFj8Ja14DDKhaB3o9A61vAUvFN0kHmmZ2REREAqx+jcRy3eft+3YfOsWQmetLzNZkZOUwZOZ6Fm46UP7xbF0A0zrBmlcBAy4aAEPXwoW3RkTQAYUdERGRgOuYUptkWyKeooEFcxamY0ptn99nT0rggzV7fNp/481zL0g6Tad1o+GD/pC9D2o1gztmw40zoFodD+8MTwo7IiIiARZjtTCuTypAiYDhuh7XJ7VEfRtv3jegY1Mysn3bf1Pac6046R+zhM+NUVg3fw6WGOg6Aoakm3t0IpDCjoiISBD0apXM9NvbYbcVX0Ky2xJLPeZd1vua1fWuns2Z+3TcPfdsy34+qzKRSXGvE5efDclt4C9LoceTEF/Vq88JR9qgLCIiEiS9WiXTI9Xuc8Xi0t6XvvOwV5/tbp+O67lrd2RgWz+N87b9C6szD+KqQvfHoOP/hfUpK29F/m8gIiISQWKsFro0933Pi6f3ufbfZGTluN23Y8GcBTpzP1Dhc/etpfOi4fDbFvOFc9Kg93NQ6yyfxxiutIwlIiISwcq7H4icbJh/H7xxtRl0qtaFfm/AwE+jKuiAwo6IiEjE83k/0E/zzePka18HDGgzEIathdY3R8xxcl9oGUtERCSE/NXN3Kv9QNkHYMGDsGWOeV0rBfq8CGdf7qffJjwp7IiIiISIv7uZe9wP5HTCurfgmycgN/v34+TD4fKHIK5KBX6DyKCwIyIiEgKeuo67qh6Xt5t5Cb9thbkjYE+6ed2wndmd3N6q4s+OEAo7IiIiQVZW13ELZtXjHqn2ci1pAVCQC98+B9/+E5z5EFcNrnocOt4L1pgKjD7yKOyIiIgEmbddx59ftI2u59R1u4+n1L0+v6TD3OFwaJt53eJq8zh5zSYB+o3Cm8KOiIiIn3i72djbruMvL93By0t3lNjH42mvz1M9G5O2b7q5PwegWj24ZjJccFNUnrLylsKOiIiIH/iy2djbruMuRffxAG72+hi0Ob6c1l+8DZZj5ktt7zDbPFR1X0ywMgnrOjsOh4OxY8eSkpJClSpVaN68OU899RSG8cf/iA3D4PHHHyc5OZkqVaqQlpbG9u3bQzhqERGJdA6nQfrOw3yxYR/pOw8X6xjujmuz8ZlLU66QsnDTgWKvl9V1/EyuT39izo88Maf4Xh87h3k17jmmx79AA8sx9liScdwxB254WUHnd2E9szN58mSmT5/OO++8wwUXXMD333/P3Xffjc1mY/jw4QBMmTKFqVOn8s4775CSksLYsWPp2bMnmzdvJjHRt+QsIiLi63Hw8mw2dlU9HjJzPRZw+153z8rIzi28tuBkYMxiHor9kBqW0+QbMcxw9OHlgr68TSu6+PJLR7mwntn57rvvuOGGG+jduzfNmjXj5ptv5uqrr2bNmjWAOavzwgsv8Nhjj3HDDTdw4YUX8u6777J//34+//zz0A5eREQijq8zNOD9ZuM1u44Ue91T1WNvtLD8yqfx45kQ9xY1LKf5wXkO1+U9zT8LbiWXeK/3BFUWYR12LrnkEhYvXsy2beZu8v/+97+sWLGCa665BoBdu3aRkZFBWlpa4XtsNhudOnUiPT09JGMWEZHIVNYMDZgzNGcuaXkbLNzd16tVMise6s4H93Zm2JXNy3xGAnmMiv2E+fFjaG/dzgkjkXH5g+iX9wRbjaaF9/m6JyjahfUy1sMPP0x2djbnn38+MTExOBwOnn76aQYOHAhARkYGAA0aNCj2vgYNGhT+zJ3c3Fxyc/+YCszOzg7A6EVEJJL4MkNTtEqxt8HC032uqscdU2rz2fp9HruXd7JsYXLCGzRjPwCLHO15PP8uDvDHWMrqcF5ZhfXMzscff8z777/PrFmzWL9+Pe+88w7/+Mc/eOeddyr03EmTJmGz2Qr/NWlSOesOiIjIH8o7Q1PWZmML5p6fsgKIp+7lSZzgmdjX+CjhKZqxn5zEegzJG8lf8keXCDrgocN5JRfWYeeBBx7g4Ycfpn///rRu3Zo77riDUaNGMWnSJADsdjsAmZmZxd6XmZlZ+DN3xowZQ1ZWVuG/vXv3Bu6XEBGRsOLppFV5Z2g8hZSi194GkOL7eAyuta5iccID9I9dat7Q/m4SR3zPDX/6K3Zb8Z5WHjucS3gvY506dQqrtXgei4mJwel0ApCSkoLdbmfx4sW0adMGMJekVq9ezZAhQzw+NyEhgYSEhICNW0REwlNpJ616pNpJtiV6XEYqbYnIFVLOfLa9HE09e7VKpkejArI+G0HtXxcDYNQ9F0ufF+GsS36/p2bZHc6lUFiHnT59+vD000/TtGlTLrjgAn744Qeee+457rnnHgAsFgsjR45kwoQJtGjRovDoecOGDenbt29oBy8iImHFm8abno6DezND06tVcsUDiNMBa18nZvGT1M47AdY4uHQ0lkvvg9ji/0e6xw7nUoLFKFqhL8wcP36csWPHMnv2bA4ePEjDhg0ZMGAAjz/+OPHx8YB5/HzcuHG8+uqrHDt2jG7duvHKK69w7rnnev052dnZ2Gw2srKySEpKCtSvIyIiIeJwGnSbvMTjBmTXrM2Kh7qzaHOGT3V2/CbzR5gzHPZ9b1436QR9XoT6LQP3mRHO2+/vsA47waKwIyIS3dJ3HmbAa6vKvO+DezvTpXkdr3tc+UV+DiyfAitfBGcBJCRB2jhofw9Yw3prbch5+/0d1stYIiIi/uDrSaugLRHtWg5zR8KRneb1+dfBtc9CUsPAf3YlorAjIiJRr6K1cPzu1BFYNBZ+mGleV7dD739Ayz7B+fxKRmFHRESinqsWTnlOWvmVYcCP/4YFD8HJ38zXOgw2l60SbYH97EpMYUdERKJeaY03va2FU+F9PMf2wvz7YPtX5nXd8+D6qdC0s6+/jvhIYUdERCqFitTC8bUTejFOB6x5FRY/BfknISYeLr0fuo0scZxcAkOnsdBpLBGRysTXGRpP9Xlc7yi1anHGRvM4+f715nXTLuZx8nrnVfj3EJ3GEhERccuXk1ZldUK3YHZC75FqLx6Y8k/DssmwcioYDvM4eY/x0O4uHScPAYUdERERD8rVCf3n/5jHyY/uMq9bXm8eJ6/huWejBJbCjoiIiAc+1ec5dQS+ehT+O8t8sUZD8zj5+b0DOELxhsKOiIiIB97V3TFIPfQVfP00nDoEWODiP8NVj0Oi9oGGA4UdERERD8qqz9PY8hvPJr5NixU/mC/Ua2keJ2/SMajjlNJpl5SIiIgHrvo88MfpK4AYHPw5Zj5fxz9IF+MH8zj5lY/B/y1X0AlDmtkRERG/CWoDzSA5sz7PBZbdTIp7jQutv29APqureZy8bovQDlQ8UtgRERG/qFDhvTDXq1UyPVokkfHFEyRveQOr4cBItGHp8RS0vcPjcfJoDH+RSGFHREQqzFPhvYysHIbMXF964b1IsHMJMXNH0ujYL+b1BTdi6TUZajTw+JZoDn+RRnt2RESkQsoqvAdm4T2Hs2IF+x1Og/Sdh/liwz7Sdx52+zxv7vHJycPw7/+D926EY79AUiMY8CHc8naZQWfIzPUlavS4wt/CTQcqNi7xiWZ2RESkQspVeM9H3syS+HUmxTDgfx/DV2Pg1GHAAp3+D7o/Bgk1Sn1ruasuS8Ao7IiISIX4VHivHLxZIgPc3nMgK4e/zlzP4K7NSEu1e7dn5sgumD8adi4xr+tfYB4nb9zBq/EGI/yJbxR2RESkQrwrvOf9fUV5M0vyxJwfAYvbe1zeWLmbN1buLn2mx1EAq16BpROh4DTEJMDlD0LXERAT5/WYAx3+xHfasyMiIhXiKrznab7Egrmc1DGlts/P9maWJCM7l4xs74KDxz0z+zfAa1fCorFm0Gl2KfwtHS6736egA4ENf1I+CjsiIlIhngrvFb0e1ye1XPtT/D37UWLDdN5Js5/Va1dCxv8gsSZc/zIMmgt1mpfrMwIZ/qR8FHZERKTCXIX37LbisxV2W2KFjp0HYvbDtWdm68rP4ZXOkP4yGE5o1Q+GrYV2d4Cl/BuHAxn+pHwshmFU8Fxe5MvOzsZms5GVlUVSkpq2iYiUl7+L6DmcBt0mL/HYm8oCNEhKACxkZru/50y1yWZs3HvcGLPSfMHWBHo/B+deXe5xuqM6O4Hn7fe3wg4KOyIi4cx1GgsoFmZcEaroaawz7ynOoJ/1Wx6Lm0ktywkMixVLp7/ClY9CQvVADF0VlANMYccHCjsiIuGtvHV2XJpaMpkY+zrdYn4EYJulGc3veZOYJu2D8wtIQCjs+EBhR0Qk/HkzS+K6Z9HmDN5cuZs4CrgnZgGjYj8l0ZJPjhHHCwX9aHvbo/S8sGmIfhPxF2+/v1VnR0REIkKM1VJmET7XPV2a1yHNto+6Sx/hXGM3ACscF/Bilb8x+Jar6Kk9M5WKwo6IiESX3BOwdCKXrJ4OhpP8+JpsbPUgMam38eHZddzumdHemuimsCMiIhHJbUDZ+Q3MGw1Ze8ybWt9KXM+JtKtez+NzdGoq+insiIhIxDkzoNQli4lV3+dq5wrzBltTuO55aJFW5nPK6rulwBP5FHZERCSiFA8oBrfELOPR2Pep6TyJw7Cw59y7SLnlaYivVupz1J288lDYERGRiFE0oDSzHGBi7BtcErMZgE3OZozJ/zOH9qSyIrYqMWU8S93JKw+FHRERiRhrdh3ht6wT/C1mPsNj/02iJZ/TRjzPF/TjDce1OIgBLwOKupNXHgo7IiISMfL2rGFu/GO0tJobkJc7WvNowT3sNRoUu8+bgKLu5JWHwo6IiIS/3OOwZAKXrf4XFqvBEaM6T+XfwWxnN0q22/QuoLi6k5fWd8uu7uRRQV3PRUQkvG1dCNM6w+oZWDBYYL2ctNx/MNt5KWcGHQvmsXFvAoq6k1ceCjsiIhKejmfCJ3fBB7dB9q9Q8yy4/d9YbvoXR0nyS0Dp1SqZ6be3w24rPhNktyXq2HkUUW8s1BtLRCSsGAasfxcWjYWcLLDEQJehcMXDhcfJ/V0IUBWUI5MagfpAYUdEJEwc2gFzR8AvvxcHTL4Irn/J/H/PoIAiagQqIiKRoyAPvnsRlj0LjlyIqwpXPgqd/gox7r+qvGkMKgIKOyIiEmp718Lc4XDQLA5I86vguuegVrOQDkuih8KOiIiERk42LH4S1r4OGFC1LvSaBK1vAYuWo8R/FHZERCT4fvoS5t8Hx/eb1xf9CXo+DVVV00b8T2FHRESC53gGfPkAbJljXtdqBte9AM2vDOWoJMop7IiISOA5nbD+HVg0DnJ/P05+yd/h8ocgvmqoR1cuOg0WOcK+qOC+ffu4/fbbqVOnDlWqVKF169Z8//33hT83DIPHH3+c5ORkqlSpQlpaGtu3bw/hiEVEpJjftsLbvWHeSDPoNGwH/7cMeoz3GHQcToP0nYf5YsM+0ncexuEMryopCzcdoNvkJQx4bRUjPtzAgNdW0W3yEhZuOhDqoYkbYT2zc/ToUbp27cqVV17JggULqFevHtu3b6dWrVqF90yZMoWpU6fyzjvvkJKSwtixY+nZsyebN28mMVHN20REQqYgF1Y8D9/+Exx5EFcNuj8Gnf4PrDEe3+bvgoH+tnDTAYbMXF+in1ZGVg5DZq5X5eUwFNZFBR9++GFWrlzJt99+6/bnhmHQsGFD7rvvPu6//34AsrKyaNCgAW+//Tb9+/f36nNUVFBExM/2rII5w+HQVvP6nB7mcfKaTUt9m6cg4VocCnWQcDgNuk1eUiyIFeVqHrrioe5a0goCb7+/w3oZa86cOXTo0IFbbrmF+vXr07ZtW1577bXCn+/atYuMjAzS0tIKX7PZbHTq1In09PRQDFlEpHLLyYJ5o+DNnmbQqVYP+r0BAz8pM+g4nAbj525224Hc9dr4uZtDuqS1ZtcRj0EHzHEeyMphza4jwRuUlCmsw87PP//M9OnTadGiBV999RVDhgxh+PDhvPPOOwBkZGQA0KBBg2Lva9CgQeHP3MnNzSU7O7vYPxERqaAtc2FaJ/j+TfO67e0wdA20vtmrujmRECQOHvc8vvLcJ8ER1nt2nE4nHTp0YOLEiQC0bduWTZs2MWPGDAYNGlTu506aNInx48f7a5giIpVb9n7zOPlP88zr2mebx8nPvtynx0RCkKhfw7u9oN7eJ8ER1jM7ycnJpKamFnutZcuW7NmzBwC73Q5AZmZmsXsyMzMLf+bOmDFjyMrKKvy3d+9eP49cRKQScDrN6sfTOplBxxoLl94HQ74rNeh4OmkVCUGiY0ptkm2JeJqnsmBupu6YouKI4SSsZ3a6du3K1q1bi722bds2zjrrLABSUlKw2+0sXryYNm3aAOZmpdWrVzNkyBCPz01ISCAhISFg4xYRiXoHfzL7We1dbV43ag99poK9ValvK+2kVY9UO8m2RDKyctzu23Ft/g1lkIixWhjXJ5UhM9djgWLjdAWgcX1StTk5zIT1zM6oUaNYtWoVEydOZMeOHcyaNYtXX32VoUOHAmCxWBg5ciQTJkxgzpw5bNy4kTvvvJOGDRvSt2/f0A5eRCQaFeTC0okwo5sZdOKrwzVTYPAir4LOkJnrS+zLcR3ZXrQ5g3F9zNn8M6NCOAWJXq2SmX57O+y24jNMdltiyE+LiXthffQcYN68eYwZM4bt27eTkpLC6NGjuffeewt/bhgG48aN49VXX+XYsWN069aNV155hXPPPdfrz9DRcxERL/zyHcwdAYe2mdfn9oLe/wRb4zLf6suR7UWbM8K6zo6LKiiHnrff32EfdoJBYUdEpBSnj8E342Dd2+Z1tfpw7RRI7evVKSuH0+Dtlbt4av6WMu8dduU5dD2nLu3PqsW6X44qSEipFHZ8oLAjIuKGYcDmL2DBg3Di94Mg7QaZbR6q1Cr9vb9zt0fHG+E4kyPhJyqKCoqISIhk/Qof/gk+GWQGnTrnwF1fwvVTfQo67vboeMO1j0e9psQfwvo0loiIBJnTAWvfgMXjIe8EWOOg2yjzSHmc90e+S6uG7A0Dcx/P+Lmb6ZFq1xKWVIjCjoiImDI3m8fJf11rXjfuaM7k1G/p86PKqobsjaIVk7s0r1OhZ0nlprAjIlLZ5efA8mdh5QvgLID4GpA2DjoMBmv5djv4s8qxWi9IRSnsiIhUZrtXmMfJD+8wr8/rDdc+C7ZGFXqsP6scq/WCVJTCjohIZXT6KHw9Fn54z7yubjdDTur1QMVryLjaKpRWDblBUgJgITM7fCsmS3RQ2BERqUwMA378Nyx4GE4eNF9rfzekPQFVagKlt3Tw9ii4N20Vnrj+AgC1XpCAU50dVGdHRKJHqTMyx/bC/Ptg+1fmdd1zzX5WZ3UpfL/ruPiZXwyuuOFrOwRvgpM/wpVUTioq6AOFHRGJBp5Cw9hrz6Xl3o9psuEfxBacwrDGYbnsfvNIeewfTZF9aengy2yLN0tiar0g5eHt97eWsUREooCnGRlb9jYafvYAKdadAKx1nss/44ZxV92r6XVG0Hl75a5Sj4uX9yh4jNVS5v3e3CNSXgo7IiIRzl0BvwTyGB77b/4SM584i4NsowqTCwYwy9Ed8qysnrm+cEnK15YOOgoukUZhR0Qkwp1ZwK+L9Ucmxr5OitXsZ7XAcTHj8u/iIH+0eXBVJ3Y6YeiskjNCpdFRcIk0CjsiIhHONdNi4wSPxr7PrbHLAMgwavF4/l187by4xHtcS1KPfbHJ66Cjo+ASqRR2REQiXP3qCfSxfsfjce9Sz5KN07DwvuMqphT05zhVS33vkZN5Xn2GjoJLJFPYERGJEG5PLGXvpfOq++gS/zUA25yNeDj/XtYb5/r1s+06Ci4RTGFHRCQCnLmJ2IqTEdUWM5QPiXWcxmmN44XcG5jhuJ48L/6r3QLUqhbHkZP5Zd47tndL7uqaohkdiVjl6/AmIiJB4zpW7go6qZbdzI5/nBGOt4h1nOZI3Q5Yh3xH6oAJ1LFVL/N5rsgy4YZWJNsS8RRhLJh1ehR0JNJpZkdEJIwVPVaeSC4jYv/NvTHzibU4yTaqMqngTyzL7sW3dVrQq56FHqn2YktdR0/m8dT84sfKiy5JWa0Wr9s1nLmM1v6sWqz75agKAUrYU9gREQljrmPlXa0bmRj7BmdZzX5W8x0deSJ/EL9RC7LzCgv9uSvO17OV3WN14l6tkpl+e7sSdXbO3KPjrhaP1QLOIglJLR4kXCnsiIiEsaOHM/hH3AxujlkOwH6jNo/n3803zvbF7iut0F9Z1Yl7tUouMSNUNBB5qs7sPOOFjKwchhQpVigSLhR2RETCkWHAxk/osfgh4mKO4DQsvOvowT8KbuWEm+PkvhT689SHyl0gcled2eOQ+aNYYY9Uu5a0JGwo7IiIhJuju2HeaNi5mDhgp6Up9+cO5gejRYlbfS3052uH8TOrM5elvP2zRAJJp7FERMKFowC+ewle6QI7F0NMAnR/jB03zmOD0aLEqSlfC/2dearLxbX8tHDTgRLvKW8fLPXPknCimR0RkXCwfwPMHQ4H/mteN7sUrnsB6p5DT2C6Nb7MTcSlKW05qrTlp/L2wVL/LAknPoedQYMGMXjwYC677LJAjEdEpHLJOwX/mQjpr4DhgEQbXD0B2t4Blj9CR1mbiMtS1nKUp+Wnjim1SbYlkpGV49W+HfXPknDkc9jJysoiLS2Ns846i7vvvptBgwbRqFGjQIxNRCS67VgM80bBsV/M6wtugl7PQI0Gbm8v61RVabxdVjrzvhirhXF9Ut3W4jmT+mdJuPJ5z87nn3/Ovn37GDJkCB999BHNmjXjmmuu4dNPPyU/v+yy4yIild7JQ/Dvv8DMm8ygk9QY/vQx3PKWx6BTUd4uK7m7z1WLx24r/rMz84zdlqhj5xKWLIZheDMz6dH69et56623eP3116levTq33347f/vb32jRouSpgXCVnZ2NzWYjKyuLpKSkUA9HRKJM4VHv7NOk/raAc36YiOX0EcACnf4K3R+FhBoBH0O3yUs8Lke5lp9WPNTd46yMKihLuPH2+7tCG5QPHDjAokWLWLRoETExMVx77bVs3LiR1NRUpkyZwqhRoyryeBGRiOc66h2b/QtPx75Ji5iNAGTbziPplunQuH0ZT/CP0pajvF1+creMpuPlEgl8XsbKz8/ns88+47rrruOss87ik08+YeTIkezfv5933nmHb775ho8//pgnn3wyEOMVEYkYCzcdYNjMtfQ58Qlfxz/EZTEbyTHimJzfn/aZj7LwWMOgjsfTcpSWnyTa+Tyzk5ycjNPpZMCAAaxZs4Y2bdqUuOfKK6+kZs2afhieiEhkcjgNPvxiLl/Ev8wFVnMD8neOVB4pGMxuIzlklYYreqpLJBL5HHaef/55brnlFhITPW92q1mzJrt27arQwEREIlbeSTJnj+WNvDeJsRocM6oxoeB2PnVchmvRKJSVhityqkskEvkcdu64445AjENEJDps/wbmjaJh1h6wwBeOS3gy/w4OY3N7uyoNiwSeKiiLiPjDid/gqzGw8RMAcqo14q9HB/IfZ5tS36ZKwyKBp95YIiIVYRjww/sw7WIz6Fis0GUYcX9fzdYanUv0s3KxYDbgVKVhkcBT2BERKa/DO+Hd6+GLv8Hpo2BvDX9eDD2fJiaxBuP6pAJUuIGniFSMwo6IiK8c+fDtczD9Eti1HGITIW083LsUGrUrvE1HvUXCg/bsiIj4Yt86mDMCMs3igJx9BVz3PNQ+2+3tOuotEnoKOyIi3sg9AUsmwJp/geGEKrWh50S4qH+x7uTueHPU+8xWDApEIv6jsCMiUpZtX8P80ZC117xufSv0mgTV6vrl8a6WEgey/jiGnmxLZFyfVC11ifiB9uyIiHhy4iB8cjfMusUMOjWbwu2fQb/X/Bp0hsxcXyzoAGRk5TBk5noWbjrgl88Rqcw0syMicibDgB/eg68fg5ys34+TD4UrxkB8Nb99jMNpMH7uZrddyA0IWUsJkWijsCMiUtThnTB3BOz+1ry2XwjXT4WGbf3+UWt2HSkxo1OUq6XE2yt3cVfXFAUekXLSMpaICJjHyZf/A17pYgad2Cpw9QTzOHkAgg543yriqflb6DZ5iZa0RMoposLOM888g8ViYeTIkYWv5eTkMHToUOrUqUP16tXp168fmZmZoRukiESeX7+Hf10GS54CRy407w5DV8Elf4eYwE2A+9IqQnt4RMovYsLO2rVr+de//sWFF15Y7PVRo0Yxd+5cPvnkE5YtW8b+/fu56aabQjRKEYkoucfhywfh9TQ4uBmq1oEbX4Xb/w21mpW43eE0SN95mC827CN952EcTsPta97qmFKbZFuix5YSRbmeOn7u5mKfUZHPF6ksImLPzokTJxg4cCCvvfYaEyZMKHw9KyuLN954g1mzZtG9e3cA3nrrLVq2bMmqVavo3LlzqIYsIuFu6wKYfx9k7zOvLxoAVz8N1dzXw3F3PLxm1TgAjp3KL3zNmyPjRWvq9L+4KS98sw0LuN2oXJRrD8+aXUfo0ryOjqyLeCkiws7QoUPp3bs3aWlpxcLOunXryM/PJy0trfC1888/n6ZNm5Kenu4x7OTm5pKbm1t4nZ2dHbjBi0h4OZ4BCx6CzZ+b17WamRWQm3f3+BbX8fAzw0jRkOPiWm7y1A7C29BUmoPHczyOqazPF6mMwn4Z68MPP2T9+vVMmjSpxM8yMjKIj4+nZs2axV5v0KABGRkZHp85adIkbDZb4b8mTZr4e9giEm6cTlj3Nrzc0Qw6lhjoOgKGpHsMOg6nwcrth3j4s41lzrq4eFpuAs81dbJO5XPsVD43t2vk1WfUrZZQ6pF1T58vUlmFddjZu3cvI0aM4P333ycx0fuNfGUZM2YMWVlZhf/27t3rt2eLSBj6bRu8c515pDw3C5LbwF+WQo8nIb6q27cs3HSAbpOXMPCN1Rw77d2Mi0vR5SYXb2rqrNhxCHuS5z08FsxlKix4dWS96OeLVGZhHXbWrVvHwYMHadeuHbGxscTGxrJs2TKmTp1KbGwsDRo0IC8vj2PHjhV7X2ZmJna73eNzExISSEpKKvZPRKJQQR78ZzLM6Aq/rIS4qmY/qz8vhuSLPL7N0wyMr4oeLfempk5Gdi4DOjYFKBF4XNfj+qRy6EQu3vD2aLtItAvrsHPVVVexceNGNmzYUPivQ4cODBw4sPD/HxcXx+LFiwvfs3XrVvbs2UOXLl1COHIRCbk9q+Ffl8J/JoIjD85Jg7+tMishl3KcvLQZGF8VPVrubfBoVrcq029vh91WfDbbbkss3Ifj7ZF1X462i0SzsN6gXKNGDVq1alXstWrVqlGnTp3C1wcPHszo0aOpXbs2SUlJ/P3vf6dLly46iSVSWeVkweInYe0bgAFV68I1k6FVvzK7k0PZMzDesGCGk44ptQtf8yWgdGlehx6pdo9d0F1H1jOyctyGMnefL1KZhXXY8cbzzz+P1WqlX79+5Obm0rNnT1555ZVQD0tEQmHLPPjyfjj+e+G9NgPNKshVvf/S98fSjwFc28oMK66Q4mtAibFa6NLc/TH4GKuFcX1SGTJzfYkj60WXu9ReQsRkMQyj0m/Xz87OxmazkZWVpf07IpEo+4AZcn6aZ17XPhuuewHOvtznR6XvPMyA11Z5da+7I+NWCxQ9BFW07o1rLxC4Dyi+HhdXnR2p7Lz9/lbYQWFHJGI5nbDuLfjmCcjNBmssXDIcLn8Q4qqU65EOp0G3yUs8zsAA1KwSx7SB7eh8tjnzsmbXERZtzuDNlbtL3HtmkPF3QClaoPDM5S6RaKew4wOFHZEIdPAn8yj5XnMW5mit1uzp+gyt2nWt8Je9rzMwroDkaa+Pa4lqxUPdibFaFFBE/MTb7++wPo0lIlJCQS4snQQzusHeVZwikSfy76T9gYe44dMsv3QH79UqucwTUUV5c6y8aN0b136cG9o0okvzOgo6IgEW8RuURaQS+SUd5g6HQ9sAWOxoy9j8u9lP3cJb/NUuoVer5FJPRBXl7aZm1b0RCQ2FHREJf6ePmfty1r0FgFGtPmNz72BmTjvOLL/nqkY8fu5meqTa3YYTb5eRSjsRVZTq3oiEN4UdEQlfhgFb5sCXD8KJ3/vdtbuTteeMYua7Wzy/jeLdwYsKxAkm1b0RCW/asyMi4SlrH8aHf4KP74QTGZxOSsFx5zy4/iUO5CV49Ygzl408tYFwLX2Vd6+Pq+4NlN7mQXtzREJDYUdEQsLhNEjfeZgvNuwjfefhPzp0O52w5jUKXroYy9YvyTdimFrQlzYHx9HtozwWbjpQrmWjshpxQsU6hfu6qVlEgkfLWCISdJ6Wkp69LI5uW56CX9cQC6x3nsPD+feyzWgC/DEDM+1PbUtdNgKzFo7TMHA4DWKsFp9OTHmzT8cdXzY1i0jwKOyISFC5lpKKhpQE8hhw8hM6fj0HLA5OUoVn8m/jfUcaziIT0K7Nx0/N38LY3qkMnVWyXYLLsdP5DHx9deF+nNwCp1fjq+iJKW83NYtI8GgZS0SCxt1SUkfLFr6MH8Pw2NnEWxz8hw5clTOF9xxXFws6Lq4ZmFrV4t0uG53JNRu0+9BJr8aoE1Mi0UczOyISNEWXkpI4wcOxH/Cn2KUAHDRq8nj+XSx0XkzJbb4lHTyeww1tGtEj1c6qnYcZOms9x07nl7jPNRv0wZo92JMSyczWiSmRykYzOyISNOYSkcG11lUsTnigMOi8X3AVabnPstDZEW+CDvwxAxNjtWC1WtwGHRcDyMjOZUDHpuDmE1xLYf0vbsK8/+0vvmFaRCKeZnZEJGgaW4/wetw/SIv5AYCdzmQezr+Xtcb5xe6rXS2eoyfzvJ6B8XafzeGTuYxMO5cP1uwhI/uP99h+717+/DfbC19T93CR6KGwI1JJBbUZpdMBa1+n3eInscScIM+IYbrjBqYV3EAecYW3uYKMp83HnmrWeLvP5t30XwCwJyUwKq0FzepWY/ehU7zwzbYSwcpfbSdEJPQUdkQqoUBUEfYo80eYMxz2fY8FOFqnLbfuH8AOo7HHINOrVTLTre1KjNHuYYxlVTAuMaTsXF74ZjvT/tSWD9fu8Vh7p6y2EyISGSyGYVT6hWlvW8SLRAN3R7/hj7Dht5mM/BxYPgVWvgjOAkhIgrRx0P4eFm7O9Cps+TL75Pq9wP1R9DNZgFrV4jhy0vNeH5cP7u2s4+QiYcjb72+FHRR2pPJwOA26TV7isbieaxlpxUPdKzaTsetbmDsCjuw0r8+/Dq59FpIaFhuLv5fR3M1Y+cOL/dtwQ5tGfn2miFSct9/fWsYSqUQCXkX41BFYNBZ+mGle10g2Q07LPiVuLU/xvbICUtEKxgs2HSjco1NRqr0jEtkUdkQqEW9PLflcRdgwYNNnsPBhOPmb+VqHweayVaLNx1G65+0+o6Ihypuw4+vJLxGJPKqzI1KJlKeBZpmO7YVZt8Fng82gU/c8uOcruO45vwYdX7uVuzYte1oYs2CGpQk3tCq8PvPnoG7lItFAYUekEvE2AHg1k+F0wKrpMK0TbP8KYuLhikfgr99C085+G3N5u5XHWC2M65MKlB5krr1Q3cpFop2WsUQqEVcAGDLT+xo2bmVsNI+T7zdPP9G0C/R5Eeqd5/cxV2SfUa9WZpAp6wi7upWLRDeFHZFKxtsA4Fb+aVg2GVZOBcMBCTboMR7aDQJrYCaKK7rPyNsgo27lItFLYUekEirXTMbP/4G5I+HoLvO65fXmSasa9oCO1R/7jBRkRCo3hR2RKORNDRuvA8CpI/DVo/DfWeZ1jYbQ+x9wfu8AjLyksqoj68SUiJRFYUckyvitFYRhwMZPzePkpw4BFrj4z3DV45DoW/HNihQQ9Ns+IxGptFRBGVVQlujht1YQR3+B+aNhxzfmdb2WcP1UaNKxXGPyR/gKaj8vEYkIahfhA4UdiQZ+aQXhKIDVM2Dp05B/yjxOftmD0HUExMb7PCZfwpc3sz9B7dQuImFP7SJEKpkKt4I48F/zOPmBDeb1Wd2gzwtQt0W5xlNWfZyiHcUXbc7wuTqyiIi3VFRQJEqU+4h23in4eiy8eqUZdBJt0GcqDJpb7qAD3oevl5fs8Lk6soiILzSzIxIlynVEe+cSmDcKju42ry+4EXpNhhoNKjweb8PXWyt3eTX7o+UqESkvhR2RKOHTEe2Th+GrR+B/H5o/TGoEvf8J513jt/F4G76Onc73+LMKd2EXEUHLWCJRw6teUNe1JGbjxzDt4t+DjgU6/RWGrvZr0AHv+nDVrBLn1bN87sIuIlKEwo5IFHG1gnDX1PLtG+vTa8NQmP0XOHUY6l8Af/4GrpkMCTX8PhZvwtfdXZt59SyfurCLiJxBy1giUaZEK4hqsXTK/AjroolQcBpiEuCKh+CS4RDj3cxKRcZSWh+uHql2Ply7V9WRRSSgVGcH1dmRwPNnfRifnrX/B/M4ecb/zOtml5rdyes0L+dvUj6ljdlViwfcV0f2uhCiiFQ6KiroA4UdCSR/Vv71+ll5J2HpRFj1ChhOSKwJV0+AtreDJfxONak6soiUh8KODxR2JFAq2r6h6IzI7kOneOGbbWU/a8c35nHyY3vMH7TqB72eger1/fAblT7GisxaqTqyiPhKFZRFQsyXCsLuvtTdzXa443rWi3PSuXrrAqwbPzZ/YGtiHic/t2cFfxPP/Dkjo+rIIhIoOo0lEiC+tG84k2tGqKyg43rSTdblzMr9uxl0LFbo/Df426qABx1VPhaRSKCZHZEAKW/7htJmhM7U1JLJxNjX6RbzIwBZSedhu206NGrv63B9UtFZKxGRYNLMjkiAlKt9A2XPCAHEUsBfY+bwdfyDdIv5kRwjjmfy+7Olz5yABx1vxljarJWISLBpZkckQHxq31BEWTNCF1p28kzc66RafwFgheMCHisYTG5SMx5oXvGeVt7wdtZq5Y7ftNFYREJOMzsiAeJV+4Y+qSWCgKcZoarkMDb2PWbHP06q9ReOGtW5L++v3JH/CL8YdrfPChRvZ61eXrqTbpOXaP+OiIRUWIedSZMmcfHFF1OjRg3q169P37592bp1a7F7cnJyGDp0KHXq1KF69er069ePzMzMEI1YpLjS2jd4OnburqfUFdYNfJ3wIINjFxBjMZjt6Epa7rN85rwMu61K0AvvldX3qihtWBaRUAvrOju9evWif//+XHzxxRQUFPDII4+wadMmNm/eTLVq1QAYMmQI8+fP5+2338ZmszFs2DCsVisrV670+nNUZ0cCzdcaMq6TTnXJYmzcu1wfkw7AXmc9Hiu4h3bdb6ZZ3WohrUfjqfKxO64luxUPddeSloj4TVQWFfztt9+oX78+y5Yt47LLLiMrK4t69eoxa9Ysbr75ZgB++uknWrZsSXp6Op07d/bquQo7EnYMg43zp3HW95NI4gQOw8Ibjmv5oOpAHro+fNoneFsLyOWDezurlo6I+E1UFhXMysoCoHZtc0PnunXryM/PJy0trfCe888/n6ZNm5YadnJzc8nNzS28zs7ODuCoRXx0eCfMHUHr3d8CcKL2BXx/4RO0btKewUGexTlzRqr9WbVY98vRwuseqXZ6pNp5ftE2Xl66o8znebuxWUTEnyIm7DidTkaOHEnXrl1p1aoVABkZGcTHx1OzZs1i9zZo0ICMjAyPz5o0aRLjx48P5HArFZX59w9Hfh6/zp9M4/9NJcaZhxFbBcuVY6jeeShXxAT/P6ruZm2sFnAWmQt2VUvuek5dr8KOtxubRUT8KWLCztChQ9m0aRMrVqyo8LPGjBnD6NGjC6+zs7Np0qRJhZ9bGamBo3+kL1tIvf88yDmGeZx8uaM1L8YP4V5bd3qFKOi46+nlPOMF1+bjaX9qW65j9iIiwRDWp7Fchg0bxrx581i6dCmNGzcufN1ut5OXl8exY8eK3Z+ZmYndbvf4vISEBJKSkor9E9+pXYAf5B5n98xhdFrSn3OMXzhiVGdU3hDuzH+Y9dk1Q/J39KWCs+uep+ZvYWxv34/Zi4gEQ1iHHcMwGDZsGLNnz2bJkiWkpKQU+3n79u2Ji4tj8eLFha9t3bqVPXv20KVLl2APt1Ipq10AmO0CHGdOBcgfti7EmNaZZjvew2ox+MzRjaty/8Fs56WAJWR/R28qOBflqpZcq1q8z8fsRUSCIayXsYYOHcqsWbP44osvqFGjRuE+HJvNRpUqVbDZbAwePJjRo0dTu3ZtkpKS+Pvf/06XLl28Pokl5eNLuwCdvjnD8UxY+BD8OBsLsMdZj0cLBvOt88ISt3rzd/T3nqnybiI+eDyHG9o0okeqXXu4RCSshHXYmT59OgBXXHFFsdffeust7rrrLgCef/55rFYr/fr1Izc3l549e/LKK68EeaSVT3mbXFZqhgHr34VFYyEnCywxbG8+iD6bLiWHhFLf6unvGIg9U+XdROx6X4zVooArImElrMOONyWAEhMTmTZtGtOmTQvCiMSlvE0uK61D22HuSPjl9w32yRfB9S9x6FRjcjatKvPt7v6OnjYRu/ZMlXfpqKyeXmfS5mMRCXdhvWdHwldZ7QIsmDMMwfgCdDgN0nce5osN+0jfeTi89gkV5MGyZ2F6VzPoxFWFq5+GPy+B5IvK/XcM5J6p0np6uRsfaPOxiIQ3hR0pl/I2ufS3hZsO0G3yEga8tooRH25gwGurwqLxpMNpsHHVIrJf7AJLJ4AjF5pfBX9Lh0uGwe/Hycv7d/Rlz1R5eOrpdeb/OLX5WEQiQVgvY0l4c30hnrlnxB6kOjuBWsapqEU/bCdr7lhucizEajE4ZCTxUtw9dGnzf/Sq1bDE/eX5OwZjz1SvVsklNhufWUFZm49FJBJEVG+sQFFvrIpxdxoICOiJHIfToNvkJR5nN0LVeHL91zNJXjmWZIs5o/Kp4zIm5A8kixoApQYwX05Vpe88zIDXyt7ro15UIhLNorI3loSnM0/fBKOqctgdfT+egfHlg7Tb8gVYYLezAY8UDOY7Z6vCWyyY+2h6pNrdhhhfTjGVtYlYm4ZFRP6gPTviV8Gqqhw2R9+dTvj+LXi5I5YtX1BgWHml4Hp65k0uFnSg4vtoigqXPVMiIpFAYUf8JphVlcPi6Ptv2+Dt3jBvJORmcbRma/rkPc2Ugv7kEu/xbf4KYJ42EWvTsIhIcVrGEr8J5tJSSJdxCnJhxfPw7T/BkQdx1eCqsfxUpx9b3lhb5tv9GcDcbSIu7/4oda8XkWilsCN+E8ylJdcyzpCZ67FAscAT0GWcPatgznA4tNW8bnE19P4n1GxKR6cRkgDmj4rF6l4vItFMy1jiN8FeWgrqMk5OFswbBW/2NINOtXpw85vwp4+hZlMgcvfRqHu9iEQ7zeyI34RiacmfyzgebZkLXz4Ax3//0m97B/R4EqqW/D1CXXvIV2XtsyrrBJmISCRQ2BG/CdXSUsAaT2bvN0POT/PM69rNoc+LkHJpqW/zNoCFwx6ZsDvCLyISAAo74leRNrPhltMJ696Eb8ZDbjZYY6HrSLjsAYjzbgmurAAW7D0ynoJV2BzhFxEJIIUd8bugLC0FysGfYO5w2LsaAKNRB/7X9kl2xzaj/p6TdExJqPDvEew2F6UFq7A4wi8iEmAKOxIQAVtaCpSCXPMo+bfPgTMf4quzOXUE925uw75PjwEbgIrPvgR7j0xZwWran9qqErOIRD2dxpKw5nAapO88zBcb9pG+8zAOp+H2tQr55TuY3hWWTTaDzrm9+E/aXHqvuoB92fnFbq3oCaVAdysvypsij0/N38LY3pF3gkxExBea2ZGw5W75pWbVOACOnfojhJR7tuX0MfhmHKx727yuVh+unYLj/BsYM2VpQGZfgrlHxttgVatafOTvsxIRKYXCjoQlT8svRUOOi897XQwDNn8BCx6EE5nma+0GQY/xUKUWa3Ye9umEki+nqoK5R8aXYHVDm0aRu89KRKQMCjsSdkpbfnHHp9mWrF/N4+RbvzSv65xjHidv1q3wFm9DwoJNB1iz6wgfrNlDRrZ3p6qCWYvI12AVcfusRES8pD07EnbKWn5xp8y9Lk4HrH4VpnUyg441Di57EP66sljQAe9Dwrvpv/D8N9uKBR0ofV9PMKssu4KVpydZMIOZNh+LSLRT2JGwU5H9Km7fm7nZbPOw4AHIOwGNO8Jfv4Xuj7qtm1NWSChLWR3eg9XmIlLbV4iI+JuWsQIkkNVxw6HybiBVZL9Ksffm58DyZ2HlC+AsgPgakDYOOgwGq+ecX1olaG+VVXk4WLWIoqLIo4hIBSnsBEAgq+NWhu7UZe1rcafEXpfdK2DuCDi8w7w+rzdc+yzYGnn1PE8hwVelzVIFa49MRBd5FBHxA4thGBUsUhL5srOzsdlsZGVlkZSUVKFneTpF5Ppa8XWZougszu5Dp3jhm21+e7avnx/ML0nX3xHKnlkp9vs3T4Svx8IP75kvVrebIadlH7D4Pm7X779g0wHeTf/F5/d/cG9nbfoVEQkQb7+/NbPjR/6ujutuFsedQHWnDuUskqeZFXd1duy2RMZd15JepMPLD8HJg+YP2t8NaU9AlZrlHkfR2Rdfwo4qD4uIhA+FHT/yZwdpTzNE/ni2N4Ldv8kdT8svQPHXap8iZsEo2LbQfGPdc83j5GddUurzfZm18mVpTZt/RUTCi8KOH/mrOq6vdWa8fba3X+7B7t9UGk/7Wro0r2MeJ1/7Onz0pHnKyhoHl94Hl46G2IRSn+vrrJUvm5a1+VdEJLwo7PiRv6rjlqfOTFnP9uXL3Z8zVAGTscnsTr5vnXndpBP0mQr1z3d7uzd7n8qatfJ4sikpgQEdm9KsbjVt/hURCUMKO37kr+q45akzU9qzfV2S8qWCMBDcL/f807BsCnw31TxOnpBkHidvf4/H4+S+7H0CeGT2Rk7nO7EnlQwuOtkkIhJ5FHb8qLSlDl/2cfhaZ6a0Z5dnScqXCsLvpv8SvKPvPy+DeSPhyM/m9fnXmSetkhp6fIuve58AjpzMZ9RHGwD3s19qqyAiEllUQdnP/FEd19cKvqU925clqfJ+fmntEfzi1BH4fCi8e70ZdGokw20zof/7pQadiux9cgn47yYiIgGnmZ0AqOhSR1kzRAYwKq2FV3tEyrNp2tcKwgHbtGwYsOkzWPAQnDpkvtZhsLlslWgr8+0V2ftUOASCuyFbRET8T2EnQCq61OGvMv/l3TTtawVhv29aPrYH5o2GHYvM67rnwfVToWlnrx9RkR5bRYXFhmwRESk3hZ0w5o/NsBXZNF30872tIFzhgOF0wOp/wZIJkH8SYuLh0vuh28gyj5OfqSI9ttzxV3gSEZHgUtgJcxWdIaropmlfKwiXFjDKrPOTsRHm/B32/2BeN73ELA5Y79wyP9ed8vTYKo2/w5OIiASHwk6UKC1I+GNJrKLH6kut83NeTfjPM/DdS2A4IMEGPcZDu0Gldicvi7d7n5rWrspT87dw9GRehUoGiIhIeFIjUPzbCDQUvC0YWNGmnp6ac5bViLS05qiXWDfxaq2ZVDu5x3wx9Qa4ZgrUsHs9Lm/GXdbfp7y/m4iIhI63398KO0R22PF3l3VvPs+XNgsOp0G3yUtKbHKuyXEei3ufm2OWA2DUaIil9z/h/Gv9NtYzx1FW0Atl41MREfGdwo4Pwins+DL74ilIuLiWX1Y81N2vR6Z9GWP6zsMMeG1VkVcMrrd+x+Nx71HXko3TsPCuowepA5+lY8tmfhtjeVV09ktERILH2+9v7dkJI77OLISqh5Uvm6aLnmBqbDnI07FvcnnM/wDY6mzMmPw/s944lxdz4/w2vopQdWQRkeijsBMmfO1fBf7rsh5I9WskEoODu2MWMjr2U6pacsk14phacCOvOq4j//f/FdRJJxERCRSFnTBQnv5V4L8u64HUMXEP86uM43zD7Ge1ytmSMfl/ZpdhBjeddBIRkUBT2AkD5V2O8leX9YDIOwn/mURM+iucbzjIMqoysWAgHzmuwLV92l2dH+2ZERERf1PYCQPlXY7yV5d1v9uxGOaNgmO/FyG84EbWpYxm+aLfoJQ6PzoNJSIigaCwEwYqshzlrx5afnHyMHw1Bv73kXmd1Biuew7O7Ul3YEU7z7M25dmzJCIi4o2oCTvTpk3j2WefJSMjg4suuoiXXnqJjh07hnpYXqnocpQ/emhViGGYAWfhGDh9xBxxp79C98cgoXrhbZ5OOpV3z5KIiIg3yl+LP4x89NFHjB49mnHjxrF+/XouuugievbsycGDB0M9NK+4lqPgj+UnF2+Xo1xB4oY2jejSvE7wQsGRXfDejTD7/8yg06AV/HkxXPNMsaBTGl/2LImIiPgqKsLOc889x7333svdd99NamoqM2bMoGrVqrz55puhHprXXMtRdlvxpSq7LTE8l3AcBbDyRXilC/y8FGIT4apx8Jf/QOP23j3CaZC+8zALNh3w6n51HRcRkfKI+GWsvLw81q1bx5gxYwpfs1qtpKWlkZ6e7vY9ubm55ObmFl5nZ2cHfJzeCPlylLf2/wBzhkOGWRyQlMvguhegTnOvH+FuM3JZVItHRETKI+LDzqFDh3A4HDRo0KDY6w0aNOCnn35y+55JkyYxfvz4YAzPZ2FdwTfvJCydCKteAcMJiTWh50Ro8yew+N5Q1Ns+JarFIyIiFREVy1i+GjNmDFlZWYX/9u7dG+ohhb8d38ArnSH9ZTPotLoZhn0PbQf6FHRK24zsTkiP0IuISFSI+JmdunXrEhMTQ2ZmZrHXMzMzsdvtbt+TkJBAQkJCMIYX+U78Zh4n3/iJeW1rah4nb9GjXI8razPymUJyhF5ERKJKxIed+Ph42rdvz+LFi+nbty8ATqeTxYsXM2zYsNAOLpIZBvz3A/jqETh9FCxW6DQErnzE61NW7ni7yfjOLmdxTavk8NyzJCIiESXiww7A6NGjGTRoEB06dKBjx4688MILnDx5krvvvjvUQ4tMh3fCvJGwa7l5bW8NfaZCo3YVfrS3m4yvaZUcvnuXREQkokRF2Lntttv47bffePzxx8nIyKBNmzYsXLiwxKZlKYMjH757CZZNhoIc8zj5FWOgy1CIifPLR4R1Py8REYlKFsMwvN0rGrWys7Ox2WxkZWWRlJQUsM8J6yaXv66DucMhc5N5ffYVcN3zUPtsv3+U6zQWuO/nFZZ1hUREJOx4+/0dFTM7kSBsm1zmHoclT8PqGYABVWpDr0lw4W0+nbLyRTD6eYV1sBQRkaDSzA6Bn9nxVFcmkDMZ3nzZO35aSMGckSScMisYO1vfirXXJKhW169jqcgYyyNsg6WIiPiVt9/fCjsENuw4nAbdJi/xeNzatUdlxUPd/TbzUOaX/YmDHPhwBMm/fgnAHmc9Hi0YzI4aHSM+EIQiWIqISGh4+/1dKYsKBlOwm1y6vuzP/MyMrByGzFzHxrkvkf9ie5J//RKHYeFfBb3pmTeZb50X/n7PehZ62asq3JTVPR3M7ukOZ6XP9yIilYrCToB5W1fGH00uS/uyb2Y5wKy4p2m97jHi8rPZ5GzG9XkTmFQwkNOYx8EjPRCoe7qIiLijDcoB5m1dGX80uXT3ZR9HAX+Jmcfw2NkkWPI5bcTzXMHNvOm4BgcxJZ5RNBBEWp2bYAZLERGJHAo7ARbMujJnfom3tWxnUtzrnG81e38td7TmkYLB/GrU9/lZkSCYwVJERCKHlrECLMZqYVyfVOCPTbIu/m5y6foSr8Zpnoh9m8/in+B8614OGzUYkfc37sx/2KugU/RZkcQVLD39JS2YG7VVsFBEpHJR2AkCV10Zu614gLDbEv16OqhjSm1uqb6RRQkPcFfs11gtBp85LiUt91m+cHbDggV7UgL2pOgMBMEMliIiEjl09JwoqaB8PBMWPAibPwfgF2d9HikYzEpna6D40WsgqisYq86OiEjloDo7PghW2AkIpxN+eBe+fhxys8ASw8/n3sM9u7qzO/uP/9Ge+WUf7YFAFZRFRKKfwo4PQhV2KvyF/Ns2szv5LyvN64Ztze7kyRd6V0FZgUBERCKYemOFuQrNrBTkwcoXYPmz4MiDuKrQ/THo9FewmsfJY6yWMo+Oe3OPiIhIpFPYCQFPLQ1cFYxL3TOzZ7XZnfy3n8zrc9Kg93NQ6yy/jE2zPSIiEm0UdoKsrJYGFswKxj1S7cVDRk42LB4Pa98w76xaF66ZDK36+a07ebTv4xERkcpJR8+DrFwtDX6aD9M6wdrXzTva3A7D1kLrm/0adDz31IrcflkiIiKa2Qkyn1oaZB+ABQ/Alrnmi7XPhutegLMv9+uYyj3bJCIiEgEUdoLMm8rEFpy0PvAZLHgWcrPBGguX/B0ufwjiqvh9TL7MNmlDs4iIRBqFnSArq1fWOZZ9/CPxTc5evcV8oWE7uP4lsLcK2JjUQFNERKKZ9uwEmaeWBvHkMzL2U76Mf5g2xhaIqwa9JsOfvwlo0AE10BQRkeimsBMCZ/bK6mD5ifnxjzAy9t/EWxzQoicMXQ2d/6ibE0hqoCkiItFMy1gh0qtVMj3OrsJvn4/Bvm0WAEa1+liumQwX3Oi3U1becM02DZm5Hgvu+2WpgaaIiEQqzeyEyuY5xLzSqTDo0O5OLMPWQKubghp0XILVmV1ERCTYNLMTbNn74csH4Kd55nWdc8zj5CmXhnRY8PtsU6pdFZRFRCSqKOwEi9MJ378B34yHvOPmcfKuI+GyByAufDb+ql+WiIhEG4WdYDi4BeYMh1/XmNeNOsD1U6HBBaEdl4iISCWgsBNI+Tnw7T9hxfPgzIf46nDVOLh4cFBOWYmIiIjCTuCcPgqv94DD283r866Fa58FW+PQjktERKSSUdgJlCq1oH5Ls93DNVMg9YaQnLISERGp7BR2Aum6582NyFVqhnokIiIilZbCTiBVqxvqEYiIiFR6KiooIiIiUU1hR0RERKKawo6IiIhENYUdERERiWoKOyIiIhLVFHZEREQkqinsiIiISFRT2BEREZGoprAjIiIiUU1hR0RERKKawo6IiIhENYUdERERiWoKOyIiIhLV1PUcMAwDgOzs7BCPRERERLzl+t52fY97orADHD9+HIAmTZqEeCQiIiLiq+PHj2Oz2Tz+3GKUFYcqAafTyf79+6lRowYWi8Vvz83OzqZJkybs3buXpKQkvz1XitPfOXj0tw4O/Z2DQ3/n4Ajk39kwDI4fP07Dhg2xWj3vzNHMDmC1WmncuHHAnp+UlKT/IAWB/s7Bo791cOjvHBz6OwdHoP7Opc3ouGiDsoiIiEQ1hR0RERGJago7AZSQkMC4ceNISEgI9VCimv7OwaO/dXDo7xwc+jsHRzj8nbVBWURERKKaZnZEREQkqinsiIiISFRT2BEREZGoprAjIiIiUU1hJ4CmTZtGs2bNSExMpFOnTqxZsybUQ4oqkyZN4uKLL6ZGjRrUr1+fvn37snXr1lAPK+o988wzWCwWRo4cGeqhRJ19+/Zx++23U6dOHapUqULr1q35/vvvQz2sqOJwOBg7diwpKSlUqVKF5s2b89RTT5XZW0nKtnz5cvr06UPDhg2xWCx8/vnnxX5uGAaPP/44ycnJVKlShbS0NLZv3x6UsSnsBMhHH33E6NGjGTduHOvXr+eiiy6iZ8+eHDx4MNRDixrLli1j6NChrFq1ikWLFpGfn8/VV1/NyZMnQz20qLV27Vr+9a9/ceGFF4Z6KFHn6NGjdO3albi4OBYsWMDmzZv55z//Sa1atUI9tKgyefJkpk+fzssvv8yWLVuYPHkyU6ZM4aWXXgr10CLeyZMnueiii5g2bZrbn0+ZMoWpU6cyY8YMVq9eTbVq1ejZsyc5OTmBH5whAdGxY0dj6NChhdcOh8No2LChMWnSpBCOKrodPHjQAIxly5aFeihR6fjx40aLFi2MRYsWGZdffrkxYsSIUA8pqjz00ENGt27dQj2MqNe7d2/jnnvuKfbaTTfdZAwcODBEI4pOgDF79uzCa6fTadjtduPZZ58tfO3YsWNGQkKC8cEHHwR8PJrZCYC8vDzWrVtHWlpa4WtWq5W0tDTS09NDOLLolpWVBUDt2rVDPJLoNHToUHr37l3sf6/Ff+bMmUOHDh245ZZbqF+/Pm3btuW1114L9bCiziWXXMLixYvZtm0bAP/9739ZsWIF11xzTYhHFt127dpFRkZGsf/+sNlsdOrUKSjfi2oEGgCHDh3C4XDQoEGDYq83aNCAn376KUSjim5Op5ORI0fStWtXWrVqFerhRJ0PP/yQ9evXs3bt2lAPJWr9/PPPTJ8+ndGjR/PII4+wdu1ahg8fTnx8PIMGDQr18KLGww8/THZ2Nueffz4xMTE4HA6efvppBg4cGOqhRbWMjAwAt9+Lrp8FksKORIWhQ4eyadMmVqxYEeqhRJ29e/cyYsQIFi1aRGJiYqiHE7WcTicdOnRg4sSJALRt25ZNmzYxY8YMhR0/+vjjj3n//feZNWsWF1xwARs2bGDkyJE0bNhQf+copmWsAKhbty4xMTFkZmYWez0zMxO73R6iUUWvYcOGMW/ePJYuXUrjxo1DPZyos27dOg4ePEi7du2IjY0lNjaWZcuWMXXqVGJjY3E4HKEeYlRITk4mNTW12GstW7Zkz549IRpRdHrggQd4+OGH6d+/P61bt+aOO+5g1KhRTJo0KdRDi2qu775QfS8q7ARAfHw87du3Z/HixYWvOZ1OFi9eTJcuXUI4suhiGAbDhg1j9uzZLFmyhJSUlFAPKSpdddVVbNy4kQ0bNhT+69ChAwMHDmTDhg3ExMSEeohRoWvXriVKJ2zbto2zzjorRCOKTqdOncJqLf7VFxMTg9PpDNGIKoeUlBTsdnux78Xs7GxWr14dlO9FLWMFyOjRoxk0aBAdOnSgY8eOvPDCC5w8eZK777471EOLGkOHDmXWrFl88cUX1KhRo3Dd12azUaVKlRCPLnrUqFGjxD6oatWqUadOHe2P8qNRo0ZxySWXMHHiRG699VbWrFnDq6++yquvvhrqoUWVPn368PTTT9O0aVMuuOACfvjhB5577jnuueeeUA8t4p04cYIdO3YUXu/atYsNGzZQu3ZtmjZtysiRI5kwYQItWrQgJSWFsWPH0rBhQ/r27Rv4wQX8vFcl9tJLLxlNmzY14uPjjY4dOxqrVq0K9ZCiCuD231tvvRXqoUU9HT0PjLlz5xqtWrUyEhISjPPPP9949dVXQz2kqJOdnW2MGDHCaNq0qZGYmGicffbZxqOPPmrk5uaGemgRb+nSpW7/O3nQoEGGYZjHz8eOHWs0aNDASEhIMK666ipj69atQRmbxTBUNlJERESil/bsiIiISFRT2BEREZGoprAjIiIiUU1hR0RERKKawo6IiIhENYUdERERiWoKOyIiIhLVFHZEREQkqinsiEjUuuKKKxg5cmSohyEiIaawIyIiIlFN7SJEJCrdddddvPPOO8Ve27VrF82aNQvNgEQkZBR2RCQqZWVlcc0119CqVSuefPJJAOrVq0dMTEyIRyYiwRYb6gGIiASCzWYjPj6eqlWrYrfbQz0cEQkh7dkRERGRqKawIyIiIlFNYUdEolZ8fDwOhyPUwxCREFPYEZGo1axZM1avXs3u3bs5dOgQTqcz1EMSkRBQ2BGRqHX//fcTExNDamoq9erVY8+ePaEekoiEgI6ei4iISFTTzI6IiIhENYUdERERiWoKOyIiIhLVFHZEREQkqinsiIiISFRT2BEREZGoprAjIiIiUU1hR0RERKKawo6IiIhENYUdERERiWoKOyIiIhLVFHZEREQkqv0/Dr8W+2BNBAwAAAAASUVORK5CYII=",
+      "text/plain": [
+       "<Figure size 640x480 with 1 Axes>"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    }
+   ],
+   "source": [
+    "plt.plot(t, y,'o')\n",
+    "plt.plot(t, y_hat)\n",
+    "plt.title('t vs y')\n",
+    "plt.xlabel('t')\n",
+    "plt.ylabel('y')"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 77,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/plain": [
+       "Text(0, 0.5, 'y')"
+      ]
+     },
+     "execution_count": 77,
+     "metadata": {},
+     "output_type": "execute_result"
+    },
+    {
+     "data": {
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjsAAAHHCAYAAABZbpmkAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAABZcklEQVR4nO3deXiTVfrG8W+SriwNlKWLgFRAsSIiIIswrkUQRFBEYEBxGf0NggqMjjIjIoIijuMCKriLIuA2iojiIKAOyCaIUkE2QRFaUJYWCi1t8v7+iAktTdqkzd77c129hiZv3p5m1Nyc85znmAzDMBARERGJUuZQD0BEREQkkBR2REREJKop7IiIiEhUU9gRERGRqKawIyIiIlFNYUdERESimsKOiIiIRDWFHREREYlqCjsiIiIS1RR2REREJKop7IhI2Pj666956KGHOHz4cKiHIiJRRGFHRMLG119/zcSJExV2RMSvFHZEREQkqinsiEhYeOihh7j33nsByMjIwGQyYTKZ2LVrl9vrR40aRZ06dTh27Fi554YMGUJqaio2mw2Ab775hp49e9KwYUMSExPJyMjglltuqXA8w4cPp2HDhhQXF5d77oorruCss87y8TcUkVBR2BGRsHDttdcyZMgQAJ566inefPNN3nzzTRo1auT2+kGDBlFQUMDChQvLPH7s2DEWLFjAddddh8ViYf/+/VxxxRXs2rWL+++/n+nTpzN06FBWrVpV4XhuuOEGDhw4wGeffVbm8dzcXJYuXcqwYcOq8duKSDDFhHoAIiIAbdu2pX379sydO5f+/fvTvHnzCq/v3r07p512Gm+//TYDBw50Pb5w4UIKCgoYNGgQ4KgDOnToEP/973/p2LGj67rJkydXeP/LLruMJk2aMHv2bK666irX43PnzsVutyvsiEQQzeyISEQymUwMHDiQTz75hKNHj7oef/vttznttNPo3r07APXq1QPg448/drsk5YnZbGbo0KF89NFHHDlyxPX4W2+9xYUXXkhGRoZ/fhERCTiFHRGJWIMGDeL48eN89NFHABw9epRPPvmEgQMHYjKZALj44osZMGAAEydOpGHDhvTr14/XXnuNoqKiSu9/4403cvz4cT744AMAtmzZwrp167jhhhsC90uJiN8p7IhIxOrSpQvNmzfnnXfeAWDBggUcP37ctYQFjhmg9957j5UrVzJq1Cj27NnDLbfcQocOHcrMCLmTmZlJhw4dmD17NgCzZ88mLi6O66+/PnC/lIj4ncKOiIQN52yML66//noWLVpEfn4+b7/9Ns2bN6dLly7lruvSpQuPPPII33zzDW+99RY//PAD8+bNq/T+N954I0uXLiUnJ4c5c+bQp08f6tev7/M4RSR0FHZEJGzUrl0bwKemgoMGDaKoqIhZs2axaNGicrMuhw4dwjCMMo+1a9cOwKulrCFDhmAymbj77rv56aefVJgsEoFMxqn/FRARCZG1a9fSqVMnevfuzeDBg4mNjaVv376uEORJq1at2LdvH0eOHGHdunW0b9/e9dzTTz/N888/zzXXXEOLFi04cuQIL730Ejk5OWzYsMGrQuO+ffvy8ccfU69ePXJzc4mPj6/27yoiwaOZHREJGxdccAGTJk3iu+++46abbmLIkCH89ttvlb5u0KBBHDlyhJYtW5YJOuAoUO7YsSPz5s3jrrvu4vHHH6dVq1YsXbrU6x1VN954I+BYMlPQEYk8mtkREanE/Pnz6d+/P1999RV/+tOfQj0cEfGRwo6ISCWuuuoqNm/ezPbt26tURC0ioaUOyiIiHsybN4/vv/+ehQsX8swzzyjoiEQozeyIiHhgMpmoU6cOgwYNYubMmcTE6O+HIpFI/+aKiHigvwuKRAftxhIREZGoprAjIiIiUU3LWIDdbmfv3r3UrVtXBYgiIiIRwjAMjhw5Qnp6Omaz5/kbhR1g7969NG3aNNTDEBERkSrYvXs3TZo08fi8wg5Qt25dwPFmJSUlhXg0IiIi4o38/HyaNm3q+hz3RGGHkyctJyUlKeyIiIhEmMpKUFSgLCIiIlFNYUdERESimsKOiIiIRDWFHREREYlqCjsiIiIS1RR2REREJKop7IiIiEhUU9gRERGRqKawIyIiIlFNHZRFREQkIGx2gzU7D7L/SCGN6ybQKSMZizn4B24r7IiIiIjfLcrOYeKCTeTkFboeS7MmMKFvJr3apAV1LFrGEhEREb9alJ3DiNnrywQdgNy8QkbMXs+i7JygjkdhR0RERPzGZjeYuGAThpvnnI9NXLAJm93dFYGhsCMiIiJ+s2bnwXIzOqUZQE5eIWt2HgzamBR2RERExG/2H/EcdKpynT8o7IiIiIjfNK6b4Nfr/EFhR0RERPymU0YyadYEPG0wN+HYldUpIzloY1LYEREREb+xmE1M6JsJUC7wOL+f0DczqP12FHZERETEr3q1SWPGsPakWssuVaVaE5gxrH3Q++yoqaCIiIj4Xa82afTITFUHZREREYleFrOJri0ahHoYWsYSERGR6KawIyIiIlFNYUdERESimsKOiIiIRDWFHREREYlqCjsiIiIS1RR2REREJKop7IiIiEjg7N0A2z4P6RAUdkRERCQwDu6Et66DuYNCGngUdkRERMT/Cn6H2QOg4DdodDY07RSyoSjsiIiIiH+dKIA518PBHWBtBkPfhYSkkA0npGHnq6++om/fvqSnp2Mymfjwww/LPG8YBg8++CBpaWkkJiaSlZXFtm3bylxz8OBBhg4dSlJSEvXq1ePWW2/l6NGjQfwtRERExMVWAu/dAnvWQWJ9GPY+JAX3lPNThTTsFBQUcN555/Hcc8+5ff7xxx9n2rRpzJw5k9WrV1O7dm169uxJYWGh65qhQ4fyww8/sHjxYj7++GO++uorbr/99mD9CiIiIuJkGLBwDGxdBDEJMORtaHRmqEeFyTAMI9SDADCZTHzwwQf0798fcMzqpKen87e//Y177rkHgLy8PFJSUnj99dcZPHgwmzdvJjMzk7Vr19KxY0cAFi1aRO/evfn1119JT0/36mfn5+djtVrJy8sjKSl002wiIiIRbdkU+PIxMJlh0Gxo3SegP87bz++wrdnZuXMnubm5ZGVluR6zWq107tyZlStXArBy5Urq1avnCjoAWVlZmM1mVq9e7fHeRUVF5Ofnl/kSERGRavjmVUfQAej9RMCDji9iQj0AT3JzcwFISUkp83hKSorrudzcXBo3blzm+ZiYGJKTk13XuDNlyhQmTpzo5xGLiIjUHDa7wZqdB9l/pJCzDn3FWV/9DRPAxffBBbeGenhlhG3YCaRx48YxduxY1/f5+fk0bdo0hCMSERGJHIuyc5i4YBM5eYV0NP3I7LgpmEx2fmk+kD1Nb2f/hj00rptAp4xkLGZTqIcbvmEnNTUVgH379pGWdrKKe9++fbRr1851zf79+8u8rqSkhIMHD7pe7058fDzx8fH+H7SIiEiUW5Sdw4jZ6zGAVqZfeSXuCRJMxSy2deCvP16N7ceTZSRp1gQm9M2kV5savBurIhkZGaSmprJkyRLXY/n5+axevZquXbsC0LVrVw4fPsy6detc1yxduhS73U7nzp2DPmYREZFoZrMbTFywCQNI4wCz4h7DajrGN/YzubN4FDYsZa7PzStkxOz1LMrOCc2A/xDSmZ2jR4+yfft21/c7d+5kw4YNJCcn06xZM0aPHs3kyZNp1aoVGRkZjB8/nvT0dNeOrbPPPptevXpx2223MXPmTIqLixk1ahSDBw/2eieWiIiIeGfNzoPk5BVSjyO8EfcY6aaDbLOfxq0n7qGQ8ismBmACJi7YRI/M1JAtaYU07HzzzTdceumlru+ddTTDhw/n9ddf5+9//zsFBQXcfvvtHD58mO7du7No0SISEhJcr3nrrbcYNWoUl19+OWazmQEDBjBt2rSg/y4iIiLRbv+RQhIp5NW4f9HKvIe9RjI3nrifPOp4fI0B5OQVsmbnQbq2aBC8wZYSNn12Qkl9dkREJBqV3jHlj4LhVdtyOfbGIC6zbOCwUZvrTkxgu9HEq9c+M7gd/dqdVuWf7Y63n99hW6AsIiIiVVd6x5STtwXDbkMSBp02TsBs2cBxI45bTtzrddABaFw3ofKLAkRhR0REJMqU3jFVmrNgeMaw9h4Dj6eQNOf0j8nYOg+7ycKoE3fxreHdMRAmINXqCEyhEra7sURERMR3pXdMncr52MQFm7DZy1/hDEmlgw7A1UffJWPrqwCYr57OwD//hVRr5TM1zgWzCX0zQ9pvRzM7IiIiUcS5Y8oTZ8HwU4u30q1lQ1cdj6eQNNDyBeNi5wLwFMPoVLcnPc5oQI/M1DJLXYcKTjBpYdkZodQw6bOjsCMiIhJF9h/xHHRKe3bZdp5dtt1Vx2NNjCsXkrLM63gs5iUAZpZcxTMlveHl1R5rf3q2SfVrQbS/aBlLREQkivhaCOys4/l8U9kzJTuZNvNc7DQsJoN3Sy7isZIh5V5zarNAi9lE1xYN6NfuNLq2aBAWQQcUdkRERKJKp4xk0qwJeBsznMtWH2zY43os07SLl+OeIP6PYyDuL7kNSt2xstqfcKOwIyIiEkUsZhMT+mYC+BR4DhYUk1w7jgxTDrPiHiPJdJzV9taMKr6z3DEQztc4mwWGO4UdERGRKNOrTRozhrX3asdUacMyY3kj9jEamfLZZD+dv5y4hyLiKnyNtzVCoaQCZRERkSjUq02aa8fUiu2/8eyyHRVeb+Uof939ELXMv7HblMaNJ+7nCLUq/TmhbBboLc3siIiIRClnwfCYHmdVWMdTi0JmJ/6bWoe3Qt000kd9yjO3XkG9xFiP9zbhaDYYymaB3lLYERERiXIV1fHEUcKM2Kc519gCCfVg2H+wNMigW6uGPDbgXExuXhMuzQK9pbAjIiJSA7ir4zFjZ0atmVxs+R5ia8HQdyEls8LXgKNZYEVHToQbnXqOTj0XEZGaw3XIZ/5xuvzwMCnb3wZzLAx9B1pcVvFrwqxZoE49FxERkXIsZhNdz0iGxQ/C9rfBZIbrXvEYdFyvadEgiKP0Ly1jiYiI1DTLn4Kvpzn+3HcaZPYL7XgCTGFHRESkJln7MiyZ6PjzFZOh/Q2hHU8QKOyIiIjUFN+/Awvvcfz5T/fAhXeGdjxBorAjIiJSE/z4CXzwV8CATrfDZQ+EekRBo7AjIiIS7X76Et69CQwbnDcEek0FU+h3UwWLwo6IiEg0+/UbmDsEbEXQ+iq4+lkw16yP/5r124qIiNQkudkwewAUF8AZl8J1r4Kl5nWdUdgRERGJRr9vgzf7Q+FhaNIJBr8FMfGhHlVIKOyIiIhEm0M/wxv9oOA3SG3rOAYirnaoRxUyCjsiIiLR5EiuI+jk74GGZ8ENH0BivVCPKqRq3sKdiIhICFX1nCmvXldwwBF0Du2EeqfDjR9C7YYBGU8kUdgREREJkkXZOUxcsImcvELXY2nWBCb0zazwBHGvXnf8MMy+Bn77EeqmwfCPICk9IOOJNFrGEhERCYJF2TmMmL2+TLAAyM0rZMTs9SzKzqn664qOwlsDIec7qNUQbpwP9ZsHZDyRSGFHREQkwGx2g4kLNmG4ec752MQFm7DZy17hzese++hbjLmD4dc1kFDPsXTV6KyAjCdSKeyIiIgE2JqdB8vNoJRmADl5hazZedCn18VSzITjj2Ha9T+Iq4tt6PusLEhn/oY9rNxxwGNYqep4IpVqdkRERAJs/xHPwaKi6yp6XQwlTIt9lkst31FiSWBd1xmMnp1HTt4q1zWe6m+qOp5IpZkdERGRAGtcN6FK13l6nRk7T8bOoJdlLUVGDAtaP8Hgz8xe199UdTyRSmFHREQkwDplJJNmTcDThm4TjlmYThnJlb7OhJ2pMS9ytWUlxYaFf8Tey9Rt6T7V31R1PJFKYUdERCTALGYTE/pmApQLGM7vJ/TNLNffpvzrDCbFvMbAmK8oMczcVXwnzbpcS26+b/U3VR1PpFLYERERCYJebdKYMaw9qdayS0Op1gRmDGvvsa+N63VJ8YyPmc2wmCXYDROTYu+i35//SvOG3h0DcWr9TVXHE4lUoCwiIhIkvdqk0SMz1eeOxb3OSeWKnBcwr/gUgJ8unMKDPf6KxWxi5Y4DXv1sd/U3VR1PpFHYERERCSKL2UTXFg18e9EXj2Fe8ZTjz72foGWn21xPOetvcvMK3dbtmHDM1niqv6nSeCKMlrFERETC2VdPwJePOf7ccwqUCjpQ8+pvqkJhR0REJFytmAZLJzn+nDURut7h9rKaVH9TFVrGEhERCSGPp46vmgmLxzsuuvQB6D66wvvUlPqbqlDYERERCRFPp46/nLmBc7592PHARffCxfd6db+aUH9TFQo7IiIiIeA8dfzUouJLjy7knG9fcXzT7W649J9BH1u0Uc2OiIhIkHk6dXyg5QsejXUEnTmWq7Fd9hCYtAxVXQo7IiIiQebu1PEB5q+YGvMSAK+W9OIfBYN46vNtFZ5eLt7RMpaIiIifeCw2PsWp3Yz7m5fzr9gXMJsMZpX04OGSGwATzy7bzrPLtrs9vdzbnyUKOyIiIn7hqdj41JACZbsZ9zMv59+xMzCbDN4quZwJJTdxascc5+nlzm3kvvws0TKWiIhIOTa7wcodB5i/YY9Xy0jOYuNTl6acIWVRdk6Zx51dj682f82TsTOwmAzmlFzKAyU3U741YNnTyz/53refJWEedmw2G+PHjycjI4PExERatGjBpEmTMIyT/9AZhsGDDz5IWloaiYmJZGVlsW3bthCOWkREItmi7By6T13KkJdWcfe8DQx5aRXdpy71GCI8FRtD2ZBSOjBZzCZmtNvFU7HPYTEZzC25lH+W3IpRwcey8/TyB+Zn+/SzJMzDztSpU5kxYwbPPvssmzdvZurUqTz++ONMnz7ddc3jjz/OtGnTmDlzJqtXr6Z27dr07NmTwkLPx92LiIi44+sMDbgvNi7NGVLW7Dx48sHs/9Buzb1YTAYLzJfzj0qCTmkHC0749rMkvGt2vv76a/r160efPn0AaN68OXPnzmXNmjWAY1bn6aef5oEHHqBfv34AvPHGG6SkpPDhhx8yePDgkI1dREQiS2UzNCYcsyY9MlPLFAKfWmzsieu67P/A+38BwwbthtK773Qa7jrMiu2/8eyyHdX+PXwZU00R1jM7F154IUuWLGHr1q0AfPfddyxfvpwrr7wSgJ07d5Kbm0tWVpbrNVarlc6dO7Ny5UqP9y0qKiI/P7/Ml4iI1GxVmqGhbLFxRRrXTYDs908GnfP+DFdPx2Kx0LVFA8b0OIs0a4Kbih0HE5BcO9b7nyUuYR127r//fgYPHkzr1q2JjY3l/PPPZ/To0QwdOhSA3NxcAFJSUsq8LiUlxfWcO1OmTMFqtbq+mjZtGrhfQkREIoLPMzR/cBYbVxRS0qwJdD72Bbx/2x8zOsOg37Ngtriu8+b08sn92nj1szplJHv1u9QUYR123nnnHd566y3mzJnD+vXrmTVrFk888QSzZs2q1n3HjRtHXl6e62v37t1+GrGIiIQ7TzutfJqhKcWbkDKj3S7M/ykVdK6eXiboOFV2ennvtumV/qwJfTPVb+cUYV2zc++997pmdwDOPfdcfv75Z6ZMmcLw4cNJTU0FYN++faSlnewrsG/fPtq1a+fxvvHx8cTHxwd07CIiEn4q6k/TIzOVNGsCuXmFbut2TDhCh7tZE2dIOfXeqdYEZrTbSbs194Jhh/OHQd/pYPY811DZ6eUV/Sz12XEvrMPOsWPHMJ/yD4TFYsFutwOQkZFBamoqS5YscYWb/Px8Vq9ezYgRI4I9XBERCWOeDt4s3bBvQt9MRsxejwnKXOfNrIm7kNL56OeYP/Q+6DhVdnp5ZYFIygrrsNO3b18eeeQRmjVrxjnnnMO3337Lk08+yS233AKAyWRi9OjRTJ48mVatWpGRkcH48eNJT0+nf//+oR28iIiEDW93Wi2/77JqzZqUCSkb5sL8OxxBp/2NcNUzXgUdb1UWiOSksA4706dPZ/z48dxxxx3s37+f9PR0/u///o8HH3zQdc3f//53CgoKuP322zl8+DDdu3dn0aJFJCSoEl1ERBx82Wnll1mTb2fD/FGOO3e4Gfo86degI74xGaXbEddQ+fn5WK1W8vLySEpKCvVwRETEz+Zv2MPd8zZUet0zg9vRr91p1fth69+Aj+4CDLjgL9D7CTBpeSkQvP38VswUEZGoV9WdVj775lX46E7AgE7/p6ATJhR2REQk6nnbC6da/WlWvwAfj3H8ucsdcOVUBZ0wobAjIiJRz5teONXqT/P1s/Dp3x1/vvAu6Pmogk4YUdgREZEaobKGfZXttPLUjJD/PQn//afjz3+6B3o8rKATZsJ6N5aIiIg/VXWnldtmhEnxvNHyC1ptmu544JJxcPF9CjphSGFHRERqFF/707hvRmgw7NgbtNo03/HtZePhonv8OUzxIy1jiYiIeOC+GaHBP2PeYmSMI+hMswzH1v1voRieeElhR0RExINTmxGasPNwzOvcFvMJAOOLb+LJgp6s2XkwVEMUL2gZS0RExIP9R04GHTN2Hol5hSExy7AbJsaV/IW3bZeWu07Cj8KOiIiIB84mgxZsPBE7k2ssK7AZJu4p/isf2P9U7joJTwo7IiIiHnTKSKZpkoV/HH+aKy1rKDYsjC4eyUJ7F8DRoye1us0IJeAUdkRERDyw2Ap5P/l5GueuociI4Y7iu1li7wD4qRmhBIXCjoiI+I3NblTvtPBwUnQU5g2hce5X2CwJ3Gv+O0uKWrueTrUmMKFvZqXNCCX0FHZERMQv3Dbei9RAcPwwzLkedq+GuDpY/vwOTzW7kCE+BrmoCn8RzGQYhlH5ZdHN2yPiRUTEPfeN904u9XhzHEPYKPgd3rwGcr+HBCsM+w806ejzbaIq/IUpbz+/1WdHRESqxX3jPQfnYxMXbDp5llQ1fo7bs6l8vKZC+TnwWm9H0KnVEG5aWOWgM2L2+jJBByA3r5ARs9ezKDvH53tK1WkZS0REquXUxnunMoCcvELW7Dzo0zENpXkzS1LtmZRDP8MbV8OhXZB0Gtw4Hxq28nmslYU/E47w1yMzVUtaQaKZHRERqRZvG+pVtfGeN7Mknq7JySvkr7PXM2nBDxXP9Py2FV7t5Qg69ZvDzZ9WKeiAb+FPgkMzOyIiUi3eNtSrSuM9b2ZJHvroB8Dk9hqnV1bs4pUVu9zP9OR8B29eC8d+h0at4YYPIanqNTWBDn/iO83siIhItXTKSCbNmoCnBRkTjuWkqjTe82aWJDe/iNx874JDuZqZX1bB630dQSetHdz0SbWCDgQ2/EnVKOyIiEi1WMwmJvTNBCgXeKrbeM/fsx9lCqa3fg5v9IeiPDi9GwxfALWrVlNUWiDDn1SNwo6IiFRbrzZpzBjWnlRr2dmKVGtCtbadB2L2wwDOO/IVpnmDoeQ4tOwBQ9+DBP+0Hglk+JOqUZ8d1GdHRMRf/N1Ez2Y36D51Kbl5hW5rckxASlI8YGJfvvtrTjXQ8gWPxbyExWTAOdfANS9CTFyVx+iJ+uwEnref3wo7KOyIiIQz504roEyYKd2wEHB7zan+YlnIA7FvAbCv5SBS/jwDzBY/j/gkdVAOLDUVFBGRqODNEpmna04yuCfmbVfQmW3pT8MhMwMadMCxpNW1RQP6tTuNri0aKOiEiGZ20MyOiEgk8GaWxHnN4k25vLpi1x+zP3YejnmdG2I+B2Bq8WDOG/KQlpKigLef3+qzIyIiEcE5S+LNNV1bNKBTRjKPfPQd9x6fztWWldgNE/+KvZ3zrhujoFPDKOyIiEhU6nVmEj2bvohpx0rsphi2d/8391w2XEtJNZDCjoiIRKQKl7WOHYQ512P6dS3E1sI86E3ObJlVtXtJxFPYERGRiFPhtu5mhuP4h982Q0I9Rw+dphdU7V5a7ooKKlBGBcoiIpHEuRX91A8vE5BhymFh/X+TeGwv1E2DGz6AxmdX6V5AtRoiSuBp67mIiESdig4GbWP6iXfjHiLx2F6M5BZwy2cVBp3KDhmFP46V8HRSukQMhR0REYkYng4G7WbeyNy4yTQwHeF7ewbfXDYX6p9epXs5GUBOXiFrdh6s7rAlxBR2REQkYrg7GLSPeRWvxT5OHVMhy23nMOTEA+wtqVOle1XnOglfKlAWEZGIcerBoDdY/svEmFmYTQYf27owtngEJ4j16gBRbw8ZDcRhpBJcmtkREZGI0SkjmTRrAiYM/hbzDpNiX8dsMnijpAd3FY+imFjSrI6t497fyz0TeH0vCW8KOyIiEjEsZhMP9TmTx2Je4s6YDwF4onggD5bchPHHR9qEvple9cixmE1M6JsJUC7wOL/39l4S3hR2REQkcpw4Rs/sexgU8wU2zNxXfBvP2q4BTGUOBvWWN4eMSuRTzY6IiESGYwdhziD4dQ3EJMCAV+kf14ULq9n1uFebNHpkpqqDchRT2BERkfB36GeYPQAObHN0Rf7z21iadaGrn27vzSGjErkUdkREJLzlfA9vXQdH90FSExj2PjRuHepRSQRR2BERkfD10xcwbxicOAKNz4Fh70FSeqhHJRFGYUdERMLT9+/ChyPAXgzN/wSD34IEa6hHJRFIu7FERCS8GAYsfwr+8xdH0DnnGsfSlYKOVJFmdkREJHzYbfDp32Hty47vu4yEKyaDWX83l6pT2BERkfBw4hi8/xfYshAwQc9HoesdoR6VRza7oe3qESLso/KePXsYNmwYDRo0IDExkXPPPZdvvvnG9bxhGDz44IOkpaWRmJhIVlYW27ZtC+GIRUTEZwW/wxtXO4KOJR7bda+zsvEg5m/Yw8odB7DZjVCPsIxF2Tl0n7qUIS+t4u55Gxjy0iq6T13KouycUA9N3AjrmZ1Dhw7RrVs3Lr30Uj799FMaNWrEtm3bqF+/vuuaxx9/nGnTpjFr1iwyMjIYP348PXv2ZNOmTSQk6PA2EZGwd2CHo4fOoZ2QUI9VXZ5jzIJEcvJWuS5JsyYwoW9mWHQ0XpSdw4jZ6zk1fuXmFTJi9np1Xg5DJsMwwisul3L//fezYsUK/ve//7l93jAM0tPT+dvf/sY999wDQF5eHikpKbz++usMHjzYq5+Tn5+P1WolLy+PpKQkv41fREQq8ctqmDsYjh+Ees34X6eZ3PjR4XJBwrk4FOogYbMbdJ+6lJy8QrfPm3AcNbH8vsu0pBUE3n5+h/Uy1kcffUTHjh0ZOHAgjRs35vzzz+ell15yPb9z505yc3PJyspyPWa1WuncuTMrV670eN+ioiLy8/PLfImISJD98CHM6usIOunnY7vlc/7+ZWG5oAO4Hpu4YFNIl7TW7DzoMeiAY5w5eYWs2XkweIOSSoV12Pnpp5+YMWMGrVq14rPPPmPEiBHcddddzJo1C4Dc3FwAUlJSyrwuJSXF9Zw7U6ZMwWq1ur6aNm0auF9CRETKMgz4ejq8exPYiuDMK+Gmhaz5LSbsg8T+I57HV5XrJDjCOuzY7Xbat2/Po48+yvnnn8/tt9/ObbfdxsyZM6t133HjxpGXl+f62r17t59GLCIiFbKVwMK/wX8fAAy44DZHs8C42hERJBrX9a4W1NvrJDjCOuykpaWRmZlZ5rGzzz6bX375BYDU1FQA9u3bV+aaffv2uZ5zJz4+nqSkpDJfIiISYEVHYN4Q+OYVDExkn3MvK8+6H9sfH0WRECQ6ZSSTZk3AUzWOCUcxdaeM5GAOSyoR1mGnW7dubNmypcxjW7du5fTTTwcgIyOD1NRUlixZ4no+Pz+f1atX07Wrv87CFRGRasvfC69dCdv+SxFx/PXE3Vy17nyGvLzatWU7EoKExWxiQt9M13hKc34/oW+mipPDTFiHnTFjxrBq1SoeffRRtm/fzpw5c3jxxRcZOXIkACaTidGjRzN58mQ++ugjNm7cyI033kh6ejr9+/cP7eBFRMQhdyO8dDnkbuR3I4lBRQ/wmb3Tyaf/2LK9eFNuRASJXm3SmDGsPanWsjNMqdaEkO8WE/fCeus5wMcff8y4cePYtm0bGRkZjB07lttuu831vGEYTJgwgRdffJHDhw/TvXt3nn/+ec4880yvf4a2nouIBMjWz+C9W+DEUXaZTmNY4b38ajQud1npLduLN+UyccGmMsXK4dRnx0kdlEPP28/vsA87waCwIyISAKtfgEX3g2EnL6Urf/r5ZvKpU+FLRl3akm4tG9Lh9Pqs+/mQgoRUyNvP77DuoCwiIhHIboNF42DNC45v2w3jlfgR5P/8c6UvfXbZdp5dtt01k9Ov3WmBHq3UAGFdsyMiIhGm6AjMHeIKOlva/I1um65h2peVB53SnHU8OmtK/EFhR0RE/OPwL/BKT9j2GcQk8G2Xp+n1TQdy8ot8vlW4dEyW6KCwIyIi1ffrN44dV/t/gNqNsQ1fyB3fNnN79IO3wqFjskQHhR0REame7Pfhtd5QsB9S2sBtS1lzIqPCox98oaMXpLoUdkREpGoMA754zLG13FYEZ/aCWxZBvaZ+DSg6ekGqS7uxRESknEp7yBQfhw/vgB/+4/i+y0i4YhKYLYD3ASW5dhyHCk64Xe5y9t7R0QtSXQo7IiJSxqLsnIqb+h3Jdey42rsezDHQ50noMLzMPZxHP+TmFVYYZMb3yWTknPWYoMx14dQxWSKflrFERMRlUXYOI2avL1dv49wKvuJ/S+DFSx1BJ7E+3PBhuaAD3p8h1butjl6QwFMHZdRBWUSih7vlJ8CrYw1sdoPuU5d6LCzubV7Nv+NmkkgRNDwL/jwPks+ocDyVzhJVMG7N6Ehl1EFZRKSGcRcs6tWKBeDwsWLXY57CxusrdroNOibs3GX5gDGx7wNwKP1i6t/4JiRYKx1TrzZp9MhMrTTIWMwmurZo4NsvLOIlzeygmR0RiXzO5Sdv/oPujBnOZSJ3IckpkUKeiJ1JH8saAF4q6U3jAVPpd34z/w1epIo0syMiUkPY7AYTF2zyuoGfgSPwTFywCbsdRs5xH5LSOMBLcf+mjXkXJwwL/yy5lXdtlzA3qbb/Bi8SBAo7IiIRbs3Ogz438HN2J35gfrbboNPBtIWZcU/RyJTP70YSfz0xmnVGa9K0FVwikMKOiEiEq04Dv4MFJ8o9NsiyjEkxrxJnsrHJfjq3nRjLXhoB2goukUlhR0QkQnjaseSvDsMxlDA+5k2GxywG4GNbZ+4t/j+Ok+C2qFkkUijsiIhEgIq2cPfITK2wgZ87JqB+7VgOFjh2adUnn+dip3GhZRMATxQP5Flbf8DE+D5nc1O3DM3oSMRSU0ERkTBXWaO/xZtyPTbwc8d5zeR+bUizJpBp+pmP4sZzoWUTR40Ebjsxlmdt12DCRJo1QUFHIp7CjohIGKtop5XzsYkLNtEjM9VtJ+J6tWJdvXacnN2Je7dNZ0a7XbwfN4Gm5t/YaU/hmhMPs9je0e1xDTa7wcodB5i/YQ8rdxzAZq/xnUskQmgZS0QkjFW208q5q2rNzoMeG/g571Om1gc7LJ5Au9VPgwlWmdpx+4k7yKcO4AhEpWt03C2jpSbFM6RTM5o3rK2uxxLWFHZERMKYtzutnNd56kRc5rHjh+D9v8D2zx3fd7ubCy59kBd+znPb5dhTw8Lc/CKe+nyb63sVMUu48nkZa/jw4Xz11VeBGIuIiJzC251WXu/I2vcDvHiJI+jEJMKAV6DHw1hiYujaogH92p1G1xYNyixdeduw0FlDtCg7x7uxiASJz2EnLy+PrKwsWrVqxaOPPsqePXsCMS4REQE6ZSSTZk3wWHhsAu8b/WX/B17OgkO7KKzdhGXdZ7Oy1qUV1t740rCwdA2R6nkknPgcdj788EP27NnDiBEjePvtt2nevDlXXnkl7733HsXFxZXfQEREvGYxmzzutHJXROyWrQQWPwjv3QzFx1hjakuXA+O5eVERQ15aRfepSz3OxvjasLB0DZFIuKjSbqxGjRoxduxYvvvuO1avXk3Lli254YYbSE9PZ8yYMWzbtq3ym4iIiFd6tUlzu9PKuauqwhqZggMw+1pY8QwAL5RcxZDj93KYuq5LKlp+qmrDwup0dRbxt2oVKOfk5LB48WIWL16MxWKhd+/ebNy4kczMTB5//HHGjBnjr3GKiNRonnZaVTijs2c9vHMj5O3GiK3FeOOvzC7sWO6y0geD9shMLXNP5zKaLw0LoeohSSQQfJ7ZKS4u5v333+eqq67i9NNP591332X06NHs3buXWbNm8fnnn/POO+/w8MMPB2K8IiI1lnOn1alFxG6tfxNe7QV5uyG5Bd/1fJ/ZR8sHHSdPy08VLaO541MNkUiQ+Dyzk5aWht1uZ8iQIaxZs4Z27dqVu+bSSy+lXr16fhieiIj4pKQIPr0P1r3m+P6s3nDNTH7+8ShwqNKXu1t+ci6jndpn51Re1xCJBJnPYeepp55i4MCBJCR4nqKsV68eO3furNbARETER4d3O5at9q4HTHDZP6H738BspnHdEq9u4Wn56dRltF2/H2Puml/IzS/VZFB9diRMmQzDqPH7A/Pz87FareTl5ZGUlBTq4YhIlPF0WrlfbV/iaBR4/CAk1odrX4ZWWWXG0H3qUo+1NyYcYWX5fZd5Pbag/F4iFfD281sdlEVEAqii08r9MgNit8P//g3LHgEMSGsH178B9U8vc5mz9mbE7PWYoEzgqeryk6duzSLhRgeBiogESGWnlVe70/CxgzB3MCybDBjQfjjc8lm5oONUrS3sIhFMMzsiIgFQ2WnlnrZ6e23Penh3OBz+BWISoPcT0P6GSl9WpS3sIhFOYUdEJAB8Oa3cp6Ugw3DstPr0PrCdgPoZjmWrtLZe30LLT1LTKOyIiASAr6eVe+VEAXw8Fr6f5/i+9VXQ7zlIrOf7AEVqEIUdEZEA8Ptp5b9tcWwr/+1HMFkgawJceBeYtPwkUhmFHRGRAKjsmAXnVm+vOg1//w4suBuKj0GdVLjuVWjezd9DFola2o0lIhIAfjmtvLjQEXL+c5sj6GRcDH/9n4KOiI8UdkREAqRaW71/3w4vZ8G61wETXHw/3PAB1Gkc0DGLRCMtY4mIBFCVtnpvfM8xo3PiKNRqCNe+CC0vD96gRaKMwo6ISIB5s9XbZjf4ZtteGq98iIxd7zgePL07DHgZktTsT6Q6FHZEREJsUXYOr81fzISiJ8gw/4zdMDErZgDpHR6ip4KOSLUp7IiIhNCi7Bz+O+cZXo19ldrmIn43khhTfAfLi9rCnO+ZYY7RMQ4i1aQCZRGRELEVHsH2nxE8GTeD2qYivrZl0rtoCv+zt3VtV5+4YBM2u7vN6yLiLYUdEZFQyM2m6PmL6WNfhs0w8WTxdQwr/gf7qe+6xHmkxOsrdirwiFSDwo6ISDAZBqx5CV66jFr5O8g16vPnEw8wzXYtdg//SZ60cDPdpy6t/inpIjWUwo6ISLAcOwhvD4NP7gFbEYdOu4TeRVNYbZxd6Utz8woZMXu9Ao9IFURU2HnssccwmUyMHj3a9VhhYSEjR46kQYMG1KlThwEDBrBv377QDVJEopbNbrByxwHmb9jDyh0HfFta+vlrmPkn+PFjMMdCzykk3fIB8dbG5Tosu6MaHpGqi5jdWGvXruWFF16gbdu2ZR4fM2YMCxcu5N1338VqtTJq1CiuvfZaVqxYEaKRikg0WpSdw8QFm8jJO3lKeZo1gfF9zqZ+7XjPDQNtJfDVv+Crx8GwQ3ILx9lW6e2w4DgyYsTs9ZjA7RlapTlreNbsPOjq22OzG741LBSpgSIi7Bw9epShQ4fy0ksvMXnyZNfjeXl5vPLKK8yZM4fLLrsMgNdee42zzz6bVatW0aVLl1ANWUSiyKLsHEbMXl8ujOTkFXLHnG/LPJZmTWBC30zHdvHDv8D7t8HuVY4nzxsCvf+FLbYOa3YccAWU5/7cnkkLywapiuw/Uugal7sA5vr5IgJESNgZOXIkffr0ISsrq0zYWbduHcXFxWRlZbkea926Nc2aNWPlypUew05RURFFRUWu7/Pz8wM3eBGJaDa7wcQFmyqddXFy1tZ8cHEu7TZMhKI8iE+CPk9C24F/BJQ1bmeIcvIKmbRwc6U/o3HdBI8BzPnzKz17S6QGCfuanXnz5rF+/XqmTJlS7rnc3Fzi4uKoV69emcdTUlLIzc31eM8pU6ZgtVpdX02bNvX3sEUkCtjsBq+v2On1jAtALY4zNeYF2q0e6wg6TS5wnFT+R9AZMXt9ufvl5hUycs63pFkTSLMmeKzhMeEIRh1Or+8xgKm2R6S8sA47u3fv5u677+att94iISGh8hd4ady4ceTl5bm+du/e7bd7i0h0WJSdQ/epS72aaXFqZ9rOwrh/cH3Ml9gNE7+eOwpu/hTqN69whsj52KSFmxnfJxOgXOBxfj+hbybrfj5UYQArXdsjImEedtatW8f+/ftp3749MTExxMTE8OWXXzJt2jRiYmJISUnhxIkTHD58uMzr9u3bR2pqqsf7xsfHk5SUVOZLRMTJ0wyMJ2bsjLJ8wHtxD9HcvI89RgMGn3iAdS3uAEssAGt2HvQqoNSvHceMYe1JtZb9C16qNcG1NOWs2amMt9eJRLuwrtm5/PLL2bhxY5nHbr75Zlq3bs19991H06ZNiY2NZcmSJQwYMACALVu28Msvv9C1a9dQDFlEIpyvNTpNTL/x79gZdDb/CMACWxf+WXwr+dRmTN2TgcWXgNKv3Wn0yEz1uMuqcV3vZrq9vU4k2oV12Klbty5t2rQp81jt2rVp0KCB6/Fbb72VsWPHkpycTFJSEnfeeSddu3bVTiwRqZLKZmBOMrjGvJyHY1+nruk4R4xEHiy+iQ/s3QET9RJjsRsGNruBxWzyOaBYzCbX9vJTdcpIJs2aQG5eodtQZsIxE9QpI9mrnykS7cI67Hjjqaeewmw2M2DAAIqKiujZsyfPP/98qIclIhHKmxkYK0d5JPYVrrKsBmCt/UzGFo9gt5Hiuubw8WKGvrzatRW8R2aq3wKKxWzy2J+ndG2P+u2IOJgMw6jx5fr5+flYrVby8vJUvyNSw63ccYAhL63y+Hw380b+HTuTVNMh7KYYXjRfz78KrsSGxe31zrgxY1h7AEbMXg+4Dyi+bhdXnx2p6bz9/FbYQWFHJJL5u4OwzW7QferScjMwCRRxX8w8bo75DACjQUtM176ILa09q3YcYOSc9Rw+Xuz2ns5Zm+X3XcbiTbl+DSjqoCw1mbef3xG/jCUiNVcgZjbcLRGda/qJp2Kfp6V5LwC/nDGYZoOfhLjaWACz2eQx6EDZreC92qRVWHxclfF6qu0REQeFHRGJSIHsINyrTRozhrVn8kcbubbgbe6M+YBYk43fqM/P3R+nY9b1Za73dSu4AopIcCnsiEjEqaxBnwlHB+EemalVnjHplXKEng2nYDrhqLH5/fTeJA98lkZ1yocUbQUXCW8KOyIScbxt0Ff6dPDSKqxzsdth9UxYMhFTSSHEW6HPEzQ8dyCY3AcnbQUXCW8KOyIScarTQbjCOp/TTsD8kbDrf44nWlwOV08H62kV/hxtBRcJbwo7IhK2PM3AVHXZyHOdz3G+mvsvshLnElNyDGJrwxWToOMtHmdzTuWs8zk1SKVqK7hIyCnsiEhYqmgGpioN+jzV+aRxgKmxL3KRZSOUgNG0C6ZrZkDyGT6P2d87rUTEP8L6IFARiV42u8HKHQeYv2EPK3ccwGY/GUM8HcTp3Gm1eFMuE/q6Px0cHMtIvds4QofzvuXrfAyus3zJZ/F/5yLLRgqNWCYVD2PVRW9WKeg4OXda9Wt3Gl1bNFDQEQkDmtkRkaCrbNbGm51Wy++7zO2ykdkEdgNeWbGLV1bsct23qMTuuiaFg0yJfZnLLBsAWG9vyT3Ff+UnI522BZ775YhIZFLYEZGgqqw/zuisVl7vtCq9bLR4Uy6vrtiF/ZQbl74vGAy0fMn4mNkkmY5RZMTwVMl1vGi7CvsfE93aHi4SfRR2RCRovOmP89qKXV7dq3SDvk4ZyYx9Z4Pb65z3XbJ6PXMTZ9LV+BaADfYW3FP8f2w3mgDaHi4SzRR2RCRovOmPU9GxC6WVnoGp+L4Ggy1LGXdiDkmm4xQZsTxZch0v23q7Du90bhcffEFTPv5+rwqLRaKMwo6IBI23/XHqJcaSd7zY651Wnu7b1LSPx2JeppvlBwB21zqHZa0nMj/bjC3/5GustWIBeOrzba7HdHq4SPTQbiwRCRpv62Fu7pYBlN9p5alB36n3NWPnZsunfBZ3P90sP3DciOPh4hu4+OA4Hvz6BGAwJqsVzwxux5isM8k7VszhY2VnlJy1Pouyc3z5FUUkDCnsiNRQFW39DhTnsQqeFodMOGZURl3WkhnD2pNqLRtiUq0Jbg/4LH3flqZfeTduIhNi36SWqYiVtkx6npjKq7YrXUXI+/KLePrzbcSaTcxb+4vHGiJw7PwKxnsjIoFjMgyjxv9bnJ+fj9VqJS8vj6SkpFAPRyTgKjwyIcDLNs7dWOD+WIXSYabCM6xO8dn3P7Pp7YmMjPmQOJONI0YiU0r+zFzbpRhu/l5nAurXjuWgF1vN597WRaeUi4Qhbz+/FXZQ2JGaxdPWb3dhI5Bj8GvY+vUbmD8KftsMwOe283mg+BZy8U9AeWZwO/q1q/h8LBEJPm8/v1WgLFKDeLP1e+KCTfTITA3oTiS/HatQdASWTII1LwIG1GqIvddUate6hCt+yOWNlT/7ZbzqvSMS2RR2RGoQb7Z+Oxv2BXrZxnmsgi9KL2udmbeC1usmYsr/1fFk20HQcwrm2g3oCmAyeRV2kmvHcajghNc7v0Qk8ijsiNQg3m799va6YHIufZXk5TIhdhZnW1YDcKx2E2pdMw1aXl7memfRcmWHhY7vk8nIOetdvXZKPw/ld36JSOTRbiyRGsTb5ZhwW7ZZlJ3DHbO/4dKjC1kSfw9XWVZTYph5oeQqOhyYxKLCzHKvsZhNHg8LLR1kerdN82nnl4hEHs3siNQg3s52hNOyjc1u8Ob8T3k37nk6mB1N/763ZzCu+DZ+MJpXWGfUq02a28NCU08phvZbDZGIhCWFHZEaxDnbMWJ2hCzbnDhGzkcP8/qJF4k12zhqJPBEyfW8YbvC1TOnsjojb4NMVWqIRCQyKOyI1DDeznaE3Nb/wid/o8nhX8AEn9k68lDxcHI8bCevqM5IQUakZlPYEYlClTXjC+tlm7w9sOh+2PwRAEW10rjz8GD+a7+gwpeFW52RiIQPhR2RKONtw76wm+2wFcPqF+CLKXDiKJgs0PUOYi66j41PrcEUQXVGIhJetBtLJIo4uyOf2ksn1IdaVnoO188r4YWL4b//dASdpp3h/76CKyZjSajr1a6qsJiVEpGwpJkdkSgRLt2RT1XhTFPzGFj8IHw3x/FEYjJkPQTn3wDmk38Xi5g6IxEJSwo7IlEinLojO3k6h2t/3jFWzJ3K5bXeJ7Y4HwB7++Gsa3kXe08k0njnociqMxKRsKawIxIlwq07sqeZpvamrUyKfY1zzD9DMRipbVl19j8Y+3UcOV9vcV0XEXVGIhIRVLMjEiXCrTvyqTNNDcnjidiZ/Cf+Ic4x/0yeUYvxxTcxvcVL/PlTe9jVGYlI9NDMjkiUCLfuyM4ZpBhKuNGymNEx75FkOg7A2yWXMLVkMAdJot7Xv4RdnZGIRBfN7IhECW/PggpWaGhcN4ELzdl8EjeOB2PfJMl0nI325lxTNJH7Sm7nIEkAHD5e7PEepeuMRESqSjM7IlEkbHYtHfqZLmv/Sde4BQAcMOryr5JBvGO7xHXMgwmwJsZWGHacwvEUdhGJHAo7IlEmpLuWThTA/56Er6djshVhN1mYVZzF0yUDyKOO6zLnSG7u1pynPt9W6W3VHVlEqkNhRyQKBX3XkmHAxndh8QQ4stfxWMZFmHs9RtpvydRasIk8NzNNPTJTmbd2d9jUGYlIdFLYEQmCys6qCtW9/OLXb2DROPh1jeP7eqdDz0eg9VVgMtErhQpnmiLqFHYRiUgmwzDc/YWqRsnPz8dqtZKXl0dSUlKohyNRxtuzqoJ9r2rL+xU+f8gxowMQWxsu+ht0GQmxvi07hdXvJSIRw9vPb4UdFHYkcDx1EHbOU8wY1r7CD/PSszi7fj/G059vrfK9/OZEAax4BlZMg5LjGJj4rcW1fNvyTpIaN6vyTFPYzViJSNhT2PGBwo4Egs1u0H3qUo9HODjrUZbfd5nbD3V3sx2eVHYvv7DbYMMcWDoZjuYCcLBhR8bmDeKLI6e5LtOMjIgEi7ef3+qzIxIgvpxVdSpPp5dX5V5+sWOZ41Tyj0Y5gk690/m2yzN0+HVMmaAD6nwsIuFHYUckQKp6VlVFp5f762d6f8PN8Nb18GZ/2LcR4q1wxWRsd6zhjm+bYpRrX3iyyHjigk3Y7DV+4lhEwoB2Y4kESFXPqqpsRsgfP7NSR3Jh2SPw7Www7GCOgQv+AhffB7WSWbPjQNidsC4i4onCjkiAVPWsqqrMzvitH03RUfh6Gnw9HYqPOR47uy9c/hA0bOnzGFds/02FxiISclrGEgmQqp5V5evsjF/60ZScgDUvwbR28OVUR9Bp0glu+QwGzS4TdHwZ47PLdtB96lLV74hISIV12JkyZQoXXHABdevWpXHjxvTv358tW7aUuaawsJCRI0fSoEED6tSpw4ABA9i3b1+IRixSlvOsqlRr2XCQak3wuFXcOSPkbWyp6F6Vstsh+314rhN8cg8U/AbJZ8DAWXDrf6FZF7cv82WMKlgWkVAL663nvXr1YvDgwVxwwQWUlJTwj3/8g+zsbDZt2kTt2rUBGDFiBAsXLuT111/HarUyatQozGYzK1as8PrnaOu5BJqvPWScu7GgfFdhAxiT1YrmDWtXrx/NjmWOpoA5Gxzf124Ml9wH7YeDJbbSl3saoztB2RovIjVOVPbZ+e2332jcuDFffvklF110EXl5eTRq1Ig5c+Zw3XXXAfDjjz9y9tlns3LlSrp0cf+30lMp7Eg4ClhX4T3r4POJsPNLx/dxdeDCu6DrSIivU/FrvRhjRebe1kUFyyLiN95+fkdUgXJeXh4AycmOIsx169ZRXFxMVlaW65rWrVvTrFkzn8KOSDjy++nlv22FpZNg80eO782xcMGt8Kd7oE6jao3xqcVbeXbZ9kqv9/vWeBERL0RM2LHb7YwePZpu3brRpk0bAHJzc4mLi6NevXplrk1JSSE3N9fjvYqKiigqKnJ9n5+fH5Ax1xRq8+8fnt7Has+EHNoFX0yF7+c5tpFjgvOGwCX3Q/3TfRpTh9Prs+7nQ+XG2K1lQ6/Cjt+2xouI+CBiws7IkSPJzs5m+fLl1b7XlClTmDhxoh9GJTrA0T8C8j7m58D/noB1s8Be7HjsrD5w2QOQklmlMZlNULpPoHOMPTJTq7TNXkQkGMJ6N5bTqFGj+Pjjj1m2bBlNmjRxPZ6amsqJEyc4fPhwmev37dtHamqqx/uNGzeOvLw819fu3bsDNfSo5ulIA+2+8Y3f38eC3+Gzfzq2ka992RF0zrgU/rIUhszxOui4G9OpDZGdY1y8KbdK2+xFRIIhrMOOYRiMGjWKDz74gKVLl5KRkVHm+Q4dOhAbG8uSJUtcj23ZsoVffvmFrl27erxvfHw8SUlJZb7ENxUdaaDjArzn1/fx2EFH4fHTbWHls1BSCE07w/CP4cYPoUmHao+pojH2yEz1eZu9iEgwhPUy1siRI5kzZw7z58+nbt26rjocq9VKYmIiVquVW2+9lbFjx5KcnExSUhJ33nknXbt2VXFygPlyyKV233jml/fx+GFYNQNWPQ9Ff9SfpbVzLFe1zAKTb7Mpvh5XUXqMfi+qFhHxg7AOOzNmzADgkksuKfP4a6+9xk033QTAU089hdlsZsCAARQVFdGzZ0+ef/75II+05qnqIZdSVrXex+OHYfVMjJXPYfoj5BTUb03CFeOxtO7jc8jxdUyeXueXomoRET8K67DjTQughIQEnnvuOZ577rkgjEicqnrIpZRVpffx+GFY/QKseg4K8zABW+xNeKbkWj7N6UTqh4lM6Jtb5WWjqv5/pv+vRSRchXXYkfBV1UMupSyf3sdjB2H1TFg1E4ocPae22k/jmZIBfGLvhPFHCZ6zaLiqdTKVjanCMYqIhKGwLlCW8FXVQy4DwWY3WLnjAPM37GHljgMRVRTtzfs4+Yo0LMsmOQqPv5wKRXkYjVrzQMxYep6YykJ7F1fQgeoXiFc0plNpp5WIRIKIOi4iUHRcRNWFus9OqH++J1U5C+vU3+OcpOM823wFGTvfhuICx4MpbeCie1kZ340hL6+pdBzVOZ7Blz472mklIqEQlWdjBYrCTvW4+2AHAr4jx9kL5tR/gJ0/JVTbnasawJzv45F9P9H251mkbH8Hk+2PTt+pbeHi++Cs3mA2M3/DHu6et6HSsTwzuB392p1W5d/F2w7KIiKhEJVnY0l4OnX3TTBmWyrrT2PiZO+XYH4Yewpg3tTRWA5up+vGpx3HOthLHA826QQX3QutepTZXRWsAnF3O6u000pEIo1qdsSvgtVV2Zf+NMFS5QaBezfAOzfCsxfAhtmOoJNxEQxfALf+F868otw2cmcRsacYZ8IRMFU0LCKisCN+FMyuyuHY58enAGYY8NOX8OY18OLFsGm+44qzesOtnzuCTsZFHnvlhFOBuIhIuNMylvhNMLsqh2OfH2+ClRk7lh/nw5I3YO+3jgdNFmgzALqP8ercKqdebdKYMax9uSXDVBUNi4iUobAjfhPM2ZZw7PNTUbCK5wQDLP/jNsvHZKzd53gwJgHOvwG6joTkDI+vrYg/j2fwdQeZiEikUNgRvwnmbItzGWfE7PWYoEzgCdUyjrsAlkw+N1gWc2PMf2lgOgKAkVgfU6fbodPtULthtX+uP45nCNct/CIi/qCaHfGbYBfNOpdxwuWU7dJ1NGeYcpgc8wpfx9/JmNj3aWA6wm57Izaf9w9Mo7Ph0n/4Jej4Q7CKykVEQkUzO+I3oZhtCatTtg2DXrW3sTrjRRrnfOF6+Ht7Bu/EXcOfrr6Fnm2bBn9cFQjXLfwiIv6ksCN+FYqi2ZCfsl1SBD98ACufhdyNNP7j4YNNLmfT6TdgyfgTE89oUC4shEONTDCLykVEQkVhR/wurGZbAunofvjmVVj7ChTsdzwWkwjnD4XOI0hu2JLuHl4a7BoZT8EqHLfwi4j4m8KOBETIZ1v8xG1IyP0OVr8A2e+B7YTjwrppcMFfoOMtUKvimqTqdFmuioqCVThu4RcR8TeFHREPSoeEGEq40ryG2vGLaWtsOXnRaR2hywjI7AeW2ErvGewamcqC1XN/Pj/stvCLiPibwo6EtVDVtThDQkMOcZdlGUNjPifFdBgMKDYs/NasF+lXjIGmF/h032DWyHgTrCYt3Mz4PpmMnBM+W/hFRPxNYUfClqfll/F9zqZ+7fiABSCbzc5H899leuxCeprXEmuyAbDfqMcc22XMKbkcy/40lp/WEYuP9w5mjYy3wap+7Th1YhaRqKawI2HJ0/JLTl4hd8z5tsxjfivsLcyD79+haMWLPF+8FWeSWWs/kzdLruBTeyeKnf/KlJp98WX2KZg1Mr4Eq37tTqsZReUiUiMp7EjYqWj5xZ1qF/bu/daxq2rje1B8jFrAMSOeD23deNPWg83G6W5f9ml2Dmt2HmTuml/IzfduV1Uwj7nwNVhFS1G5iMipFHYk7FS2/HKqKhX2Fh2B7Pdh3esnD+QEaNSanadfT7/lTcmndoW3eGPlz24fryh8BbPxYjieHyYiEgo6LkLCTlXqVUoX9nq+yIA962HB3fDv1o7/3fstWOLg3IFw86dwxyqa9R5LbWsDj8deeDMWcIQvm718zAjWMRelj6849XdR8bGI1CSa2QmQQO4iCofOu4FUnXoVt0Hp+CH4/l1Y/wbs23jy8QYtocNNcN6QMudUWUx4nH3xVmW7qoLVeDEUHa1FRMKNwk4ABLI7bk04nbqy5ZeKuIKS3Q67voL1b8LmBWArcjxuiYfMqx0h5/RuYHIfLjyFBF9VNEsVrBqZGtPRWkTEA5NhGFX5i2tUyc/Px2q1kpeXR1JSUrXu5WkXkfNjxddlitKzOLt+P8bTn2/1273DmfN9BO9mVpz1J8tvPwPL93Nhw1zI++XkBSnnQvsbHMtVlXQ4Ls35/n+aneOxRqcic2/roqJfEZEA8fbzWzM7fuTv7rjuZnHcCeTp1KFaMvNlZqUOx+llWcO4pG+xTF9z8ol4K5w7ANrfCGntPM7iVKT07IsvYUfFvyIi4UNhx4/82R3X0wyRP+7trVAvmblbfjlUcIJJCzexL+8YF5p/4FrL/7jSspZEiuA3ABOccQmcPwxa94HYRI/39yXI+bK0puJfEZHworDjR/7qjutrnxlv7+3Lh3uwD6v0pFxdS242vdp/QfF37xB/LPfk48ktHIXG5w2Gek0rva+vQa6iLeOnUvGviEh4UdjxI391x/W1z4w39/blwz3Yh1VWKu9XR8O/79+B/T9gBuIBEupBmwGOkNOko9fLVFUNch53NiXFM6RTM5o3rK3iXxGRMKSw40f+auJWlT4zFd3b1w/3YB5W6VHBAdj0oSPk/PL1ycctcXBmT2g7CFpdATHxXt3OOauVm3ecSQs3VznIaWeTiEjkUdjxI391x/W1z0xF967KLI23YevT7BwA/33YF+bDlk8cnY13LAV7ycnnml0I5w2CzH6QWN+n23pb6A0ng9xTi7fSrWVDt7+bjlUQEYksCjt+5o8mbr72mano3lWZpfE2bL2x8mfeWPlz9YqWi47C1kXwwwewbfHJfjgAqW0dW8XbXAvWJr7fG98LvZ2eXbadZ5dtj7oeRiIiNZHCTgBUd6mjshkiAxiT1cqrGpGqFE37GrZ8Llp2BpxNH8K2z6Hk+MnnGp4J51wL514HDVt5NXZPqlPo7RTsgmwREfE/hZ0Aqe5Sh7/a/FelaNqXnUfgZdHy8cOw7b+waT5s/xxKSoWw+hmO2ZtzroWUc6rUD8ed6hR6O4WkIFtERPxKYSeM+aMYtqpF074el+C2aPnofkcNzuYF8NOXYC8++YLkMyCzP5zT37Fc5aeAU1pVCr3dCUpBtoiIBIzCTpir7gxRdYqmS4ctb49LKMj5EXLXwY8LYffqsj+xUWs4u6+jyDilTUACTmnVOVDUHX+FJxERCS6FnRqgOktilR2XYMZOO9N2LresJ8u8nrM+/7XsBWntHAdvtu6LrUErxyxVbiGNCw4GfMu2N7VHybVjGdKpGc8t21Hp/fwdnkREJDgUdqJEZd2Rq7skVjo41OI43c0byTKv51LLBhqa8l3XGeYYTM27Q+ur4KwrXbuoFmXnMPHlpUE9esKbWa1HrzmXHpmp/Gf9nmr3RxIRkfCkU8/x76nnoRCUM6x+387mr97lwLcL6GT+kTiTzfVUvlGLL+1tadZ1AOddOrBcHxx/nwTvK2/eH0+nrEfjifIiItHC289vhR3CK+z4esp4wILEiQLYtdyxc2rbYji0s8zTO+0pLLW353N7e3bXOY8Hrm7r9ufY7Abdpy71WOTsnDVZft9lAV3S8uZ9DfXBpyIi4htvP7+1jBVGfP2w9esZVoYB+35wdC7esRR+/rpsgz9zLDTvBq2uwNbyCnLzk2l4pJC7KglkYXH0BN4VeusoCBGR6KSwEyaqcjhltYNE/l7HlvCfvoCflsHRfWWftzaDVlnQsgdk/Ani6wJgAbo28u738tdJ8MGioyBERKKPwk4YqOoMjc9B4vghx4zNT184Qs7vW8peGFsLmneHFpdBi8sdHYyruT3cXyfBi4iIVJXCThio6gxNZQGhDse4wLyFC7YuhdVrIed7ypXfpreDMy5xBJymnb0+Rdxb/joJXkREpKoUdsJAVZd6Tg0SVo5ygXkLnc2b6WzezDmmXVhMBmwq9aIGLSHjYkfAad4dagU2ZPja1NDXAm0REZHKKOyEgaou9Vjyf2XmeTv4/uvPuMC8hdbm3eVec6x2U2qdeQlkXOQIN0np/hiyT7xtaqjdUCIiEgjaek7ot547t2dXtNRzWlIMX97YCMuvaxzHMOxeA/m/lrt2hz2N1faz2ZLQlkuu6M+lnc4P+Pi9VdGsTah78YiISOSpcX12nnvuOf71r3+Rm5vLeeedx/Tp0+nUqZNXrw112IHyTe1SOcD55u2cb95OO/MOOsTuwmI7ZbnLZIG0ttDsQmxNO/OtcRZ7SupG3PJPuPTiERGRyFKj+uy8/fbbjB07lpkzZ9K5c2eefvppevbsyZYtW2jcuHGoh1e5YwfplfADn3VYTc7mlZxl20aq6VDZa2xAQj1o2umPr86Q3h7i6wCO7eAd//iKNOHSi0dERKJTVISdJ598kttuu42bb74ZgJkzZ7Jw4UJeffVV7r///hCP7hQFByDnW8j5DvZucPzvYccBm2f+8YUJ7CYLx+q3plZGZ8xNOkKTjtCgFZjNIRy8fzmXtT7NzvHq+nDpxSMiIpEl4sPOiRMnWLduHePGjXM9ZjabycrKYuXKlW5fU1RURFHRye7A+fn5bq+rFsOAgz9B7sayX0f2ur8++QzHTM1p7SG9Pea0ttSJq+3/cYUJd8XIlVEvHhERqYqIDzu///47NpuNlJSUMo+npKTw448/un3NlClTmDhxYuAH93IWHD9Y/vHkFo7+NmnnQVo7R93NKYdnRjNPxcieqBePiIhUR8SHnaoYN24cY8eOdX2fn59P06ZN/ftDTCZHbc3RfZB6LqS2dfxvyjmuYxdqooq6RbvjrhePiIiILyI+7DRs2BCLxcK+fWXPddq3bx+pqaluXxMfH098vH87Bbs1ZF61j1uINpUVI5/q1F48IiIivor4ate4uDg6dOjAkiVLXI/Z7XaWLFlC165dQzgyFHTc8LbI+MaupzP3ti4sv+8yBR0REamWiJ/ZARg7dizDhw+nY8eOdOrUiaeffpqCggLX7iwJH94WGV/ZJk3bzEVExC+iIuwMGjSI3377jQcffJDc3FzatWvHokWLyhUtS+jpYFAREQm2qOmgXB3h0EG5Jjm1W7STjoYQERFfePv5HfE1O5HEZjdYueMA8zfsYeWOA9jsNTNnOg8GTbWWXdJKtSb4LejovRYREaeoWMaKBOF4ondFB3MGWq82afTITA3Izw/H91pEREJHy1gEfhkrFCd6VxZkojUQ6PR0EZGao8adel4dgQw7oTjRu7IgE62BQKeni4jULKrZCRO+nOjtD84gc+rPzM0rZMTs9Xzy/V6PHYydj01csCkia1yC/V6LiEhkUNgJMG+b6PnjRO+KjmJwPvbA/OyoDQTBfK9FRCRyKOwEmLdN9Pxxorc3MxsHC4q9ulckBoJgvtciIhI5FHYCzNlEz1OFiAlHPY0/muj5M6BEYiAI5nstIiKRQ2EnwCxmExP6ZgKU+xD294ne3gaU5NpxURkIgvlei4hI5FDYCYJgNNED72c2Jvdr4/r+1OchsgNBsN5rERGJHNp6TvCOiwhGEz9vj2KI1j47TqFsmCgiIsGhPjs+CNXZWIH6QPY2yCgQiIhIJFPY8UEowk6gZ1YUZEREJNop7Pgg2GEnWjsYi4iIBJM6KIcpbxr/hbKDsU4LFxGRaKNTz4PMlyMNurZoELyBEb2Hg4qISM2mmZ0gC9cjDSo7U2tRdk5QxyMiIuIvCjtBFo5HGoT70pqIiEh1KOwEWTgeaaDTwkVEJJop7ARZOB5pEK5LayIiIv6gsBMC4XakQTgurYmIiPiLdmOFSK82afTITA2Lxn/OpbXcvEK3dTsmHEEsEg8HFRERUdgJIYvZFPTt5Z7GMaFvJiNmr8eE+zO1IvlwUBERqdm0jCVA+C2tiYiI+ItmdsQlnJbWRERE/EVhR8oIl6U1ERERf9EyloiIiEQ1hR0RERGJago7IiIiEtUUdkRERCSqKeyIiIhIVFPYERERkaimsCMiIiJRTWFHREREoprCjoiIiEQ1dVAGDMNx9GV+fn6IRyIiIiLecn5uOz/HPVHYAY4cOQJA06ZNQzwSERER8dWRI0ewWq0enzcZlcWhGsBut7N3717q1q2LyeS/Qy/z8/Np2rQpu3fvJikpyW/3lbL0PgeP3uvg0PscHHqfgyOQ77NhGBw5coT09HTMZs+VOZrZAcxmM02aNAnY/ZOSkvQvUhDofQ4evdfBofc5OPQ+B0eg3ueKZnScVKAsIiIiUU1hR0RERKKawk4AxcfHM2HCBOLj40M9lKim9zl49F4Hh97n4ND7HBzh8D6rQFlERESimmZ2REREJKop7IiIiEhUU9gRERGRqKawIyIiIlFNYSeAnnvuOZo3b05CQgKdO3dmzZo1oR5SVJkyZQoXXHABdevWpXHjxvTv358tW7aEelhR77HHHsNkMjF69OhQDyXq7Nmzh2HDhtGgQQMSExM599xz+eabb0I9rKhis9kYP348GRkZJCYm0qJFCyZNmlTp2UpSua+++oq+ffuSnp6OyWTiww8/LPO8YRg8+OCDpKWlkZiYSFZWFtu2bQvK2BR2AuTtt99m7NixTJgwgfXr13PeeefRs2dP9u/fH+qhRY0vv/ySkSNHsmrVKhYvXkxxcTFXXHEFBQUFoR5a1Fq7di0vvPACbdu2DfVQos6hQ4fo1q0bsbGxfPrpp2zatIl///vf1K9fP9RDiypTp05lxowZPPvss2zevJmpU6fy+OOPM3369FAPLeIVFBRw3nnn8dxzz7l9/vHHH2fatGnMnDmT1atXU7t2bXr27ElhYWHgB2dIQHTq1MkYOXKk63ubzWakp6cbU6ZMCeGootv+/fsNwPjyyy9DPZSodOTIEaNVq1bG4sWLjYsvvti4++67Qz2kqHLfffcZ3bt3D/Uwol6fPn2MW265pcxj1157rTF06NAQjSg6AcYHH3zg+t5utxupqanGv/71L9djhw8fNuLj4425c+cGfDya2QmAEydOsG7dOrKyslyPmc1msrKyWLlyZQhHFt3y8vIASE5ODvFIotPIkSPp06dPmX+uxX8++ugjOnbsyMCBA2ncuDHnn38+L730UqiHFXUuvPBClixZwtatWwH47rvvWL58OVdeeWWIRxbddu7cSW5ubpn/flitVjp37hyUz0UdBBoAv//+OzabjZSUlDKPp6Sk8OOPP4ZoVNHNbrczevRounXrRps2bUI9nKgzb9481q9fz9q1a0M9lKj1008/MWPGDMaOHcs//vEP1q5dy1133UVcXBzDhw8P9fCixv33309+fj6tW7fGYrFgs9l45JFHGDp0aKiHFtVyc3MB3H4uOp8LJIUdiQojR44kOzub5cuXh3ooUWf37t3cfffdLF68mISEhFAPJ2rZ7XY6duzIo48+CsD5559PdnY2M2fOVNjxo3feeYe33nqLOXPmcM4557BhwwZGjx5Nenq63ucopmWsAGjYsCEWi4V9+/aVeXzfvn2kpqaGaFTRa9SoUXz88ccsW7aMJk2ahHo4UWfdunXs37+f9u3bExMTQ0xMDF9++SXTpk0jJiYGm80W6iFGhbS0NDIzM8s8dvbZZ/PLL7+EaETR6d577+X+++9n8ODBnHvuudxwww2MGTOGKVOmhHpoUc352Reqz0WFnQCIi4ujQ4cOLFmyxPWY3W5nyZIldO3aNYQjiy6GYTBq1Cg++OADli5dSkZGRqiHFJUuv/xyNm7cyIYNG1xfHTt2ZOjQoWzYsAGLxRLqIUaFbt26lWudsHXrVk4//fQQjSg6HTt2DLO57EefxWLBbreHaEQ1Q0ZGBqmpqWU+F/Pz81m9enVQPhe1jBUgY8eOZfjw4XTs2JFOnTrx9NNPU1BQwM033xzqoUWNkSNHMmfOHObPn0/dunVd675Wq5XExMQQjy561K1bt1wdVO3atWnQoIHqo/xozJgxXHjhhTz66KNcf/31rFmzhhdffJEXX3wx1EOLKn379uWRRx6hWbNmnHPOOXz77bc8+eST3HLLLaEeWsQ7evQo27dvd32/c+dONmzYQHJyMs2aNWP06NFMnjyZVq1akZGRwfjx40lPT6d///6BH1zA93vVYNOnTzeaNWtmxMXFGZ06dTJWrVoV6iFFFcDt12uvvRbqoUU9bT0PjAULFhht2rQx4uPjjdatWxsvvvhiqIcUdfLz8427777baNasmZGQkGCcccYZxj//+U+jqKgo1EOLeMuWLXP73+Thw4cbhuHYfj5+/HgjJSXFiI+PNy6//HJjy5YtQRmbyTDUNlJERESil2p2REREJKop7IiIiEhUU9gRERGRqKawIyIiIlFNYUdERESimsKOiIiIRDWFHREREYlqCjsiIiIS1RR2RCRqXXLJJYwePTrUwxCREFPYERERkaim4yJEJCrddNNNzJo1q8xjO3fupHnz5qEZkIiEjMKOiESlvLw8rrzyStq0acPDDz8MQKNGjbBYLCEemYgEW0yoByAiEghWq5W4uDhq1apFampqqIcjIiGkmh0RERGJago7IiIiEtUUdkQkasXFxWGz2UI9DBEJMYUdEYlazZs3Z/Xq1ezatYvff/8du90e6iGJSAgo7IhI1LrnnnuwWCxkZmbSqFEjfvnll1APSURCQFvPRUREJKppZkdERESimsKOiIiIRDWFHREREYlqCjsiIiIS1RR2REREJKop7IiIiEhUU9gRERGRqKawIyIiIlFNYUdERESimsKOiIiIRDWFHREREYlqCjsiIiIS1f4fATw1KXzpHGMAAAAASUVORK5CYII=",
+      "text/plain": [
+       "<Figure size 640x480 with 1 Axes>"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    }
+   ],
+   "source": [
+    "A_new = np.column_stack((one_vector, t, t**2))\n",
+    "x_hat_new = gen_xhat(A_new, y)\n",
+    "y_hat_new = A_new @ x_hat_new\n",
+    "plt.plot(t, y,'o')\n",
+    "plt.plot(t, y_hat_new)\n",
+    "plt.title('t vs y')\n",
+    "plt.xlabel('t')\n",
+    "plt.ylabel('y')"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 73,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "e_hat = y-y_hat_new"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**End of notebook.**\n",
+    "<h2 style=\"height: 60px\">\n",
+    "</h2>\n",
+    "<h3 style=\"position: absolute; display: flex; flex-grow: 0; flex-shrink: 0; flex-direction: row-reverse; bottom: 60px; right: 50px; margin: 0; border: 0\">\n",
+    "    <style>\n",
+    "        .markdown {width:100%; position: relative}\n",
+    "        article { position: relative }\n",
+    "    </style>\n",
+    "    <a rel=\"license\" href=\"http://creativecommons.org/licenses/by/4.0/\">\n",
+    "      <img alt=\"Creative Commons License\" style=\"border-width:; width:88px; height:auto; padding-top:10px\" src=\"https://i.creativecommons.org/l/by/4.0/88x31.png\" />\n",
+    "    </a>\n",
+    "    <a rel=\"TU Delft\" href=\"https://www.tudelft.nl/en/ceg\">\n",
+    "      <img alt=\"TU Delft\" style=\"border-width:0; width:100px; height:auto; padding-bottom:0px\" src=\"https://gitlab.tudelft.nl/mude/public/-/raw/main/tu-logo/TU_P1_full-color.png\" />\n",
+    "    </a>\n",
+    "    <a rel=\"MUDE\" href=\"http://mude.citg.tudelft.nl/\">\n",
+    "      <img alt=\"MUDE\" style=\"border-width:0; width:100px; height:auto; padding-bottom:0px\" src=\"https://gitlab.tudelft.nl/mude/public/-/raw/main/mude-logo/MUDE_Logo-small.png\" />\n",
+    "    </a>\n",
+    "    \n",
+    "</h3>\n",
+    "<span style=\"font-size: 75%\">\n",
+    "&copy; Copyright 2024 <a rel=\"MUDE\" href=\"http://mude.citg.tudelft.nl/\">MUDE</a> TU Delft. This work is licensed under a <a rel=\"license\" href=\"http://creativecommons.org/licenses/by/4.0/\">CC BY 4.0 License</a>."
+   ]
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "TAMude",
+   "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.12.5"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 2
+}
diff --git a/src/students/Week_1_3/data.csv b/content/tutorials/Week_1_3/data.csv
similarity index 100%
rename from src/students/Week_1_3/data.csv
rename to content/tutorials/Week_1_3/data.csv
-- 
GitLab