{ "cells": [ { "cell_type": "code", "execution_count": 1, "id": "2942bc59", "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import math as mth\n", "import matplotlib.pyplot as plt\n", "import sklearn.tree as tr\n", "import sklearn.model_selection as ms\n", "import sklearn.metrics as met\n", "import scipy.stats as stt" ] }, { "cell_type": "code", "execution_count": 2, "id": "2915f23f", "metadata": {}, "outputs": [], "source": [ "#Creating Artificial Data:\n", "np.random.seed(25)\n", "Min = 0\n", "Max = 2 * mth.pi\n", "N = 2500" ] }, { "cell_type": "code", "execution_count": 8, "id": "dbbad9b0", "metadata": {}, "outputs": [], "source": [ "t = np.random.uniform(Min, Max, N)\n", "y = np.sin(t) + np.random.rand(N)/35" ] }, { "cell_type": "code", "execution_count": 9, "id": "4f768726", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAfIAAAD4CAYAAAAXfWQCAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAA3GUlEQVR4nO3deXiU9b3//+c7yyRhwpJg2CKbEUFEwDSC4oZVzhGrYm2xYK1YtUhdip7WHvXbq9ae/qy/+m2VaiviBrQqlVYLrVuVIxZElhABkUUJENmEyIQlA8lk+Xz/mMmQQFgTcs9kXo/rypW5t8w7UfLK574/iznnEBERkfiU5HUBIiIicuIU5CIiInFMQS4iIhLHFOQiIiJxTEEuIiISx1K8LuBEnHLKKa5Xr15elyEiItIili5d+pVzLqexY3EZ5L169aKwsNDrMkRERFqEmZUc7phurYuIiMQxBbmIiEgcU5CLiIjEMQW5iIhIHFOQi4iIxDEFuYiISBxTkIuIiMQxBbnEjEAwxDMfFBMIhigqKeOy385l7pod0X0iInIoa471yM3sBeAqYIdzbkAjxw2YBFwJ7ANuds4VRY5dETmWDDznnHv0aO9XUFDgNCFMfCouLedX/1zF3V/vwwef7WDbrgrmflZKsLKKtNRkAsEqstukUravCgekJhtVNY78Hh342Tf68/+/vYb9oRrAkeFL4ZHrziYvJ9Prb0tE5KQys6XOuYJGjzVTkF8MlAPTDxPkVwJ3Ew7yocAk59xQM0sGPgNGAJuBJcBY59yqI72fgjz+BIIhpi3YwGtFW9hUtj8a0McjyaDWHbpv+BmnsGt/NanJSQp2EWmVjhTkzTJFq3Pu32bW6winjCIc8g5YaGYdzKwr0AtY55xbHyl0RuTcIwa5xL5AMMTkucWs2rab/xrRl9/+ay0fFu+MHj/eEIdDQ7xu3/+u/Sq6PX56ITMnDANgZuEmRhd0J9vvO/5vQEQkTrTUXOu5wKZ625sj+xrbP7SxL2Bm44HxAD169Dg5VUqTBIIhHp69kndWbWfYaR3537WlACxav4Cq2uZ5DwMO9yeAAcWlQS58dA6nd/KzYstefvfuZ4wa1I37rzxTgS4irVJLBbk1ss8dYf+hO52bAkyB8K315itNmkNxaTk3Pb+ILbsqAKIhDhwS4klALeBLhozUFPZXVfPT/+zHh8U7ufvrfXh75TYKS8oIlFfyRWA/2f4UHrp6AK99vIXrzsnl0bfXsLO8gopqSEmCDhkpJCcls31vJQD7qmpZsWUvAJXVtby6dDNrvtxDu4xUHh41QLfeRaRVaakg3wx0r7d9KrAV8B1mv8SJQDDEE++u5S+Fm6msPnKz24DRBady/8jGW8e3XZwHQH7PrEavv3pwbvRzIBhqcOv88Xc/Y9Kczw/bYl+xZQ8AVzzxb8ae24N7RpyhFrqItAotFeSzgbsiz8CHArudc9vMrBToY2a9gS3AGOCGFqpJTlDd8+/CkgBf7q5g6+6KRs+74qzO9O3SlrJgFcWl5fzy2uZrDWf7fdx+SV50e9ywXrTxJXN5/87c/MJiNpXtByAt2ais9zy+qsYxfWEJf1pYwiPfHMDYoT2bpR4REa80V6/1V4DhwCnAduAhIBXAOTc5MvzsKeAKwsPPvu+cK4xceyXwBOHhZy845/6/o72feq17JxAMMeFPhSzeWHbE84b2zubpG7/mSau3uLSch2atpH/X9gzL68it05ZwuL51WW1Sefz6wQzv16llixQROQ4nffhZS1OQe6OopIzbpi8hEKyK7ksC+nXJZNWX5QztncWgUzuQ4Uth3LBeMXPruqikjHEvLGJvZc1hz5l687kKcxGJWQpyOWGBYIhH31rNW59so7yyptHnzxeefgoX9Tklpod6hcexb2TjV3v5x4ovDxnKlpwEM28fxpKNgZj+PkQkMSnI5YQEgiF++OelLNoQOORYZloyAKed4ufxMefEVU/w4tJy7pu5nBWbd1G/f57fl0wwVEPntmm8PP68uPqeRKR1U5DLcSkuLefB1z5hc9m+6HCy+rpnZTD1liFxH3R1gb52+x6ClQ173KclG99R73YRiREKcjlmxaXlXDXpA/ZXN368e1YGs+66sFWFWyAY4uu/fZ9d+w79pvN7dOCiPjkx9cxfRBLPkYJcq59JVFFJGVdOmndIiLfPSGHUoK5ckNeRqbcMaXWBlu338cK4IWS1ST3kWNEXu5g053Mmzy32oDIRkaNrqXHkEuPmrtnB96cuOaQzW9u0ZN7/yaWtLrwPlt8zizk/Hs6PXvmY+eu+OuT4ix+u54oBXQ47WY2IiFfUIk9wgWCIR95Yza3TCxuEeHabFHI7pDPtlqGtPsTrZPt9/H7sOUy8rA95p7RpcKyqFr773EKKS8s9qk5EpHF6Rp7AAsEQd79c1GBVsjoPjOzXYOa0RFM3XO3lhSWUBkPR/V3apfGdc7szbljvhPkDR0S8p2fkcoi6oWUHh3j79BTGX3waowu6H+bKxJDt93HviDOYMeF8LsjrSNvIcLsv91Qyac46fvRKEYF6AS8i4hW1yBNQIBjix68u4/16K5QBpCYbb99zcdwPKzsZikrK+NbTCw7pQ3BhXkd+f0O+WuciclKpRS4NTP6g+JAQ75CRyl/Gn68QP4z8nln8fszgQ9bdnV+8k8lz13lSk4gIKMgTSrhj2ypeXlTSYH/3rAz+9yfD1SP7KK4enMvffjiMUztkkFov0V9Z/IU6wYmIZxTkCaLudvqUeRsor7d4SOe2aa1ybPjJkt8zi++d35MqB0mRMN9bWcNNzy/i8XfX6rm5iLQ4jSNPAPVX/zrFn8pXwSq6tU/nqoFdmTD8dIX4carrCFhcWs6rhZsB2LKrgklz1rFtdwW/+fYgL8sTkQSjIG/lAsEQt01fEl3Cs6K6lgdG9tMKX02Q7fdx+yV5BIIhOmT4mP7RRioiq6/8delmbr8kT30NRKTFNEuQm9kVwCQgGXjOOffoQcfvA75b7z3PBHKccwEz2wjsBWqA6sP1ypPjEwiGmFm4iZ3lIQLBKlKTjaoax5Vnd03o8eHNKdvv48FvnMkVA7pw3dMLAKh1cNlvPyAt2fi/owdx9eBcj6sUkdauyUFuZsnAH4ARwGZgiZnNds6tqjvHOfcY8Fjk/KuBe51z9dfGvNQ5d+i8mHLCZhZu4tdvrYluf3NwLnmdMhN+fPjJkN8zi4G57VixZU90X2WNY+JflnFBnxzd+RCRk6o5OrsNAdY559Y750LADGDUEc4fC7zSDO8rR9C3c9sG21t3V3D7JXkKlZPk8THncE739g321TqYtmCjNwWJSMJojiDPBTbV294c2XcIM2sDXAH8rd5uB/zLzJaa2fjDvYmZjTezQjMrLC0tPdxpCS0QDPHIm6u5fvIC7nx5aXR/386ZPDzqLA8ra/3ycjJ5/c4LGZjbrsH+38/5nAl/KlRvdhE5aZojyA+eIwM4ZAKsOlcDHx50W/0C51w+MBK408wubuxC59wU51yBc64gJyenaRW3UjMLNzHl3+tZvLGMYCjc+Sq3Qzrv3HuJOl+1kMfHnMOlfXMYeVZnIPwP4e1Pt/PoW6u9LUxEWq3mCPLNQP0Hr6cCWw9z7hgOuq3unNsa+bwDeJ3wrXo5ToFgiJ3llSQnHfi7KsngybH5HlaVePJyMnnx+0NY/eXeBvtnL9uqVrmInBTNEeRLgD5m1tvMfITDevbBJ5lZe+ASYFa9fX4za1v3GvgPYGUz1JRQikrKuPx3c5kybwM1tY6UJMhMS+aFcedqtjaP/O76wWSkHvjnVVFdy80vaNIYEWl+Te617pyrNrO7gHcIDz97wTn3qZlNiByfHDn1m8C/nHPBepd3Bl43s7paXnbOvd3UmhLJP5Zt4UczluEIP+NwwKlZbZh736UeV5bY8ntm8c8fXcRDsz5lxeZd7KmoZsWWPdGe7feO6OtxhSLSWmj1szgWCIY491fvUhP5Tzi0VxZfBUM89u1BaonHkKKSMkZPXhD97zQwtx1TbxmqEQQicsy0+lkrNXlucTQcALbvrWTOj7X4SazJ75nFd4f2jG6v2LJHw9JEpNkoyONQIBji57NWHhIGv7t+sCf1yNHdM+IM8nt0iG7/dekmvvPMR1o1TUSaTEEeh5549zOmf1RCZU1tdN/A3PZqicewbL+P58ady8TL+tA9K4MtuypYtCHATc8vUuc3EWkSBXkcCQRDPPNBMe+v3RHd179LJhee3pHHxwz2rjA5Jtl+H/eOOINL+3aK7tuyq4JH39QYcxE5cVr9LI5MnruOKfM2kN0mhbSUJK44qzMPXTNAnabiTJY/tcH2zKLNjBnSQ3dUROSEqEUeRxauD0+IF9hXTWV1LXsqqhXicWjcsN4M7X0gtJ2Dsc9+xCNvrNJtdhE5bgryOFK6tzL6untWBj+7qr+H1ciJyvb7ePrGAsZfdBpt05IBqKx2TJm3QbfZReS4KcjjwNw1O8j/n3/Rv1t4RbP8Hh2YddeFmj89jtWtZf73uy6kS7u06P7XPt6inuwiclwU5DFu7pod3Dx1CYFgFYs3BnhgZD+eG3eubqm3Enk5mZze6cAfZNW1jm89/aHCXESOmYI8hgWCISb8+cAMdr07+rWmeCv08KgBXHj6KdElUHftq+bB11Z4XJWIxAsFeYwqLi1n1FPzqagOT93m9yXz+JhzPK5KToa8nEz+fNtQpt4ylNwO6QCU7NynBVZE5JgoyGPUT15dxqay/dHt7w7tqWfirVy238dlZ4bXMf9yTyWT5qxjZuEmj6sSkVinII9BgWCItfXWsx7aO4sJw/M8rEhaSlabhmPMu7VP96gSEYkXCvIYUjdz2xPvrmVfVXj61STg6RsL9Fw8QYwb1pv+XdtGt++esYyJrxTpFruIHJaCPIZM/qCYX7+1hukLv4jumzRmsEI8gWT7ffz5tvPo3PbAkLRZy7dx47MfKcxFpFEK8hhRVFLG9AUbGuwb2juLqwfnelSReCXb7+Pl8eeRmnxg36ovy7X0qYg0qlmC3MyuMLO1ZrbOzO5v5PhwM9ttZssiHz8/1msTQSAY4nvPL4z2UAfo0i6NR64b6GFV4qW8nEz+Mn5Yg31LSwJqlYvIIZoc5GaWDPwBGAn0B8aaWWNzh85zzg2OfPzyOK9t1SbPLSYYOrAkaa+ObXhz4sXqpZ7g8ntmMfXmc8lMSybVYP66nYx6ar4mixGRBpqjRT4EWOecW++cCwEzgFEtcG2rMHfNDp6bvz66nZoEz9+smdskbHi/TnT0p1EVuVmzqWw/N7+wWC1zEYlqjiDPBeoPdt0c2Xew881suZm9ZWZnHee1mNl4Mys0s8LS0tJmKNt7xaXl3DptCbUH7qgzdkgPtcSlgUv65jTY3lS2n4dnr/SoGhGJNc0R5NbIPnfQdhHQ0zk3CHgS+PtxXBve6dwU51yBc64gJyensVPiSiAY4uYXFlNT77vt0i6Ne0b09a4oiUn3XH4GEy87nZvO6xn9B/OPFdt0i11EgOYJ8s1A93rbpwJb65/gnNvjnCuPvH4TSDWzU47l2tbqiXc/azBzW7f26bz0g/N0S10Oke33ce+Ivvzy2gH851nhmd9qHYyfXqhb7CLSLEG+BOhjZr3NzAeMAWbXP8HMupiZRV4PibzvzmO5tjUqLi1nxpIDY8W7Z2Xwzx9dpFvqclR9u4Qni0lLSaK4NMgP/6wwF0l0TQ5y51w1cBfwDrAaeNU596mZTTCzCZHTvg2sNLPlwO+BMS6s0WubWlMsKy4t59qnPiQUuaeelgxTbxmilrgck3HDenNp3xwqq8OjHBZtKOPHry5TmIskMHOu0UfSMa2goMAVFhYe/cQYU1RSxthnP6IyMl48Jcl49fbzye+Z5XFlEk8CwRDnP/IelfU6WEy87HTuVf8KkVbLzJY65woaO6aZ3VrQxBkfR0M8CRTickKy/T5GndNwcMe23RUeVSMiXlOQt5Di0nJ276+Kbndpn64QlxN2/8gzaZ+REt2evWyLerGLJCgFeQsIBEPc8MxC9lRUk5IU7qH+1A35XpclcSzb7+PFm4eQ2yEdXxJUVDuu++OHCnORBKQgP8kCwRB3v1zE9vJKAKprYdywXmqNS5Pl98ziw/svY3CP8P9Lu/dXM3bKQnV8E0kwCvKTbPLcYj4s3hndPrtbW0YXdD/CFSLH59ffOrC4zo69lVz95Dy1zEUSiIL8JFu1bXf09dDeWUy7VZO+SPPKy8nEn3ZgzdMtuyr41T9XeViRiLQkBflJEgiGeOaDYv5rRF96dWwDwHmndVSIy0nxh7H5pEb+NWekJrFrX0itcpEEkXL0U+R4BYIhJvypkMUby2ifkcKk75zD2u17dUtdTprh/Tox9LRTmL/uK/ZX1fLxpt08+NoK/nL7sKNfLCJxTS3yk2Dy3GIWbywDwh2QHpr9KbdfkqfWuJxUD486iwvyOnJKZioAn27do1a5SAJQkDez+vOo161UdfAylCInQ15OJi/94LzonP3llTV6Vi6SABTkzahuadI9FdVAeD3WvBw/91x+hreFSUJ55LqB5PfoQLf26XTPaqPhaCKtnIK8GU1bsKHB0qTt0lOYclOBbqlLi8rLyeQ/z+rC1t0VTF9YwgWPzuGnM5cr0EVaKQV5MykqKeP5+Rsa7BtzbnctTSqeGF3QnYmX9SE9NYn9VbW8unQzD89e6XVZInISKMibyR0vLaW8sgYIrxU9/qLeTBh+usdVSaLK9vu4d8QZ/GTEgcc6b3+63cOKRORk0fCzZhAIhvgqMgUrwP/99kCuHpx7hCtEWsYrSzZFX/t9SQSCIT3qEWllmqVFbmZXmNlaM1tnZvc3cvy7ZrYi8rHAzAbVO7bRzD4xs2VmFn+LjBMeblZde2D7tY+3eFeMSD2PfXsQ6Snhf+aBfdVMnlvscUUi0tyaHORmlgz8ARgJ9AfGmln/g07bAFzinBsI/A8w5aDjlzrnBh9u0fRYFQiGePzdz5i9PBzcqclGfo8O/Oyqg799EW/k98zijYkX0SEjPLZ82kcbtUqaSCvTHC3yIcA659x651wImAGMqn+Cc26Bc64ssrkQOLUZ3tdzkz8oZtKcz/lyT/i2en6PDrx2xwXq4CYxJS8nk7/dMYy0FKOyupaiL3Zx38xlXpclIs2kOYI8F9hUb3tzZN/h3Aq8VW/bAf8ys6VmNv5wF5nZeDMrNLPC0tLSJhXcHALBEP9cvjW6ndshnUeuG3iEK0S8k5eTyVnd2ke3P9teruFoIq1EcwS5NbLPNXqi2aWEg/y/6+2+wDmXT/jW/J1mdnFj1zrnpjjnCpxzBTk53s+UNnnuOrburohu9+roV0tcYtpjowfRNj28Slp5ZQ0/fnWZwlykFWiOIN8M1F8N5FRg68EnmdlA4DlglHMuukC3c25r5PMO4HXCt+pjWiAY4o1PtgEHnov/8toBHlclcmR5OZn8/c4LGZjbDl+y8f7aUqYt2Oh1WSLSRM0R5EuAPmbW28x8wBhgdv0TzKwH8BrwPefcZ/X2+82sbd1r4D+AmJ61IhAMcffLRWzZFW6NV9U4/vOsLmqNS1zIy8kkGKohVBO+afb+mu1qlYvEuSYHuXOuGrgLeAdYDbzqnPvUzCaY2YTIaT8HOgJ/PGiYWWdgvpktBxYDbzjn3m5qTSdLIBjiR698zIfF0RsKDO2dreVJJa489u1BpCSFn4it2LKHi3/zvxSVlB3lKhGJVeZco4+zY1pBQYErLGz5IeePv/sZk+Z8DsDA3HZc2q8T44b11gQbEnd+Pmsl0z8qiW736tiGufdd6mFFInIkZrb0cEO0NUXrMSouLecvkeVJ69w7oq9CXOLSPZefQZd2adHtiqoajS0XiVMK8mP0k1eXRceLA5Ttq/KwGpGmyfb7eOkH59GrYxsAvtxTyQN/W+FxVSJyIhTkx6C4tJxV2/ZEt9NSkpg05hwPKxJpurycTF674wJyO6QDUFFVyzMfFKvzm0icUZAfRVFJGVdOmkdldbgvQadMH29OvIj8nlkeVybSdNl+H984uysAKcnGr99aoyFpInFGq58dxX1/XU5lvRVR+nRuq6Fm0qpMGH46HTPT2FK2n6IvdrFw/U6tkiYSR9QiP4KikjJ27ztwm7FdeoomfpFWJ9vv4/ZL8khPDf86WLQhwKNvrfa4KhE5VgrywyguLeeGZz/iq2C4U1tGahKv36kFUaT1yvAduEE3s3CzxpaLxAkFeSMCwRA3v7CYishz8XbpKbx023kKcWnVxg3rxcDc8MIqDhj34iJ1fBOJAwryRkyeW8ymsv1AeFWzufddqs5t0upl+31MvWUIbdPCC6vsrajh5hcXK8xFYpyC/CBz1+zg+Q83RLe/cXZXdfqRhJHt9/HNc06Nbq/YvJvJc4s9rEhEjkZBXk8gGOKHLy2lpjZ8S71LuzQmDD/d46pEWtY9I86gf9e20e1V23Z7WI2IHI2CvJ7JHxSzv+rAULOeHf1qjUvCyfb7GNG/MwDJSbBjTwWPvLlat9hFYpSCvJ7CjQdWNcv2p/LIdWd7WI2Id8YN601Wm1RqauGzHUGm/Hs9kz/QLXaRWKQgJ3xL/aczl7F8U/gWYvuMFGZOGKZe6pKwsv0+rh7YtcG+VxaVaGEVkRikIAdmFm7i1aVbqHGQkmS8dofGi4vcM6IvN53Xg8jS5eytrOG+mcu9LUpEDtEsQW5mV5jZWjNbZ2b3N3LczOz3keMrzCz/WK9tCef2yiYtOfzbqn9XTcEqAuFW+S+vPZsbh/aM7vtk0y61ykWOIhAMtegCRE0OcjNLBv4AjAT6A2PNrP9Bp40E+kQ+xgNPH8e1J1UgGGLijI+prHG0z0jhca1qJtLAPSMOrF1e5WD89EJ1fBM5gslzi/n1W2tabOhmc7TIhwDrnHPrnXMhYAYw6qBzRgHTXdhCoIOZdT3Ga0+qmYWbopO/9NWCKCKHyPb7uGZQLgB+XxLFpUGeeO8zj6sSiV2FJTsbfD7ZmiPIc4FN9bY3R/YdyznHci0AZjbezArNrLC0tLTJRdcZXdCd8RefxgV5Hfn1twY229cVaU0mDM/jgZH9yExLBeDPC0s0F7vIQYpLy7nmyfks+6Ju7gVrkfdtjiBvrFJ3jOccy7Xhnc5Ncc4VOOcKcnJyjrPEw8v2+3jwyjN56QeaS13kcOpWSMvNygCg1mkudpGDPfjaJ6zYspu62UgKWmhq7+YI8s1A93rbpwJbj/GcY7lWRGLEY6MHkZYS/rWxt6JG07eK1PPV3v3R1wNz27XYzKDNEeRLgD5m1tvMfMAYYPZB58wGbor0Xj8P2O2c23aM14pIjMjLyeSVH5xHZOly5n++w9uCRGJEcWk5G3aGgzwlyZh6y9AWmxm0yUHunKsG7gLeAVYDrzrnPjWzCWY2IXLam8B6YB3wLHDHka5tak0icvLk98yKrl2+cec+LvvtXD0vl4RWXFrON//wIbUOzODx6we16PTe5lyjj6RjWkFBgSssLPS6DJGE9dOZy3h16Zbodtd26Xz04GUeViTijUAwxBWPf8CO8nB/ke5ZGcz77683+/uY2VLnXEFjxzSzm4gct/uv7E9ejj+6/VWwUhPFSEJ6ePbKaIgDTPJgLhIFuYgct2y/jyk3FeBPC/8KqapxXPPUPIW5JJSikjJmLd8W3R41qCv5LdRTvT4FuYickLycTL475MD0rcHKWh6atdLDikRaTiAY4vtTF0e3+3dty0PXDPCkFgW5iJywCcNPx1fvt8iC4p38Y9mWw18g0krMLNzE7v3VAHTK9PHn285r0Q5u9SnIReSEhW+xn0tKZIm0Wgc/+atWSJPWKxAM8cgbq5i1bGs0QHvnZHoW4qAgF5EmGt6vE+/cezEZqeEwD1U7rnny33peLq3So2+uYsq8DazatodaICM1iUeuO9vTmhTkItJkeTmZnJ3bAQjPsbxiy14e+NsKT2sSORne/nR79HXntmm8dJv303sryEWkWfz6WwPJqXd7cfW2PWqVS6tSVFLGvsrq6Hb/bu086aV+MAW5iDSLvJxMbjjvQC/2vZU1apVLqzF3zQ6ue3oB1ZE51Nqlp/Czq/p7W1SEglxEms24Yb1on5ES3VarXFqDQDDE7X9a0mDfmHN7eH5LvY6CXESaTbbfx4s3DyEtOby9t7KG7z2n5U4lvk1bsIHKmgPb+T06MGF4nncFHURBLiLNKr9nFm/ecwm+5HAv9q27K7TcqcStQDDE0pJd0e1u7dN5bty5ng43O5iCXESaXV5OJgNy20e3p3+0USukSdwJBEPc/MJC5q/7ihSDgbnt+dNtLbc86bFSkIvISfHY6EHkdkgHoKK6lnEv6ha7xJeZhZtYsWUvANUOgqHqmHkuXp+CXEROirycTL79te7R7b0VNdz8gsJc4kMgGGJfqIb+XcKr/Pl9STz27UEeV9U4BbmInDTjhvViYG676PaKLXv0vFziwhPvfsakOZ+TkpzCxMv6MO+/L4uJMeONaVKQm1m2mb1rZp9HPh/yXZpZdzN738xWm9mnZjax3rFfmNkWM1sW+biyKfWISGzJ9vu4tF+nBvveWrlNrXKJaYFgiL9HFv9ZsWU3bXzJMfdcvL6mtsjvB+Y45/oAcyLbB6sGfuycOxM4D7jTzOqPon/cOTc48vFmE+sRkRgzblhvxl98Gl3bhZ+Xbyrbz6in5mt8ucSkQDDEbdOWsKciPINbbod0Rhd0P8pV3mpqkI8CpkVeTwOuPfgE59w251xR5PVeYDWQ28T3FZE4ke338eCVZ/LGxIvo1v5AmGt8ucSauhAv+mIXAGnJxvRbY6+X+sGaGuSdnXPbIBzYQKcjnWxmvYBzgEX1dt9lZivM7IXGbs3Xu3a8mRWaWWFpaWkTyxaRlpbt99E9OyO6HR5fvs7DikQaeuK9z6IhDnBWbruY7KV+sKMGuZm9Z2YrG/kYdTxvZGaZwN+Ae5xzeyK7nwbygMHANuC3h7veOTfFOVfgnCvIyck5nrcWkRjxyHUDSUux6Pbyzbs9rEbkgKKSMl5eVBLdbpeewmOjB3tX0HE4apA75y53zg1o5GMWsN3MugJEPu9o7GuYWSrhEH/JOfdava+93TlX45yrBZ4FhjTHNyUisSkvJ5NXfnA+7dPD87Gv3raHn89aqVvs4rmJMz6mujb8ulNbH6/feUFctMah6bfWZwPjIq/HAbMOPsHMDHgeWO2c+91Bx7rW2/wmsLKJ9YhIjMvvmcX7911Ku/QU9lRUM/2jEqYt2Oh1WZLAAsEQNbUuun1G5/i4pV6nqUH+KDDCzD4HRkS2MbNuZlbXA/0C4HvA1xsZZvYbM/vEzFYAlwL3NrEeEYkD2X4fp3c68Ity3uelapWLJ+o6uG3dXUEbn5HfowMPjzrL67KOiznnjn5WjCkoKHCFhYVelyEiTVBcWs51f/yQ3fvDw3yu/9qp/GZ0bM6cJa1TcWk5N7+wmE1l+6P7HhjZj9sviZ2VzeqY2VLnXEFjxzSzm4h4Ii8nk+8U9Ihuz1y6WQurSIu6d8bHDUI8v0eHmB8z3hgFuYh4ZsLwPDIji5c74LbpS3SLXVpMSWBf9HX7jJSYW570WCnIRcQz2X4f028ZSrf26bRLTyEQrGJm4Savy5JWrqikjPMeeY9QVQ0ZKUbntmm8ePOQuAxxgBSvCxCRxJbfM4sFD1xGcWk5D772CXNWb+fy/p3jqtewxJeJMz7myz2VAGT7U1n0fy73uKKmUYtcRGLCe6u2s2hDgMUbyxgz+SPdYpeTIhAMEYoMGE8Cfhcnk74ciYJcRGLC6ILuJCeFZ30rjQwJUphLcwoEQ9z9chHb94Zb4zee35Ph/Y44s3hcUJCLSEzI9vu46uwu0e2iL3bxxLufeViRtCZFJWVc+n/f58PindF9WW3i85n4wfSMXERixkPXDGBT2T6KvgjPwf6nhSVUVFVz/5X947YjknivuLSc7zyzgKrIFKzd2oeXJh03rJendTUXtchFJGZk+308N25IdLlTB7y6dIumcJUTFgiGuHXqkmiIA1w1sBv3jjij1fxxqCAXkZiS7ffxp9uG0ib1wK+nlxZupLi03MOqJF5N/qCYjTsPjBdPTTYmDI+9mduaQkEuIjEnLyeTs0/tEN3+KljFg6994l1BEpcCwRD/XL4VgIyU8KQvz36voNW0xOsoyEUkJj1y3dkkHVi6nKUlAU3hKsfl0bdWs3V3BQD7q+GO4ae3il7qB1OQi0hMysvJ5IVx50a3q2vhW08vYO6aHR5WJfGiqKSMmYWbo9vdszLich71Y6EgF5GYNbxfJ6befCDMw/OxF2p8uRxRcWk5331uIXVre3bK9DH1lvidgvVoFOQiEtOG9+tE/y4HpmutrnXc/XKRwlwO66FZK9lfr5v62KE9W/WUv00KcjPLNrN3zezzyOesw5y30cw+MbNlZlZ4vNeLSGK7sE/D55ofFu/k0TdXe1SNxLKikjKWbdoFQP8ufiZe1qfVjBc/nKa2yO8H5jjn+gBzItuHc6lzbvBBC6Mfz/UikqAmDM/jgZH9mPPjS0hNDveA+/vHm9UqlwaKS8u54dmFlFfWAJCdmd6qxosfTlODfBQwLfJ6GnBtC18vIgkg2+/j9kvyyMvJpG16eELKUC1cNenfCnMBwkPNbpiykIrIgiht05N5eNQAj6tqGU0N8s7OuW0Akc+H69fvgH+Z2VIzG38C12Nm482s0MwKS0tLm1i2iMSr9JTk6OuteyqZPLfYw2okFgSCIX70yoHFUJIN/n7nha36uXh9Rw1yM3vPzFY28jHqON7nAudcPjASuNPMLj7eQp1zU5xzBc65gpycnOO9XERaiaduyCcj9cD2O59+qVZ5AisqKePy381l/roDi6F862unJkyIwzEEuXPucufcgEY+ZgHbzawrQORzowM8nXNbI593AK8DQyKHjul6EZE6+T2z6NbBD0BKklES2Mc3//ChpnBNQEUlZVz/zEcEglVYZPKgC0/vyP0jz/S2sBbW1Fvrs4FxkdfjgFkHn2BmfjNrW/ca+A9g5bFeLyJysMe+PYi8HD/XnZMLQElgH7dO1frliWbijI+prg2PFncO8nL8/H5sfqvv3Hawpgb5o8AIM/scGBHZxsy6mdmbkXM6A/PNbDmwGHjDOff2ka4XETmS/J5ZzPnxcO6/8kzye7QHYOPOffzwz5osJlEUl5ZTti/837pNahJDe2cz5abWN4/6sWjSeuTOuZ3AZY3s3wpcGXm9Hhh0PNeLiByLbL+Pi/p0iq5fvmhDGTe/sIiptwxNyF/oiaKopIzvTF5AVWTqtn1VtXy9X6eEei5en2Z2E5G4Nm5Yr+j65QArtuxhZuEmDyuSk+2ul4uiIQ5wQV7HVjuP+rFQkItIXMv2+7hqYLcG+y7v39mjauRkCgRDPPLGKkojw8wA/GlJPHlD4j0Xr09BLiJxb8LwPE7xHxiT9t1nF/KdZz5ST/ZWZvLcdUyZt4GqSAc3X7Lxp1vOS+gQBwW5iLQC2X4ff5kwjA6RAeZf7qlk0YYA//n4v7XsaSsRCIZ445NtDfbdPKw3+T21RIeCXERahbycTP52xzC6Z2UAYIRXSrv31WWe1iVNE76dvpr/+N1ctuyqAKCNzxh/8WlMGJ7ncXWxQUEuIq1GXk4ms+66kAdG9uOCvI4A7N1fpVZ5HJs8t5gp89bzVbAquu/s3CwevPLMhL+lXkdBLiKtSt0CK5t37Qeg2sGdLy/1uCo5UYUlgQbbnTJ9PHLd2R5VE5sU5CLSKv3u+sHR18FQLdc8OY9H3litCWPiSHFpOZ9t3xvdzu/RgbfvvSRhx4sfjoJcRFql/J5ZPDlmcHR7xZY9TJm3XmPM40AgGOLnf/+EkY9/EF1bvFv7dJ4bd65upzdCQS4irdbVg3PJ79Ehut2rYwbvfPqlWuYxrG5J0ukLvyBUe2D/VYO6KcQPQ0EuIq3aY6MHkZIUXhpr4879FH2xiynz1jNtwQaPK5ODBYIhfvjnpQ2WJE0yuOm8nky4RD3UD0dBLiKtWl5OJs/dVEBGasNfd0tLdqlVHmMefWs1izY07Nx249Ce/PLaAWqNH4GCXERaveH9OvHh/ZeRlmLRffPXfcX5v35PQ9NiQN0z8ZlLNzfYP6RXFveMOMOjquKHglxEEkK238eoQbkN9lVWO74/dQlFJWUeVSVFJWUMf+x9pi/8AldvIZTcDulM/l5iLkt6vBTkIpIw7r/yTG46vydJBxrmOGD05AUKc4/c99fl7KmobrAvq00q02/VUrTHSkEuIgkj2+/jl6MGcON5PRvsr3Fww7MLtchKCwoEQzz+7mcM6NauQRC1TUvmrz8cprHix6FJQW5m2Wb2rpl9Hvl8yOz1ZtbXzJbV+9hjZvdEjv3CzLbUO3ZlU+oRETkW91x+BkN7ZzfYV1Fdy/eeW6QOcC1k2oKNTJrzObOWb6NulFm7jGT+fteFCvHj1NQW+f3AHOdcH2BOZLsB59xa59xg59xg4GvAPuD1eqc8XnfcOfdmE+sRETmqbL+Pp2/8Gjed35P2GSnR/Vt3VzDqqflqmZ9ExaXlXPfH+Uz54PMG+1OSjNfvUIifiKYG+ShgWuT1NODao5x/GVDsnCtp4vuKiDRJ3W32939yaYNJYzaV7eehWSu9K6yV+/nfV1L0xW72Rx6LG9A2PZnnbipQiJ+gpgZ5Z+fcNoDI505HOX8M8MpB++4ysxVm9kJjt+brmNl4Mys0s8LS0tKmVS0iEpHt9/HcuHPp1j49um/xhgDXPDlPLfNmVlxazuc7Gv5MRxecyie/uILh/Y4WH3I4Rw1yM3vPzFY28jHqeN7IzHzANcDMerufBvKAwcA24LeHu945N8U5V+CcK8jJyTmetxYROaJsv4+rBnUDIDnJCNU4VmzZw+W//UDjzJtBcWk53/zDfK544t/s2FsZ3X9p3xzuH3mmh5W1DilHO8E5d/nhjpnZdjPr6pzbZmZdgSP9Hz8SKHLOba/3taOvzexZ4J/HVraISPOacEkeHf0+Vm3dzazl24Dw0LTxfypk4YOXayjUCQoEQ9z8wmI2le1vsD+3Qzq/vX6wfq7NoKm31mcD4yKvxwGzjnDuWA66rR4J/zrfBPRgSkQ8UbeO+UPXDKCjPzW6P1TjGPXUfB55UwutnIgn3v2sQYhnpBj5PTponHgzamqQPwqMMLPPgRGRbcysm5lFe6CbWZvI8dcOuv43ZvaJma0ALgXubWI9IiJNku338eqEYbRLP3DDclPZfqb8ez3X/fFDPTc/DnPX7GD6woZ9m/N7ZvPaHReoY1szMld/Trw4UVBQ4AoLC70uQ0RaseLScm6duoSNO/c12J/bIZ1/3H2RWpNHUFRSxt2vFLFlV0WD/Tl+HzMmnK8QPwFmttQ5V9DYMc3sJiLSiLycTF674wKu/9qpDaZ03bKrgtumLeHxd9fqVvtBAsEQP5+1ktHPLDgkxEcN6so7/3WJQvwkOGpnNxGRRJXt9/Gb0YMoCexrsLxm0Re7KPpiF/tDtTz4DfW6hvAdjPHTCykuDR5y7KbzevDLa8/2oKrEoBa5iMhRPHLd2QztncUp9TrBAUyZt578//mXFlwBHnztE4pLg6TWv30BXHFWZ+4Z0dejqhKDglxE5CjycjL5y+3D6Nul3SHHAsEqrn/mo4QM8+LScr7/4mJeWVRC4cbwHYuq2nC/q6w2qbz2w2FairQFKMhFRI7RL68dQH6PDnRrn86gUw+EenWt49uTFyTc5DH3zVzO+2tL+dnfV1JTr9/0had3ZM6Ph5Pf87CTdUozUpCLiByjug5wCx64jCvP7tbgWK2D709dwk9nLmv1neCKS8u5ctI8ir7YBdAgxIf2zub3Y/PVCm9B6uwmInICRhd0Z97npcxftzO6zwGvLt3Cll0VPHlD6wuzuWt2cM+rH5OalERp+aF/rFx4ekeFuAfUIhcROQHZfh+/H5vPBXkdDzn2YfFObp26mKKSMr7/4uK4n0QmEAzx+LtruXXaEnbtqz4kxDPTkhl/UW+FuEc0IYyISBMEgiEmz13H7OVb+XJPZYNjRriV3j0rg1l3XRiXIVdUUsat05ZQtq/qkGOZacmc0TmTx0YP1vjwk+xIE8IoyEVEmkEgGGLUU/MPWRykTkqS8Y2zu/DQNQNiPtCLS8u5d8bHrP8qSKi6llBNw5xIAm48ryf3jDgj5r+X1uJIQa5n5CIizSDb72PqLUO4d8YyVm7dTe1BbaTqWses5dvo3C6dz3eU87Or+sdkK7a4tJzr/vghu/dXH3KsU6aPtNRkJo05Rz3SY4ha5CIizay4tJxf/XMVa7/cy9bdDacqbZOazL6qGlKTjU5t03hybL6noRgIhnji3c94b/V22qansOGrckI1Dc/JSIH91TDxstO5V5O7eEJzrYuItKC8nExe/P4Qnrohn14d23DFWZ2pm+9sX1U4JatqHFt2VXDDswv5zjMftWiHuEAwxDMfFBMIhpi2YCPTF5awdXcFa7cfGuIAqSkpPDCyH+OG9W6xGuXYqUUuItIC5q7ZwcS/fEzndmlsDuxnX1VttDMchJ+hX9znFNZu34vfl8y+UC2/unYAa7fvZXRB9yY/iy4qKWPijI9pm5bM7opqtuyqIMfvY1dFFVU1jeeALxnSUpJ5cmw+w/t1atL7S9Oos5uISIx45oNifv3WGgC6tEs7pKd7YzLTkjntlDYM7pFNVhsf44b1ItvvIxAMMbNwE5f378x7q7YzuqA7ZftCPDRrJf27tueKAV341Rur2LprP6V7KzlMXgOQlpJEZXUt+T06UNAzmwxfcvR9xHsnLcjNbDTwC+BMYIhzrtF0NbMrgElAMvCcc+7RyP5s4C9AL2AjcL1z7qgTFivIRSRehW9nbwCMS87IOezQriMZmNuOX1wzgB+98jGbd+2P/kGQbNC3cyarvgzfpk82jhjedTLTkimvrOHSvjn89vrBCu8YdDKD/EygFngG+EljQW5mycBnwAhgM7AEGOucW2VmvwECzrlHzex+IMs5999He18FuYi0FsWl5dw3czlbd+0nVF1DYN+hvcUb40uCUO2h+1OSoLqR/fUlWXhK2WSDvl3a8qtrz2bJxkCz3MKXk+OkDT9zzq2OvMGRThsCrHPOrY+cOwMYBayKfB4eOW8aMBc4apCLiLQWdfO3w4HJZd5a+eVhx6PXaSzEAWqOEuLdszKYesuQ6K34uuDWcLL41RLjyHOBTfW2NwNDI687O+e2ATjntpnZYXtTmNl4YDxAjx49TlKpIiLeyfb7ePAb/Zkw/HSmLdgIOK4ZnMt9M5dT9MUuUpOgbXrKIa32uhY2hDvP1X/efXqOn4/WB8jv0Z6dwWoeHnUWeTmZ5F0Se2PY5cQcNcjN7D2gSyOH/o9zbtYxvEdjzfXjvp/vnJsCTIHwrfXjvV5EJF5k+33cO+KM6PZz485lZuEmRhd0B2Dagg3sr6oF58jwpXBO9w784h+fMqR3Nl3bp3PN4NxDWtzSeh01yJ1zlzfxPTYD3ettnwpsjbzebmZdI63xrkBiLeYrInIMsv0+br8kL7rd2KQscw8aHqYWd+JoiQlhlgB9zKy3mfmAMcDsyLHZwLjI63HAsbTwRUREJKJJQW5m3zSzzcD5wBtm9k5kfzczexPAOVcN3AW8A6wGXnXOfRr5Eo8CI8zsc8K92h9tSj0iIiKJRhPCiIiIxDjNtS4iItJKKchFRETimIJcREQkjinIRURE4lhcdnYzs1KgpBm/5CnAV8349eJRov8MEv37B/0MQD+DRP/+IXZ/Bj2dczmNHYjLIG9uZlZ4uN6AiSLRfwaJ/v2Dfgagn0Gif/8Qnz8D3VoXERGJYwpyERGROKYgD5vidQExINF/Bon+/YN+BqCfQaJ//xCHPwM9IxcREYljapGLiIjEMQW5iIhIHEv4IDezK8xsrZmtM7P7va6npZnZC2a2w8xWel2LF8ysu5m9b2arzexTM5vodU0tzczSzWyxmS2P/Awe9romL5hZspl9bGb/9LoWL5jZRjP7xMyWmVlCrkplZh3M7K9mtibyO+F8r2s6Fgn9jNzMkoHPCC+hupnw2uljnXOrPC2sBZnZxUA5MN05N8DrelqamXUFujrnisysLbAUuDbB/h8wwO+cKzezVGA+MNE5t9Dj0lqUmf0XUAC0c85d5XU9Lc3MNgIFzrlYnAylRZjZNGCec+45M/MBbZxzuzwu66gSvUU+BFjnnFvvnAsBM4BRHtfUopxz/wYCXtfhFefcNudcUeT1XmA1kOttVS3LhZVHNlMjHwn1F76ZnQp8A3jO61rEG2bWDrgYeB7AOReKhxAHBXkusKne9mYS7Je4HGBmvYBzgEUel9LiIreVlwE7gHedc4n2M3gC+ClQ63EdXnLAv8xsqZmN97oYD5wGlAIvRh6xPGdmfq+LOhaJHuTWyL6EaolImJllAn8D7nHO7fG6npbmnKtxzg0GTgWGmFnCPGYxs6uAHc65pV7X4rELnHP5wEjgzshjt0SSAuQDTzvnzgGCQFz0m0r0IN8MdK+3fSqw1aNaxCOR58J/A15yzr3mdT1eitxKnAtc4W0lLeoC4JrIM+IZwNfN7M/eltTynHNbI593AK8TfvSYSDYDm+vdjfor4WCPeYke5EuAPmbWO9KxYQww2+OapAVFOno9D6x2zv3O63q8YGY5ZtYh8joDuBxY42lRLcg594Bz7lTnXC/CvwP+1zl3o8dltSgz80c6exK5nfwfQEKNZHHOfQlsMrO+kV2XAXHR6TXF6wK85JyrNrO7gHeAZOAF59ynHpfVoszsFWA4cIqZbQYecs49721VLeoC4HvAJ5FnxAAPOufe9K6kFtcVmBYZxZEEvOqcS8ghWAmsM/B6+O9aUoCXnXNve1uSJ+4GXoo07NYD3/e4nmOS0MPPRERE4l2i31oXERGJawpyERGROKYgFxERiWMKchERkTimIBcREYljCnIREZE4piAXERGJY/8PcdJMIwssco0AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.figure(figsize = (8, 4))\n", "plt.scatter(t, y, s = 1)" ] }, { "cell_type": "code", "execution_count": 11, "id": "958f87d7", "metadata": {}, "outputs": [], "source": [ "X = np.reshape(t, (-1, +1))\n", "Y = np.reshape(y, (-1, +1))\n", "Y0 = np.sin(t)" ] }, { "cell_type": "code", "execution_count": 12, "id": "8cf911a0", "metadata": {}, "outputs": [], "source": [ "xtr, xte, ytr, yte = ms.train_test_split(X, Y, train_size = 0.65, random_state = 22)" ] }, { "cell_type": "code", "execution_count": 13, "id": "0cb826ec", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "DecisionTreeRegressor(max_depth=10, random_state=45)" ] }, "execution_count": 13, "metadata": {}, "output_type": "execute_result" } ], "source": [ "DTR = tr.DecisionTreeRegressor(max_depth = 10, random_state = 45)\n", "DTR .fit(xtr, ytr)" ] }, { "cell_type": "code", "execution_count": 16, "id": "8fe8e247", "metadata": {}, "outputs": [], "source": [ "trPred = DTR.predict(xtr)\n", "tePred = DTR.predict(xte)\n", "\n", "#Mean squared error:\n", "trMSE = met.mean_squared_error(ytr, trPred)\n", "teMSE = met.mean_squared_error(yte, tePred)\n", "\n", "#Pearson correlation coefficient:\n", "trPCC = stt.pearsonr(ytr[:, 0] , trPred)\n", "tePCC = stt.pearsonr(yte[:, 0] , tePred)" ] }, { "cell_type": "code", "execution_count": 17, "id": "c35c850c", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "train MSE 3.7416187312649124e-05\n", "test MSE 0.00012556610215866233\n", "-----------------------------------\n", "Train PCC (0.9999631250661607, 0.0)\n", "Test PCC (0.9998749456077313, 0.0)\n" ] } ], "source": [ "print(\"train MSE\", trMSE)\n", "print(\"test MSE\", teMSE)\n", "print(\"-\" * 35)\n", "print(\"Train PCC\", trPCC)\n", "print(\"Test PCC\", tePCC)" ] }, { "cell_type": "code", "execution_count": 28, "id": "d297f4de", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 28, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmkAAAE9CAYAAABQhvWuAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABDOUlEQVR4nO3deXiU9dn28e81M5lJWJSwRJQloBLBBVASQBDccAHc6griUtH62NZqK221r7ba2qe1PrV1q3VBUKuIS7XahrqgRUBBCAjKLioogoLsa5KZud4/ZsAUIQRJck+S83McOWbmnnvuudLphNPfau6OiIiIiGSWUNAFiIiIiMg3KaSJiIiIZCCFNBEREZEMpJAmIiIikoEU0kREREQykEKaiIiISAaKBF1ATWjZsqV36NAh6DJERERE9mjGjBlfuXurnY/Xy5DWoUMHSkpKgi5DREREZI/MbOmujqu7U0RERCQDKaSJiIiIZCCFNBEREZEMVC/HpO1KeXk5y5YtY9u2bUGXUmuys7Np27YtWVlZQZciIiIie6nBhLRly5bRtGlTOnTogJkFXU6Nc3dWr17NsmXL6NixY9DliIiIyF5qMN2d27Zto0WLFg0ioAGYGS1atGhQLYciIiL1SYMJaUCDCWjbNbTfV0REpD5pUCEtKKtXr6Z79+50796d1q1b06ZNmx2Py8rKKn1tSUkJ1113XS1VKiIiIpmiwYxJC1KLFi2YNWsWALfddhtNmjThpz/96Y7n4/E4kciuP4rCwkIKCwtro0wRERHJIGpJC8h3v/tdbrjhBk488URuvPFGpk2bRp8+fTj66KPp06cPCxcuBGDChAmcccYZQCrgDR8+nBNOOIGDDz6Ye++9N8hfQUREpP5atAhGjgy0BLWkBWjRokWMHz+ecDjMhg0bmDhxIpFIhPHjx/P//t//4+9///s3XrNgwQL+85//sHHjRg477DC+//3va4kNERGR6jJnDvzud/DMM9C4MVxwAey/fyClKKRVYsbStdwzfhHXDyigR35utV//ggsuIBwOA7B+/Xouv/xyPvzwQ8yM8vLyXb5m8ODBxGIxYrEYeXl5fPnll7Rt27baaxMREWlQ3nsPfvtbeOEFaNIEfvYzuOGGwAIaqLuzUveMX8TED7/invGLauT6jRs33nH/l7/8JSeeeCJz5szhn//8526XzojFYjvuh8Nh4vF4jdQmIiLSILz7Lpx5JhxzDLzxBvzyl7BkCdxxB+TlBVqaWtIqcf2Agv+6rUnr16+nTZs2ADz22GM1/n4iIiIN2sSJqZaz11+H5s1T96+9NtCWs52pJa0SPfJzeeLKXjXS1bmzn//85/ziF7+gb9++JBKJGn8/ERGRBscdxo+H449P/cyeDXfeCUuXws03Z1RAAzB3D7qGaldYWOglJSX/dWz+/Pl06dIloIqC01B/bxERkR3cYdy4VGvZ1KnQpg38/Ofwve9BTk7Q1WFmM9z9G+ttqbtTRERE6qdkEl56KRXOZs6E/Hx48EH47nehwhjvTKXuThEREalfEgkYOxa6dYNzz4UNG2DUKPjwQ/if/6kTAQ0CDmlmNsrMVprZnN08b2Z2r5ktNrP3zeyY2q5RRERE6oh4HB5/HA4/HIYOTbWkPfUUzJ8PV1wBdWxd0aBb0h4DTq/k+YFAp/TP1cBfa6EmERERqUvKyuCRR6CgINWVmZMDzz0HH3wAF18Mu9l6MdMFWrW7TzSzDpWccjbwhKdmN0w1s2ZmdqC7r6idCkVERCRQxSPw6Y/iOE/GT6FnaD6H2bLUc3GHmeXwTim2wfGDQnBRDhR8DB9cSWLuNURu/SrY+vdBpkfLNsBnFR4vSx9TSBMREWkISkZjOAZcHH6DMEms3GFGGbxTBpsc2ofhrBh2cBjMdrw07LvevaeuyPSQZrs4tss1Q8zsalJdorRv374ma9prq1ev5uSTTwbgiy++IBwO06pVKwCmTZtGNBqt9PUTJkwgGo3Sp0+fGq9VREQkoxResaMl7bkt/Tl5xlRavbsK2+J4hzCcG4P8CuGsQkpIhLIyPuhUJtNrXwa0q/C4LbB8Vye6+8PAw5BaJ63mS6u6Fi1aMGvWLABuu+02mjRpwk9/+tMqv37ChAk0adJEIU1EROq/4hF4yWiWhtrSLvEZX0bzab3FsXe3MWTay9g28EMjeL8otIswJnkqNvguLu71zQaaTA85exL0xIE9eRm4LD3Lszewvr6MR5sxYwbHH388PXr04LTTTmPFitSvde+993L44YfTtWtXhgwZwpIlS3jwwQf585//TPfu3Zk0aVLAlYuIiNSA4hFwWzN8+kjME+THlxLeHOfAcfMJ3bMBm1iGdYjgVzXGhjXC2kcwg4tC4/m/VxcEXX2NCDRkmtnTwAlASzNbBtwKZAG4+4PAOGAQsBjYAlwRTKXVy9350Y9+xEsvvUSrVq145plnuPnmmxk1ahR33HEHn3zyCbFYjHXr1tGsWTOuueaavW59ExERqVNKRkN67JlvSGJTymBGGVYOHBGBfjE4ILyjO3P7hknPJAfws0GdAyq6ZgU9u3PoHp534Ie1VM43FY9I/Z+m8AoYfFe1Xba0tJQ5c+ZwyimnAJBIJDjwwAMB6Nq1K8OGDeOcc87hnHPOqbb3FBERyURr/9iDZhsXU0qE2PokvFOamrGZBD8qC46LQsswkIpwi2nDobaCUFHq3+ZhAddfk+p6d23NKhkNnkjdVmNIc3eOOOIIpkyZ8o3niouLmThxIi+//DK33347c+fOrbb3FRERyTTNNi7G1iWJTd6IzSrHAeuWBcfFIDeEGSQIMeuKj+mRn0tB0AXXokwfkxaswivAwqnbahSLxVi1atWOkFZeXs7cuXNJJpN89tlnnHjiidx5552sW7eOTZs20bRpUzZu3FitNYiIiARl5dhr8dv2x3/YBH9pK37fJphdDsdkYdc1wc/MgeapgOZAuGg4PfJzgy671qklrTKD76rWFrTtQqEQzz//PNdddx3r168nHo/z4x//mIKCAi655BLWr1+Pu/OTn/yEZs2aceaZZ3L++efz0ksvcd9999GvX79qr0lERKS2tJj4BDZpK8yJQwSsVxQ/Ngr7hVItaaTHnBlY0VU18m9xXWDuGbVaRbUoLCz0kpKS/zo2f/58unTpElBFwWmov7eIiGSW9be1Yb/l62FSKbYgjkeBwigcG8Ubh1jkbSkIryDUqgBWLar28eCZzMxmuHvhzsfVkiYiIiLVr8LkuwVlp3HYmJXYh3GIAf2jWO8Y5KQWoE0QovOvNQZ7ZwppIiIiUq1W/uFoWm35GPs0jj9+D50/+ROeY/iJMSiKYjmW6s701IzNj/MvbFATAqpKIU1ERESqxx/yYctaWn2UwCaVwqcJrLHBKTGsMIpn2Y7dm8yA29ZjoIC2GwppIiIisu/cYfYqmFiGfZ6ApganZ8MxWXgkHc7SLWeGQ9FVQVec8RTSREREpOqKR8D0kTseujssiOMTywh9kcCbGQzOhm5ZEDEww/I6w6pFWFHDmQxQHRTSREREpOpKRqdukw5z46luzVVJaB6Cs7Oxo7IgbF+fbyHI7ws/mBpMvXWYQpqIiIjsWboFzRMOH5TD5DJsdRJvFcLOzYHDI3goHc62r3GG1cjOPQ2FQpqIiIhUrngETHkEZpdjk0thnbP6gP25+ZwfcOzhCxkWeZOnEycTHvQnLu7V/r9ft30PbNlr2haqlkyfPp2uXbuybds2Nm/ezBFHHMGcOXO+9fU2btxIx44dKS8vB2DDhg106NBhx2MREZF9UjwitXXT/9sPv/8BuG8T/GsbNDYYmkPu1Un+0u2vDDzqQIa3f4XDv/fIfwc0SLWe3bpGrWjfklrSaklRURFnnXUWt9xyC1u3buWSSy7hyCOP/MZ5/fr12+U+nX/84x8ZMGDAjsdNmzblhBNOoLi4mHPOOYexY8dy3nnnkZWVVaO/h4iINBBvj8JKSmFKGWxyaB/Gz4xhh4RxDDDCJMlbOIYnbr0/6GrrpYYZ0n78Y5g1q3qv2b073H13paf86le/oqioiOzsbO69995dnjNp0qQqv+VVV13FnXfeyTnnnMPo0aN55JFH9qJgERGRr81YupZ7xi/iV9se5OAXnsamlGFbHTqG4bwYnh9JbXju8HJ0IMce0pK8hWPUlVmDGmZIC8iaNWvYtGkT5eXlbNu2jcaNG3/jnKq2pAH07duXJUuW8NZbb5FIJHbZMiciIrInM5au5Ya/vsE97/yKQ6bPx7YBnSLQLwrt0lEhvdW3HdCFs38wNv1KtaDVpIYZ0vbQ4lVTrr76am6//XY++eQTbrzxRu6//5v/596bljSAyy67jKFDh/LLX/6yusoUEZGGongEPuFhjplSxoTpZVgZ0DmCHxeDg8Kpc9IL0C7Ov4iC4Q8FWm5D0zBDWgCeeOIJIpEIF198MYlEgj59+vDmm29y0kkn7dN1hw0bxi233MLQoUOrqVIREanX0jMu1yXasv/rC7EZZVAOHBGBfjE4IBXOLm/3GtcPKKBHfq62bgqIQlotueyyy7jssssACIfDvPvuu9Vy3cmTJ3P++efTrFmzarmeiIjUU9uXw1hXDm+Xsv/MD7Ak+FFZWL8otAxv79HE8rrwxJW9Ai1XFNLqtB/96Ef8+9//Zty4cUGXIiIime61kTBpKz67HAPolkW8bzaLc9tREPocw7Giq7RcRgZRSKvD7rvvvqBLEBGRTFNxAdnBdxH/US7hiZtTuwSEwHpkQd8Y7Beia/IZbhl8OJ13Xt9MMoJCmoiISH1SMho8QbJ4JJ/9/Fnaz1uHZQG9otAnijcJYQYGzLv99KCrlUo0qJDm7pjZnk+sJ9x9zyeJiEi9sHLstbRc8DdseQImlRJaEKd9dAPWN4r3jmJNQriDG2BhTOubZbwGE9Kys7NZvXo1LVq0aBBBzd1ZvXo12dnZQZciIiI1ZMbStax99kecvPmftPo0jk0ugw/jeDZ81b8lrxceywVNJrHYD6ITy1ndZRh5Q7S2WV3RYEJa27ZtWbZsGatWrQq6lFqTnZ1N27Ztgy5DRERqyD3jF/HYnH9gk7bBJwnIMfzEGNYzSotYGe92+TX/3lzG9QMKCOfnkhd0wbJXGkxIy8rKomPHjkGXISIism+KR8D0UWz8vCWPv/YJ9mkCGht+Sgx6RCGaHnAG3D3k6GBrlX3SYEKaiIhInVU8AqaPTI0p+7Cc0MQymn6+DvYzGJiNd8/i2ejpJAb+kYu/ugdKRmvMWT2gkCYiIpJpHugNK+dDTjN820ZIxrH5cWxSKfZFEm9m2BnZ0C0LDxsf5g/hoh1bNt2ltc7qCYU0ERGRTLNyfup281psbhyfVAqrkniLECvOOojWR22EsOEGX3W+lAJNBqiXFNJEREQyyMqx19Iq4dgH5fikMmxNEmsVws/NYWGXDhwaXkHIDLcwoVvXaDJAPaaQJiIikilKS2n55Cjs7W2wzqF1iOQFjaBzmCeTp7K8z+0Uzv1fTt4yTmPOGgCFNBERkYA989YiDv79RRROmU1og+NtwjAwypiDB3JLfDhtGudwz5CjuSw/FwY9GXS5UksU0kRERGrZjKVruf2fc7lq5d0MmlrMee+UE9mcgPZh/MwYdkiYhIUZdtvfGRZ0sRIYhTQREZHakt78vPm21rzw7sfY1DJsq2MHh6F/I7x9BDNwh1ezBzIo6HolUAppIiIiteGB3rB0HkwtpcO0+dg2oFME+kWhXQQntQatA8uyOnDA0L8EW68ETiFNRESkpq1ciY95D5teBmVA5wj0j+Gtwzt2B9i+q7QVXUW7wXfRLqhaJWOEgi5ARESkvrrvyYnM69UJb3MAvF2Gd4rg1zRm4QUFxFtn8WH+EOy29VjRVWBhKLpKC9HKDmpJExERqWb/eGkKm3/zv/xg9r8JJZNY1yySfaMkW2bxrA+AwXfRuVd7Cra/YLB2CZBvUkgTERGpBjOWruWR0a9z9r8f56ySf2M4dM/CjmuE54YAI0KSi0NvQK/2QZcrdYBCmoiIyL5asIAWF53OXz9YSjIUgh5ZWN8oyf3CJM0IFw1PjTkrGQ1ahFaqSCFNRESkqtJLaFB4BQy+i9V3DiX3ib9j88rJj4D1imLHRhnb5HQuCo3n1eyBHDD0L/TIz029Xl2ashcU0kRERKqieARMHwlA8uVHmPDTqZy04B2IQrJPDHpHsSYGGIdd8TDh/Fytcyb7RCFNRESkKkpGw2dxfGIZocVx+mZPZ8vxTcjpaViOsdDbUmArCBVe8XXLmcg+UEgTERHZWYVuzTEtruPf94/lwbdLabxkC9bI4KQY0aIoZNuO9c0OC63Abl0TaNlSvyikiYiIVLS9W9OdxNiHOHTCK/xt2Ty8icGpMfyYKBb7OpxtZ5oQINVMIU1ERKQCnz4KW1iOTywlvDxJ0X7z+eL01ozv2psh2ROIZzUhO77h6xfctj64YqVeU0gTEREpHoFPG8WU+V04ZuJWsr8shWYGZ2Rj3bM4ILSFhxr9iC5DRqXGm1Wc5SlSQ8zdg3tzs9OBe4AwMNLd79jp+ROAl4BP0odecPff7Om6hYWFXlJSUr3FiohI/RSP4xc1wyZthVVJvEUIOy4KR2XhoVSn5rqmh5L70xkBFyr1lZnNcPfCnY8H1pJmZmHgL8ApwDJgupm97O7zdjp1krufUesFiohI/bO9BaxVAayYR/L9OJsmhdhv7Wa8VQg/N4dFXTpwWPhzAJZldaDdLbPRXE0JQpDdnT2Bxe7+MYCZjQXOBnYOaSIiItVj+qMQT+LjZmNvlxJa5zRtHWLjhbk0OSxOKGR0YgUzr1hCj/xc2gVdrzRoQYa0NsBnFR4vA3rt4rxjzWw2sBz4qbvPrY3iRESkntm6Fd7dBm+XYRsdbxPGBkaxThGakMAs1bUZLhqudc4kIwQZ0naevQyw8wC5mUC+u28ys0HAP4BOu7yY2dXA1QDt22vjWhERSXnvT5fS9d/PY1PKCG1O4u3DcHaM0o4xYpYAA8tpBts27tjuSSQTBBnSlsF/tSS3JdVatoO7b6hwf5yZPWBmLd39q50v5u4PAw9DauJAzZQsIiIZLz3ubFW783jxmQTfm/QUttXxg8PQvxGWn/qnL0YS0/IZksGCDGnTgU5m1hH4HBgCXFzxBDNrDXzp7m5mPYEQsLrWKxURkTrDJz6KTd1Ky3dHcnUpeKcI3i+KtUv9k+ekunLsG503IpklsJDm7nEzuxZ4ldQSHKPcfa6ZXZN+/kHgfOD7ZhYHtgJDPMg1Q0REJHOtXMncn95GwTObySqL8+lhB9Gm/wYiB4VwrxDO8rrAqkVa40wyXqDrpNUUrZMmItJwvPiv6XS+/WI6z1wMCeCICC8fewLr8/bnksjrGIYVXamxZpKxMm6dNBERkX3x1X2X0fyJZzjrvTihZBLrmgX9otAizBmUEDLDHLCQAprUSQppIiJSp7z8/ETKf/c7zp39KgZ49yzo2wjPDe1YNyBcNDx1R1s3SR2mkCYiInXCyh93ptWbH3PmnPLUNLIeWXifGDQLpcIafHO2plrQpA5TSBMRkcx2/2X4X8fSal45lgX0jsKxUWga2jEZANITAkTqEYU0ERHJSEt/dS7tnv0noYVxLAp+XBTvFYXGoR3naBFaqc8U0kREJHgVNz6fMQefuI38xQnIBj8+RqJnjKdjpxA+409c3Eu7ykjDoJAmIiK1b3so2z6of9ojsCSBP1aCLUlgjQxOjkFRFI8a5+eN45dnHK49NaVBUUgTEZHaVzIaPAHTHsE/SmATS+GzBDQxkqdkYz2yIJraE2Bx/hD+Mbxv0BWL1DqFNBERqX09vgtPPQgTS7HlSXw/g4HZLOregdOSfwSHa3odzE2DulAQdK0iAVFIExGRmvONbs1RLF9ZyIYXP6bzl1vxXKPsonxO6PwAK7Y553Q9iCVDjg62ZpEMoZAmIiI1Z3u35ruj8Dml2KRSDvpqPFubt+WGwT9hzVnn8djVfZkSdJ0iGWiPIc3MDgGWuXupmZ0AdAWecPd1NVuaiIjUacUjIJ6A98uIT4sR+WIr5IXw83I45eC/cEjr/bnjtMODrlIkY1WlJe3vQKGZHQo8CrwMjAEG1WRhIiJSh5WW4n99EJu8FdY7C1q3pdkFTWnTeT0LacdZh7fjbnVrilSqKiEt6e5xM/sOcLe732dm79V0YSIiUodsH3t25CVsfOwFmry1HNvoeJswDIoSOSSLM7P+wsjLi+iRn8vdQdcrUgdUJaSVm9lQ4HLgzPSxrJorSURE6hp/exQ2fSv+f/fTdLNDfhjOiUGHMBYyDuPzHQFNRKqmKiHtCuAa4H/d/RMz6wg8WbNliYhInbB+PZ/ffictH9hCbGsZyw/OI94vh/b5awEozdqP7PgGLK+LAprIXtpjSHP3eWZ2I9A+/fgT4I6aLkxERDJMheU07t54Ppf85SxaTP+KNqXgnSKs79eCA9qWURw9ia8ueYAe+blkB12zSB1m7l75CWZnAn8Eou7e0cy6A79x97Nqob5vpbCw0EtKSoIuQ0SkfrmtGWxK4FPKoKQMKwO6RKBfDA4M4w5mgIXh1jUBFytSd5jZDHcv3Pl4Vbo7bwN6AhMA3H1WustTRETqu3Tr2cdNBtHxla3YjDJIgB0RwY+LQV44FcyoENC2L1wrIvukKiEt7u7rbfu3MKXy5jcREakfxo+ESVvpOOspSALdsqBvlHiLLEIkCVnqHwTDoOhKGHxX0BWL1BuhKpwzx8wuBsJm1snM7gPeqeG6REQkQHPemsHGHs3xezfg75VDtyw+/uGhxM9qzMsHnsHs4R8T6nkVWBgrugpuW6eAJlLNqjImrRFwM3AqYMCrwO3uvq3my/t2NCZNRGQvFI+A6SMB8FUJyiYlyJqzDQuDHROFvlHYL0SCELOu+FizNEWq2bcek+buW0iFtJtrojAREQnY9JHwRQImlWLz4kSzwHvHsD5ZeOMQWKo7M1w0XAFNpBZVZe/O/7CLMWjuflKNVCQiIjVj+xIarQpg1aLU7ax5+MSt2MI4HgWOi7Kudy4vx45jWOQN1nQeRt6Q+4OuXKRBqsrEgZ9WuJ8NnAfEa6YcERGpdtvDmSdSj1fOh0/j+JMl2OIEZIMfH+PDog6c1+huHhvei8vTLWZ5AZYt0tBVpbtzxk6H3jazt2qoHhERqW7bA5qDL4ljE0thSQJvFCJ5cjbvDzyP878aykHNcnhsyNHq0hTJEFXp7mxe4WEI6AG0rrGKRERk31TYGYDBd0GP75Ic+whrJsZo+dkGvIlhp8ZIHpPNwFYv8foNx/NR0DWLyDdUpbtzBtuXwUl1c34CXFmTRYmIyD5Iz9Rk2iMkn3yItROjtFi+idKmOfzylO9z2DErGJI9gdeyB3LHeV2DrVVEdqsq3Z3aXUBEpK4oHgHuMC+OTyol9GWS5rlbufH0H/HCkSdRHs6if6eWXHLliwwKulYRqdRuQ5qZnVvZC939heovR0REqmTnLk3g6bc/5sKRDxKevA2+SpJsEcbOyYYjo/SKfsnMRs1onJ3F9QMKAi5eRKpit4vZmtnoSl7n7j68Zkrad1rMVkTqvV83T00GsDD84gvG/fxODn/8ATqsW4HnhbD+MZKdI1jIduytyW3rAy1ZRHZtrxezdXftkCsikqkKr4Apo9g4tzFN8rIZtN75qnUzrv7OzfTtMp9hkTcJm2M7lrm0Si8nIpmnKhMHMLPBwBGk1kkDwN1/U1NFiYhIJbZsYcb4JrR7NEzexmXQNgyDozQ7xHm97FhOOOtqIr3aV9juKb35uYjUKVVZguNBoBFwIjASOB+YVsN1iYgIpIPWo4BDKWx5vwk5b62gx2bH88NsPKc5jTvEIQSvZQ/i+Sv7fL3O2eC7tOm5SB1WlZa0Pu7e1czed/dfm9ldgCYNiIjUhpLRsC0J75bBu2U02roBDglD/xjWPkIT4lh6rJlma4rUL1UJaVvTt1vM7CBgNaBlOUREakq6m9K3JGFKGUwvw0rBCyLQLwZtwjtOtbwuARYqIjWpKiHtX2bWDPg/YCaphW0fqcmiREQatAmPwjvbsOllUA50ieD9YnBgGCO9urhmaorUe5Wtk1YMjAH+5O6bgb+b2b+AbHfXXwcRkX1VYa2zGUfewhPPTebSSc/SfdxGwokkHBGB42KQF2ZpJJ+8pjEarVuk1jORBqKylrSHgSHAn83sP8DTwDgFNBGRarJ9QsDrj/DRbQu5871XCbnzyZHt6NhvDU81O5U/hr/HLwZ14eJe7YOuVkRq2W4Xs91xglkOcBapwHYsMA542t1fr/nyvh0tZisiGatC65m/+jA2qRR/v5wyi/D8UQP4a+8LmJA3goglSRIidNvaoCsWkRq214vZbufuW4FngGfMrCvwOHA5EK70hSIi8k3TR8KqBP6L+2BOeeovaWGUrD4xhu73Dgc0aUHk8OFQMppQodYUF2nIqrJO2gHAhaRa0g4EngP0l0NEZG/Nng3PbYF5cciCNb1b0PzYUmgSIpTeEGDAlnEweI3WNxORSicOfA8YChxGal20n7v727VVmIhInfdAb1g5H/88waaJRtNFm0jEQoT6xaBXFrmN4nxs+RzMZ1+/Rq1nIpJWWUtaH+AOYLy7J2upHhGR+qPkA5hYin2UoEk2cEIM7xkjmRMiTBLDOWR7QLMw3Lom0HJFJLOEdveEu1/h7q8poImI7AV3Fj39EmsOaQ2jt+ArkvjJMT687jDi/XN4qcnprO48LH2yQV6XVEBTC5qI7KRKG6yLiEglikfA9FEkPkoQemsLBZ8l8CYGp8WwHlHikTAv9HqRmwZ14bwdL7o/wIJFpC5QSBMR2Rf/vAF/6kFsUinh5UnY32BQNnZ0Fh4xHOPj9hdy0yAtQCsie6eyiQPNK3uhu+/z4AkzOx24h9Qk9JHufsdOz1v6+UHAFuC77j5zX99XRGSfJRLw97/j192LfZnAc434mY2JdA1B2MDAiq7CBt9FQdC1ikidVFlL2gzSW8QB7YG16fvNgE/Zx03WzSwM/AU4BVgGTDezl919XoXTBgKd0j+9gL+mb0VEghGPs+KaU2j90kTsqyS0DOHnZLPwiI5c1fR+Jm07H/NEapyZltEQkX1Q2cSBju5+MPAqcKa7t3T3FsAZpJbk2Fc9gcXu/rG7lwFjgbN3Ouds4AlPmQo0M7MDq+G9RUT2TlkZG/7nJDwvxoGPTsDC4OfnYN9vjHWLUpD1BZNvPAkrvEITAUSkWlRlTFqRu1+z/YG7/9vMbq+G924DFRcHYhnfbCXb1TltgBXV8P4iIpUrHgFTRvHp4s6E//0RbTasggND+EU5eEEEs1S3JvD17gCD71ILmohUi6qEtK/M7BbgSVLdn5cAq6vhvW0Xx3beSLQq56RONLsauBqgfXttRCwi+6B4BLzzCD6jDHunjPYbp1LSpgsHnbEZOzQMZpiFU92aAEVXKZiJSLXbbXdnBUOBVsCL6Z9W6WP7ahnQrsLjtsDyb3EOAO7+sLsXunthq1atqqE8EWmInn1zHok/PwB3b8JeLYUWIZKXNuL8YXfih0ZSAQ2+7tZUQBORGlKVDdbXANebWRN331SN7z0d6GRmHYHPSe0NevFO57wMXGtmY0l1ha53d3V1isi+Kx4BJaNTY8cG38Ws95cw7xe3M+jNZwlv24YfEsb6x/B2Ef6WOIXcaBZfFVxK3sIxO16jcCYiNakqG6z3AUYCTYD2ZtYN+B93/8G+vLG7x83sWlITE8LAKHefa2bXpJ9/EBhHavmNxaSW4NBIXBGpHiWjwRPw1qO8/bfNHPXi3+hetoXXD+3F/cdeyLnt3+GSyBus7jyMy4fcz+UAnIYWoRWR2lKVMWl/JvWX6WUAd59tZv2r483dfRypIFbx2IMV7jvww+p4LxGRHYpH4Bvj8E4pXlJOn/KHoEuEl/uczK3tf8baLXE6Hnkilw95nrygaxWRBqtKOw64+2epdWV3SNRMOSIiNWfG0rWsvO8yTp/0GjazDE8AR2Zh/aLQKswZzOCsX50WdJkiIkDVQtpn6S5PN7MocB0wv2bLEhGpRsUj8NdHst+kFhwz+xPMwbtm8eKxJ7KxeVMuDb8OGKGi4UFXKiKyQ1VC2jWktmZqQ2q25WvAPo1HExGpcemJAav3G0TuE88Rer+MQ20Di7p15JDjvuKppqdyW/wKzj7yIN7Y/DOuH1BAj/zcoKsWEdmhKiHtMHcfVvGAmfUF3q6ZkkRE9l2i+BHCk7fSfM5TeNjwoigcG6Vgv9VgIS4vyufywYODLlNEZLeqEtLuA46pwjERkWA80BtWzoecZvgna/ni7WYcOGcjngVzexXQuc8XRJpWfIGnZndqCQ0RyWC7DWlmdizQB2hlZjdUeGo/UktmiIhkhpXz4fMEPnE5tihO69hGZvftzBHHLuegRhsJk9qq5OvpT6a9NUUk41XWkhYltTZaBKj436AbgPNrsigRkSqbPBmeNZi/GcsxODGGF0U5IvsLIgbN2fz1ubetD65OEZG9tNuQ5u5vAW+Z2WPuvrQWaxIR+aaKOwS4488+DBO3YUsTrG2Sy2cDjuLIHkshBmP9VArzm1Pw2XOQ3RS2roO8LkH/BiIie6UqY9JGmtkF7r4OwMxygbHursWERKR2FI+A6SPBncSYhwhN3IYtS0BTI3lajOJu/Rma8xZvNj6T3Avv4+IdszQfCrRsEZF9UZWQ1nJ7QANw97VmpkW4RaTmbW89SyZgYTk+sZTwiiSb9m9E40FJ7OgsCIcYwluESTJgyzjIfzLoqkVEqkVVQlrSzNq7+6cAZpZPagyuiEjN2B7OEnGYF8cnlWIrk5BrrDjzQFp13cLCUBs6hZbzavZAAAaWvoJpMoCI1CNVCWk3A5PN7K304/7A1TVXkog0eO+Ogg+24ZPLsK+SWMsQ/p0c7MgIrW0zZnCYrcBuXcugoGsVEakhewxp7v6KmR0D9CY1g/0n7v5VjVcmIvVfhckAWz6cSM7qhZS+HyY2aT221uGAEGXnNyHUJcxiUi1npc0OpdH6j9RqJiL1XmXrpHV29wXpgAawPH3bPt39ObPmyxOReq1kNHgCpowiZ+Y27O1Ssjc4HBTCT83GDouQhfGTLv/h7iFHA9Ao4JJFRGpLZS1pI4DvAbtaktuBk2qkIhFpOI68hOTDj7DtHWi0aRveNgxnxLBDw2CWWoA2p9mOgCYi0pBUtk7a99K3J9ZeOSLSEDz75jwK/vciuk2dS2iL4x0aMfSM33FOwUwutNeA9O4AFoYbtUyjiDRMlXV3nlvZC939heovR0Tqs1nvf8JHv/o9p736NPtv2wSHhqFfjEbtIhzQ5XQuWvh/X88dt7C2bhKRBq2y7s4z07d5pPbwfDP9+ERgAqCQJiJVMuu9xXx40284/a2/0710C6936k2n/l+S32Y1lt5Q8+6FJ0OrAli1KBXOtPm5iDRwlXV3XgFgZv8CDnf3FenHBwJ/qZ3yRKTOeqA3/vFcSt9xupVspls5LO3Shgt738H8vIP5OPviHRueG5aaQLBqEdy6JtCyRUQyRVXWSeuwPaClfQkU1FA9IlIfLFuGP/4eNrOMWALsqCw4Lkr7lhuZX3ow53Q/iNACI9W3aVB05df7coqICFC1kDbBzF4Fnib1F3UI8J8arUpE6pyVY6+lxZQn+HxqHm1nfIJ5Eu+ahfWL4bkhzFItZkvuGJx6QXGFYDb4LnVviojspCqL2V5rZt8htdMAwMPu/mLNliUiGa3CIrQMvout17em1VtrsPfLaRfaiB8dxfpGSe4fImmhXS9Aq2AmIlKpqrSkAcwENrr7eDNrZGZN3X1jTRYmIhksvQitFz/EhhtHst+8DVgYvGeU53sOwFs05tzk68w+4Dv0+MEoLUArIvIt7DGkmdn3SO3V2Rw4BGgDPAicXLOliUim2rDuAJq+shhbEGe/LODYKN47ytomTbmn0QjuGXI0kfxcegRdqIhIHVaVlrQfAj2BdwHc/UMzy6vRqkQkI7366EvE/vA7TvhwAcTA+0WhVxRrHCLuIS7a/2km33B80GWKiNQLoSqcU+ruZdsfmFmEr5ebFJH6qHgE/Lp56hZY+Nw4Vh7amtOuOofey2ax6oRW+PVN+fDEQ/mwy8W4hZnQZDB3nNc14MJFROoPc688b5nZncA64DLgR8APgHnufnONV/ctFRYWeklJSdBliNRdv24OyTh84mxc2JWm097BGxt2bBQvjOIxS/0XnoW1rpmIyD4ysxnuXrjz8ap0d94IXAV8APwPMA4YWb3liUhGKB4B00bii8vhrTLs8wTh/d4jeVo2dkwWnmU7ltLAQlrXTESkBlUa0swsBLzv7kcCj9ROSSISiGQS/vYQTNyKrUji+xtTTj+a3scsxrJSewN4xXCm5TNERGpUpSHN3ZNmNtvM2rv7p7VVlIjUguIRMH0kJB2fF4dJpdjKJN48hJ2VDUdl0Sv8EbZ9c03Aiq5UOBMRqSVV6e48EJhrZtOAzdsPuvtZNVaViNQ4f3cU9kEZTCrDViehZQi+kwNHRiBk6X011XImIhKUqoS0X9d4FSJSe8rK4PHH2XJ/OY3XbsMPCMH5OdAlFc7i0f3JKlsPmFrOREQCtNuQZmbZwDXAoaQmDTzq7vHaKkxEqsn2LZyaHUrytTlsnQyNN2zlo4MKuOeki3jjkJ40ioa5ZfARXNyrPVlB1ysiIkDlLWmPA+XAJGAgcDhwfW0UJSLV6J1RULIVf6eE0CanUbswl532a87rOo+H4vfxcfsLKRj+UNBViojITioLaYe7+1EAZvYoMK12ShKR6vDenKXMvPF/GTZhC9lbytjSoRHZ5xpPtj2NFS2O56yNd2MkKfjsOUAhTUQk01QW0sq333H3eMUZXiKSgf6QD1vXUZ5oyqi3T+KiKS9y9LZN+KERlh7XgRNb389vsx7j0vAbXH5YB+CKVDeo1joTEclIu91xwMwSfD2b04AcYEv6vrv7frVS4begHQekQfpZU5hahk8rw8ogeVgW1i+KtQnjwMwrltDjsUPAE9opQEQkg+z1jgPuHq7ZkkRkn6TXOfONSda+EyV3xmYod8oOzybcL4tIa9uxya7ldaFHfm6q1UytZyIidcIe9+6si9SSJg2B39AMm7wFZpbjSfj4yPZc3fMXLDugPbfYKIZF3iBUNFxLaIiIZLh92btTRDJB8Qi8ZDSfrWlF28lLsFllqZayblnYcTE65q5nQ/QQxlxSSI/8QUFXKyIi+0ghTSSTbd+6CcO/imOTy2j3/losBH5MFtY3Bs1CkDqD6TefEmy9IiJSbRTSRDJNxWCGYysTMKkU5sbxMFjPKPSNQpPQ12POACu6KsCiRUSkuimkiWSa6Y+mblfEYWIpLIjjUZjTu4B5vTpxwX6TwZyvOl9K3pD7g61VRERqjEKaSCYpHoEvK8cmlsGHcYjBrOO6cGTvz5kZO4q/5FzDoZf8jR75ueQFXauIiNQohTSRDPH6g89x8p33EvokjucYnBhjXVEzjsxeQcTgMnuTy29+PugyRUSkliikiQSleARMH8UHXxZyxLgJnPJpHG9s+IAYFEaxmNHMt/Jh/hAKPnsO09pmIiINikKaSBDc8TEPY29t5ajPx+NNDU6PYcdE8Yixrumh5G5ajB3QJb35ufbWFBFpaAIJaWbWHHgG6AAsAS5097W7OG8JsBFIAPFdLfQmUicUj4Dpj+KeZOZH3Tj4jTnkfrEJb2ZMHXg0Rx+9hFgkAUDcssj96YyACxYRkaAF1ZJ2E/CGu99hZjelH9+4m3NPdPevaq80kRowbRTMLcMmldJj5WS8eQjOziZxZIyrw79nY2mc/ge35Ikre5EVdK0iIpIRggppZwMnpO8/Dkxg9yFNpG6osL4ZRVcC4FMfwT+IUz4pQWxNGbQK4d/J4W+dT+firP8w+4Dv8NjgntwzfhHXDygItn4REckogezdaWbr3L1Zhcdr3T13F+d9AqwFHHjI3R+uyvW1d6cE4tfNwVNdliRC+Kxt2ORSWOf4ASHoH8O6REgS4oYu/+HuIUcHW6+IiGSEWt+708zGA6138dTNe3GZvu6+3MzygNfNbIG7T9zN+10NXA3Qvn37va5XZJ8VXgHvPIK/V862d0LkrN8GB4Xg9GysIAJmAISKruTuwQpoIiJSuRoLae4+YHfPmdmXZnagu68wswOBlbu5xvL07UozexHoCewypKVb2R6GVEvavtYvslc2b+bNF8McOTZC3uaNfND2cArPWELokDCO7TjNcprB4LuCq1NEROqMoMakvQxcDtyRvn1p5xPMrDEQcveN6funAr+p1SpFdqV4BJSMhlYF+GcL+WBaWwqmLuKkLWUs79CKi876PfM6dee30cc5o/wVVncepu2bRERkrwUV0u4AnjWzK4FPgQsAzOwgYKS7DwIOAF60VBdRBBjj7q8EVK/I10pGw5Y4/ux72LtldN22juShEejfiLy25cyKd2XM8F70yD8dQNs3iYjItxJISHP31cDJuzi+HBiUvv8x0K2WSxOp3KpVbJ7ZkUavzMTKIHlYFsV9+7HmwFyGhd+gOHoaY4YfS4/8b8yDERER2SvacUCkKv72PyQfeJzEjDg55Un+1bkfTQY0pn+LqWzy5mQNuotIr/acHXSdIiJSbyikiezGlruLyFm6AN4uhZnlhJJgR2UxoOf9fNSiHR9lX0KYJENDb0AvzSgWEZHqpZAmsrPiEfhrI8mZtBWbXZ461i0LjotBbohX+QWv5QxkdYdh5C0ck1p6Q0REpJoppIlUtHAhfvP92PtleAj8mCzoG8OahXacEiHJoNJXYMgaQLM2RUSkZiikScOU3vA8tZkFlK80IhM3w9w4RMB7RqFPFNsvxPZNORywvC6wapFaz0REpMYppEnDVDIacFiRgImlZC2IQxS8bxR6R7EmqXCWIET4gMNg1SKs8AotRCsiIrVGIU0anj/k45+VwsQy7MM4HgP6R6FXjEROiHhWU2LxDWwKNWXRdz/QchoiIhIIhTSp/9I7BCxqdwG3vtmRMe98jn2SwHMMTopBYRTLSW3dFMGJ3PIZAE2BHgGWLSIiDZtCmtR/00fBR2V0GvUoT3+aINk4BKfEoEcUYoaRHm8GkNcl2FpFRETSFNKk/vrXDfiYB1Pdmp8nsKZG4vQcnu56GkNyJmChMOFkaokNK7pK481ERCSjhPZ8ikgdUDwCft08dZtM8tGDj+Pfuxd7eiu2OQlnZOM/asKTPQYxNGcCEZI7AhoWVkATEZGMo5Am9UPJaEjESY5+iNLWjTjk+98lWQbJs7LxHzYheUyMF7IH8mLrH7O687BUMMvrkrrVchoiIpKB1N0pdd7Tb39Ex5mH0GvyTEJrkkRbheDcHEKHR0iEwoxNnkyX4Y9wXn4u5wFwHFqEVkREMp1CmtQN6RmapNcq23J3ETlfLaT0/TBDJm3A1iVZfcD+5F5QhnWO4JaerWlJhkXexLSMhoiI1DEKaVI3lIwGT+Alo/nP7M85cdIc7J1Ssjc4tAmROL0xuZ2ShEJZQGoR2tWdU3trmrozRUSkDlJIk7qh8AqSb49i5rS2nDj1WWyz4+3DlJ/ZiNDBIZ5KDuBSG8/2bZ7CRcPJG3wX6tYUEZG6SiFNMtaYdz+F4hEMKXsVm1aGTS2jcOs6vGMY+sewDhFCHqJr8mkKDmzKwBatyVs4ZkeXqIiISF2mkCaZ6YHeDF0yD94txaaVwTagUwT6RbF2EdxTbWav5Qzkb0N7pbdu6otazkREpL5QSJOMsHLstbRa8DdwKNsSIjplCza9DMqAzhG8Xww7KIw74LA1t4BGP57OoKALFxERqSEKaZIRWix4CtuYhHfKiJaUYXHgiAj0i+F5YcxSLWdvNDmT3Avv06bnIiJS7ymkSeD+8dIUji82mr23CZLgR2XBcVFoGQZgobelwJYTKhrOAI01ExGRBkIhTWrHTuucrRx7LS2nPI5NKuXs2en9M7tn4X1jWG5oR8uZ3baezsFWLiIiEgiFNKkd6XXOmD6ST2Yto8OzL2IflKc2JuuRhfeJ4ft/Hc4ADAuyYhERkUAppEnNKx4BnsS/TGCTSukw91nIAnpFoU8UmoYw0uHMwlirAli1SHtqiohIg6aQJjXngd74yvmwPBXObEEcokDfKIt7deDQJl+wvbHMACu6SuubiYiIpCmkSY1YOfZaWs34AJtUBh/GIRu8fxR6xSjNyeJgW/l1Z6aF4dY1QZYrIiKScRTSpFrNWLKGx373GPeOfxD7JAGNDE6KQc8oRA0zyCJJuGg4TH8UcHVrioiI7IJCmlQPdxY9+QJdbr6U+z7bijc2OCUGhalw5ukFaGPrFqc2Ph98l7o2RUREKqGQJvtkxpI1HPiLozlw4jIKlifx/QwGZkP3LDwr1aFpAGY0+vF0APKCK1dERKTOUEiTb+efN+BPPszBE7PI/WIdNDM4IxvrlgUR2zFb843GZ3LylnGYujRFRET2ikKaVMmWu4vIWbeIRBLCc8phcim2KkmzFiH87Gw4MguLpFrOnPRszZxmDPjZk4HWLSIiUlcppEnl0jsF5MTj2AflhCeVYWuS0CoE5+VAlwgWTo052xHO8rrAD6YGXbmIiEidppAmu/dAb1g+D2aVw9ulsM6hdQi/IAc6R8CMUiJEPMnH+RdSMPyhoCsWERGpNxTS5L8Vj8BLRvNC+YmcWzILe6cMNjq0CcPAKNYpQpzUxucRS5JtDreupSDgskVEROobhTRJdWlOHwkYlDlWUsq577yIbXa8fRg7J4Z1DOPp5WefYwA98ptT8NlzWuNMRESkhiikNVDbJwJ8Yu3p4J8S2uYwrRSfWoZtdegYJtk/RqhDBE8PNnuj8ZnkXngfQ/Nz01dR96aIiEhNUUhriIpHkLNuEQZ03LwUppbCtDIoBTpFoH8UaxvB3Yi77RhvNiDoukVERBoQhbSGqGQ0tjkJU8qw6WVQBnSO4P1iPJF3OpeGXwc3/hkdSNtLHqDHjpYzERERqS0KaQ3EmHc/5bf/msvLm67nkCmbYUYZxIEjItAvBnlhSrP248hLH+G74xdx/YACzlY4ExERCYy5e9A1VLvCwkIvKSkJuozgpNc2W3nYxWz5cCL5iaWwNgnvlMJ75VgSvGsWdlwUbxHGUvs2wW3rAi5cRESk4TGzGe5euPNxtaTVRyWjwRO0WPAUrVbHsbdLYXZ56rnuWXjfGNY8BKT31QRSS9GKiIhIplBIq2+KR+CegFUJmLwNPijDQ0CPLOgTg/1DJDC+XlAjreiqgAoWERGRXVFIqw/SC9CWuRH9ohSbXApz44SywHpH4dgoNA2RdKMbz/DBr08LumIRERHZA4W0+mD6SGx5gujEUmxhHKLgfaPYsVFoHKpwotN6v1hgZYqIiEjVKaTVYT8e+x6nvnIDA9/egi2OY9ngx8egZxRyDIwdC9ECvNnkTO44v1ugNYuIiEjVKKTVRe7w1ltccNP19F36PjQy/KQY1jMKUaNj6RgeyHmKQaWvYEVXwOC7ALQYrYiISB2ikJbJikfA9EcBxwFPwuaPIjSZtB77LEGfJoafEoMeqXDmQCkR7mn6JANLX03tq5kOaCIiIlK3hPZ8SvUzswvMbK6ZJc3sG+uCVDjvdDNbaGaLzeym2qwxI5SMBhzcsYXlhB7dRNMxa7D1SRiYjV3XBOsTw2Kp/kwziJDkjPJXME+kXy8iIiJ1UVAtaXOAc6lkh24zCwN/AU4BlgHTzexld59XOyUGY8bStRw86giasTnVrTk/DpNKsS+TeDODM7KhWxaE0wPN0mPOEmbgxouhUzn+sFbkLRyTakkTERGROimQkObu8wHMrLLTegKL3f3j9LljgbOBeh3S7hm/iMeTm7C5qXDGV0loEYJzsrGjsiBUYflZh+3/Ew5v9yoTP/yK/p1acsGQXsD9wfwCIiIiUi0yeUxaG+CzCo+XAb0CqqV2lJfzmzXTST6whfCaBJ4XgvNyoEvk63C208YADlheF64fUACw41ZERETqthoLaWY2Hmi9i6dudveXqnKJXRzb7d5FZnY1cDVA+/btq1RjxigthdGj4Y476LB0KRx9NIz8JXb22RCq2rDBHsATV9bvDCsiItKQ1FhIc/d9XfFhGdCuwuO2wPJK3u9h4GFIbbC+j+9dO7ZsgUcegTvvhOXLoXdveOABGDjw635MERERaZAyubtzOtDJzDoCnwNDgIuDLamabNoEf/0r/PGPsHIl9O8Pjz8OJ5+scCYiIiJAcEtwfMfMlgHHAsVm9mr6+EFmNg7A3ePAtcCrwHzgWXefG0S91Wb9evjtbyE/H37+c+jWDd56K/UzYIACmoiIiOwQ1OzOF4EXd3F8OTCowuNxwLhaLK1mrF4N99wD996bCmpnnAG33AK9NIZMREREdi2TuzvrvpUr4a67UuPMNm2Cc89NhbOjjw66MhEREclwCmk14fPPU+PNHnooNXPzoovg5pvhiCOCrkxERETqCIW06rR0KfzhD/Doo5BIwKWXwi9+AQVau0xERET2jkJadVi8GH7/e3jiidTg/yuugJtugo4dg65MRERE6iiFtH0xfz787ncwZgxEo/D976dmbbZtG3RlIiIiUscppH0b77+fWkrj+echJwduuAFGjIDWu9pgQURERGTvKaTtrWQSzj8fvvgiNd7sJz+Bli2DrkpERETqGYW0vRUKwTPPQIcOkJsbdDUiIiJSTymkfRta50xERERqWCDbQomIiIhI5RTSRERERDKQQpqIiIhIBlJIExEREclACmkiIiIiGUghTURERCQDKaSJiIiIZCCFNBEREZEMpJAmIiIikoEU0kREREQykLl70DVUOzNbBSwNuo590BL4KugipFros6xf9HnWH/os64/68Fnmu3urnQ/Wy5BW15lZibsXBl2H7Dt9lvWLPs/6Q59l/VGfP0t1d4qIiIhkIIU0ERERkQykkJaZHg66AKk2+izrF32e9Yc+y/qj3n6WGpMmIiIikoHUkiYiIiKSgRTSMoCZXWBmc80saWa7naFiZqeb2UIzW2xmN9VmjVI1ZtbczF43sw/Tt7m7OW+JmX1gZrPMrKS265Td29P3zFLuTT//vpkdE0SdsmdV+CxPMLP16e/hLDP7VRB1yp6Z2SgzW2lmc3bzfL38XiqkZYY5wLnAxN2dYGZh4C/AQOBwYKiZHV475cleuAl4w907AW+kH+/Oie7evb5OHa+Lqvg9Gwh0Sv9cDfy1VouUKtmLv5mT0t/D7u7+m1otUvbGY8DplTxfL7+XCmkZwN3nu/vCPZzWE1js7h+7exkwFji75quTvXQ28Hj6/uPAOcGVIt9CVb5nZwNPeMpUoJmZHVjbhcoe6W9mPeLuE4E1lZxSL7+XCml1RxvgswqPl6WPSWY5wN1XAKRv83ZzngOvmdkMM7u61qqTPanK90zfxbqhqp/TsWY228z+bWZH1E5pUgPq5fcyEnQBDYWZjQda7+Kpm939papcYhfHNDU3AJV9lntxmb7uvtzM8oDXzWxB+r8UJVhV+Z7pu1g3VOVzmklqO55NZjYI+Aep7jKpe+rl91IhrZa4+4B9vMQyoF2Fx22B5ft4TfkWKvsszexLMzvQ3Vekm9pX7uYay9O3K83sRVJdMwppwavK90zfxbphj5+Tu2+ocH+cmT1gZi3dva7vA9kQ1cvvpbo7647pQCcz62hmUWAI8HLANck3vQxcnr5/OfCNVlIza2xmTbffB04lNXlEgleV79nLwGXp2WS9gfXbu7glo+zxszSz1mZm6fs9Sf2buLrWK5XqUC+/l2pJywBm9h3gPqAVUGxms9z9NDM7CBjp7oPcPW5m1wKvAmFglLvPDbBs2bU7gGfN7ErgU+ACgIqfJXAA8GL634YIMMbdXwmoXqlgd98zM7sm/fyDwDhgELAY2AJcEVS9sntV/CzPB75vZnFgKzDEtcJ7RjKzp4ETgJZmtgy4FciC+v291I4DIiIiIhlI3Z0iIiIiGUghTURERCQDKaSJiIiIZCCFNBEREZEMpJAmIiIikoEU0kQk45lZCzOblf75wsw+r/A4Ws3v1czMfrCb5yaY2Wk7HfuxmT1QyfUmmFlhddYoIg2DQpqIZDx3X+3u3d29O/Ag8Oftj9ObZ++SmX2btSCbAbsMacDTpBZFrWhI+riISLVSSBOROsnMvmdm09ObY//dzBqljz9mZn8ys/8AfzCzQ8xsavrc35jZpgrX+Fn6+Ptm9uv04TuAQ9KtdP+309s+D5xhZrH06zsABwGTzeyvZlZiZnMrXGvnmiu+9/lm9lj6fqv07zA9/dM3ffz4Ci2G723fqUJEGgaFNBGpq15w9yJ37wbMB66s8FwBMMDdRwD3APe4exEV9vIzs1NJbabdE+gO9DCz/sBNwEfpVrqfVXxDd18NTANOTx8aAjyTXqX+ZncvBLoCx5tZ1734Xe4h1TpYBJwHjEwf/ynww3QLYj9Sq+KLSAOhkCYiddWRZjbJzD4AhgFHVHjuOXdPpO8fCzyXvj+mwjmnpn/eA2YCnUmFtj2p2OVZsavzQjObmb7eEcDhe/G7DADuN7NZpPYg3C/davY28Cczuw5o5u7xvbimiNRx2rtTROqqx4Bz3H22mX2X1L5+222uwusN+L27P/RfB1NdmJX5B6ngdAyQ4+4zzawjqVavIndfm+7GzN7Fayvuw1fx+RBwrLvv3FJ2h5kVk9qTcKqZDXD3BXuoT0TqCbWkiUhd1RRYYWZZpFrSdmcqqS5E+O9B/68Cw82sCYCZtTGzPGBj+tq75O6bgAnAKL5uRduPVDBcb2YHAAN38/IvzayLmYWA71Q4/hpw7fYHZtY9fXuIu3/g7n8ASki19olIA6GQJiJ11S+Bd4HXgcpal34M3GBm04ADgfUA7v4aqe7PKeku0+eBpulxZ2+b2ZxdTBzY7mmgGzA2fa3ZpLo555IKb2/v5nU3Af8C3gRWVDh+HVCYnsAwD7hme+3pOmaTGo/270p+TxGpZyw13lVEpH5Kz/rc6u5uZkOAoe5+dtB1iYjsicakiUh914PUoHwD1gHDgy1HRKRq1JImIiIikoE0Jk1EREQkAymkiYiIiGQghTQRERGRDKSQJiIiIpKBFNJEREREMpBCmoiIiEgG+v+ZHB3a72q90AAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.figure(figsize = (10, 5))\n", "plt.scatter(ytr[:, 0], trPred, label = \"Train\", s = 2)\n", "plt.scatter(yte[:, 0], tePred, label = \"Test\", s = 4)\n", "plt.plot([-1.1, 1.1], [-1.1, 1.1], label = \"x = y\", c = \"r\")\n", "plt.xlabel(\"Target Values\")\n", "plt.ylabel(\"Predicted Values\")\n", "plt.legend()" ] }, { "cell_type": "code", "execution_count": 29, "id": "66b25043", "metadata": {}, "outputs": [], "source": [ "allpred = DTR.predict(X)" ] }, { "cell_type": "code", "execution_count": 30, "id": "2f157846", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 30, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnAAAAE9CAYAAACLPV+MAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABMo0lEQVR4nO3deXiU5fn+//c1kwU0irKoiGzyEayNCDSA/lDADVxQAraVRUSRIooWaatlMRKMLF/UpmpxiUCLFkxtgUgBK2oLWqxIBMWogAuyWKosigZDkpm5f3/MZEwgIEuSJzM5X8cxR2aeeZ6ZK8OSM/dqzjlEREREJHb4vC5ARERERI6MApyIiIhIjFGAExEREYkxCnAiIiIiMUYBTkRERCTGKMCJiIiIxJgErwuoSY0bN3atWrXyugwRERGRH/T222/vdM41qey5OhXgWrVqRX5+vtdliIiIiPwgM9t8sOfUhSoiIiISYxTgRERERGKMApyIiIhIjKlTY+BERERiWWlpKdu2bWPfvn1elyJVqF69epxxxhkkJiYe9jUKcCIiIjFi27ZtnHDCCbRq1Qoz87ocqQLOOXbt2sW2bdto3br1YV+nLlQREZEYsW/fPho1aqTwFkfMjEaNGh1xq6oCnIiISAxReIs/R/NnqgAnIiIiR2ThwoWYGevXr48e++yzz0hNTQVg+fLl9OnT54Drli9fToMGDejYsSPt2rWje/fuLF68+Affb/ny5bzxxhtV9w3EAQU4EREROSLPPfccF154Ibm5uUd87UUXXcTatWvZsGEDjz76KHfccQevvvrqIa9RgDuQApxIFcjIK6DV2CW0HruEjLwCMvIKaDNuKRl5BV6XJiJSpQoLC1m5ciWzZs06qgBXXocOHbjvvvv4wx/+AMDf//53unbtSseOHbnsssv44osv+Oyzz3jyySfJzs6mQ4cOvP7665WeV9d4OgvVzGYDfYAvnXOplTxvwCPAVcB3wE3OuTWR566IPOcHZjrnptVY4VLn9MpewcYvCkn0G6VBB1DhfhkHPPvm9zuf/PnNzcxbtYVBXVuw6N3P2VMUINFvhEIwqGsLstIP+GsvIlKr5eXlccUVV9C2bVsaNmzImjVr6NSp01G/XqdOnXjwwQcBuPDCC3nzzTcxM2bOnMn06dN5+OGHGTlyJCkpKfzmN78B4Kuvvqr0vLrE62VE/gT8AXjmIM9fCZwVuXUFngC6mpkfmAFcDmwDVpvZIufcB9VescS9jLwC5q3aQko9P3uKAhWeKx/Y9g9vlXFA0LkKoa7sumff3Myzb27GgBvOb6kwJyLVouz/tKr6pfG5557jrrvuAmDAgAE899xzxxTgnPv+/9Jt27Zx/fXXs337dkpKSg66rMbhnhfPPA1wzrnXzKzVIU7pCzzjwn+6b5rZSWbWFGgFfOyc+xTAzHIj5yrAyTFZMeN2Mr+cy/1JEAqBL/nYXs8BVu5rZc+HAP9aCK4Fn1V+HmdeAjcuPLZiRKROmrdqC0HnmLdqyzEHuF27dvHPf/6TgoICzIxgMIiZMX369KN+zbVr1/KjH/0IgDvvvJNf/epXXHvttSxfvpzMzMxKrznc8+KZ1y1wP6QZsLXc422RY5Ud71qDdUmMy8griLaKtT01BYCNXxSyMTkXfyRB+avgfWy/r5U9XzYQ9ZDv9+k/q6AaEamLBnVtEW2BO1Z/+9vfuPHGG3nqqaeix3r06MG///1vmjdvfsSvt27dOrKyspg5cyYAe/bsoVmzZgDMmTMnet4JJ5zAN998E318sPPqkto+ieFgjRYHO37gC5iNMLN8M8vfsWNHlRYnseesCUtpNXZJhS7NjV8UsvGLQgBmBq4k6MA5ol+P5Rba72tlzwcqeb/9j/2b87z6yEQkxmWlp/LJ1KuqrPu0X79+FY5dd911zJs377Bf4/XXX48uIzJq1CgeffRRLr30UgAyMzP52c9+xkUXXUTjxo2j11xzzTUsXLgwOonhYOfVJVa+79mTAsJdqIsPMonhKWC5c+65yOMNQE/CXaiZzrnekePjAJxzUw/1XmlpaS4/P78qy5fabutb8Gx/KPmWEGCRv+7lfwv4ziWQbEH8OIJAMUlklQ4hN3RppS+Z6Dc+mnwV8H1LXtk4ttzVW6Jj3Iac35JVm3ZFw+GxalA/gcJ9QU1+EKnDPvzww2h3o8SXyv5szext51xaZefX9gB3NXAH4VmoXYFHnXNdzCwB2AhcCnwOrAYGOefeP9R7KcDVLRl5BVz69m309K874mt3uxSyO75U5UGpssHEbcYtJXgU/w7bnprCsjE9qrQ+EandFODi15EGOK+XEXmOcItaYzPbBkwEEgGcc08CSwmHt48JLyNyc+S5gJndAbxEeOjQ7B8Kb1L3zFu1hffpT0f7iBOtqEKrW/n7pZaIjwB+5wgaJCQcR8MrJpOVVvWtXFnpqQeEwrLxKZXNeh3ge5VJiX8iieABr7Vvt5/+4zIoPKWTgpyISB3jeQtcTVILXHwrPzGhk23k3oRnacQeTvLtZWrpoEq7RGtjK1bZmnMAa5JH0NAO3gW7PNiem0rHMkTLkIjUCWqBi18x1QInUlUy8gpounoqnyQvjs7MKb838D2JfyG3uGKAq62hpyxQZuQVMP2t6w/aAgfQw7eOTcmDILIMyV+Tr2PA+Nk1Wa6IiHhAAU5iXlnL28bkF6NLgJQJOiikPtNLrwdqb2irTFZ6KqT/DvhdhdZFgE3JgyoEVAiPJehfnAcowImIxDsFOIlJ5016iT1FAW6ut5z73EzuTw4RIrwsR1muKcHHxNKbSewyjL+u2sKQ82N39mb5sXOtxy5hTfBMOvk/PeC8RIKUTGzI0ua/Jn34hJouU0REakhtXwdO5AAZeQXRwf53unkkWggz8BsE8NO6eB6ti+fRrvjPrGnSt0rXQKoNbji/JT8NPBD9PlsXh9dfMgvfkixI962Pk5FXQJtxS8nIK/C4YhGJJ36/nw4dOpCamsrPfvYzvvvuu6N+rZtuuom//e1vAAwfPpwPPjj4hkrLly/njTfeOOL3aNWqFTt37jzgfcsvRgzhPV6vuuqqw6q1NlCAk5iSPW08E9Z2Z1PyIDYlD+JkCissfjszcCVDzm/JZ9Ou5rNpV9e6CQpVISs9lU2R769sF4nlwXMrLBB8IoVMWtuNj5MGckn+SHplr/C4ahGJF/Xr1+edd96hoKCApKQknnzyyQrPB4OVj9n9ITNnzuScc8456PNHG+AqM3DgQHJzcyscy83NZeDAgVXy+jVBAU5iQkZeAf3HZXNn0QzqWTDa2lR2WxFqz/8Vz2N753Fx09J2OJaN6cFn067mn2lPcmbxPErxYwYJFtlX1aCn/z02flFIq7FLaD12iVrkRKTKXHTRRXz88ccsX76ciy++mEGDBnHuuecSDAa5++676dy5M+3bt4+2djnnuOOOOzjnnHO4+uqr+fLLL6Ov1bNnT8pWivjHP/5Bp06dOO+887j00kv57LPPePLJJ8nOzo7uxrBjxw6uu+46OnfuTOfOnVm5ciUQ3q+1V69edOzYkVtvvZXKVtu47LLLWL9+Pdu3bwfgu+++45VXXiE9PZ3777+fzp07k5qayogRIyq9vnyrXn5+Pj179gRg7969DBs2jM6dO9OxY0deeOEFAN5//326dOlChw4daN++PR999NExf/YKcFLrnTfpJZ59czO/TFhAQmSAW/nWpq+D9fhD8LqYmqBQ1cpa5f5zygACzihxvuj2XRCe9PBx8iDu9s9j3qot3hYrInEhEAjw4osvcu655wLw1ltvMXnyZD744ANmzZpFgwYNWL16NatXr+bpp59m06ZNLFy4kA0bNvDee+/x9NNPV9qitmPHDn7xi18wf/583n33Xf7617/SqlUrRo4cyZgxY3jnnXe46KKLGD16NGPGjGH16tXMnz+f4cOHAzBp0iQuvPBC1q5dy7XXXsuWLQf+n+f3++nfvz/PP/88AIsWLeLiiy/mhBNO4I477mD16tUUFBRQVFTE4sWLD/szmTx5MpdccgmrV6/mX//6F3fffTd79+7lySefZPTo0bzzzjvk5+dzxhlnHM1HXoEmMUit9uS9N7DG/3d8yeHHzkEQyAn0YXpwEFA713LzSo9Rj9Mr+3o2flHInxKnVdiFIgG4LWExI1hM4cQkXml+lyY6iNQFW9+C5dOg51ho3uWYX66oqIgOHToA4Ra4W265hTfeeIMuXbrQunVrAJYtW8a6deuiY8b27NnDRx99xGuvvcbAgQPx+/2cfvrpXHLJJQe8/ptvvkn37t2jr9WwYcNK63jllVcqjJn75ptv+Pbbb3nttddYsGABAFdffTUnn3xypdcPHDiQu+++m9GjR5Obm8uNN94IwL/+9S+mT5/Od999x+7du/nxj3/MNddcc1ifzbJly1i0aBEPPfQQAPv27WPLli1ccMEFTJ48mW3bttG/f3/OOuusw3q9Q1GAk1qpbNmMj5KXHLA0iM/Bw6HBfDbt4INN67KyMJsz91vO3fDr6ELAZd3NCUAKJZGJDn3rbKulSJ2xfBp88mr4/pAFx/xyZWPg9nf88cdH7zvneOyxx+jdu3eFc5YuXYrtvwbSfpxzP3gOQCgU4j//+Q/169c/4LnDub5bt25s376dd999lzfeeIPc3Fz27dvH7bffTn5+Ps2bNyczM5N9+/YdcG1CQgKhUAigwvPOOebPn0+7du0qnP+jH/2Irl27smTJEnr37s3MmTMrDa9HQl2oUrvM/wWhzAbcv7Ybm5IHkUD4H0j5LtOFwW4M6trC40JrvxGDB9Do/s+5r+NKFgS7RT+/sq7VBhQyce2FPD5hsMbFicSznmOhzaXhrzWkd+/ePPHEE5SWlgKwceNG9u7dS/fu3cnNzSUYDLJ9+3b+9a9/HXDtBRdcwIoVK9i0aRMAu3fvBuCEE07g22+/jZ7Xq1cv/vCHP0Qfl4XK7t27M3fuXABefPFFvvrqq0prNDN+/vOfM3ToUK666irq1asXDWONGzemsLDwoLNOW7Vqxdtvvw3A/PnzK3zfjz32WHTc3Nq1awH49NNPOfPMM/nlL3/Jtddey7p1R75H9/4U4KTWyMgrILTueXxUnKDgHMy2vvxfyXPc13El1z2wVK1GRyArPZXrHlhK75MW0bp4HgbRZVcSzDE84UWNixOJZ827hFveqqD79HANHz6cc845h06dOpGamsqtt95KIBCgX79+nHXWWZx77rncdttt9Ohx4PCXJk2akJOTQ//+/TnvvPO4/vrwQuzXXHMNCxcujE5iePTRR8nPz6d9+/acc8450dmwEydO5LXXXqNTp04sW7aMFi0O/gv/wIEDeffddxkwYAAAJ510Er/4xS8499xzSU9Pp3PnzpVeN3HiREaPHs1FF12E3++PHs/IyKC0tJT27duTmppKRkYGAH/5y19ITU2lQ4cOrF+/Ptpdeyy0F6rUGk/cewO3+v9O+Ybvslbw3S6FhpM+96SueJKRV8Dg/J/Rzv9fQoR/gwsBTwev4cHgIAZ1jd3FjkXqAu2FGr+OdC9UtcCJ5zLyCnh8wmBG+v8eXfrCAdeVZLIh2Iw9rj6vNb/d6zLjQlZ6Kr9s9CSti+exh5RoS9ww/1KCzqklTkQkRmgSg3juz29uZmPy0mh3qSM8zm2Na8uff/I8WemppHtdZBwpm+Qwdvz1TEx8hiRKSSTIp8mDCGIEJjpKSeSl5mM0S1VEpJZSC5x4JiOvgFZjlzA7YWp0soIZvBZqz68Do+r0um41IbHLMFJLnomOOfQZJJojwaC+lUa34xIRkdpHAU488+c3NwPhnQLKxrqVOh+PBvorvNWAsj1i19n/VZidCuGv00uv59nIn5GI1B51aex6XXE0f6YKcOKJ7GnjeS/5Jj5ODi/G6xyUOB8ZpTfz466XK7zVoPMy36b3SYs4s3geG4Knh8McMDVxFn9MmKpWOJFapF69euzatUshLo4459i1axf16tU7ous0Bk5qVEZeAX9+czNvJ88hxUqix52DtsV/Du+qoPBW48rGxbUaC2v8I6KL//b0v8eHq6fSZtVgzVAVqQXOOOMMtm3bxo4dO7wuRapQvXr1jnh7LQU4qREZeQXMW7WF6f4/MCl5JXtdEoHIc35gpZ3HZ9Ou9rJEAYac35Lpb13P5MRZ+CMzgm9LWMwtvMik1TdD+sNelyhSpyUmJka3mJK6TV2oUiPeX/Uy/0y4k/7+lfgMTvCVkGCwMtQem7SHCzNf87pEITwubtqU3/GrlOlsDjbGuXCIS7YgkxNmkj1tvNcliogICnBSA3Lm5vKnpP9HS//O6FIhO0PHszbUhg1nj/K6PKnEI3ffysyfvBDdggvCQe6XRTPCG2OLiIinPA1wZnaFmW0ws4/N7IBN2szsbjN7J3IrMLOgmTWMPPeZmb0XeU7bK9RibdfP4EQrAsIzHZ8I9GHQyc/R8f41jBg8wOPq5GCy0lN5qtFvK4Q4v4GbeTnz771KkxtERDzk2VZaZuYHNgKXA9uA1cBA59wHBzn/GmCMc+6SyOPPgDTn3M7DfU9tpVWzcubm0nv9BJr5drEjdCLbaUynEU/W6H58cuxajV3CHxOm0tP/HvD9/rQA6+z/OC/zbQ+rExGJX7V1K60uwMfOuU+dcyVALtD3EOcPBJ6rkcqkSrRbP4OW/p0kmOM03x7W+lMV3mLQkPNbcnNgHKX4o+HNIhMc2ruP1RInIuIBLwNcM2BrucfbIscOYGbHAVcA88sddsAyM3vbzEZUW5Vy1Ha6EyqMnxriFntbkByVrPRUhpzfkpmBK6PhDcJBbk3wTO2fKiLiAS8DnFVy7GD9udcAK51zu8sd6+ac6wRcCYwys+6VvonZCDPLN7N8rZtT/UY/+BRrMjqx675m9PevjLbYBJ2xIDnd6/LkKGWlp7K98zjGld5CoUui0CXzv2ADOvk/5d3EofQfl62WOBGRGuTlGLgLgEznXO/I43EAzrmplZy7EPirc27eQV4rEyh0zj10qPfUGLjqlT1tPHcWzSChXDR3Dv6SfB0Dxs/2rjCpUhl5BTz75mY2JQ+q0Bq3PHguFz/wb2+LExGJI7V1DNxq4Cwza21mScAAYNH+J5lZA6AH8EK5Y8eb2Qll94FegH7999jQojkHhLe8YDeFtzhT1qVaGEqKHjML79qQN3Oyh5WJiNQdngU451wAuAN4CfgQeN45976ZjTSzkeVO7Qcsc87tLXfsVODfZvYu8BawxDn3j5qqXSrKyCtg7PhfcRxF0TFvIQdPBa+h3wNLvS1OqkVWeirD3IToYr8QDnHdtz6urlQRkRrgWReqF9SFWvVGP/gUQ795inN9m0i0EAClzsdvUqbxyN23elydVLde2StI+XIN2Ql/4AzfTnyEu1L/mfak9k0VETlGtbULVWLcu5k/4feF99DJ/wmJFiIU+V3g01BThbc6YtmYHvy46+Vs4vTo3qk9/e/x7Jub1RInIlKNFODkqOTMzaW9+7jCIPbnk69jRbA9K86+z9vipEZlpaey8MQb2Bk8HufCU8nXJw9h+Jq+bJ1+gbbeEhGpBupClaOyPaMFTf17gO8Xdl0RbE+PrNc9rky80mrsEjYmDyHJghWfaHMpDFngTVEiIjFMXahSpTLyCjjVtyf6eC/JLA+218b0dVz5xX7LfOPqk+P7uXdFiYjEKbXAyREpm7RwJp/TwLePUnzcV3oziV2GadC6ADAp827ucs/yv1BDxgV+wRrXFggHPP0dERE5fGqBkyqRkVfAr76ZTif/J5zk34cZvBFK5cXk3vrBLFETMx+kwaT/cWfDJ1jj2jLA9yprkkdQ+pbWAxQRqSoKcHLY5q3awum+73czK3U+Fp54A+9O7O1hVVJbnbhjLe8m3cLUxFk0tEKmJs5iVuaNXpclIhIXErwuQGq/WZk3MtQt4v0kH3vccZxMIcUksqz5GB4ZruVCpHJTGy2lwbdF0cdmcJN74RBXiIjI4VILnBxSztxcbnYvkGCOehaksa+QbaHGHDdpJ+nDJ3hdntRiZ/38AfaE6oeXFokMtfUBLz/Q19O6RETigQKcHFLa+un4yq31BlToRhU5qOZdeKjTK7QunscTgT7R5WYuK11O/3HZ9Mpe4XWFIiIxSwFODul49gGw1yWxINiNEudnfnI/j6uSWJGVnspn065me+dxLAh2i4a4+UmZzPlqiBb5FRE5SlpGRCqVO2UY/Yvz2O1O4DTf16wJtqFT1hqvy5IY1it7BZ12vMDUxFnRHTyWB9vz6k+e0CxmEZFKaBkROWI/LV5AkgVpbF+zItie/LPv8bokiXHLxvTglkYFFbZfO9/3vpYXERE5CgpwcoCF916Fn3DLbAA/PbJeZ8TgAR5XJfFg3K6r2BxsHO1KrWdB7kn8iza+FxE5Qgpw8r38P1E86VTS/SujrSRfuJO9rUniyo+6XMYlgcd4ItCHgDOKnY8kStUKJyJyhBTgBAjvsrD77xNIduEdFsqGRu5yDbwtTOJKVnoqn0y9iu2dx9GuZB57OY4UK1YrnIjIEVKAEwDeX/UyO1wDQpHgZga7XYrGvkm1KAty00uvZ7dL4V/B8xiztjd5Myd7XZqISEzQTgwCQE7iQzT2FVIYSiKIn/+FGtJuxJ8Y0byL16VJHDv/zEacsHUv/SPd9t23Pg5ogWgRkR+iFjghd8owGlkhAMdbCXMDl/LnnzwPCm9SzdJ3zSTRXLTb/ni+o3DiKWqJExH5AQpwdVjO3FxWZFzET4sXRCctmMHwhBe1LpfUjEsnRsdbmkGyhUixYrpvfVxj4kREDkEBrg7rveFeevjXRZcMKXVGqfOxIDnd28Kk7ki7CbvmEUrd9/8VOQcnU8iktd3InTLMw+JERGovBbg6KmduLs1tB0C0++r6konc3/F1BozXkg5Sg9Ju4o9tn2CPqw+E/z6agc+gf3Get7WJiNRSngY4M7vCzDaY2cdmNraS53ua2R4zeydyu+9wr5WD65W9IrpJvXPh24JgN37c9XJ1nYonRgwewF/a/p6vQ/WiXarOwczAld4WJiJSS3kW4MzMD8wArgTOAQaa2TmVnPq6c65D5Hb/EV4r+8mbOZmlX19LR98nQLilY0WoPe+kTVd4E0+NGDyAk+7/gnGlt7DH1WdjqBmvhNLolb3C69JERGodL1vgugAfO+c+dc6VALlA3xq4tk7rvvVxEiJdVABFLpGNZ49SeJNaI7HLMD4JnU47/+fMT8rkhd39KJ14ssbDiYiU42WAawZsLfd4W+TY/i4ws3fN7EUz+/ERXivljH7wKXaEGkS7qABKSdA+p1KrlP9lwgzq+4MkWkjj4UREyvEywFklx9x+j9cALZ1z5wGPAXlHcG34RLMRZpZvZvk7duw42lpj3ugHn2Jy4b20838OhMcX7XN+/tn8To8rEzlQ/tn3RMfDFQd9hBzsCqUw+sGnvC5NRKRW8DLAbQOal3t8BvDf8ic4575xzhVG7i8FEs2s8eFcW+41cpxzac65tCZNmlRl/TEjHN4ySLES4Pvu08kdXyN9uFa9l9pnxOABPNzpVc4snsd/SMVn0NS/h+mF48iZm+t1eSIinvMywK0GzjKz1maWBAwAFpU/wcxOMwvHDTPrQrjeXYdzrYTlzM0lq3AiKVZc4fgud7zGvUmtlpWeyqZpV9NzxMPRbv9kC9Ju/QxvCxMRqQU8C3DOuQBwB/AS8CHwvHPufTMbaWYjI6f9FCgws3eBR4EBLqzSa2v+u6j9eqy/nxOtKPrYOfgvp9D4F3neFSVyJJp34X+hBgAEHbwY7EybcUu1U4OI1GnmXKVDx+JSWlqay8/P97qMGpMzN5dbNo7Eb+E/45CDJwN9uH3yXI8rEzkyOXNzGbTxTlKshM3BxmzidB4N9GfB1DFelyYiUm3M7G3nXFplz2knhjiVN3MywzbeWiG8jS+9he2dx3lcmciRGzF4AEH8ALTw7aSnfx1TE5/2uCoREe8owMWpsvXeINxt6jO40r9a494kZv2r+Z049/0knJb2P3WjikidpQAXh2Zl3kgDCqMDv//nGrAi2J4NZ4/ytjCRY5A+fAKvnTL4+wkNBJm49kJWzLjd28JERDygABeHhrpF+CO7LYQc/L3ddHpkva4FeyXm9Rj1eDTEmUGCOS74UsuKiEjdowAXZ7KnjcdXbk3jfSQquElc6THq8eh958DhKJzYhLyZkz2sSkSkZinAxZlbi57CZ+EfbN+5RJ6qP8LrkkSq3Eo7D+egBB/JFiLFSui7dbq6U0WkzlCAiyPZ08ZTnwAQ7l5aHfoRY8ZO8bgqkap3YeZr2KQ9TCy9mYALz2owQ92pIlJnKMDFidwpw7iraEZ0ht53LkGTFiTuJXYZRk7gakIu3OrsJ6iuVBGpExTgYlzO3FzWZHTi58XzsUjX6YZgM261TI19k7iXlZ7KoMQV+CKTdvwGV2990OuyRESqXYLXBcixabd+Bp38n0Qfb3FNaJf1Ac96WJNITZpeej2TEmaRgAvPTKXu7C4jInWXWuBi2IoZt/P/+QrYHjqJDcFmrAm2YUypuk2lbvF3vpm2xXNZEzwz0gJ9OisyLiJnrsbDiUj8UoCLYRd8mUuihTjNvibXfzU/CzzAj7te7nVZIjUqKz2VTdOuZnLwRlaE2lPfSujhX8fgjXcqxIlI3FKAi1GzMm/ETzC6oOmdbh6fTL1KW2VJnfXLhAX09K/jdN9uAI63Etqtn+FxVSIi1UMBLkYNcYujuy2UOD9z6g/1uiQRT208exQrgu35d/DHOAcBBx/7W3tdlohItdAkhhiTPW08txblUC/S+vY/14BRpWNYMGmM16WJeCo863oAuyc2i0xmgH7uVa/LEhGpFmqBizFDi+ZwnJVGl00oIYm1rq3XZYnUGnPqD6XQJbPP+UmilOxp470uSUSkyinAxZCMvAKml16PK7dKwum2ixvOb+ldUSK1zJixU0iZ9CXfUZ8UK+YXRU+ze2IzBTkRiSsKcDHk2Tc3s9E1J0h4uwXn4I1TBmrigkgl5tQfym6Xgh9HQytkWNFszpv0ktdliYhUCQW4GNF/XDb/SbqN+UmZJJijxPm530bQY9TjXpcmUiuNGTuFhpM+Z0voFABOpIj/BAdpaRERiQuaxBADZmXeyPykF6L7nDoHf2r7OBO1VZbID/qva0hb9zlmcJwFIkuL6N+OiMQ2BbgYcKP7e4Xw9kpiT+1zKnKYevjfi+4TDPCNv4G3BYmIVAFPu1DN7Aoz22BmH5vZ2EqeH2xm6yK3N8zsvHLPfWZm75nZO2aWX7OV15yMvAJ2hk6ocOzye1/wqBqR2LPELiLkIER45nZv94bXJYmIHDPPApyZ+YEZwJXAOcBAMztnv9M2AT2cc+2BLCBnv+cvds51cM6lVXvBHumQfw+n+fYA4RaEvGA3jysSiS3XZC5mYseVPBXoQ6nzsdudwJqMThoLJyIxzcsWuC7Ax865T51zJUAu0Lf8Cc65N5xzX0UevgmcUcM1eiZ72nj2TWxIf//KaPfP022fot8DS70uTSTmZKWn8nBoMCtDqZzm+5pO/k8YvvFWcqcM87o0EZGj4mWAawZsLfd4W+TYwdwCvFjusQOWmdnbZjaiGurzTM7cXEYVPUk9C0bD22K7SOPeRI7BoK4teDTQn12hFJwDn8F1xQu9LktE5Kh4GeCskmOukmOY2cWEA9xvyx3u5pzrRLgLdpSZdT/ItSPMLN/M8nfs2HGsNdeIa9bfQ5KFt8pyDhYEu3FN5mKvyxKJaVnpqSyYOoaTrfD7SUE4daWKSEzyMsBtA5qXe3wG8N/9TzKz9sBMoK9zblfZcefcfyNfvwQWEu6SPYBzLsc5l+acS2vSpEkVll99To2MeTODr0jh14FRHlckEj9eTewZnZGaZC6yrIiISGzxMsCtBs4ys9ZmlkR4YaZF5U8wsxbAAmCIc25juePHm9kJZfeBXkBBjVVejfJmTqYEP85BsfMzvfR62p6a4nVZInHj8ntf4IlAH4LOCDjo7PuQvJmTvS5LROSIeLYOnHMuYGZ3AC8BfmC2c+59MxsZef5J4D6gEfC4hfs8ApEZp6cCCyPHEoB5zrl/ePBtVKmcubnctPVhkixIifMzsCSDBVPHeF2WSNzZ3nkce9Yup6EVkkApV299iJy5bTTOVERihjlX6bCzuJSWluby82vvknEbMs6hnf9znAt3ny4Ptqdn1utelyUSl8aO/xUPJM4mwcL/B24ONmbmT17Q3sIiUmuY2dsHWypNe6HWIqf5dgPh8LbbpbDwxBs8rkgkfiV2GcbPSyYSivwO28K3kw7592jDexGJCQpwtUDezMnsmXga34WSKHRJbAg1429tH+aRu2/1ujSRuJWVnsqPu17OiuC50Vbvfv6V7CkKeF2aiMgPUoDzWO6UYfTdOp0GVkRT/x52hU7kqtKHNBZHpAZkpadynn9TdL3FID5WJP6StdqpQURqOQU4D+XMzeXnxfMrbLTdzLeLQV1beFuYSB0yp/5QSpwfM0i0EC39O+no/4S09dO9Lk1E5KAU4DzUbv0MfJHwVraw6OeusQZRi9SgMWOncHfK1OhYuDo0r0tEYpgCnEcy8gp4Mdg5Gt6cgzXBNvzuhLu9Lk2kznnk7lt5/ZTBBJxRQgIbgs3IP/ser8sSETkoLSPigbyZk7l660MkEIqGt9eC59LjgX97XZpInbZ7YjMaWiF7XH3eCZ3FhrNHaTyqiHhGy4jUMt23Pk6ihaLdpgF8Cm8itcD00uvZ7VL4X6ghPfzraLt+Bhl5cbHJi4jEGQU4D0wvvb7CeJunA1d5W5CIAOG14TqXPM3jJ45mc6gJ3XwFNF09VSFORGoddaHWsOxp4xlZlEMipXzumjCmdBQ/7nq5Ji6I1DKlE08m0UI4B8uD53KxWslFpIapC7WWmJV5I3cVzaC+lZJg0NR2K7yJ1FL/dY2A8CSjnv73GP3gUx5XJCLyPQW4GjTELY6Oews5eNb6KLyJ1FK/O+FudoVSojPFHyocS97MyV6XJSICKMDVmCGZj7HLnUDIwfZgA2a2fYpbMp/xuiwROYhH7r6VgSfP4wt3EhBe5Pfyrb/3tCYRkTIKcDXg5Qf68oy7l6a+r/EZbKSlliYQiQHLxvSgkX0TfVyPEnWlikitoABXzWZl3shlpcujXaffuQQ2nD3K26JE5LA9Y9dEd2fwG/T75s/eFiQiggJctSs/7g3A/MlqfROJIbdkPsMTgT6UOh+bQ415NNDf65JERBTgqkvezMkUTWyMnyDOQYkzCl0S9a+a4nVpInKEtncex6LgBbSwnfwtKZMvMlvC1re8LktE6jAFuGrSfevj1LdS/BaewVYQOpP+J/0N0m7yujQROUJZ6an086/EDHwGp/I1e2ZeS87cXK9LE5E6SgGumsypP5QilxjdcQHCA6JFJDa9e3IvnCM6Hq6BFdFjw/3eFiUidZYCXDXImzmZ24ueIplSVgTPZXmwPfln3+N1WSJyDDre9Vfu67iSQLn/Ns+yz8meNt7DqkSkrlKAq2I5c3O5auvDJFsAn0EP/3u8+pMnNHFBJA5kpacyP7lftBXOZzCq6El1pYpIjTtogDOzpWbWqgZriQtt188gyYIEXXi3hYXBbtptQSSODBg/m9nWl0CkOzXJgvReP8HrskSkjjlUC9yfgGVmNsHMEqvjzc3sCjPbYGYfm9nYSp43M3s08vw6M+t0uNd65dFAf5YH2/OzkkzOLJ7HbwJa800k3tyS+Qw+iC4R1My3i4y8Ak9rEpGakTM3lxUZF9F/XLan/+4PGuCcc88DHYETgXwz+42Z/arsdqxvbGZ+YAZwJXAOMNDMztnvtCuBsyK3EcATR3BtjcuZm8uMxN/T3beOwf6X8Ztxw/ktvS5LRKrBErso2pX6uWvMvFVbvC1IRGpE2vrp9PCv496EZz39d/9DY+BKgb1AMnDCfrdj1QX42Dn3qXOuBMgF+u53Tl/gGRf2JnCSmTU9zGtrXLv1M6LbZfXzr+STqVep+1QkTl2TuZin2z7FimB7xpSOIugcvbJXeF2WiFSzRuwB4GzfZn6b+s0PnF19DjUG7grgHeA4oJNzbqJzblLZrQreuxmwtdzjbZFjh3PO4Vxb4zacPYrtoZMIOXg1safX5YhINRsxeABftbyCZ5KmsCF5CAu+/plmpYrEqbIF+lv4dgJwnAUYEXres3oSDvHcBOBnzrn3q+m9rZJj7jDPOZxrwy9gNoJw9ystWrQ4kvqOWHimaXi26eXV+k4iUlt03/o4KVYCQDJBfln0OGxNh+ZdvC1MRKrUxVsfo76VRh8HHCT09G4I/qHGwF1UjeENwq1mzcs9PgP472GeczjXAuCcy3HOpTnn0po0aXLMRYuIlPda89srLNjtN8cHf9TkJZF4U+rCkals7OsnoWae/qLm5Tpwq4GzzKy1mSURbrpatN85i4AbI7NRzwf2OOe2H+a1IiLVLn34BBY1v4cilxj9j71FcLPWhhOJIzlzc2lkewEIAcuD7Vlx9n2e1uRZgHPOBYA7gJeAD4HnnXPvm9lIMxsZOW0p8CnwMfA0cPuhrq3hb0FEBAiHuCkdl7Ml1BiAFCum9watDScSL9qtnxFdNqiYRHpmve75Av2HGgNX7ZxzSwmHtPLHnix33wGV9kVUdq2IiFey0lPZs3Zv9PHpttvDakSkKmTkFVD61mzGJ37E/0IncZwVM6v+MMZ4XRjaSktEpMrMrj+M71wiQQd+QoQmNuDvmX28LktEjlKH/HuYmjiLE62IU+1r/mq9GDN2itdlAQpwIiJVZszYKZxb8gxB/PgsvFfq1e517dIgEoNmZd5If/9KzMITF8xgiFvsdVlRCnAiIlVoUNcWzAxcGZ3QsDeURMba7qyYcbu3hYnIEbnJLaoQ3pyDBcnpXpcVpQAnIlKFstJT2d55HE8E+hBwRoqvhCQL0u3LeV6XJiKHISOvgLHjf4UvsrxsCAg5eCWxJwPGz/a2uHI8ncQgIhKPstJT2b12OQn2/QJxJSToP1yRGPDnNzfzbvK86KzTba4JLe//uNYt0K8WOBGRavBa89uj3ajOgY8QeTMne1uUiPygu/3zSKEIgEKXzEvtHvC4osopwImIVIP04RN4J3hmdPxMPQvSd+t0ZmXe6HVpInIQOXNzGZmwGF9kzNszgcs9X+/tYBTgRESqSccH1rLRNYs+NoOb3Quw9S0PqxKRyvTKXkHb9TPwRbpOzWBAwnJPazoUBTgRkWq0ot19bAg2i+6X6jMomXkF2dPGe1uYiETlzM3l6a9u4ULfe+wMpfCdS6DQJfNa89o7e1wBTkSkGo0YPIB2WR9QSH0g3C2TZEGGFs3xuDIRKZO2fjotfTtIMEcjK+T+0qH0P+mvpA+vvVviKcCJiNSAWfWHVVhP6kT2kjtlmNdlidR5GXkFNGVn9LEZ3JP4F5aN6eFhVT9MAU5EpAaMGTuFvGC3aFdqgjl+WrxAuzSIeKxD/j2c5tsDhH+52ucSanXXaRkFOBGRGvKb4B0E8EfXlzIc76962duiROqonLm5bMg4p8J2WQuC3bj2pAW1uuu0jAKciEgN2X+bLb9BduIMtcKJeKDH+vtp5//8+1+oDBrbt7W+67SMApyISA3JSk/l9slzeSLQJxrimtkO5q3a4m1hInVMRl4Bp/l2A+GWt+9CCawJtmHD2aM8ruzwKcCJiNSw7Z3HsY9EAAIkMCthKqMffMrjqkTqjnmrtvBKsFN0YlF9C9Apa02tXbS3MgpwIiI1LCs9lUmlN1LkEkkiQE//OoZ+owAnUhN6Za9gpn9KhbFvqxM6el3WEVOAExHxQGKXYfgJRVd9b+vbxpDMx7wtSiTO5czN5bHdt9HT/140vI0rvYUuGcu9Lu2IKcCJiHggKz2VBcnp0WVFUqyYP7l7tUODSDVqt35GdOKCc7A8eC6JXWJzPUYFOBERjwwYP5vPjz8n+thvMKxotocVicSvjLwCHgn0J+DCzd57SeLiB/5NVnqqx5UdHQU4EREPNR/4CEH3/eMTKeLvmX28K0gkTr2/6mV+mbCAnMDV7HYpvNL8Lq9LOiaeBDgza2hmL5vZR5GvJ1dyTnMz+5eZfWhm75vZ6HLPZZrZ52b2TuR2Vc1+ByIiVaR5Fx6tPyq6rIgZXO1e97YmkTiTkVfAjIRsevrX0c/3Og0nfR4Ti/UeilctcGOBV51zZwGvRh7vLwD82jn3I+B8YJSZnVPu+WznXIfIbWn1lywiUj3GjJ3CgmC3aIgrxceajE7kzM31tjCROFH61uzodlmnRr7GOq8CXF9gTuT+HCB9/xOcc9udc2si978FPgSa1VSBIiI16Z206VxXkslul0Kyhejk/4RfbLxV3akixyAjr4DHJwxmauKs6MSFVxN7el1WlfAqwJ3qnNsO4aAGnHKok82sFdARWFXu8B1mts7MZlfWBSsiEkuy0lNZMHUMc+oPjQ6yNoM+6k4VOWpNV0/ltoTF0e2yghiX3/uCt0VVkWoLcGb2ipkVVHLre4SvkwLMB+5yzn0TOfwE0AboAGwHHj7E9SPMLN/M8nfs2HF034yISA3ptPd1EsxVODYr80aPqhGJbbcmLIm2vAWd8bfk/l6XVGWqLcA55y5zzqVWcnsB+MLMmgJEvn5Z2WuYWSLh8DbXObeg3Gt/4ZwLOudCwNNAl0PUkeOcS3POpTVp0qQqv0URkSq34exRrAm2YWfw+Og2Pze7FzQeTuQIZeQVAOFfhkLAlSe9wIDx8bNMj1ddqIuAoZH7Q4ED2jPNzIBZwIfOud/t91zTcg/7AQXVVKeISI0aMXgAnbLWkFb6dHRig8/CC5CKyA/LmzmZPRNP44a3f17uqLFsTA/PaqoOXgW4acDlZvYRcHnkMWZ2upmVzSjtBgwBLqlkuZDpZvaema0DLgbG1HD9IiLVasj5LWlk32IGJc7Pi8HO9Mpe4XVZIrVe962P08CKaOf/HB8QdPAnu9brsqqcOed++Kw4kZaW5vLz870uQ0TksOTMzeWnG39NQysk5ODJQB+2dx4XsyvHi9SElx/oy2Wly6MTF0IOfJNic+kQM3vbOZdW2XPaiUFEpJYaMXgAf2v7cPgHkMHIhMW8v+rlyNgeEdlf3szJXBoJb86Fw1u8LBuyP7XAiYjUcpvva0NL304g/EOpiESOu+YhSLvJ28JEapndE5vR0AoB2BBqRrv7P/C4omOjFjgRkRj2UrvJhMpttXWclcKrk7wtSqQWyZmby/b7WnIShTgHO0MpPOC7zeuyqpUCnIhILTdi8ACeT74uutWWc/D3oh97W5RILdJj/f009X2Nz8K/5JxoRTybeafXZVUrBTgRkRgwYPzs6LIiZtDdaTiICITXe2tWbohB0MF/ThngcVXVTwFORCRGvJM2nY2h8JbQ34WS2D2xGXkzJ3tclYg3yvY5zVjbnX0uEYCNoWZcedIieox63OPqqp8CnIhIjMhKT+XxE0ezPNie43wlNLRCum+N/x9UIpV5f9XLjExYTJIFaWSFbA41YcXZ98Xdgr0HowAnIhJDHrn7Vnpmvc7s+sP4ziVyInvDrRBaWkTqmF8mLMAXWevNDJrabkYMjv+u0zIKcCIiMegPey4kiQAJ5hiRsIR5q7Z4XZJIzdn6Fu1sMyEH34aSKHU+FiSne11VjVKAExGJQYO6tiBEuPnBj+ODpBvInjbe46pEasbW50bT1LcHn0F9KyVx0ldxtVH94VCAExGJQVnpqZQQHrhtBskW5I6iJ8iZm+txZSLVJyOvgDbjlvLdt99vjVVKgocVeUcBTkQkRr3SfDT7nD+6PlyihRi2caRa4iQuZeQV0CH/Hj5KGshZvs8BKHU+Xmo+xuPKvKEAJyISo9KHT+DakxZyXUkmgUiISzDH6KIZWl5E4s6zb26mn38lPgvvDVzi/PwmZRrpwyd4XZonFOBERGLYsjE9WDB1TLQ7FcI/3LS8iMSTnLm5/CfxNozwYr0BB/eV3sQjd9/qdWmeUYATEYkDW90p0fvOwTfBeh5WI1J1MvIKSFs/nab+PVhkqyw/kNhlmNeleUoBTkQkDqxodx8BF56VagYt/TvZc99psPUtjysTOTZNV0+lg+8TIPzLiXPwTvBMstJTPa7MWwpwIiJxYMTgATxW//bohAYzaOArguXTvC1M5BiMfvApRiYsxmfh4Dbb+mKT9tDxgbVel+Y5BTgRkTgxZuwUxpXeQiDSSlEc8rHn4zc0oUFi1thvp1TYbaGfe9XbgmoRBTgRkTiS2GUY35CCGSRZiAZWxMVbH6NX9gqvSxM5bDlzc1mb0YnT7OvosYAz5tQf6l1RtYwCnIhIHMlKT+UxG4Rz4RYLgBMpYsnXfdUSJzEhZ24uP934azr6P8EiXac7g8dzs2UxZuwUr8urNRTgRETizMTMBykqtzq9GSSa09IiEhPS1k+noRUSKjee833a8Gzmnd4WVssowImIxKFlzX8VndDgHIQcnEwh/87s7m1hIoeQO2VYdMbpF+4kvnH12RBqxoazR3lcWe3jyQZiZtYQ+AvQCvgM+Llz7qtKzvsM+BYIAgHnXNqRXC8iUlelD59A9rS9DC2aw2M2iPtcDmbQzb3rdWkiB/XT4gXRSQslJNB00mZOBNp5WlXt5FUL3FjgVefcWcCrkccHc7FzrkNZeDuK60VE6qQxY6fQcNLnTMx8kBL8AJTgJyOvwOPKRCrKyCug9dgllEb+ngKcZHs9rKj28yrA9QXmRO7PAdJr+HoRkTqlINgS5yCRIPev7cZbWT29Lkkk6pL8kXyaPIhEAuxzfgpdMrPr1+2dFn6IVwHuVOfcdoDI11MOcp4DlpnZ22Y24iiuFxERoKP/UyyyCbgZdA6sJXvaeK/LEiFnbi49/e9hBgkGAfz8v47/1IzTH1BtAc7MXjGzgkpufY/gZbo55zoBVwKjzOyIR9+a2Qgzyzez/B07dhzp5SIicWGJXRSd1ADhEDe0aM7BLxCpIe3Wz4guFxJmdX6brMNRbQHOOXeZcy61ktsLwBdm1hQg8vXLg7zGfyNfvwQWAl0iTx3W9ZFrc5xzac65tCZNmlTdNygiEkOuyVzMxhO6RH9IhhwstEu9LUrqtCGZj7E84yKWuS6sCbZhS6gxe1x9Xmk+2uvSYoJXXaiLgLLllIcCL+x/gpkdb2YnlN0HegEFh3u9iIhU1O43L3NdSSa7XQo+g+vdixRlNqZwYhN1p0qNypmbS47LpKd/HRP9sznXt4mW/p28EzqL9OETvC4vJngV4KYBl5vZR8DlkceY2elmtjRyzqnAv83sXeAtYIlz7h+Hul5ERA6t8JRODC/5DbtdCilWQn1KSbESdadKjeq9fgL1LQCEt3xLtBAlzq/13o6AJ+vAOed2AQe03Ue6TK+K3P8UOO9IrhcRkUNbNqYH0IOcuU3psf5+Wvi+JIgxvfR6EvMKNPZIql1GXgH3+XYD4XFv37j6mMGs+sMYM3iAx9XFDu3EICJSB40YPIArSh/kj4HeJBPgNv8LvL/qZa/LkjqgQ/49+AlF9+t9x53FiZP+p1mnR0gBTkSkjrrh/JYMT3iRRAvR0r+TvyTdz9jxv9JCv1ItMvIKmH/vVfT3r4wuZ7Pbpajb9CgpwImI1FFZ6aksSE6PzkxNtBD3JT7DvFVbvC1M4tL7q16mv38lFtkqyzn4W9uHGaFu06OiACciUocNGD+7wuP6lPLb1G88qkbiUUZeAa3GLmFqwtMVwtvy4LkKb8dAAU5EpI7bEDw92gpnBiNCz3tbkMSVS/NHsil5EG19nwMQdHBdSSZTG032uLLYpgAnIlLHnf3Ah8y2vgScUeL8FH28gs33tWFNRidy5uZ6XZ7EsIy8AnpEtskqa30rIokFU8dEZkTL0VKAExERbsl8hpWhc0myIPUtQEvfTjr5P+HmjbcxKfNur8uTGJQzN5fRb/cCwl2mzsE+5+eV5nd5W1icUIATEREANpw9ikKXXOFYooW4083zqCKJVRl5BfReP4HG/r3RlrfWxfO49qSF2mmhiijAiYgIEF4bbl7bR1kbbMNWmlDsfIQc/CtY6ZrqIpUqa3lr4dsJfD9hAVC3aRVSgBMRkagRgwfQMWsNnwabkmwhfAb9/SvZfm8Lr0uTGNF7Q8WWtyBwc2AcbU9N8bSueKMAJyIiB9hw9igCLvwT2AxO8++hzbilWuRXftDp9v02WSEHOYE+DDm/pVrfqpgCnIiIHGDE4AEkXPN7ApHB50EHHyUN5PY1VzH6wae8Lk9qqV7ZK1gUvICQgwXBbrQpnsftk+dqj91q4Mlm9iIiEgPSbuInS5ryfGAM7fzhNbya2h76ffNn4FZva5NapVf2ClK+XMO8xIdo5C/EDFrb/7jh/JZelxa31AInIiIH9e7E3pzm+75LLHxz6k6VqIy8AlK+XMNzSVk09hVGx74BanmrRgpwIiJySLPrD2O3S8ERHg/X0/8eHyUNJHNtN3KnDPO6PPFQRl4Bz765mRkJ2SRbEAiH/J2hFPLPvsfj6uKbApyIiBzSmLFTyO74EguD3QhFttzyGfgN+hfneVqbeCcjr4AO+fewKXkQp/n2AOHwtiDYjcb3f659TquZApyIiPygrPRUrntgKV+TUmFD8metj7pT66imq6fS378yuk1WWXh7qtFvvS6tTlCAExGRw/Za89ujG9/vI5GsfddzHhu49O3bNDu1jsjIK2Ds+F9xW8LiCmF+efBcrntgqZYLqSEKcCIictjSh09gXOktFLokfIT4JHkQ85My6elfx+8L72FW5o1elyjVqGzM28TEZ6KtbkEHTwT6MLXRZK/Lq1MU4ERE5IgkdhlGCUkkWxB/pPsMwl+HuMX0yl7hbYFSLTLyCih9azYfJQ+mHqVA+M/89VB78hqPUMtbDVOAExGRI5KVnkrDaybzHckEXfnlRSCBIOk7cxTi4tDp+dOYmjiLRHPfb5PljI1nj1J484AnAc7MGprZy2b2UeTryZWc087M3il3+8bM7oo8l2lmn5d77qoa/yZEROqytJs4LvNLrjxpEd9SPzqQ3WcwMmEx03aPIWdurtdVShUoG/M20v/3CmPeSp2PB+wXmm3qEa9a4MYCrzrnzgJejTyuwDm3wTnXwTnXAfgJ8B2wsNwp2WXPO+eW1kTRIiJS0bIxPdjuGgLft8L5DDr5P6Hd+hkeVydVofSt2UxJnBUd8+YcbAiezv0dX2di5oNel1dneRXg+gJzIvfnAOk/cP6lwCfOuc3VWZSIiBy5Fe3uY0WwPdeVZEaPOQePBPrTauwSLTESo3Lm5rI5ow1TE2fhKzfOcUWoPVcEHtIuCx7zKsCd6pzbDhD5esoPnD8AeG6/Y3eY2Tozm11ZF6yIiNSMEYMH0CPrdRZMHUMx/ujx+UmZbEgeQulbsz2sTo5W2/UzaOnfWaHbdEOwGY8G+jNEe5x6zlzZgj5V/cJmrwCnVfLUBGCOc+6kcud+5ZyrNISZWRLwX+DHzrkvIsdOBXYCDsgCmjrnKt3PxcxGACMAWrRo8ZPNm9WIJyJSXbKnjeeuohkV9sN0Lrxm3EvNx5A+fIJ3xclhy5s5mYu3PkZpyMfJvr34CK/zdvED//a6tDrFzN52zqVV9ly1tcA55y5zzqVWcnsB+MLMmkaKawp8eYiXuhJYUxbeIq/9hXMu6JwLAU8DXQ5RR45zLs05l9akSZOq+eZERKRSY8ZO4S/J1xEoNzPVDOpbKVdv1XipWNArewUXb32MBlZEI99e8oLdaF08j3+mPel1aVKOV12oi4ChkftDgRcOce5A9us+LQt/Ef0ADbAQEaklBoyfzaSOK/m/kud4ItAnunNDAo75917FeZNe8rZAOdDWt/jw/s6syejE1F13cSJFQDh89/OvxG+mMW+1TIJH7zsNeN7MbgG2AD8DMLPTgZnOuasij48DLgdu3e/66WbWgXAX6meVPC8iIh7KSk8lKz2V1mMdXe0DOvk/xQz6+1eyuvhsoLfXJUo5Hz1/Lz8KbQT/962mEL6/MNiNQV1beFugHKDaxsDVRmlpaS4/P9/rMkRE6oyyrZc+Th5EQiQUFDsfD3R8XS06tUSv7BWkfLmGOQlTSPGVRJcL2eWOZ0Tp3fy46+X6s/LIocbAedUCJyIidUBWeiqrNu0iZ2ef6ObnSYSYuLYbj6/uw1NJN/LuRLXGeSUjr4CNXxQCbaPhDWAfCTzSaRkLFNxqLQU4ERGpVsvG9KBXNvxo12Z6+t/DLPzD57aExWwpPpXWYwPccH5LtfLUsIy8ApqunsqnyYspmzTsHAQwnqx/q/48ajl1oYqISI15fMJgRiYsji4Mu8/5Obv4WQA+m3a1h5XVHb2yV0Ra3WBj8hCSLBh9zjm4r+NKhbdawpNlRERERPa3vfM4flqSGZ2ZmkyQTcmDeDfpFoZkPuZtcXEuI6+ANuOW8uiukWxKHsSm5EEkEowu9+IcbDnuRwpvMUItcCIiUuNefqAvl5UuP2DBXwh34S1pfrcW/a1iPx3/eyYn5NDWPq/wuQOEHPgm7fGmMDkotcCJiEitcvm9L/CX5OsIue+Dm1n4lmiO7lsf1x6qVSQjr4DWY5dwh38+7XyfV9iUvuz2amJPr8uUI6QAJyIinhgwfjYTO67kiUAfgvsFipMp5JL8kVr09xiNfvAp7lrTm0+TB3Ghbx0hB0UugXGlt9C6eB6ti+dxX8eVXH7vodbTl9pIAU5ERDyTlZ7K9s7jaBMJE1+REm2J6+l/j/zQ9WRPG+91mTEpI6+Aft/8mUa+wvDMXwOfQRH1yA1digFDNPs3ZinAiYiIp7LSU2l7agoA00uvp9RZdDeARHMMLZpDq7FL6JW9wuNKY0f2tPFkrL2IHr510VbNgIMil8j00usZcn5LNk27WuEthmkdOBER8dyyMT0i965mRcZqevgjwYNwd+qm5EHs2n08bM2D5l08rLR2K9v5Yk3yHJIsFD3uHPxf8TwMtOZenFCAExGRWmXD2aNw62fQ1HbRzvd59Hhj/16+zrmGjqWzFUIOovzCvOUXmVgePFfdpXFGAU5ERGqVEYMHkJGXSr+3bwQqzlJt4NuHA+at2qIwEtF67BLKstrG5BejiyRDeHmQM4vn0aB+Au/q84orCnAiIlLrZKWnkrP3Hr5ZP4NHA/2Z4H+GTv5PWRM8EwCf7/vgUpdbljLyClibMIwG/n0AlIR8uHKj2xcGu9XpzyeeaSFfERGptcrGdAH8MWEqPf3vARAkPAvPgF3ueEaU3k3hKZ3KjaWLX6MffIpffTudM2wnAXwkEap0Yd42xfPU1RzjDrWQrwKciIjUehl5Bdy/ttsBQaXM8mB7biodC8Rvi1xGXgHzVm1hVsJUevrXHfC8c1AYSuJ4XwnvntyLjnf91YMqpSppJwYREYlpWemprE7oWGFJjFC5hX97+NaxKXkQHyUPpvSt2V6XW+XOmrCUZ9/cTNA5Hg30Z3OoMUEHxc5HifMTcMYTgT50CMxhYseVCm91gFrgREQkZpS1QgUjP7s+TR5UYdA+wG6XQnbHl/jzm5tjfoxc9rTx3FqUQz1KKf9tOsLj234dGBU9Fsvfp1ROXagRCnAiIvGhLMhN9/+Bfv6V0XBT9hMtgI8SEnmg9AZyQ5cChGdiTuztSb2HIyOvIBo6y6xJHkFDK6z0/JCDs0qe45OpV9VMgVLjFOAiFOBEROJLWZDz+WBVwi8OCDu7XQqdinOijxvUT6BwX5BBXVvUutaqNuOWRgNpCPCXe85BpS1w76RNr3Xfh1QdBbgIBTgRkTiW/yeK/v5rkggQrKQF7lC86H4s3+LW9tQUNn5RWGmXMFQMotpNoe5QgItQgBMRiW9lLXKDurbgkvyR9PS/V6FLcv+Wrf0ttou4q/j2Km2hK19T7uotlAYP/nP34YQZB7TAFeMns/QmErsMU2irYxTgIhTgRETqDpfZgIOsOnJQZTsXlDG+H1fX9tSU6DpzZaGszSnH88mXe6OBr1f2CjZ+URjtqk2p52dPUeCov4fa3OUr1a/WBTgz+xmQCfwI6OKcqzRVmdkVwCOEfxGZ6ZybFjneEPgL0Ar4DPi5c+6rH3pfBTgRkTrkmX7w6T+BcAhz7odb4MqUjTkrP/bMEX7g3PdrcAX3e719zk89C0av2//58gLOSLCKP4ODGPtI5JXmd5E+fMJhVCrxrDYGuB8R/nf0FPCbygKcmfmBjcDlwDZgNTDQOfeBmU0HdjvnppnZWOBk59xvf+h9FeBEROqu8rs67O9PidMqXRzXM/Ubwm83eV2FeOxQAc6TvVCdcx8C2MGW1A7rAnzsnPs0cm4u0Bf4IPK1Z+S8OcBy4AcDnIiI1F1Z6akVuiHLB7pHA/3paB9xohUBh2iBi3ytiha4oPnwE4q+Sfh9fJBYDy6deCzfqtQBtXkz+2bA1nKPtwFdI/dPdc5tB3DObTezU2q6OBERiW3lA11GXkvOe7Mt8P2M0OqkRXflWFVbgDOzV4DTKnlqgnPuhcN5iUqOHXF/r5mNAEYAtGjR4kgvFxGROmD/1rmyyQgQ/mF01qkpfPRF4QE/hMpPcijT9tSU6MQGIDoDVYFNqlK1BTjn3GXH+BLbgOblHp8B/Ddy/wszaxppfWsKfHmIOnKAHAiPgTvGmkREpA4om21aFRTcpDrU5s3sVwNnmVlrM0sCBgCLIs8tAoZG7g8FDqdFT0RERCQueBLgzKyfmW0DLgCWmNlLkeOnm9lSAOdcALgDeAn4EHjeOfd+5CWmAZeb2UeEZ6lOq+nvQURERMQrWshXREREpBY61DIitbkLVUREREQqoQAnIiIiEmMU4ERERERijAKciIiISIxRgBMRERGJMQpwIiIiIjFGAU5EREQkxtSpdeDMbAewuZrfpjGws5rfI97pMzx2+gyPnT7DY6PP79jpMzx2sf4ZtnTONansiToV4GqCmeUfbNE9OTz6DI+dPsNjp8/w2OjzO3b6DI9dPH+G6kIVERERiTEKcCIiIiIxRgGu6uV4XUAc0Gd47PQZHjt9hsdGn9+x02d47OL2M9QYOBEREZEYoxY4ERERkRijAFeFzOwKM9tgZh+b2Viv64k1ZjbbzL40swKva4lFZtbczP5lZh+a2ftmNtrrmmKNmdUzs7fM7N3IZzjJ65pilZn5zWytmS32upZYZGafmdl7ZvaOmeV7XU8sMrOTzOxvZrY+8v/iBV7XVJXUhVpFzMwPbAQuB7YBq4GBzrkPPC0shphZd6AQeMY5l+p1PbHGzJoCTZ1za8zsBOBtIF1/Bw+fmRlwvHOu0MwSgX8Do51zb3pcWswxs18BacCJzrk+XtcTa8zsMyDNORfLa5h5yszmAK8752aaWRJwnHPua4/LqjJqgas6XYCPnXOfOudKgFygr8c1xRTn3GvAbq/riFXOue3OuTWR+98CHwLNvK0qtriwwsjDxMhNv+UeITM7A7gamOl1LVI3mdmJQHdgFoBzriSewhsowFWlZsDWco+3oR+e4hEzawV0BFZ5XErMiXT9vQN8CbzsnNNneOR+D9wDhDyuI5Y5YJmZvW1mI7wuJgadCewA/hjpyp9pZsd7XVRVUoCrOlbJMf3mLjXOzFKA+cBdzrlvvK4n1jjngs65DsAZQBczU3f+ETCzPsCXzrm3va4lxnVzznUCrgRGRYaYyOFLADoBTzjnOgJ7gbgam64AV3W2Ac3LPT4D+K9HtUgdFRm3NR+Y65xb4HU9sSzS3bIcuMLbSmJON+DayBiuXOASM/uztyXFHufcfyNfvwQWEh6mI4dvG7CtXAv63wgHurihAFd1VgNnmVnryGDJAcAij2uSOiQyAH8W8KFz7nde1xOLzKyJmZ0UuV8fuAxY72lRMcY5N845d4ZzrhXh/wf/6Zy7weOyYoqZHR+ZiESk268XoNn5R8A59z9gq5m1ixy6FIirCV0JXhcQL5xzATO7A3gJ8AOznXPve1xWTDGz54CeQGMz2wZMdM7N8raqmNINGAK8FxnDBTDeObfUu5JiTlNgTmRWuQ943jmnZTCkpp0KLAz/TkYCMM859w9vS4pJdwJzI40qnwI3e1xPldIyIiIiIiIxRl2oIiIiIjFGAU5EREQkxijAiYiIiMQYBTgRERGRGKMAJyIiIhJjtIyIiMhRMLNGwKuRh6cBQcJb9wB0ieyJLCJSLbSMiIjIMTKzTKDQOfeQ17WISN2gLlQRERGRGKMAJyIiIhJjFOBEREREYowCnIiIiEiMUYATERERiTEKcCIiIiIxRsuIiIiIiMQYtcCJiIiIxBgFOBEREZEYowAnIiIiEmMU4ERERERijAKciIiISIxRgBMRERGJMQpwIiIiIjFGAU5EREQkxvz/leYyco4zqzAAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.figure(figsize = (10, 5))\n", "plt.scatter(t, y, label = \"All Data\", s = 4)\n", "plt.scatter(X[:, 0] , allpred, label = \"Predicted Values\", s = 2)\n", "plt.xlabel(\"T\")\n", "plt.ylabel(\"Y\")\n", "plt.legend()" ] }, { "cell_type": "code", "execution_count": 35, "id": "c05ebeb1", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 35, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnAAAAE9CAYAAACLPV+MAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABat0lEQVR4nO3deXhU1f3H8feZO5lkkkAWCFtYAkgwhLCLBaJsokFFLEhFqEVFMSi44AJq+8O6IlqrUUtKVVAqpeKuKKhUioIaQowCIqIIhH2y75nMnfP7Y8YUNSBLkptJvq/nGWfuNvPJPAjfnHPPOUprjRBCCCGECBw2qwMIIYQQQoiTIwWcEEIIIUSAkQJOCCGEECLASAEnhBBCCBFgpIATQgghhAgwUsAJIYQQQgQYu9UBGlLr1q11XFyc1TGEEEIIIX7V5s2bc7XWMbUda1YFXFxcHJmZmVbHEEIIIYT4VUqpPcc6Jl2oQgghhBABRgo4IYQQQogAIwWcEEIIIUSAaVb3wAkhhBAiMFRXV7Nv3z4qKyutjlLvQkJC6NixI0FBQSd8jRRwQgghhGh09u3bR4sWLYiLi0MpZXWceqO1Ji8vj3379tG1a9cTvk66UIUQQgjR6FRWVtKqVasmXbwBKKVo1arVSbc0SgEnhBBCiEapqRdvPzqVn9PSAk4p9bxS6ohSausxjiulVJpS6jul1FdKqQFHHUtRSu3wH5vXcKmFEEII0dRprUlOTua9996r2ffyyy+TkpLC6tWr6dmzJ2eccQYLFiywJJ/VLXBLgZTjHB8L9PA/ZgCLAJRSBvCM/3gv4AqlVK96TSqEEEKIZkMpRXp6OnPmzKGyspKysjLuuecennnmGW688Ubee+89vv76a/71r3/x9ddfN3g+SwcxaK3XK6XijnPKeOBFrbUGPlNKRSql2gNxwHda610ASqkV/nMb/hsUTcsH81n51bM8HhWBW2vchkGYaXJbQRHPRbTkUJCdlJIy9jqCKFOKKpvisM2Gx/a/34VsXi82FEqBBmxaE2N6qVKKfLuB0zQpNwwcWjM3r4BJpeVkBwcxN6Y1h+wG0aZJS9NLlU2RbzNo7/Hw58j+9PvdCuu+FyGEaIZ69+7NuHHjeOSRRygrK+MPf/gDR44c4YwzzqBbt24ATJ48mTfffJNevRq2Hamxj0KNBXKO2t7n31fb/rMbMJdoqj59mrSObSk1jJpdZXY7adGRFPr3rWoRBse5X8FrGHh/tm//z94PoEop0qIjmVRaRnpkBAeCfPtz7XZyj/o/c5fhID0vk/TT/NGEEKKpq3CbfO8qpXtMOE6H8esXnID58+czYMAAHA4HmZmZvP3223Tq1KnmeMeOHfn888/r5LNORmMv4Gr7V1IfZ/8v30CpGfi6X+ncuXPdJRNNyubd+SxYvZ3rIy5jVv47PBH90xa4m+qpBW52fiGmhtTCIn6wB3EoqPYWuNRWg6z7coQQIgBUuE3Of+K/5Ja4ad3Cwfu3DK+TIi4sLIzLL7+c8PBwgoOD8XUK/pQVgy0aewG3D+h01HZH4ADgOMb+X9BaLwYWAwwaNKjWIk80YTkZPL5qOstCNIPLK/giJBivUrTymIRoTZjWXFzq5j9OB78rL2VFaBhXlymGf5fCCu9oAEqAOwByfW95Kh2ZeUe9Ljrq9Z3+B1X/O1AETOrfjrF9OjGke6s6+y1SCCGasu9dpeSWuKmoNsktcfO9q5TesRF18t42mw2b/xf1jh07kpPzv07Affv20aFDhzr5nJPR2Au4t4BZ/nvczgaKtNYHlVIuoIdSqiuwH5gMTLEwp2is1i1gWYjGoxQbQ501XZ8HHP9rMcsJslNoGOwI8T3bleahqn+zomp0rW/Zp10YpdVe3B4vUU4HYcEGuaVVlFab2G2+3wK19tIhMozOrUN9f6FUVrPlUNkJx175xSFWfnGIEDucFdeKeWMTSKyjv4iEEKIp6h4TTusWjpoWuO4x4fXyOWeddRY7d+7khx9+IDY2lhUrVrB8+fJ6+azjsbSAU0r9CxgBtFZK7QPmA0EAWut04F3gQuA7oBy42n/Mo5SaBawBDOB5rfW2Bv8BROM3Yh5X/koL3NjiKv4bFswF5aWsCQ3jysJKFlZfDkBEiMJht+O027i4byzXntOd6HDHKUXZuNPF3a9toaSiAsPu4Eip+1ev6Rr0CeWeN7lnlUmRYSP3qHvpWiuDv459gX5t+p1SHiGEaEqcDoP3bxle5/fA/Zzdbufpp5/mggsuwDRNrrnmGhITE+vls45H1daX21QNGjRIZ2ZmWh1DNIScDFh9F5nuQp6yezhy6EK2VyT/6mUT+rbDVVbN2V1bMeXsLqdcrJ2I/QUVvPTZHlBevj1UyoffuH5xzuDOf2R7mOeY7zGswzDSx8jwBiFE07N9+3YSEhKsjtFgavt5lVKbtda13gTd2LtQhTgllR8+RMj+TJ5tG0NWiJOEmNWwt/YC7sbh3ejSOozzEtrVa8H2c7FRTu4ce2bN9v6CCpZs2MWRokre2nIIgBJXCvG2NzFtv2yBA9iw/xP6LUnEoTV35Bcxqe+1MObPDfYzCCGEsIYUcKLp+GA+L3/5Dx6LjgI0t4eHMSnfpFgrCly/nC/6goQ23Dk2ge5t6uc+iZMVG+Xkjxf7muHnFlSQ/t+dLPssGXb/svDsFH9HzbQmplJUAGnREUz69Gkp4IQQohmQAk4EtJU7VvL4Z/cTUe0m2vTyQ3QkFYZvgMJT0ZG8uaeIjXsfrjnfBswedQbThnZt0Na2kxUb5eT+S/tw63ln8srmvXyw/TCbdhfWHE/Obc0HMXlUKYWhNYaGUNPLlLatuPORDvQb/RAMusqy/EIIIeqXFHAioKV9kUYpmtKgIPYHQbcqNyYK0MzML2Fh9eU4FJzZoSVX/qYLF/eNDahpOaLDHcwYfgYzhp/B5t353PbyF+zOr+Rf+XMh/3/nXdDlZjaGOjkApLf08vQH92KXAk4IIZosKeBEYMlcysqP7+XRiFC8ShHi9aJsBq1Mk3amScvDyXxZPhGAu4H5FyXw5dldAqpoO5aBcdGsu3M0G3e6uHnFZlxlZs2xhLwulNhyKFOKg4bBb1vZcf8jgWJDMSd6IJMuXWZhciGEEHVNCjgRUFZ+ch/3RYXXzOdW5Z9Y8VB5Aj/kTK85r1/Hljz42z5Ncu60oT1i2PSnFFZvOcC8V7+ksNLLE6VzoBRGdZnDptCfdg2n5W9mkkVZhRBC1A8p4ETgyMkgraV/Ml6tCdaaVh6TIDOYr3PPAyDYgHdvHt5oBibUp5SkDqQkdWDz7nwuS/8UDQS5hpHQdgPlyqDUpsm3G5SjWPlYByaVlsOwm2WQgxBCnIC8vDxGj/ZN6H7o0CEMwyAmJgaAjIwMHI5j30edmZnJiy++SFpaWr3lk3ngRKO28rOFPP71C1Qr3+LvCgj2emlx5Fx2FY77ybl3p/TkyqFdm0R36cmqcJu88+V+5r26hR87Vp2dnsMevhMAh+ml2qbQQGtT097m4M6z7qBfb1nARAjRODWmeeDuvfdewsPDuf3222v2eTwe7Pa6awc72XngbLXtFKKxSPvmn5QaNl9XqVJopQiGnxRvwQasnTOcGSPOaJbFG/hmIJ90Vme23pfCk5f3BaAq9zyiqoIJN008Nt93h1Lk2m1ssXlI3/ykxamFECKwXHXVVcyZM4eRI0cyd+5cMjIyGDp0KP3792fo0KHs2LEDgHXr1nHxxRcDvuLvmmuuYcSIEXTr1q3OWuWkC1U0Ko9/dAdL9qwGfC3DEaaJUym8SlGlFDYNnY/0Zp///Obc6lYbp8NgfP+OxIQHM+U52LvL113aKeYlClt9BUCQaVJtGBRXFpCa3pPUXtPod+7dVsYWQoiA8e233/Lhhx9iGAbFxcWsX78eu93Ohx9+yN13382rr776i2u++eYbPvroI0pKSujZsyczZ84kKCjotHJIAScalWV7V4MC/38ostuJNE165ZzNGv/o0o3+c9fOaR73up2KoT1iWDtnOH9+awvrv8snxzUVXFMBaHHmXFCKLSEO3/2E3ywjXQo4IURTkJMB6xbAiHnQaXC9fMSkSZMw/BOpFxUVMW3aNHbu3IlSiurq6lqvueiiiwgODiY4OJg2bdpw+PBhOnbseFo5pAtVNCpXhvUArWseDtNLoWFQ0PqLmnMu6dOWrD+OkeLtV3RvE86L1w4hfWr/n+zvURwJWpNY6SbSNBldXAwfzLcmpBBC1KV1C+D7tb7nehIWFlbz+k9/+hMjR45k69atvP3221RWVtZ6TXBwcM1rwzDweI69xvWJkgJONBrZ6x8iM28rSVVu0g8U0mPHdOYeDCKpTFHiXwrrsYlJpE0Z1KhXUWhsUpI68GrqkJrm9qwDd1HyzSN4vWEUGgZrw5x4Nz5taUYhhKgTI+ZB99G+5wZQVFREbGwsAEuXLm2Qz/yRFHCiUche/xCzv/snW0KC2RISTHpkJFk6njtLHmTj3ofZXpHM2jnDueyszlZHDUgD46LZcl8Kt4w6o2ZfiSuFAWUm8VVuhnVqy0VLBpJ9JNu6kEIIcbo6DYYrX6u37tOfu/POO7nrrrsYNmwYpmn++gV1SKYREY3CjPSefOp0EG6adHBD+eHxbK/wLeL+uwEdmHdhorS61ZFXNu3l9le31Gx3ir+DQv/9HB3cHlppLxjB3Hn2PJlmRAhhmcY0jUhDkGlERMB55sXfscVh0L3KzaLDLmbvb1VTvC2fPpiFv+svxVsduuyszmyYO4pwh2+gSHJua989h8CBIIMtwQ622LVMMyKEEI2YFHDCGplL0Q+0J+vhtjzr/ZpSwyDPbpBY6SHNMwGAV1OHMLRHjMVBm6bYKCeb/ngBd6f05F/5c6k4NIEgz1FD2rXmU13G4+9cY11IIYQQxyQFnLCE98N7UZ5yFke2wKMUdq2ZmV/Cn6qvZgvxrJ0znIFx0VbHbNKcDoMZI84gfWp/PIVnk7/z/poRqvjn3nsxN0NGqAohRCMkBZxocBVukxcrz8WjFdcVFDO0vIJp+9pyd2467waN4fN7ZIqQhpSS1IG1c4aj8I1Q7bhnHEGmr0u1XbVHRqgKIUQjJAWcaFCZb93I9CUJvN0uk60hQZRU9GTNnid5onQOd5zfg8/vOV/ud7NA9zbhLPLPF7e9IpnCnBvoVBZOn0o3/Tu35/zn+5C9dbnFKYUQQvxIVmIQ9S8ng5X/mceT5NO6qoLvQ3wTGj7k7YineAJBhmLFdb+RLlOLpSR1IO1yLzf9+0u8FV34eu8f2d9zHl4bHEQz99P76Jr9N1LPf5p+bfpZHVcIIZo1KeBE/Vu3gDTzEEWGgWk3SKqs4rCOJNN1FVp34ZPbRxIb5bQ6pQAu6d+R2KhQJqZ/CkDL/CTfGqpKcSDIzgGzCL5MJ31MusVJhRCifuXl5TF69GgADh06hGEYxMT4BtZlZGTgcBy/t2jdunU4HA6GDh1aL/ks7UJVSqUopXYopb5TSv1i2mSl1B1KqWz/Y6tSylRKRfuP7VZKbfEfk8ndGqOcDFg2gYIuKaTkheP0KG7OLyZ571l8v3s+juqufDJ3lBRvjczAuGjWzhnO2V2jyHFNJelQL+z+gQ1O0yRz38f0eT6RpKW9uf6ty62OK4QQ9aJVq1ZkZ2eTnZ1Namoqt956a832rxVv4CvgNm7c+KvnnSrLCjillAE8A4wFegFXKKV6HX2O1vpRrXU/rXU/4C7gv1rr/KNOGek/Xuskd8Ji6xaw8vBnXPDdk7wWXUG4ayjzcv/OQtM3Oeyqm86R4q2R6t4mnKVXn00Lh42NhdMYvHcwvSrduG02qmw2tGEDpdhY8LXVUYUQosFs3ryZ4cOHM3DgQC644AIOHjwIQFpaGr169aJPnz5MnjyZ3bt3k56ezl//+lf69evHxx9/XOdZrOxCHQx8p7XeBaCUWgGMB471L8IVwL8aKJuoA+udbXiwdTSm8k0YWx6zAQrHoYB3ZifLSNNGzukwWH3rCIY98h/WlE9klHcDplLYtEZ7NdqmGBqdaHVMIYRoEFprZs+ezZtvvklMTAz//ve/ueeee3j++edZsGABP/zwA8HBwRQWFhIZGUlqairh4eHcfvvt9ZLHyi7UWCDnqO19/n2/oJQKBVKAV4/arYH3lVKblVIz6i2lOCXfHyllqeujmuJNaU2oaxgj41uz+Y9jSIyNsDihOBGxUc6aKUaCXMMYWl7BCwcP89XefdydW8DXuVtYmdbD110uhBAWyz6STeoHqfWyrnNVVRVbt25lzJgx9OvXjwceeIB9+/YB0KdPH6ZOnco///lP7PaGaRuzsgVO1bLvWAuzjgM2/Kz7dJjW+oBSqg3wgVLqG631+l98iK+4mwHQubMshF6fsrcuJ33TY7So9rDG4cUeYqeDu5oorybi8DDWlI9j2fV9ZJqQANO9TTifzB3F6Mfg0pxIzgj6J2U4+Ft0BIWGQVo4TFq3wLeAtBBCWCj9y3Q2HNjge13Hg6201iQmJvLpp5/+4tiqVatYv349b731Fvfffz/btm2r08+ujZUtcPuATkdtdwQOHOPcyfys+1RrfcD/fAR4HV+X7C9orRdrrQdprQf9OHpE1I/0zU+ywVbNaocXrRTVNhulhsFnPzzOmvKJvJo6RO55C1CxUU423jWGNxhN76olJFYtZXpeJeGmSZjXy5TyrWS/PNnqmEKIZi61byrDOgwjtW9qnb93cHAwLperpoCrrq5m27ZteL1ecnJyGDlyJAsXLqSwsJDS0lJatGhBSUlJnef4kZUF3Cagh1Kqq1LKga9Ie+vnJymlIoDhwJtH7QtTSrX48TVwPrC1QVKLY7qk2wwM/2hFtMbh9dLliO8eqVdTh8g8bwEuOtzBxrvHEBbkazx/K382SVXV7A8KYktIMLNLv5LJfoUQlurXph/pY9LrZa5Km83GK6+8wty5c+nbty/9+vVj48aNmKbJ73//e5KSkujfvz+33norkZGRjBs3jtdff73eBjEorY/Va1n/lFIXAk8ABvC81vpBpVQqgNY63X/OVUCK1nryUdd1w9fqBr5u4OVa6wd/7fMGDRqkMzNlxpE6lZMBq+9iuS7lcaOKEeWlfBbmoPOR3mwsnIbDgPduHi4DFpqQCrfJ1c9/zme7C7gg9FVK2n7CD0F2Sg2DwWYIz12zyeqIQogmYPv27SQkJFgdo8HU9vMqpTYfa6YNSwu4hiYFXD1YNgG+X8s5nWMpNAwiTZN39uTTt+o5FMg8b03U/oIKhj3yHwAGqG+5Kvwf/KWtJt9ucF6rwfx13PMWJxRCBDop4I5fwMlaqOKUPf7KBPqa3zKwc0fOKvPQwqO5Kb+QQ15fV+k7s5OleGuiYqOcrJqdDECWjqdlZSvy7QYoxYd5GUx5TtZOFUKI+iQFnDg1ORm8WPotXqVwGzY2hdk59/vziSnqzl2e61g1O1mmCmniEmMj2DB3FCEGpHkmMLTEd+8jSrHFrln46QOkvj25XobzCyFEcycFnDgl2W9Mx+bvflfa1/KWYmziRnUPT92ZKsVbM/Hj6NSttnjW7F/A0L2DSaysIqmyCrweNuRvI/1LWTdVCHFqmsttXqfyc0oBJ07JInsF1TbfH5+ulYqYou6keSbw/q3Dpdu0mYkOd/D6Db7u1DXlE/nbgTKWHzzMnfkF9ClXXJ0o82wLIU5eSEgIeXl5Tb6I01qTl5dHSEjISV1n5US+IgBlr5rFwn3vU24YdKty49SaosMTuao6mSXTBkrx1kwlxkawfPpgpjyXwcLqy7kraDnO8mgqiq+mZdbHsGwijJ4Pg66yOqoQIkB07NiRffv24XK5rI5S70JCQujYseNJXSOjUMVJuX5RDzaG+n5LGFJeSe7eG8jS8VwzpDP/Nz7J4nTCapt35zMx3TfJpc25h7j2/6DUXkUHj8m9xVX0u323tQGFECKAyChUcXpyMlj5wiiSl51FhyobvSurSKx0o1znkKXjiXDAHWN7WZ1SNAID46JZO2c4IXYbwa0/xBXsocIw+D7YwR9ah7Pys4VWRxRCiCZBCjjxq7LX3M5D+ghF3kreiAziilx7zfJY153Thc/+mILTYVgdUzQS3duEs/a2EVTlnkdMld03MhXQSvHk1y9anE4IIZoGKeDE8WUuZVH1ATz+5bE8SvFmtAeAcIeNOWMSpHgTvxAb5WT2kNHs2vUAFYcm4DS9GKaXEhs8sniQbwUPIYQQp0wKOHF8a//MzMIihpVXcKPLzVnlbuwu34jDN2edI8WbOKZpQ7sC4Ck8m6u+74eyKbxKscJRifnRwxanE0KIwCYFnDiufXGXkVjp4a4DmvfzZ/KfPb6u01dTh8j6puK4osMdrJ0zHICF5hTC8/pgeOHcEjc3lG0h+6kkaYkTQohTJAWcqN0H89F/boW59TWClJcf6ECWjifIUKydM5yBcdFWJxQBoHubcDbMHUWQDXJcU+n27TXsDlZsDHUyN6Sa1NXTZaUGIYQ4BVLAiVrpDWko7aGjLY91Zh/SPBMAeOOGYdLyJk5KbJST1bf4WuJusr/GEbuv2/1wkJ0NDkj/fIGV8YQQIiBJASd+Ye2/Z5DathXZwQ5MFFdVzyNLx7N8+mBZIkucku5twlk0tT9pnglMyguipcfLHwqLSaqsovjgl7LwvRBCnCRZiUHUyN66nEcyFrIbN6WhThTw2KEizmwXxnPTzpZVFsRpGdGzLbfb48nKewjy4BP1LZGd09kUamfR5if4e+8pVkcUQoiAIS1wokb65ifZapiUGgaRpsk1BaU8UH0lQ7q2luJNnDanw+Ct2cNrtrN0PAePTCChzE7P4AGkvnC2tMQJIcQJkgJO+ORkcFW5Qe/KKpIqq3jqsIvyih6s8I7m2nO7W51ONBHd24SzanZyzfbO4CB2BHtYUvkJGygnffOTFqYTQojAIQWcAGDhquuY6SymV7lm0YFinOUxpHkm8GrqEGl9E3UqMTaCDXNH4TAgOGYNXjvgnyh6TMdrrI4nhBABQQo4wUer72VZiIlHKV6OCuZO40IuqH6U26/5vUwXIupFbJST129Ipsp1ATYPvuKtpJIOn9/HjGf7SleqEEL8CingmrOcDFg2geU/LK9pAUEpqmM2MDaxDUN7xFidUDRhibERTImfRM8qOyhFuR0WtA7n0yAv8z97wOp4QgjRqEkB15y9ci3Z+z6mzAZJlVXclVvI2eVVBLmGcfsFCVanE83ALef3pMSVwuDyKlILi2rmiNtltzF4SW9WfrbQ4oRCCNE4SQHXXGUuRRftIT0ygi0hwbTwenkl7zY+3PNXpl3xiEzWKxpEdLiDx66dR8HemTgrYrgpv6imJbjCpkjb8ZLVEYUQolGytIBTSqUopXYopb5TSs2r5fgIpVSRUirb//i/E71WHJ/ng3tRQGpBEUPLK2ib27tmsl7pOhUNKTE2gnuun8YF7ke5K3cxd+UW4jRNnKaXiTG/tTqeEEI0SpYVcEopA3gGGAv0Aq5QSvWq5dSPtdb9/I/7TvJa8TPZW5eT8lwig9qH83hkBP3cbq440J4Xy65jybSBUrwJSwyMi+axiUkAfJX/O/6zx8V/9rg49Ekx+wsqLE4nhBCNj5UtcIOB77TWu7TWbmAFML4Brm2+cjJI33g/++02TKVYEtmSjY4WpHkmcOOIboxMaGd1QtGMXXZWZwZ1jmSFdzRuHLzXws5/z1jNH14bzJR/JJC9/iGrIwohRKNhZQEXC+Qctb3Pv+/nhiilvlRKvaeUSjzJa8XRVt9FZFWl7x4jAKV4OjKaLB3P9GSZrFdY75HL+gKwsPpy0qIjKTUMjthhi8NO+jfLLE4nhBCNh5UFnKpln/7ZdhbQRWvdF3gKeOMkrvWdqNQMpVSmUirT5XKdatbAl5OB59DXrGkRVjNlSMcqTYHrIl5NHUJ0uMPqhELQvU04r6YOYYV3NMm5rQk3TSKqwWl6OaiQVjghhPCzsoDbB3Q6arsjcODoE7TWxVrrUv/rd4EgpVTrE7n2qPdYrLUepLUeFBPTTO/vysng5Vd+xzmxUQSbJmhNh6L2bN/1CAnxl8lkvaJRGRgXzZJpA/lX/lwOfvsoeVU9qDBs7HLYWbh9qW/+QiGEaOasLOA2AT2UUl2VUg5gMvDW0ScopdoppZT/9WB8efNO5Frhl5NB9soreDgqnFLDoMzumzR1r70FAH8cl2RxQCF+aWRCOx71D2qoyj2PUPN/DezmRw9bFUsIIRoNywo4rbUHmAWsAbYDL2uttymlUpVSqf7TLgO2KqW+BNKAydqn1msb/qcIAG/NJt2p8CiFoTUd3NVElLekKvc8Vs1Olq5T0Whd3DeW1mFBeCu6cNtBO8PKK5iTV8SnjqFWRxNCCMsprWu9daxJGjRokM7MzLQ6RsPJySD7n+NYGB0BQGpeFU+U3CLzvYmAkV/qZvCDH9CHb7mtxV9YFhXChfkOzu8QT8h5d0OnwVZHFEKIeqOU2qy1HlTbMVmJoYnK3rqc61f9noXRvpUWwr1e3i2dSJaO57rkOCneRECIDnfw7+uHkKXjeSYymg2hTh5uD9e4vybzvTlWxxNCCMtIAddEpW9+ko3OYACGlldwQ2ERKcYmAGaO6GFlNCFOysC4aOZflECB6yIiTJNSw2BLSDB/opD8UrfV8YQQwhJSwDVFH8xntGsfkabJb0tKufdgNYXlZ5LmmSD3vYmANPnsLuzzjiBx328w/Ld97HfA2JcHyoL3QohmSQq4Jih7899Ji46g0DD4ICyUG6tv4arqedxz/TQSYyOsjifESXM6DN64cRhryicSfSgZtEYrRbkBT34jC94LIZofKeCamOyty5nVJppCwyDSNLm6oJwsHc/No86Q+d5EQOveJpzl0wezq3Dc//7i0pooj4fsx+Igc6l14YQQooFJAdfE/OnzRygyDJymycJDhbxZOhmAaUO7WpxMiNM3tEcMZ8dFEpcXh6E1MdWwO8jG/JbBnPPVQulOFUI0G1LANRHZW5dzxfP92GOYABhAVUV3VnhHy1JZokn5v3G9+dI1k8JvHiH4wDjOLq/iiN2g0DBI2yHdqUKI5kEKuCZiUcZjbDVMtFLYteaG/DLSPBO4O+VM6ToVTUpibARLpg0EYHtFMoYrmTYeE6fpJdrtJnvrcosTCiFE/ZMCLtDlZMA/RjPjyCGSKqtIrHRz1347r+fdwlYVz5VD46xOKESdG5nQjlWzkwEoaP0Fu4IdmGh2OezM/VwWvBdCNH12qwOI07RuAezPZCCw/GAZWWZ3JlTfD8CqWck4HYa1+YSoJ4mxEdxxfg/e+TiFwbxKZqjvNoHD8kdeCNEMSAtcAFv52ULOMb/jucjO7DBjyTK784DnSgwFq2Yny5Qhosm7YnAc2yuSWbv3r/TJj8HQmt8UtKRyyaW+1mkhhGiipIALYGk7XqJQaZ6K1NxiH8+E6vvJ0vE8PaW/FG+iWYgOd5A+tT8AZa5LeOL7KO7KP0zIno/gxfFSxAkhmiwp4AJU9vqHCK12g9aYStEiZjUAMWF2RvRsa3E6IRpOSlIHerQJ4yb7a4wwvqKDLd93oLrcd4uBEEI0QVLABahF25dxIMgOShHmgRJXCi2CDdbcOlLuexPNzoIJfUjzTGCd2YdPzES8GkwN+S16Wh1NCCHqhRRwASZ763IuXdKPrCAbHdzVtKlwcmTfTLZXJLP6luEy35tolgbGRXPP9dO4O2w+fY0fsCkwFPDFMircptXxhBCizkkBF2DSNz/J9zaTSsPGgSA7F5VU4a3owl1jexIb5bQ6nhCWGRgXzXXJ3VhYfTmlOphyHYTL25LsjWusjiaEEHVOCrgAMznhBpym17ehFCujfa8nDexsYSohGofx/Tryb+9oelct4RtvZ3oa+4n9z014F8TJWqlCiCZFCrgAs3W7wd8O5dK9yk2YaeJ0DePmUd2l61QIfKNSP5wz/Cf7OqhcbJUF8MGfLEolhBB1Twq4QJGTwYLFg3hWL2B9aDAv73cx4rsUdhWOY9rQblanE6LR6N4mnLVzhvOgeSXrzD7s97ZiXqtoktpHMC89QaYWEUI0CUprbXWGBjNo0CCdmZlpdYyTlr3+Ie795gW+dwSBUhha023HdLJ0PK+mDpG1ToWoRX6pm5WZe6n45Hqe63gYlAKt2WLEw5WvWR1PCCF+lVJqs9Z6UG3HpAUuACza/iLfBztq/gFqkdeHLB3P8umDpXgT4hiiwx3sOFzC16321vy/E6Q1z0fEWh1NCCFOmxRwAeDKElUzcKGru5oc11SWTx/M0B4xFicTonFrGWKnXW4iQ8orCTWh2mZjsWudTC0ihAh4lhZwSqkUpdQOpdR3Sql5tRyfqpT6yv/YqJTqe9Sx3UqpLUqpbKVU4PWLnqBt+4sILXKy+PARhpVXcK3LoFNkiBRvQpyA6849gxfLruP9PU8Q5komxAPKE0bm0hS5F04IEdAsK+CUUgbwDDAW6AVcoZTq9bPTfgCGa637APcDi392fKTWut+x+ocDWk4G5ou/5U8vT+WabtW8EhLDjAMtWVp6HfeNT7Q6nRABITbKyZJpAwHYVTiO0soelIaUsUjvIXXVVLLXP2RxQiGEODVWtsANBr7TWu/SWruBFcD4o0/QWm/UWhf4Nz8DOjZwRmvkZMCyCTyZv5mdUbsxbfB2pMGE6vsZcm4KIxPaWZ1QiIAxMqEdD/p/6anKPY+oqmC+DnawwRlC+vYXLU4nhBCnxsoCLhbIOWp7n3/fsUwH3jtqWwPvK6U2K6VmHOsipdQMpVSmUirT5XKdVuCG8vj7s+nfIZKlkS1rbr5umZ9EWJCNWaPjrY4nRMCZMLATrcMdeCu6YNorMP3/X/WorJKuVCFEQLKygFO17Kt1ThOl1Eh8Bdzco3YP01oPwNcFe6NS6tzartVaL9ZaD9JaD4qJCYz7xpYZ5XiU7+sxtCYyrw85rqm8NfscWaheiFPgdBi8f8twgu02uhxJBK1BKd5o4aRkzQNWxxNCiJNmZQG3D+h01HZH4MDPT1JK9QGeBcZrrfN+3K+1PuB/PgK8jq9LNuBlb11OaxNsWnNVYTHPf+8gxzWV28fE071NuNXxhAhY0eEOll97NhsLpxFzKJlQUxNiaq72/kD21uVWxxNCiJNiZQG3CeihlOqqlHIAk4G3jj5BKdUZeA24Umv97VH7w5RSLX58DZwPbG2w5PUlJ4O/bbifQ3YYVOFmRF5rHvBciQKmnN3F6nRCBLyBcdFcObgjuwrH0aUyiEMOOzuC4W+f3CddqUKIgGJZAae19gCzgDXAduBlrfU2pVSqUirVf9r/Aa2Av/1supC2wCdKqS+BDGCV1np1A/8IdS5r9W2UYpJUWcXswgKKCSNLx/PO7GRZ61SIOnLr+QkAlLhSSKh0k1RZxQ0FeVR+KCNShRCBw9J54LTW72qt47XW3bXWD/r3pWut0/2vr9VaR/mnCqmZLsQ/crWv/5H447WB7hlVypaQYFp6vXSudJDmmcANw7uRGBthdTQhmozocAerZiezvSIZz+5U/rKvkm+DgjhPf8uKNTdbHU8IIU6IrMTQCGRvXc71zw9gQGEV/co0o/LDudZ9O1/oeK49p7vV8YRochJjI3jy8r5k6XhibMU8FR1JkWHw9P4PrY4mhBAnRAo4q2UuJf2T+Ww0qtkS5mb6wSjuLHmQL3Q8H84ZLl2nQtST8xPb47CBG4Ob8gsJN00iTdiU/U+rowkhxK+SAs5KORmwag6pBYUkVVZRbLORHewr2P50US8ZdSpEPXI6DP41Ywj3VU9jfEklbTwmexw2ZmYtIPvx7pC51OqIQghxTFLAWWndArIdBumRvnvctoQE83WrvQBc2v94cxoLIerCwLhocuImMdn9Jw7b7QBUGYp0p4K1f7Y4nRBCHJsUcBba1/cmFkVGsCHUCcCQ8gra5SayfPpg6ToVooHcfVEvsnQ8oUfOIdiEYFMzrMxDWfI9VkcTQohjkgLOAtlbl5P6wtm8+VE6MwqKGFZewR9yNe/veZKRl6UztEdgrBghRFOQGBvB/IsT2FU4jm6VdqoMxXuhLfnygxep2vWp1fGEEKJWUsBZIH3zk2ygnPSIrdzSNoZrC0pYXHoDY85sLQvVC2GByYO7EOG0++aGK7NzaYGNoXwpc8MJIRotKeAsMDnhhpq1GAvtBk9GxpCl45l3YaLV0YRolpwOg3dvOpftFclk7H2AFaXTeDwsjouCdrPys4VWxxNCiF+QAq6BuT9/jpz1/4ehNWhNC48m1zWO+RclyKhTISwUG+Vkw9xRRDrtZOl4Xo72UmSD+755kcffucbqeEII8RNSwDUw9cF8FkeHYdpshJmKAzsfYVf1uUyWtU6FsFxslJNnp52Fw26jmyuhpqX8xVxZJ1UI0bhIAdeQPpiP3VPCTfmFRJgmoa5hACy/9mycDsPicEIIgMQOEYTYbcwoO8jVhcUYWjOmpJwZz/Yle+tyq+MJIQQgBVzDyclAb3gCBUwsKaPNzuvYVejrOh0YF211OiGEn9Nh8MjEPqR5JjAgrwt/P5DHJ2EhfBrkZeGmR62OJ4QQgBRwDSJ763JmrPoDX/pXWcjxtiZLxxMebJeuUyEaoRE927Aj6Eyuqp7HP6JaUGr4WsjzzErfCipCCGExKeAawCOfP8ynziAWRkdRpJ3c6pkFwJs3DpOuUyEaIafD4K1Z5wBguJKxaQ1AsU3J1CJCiEZBCrj69sF8lKcKgK+DHfzWeSlZOp5HJ/aRUadCNGLd24SzZNpA1pRPpNWhZMI9MCkviOm7hlPhNq2OJ4Ro5qSAq0/++97uzC/ArjWmUhyO+YIop8HFfTtYnU4I8StGJrRjcJdIdhWOw16YxAutveyN3MSn3+daHU0I0cxJAVdPsrcuZ8p705javi0Ad+QWYfMEU+W6gKenDJSuUyECxPxLegNQHL0Frw0KW33FoXcH4/78OYuTCSGaMyng6kn65ifZEmxnS0gwf4uM4Kv8SRTt/DMtqoYyoIuMOhUiUCTGRjD/ogRa5ifVzAu3KCoM+3tzpIgTQlhGCrh6MjE+lW4V0KuymhsKi0gxNhFih7W3jZLWNyECzOSzu5BX+AeSDvUixDSpUIqxse25KvsRmRtOCGEJKeDqWk4GmYuTeXjrX9gdDJWlZ1BYfiZpngmsumk40eEOqxMKIU6S02GwMnUIGwunYWiDKpuNA44gtoQEs2jD/TK1iBCiwUkBV9fWLeBZ8xAuh8Jrg93Re7iqeh7TfjdJRp0KEcASYyNYMm0goa5hhJleOrir6V7lpgST7LdvsDqeEKKZsbSAU0qlKKV2KKW+U0rNq+W4Ukql+Y9/pZQacKLXWsU18FZG5YcT49bYvNAyP4kWDsX5ie2tjiaEOE0jE9rRo/01dNh5Lfe6ismzG777XG2lVkcTQjQzxyzglFLvKqXi6uuDlVIG8AwwFugFXKGU6vWz08YCPfyPGcCik7i2wVW4TX73/iIe7GBSXdyHoh0LyHFN5dYxCXLfmxBNxJ1je5Gl43k2KpxCwyDSNGnnOpP9BRVWRxNCNITMpfBIV9+zhY7XArcUeF8pdY9SKqgePnsw8J3WepfW2g2sAMb/7JzxwIva5zMgUinV/gSvbXCbPn6Pgsiv8Np8Uw4AKODS/rHWBhNC1JnubcK5O6UndlcyQ8oreeqwi1HuPF76bLfV0YQQDSB73b1MiQxiStYjZB/JtizHMQs4rfXLQH+gJZCplLpdKTXnx0cdfHYskHPU9j7/vhM550SubVAVbhPb+oX8obAYQ2ta5icxsHMkm/84RgYuCNHEXDm0KxvN35G79wa+9caSFlvA5zmXs3HzMqujCSHq2d9ahrIlJJgtwXae+WKRZTl+7R64aqAMCAZa/OxxulQt+/QJnnMi1/reQKkZSqlMpVSmy+U6yYgn7ntXKU+ZExiQ14VuO6aT45rKwsv6SvEmRBPkdBisvmU4WTqeV6IUO52Kb4M1M7c8wiXLhlr6W7kQop5kLmXzI50o1R66VblJrHSTveUsy5bWO949cClANhAKDNBaz9da//nHRx189j6g01HbHYEDJ3jOiVwLgNZ6sdZ6kNZ6UExMzGmHPpbuMeEcaNmH6/VdfOvoxarZyTLqVIgmLDbKyQPjEylxpRBmegHwKsUP3hLSv0y3OJ0Qos69/0f+0TKYLSHBtDdN/rzPID+vPd+7rBnEZD/OsXuASVrrbfX02ZuAHkqprsB+YDIw5WfnvAXMUkqtAM4GirTWB5VSrhO4tkE5HQbv3zKc712ldI8Jl0ELQjQDFyZ14I9vJpOQA9Ht38Rlhw4eL9N3bIIzM6DTYKsjCiHqQk4GX9g8FNtsJFVWcWG+g7s81+EMstE9xprGmmMWcFrrc+rzg7XWHqXULGANYADPa623KaVS/cfTgXeBC4HvgHLg6uNdW595T4TTYdA7NsLqGEKIBhId7mDtnOGMfhzYlUxW8Az2hrhJj/RgrL6NAdd9bHVEIUQdqPzwIf4eEc6WkGCGlFdwU8mjAKyddY5lDTbHa4Grd1rrd/EVaUfvSz/qtQZuPNFrhRCioXVvE84d5/fg0fd38pHZl/cjv2VjqBN3VSHPWx1OCHHaKtwm1+0exezgzwG4NF/xPvBq6hBLb5WSlRiEEOI0XTE4DoCRxpfMLCwiqbKKXKq5YnEC2esfsjacEOK0/PvDv1PefjF/bdWSC/MdLC69gT6xEQyMi7Y0lxRwQghxmqLDHbyaOoSF1ZfTtdIgzKv5ITiIrcF2Fm1/0ep4QohTkZOB+5Ez+HTP4+x0KraFOHgj2iRLxzN/nOVrB1jbhSqEEE3FwLhoiq+8k74vjGZh/j2U2coBuLpMBjQJEYiy3pjO4pZwXlk5JTZfe1eQK5kbhne1vPUNpAVOCCHqzMiEdtx+fjwrSqex6EAxyw8epqyiiusX9SB71Syr4wkhTsLf7RVsCHXyYVgoLx04zMScM1hTPpFrzznD6miAtMAJIUSdmjK4C4+9H8+dxoX80H4jB4PsoBTsXcPfrQ4nhDghxW/fzcyCIhQwKd+ka9VyANbOGd5oJuiXFjghhKhD0eEOxia2oaD1Fxx0BPmKN62ZWiJdqUIEgpWfLeQi1xvsdASRfthFVGUUQKOboF8KOCGEqGO3X5BAiSuF9u5qbFrT2mNyX5SN3yxJYuVnC62OJ4Q4jke/+SeFhsGj0VHkmmE84LmSlF5tSGxk87xKASeEEHWse5twZl5yI+G7rufsiipyg+wcdhiU2SBtx0tWxxNCHEPxtxswzR+XVtfYbIosHc8dKQmW5qqNFHBCCFEPLunfkZgzk1Guc0iqrCK2uhqnadICpyx2L0Rjk5PByueGcuEnMxhTXkoLj+aG/EoWVl/OFWd1bFRdpz+SAk4IIerJ/13SmzXlE+m852LezTlIB49JDmXc++YVkLnU6nhCiB+tvovHVSFFho3/hjn5/a6+3Jv7FCu8o5k1Kt7qdLWSAk4IIepJbJST5dMHs8I7mp06lsN230CGA4bi+k0PyCoNQjQSVR4vbTwmAG08JpPt6wDfclmxUU4Lkx2bFHBCCFGPhvaI4brkOO6qvo6peRBhmrT3mGwMdZL+zTKr4wnR7FW4TVYfCuPe3HyGlldwd24xC6sv566xZzaKCXuPRQo4IYSoZzNH9CBLx/NY3qMk7vsNAE7T5GCQQ+6HE8Ji/1k6gXfbfYNS8PShPC4vfpYV3tFMGtjJ6mjHJQWcEELUs+hwB/Mv9o1iq47ZwK5gBxWGwS6bl/Qv0y1OJ0QzlZNBxRODeMvc5msRj4zggPa1uC2fPrjRTNh7LFLACSFEAxjfryMKCHINI7GyigiPB7RGHz5odTQhmqWs1bdxa3Ah55WVM7S8gkn5JrdWz+K65DiG9oixOt6vkgJOCCEaQHS4g3dmJ7OmfCKf/fBXSg0DlOJzz/e4l14KORlWRxSiWXmaEjaEOrmvdTSqsg3TS54mS8czc0QPq6OdECnghBCigSTGRvDYxCQAuuTFYdOaNh4PXx/cAKvvsjidEM2H6+uPGVDoBq1BKT6LKsYGbJg7qtF3nf5ICjghhGhAl53VmWHdo/jSNZPEchsHg4JIj4ygsqzI6mhCNBvPfHQTf4+x1axV3CUvjjljejTaKUNqIwWcEEI0sPvG9wHg0gIbw8orSC0soqpgH/u+WmdtMCGaupwMPk8fxqth1TXFW8c94/jSNZMpZ8dZne6kSAEnhBANrHubcN8Ev6XTGJUfziPRUczs0JLD70yibMOzVscToulat4AnlaumeIvM68P2imTWzhkeMF2nP5ICTgghLDC0Rwyjx1zE/MhubA0JZktIMFd3aMvIHY+z8rOFVscToknalTgLDwqAbpWKHNdUHp2Y1CjXOv01UsAJIYRFrhgcR1XueYSaGgCtFBWGwZPf/NPiZEI0PR+tvpcHPp9JTGFXEsrsuA9fjNMOF/eNtTraKbGkgFNKRSulPlBK7fQ/R9VyTiel1EdKqe1KqW1KqZuPOnavUmq/Uirb/7iwYX8CIYQ4fdHhDv45ZRLROeOIq6om2OvFaZrckFuAuXiUTC0iRB1Zl/ECdx94mU2hDnIjd5Ox9wG2VyTzysxknA7D6ninxKoWuHnAWq11D2Ctf/vnPMBtWusE4DfAjUqpXkcd/6vWup//8W79RxZCiLo3tEcM4W0uIuiH61m7O4+MvfuZUlqMcWAzLL9cijgh6sAz2Y/55l4EqrUduw1WzU4mMTbC4mSnzqoCbjzwgv/1C8ClPz9Ba31Qa53lf10CbAcCs51TCCGO495LksjS8VztnkuW2Z09ZmuygoNJbWmQ/fYNVscTImBlr3+I8f9I4FuH7zaFcNOk4vA47rigZ0AXb2BdAddWa30QfIUa0OZ4Jyul4oD+wOdH7Z6llPpKKfV8bV2wQggRKBJjI0if2p8sHc8DniuJspXxWHQkG0KdzAwpJ3vrcqsjChGQFm1fxi6HHa9SGFrTZ99v2F6RzKSBna2OdtrqrYBTSn2olNpay2P8Sb5POPAqcIvWuti/exHQHegHHAT+cpzrZyilMpVSmS6X69R+GCGEqGcpSR24a2xPbrK/RktVUbO/1DBI3/SYhcmECFyXVLalW5WbMNNL9CHfUnaBsFD9iai3Ak5rfZ7WunctjzeBw0qp9gD+5yO1vYdSKghf8faS1vq1o977sNba1Fp7gX8Ag4+TY7HWepDWelBMTONfnFYI0XxNGtiZp8wJZJnduTbXQ1JlFUmVVVzrypN74YQ4SZt351NSvZN8u8Et+YXsKhzHg+N7B8RC9SfCbtHnvgVMAxb4n9/8+QlKKQU8B2zXWj/+s2Ptf+yCBX4LbK3fuEIIUf+iwx0svOU6Rj8eD9UwuWwt9wUtxaFMzI8exvjD61ZHFKLRy966nPRNj9Gu2M2rraNAKdKiorDnw4SBHa2OV2esugduATBGKbUTGOPfRinVQSn144jSYcCVwKhapgtZqJTaopT6ChgJ3NrA+YUQol50bxPO3WN7ApBibMKhTNza4OuIcy1OJkRgSN/8JBts1bwaQc2KCyGuZH7bt0PAThlSG0sKOK11ntZ6tNa6h/8537//gNb6Qv/rT7TWSmvd5+fThWitr9RaJ/mPXXJUa5wQQgS8K4d0JTTIRppnAvk6nDdbhDAz7x+8/OwQ6UoV4lcMjhqMXWtQCrvWXO0y2FU4jlvOP9PqaHVKVmIQQohGxukweHv2OWTpeK51386T0ZEUGQb320t4ec3Nv/4GQjRTGzcv4/nctXiUIsI0+fN+L5/kX8Xy6YOJjXJaHa9OSQEnhBCNUPc24SyaOoAsHc+EvCDwtyjcby/hrGUDWbljpdURhWh00rOeoMgwiDRNbjro4KaSRzln5NgmM3DhaFLACSFEIzWiZxvCHAaf5F/F1S6jpoir9LpJ+yLN6nhCNB45Gbz0TD+2GJWgNb3KDFaUTgNg2tBuFoerH1LACSFEI+V0GLx8/RCydDxpeQ8R4vXNJo/W/CGkj7XhhGhEyl+5gXSnG49hA6XYEG4jS8ezanZyk5jzrTZSwAkhRCOWGBvB8um+qS57HEkk0jT5v9x8fvfFq+SXui1OJ4T19hdU4C3cw035hThML2hw5yfz4PjEgF8u63ikgBNCiEZuaI8YLkhow8bCaTz7A0wqLeNTRwvmrBgqy2yJZi17/UPM/3d/MhwtmFRaxvLdmpJvFhBWdgkTBnayOl69kgJOCCECwJ1jEwC4y3Md68w+vBjlZHNQFembn7Q4mRAWyclg0bYlfO4MYn77EN4KiuWu6uto6YC1c0Y2qTnfaiMFnBBCBIDubcJ5dGISWTqeq6rnUeC6iMGVXlIP7oEP5lsdT4gG5177EDMLC4k0TQoNgzeiTbbb4ll355gme9/b0aSAE0KIAHFx31g6RIZgt8H2imT+fmA//Sor8G6QEamieck+ks1l3iM8Eh3F9LxKzi6vIsg1jHduHt4sijeQAk4IIQKG02Gwds4Ibh4dD8AXwSFcENueAXGx/OWpMyBzqbUBhWggj61/gB+MEraGBPNeSwcf7vkrfZLn0b1NuNXRGowUcEIIEUCcDoOpZ3cBIC2yNQccQZhK8c8WIXjfmSNLbYkmr8JtcsmO7wg3TQCqtR2AKwbHWZiq4UkBJ4QQASY63MEVZ3Uk1zWO1m6wac0FJWXc0Daa7JcukZY40WRlr5rF9CUJvNa6nDn5hQwur6Li8DgeHJ/YbLpOf2S3OoAQQoiTN2tUPMM2JcP3yQxQ3/JDXDpfhzg4aBi03/QgqSEO+vWeYnVMIerUor1r2BIaAkCoV1G4N5X9QfFNfsqQ2kgLnBBCBKDYKCerZicDkKXjqdZBAByxG2wIDZHpRUTTkpPBw+kD2egMBq2JrvZw8MgEsonnvVtHNfkpQ2ojBZwQQgSoxNgINswdRbChmJxnJ6myCkNr0JqWthCr4wlRZ6rfncfLIVWgFChFkd1ge0Uyz/5hELFRTqvjWUIKOCGECGCxUU6mDO7MitJphHqhyG4HpXjPzCP75clWxxPitGWvf4hZ3r1cUFIGWqO0Jik/htAgG0O6t7Y6nmWkgBNCiAB37bndydLx7D8ykQ7uatAalGJ26VdkH8m2Op4Qpy4ng4Xbl7Ix1MleRxAzd/al+JtH+PjI7bw9+5xm2XX6IynghBAiwMVGOVkybSDbK5Lp/8OFvHjgcM3s9LPenSbrpYqAlfnaNfwQ9L/xlpPt6wBYNTu5Wc35Vhsp4IQQogkYmdCOIXFRrPCOpnuVwW+LS0FripSXRRmPWR1PiJNW4TZZ7Kik1DCINE3m5BWxsPpyXk0dQmJshNXxLCcFnBBCNBEPTOgDwEPVU3itZbjvhm+t2Whz0/+FvqzcsdLihEKcoJwMvE8kkVpQxNDyCu49WMlDJfPw9J3KwLhoq9M1ClLACSFEE9G9TTjpU/uzwjua6/PLsPnvhUMpPHhJ+0LWTBUBIHMpeslYwsr3M8Dt5u+HXRiVsWTpeG49v5fV6RoNKeCEEKIJSUnqwN0pPdmWP5FQr9e3U2vwQnFlAY+/c421AYU4ng/mo9+5GeX14NHg0Yo93hjSPBNYMm1gs50ypDaWFHBKqWil1AdKqZ3+56hjnLdbKbVFKZWtlMo82euFEKI5unJoV1Z6RzPA1ZkI0+Tcwx0xlMarFMtyZa1U0XjpjU+j/K9dOoIzql5iuPtJep89mpEJ7SzN1thY1QI3D1irte4BrPVvH8tIrXU/rfWgU7xeCCGaFafD4OXUIawqmM2+bx9lVcFsrigsw641U4vKrY4nxDG9GTOU69vGkB3sIFS5a/ZfP7yHhakaJ6sKuPHAC/7XLwCXNvD1QgjRpA2Mi2bJtIE124XuToR4vbwRHkLfpb2lK1U0KivfuJIhz/fiwZDv2Rjq5G+RkTxcPYWBnXyrjUjX6S9ZVcC11VofBPA/tznGeRp4Xym1WSk14xSuF0KIZmtkQjvWzhkOwPrWeZQaBkV2A69SLMnNkPnhRKORlr+ZUsOg0jAIM01yjlzGCu9o0qbIfW/HUm8FnFLqQ6XU1loe40/ibYZprQcAY4EblVLnnkKOGUqpTKVUpsvlOtnLhRAioHVvE87gLpF0OZLoG8wANSs1LNr8hKXZhABwff0xU/M8vnV8gXZu3zqnN4/qLsXbcdh//ZRTo7U+71jHlFKHlVLttdYHlVLtgSPHeI8D/ucjSqnXgcHAeuCErvdfuxhYDDBo0CB96j+REEIEpvmX9Oaip6bR3rmY0ojvGVpegVaKcyPHWh1NNGc5GbBuAZUHD5Bafpik6hY8ExlNgesiAKYN7WZxwMbNqi7Ut4Bp/tfTgDd/foJSKkwp1eLH18D5wNYTvV4IIYRPYmwE916UwMGDM+ixYzoPH65k8WEX1ftfIfmFAaxM6wGZS62OKZqTnAx4YRx8v5bDnt3MaBtDmKpixsEItlcks3bOcKLDHVanbNSsKuAWAGOUUjuBMf5tlFIdlFLv+s9pC3yilPoSyABWaa1XH+96IYQQtbvqnG7MSI4jS8dzrft28nU4z0WHUEQ1aeFBsPbPVkcUzcmr14KnEoBno8L5NNRJWmQ0aZ4Jss7pCaq3LtTj0VrnAaNr2X8AuND/ehfQ92SuF0IIcWy3nn8mK7P2k1XuK+JS8v7F69FubiypgNHzrY4nmhFv4T5s+G7HvCK/mmIdykHXBFJ/N0nWOT1BshKDEEI0E06Hwdo5I2gVHkSWjucfefMZuH8AHwbb+OzTJ33dWkLUs8y3bmRm22iyHQ7fcr2Vndm492FcxkjOT2xvdbyAIQWcEEI0I9HhDj65czQ3jToDgOqYDXweGsxSewne51PkXjhR7+4/9CEbQ53Mbx1Nvg4nzTOBiBCD928ZgdNhWB0vYEgBJ4QQzYzTYXDV0K4ABLmGkVRZxUHD4PftWvPFf/5kcTrRVGUfyWbsC2ezyxEEwEG7nWvdt7PVFs9Ht4+SQQsnSQo4IYRohqLDHSyfPpg15RNp4fWyK9jBlpBgbmjdgs8PZFkdTzRBD63/E/soB+Vb7TTaYyNLx/P6DclSvJ0CKeCEEKKZGtojhtgIB6Pywgk3TQBKDYP7Pn7C2mCiSZm3+lqSliaxvfQHAJTWdKoyUQcv5sbhXWXQwimSAk4IIZqxBy5NYm7pg/TZ9xt6VbrpVOUluCCL5X8/m9S3J5N9JNvqiCLArTr8OShqWt66Viq+3vUoP7iTmTW6p7XhApgUcEII0YyNTGjHHWN6sKZ8IjP3taGzWcVOp+LR4DI25G8jfdU1MjpVnLLHP7rDt6q51hheLx2rNO7DFwOw6ubhMmjhNFgyD5wQQojG48bR8XRvE07a8gncVvAXtgU7KDQMIk2T1MOHYPVdcN1aq2OKALPy/TksOfA+KIXSisIdj1DoP7Z8+mCZrPc0SQucEEIIUpI6cP+sa/hLyW1cczCCTmXhzMovIj0ygk1VhVbHE4EkJ4PsJ87kAX/xBmBzh9UcXjJtIEN7xFiVrsmQAk4IIQTgWzP1uimX8+eSP3LnQTsfhTnZEOpkTmgF61bOsDqeCBTrFpAe7MHrL96U1hQf/AMAN47oysiEdlamazKkgBNCCFFjRM+2RDkN0jwTiKwMBa0pNAwW568jv9RtdTwRAFa06cSXwQ5aV3sIM720PpSMt6ILANOTz7A4XdMhBZwQQogaTofBOzcN5wsdz9stnTVdYGVKMW9ZL7KfSpJBDeKYKtwmT+7/D6WGQaWyc+jbhewqHEff2JZk/XGMzPdWh6SAE0II8ROxUU7emZ1MlesCHB6IdVdzxG7waaiTdHsFrFtgdUTRSK1dOoE2HjfhpsnUfLNm/99+P0iKtzomBZwQQohfSIyN4PELZpC3cwHR1Q5KDYNw06TYZuMDe7TV8URjk5OB5+E43ja3sSvYQd8qN0OLQwDfoIXYKKfFAZsemUZECCFErS7p3xGARW9dRAKrUYabLSHBPF+8niV/jweluPPMafQ7926Lkwqr/evt6/lb21B+W1wKwO/yvTzguZK0y/vKoIV6Ii1wQgghjumS/h255sIbyNj7AFNzDZIqq/g+yM6WkGC2BDtI/2aZ1RGFxfYXVPBMSBWFhsHrLcOZvL8d15Q8xcBh59f8EiDqnhRwQgghjuvivh2IchosLb2OEK+iwvDNnh9umlzVbYrF6YRVstc/xKX/OJPfvj6QIK0JM006HenN1Z67AJg5oofFCZs2KeCEEEIcl9NhsPa2UXxj78nBIxM4s8JDUmUVlx9oz2t5F1odT1ghcymLtj3P944gKgyDXLudvlVuNhZOA+DV1CEyaKGeSQEnhBDiV0WHO/hgzki2VyTz94PF3JlfwNet9tDm4G8Z/Hwig5/vxco3rrQ6pmggWf/5EyU2G7HuapymSfcqNwl5vrneXk0dwsA4GehS35TW2uoMDWbQoEE6MzPT6hhCCBGwVm85wLp/PUZ+p7f5NDQEu9Z4/HPF2bVmSbcpMqihCcveupxbMx4m16ZBKYaVVzD5QHuuqp4H+NY4lWWy6o5SarPWelBtx6QFTgghxAlLSepAn3E3kXPkMvqVaaYWFuM0veAv5GbtXAYfzLc6pqgHKz9byLTMh8g18E3wrDXXFRST5pkAyBqnDU0KOCGEECdlypA4/jRlDh/vfYSnDy1i2Q//68kpMgz49GkL04n68uQ3L9Wsb4rWXFBSycrS35Ol43k1dYhMF9LALCnglFLRSqkPlFI7/c9RtZzTUymVfdSjWCl1i//YvUqp/Ucdk7tohRCiAQ3tEUP61P4AtLPlM7S8ArRmSHklZUYk/GO0LLnVhBz5KJ3rcwtxml6cpsnlR5y8sv9JVnhHs2TaQLnnzQJWtcDNA9ZqrXsAa/3bP6G13qG17qe17gcMBMqB14865a8/Htdav9sQoYUQQvxPSlIH7hjTg4erp/Dw4Uou+vZ8JhQobotWZOd+JUtuNREbNy/j5m8f570WwfztUC4dd17Hs/n3AshEvRayaiWG8cAI/+sXgHXA3OOcPxr4Xmu9p35jCSGEOBk3jo7nUfe1DPjvaAD2tH6HrSEhFNts3NRpOr+xOJ84DTkZrHz/Fh4wivGGBAPwbFQLbqp4jauq53HvRQkyUa+FrGqBa6u1Pgjgf27zK+dPBv71s32zlFJfKaWer60LVgghRMO4Y2wvHhqfCIBb+9oFtgY7uPGHO7nkuSSyty63Mp44VesWkKYKa+57M7RmVH44aZ4JpF3el6vO6WZxwOat3qYRUUp9CNTWrnoP8ILWOvKocwu01rUWYUopB3AASNRaH/bvawvkAhq4H2ivtb7mGNfPAGYAdO7ceeCePdKIJ4QQ9WHxuu94/b/PcbDL25g/3uwOoKGDMnhk4J306y0rNwSEzKVUr76H54LDeK5VEAo4y9WJVQWzWTtnON3bhFudsFk43jQilswDp5TaAYzQWh9USrUH1mmtex7j3PHAjVrr849xPA54R2vd+9c+V+aBE0KI+lPhNvnNg+/TyvkmxTGfoG0at1K+KSeAWA+snr7F4pTi16x840oez99MiFdTaVN0PdK7ZoWF9Kn9SUnqYHHC5qMxzgP3FjDN/3oa8OZxzr2Cn3Wf+ou+H/0W2Fqn6YQQQpw0p8PgoztG4yq/hNydj5C3YyFVecPB31Cw39BkvnWjxSlFbbKPZJP69mReTu/PQ4VfUGoY5AbZKTUM9rTZBvgGLEjx1nhYVcAtAMYopXYCY/zbKKU6KKVqRpQqpUL9x1/72fULlVJblFJfASOBWxsmthBCiOOJDnfw3zvPo1Wo7144t+tC+uXH+Io4pbg2/7+89NwIa0OKX0j/Mp0N+dt4ylGFRykMrYmu9hBmeulyJJEl0wbKgIVGRpbSEkIIUecq3CbvbzvEzf/OBuA3XW9lm38kY6RpsvTiDLmPqhFZvyeT+9dcQ77dxEBzlqszqwpmA9JtaqXG2IUqhBCiCXM6DMb3j+WxiUkARBxOpluVG6dpYvdqprz9G5Z9LPPEWS17/UOkPJvAjR9dxSGHF7dho8IwcLU4DEjx1phJASeEEKLeXHZWZ+5O6cma8ol8uetxeleZ5AbZKbcrFn7/T+asnGh1xGbtb9tfZH+Q3TfQxL++aZcqDyWuFCneGjmrJvIVQgjRTMwYcQbtIkK46d9fEuQaRkjHz6k0bKAUH5Tt4PpFZzCzwku/SSug02Cr4zZ5K18YyQLzMG6bDYLt4PWCUiigb34MHx+5XQYsBAAp4IQQQtS7S/p3JDYqlInpYMsZRFTsc7jtVbT2mGwMdaKoIH3dArjy52PWRF1LMw/jNoz/7VCKIeWVvL/nCT4GHhqfKAMWAoB0oQohhGgQA+OiWTU7GW9FF/K+u4+ZO/vylyO5JFVWsdURRJL5LUlLe3OBrN5QP3IyyH72XKJMjd30+kYGa02Ex0O7XN9KGulT+zNlSJy1OcUJkVGoQgghGtTGnS6mPJdRsz2qyxw2hTp+ck64abIovA/9freioeM1We6ll3JTxVY2hDpJKLOTsfeBnxxfNTuZxNgIi9KJ2sgoVCGEEI3G0B4xrJqdjMP/L1CQaxitqj01LUIApYbBVeVbufTZRGmNOx05GWT/9Uym/D2eaVVfM6qsggFlJiWulJpTFLB2znAp3gKMtMAJIYSwRIXb5PmPv+fRD3bW7EtwfkJ42zfYGWLUrKc6jFDSp31uVczAtmwCqaVfsSHUCcCw8goeOlTBgKrFAFzSpx33XpJEdLjjeO8iLCItcEIIIRodp8PgxtHxpE/tX7Nve0Uy5p4ZzHZVEm6axLqrOVRVyPjFPZmy7DdkH8m2LnAgyVzKyr924Zzqb4ivcpNUWUXvyiqmFVSwsPpywHe/W9qUgVK8BShpgRNCCGG5t77Yx03//hKApUELGGF8RZF2MrddeE3rEUBrI5SPfi+tcceSvf4hbt25jFzDQAFaKSJNk/V79tO1ytcVHe6wsebWEcRGOY/7XsJ6x2uBk2lEhBBCWO6S/h1pHR7MlOcySPNMAKA9eaQWuii22dgS7AClyPWUWZy0cVv0zTJyg33/tGvArjU35ReyzvStiBHlNFh72yhpdWsCpIATQgjRKAztEUPWH8fw8Kr2XPVFPAPUt9xU/hqVxRPo0O5dDkQcpFNJe/JL3VKAHGXlU2fyeJiNNh6TM6pMcPjmeGttmpxzoCd3lF0HwC2jzuD6EWfgdBjHezsRIKQLVQghRKOzessBUl/64pjHZZmn/znn+V4U+ifmtWuNRylaml72f7uw5pzHJiZx2VmdrYooTpF0oQohhAgoKUkdWDU7jIue+qRmX9+YRexqtRuAeZs0X68tYUewg5mVdvpd9mLzWYYrJ4Pst2/gXkcZh4Ic9KmsIiPUSdtqD8mlXtZE2Ol8pDf7/acvmTaQkQntLI0s6p60wAkhhGi08kvdPPbeNpZvPkDkmXNrphaB/7U2dXBXcyDI1x7R2qv46+C76Nd7ilWR69/Ppgb58Xs4q9zNf/Y8XnPamDNbc+/4PjJYIYBJC5wQQoiAFB3u4KFJ/fnNGTE895+4mhY4w6u5stjXAve5MwT8hV2uAembnyS9qRRwmUvJ/uj/WNgyhL12g6Ifu0qDg2lV7aFE2bm8pIxvg4MwXMm+YwpevGYwQ3vEWJlc1DNpgRNCCBEQvj9Syt2vfcnnuwtr9v3F/gx7Wm9lSWRLAAytaevx0grFnUXl9BtxLwy6ypK8deKRrqS2NH4ylcqPfn6fG8A1Qzpzx9heMlChiTheC5wUcEIIIQLK/oIKHlu9jde/PPyT/YM7/5HtYZ6a7WHlFaQXmzD3h4aOWGfKNjzL1x/P5/GoMHKOaoFzeDU9jySysXAaAC1DbPzfRYkyUKGJkS5UIYQQTUZslJO/XjGISYNcTHkuo2Z/iSuFHra38dg8BHsVVxZUsjjkatxrv2XK2XGBMfXIq9eR/e0bPNKqFQU2KLTZmGU36LLnYj71jv7JqRuBcIdiza0j5T63Zkha4IQQQgSs/FI3T7z/DS9m5PzquWN7t+X288+ke5vwBkh2avS9EcxsG/OTLtNI0+TNPUU165f+6O6Unlw5tKt0lzZh0gInhBCiSYoOd3DfhD6MH9CRe9/aypYDJTXHLop6iq9i9vLb4lJWh4eysdrOw69Usi3EQVuPydTiEp6LbMkBux2lNSjFtC5jmTPy0Qb9GTbvzue+t7ZR7jZJ9Qzj+oJNFCiDQruiyGbj+vzSmvVLAQZ3juThy/o26kJU1D9pgRNCCNFk5Je6eWXzXh5+bwex8XdQZBg102wA4C/UwNey9eMEuD+yA19M21KvGSvcJtsOFLHjYDHPfbyLXfkVJ3RdILQgirrV6AYxKKUmAfcCCcBgrXWtVZVSKgV4EjCAZ7XWC/z7o4F/A3HAbuB3WuuCX/tcKeCEEKJ5yC918+yKybxt+6amBe6wYeesimO3wHn9hZ3dq0lwu9kSEozT66XCZgN8S1O195j8tqSUl1q24IDdIExDsVIYaO7IL6SHqZgfFc6+IDtoaKl1zfFRZRWsCXPisdmIME2UtlFsVzhMk0rDoLXH5PrCYlaFRpHrGsf2imScQTZuHNmdKYMD5B4+UacaYwGXAHiBvwO311bAKaUM4FtgDLAP2ARcobX+Wim1EMjXWi9QSs0DorTWc3/tc6WAE0KI5qXCbfLp97m8t+UgK7P2H/O8pUELmN29oKZ17uiWup+rreXux/2JVe5ap/z4tff8+Xv3dzu4fsQaBnSJknvcmrHjFXC2hg4DoLXerrXe8SunDQa+01rv0lq7gRXAeP+x8cAL/tcvAJfWS1AhhBABzekwGJXQlkd/14+sP47h9jE9SGz3yy7INM8EzinxgtbYTS9JlVWgNU7T9BVeWtO62kPvyipm5xfSvcpNiGnSqtrEYXpxmiY35RdyVX4p3arcOEwvDtP7k+MXlZRh979fhMdDZLUXm9aEeDw17z8rv5gh1TbmDLmNYT1aS/EmjsnSe+CUUus4dgvcZUCK1vpa//aVwNla61lKqUKtdeRR5xZoraOO8RkzgBkAnTt3Hrhnz566/0GEEEIElG37i3ho1dfkFJSxt6DKshw9WocRHe7g/8YlkhgbYVkO0ThZMgpVKfUhUNvqufdord88kbeoZd9JV5ta68XAYvB1oZ7s9UIIIZqexNgIXpoxBICPth/iT29uIcwRRKXHZE9+Zb19bstgG5FhdhLaRnDn2F4yIEGcsnor4LTW553mW+wDOh213RE44H99WCnVXmt9UCnVHjhymp8lhBCimRqZ0I5PEv7X3pBf6ubp/3zDO18dIiLEzqHiCsKCg+jfOZrI0CA8HpPCimpsykablg4+/S6f/PIqEtq3ZHh8Gz7+LpfwEIPiCg+dI0P45nAZ5/aMYdqQrjIQQdSZxtyFasc3iGE0sB/fIIYpWuttSqlHgbyjBjFEa63v/LXPk0EMQgghhAgUjW4Qg1Lqt0qpfcAQYJVSao1/fwel1LsAWmsPMAtYA2wHXtZab/O/xQJgjFJqJ75Rqgsa+mcQQgghhLCKTOQrhBBCCNEINboWOCGEEEIIceqkgBNCCCGECDBSwAkhhBBCBBgp4IQQQgghAowUcEIIIYQQAUYKOCGEEEKIACMFnBBCCCFEgGlW88AppVxAfa9m3xrIrefPaOrkOzx98h2ePvkOT498f6dPvsPTF+jfYRetdUxtB5pVAdcQlFKZx5p0T5wY+Q5Pn3yHp0++w9Mj39/pk+/w9DXl71C6UIUQQgghAowUcEIIIYQQAUYKuLq32OoATYB8h6dPvsPTJ9/h6ZHv7/TJd3j6mux3KPfACSGEEEIEGGmBE0IIIYQIMFLA1SGlVIpSaodS6jul1Dyr8wQapdTzSqkjSqmtVmcJREqpTkqpj5RS25VS25RSN1udKdAopUKUUhlKqS/93+Gfrc4UqJRShlLqC6XUO1ZnCURKqd1KqS1KqWylVKbVeQKRUipSKfWKUuob/9+LQ6zOVJekC7WOKKUM4FtgDLAP2ARcobX+2tJgAUQpdS5QCryote5tdZ5Ao5RqD7TXWmcppVoAm4FL5c/giVNKKSBMa12qlAoCPgFu1lp/ZnG0gKOUmgMMAlpqrS+2Ok+gUUrtBgZprQN5DjNLKaVeAD7WWj+rlHIAoVrrQotj1Rlpgas7g4HvtNa7tNZuYAUw3uJMAUVrvR7ItzpHoNJaH9RaZ/lflwDbgVhrUwUW7VPq3wzyP+S33JOklOoIXAQ8a3UW0TwppVoC5wLPAWit3U2peAMp4OpSLJBz1PY+5B9PYRGlVBzQH/jc4igBx9/1lw0cAT7QWst3ePKeAO4EvBbnCGQaeF8ptVkpNcPqMAGoG+AClvi78p9VSoVZHaouSQFXd1Qt++Q3d9HglFLhwKvALVrrYqvzBBqttam17gd0BAYrpaQ7/yQopS4GjmitN1udJcAN01oPAMYCN/pvMREnzg4MABZprfsDZUCTujddCri6sw/odNR2R+CARVlEM+W/b+tV4CWt9WtW5wlk/u6WdUCKtUkCzjDgEv89XCuAUUqpf1obKfBorQ/4n48Ar+O7TUecuH3AvqNa0F/BV9A1GVLA1Z1NQA+lVFf/zZKTgbcsziSaEf8N+M8B27XWj1udJxAppWKUUpH+107gPOAbS0MFGK31XVrrjlrrOHx/D/5Ha/17i2MFFKVUmH8gEv5uv/MBGZ1/ErTWh4AcpVRP/67RQJMa0GW3OkBTobX2KKVmAWsAA3hea73N4lgBRSn1L2AE0FoptQ+Yr7V+ztpUAWUYcCWwxX8PF8DdWut3rYsUcNoDL/hHlduAl7XWMg2GaGhtgdd9v5NhB5ZrrVdbGykgzQZe8jeq7AKutjhPnZJpRIQQQgghAox0oQohhBBCBBgp4IQQQgghAowUcEIIIYQQAUYKOCGEEEKIACMFnBBCCCFEgJFpRIQQ4hQopVoBa/2b7QAT39I9AIP9ayILIUS9kGlEhBDiNCml7gVKtdaPWZ1FCNE8SBeqEEIIIUSAkQJOCCGEECLASAEnhBBCCBFgpIATQgghhAgwUsAJIYQQQgQYKeCEEEIIIQKMTCMihBBCCBFgpAVOCCGEECLASAEnhBBCCBFgpIATQgghhAgwUsAJIYQQQgQYKeCEEEIIIQKMFHBCCCGEEAFGCjghhBBCiAAjBZwQQgghRID5f6qoFtZFm/n7AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.figure(figsize = (10, 5))\n", "plt.scatter(t, Y0, label = \"Y0\", s = 6)\n", "plt.scatter(xtr[:, 0] ,trPred, label = \"Train\", s = 2)\n", "plt.scatter(xte[:, 0] ,tePred, label = \"Test\", s = 2)\n", "plt.xlabel(\"T\")\n", "plt.ylabel(\"Y\")\n", "plt.legend()" ] }, { "cell_type": "code", "execution_count": null, "id": "d3e6df7c", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.9.12" } }, "nbformat": 4, "nbformat_minor": 5 }