Skip to content

Instantly share code, notes, and snippets.

@kspicer80
Last active May 31, 2022 18:30
Show Gist options
  • Save kspicer80/5faed2dd26d1cd09df75fc2901bdfda3 to your computer and use it in GitHub Desktop.
Save kspicer80/5faed2dd26d1cd09df75fc2901bdfda3 to your computer and use it in GitHub Desktop.
date metric_tons
1870 13
1871 12.9
1872 15.4
1873 17
1874 16.4
1875 16.5
1876 16.6
1877 16.3
1878 16.4
1879 16.6
1880 18.8
1881 19.2
1882 20
1883 20.8
1884 19.5
1885 19.1
1886 19.5
1887 20.8
1888 22.2
1889 23.9
1890 25.6
1891 25.5
1892 25.7
1893 25.2
1894 27
1895 27.6
1896 28.8
1897 30.3
1898 31.8
1899 32.3
1900 32.7
1901 31.6
1902 29.4
1903 34.2
1904 33.5
1905 35.2
1906 33.5
1907 36
1908 36.6
1909 37.1
1910 37.6
1911 38.5
1912 40.4
1913 40.1
1914 26.8
1915 18.9
1916 20.5
1917 27.8
1918 24.9
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"id": "c928f2cb",
"metadata": {},
"outputs": [],
"source": [
"import pandas as pd\n",
"import matplotlib.pyplot as plt\n",
"import numpy as np\n",
"from sklearn import linear_model"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "426f74bc",
"metadata": {},
"outputs": [],
"source": [
"df = pd.read_csv('french_coal_production.csv')"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "1dc7db0a",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"date int64\n",
"metric_tons float64\n",
"dtype: object"
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"df.head()\n",
"df.tail()\n",
"df.dtypes"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "38d53a2e",
"metadata": {},
"outputs": [],
"source": [
"X = df['date']\n",
"X = X.values.reshape(-1, 1)"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "cb3433a1",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"0 13.0\n",
"1 12.9\n",
"2 15.4\n",
"3 17.0\n",
"4 16.4\n",
"Name: metric_tons, dtype: float64"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"y = df['metric_tons']\n",
"y[0:5]"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "a2f97fc8",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"Text(0, 0.5, 'Coal Production in Metric Tons')"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAEWCAYAAACT7WsrAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAAAxdElEQVR4nO3debgcZZn+8e8NBsMWEiBAgIQIwiAwiBKRgRlFHBRREIdxjWyi6OACDm7IqCiioOP2G5cxChg1oLgzOICIgIAjmChhkV0TWQJEFoEAMZDn98db7an0qepT3af79HZ/rquv7q6uqn7fc7r7qXdXRGBmZlZkrW4nwMzMepeDhJmZlXKQMDOzUg4SZmZWykHCzMxKOUiYmVkpBwmrRNKlkt7cA+mYLSkkPa0L792Rv4GkD0r6ervPa9YODhJ9QtISSY9LejR327Lb6SojaYqkz0v6U5bW27Lnm3b4ffN/p3slnSlpg06+ZzMk7SPpzvy2iPhERHQi+MyQdK6ku7PAOrvu9Y0lfVfSn7PbAklTstf+qe6z9mh2jkNyx79b0j2S/iLpDElPb5CWXSRdmL3PqMFZkp4l6RfZuW6T9Krca+tI+n72vw1J+9QdO1XSfEn3ZbeTWv2b2WgOEv3lwIjYIHe7O/9iN66ui0haB7gY2BnYH5gC7AXcD+wxAUk4MCI2AJ4LPA/4j4I09sTfqsNWAxcAh5S8/nFgGrAtsB2wOXASQERcnv+sAa8AHs3Oh6SXAh8AXgzMzs7x0QZpWQWcAxxV/0L2v/gJcB6wMXA08G1JO+R2uwJ4I3BPwbk/B6yXpWMP4FBJRzZIizUjInzrgxuwBPjngu0BvB24Ffhjtu0VwDXAQ8CvgF3rzvMe4FrgL8B3gcm511+ZHfswcDuwf7b9UuBk4ErgEeBnwKYlaX0zcC+wQYP8PCs750PADcBBuddeDvwuS8MdwEm512ZneX5alb8T8GngvAZ/q7cAtwEPAOcCW+aO3Q+4Kfs7fRG4DHhz9tpJwLfL0kX6sTsTuBt4EPgxsD7wOOnH+9HstmXBuQ7K/iYPZX+jZ1X9/5X8TZ6WpW123fbzgWNyz98OXFhyjjOBM3PPzwI+kXv+YuCeCp/jZwJRt22X7G+h3LafAScXHH8nsE/dtj8Dz8s9/yBwebe/s4Nyc0liMBwMPB/YSdJzgTOAtwKbAF8Fzq2rCngN6Qr/GcCuwBEAkvYAvgm8F5gKvID0o1TzBuBIYDNgHdKPVZF/Bi6IiEeLXpQ0Cfgf0g/BZsA7gQWS/i7bZQVwWJaGlwP/Jungxn+CwveZCRxACjg1BzPyt9oX+CTp7zEDWAp8Jzt2U+AHpFLIpqSAuXcTb/8t0tXtzqQ8fi4iVgAvA+6O8tLgDsDZwHHAdOB/gf/JSmc1hf+/FnwJeIWkaZKmkUoc59fvJGk94F+B+bnNOwOLc88XA5tL2qSFdKhk2y4tnqPZY60BB4n+8mNJD2W3H+e2fzIiHoiIx0lXxl+NiKsi4qmImA+sBPbM7f//IuLuiHiA9GO9W7b9KOCMiLgoIlZHxF0RcVPuuDMj4pbsfc7JHVdvE2BZg3zsCWwAnBoRf42IX5CqGl4PEBGXRsR1WRquJf1ovnCMv03ejyU9RKqiuAz4RO61/N9qbpbf30bESuAE4B+yuvsDgN9HxPcjYhXweYqrOkaRNIMUDN4WEQ9GxKqIuKxi2l8L/DT7H6wC/hNYl1RdV1P2/2vWb0nB/v7s9hTw5YL9DiFdrefzsAGpJFNTe7xhC+m4CbgPeK+kSZJeQvp/r1fx+AuAD0jaUNIzgTc1cayNwUGivxwcEVOz28G57XfkHm8DHJ8LJg8BM0nVGjX5H7vHSF94sv1ub/D+ZcfVu590ZV5mS+COiFid27YU2ApA0vMlXSJpuaS/AG8jXc1XVfs7bRMRx2QBoSb/t9oye18AspLP/Vk6tszvGxFRd2wjM4EHIuLBJtJclqbV2ftuldun6v9hLN8DbiH9sE8h/e+/XbDf4cA3s79BzaPZMTW1x49Imptr7B5VMqmXBcODSaXGe4DjSRchdzY4LO9dpGq8W0ltG2c3cayNwUFiMOS/vHcAp+SCydSIWC8izq5wnjtIDZjj9XPgpZLWL3n9bmCmpPznbxZwV/b4LFL7wMyI2Aj4b4qrJFqR/1vdTQqqAGTp3SRLxzLSj33tNeWfk6rE8lerW+Qe3wFsLGnqGO9fpD5Ntfe9q/SI1j2bVOpckQXI/yaVoP4mq7Lbh1QNmXdDdnz+XPdGxP0RsSBXnfayKgmJiGsj4oURsUlEvJTUEH51xWMfiIi5EbFFROxM+l2rdKyNzUFi8HwNeFt2NS5J60t6uaQq1QCnA0dKerGktSRtJWnHFtLwLdIP5Q8k7ZidaxOl8QAHAFeRfmTfl1Uv7AMcSNYeQLqyfSAinsjaSd7QQhqqOIuU392yNptPAFdFxBLgp8DOkv4l633zLtYMBNcAL5A0S9JGpKoqACJiGalu/8tZff8kSS/IXr4X2CQ7psg5wMuz/8Ek0lX1SlIHhKZJmgzU2qOenj2v+Q3wZknrSlqX1Ktocd0pDgV+FRH1JcxvAkdJ2ilrz/gP4BsN0qHsvdeppSvfTiZp12zbepLeQyqJfiP3ej7t62T7Knttu+zztbakl2X5+PhYfxurxkFiwETEQlK7xBdJvWpuo2LDZkRcTWqY/hypjvkycle1TaRhJanx+ibgIlIvpatJVUZXRcRfST14Xkaq6/4ycFiu/eMY4GOSHgE+TPrhbLuIuBj4EKmBehmpFPW67LU/A68GTiVVQW1P6tlVO/YiUs+ia4FFpDaVvENJ3T5r9e3HZcfdRKoO+UNWHbjGWJeIuJnU1fO/SH+bA0ldev/aYjYfJ1UNkaUlX/X2JlKvrDtJJZVtGf1ZOYw1G6xr6bwA+BRwCal6bCnwkQbp2CZ77xty6bo59/qhpP/BfaSeUvtln6Oam7NjtgIuzB7XPpu7A9eRet19EpgbETdgbaE1qxnNzMxGuCRhZmalHCTMzKyUg4SZmZVykDAzs1J9McnZpptuGrNnz+52MszM+sqiRYv+HBHTx3OOvggSs2fPZuHChd1OhplZX5G0dOy9GnN1k5mZlXKQMDOzUg4SZmZWykHCzMxKOUiYmVkpBwkzG1oLSDMcrpXdL+hmYnpUX3SBNTNrtwWkOcUfy54vzZ5DWrLQEpckzGwonchIgKh5LNsOLmXUdDxIZAuB/E7SednzjSVdJOnW7H5ap9NgZlbvTyXbl5JWnzo6exyMlDKGMVBMREniWODG3PMPABdHxPbAxdlzM7MJNavBa2OVMoZJR4OEpK1Ji5t/Pbf5lYysdDWftAC6mdmEOoVsLdWc9YCvNDimrPQxyDpdkvg88D5gdW7b5tkawLW1gDcrOlDS0ZIWSlq4fPnyDifTzIbN64CNgUmASGuhzgPeRvmavY1KH4OqY0FC0iuA+yJiUSvHR8S8iJgTEXOmTx/XJIZmZqP8BLgH+DbpKnYJI72aTiGVKvLWy7YPm052gd0bOEjSAcBkYIqkbwP3SpoREcskzSAtfG5mNmECOA3YFjik4PVasHgz8ASpZHEKw9k1tmMliYg4ISK2jojZpJLdLyLijcC5wOHZboeTArqZ2YS5DLgaeA+wdsk+c0kNqDuwZilj2HRjnMSpwH6SbgX2y56bmU2Y00iNoUeMsd804MGOp6a3TciI64i4FLg0e3w/8OKJeF8zs3rXAhcAHwfWHWPfWpAIUuP2MPKIazMbKp8CNgCOqbDvNOBJYEVHU9TbHCTMbGgsAb5DGj1dZaqH2j7DXOXkIGFmA682D9MzgKeyx1U4SHgWWDMbcPWzvUKaC2hjxu6x5CDhkoSZDbjxzMPkIOEgYWYDrmy+pSrzMDlIOEiY2YDbqGR7lXmYHCQcJMxsgH0aeIjRo6qrzsO0IelH0kHCzGwA5FeT25g0BfXrgDNJ8y/lZ3utMs3GWsBUhjtIuHeTmQ2E+l5MD5JKEAcAh2a3Vgz71BwuSZjZQCjqxfQU8KFxntdBwsxsAIynF1MjDhJmZn0k3+4wG/giaTbXKNl/vKvJDXuQcJuEmfWN+naHpcA7SW0PLwd+ATye278dq8kNe5BwScLM+kZRuwPA5sB5wNdorRdTI/npwoeRSxJm1jfK2heWZfdzaf8KctOAVaTgtH6bz90PXJIws56Ub3vYGvgHOtfu0Miwj7puKkhImiZp104lxswMRtoelpICw13Ar4E5jF5Nrh3tDo04SIxB0qWSpkjaGFgMnCnps51PmpkNq7K2h+V0pt2hEQeJsW0UEQ8D/wKcGRG7A//c2WSZ2bB6jFSCKPInUkBYAqzO7jsZIMBBokqQeJqkGcBrSB0IzMzaJt/2sAWwbYN9O9n2UMZBYmwfAy4EbouI30jaFri1s8kys2FQ3/ZwL3AfcBCprSGv020PZRwkxhAR34uIXSPimOz5HyLikM4nzcwGXVHbQ5AaP+cxsW0PZaZkaRjWIDHmOAlJ04G3kEqEf9s/It7UuWSZ2TBoNN9SJ8Y8tGLYpwuvMpjuJ8DlwM9JkyqambXFLIobqbvR9tDIME/NUSVIrBcR7+94Ssxs6HyUNDlfXrfaHhoZ5iBRpeH6PEkHdDwlZtb36mdoXTDG/ptk99PpfttDIw4SjR1LChRPSHokuz081kGSJku6WtJiSTdI+mi2/SRJd0m6Jrs5AJkNgPqeSkuz540CxVmkQHEXEzfuoRXDHCTGrG6KiA1bPPdKYN+IeFTSJOAKSednr30uIv6zxfOaWQ8q6qn0WLa96If/UVKD52HApM4mbdwcJMYg6SDgBdnTSyNizEF1ERGkzwGkz8Akhne2XbOB1+zKcD8hBZFeLDnUy08Xri6nZaJVmbvpVFKV0++z27HZtjFJWlvSNaTxMRdFxFXZS++QdK2kMyRNKzn2aEkLJS1cvnx5lbczsy4q65FUtv2s7LW9OpOctpoG/JU1FzQaFlXaJA4A9ouIMyLiDGD/bNuYIuKpiNiNNNPvHpJ2Ab4CbAfsRpoG/jMlx86LiDkRMWf69OlV3s7MuujQgm2TKO6ptJw0jcPr6Y/1CoZ51HXV/8/U3OONmn2TiHgIuBTYPyLuzYLHatKEjns0ez4z6y1B+oJvBMwkVclMzm6vLNj/e6RBV/1Q1QQOEoUk/Sx7+Engd5K+IWk+sAj4xFgnljRd0tTs8bqkmWNvyiYLrHkVcH2LaTezDmqmO+uFwBWkH4Y/kXoqXQI8Avy/gv3PAnYB/r5tqe2sYQ4SjRqupwNExNmSLgWeR7pAeH9E3FPh3DOA+ZLWJn3OzomI8yR9S9JupIuPJcBbW0++mXVCrTtrrbdSrTsrjL76D+A/SIHkzbntewKvAD4NHMNIdcQS4EoqXGn2EAeJYhtJ+peC7XtJIiJ+2OjEEXEt8JyC7UVVl2bWQ5rpzvojUvXCmcA6da+dTPoR+Ez2GODs7P717UrsBHCQKLYR6UKgqMdXAA2DhJn1r6rdWZ8CPgz8HfDGgv13Iy1E83ngXaTqibOAvUklj34xNbt3kFjTUs/0ajacyibemwz8hZHeK98Bbsjuy35MPgp8HzgVOJzUCPmldiZ2AtTyO4xBolHvpmEbM2JmmVMYvejPJOAJUuPkaaR5lt6YbV/V4Fw7kkoOnwWenW3rh26veWuTAoWDxJrcdmA2pOYCX8w934bU5nAZafW4DzBS9bSK1PukrPfTAmBh3bbjG+zfq4Z1ao7SIBER7ppqNsR2zu5/yMjEe/8EFE3mVmvULnIio0cqN9q/VzlImJnlLM7un123/e6S/Zudu6lse69ykCghaX1Ja+WeryWpvrrSzAbMYlKpYXbd9mbnaGp2e69ykCh3MWu2Ya1HWsrUzAbYYmBXRv9IFDVqN1pNrtn9e5WDRLnJEVGb8pvssUsSZj2o2ZXhygRwLaOrmiC1TcwjNWZXWU2u2f171TTgoW4noguqrCexQtJzI+K3AJJ2ZzhnzDXrac1MpTGWJcDDFAeJ2vmaOWez+/eiaaQuwE+QxosMiyolieOA70m6XNLlwHeBd3Q0VWbWUFGJ4YOUT6XRrLJG62E2rFNzVFm+9DeSdiSNvBdwU0Q0GjtjZh1UVGI4kvIBba30IlpM+rLv0sKxgyofJGY02nHAlAYJSftGxC8KJvnbvsoEf2bWGUWT7zW6atu6hfdYDGwPrN/CsYPKJYnRXgj8Ajiw4DVP8GfWJY1KBusxOoBsRlrfoZlBUYuB3ZtM16Ab1iDRaMT1R7LxEedHxJF1N0/8ZzYB8m0P2wBvoXxStVqvoXwvojeSpvFuprvpw8AfcHtEPQeJAtkSo26kNuuCWtvDUlLR/U/A14FNGd27pjbuYC6pZ9Lq7P6bwGGk6byPp1r32OuyeweJNTlIlLtI0nskzZS0ce3W8ZSZDbmitgdIAeLrVBt3IOC/SUHhs4wEnFr32KJA4Z5NxaZm98MWJKqMk6hVLb09ty2AbdufHDOrKWt7uIPmxh2sS3HDdtlKc4tJV82tNHgPsrWBKThIFHlWRDyR3yBpmMaSmHXFVsCdBdtbmfOomUn5riGVIrygzGjDODVHleqmX1XcZmZttEPBtlbnPKo6yd5TpDYJVzUVG8Yg0WicxBaki5l1JT2HkQuLKXjuJrOO+i1wCbA/cCPpin8WI43TzTqFNQfgQXHAuY00546DRDEHiTW9FDiCVDX52dz2h0kzAJhZBwTwTlIvprMZaTAdj1pg+SAp4Dyd4sZuN1o3Ng24uduJmGClQSIi5gPzJR0SET+YwDSZDbUFpPrc02lPgKipNXafCpwAzCnYZzGpgXanNr7vIBnGkkSVNokrJZ0u6XwASTtJOqrD6TIbSo8A7wOeRyrGd8IRpKvDeQWvLQZ2ZLhmOW2Gg0SxM4ELgS2z57eQZoY1szb7OLAM+C86t7bwFsDBwDdI017nLcZVTY1MI7XZrOx2QiZQlc/hphFxDmkQJxHxJKkThJm1QX7qjU8B/wQ8v8Pv+VbgAdacgO0BUpdbB4lywzjqukqQWCFpE1J7GpL2BP7S0VSZDYn6qTcAFtL6inJV7UsaDfvV3DY3Wo/NQaLYvwPnAttJupI0Hcw7xzpI0mRJV0taLOkGSR/Ntm8s6SJJt2b308Y6l9mgKpp643FaWyioGWuRgtMvgZuybQ4SY3OQKJAtW/pCYC9SKXXniLi2wrlXAvtGxLOB3YD9s1LIB4CLI2J74OLsudlQKpt6o5WFgpp1JDCJkQbsxaRpxbeYgPfuV8MYJBoNpqtfbKhmhyqLDkVEAI9mTydltwBeCeyTbZ8PXAq8v3qSzQbHLFJVU9H2TtuM1IA9H/gEbrSuwkFiTd8nTeNyTfY8P5VLpUWHJK1Nms7+mcCXIuIqSZtHxDKAiFgmabOSY48mW8d91qyJ+MqYTbzjgXfVbWt16o1WvBX4HvAd4IaCtNiaHCTWdAjwWmBX4CfA2RFxWzMnj4ingN0kTQV+JKnykrkRMY+sJDxnzpwYY3ezvvR/pC/h5qRJ+MYz9UYrXkQqUbyZ1GXxTFLd8ES9f7+Zmt0PU5BotDLdjyLidaT2iNuBz0i6QtILm32TiHiIVK20P3CvpBkA2f19LaTbrO/9kjTtxgdJXU9rCwVN5A/02aQfvFqf9vspX2fCUkDfEAeJek+Qurw+TFoXvdJgTEnTsxIEktYF/pnUkeJc4PBst8NJpRSzofIkqYvgLLrbIHcio9eaqK0zYcWGbdR1o4brFwGvB/YAfg58ISIWNnHuGaS5n9YmBaNzIuI8Sf8HnJNN7fEn4NUtp96sT80DriW1B3RzSuVu9q7qV8MWJJQ6IRW8IK0mfY6vIDVUr7FjRExYG9ecOXNi4cJm4pNZ7/ozaa2I55Cuvrq5uM9sintXbUOq+rLRXkSqnvtltxNSgaRFEVE0l2NljRqujxzPic1sTQtI1Ti1H+WX0P3V36quM2EjpgG3djsRE2isqcLNrA1q02/kf4w/RlqspZs9iWrvfSLjX9hoWAxbdVOnJpo0s5yi6Td6pYF4LqlqqRu9q/qRg4SZtZ0biAfHNFKA/2u3EzJBHCTMJsD0ku2eS6D/DNuo60YN10Aa7wC8hdQR4m/7R8SbOpcss8FxN+nKU6zZRdANxP0pHyQ272ZCJsiYQYI02O1yUm89LzZk1oSngDeS6vtPA76EG4j7nUsSo60XEZ6l1awFnwQuAc4g9Sl/b3eTY23QziBR6xbdyxcOVdokzpN0QMdTYjYg8suRfgj4B+CILqbH2qtdQaJ+VcKl9Oa8WVWCxLGkQPG4pIclPSLp4U4nzKwfFS1Heg1wVrcSZG3XriDRy92i86qsTLdhRKwVEetGxJTs+ZSJSJxZv+nWcqQ2caZm9+MNEv3SLbrRBH87RsRNkp5b9Hq2rKmZ5fTLF99aN4k0HfZ4g0Q3VyVsRqOG638nlZw/U/BaAPt2JEVmfWgVcCp1s2Dm9NoX38anHaOu+2XerEZzNx2d3b9o4pJj1h/yvVK2IC2y8kdSI/U1pCqmml784tv4tCNIzCWtK3JEbttp9GfvJjPLqW+cXkYKEO8EfgV8jTTVtrL7efTeF99atwC4hTSAbDbj6430j9n98dn9pHGcq1McJMyaVNQ4DWnJRfCEeYOsdoGwMns+3m6rtTaJA4BnAP8zrtR1hoOEWZPcOD282t1ttRYktgEOIk1rsaLFc3VKpSAhaStJe0l6Qe3W6YSZ9aqyRmg3Tg++dl8gLCVVS84kBYmVwEUtnqtTqkzwdxrwWuD3jMzdFPTH6n1mbXcSo5dtdOP0cGh3t9WlwAxgHeCfgI1I1ZYHt3i+Tqgyd9PBwN9FxMqxdjQbBrXqgM2A5fTunDvWfu3utvonUlUTpEbrA4DzSFfja7d4znarUt30B3qz0d1swq0kjYfYG7gHN04Pm7mk3mpbZM+nM77ea0sZCRKQqpyWA1e1msAOqFKSeAy4RtLFjDTqExHv6liqzHrUfOBO4OukumQbPnOB/UhrSXyY1gPEauAO4F9z2/Yn/SifC+w1jjS2U5WSxLnAyaQu4ItyN7OBkZ+5dTbFXRpXkab+3gN4yUQlzHrSdFI7wh3jOMc9pCVQ8yWJqcAL6K2usFUm+JsPnM1IcDgr22bWk8p+8BttrzJl8wJS1dKHcCli2AnYmvEFiVoDeH2j90GkXkK3jePc7TRmkJC0D3AraVGtLwO3uAus9aqyH/xjSrZ/hjQX/lh9358kNU4+B3h555JvfWQmqeqxVfkxEnkHZve9Upqo0ibxGeAlEXEzgKQdSCWL3TuZMLNWlA12+krBvo8B72lwrnzf9++Srux+iEsRlmwNXDmO48uCxLbALqR6/neP4/ztUqVNYlItQABExC24t5P1qFYGNc0o2T6JNCZiG9I61ZPovdGw1j0zgbtIDdCtWEqaKHDDgtcOAi4HHmjx3O1UJUgslHS6pH2y29eo0HAtaaakSyTdKOkGScdm20+SdJeka7Kbl0a1tplZsr2sz/k2wKdJfd3z1iF9OT7KSOBZBbyV3lte0rpja9Jn4r4Wj6/v/pp3EGmsxPktnrudqgSJfwNuAN5Fqr79PfC2Csc9CRwfEc8C9gTeLmmn7LXPRcRu2e1/W0i3WaGi9oL1SO0P9YGgNgiq1vc9P3PrGaQeLPV6cXlJ647aBUmr7RKNgsTzgCmkz22jHncTYcw2iWyk9WezW2URsYw0izIR8YikG4GtWkmkWRWrSPPezCT92N/BmqOh92ZkDYj6UdJzGd3f/dCS9/FEfgYjQeIOYE6Tx9Y6T7y45PWzSVWbtXmQah0tYOIHbpaWJCSdk91fJ+na+lszbyJpNqljSG0g4Tuy85whaVrJMUdLWihp4fLly5t5OxtS80mNy18ifanqR0M3O4W3J/KzRrbO7lvpBvsg8CjlJYkTGQkQNd0qxTaqbjo2u38FqVdW/a0SSRsAPwCOi4iHSR1NtgN2I5U0ipZHJSLmRcSciJgzfXpRwd9sxErgY8DzSR/YdjiF8ioqs02Bp9NadVNZz6aaXpqOvjRIZNVFAMdExNL8jdTtfEySJpECxIKI+GF23nsj4qmIWE1axGuP8WXBDL5KuqL7OO3rolrUVuFV5qxmPAPqxgoSvVSKrdJwvV/BtpeNdZAkAacDN0bEZ3Pb8z0OXwVcXyENZqVWkK7u96G8jrdVXmXOGml1QN1YQaKXSrGlDdeS/o1UYtiurg1iQ9I8TmPZm9T2d52ka7JtHwReL2k3UtvNElKvQrOW/RepG+KP8EA3m1gzaW1hnaXAuqQqqyK1i5GyjhYTqVHvprNI3XQ/CXwgt/2RiBhzjEdEXEHxd9ZdXq0tFpA+mHcCk4E/0jszZ9pw2Jo0oK7Z9R9q3V8bXdQU9bjrhkZtEn+JiCXAF4AHcu0RqyQ9f6ISaFakNkdTraj/BONbkN6sFTNJA8KaHVCXX2yo11Vpk/gKqbdWzQqKp8IxmzDtXpDerBWtdoNtNJCu11QJEoqIqD3JeiVVmRjQrGN6qYugDa9WRl0/xsiyt/2g0vKlkt4laVJ2O5a0pKlZ12xZsr1fvng2GPKjrquqXcgMUknibaT2wLtIAfP5jIwQN+uKZxVs80A3m2gbkzpNNBMkxur+2muqzN10H/C6CUiLWSVLgMuAfYHb6X4XQRteovmxEgMXJCSdSRrTsIaIeFNHUmQ2ho+RisDzGWk4NOuWZkddLyV1ly2rMu01VRqgz8s9nkwaJX13Z5Jj1tjNpODwLhwgrDfMBC5tYv+lpM9uv/T+qVLd9IP8c0lnAz/vWIrMGvgIaaTqCd1OiFmmtkJd1QF1/dT9Fao1XNfbHncisS5YTFpr+lhgsy6nxaxma1KAuLfi/v0WJKq0STxCapNQdn8P8P4Op8tslA8BGwHv6XZCzHLy3WDHamdYRSp1DFSQiIiidbrNJsQCRiY5C+BfSYvHm/WK/KjrseYruos0o/BABAlJz210YET8tv3JMRtRm58pP/3GT7Pt7upqvaKZUdf91v0VGpckaivGTSYt4bqYVOW0K2kZ0n/sbNJs2BXNz/R4tt1BwnrFNFJniirdYPsxSDSaBfZFEfEiUr6emy0lujtprerbJiqBNrw8P5P1g2YG1NWCRD/1/KnSu2nHiLiu9iQirietT23WUb20hKNZIzOpXpLYnFQ90y+qBIkbJX1d0j6SXijpa8CNnU6Y2TsKtnl+JutFVUdd99M6EjVVgsSRwA2k7unHAb/Ptpl1TJAaqdclfQFF+nLNw+0R1ntmAstICxA1spT+KwlX6QL7hKQvkUZZB3BzRKzqeMpsqH2DNNXBV/GUw9b7agPq7qF8upgglSQOnKhEtUmVwXT7kKbLWULWRiPp8IhoZf1vszHdRxow94/Am7ucFrMq8t1gy4LEfaRldvutuqnKHFOfAV4SETcDSNoBOBvYvZMJs+H178AjpKqlVuaNMZto+VHXe5bs04/dX6Had3BSLUAARMQtwKTOJcl6zQJgNunDMjt73sn3WAC8guKFhcx6UZW1rvs1SFQpSSySdDrwrez5XGBR55JkvaR+1PNSRtoI2tWAXDSy+kI8str6x1RgfRqPlejXIFF1+dIbSFP4H0vq3fS2TibKekfRqOfHsu399B5mnSQad4NdQFosC9Igs06UxjulYUlC0lrAoojYBfjsxCTJeslEjHr2yGobBGWjrieiNN5JDUsSEbEaWCyp37r2WhvcRPkHpJ0fiLLeIP7QWT8pG3Xd7yXlKm0SM4AbJF0NrKhtjIiDOpYq67qrgQNII5xXkbru1UymvaOed2b0l8sjq63fbM3IgLr8D2u/l5SrBImPtnJiSTOBbwJbkKZQnxcRX5C0MWmBsdmksReviYgHW3kPa5/8ug3TgQdJV0YXkqb8rb0mYCvgtW16319m77EfcEv2HrNIAaIfiuJmNTNJP3TLGOkSC+nzvLRg/34pKZdWN0maLOk44NXAjsCVEXFZ7Vbh3E8Cx0fEs0hdh98uaSfgA8DFEbE9cHH23LqoVme6lDQq9D6yfx7wTNKP9RLSF2ABcDvwn21430eBI4BtgR/l3mMJDhDWf8q6wR5WsG8/lZQbtUnMJ60jcR3wMkbWl6gkIpbVFiaKiEdIkwJuBbwyO3ftPQ5uLsnWbkV1pgF8qmDf15KuGj5M+mCMx3tJAeEbpO6DZv2saPGhAH5G6iI7k/6cg6xRddNOEfH3ANk4iatbfRNJs0nrUFwFbB4RyyAFEkmFa9pLOpqsE8CsWf1SMOtPzdSZCvgycBlpDpogXTlVrSKqX470ALx6lQ2GopLE90g/eqcDb5rwFLVHo5LE3ybxi4ixJjcsJWkD4AfAcRHxcNXjImJettDRnOnTp7f69jaGVaSib5Gy0LwpKRgsZeTHvtatr1H/7/pqLUiT+PVTn3GzMhsBGzBSklgJnAD8PXB4txLVBo2CxLMlPZzdHgF2rT2WVOnHXtIkUoBYEBE/zDbfK2lG9voMUhW4dcEKUt3fCkbPszJWnekPC7aN1a2v37sCmjVSW6GuVpL4CvAHUrXt2t1KVBs0Wr507YiYkt02jIin5R5PGevEkkQqZd0YEfmBeOcyElgPB34yngxYc/JzJG0MnE+qHz2TVFdatc60URVV0VxPT1Dcw6PRucz6TW3U9UPAyaReey/tZoLaoEoX2FbtDRwKXCfpmmzbB4FTgXMkHUX6fXh1B9NgOfUjP/8KPJ1UaphLcw1pZd361gaOIhW1yfY5EjhmjHOZDYKZpDmMPkHqRv5p0oVXP+tYkIiIKyj/+7y4U+9r5d7P6OqelaTqnmZ7WpzC6En5JpG6ztY3YK0iBY8TgC/UHdNPXQHNxvIgcDcpOKwPXA88u6spGj9P1z+g8lU+M4C9gLtK9m2lumcuqUoqX0V1ZoP9V5KuruqP6aeugGaNLCAtuVuzgrE7c/QDRcTYe3XZnDlzYuHChd1ORt8omnob0nQaT4zenW1I4xXaYTbF1VDtfA+zXjSb3vvsS1oUEXPGc45GI64fyfVuyt8q926y7ijqRQSp+Fvf3bXd1T2nTMB7mPWifp+jqUyj3k0b5no3Tanr6TRm7ybrnrIP5QN0vrqnqBrKVUo2DMo6YPR7x4zKDdfZyOjJtecR0e8BcmBtRfG89rNovhdTKybiPcx6TVFnjkEoRY/ZcC3pIEm3An8kzcawhNS93nrUPgXbBuHDatbLBrUUXaV308mkWVxviYhnkLqvXtnRVFnLgrQA+XYM3ofVrNflZ0xewmB856pUN62KiPslrSVprYi4RNJpHU+ZteRK0nS7/TyhmJn1jipB4qFskr7LgQWSassNWA/6KjCF9i0KZGbDrUp10yuBx4HjgAtIa84c2ME0WYvuJ01N/Ea8PoOZtceYJYmIWCFpc+B5pN+h8yPi/o6nzJr2LdLI5rd2OyFmNjCq9G56DWnBoVcDrwGukvSvnU6YNSdIVU17Art2OS1mNjiqVDedCDwvIg6PiMOAPYAPdTZZvaVo6utecwVwE9lSfmZmbVKl4XqtiMgvDHQ/QzQxYP08SLUV2KC3urd9lbQylhuszaydqvzYXyDpQklHSDqCNNHh0Aym64fV1O4Hvk9avKNsKVIzs1aMGSQi4r2kC9VdSVOjz4uI93U6Yb2i3ZN2tbPqqnauTUkN1luN41xmZkUazQL7TEl7A0TEDyPi3yPi3cD9krabsBR2WdnkXJuNcVxRMKhVXS0lNTTXqq5aCRT5c9Wc3OK5zMzKNCpJfB54pGD7Y9lrQ+HjjF5eT8C9wAeA+VQPBu+kcdVVWSmjaHs/VIOZWf9r1HA9OyKurd8YEQslze5cknrLWqQf+k1IU23PInXtugo4jRQwass2LSVNhfE0in/Ai9Z4qB33deBYRjeQX0kKRPnthwNPlZzLU/OaWTs1ChKTG7y2brsT0oueIK3L/FzgN6xZ7DoK+DGwvO6Yv2a3Zr2lYNtjwFcKtpcFCOj/uevNrLc0qm76jaRRv12SjiJNNDrw/ot0Zf5piv9Qf27yfJtQvGrbh5tP2t+OrX/u6cDNrJ0alSSOA34kaS4jQWEOsA7wqg6nq+vuJ/3gHgDsW7LPLIrXtN2ENNlV/eIjX8gen0gKPrOy95hLqlIqOtfaFJcctsmOLTqXmVm7lAaJiLgX2EvSi4Bdss0/jYhfTEjKuuzjpFb7TzXYp2wlqkbBAIp/yMvOdThrtknUttfO56BgZp1UZYK/S4BLJiAtPeN24EukRuidG+xX+4FuJhi0cq69G7yHmVknVV7jehjUupYuJfVa2q3CMe28mi87l0sMZtYtDhKZ+jmaAngfMBX/QJvZ8BqaifrG4sFpZmajdSxISDpD0n2Srs9tO0nSXZKuyW4HdOr9m9XuOZrMzAZBJ0sS3wD2L9j+uYjYLbv9bwffvyllg9A8OM3MhlnHgkRE/JI0k0VfOIXRczR5cJqZDbtutEm8Q9K1WXXUtLKdJB0taaGkhcuX109+0X47khqrNyYFi22AebjR2syG20QHia8A25F6ly4DPlO2Y0TMi4g5ETFn+vTpHU/YAtJQ8tuA1cASHCDMzCY0SETEvRHxVESsBr5GWi+7654CvkOagqO0aGNmNoQmNEhImpF7+irg+rJ9J9JlpGKNSw5mZmvq2GA6SWcD+wCbSroT+Aiwj6TdSNX/S4C3dur9m7EA2BB4ebcTYmbWYzoWJCLi9QWbT+/U+7XqCeAHwL8wJItkmJk1YehHXP8v8Bdc1WRmVmTog8RZwObAi7qdEDOzHjTUQeIvwHnAa/FMh2ZmRYY6SPwQWImrmszMygx1kDiLNLLved1OiJlZjxrYILEAmE3K4Ozsed4y4BfAGxg9Z5OZmSUDGSRqCwgtJQ3IWJo9X5B7fWfS9BunMzqAmJlZMpDttWULCL0XWAG8O/f63aQAAm6bMDOrN5AlibKFgpaRhnh7BTozs2oGMkiULRS0aYNjvAKdmdloAxkkTiEtGJS3HvB50joRRbwCnZnZaAMZJOaSFgzahtELCJUFEK9AZ2Y22kA2XEMKCEUN0bVtJ5KqmGaRAoQbrc3MRhvYINFIWQAxM7M1DWR1k5mZtYeDhJmZlXKQMDOzUg4SZmZWykHCzMxKKSK6nYYxSVpOmqevFZsCf25jcvrNMOffeR9ew5z/fN63iYjp4zlZXwSJ8ZC0MCLmdDsd3TLM+XfehzPvMNz5b3feXd1kZmalHCTMzKzUMASJed1OQJcNc/6d9+E1zPlva94Hvk3CzMxaNwwlCTMza5GDhJmZlerLICHpDEn3Sbo+t203Sb+WdI2khZL2yLbPzbbVbqsl7Za9truk6yTdJun/SVKXslRZk3mfJGl+lscbJZ2QO2bQ876OpDOzPC6WtE/umL7LO5Tm/9mS/i/Lz/9ImpJ77YQsjzdLemlue9/lv5m8S9pE0iWSHpX0xbrz9F3eoen87ydpUbZ9kaR9c8c0n/+I6Lsb8ALgucD1uW0/A16WPT4AuLTguL8H/pB7fjXwD6S1ic6vHd/Lt2byDrwB+E72eD1gCTB7SPL+duDM7PFmwCJgrX7Ne4P8/wZ4Yfb4TcDJ2eOdgMXA04FnALcDa/dr/pvM+/rAPwJvA75Yd56+y3sL+X8OsGX2eBfgrvHkvy9LEhHxS+CB+s1A7SpqI+DugkNfD5wNIGkGMCUi/i/SX++bwMEdSXAbNZn3ANaX9DRgXeCvwMNDkvedgIuz4+4DHgLm9GveoTT/fwf8Mnt8EXBI9viVpAuElRHxR+A2YI9+zX8zeY+IFRFxBfBEfud+zTs0nf/fRUTte3ADMFnS01vN/yAtOnQccKGk/yRVo+1VsM9rSV8egK2AO3Ov3Zlt60fHUZz375Pyu4xUknh3RDwgaQ6Dn/fFwCslfQeYCeye3a9mcPIOcD1wEPAT4NWkPELK069z+9XyuYrByX9Z3ssM0ncequX/EOB3EbFSUkv578uSRIl/I/0IzgTeDZyef1HS84HHIqJWp1dUF9ev/YHL8r4H8BSwJanK4XhJ2zIceT+D9CVYCHwe+BXwJIOVd0jVDG+XtAjYkFRahPJ8DlL+y/JeZpDyDmPkX9LOwGnAW2ubCs4xZv4HqSRxOHBs9vh7wNfrXn8dWVVT5k5g69zzrSmuouoHZXl/A3BBRKwC7pN0JTAHuJwBz3tEPEkKGgBI+hVwK/Agg5N3IuIm4CUAknYAXp69dCdrXlnW8jkwn/sGeS8zMHmHxvmXtDXwI+CwiLg929xS/gepJHE38MLs8b6kHwQAJK1FKo59p7YtIpYBj0jaM2vhP4xUbOtHZXn/E7CvkvWBPYGbhiHvktbL8oyk/YAnI+L3A5Z3JG2W3a8F/Afw39lL5wKvy+qinwFsD1w9SPlvkPdCg5R3KM+/pKnAT4ETIuLK2v4t57/brfYttvSfTapnr9WvHkXqzbCIVBd9FbB7bv99gF8XnGcOqV7vduCLZCPQe/nWTN6BDUhX1zcAvwfeO0R5nw3cDNwI/Jw0ZXLf5r1B/o8Fbslup+bzApyY5fFmcr1Y+jH/LeR9Camh99Fs/536Ne/N5p8UMFYA1+Rum7Waf0/LYWZmpQapusnMzNrMQcLMzEo5SJiZWSkHCTMzK+UgYWZmpRwkzApkY0uukPSy3LbXSLqgm+kym2juAmtWQtIupHEmzwHWJvU33z9GRrA2c661I+Kp9qbQrPMcJMwakPQp0sCk9bP7bUhTzj8NOCkifiJpNvCtbB+Ad0TEr5TWsPgIaRDUbhGx08Sm3mz8HCTMGsim9vgtafK084AbIuLb2dQHV5NKGQGsjognJG0PnB0Rc7Ig8VNgl0jTdZv1nUGa4M+s7SJihaTvkqZ3eA1woKT3ZC9PBmaR5o/6otKKh08BO+ROcbUDhPUzBwmzsa3ObgIOiYib8y9KOgm4F3g2qTNIfrGbFROURrOOcO8ms+ouBN5ZWxdY0nOy7RsByyJiNXAoqZHbbCA4SJhVdzIwCbg2W5D+5Gz7l4HDJf2aVNXk0oMNDDdcm5lZKZckzMyslIOEmZmVcpAwM7NSDhJmZlbKQcLMzEo5SJiZWSkHCTMzK/X/ASHqNObVuQfBAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.plot(X, y, color='cyan', marker='o')\n",
"plt.title(\"French Coal Production 1870-1919\")\n",
"plt.xlabel(\"Year\")\n",
"plt.ylabel(\"Coal Production in Metric Tons\")"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "20b9a3a7",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"(array([0.43477551]), -797.8321632653062)"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"regr = linear_model.LinearRegression()\n",
"regr.fit(X, y)\n",
"regr.coef_, regr.intercept_"
]
},
{
"cell_type": "code",
"execution_count": 8,
"id": "3bc79e5e",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"Text(0, 0.5, 'Coal Production in Metric Tons')"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAEWCAYAAACT7WsrAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAAAupUlEQVR4nO3dd5gUVdbH8e8BUcRAUEyoYFwJKiiLaVdQV9ccV901R1RMrDkuIubXAK4YQMwuGFlUVBQEEUmCBAkiQVxFVAQEBCWe949boy1O9/SE6pru/n2ep5/prurqOrenZ07fUPeauyMiIlKaGkkHICIi1ZeShIiIpKUkISIiaSlJiIhIWkoSIiKSlpKEiIikpSRRgMxsiJmdVw3iaGJmbmbrJHDuWN4DM7vBzB6v6tct45xvmdmZuTxnITKzU83snaTjyDdKEjEws9lm9pOZ/Zhy2yrpuNIxs43NrKuZ/S+KdUb0eNOYz5v6Pn1rZk+a2YZxnrM8zKydmX2Vus3d73D3OJLPWWY2rLR97n6Yuz9d1eesiCjpL41+Z3PM7H4zq5l0XNlw9+fd/ZCk48g3ShLxOcrdN0y5fZ26M4lv16Uxs3WBQUBz4FBgY2BfYD7QJgchHOXuGwJ7AH8EbiolxmrxXhWTMt7z3aPfWVvgZOCcHJ9fckhJIoeib2EXm9l0YHq07UgzG29mP5jZcDPbLeX5s83sKjObaGaLzOwFM6udsv+Y6NjFZjbTzA5NOV1jM/vQzJaY2TsZagVnANsCx7n7FHdf4+7fuXsXd38zOk/TqPnmBzObbGZHp8RwhJmNi2L40sxuqch74+5zgLeAFhneq/OjWs4CM3sttXZmZgeb2afR+/QQYCn7bjGz51Ie/6YZzMwaRLWYr81soZn918w2iOLZKrU2WMprHR29Jz9E71HTlH0Zf3/ZSm06K6lxmNm9Uayfm9lhKc+ta2a9zGxu9E3/tpJv+ma2g5m9Z2bzzex7M3vezOqtFe+1ZjYRWFrWP2p3nwF8CLRMeY1Mn+c9os/KEjN7KXo/bov2tTOzr6LzfwM8aWY1zOy66LM938xeNLMG0fNrm9lz0fYfzOwjM9s85T2aFZ3nczM7NfW9S4ln3+i4RdHPfdd6z7tk+TdU0JQkcu9YYC+gmZntATwBXABsAjwGvGZm66U8/yTCN/ztgN2AswDMrA3wDHA1UA/YH5idctwpwNnAZsC6wFVp4vkL8La7/1jaTjOrBbwOvBO91qXA82b2h+gpSwmJph5wBHCRmR2b+S0o9TzbAIcD41I2H8uv79WBwJ2E92NL4AugT3TspsArhFrIpsBMYL9ynP5ZoA6hNrUZ8IC7LwUOA77OUBvcGegNdAQaAm8Cr1uonZUo9fdXSXsB0whlvQfoZWYlSfFpYBWwI9AKOAQoaR4zwnu4FdAU2Aa4Za3X/gfh91jP3VdlCsLMdgH+DMyIHqf9PEfvSV/gKaAB4X07bq2X3CLa1xhoD1xG+Ay0jWJeCHSPnnsmUDcqwybAhcBPUXJ/EDjM3Tci1IrHlxJ7A6B/9NxNgPuB/ma2ScrTsv0bKmzurlsV3wj/rH8Efohu/422O3BgyvMeAbqsdew0oG3K65yWsu8e4NHo/mOEf2alnX8IcFPK4w6ERFDac98F7spQlj8D3wA1Urb1Bm5J8/yuJXEBTaIyr5PF+/QF8DCwfpr3qhdwT8rjDYGV0TnOAEam7DPgK+C86PEtwHMp+3+Ji5Bw1gD1S4mvHfDVWtt+eS3gZuDFlH01gDlAu7J+f6Wc6yxgWIbf53kpz5uRsq9OVJYtgM2B5SXvYbT/H8DgNK97LDBurd/HOWV8th1YTPhy4NFnYb2yPs+ELzFzAEvZNwy4LeW9XgHUTtk/FTgo5fGW0e98HUIT13Bgt7XOt0H0eToh9X1Y+z0GTgdGr7V/BHBWef+GCv2mmkR8jnX3etHt2JTtX6bcbwxcGVWXfzCzHwjfjFI7ub9Jub+M8M+R6HkzM5w/3XFrm0/440tnK+BLd1+Tsu0LoBGAme1lZoPNbJ6ZLSJ8oytPtbzkfWrs7h3c/aeUfanv1VbReQHwUPOZH8WxVepzPfxVpx6byTbAAndfWI6Y08W0Jjpvo5TnZPt7KI9fXtPdl0V3NyR8nmoBc1M+T48RvgljZpuZWZ+oGWox8By//11l877tEZ3vZEKtZoNoe6bP81bAnOh3k+5c89z955THjYG+Ka81FVhNSIbPAgOAPhaaCe8xs1oeaoAnEz6Hc82sf1TjWdtvfneRXz7XkTh+d3lHSSL31v4juT0lmdRz9zru3juL1/kS2KEK4hkI/DWqppfma2AbM0v9rGxL+FYI8B/gNWAbd68LPEpKf0Alpb5XXxP+aQAQxbtJFMdcwj+jkn2W+pjwrbdOyuMtUu5/CTRIbZtPc/7SrB1TyXnnpD0iXl8SahKbpnyeNnb35tH+Owll2s3dNwZO4/e/q6ymhfbgRcK373+lnD/d53ku0CilWQx++zsq7dxfEpqNUl+vtrvPcfeV7t7Z3ZsRmpSOJNQocfcB7n4w4cvPp0DPUorwm99dJPVzLREliWT1BC6Mvo2bmW1goSN4oyyO7QWcbWYHRR18jdJ8YyrLs4Q/xlfMbJfotTaxcD3A4cAowj/Za8yslpm1A44i6g8ANiJ8E/856ic5pQIxZOM/hPK2jPps7gBGuftsQttyczM7PupsvYzfJoLxwP5mtq2Z1QWuL9nh7nMJHdQPm1n9qIz7R7u/BTaJjinNi8AR0e+gFnAl4Z/08AqW0aIO2V9u5Tk4Kss7wH0WhjXXsNBZ3TZ6ykZEzXtm1ojQn1VZdwHtzWwLMn+eRxBqAZeY2Tpmdgxlj557FLjdzBoDmFnD6DjM7AAz29VCp/xiQjPUajPb3MJggg0Iv4sfo/Ou7U1gZzM7JYrnZKAZ8EYl34+CoySRIHcfA5wPPETolJtBlh2b7j6a0Kn2ALAIeJ/ffzPK5nWWEzqvPyX0TywGRhOaIUa5+wrgaEIn7veEfoMz3P3T6CU6ALea2RLCN8oXyxtDlnEOIvQBvEL4VroD8Pdo3/fAiYR/WPOBnQijbkqOfRd4AZgIjOX3/whOJ/yT+RT4jtARTVTG3sCsqMnjN9e6uPs0wrfxfxPem6MIQ3pXVLCY+wI/pd6s/ENBzyB0sk4hfKZe5tfmxM6EpqJFhMT6agXj/IW7f0L47F2d6fMcvSfHA+cS+gxOI/welmd4+W6EWuo70edrJKF5C8KXgJcJn9epUQzPEf6nXUmoKSwg9Id0KCXu+YTax5WEz8w1wJHRZ0lS2G+bCEVEcsPMRhE68p9MOhZJTzUJEckJM2trZltEzTtnEoYEv510XJKZrmoUkVz5A6E5ckPCyLy/Rf0oUo2puUlERNJSc5OIiKSVF81Nm266qTdp0iTpMERE8srYsWO/d/eGlXmNvEgSTZo0YcyYMUmHISKSV8xs7avKy03NTSIikpaShIiIpKUkISIiaSlJiIhIWkoSIiKSlpKEiIikpSQhIiJpKUmIiFRH8+dDx46waFGiYShJiIhUJ+7w0kvQrBl07w5DhyYajpKEiEh18fXXcPzxcNJJsM02MHYsHHVUoiEpSYiIJM0devUKtYe334Z77oGRI2G33ZKOLD/mbhIRKViffw7nnw+DBsH++8Pjj8NOOyUd1S9UkxARScLq1dC1K7RoAaNHwyOPwODB1SpBgGoSIiK5N2UKnHtuaFI6/HB49NHQB1ENqSYhIpIrK1ZAly7QqhVMnw7PPw9vvFFtEwTEmCTMrLaZjTazCWY22cw6R9tvMbM5ZjY+uh0eVwwiItXGmDHQujX8619hBNPUqXDKKWCWdGQZxdnctBw40N1/NLNawDAzeyva94C73xvjuUVEqodly+CWW+C++2CLLaBfPzj66KSjylpsScLdHfgxelgrunlc5xMRqXbefx/OOw9mzID27cPQ1rp1k46qXGLtkzCzmmY2HvgOeNfdR0W7LjGziWb2hJnVT3NsezMbY2Zj5s2bF2eYIiJVa9EiuPBCaNcuXAPx3nvw2GN5lyAg5iTh7qvdvSWwNdDGzFoAjwA7AC2BucB9aY7t4e6t3b11w4aVWsdbRCR3+veH5s2hZ0+48kqYOBEOOCDpqCosJ6Ob3P0HYAhwqLt/GyWPNUBPoE0uYhARidW8eXDqqXDkkVC/PowYAffeC3XqJB1ZpcQ5uqmhmdWL7q8P/AX41My2THnaccCkuGIQEYmdO/TpE6bUeOkl6NQpzLnUpjC+/8Y5umlL4Gkzq0lIRi+6+xtm9qyZtSR0Ys8GLogxBhGR+MyZAxddBK+/HpJCr17hCuoCEufopolAq1K2nx7XOUVEcsI9zLF01VWwciXcfz9cdhnUrJl0ZFVO03KIiJTHzJlhQr7Bg0OHdM+esMMOSUcVG03LISKSjdWrQ41h111Dn0OPHmHm1gJOEKCahIhI2SZNChPyjR4dFgF65BFo1CjpqHJCNQkRkXRWrIDOnWGPPcK6D336hGk1iiRBgGoSIiKlGz061B4mTQrXP3TtCptumnRUOaeahIhIqmXLwpXS++wDCxeG4a3PPVeUCQJUkxAR+dXgwWFCvlmzwtxLd98NG2+cdFSJUk1CRGTRojBL64EHQo0aMGRI6Jwu8gQBShIiUuxeey1MqdGrF1xzTZiQr23bpKOqNpQkRKQ4ffcd/P3vcMwxob9h1KjQvLT++klHVq2UK0mYWX0z2y2uYEREYuce1pZu1gz69g1rTpcsLSq/U2aSMLMhZraxmTUAJgBPmtn98YcmIlLFvvwyTOV92mmw004wbhzcdBPUqpV0ZNVWNjWJuu6+GDgeeNLd9yRM+y0ikh/WrAkd0c2bh07prl1h2LBQm5CMskkS60RrQJwEvBFzPCIiVWv69DARX4cOYTrvSZPg8ssLcsbWOGSTJG4FBgAz3P0jM9semB5vWCIilbRqFdxzD+y2G0yYEEYvvfsubLdd0pHllTIvpnP3l4CXUh7PAk6IMygRkUqZMCFMqTF2LBx7LHTvDlttlXRUeanMJGFmDYHzgSapz3f3c+ILS0SkApYvh9tug7vuggYNwnKiJ5wAZklHlreymZajH/ABMBBYHW84IiIVNGJEqD1MnQqnnw4PPACbbJJ0VHkvmyRRx92vjT0SEZGKWLoUbrwRHnwQttkG3noLDj006agKRjYd12+Y2eGxRyIiUl4DB0KLFtCtWxi9NGmSEkQVyyZJXE5IFD+b2ZLotjjuwERE0lq4MDQtHXwwrLsuDB0KDz0EG22UdGQFJ5vRTXrXRaT66Ns31BrmzYPrroNOnaB27aSjKlhZrSdhZkcD+0cPh7i7LqoTkdz69lu49NIwYqllS+jfPywrKrHKZu6muwhNTlOi2+XRNhGR+LnDM89A06ZhWu877ghLiypB5EQ2NYnDgZbuvgbAzJ4GxgHXxRmYiAhffAEXXAADBsC++4arpnfZJemoikq2U4XXS7lfN4Y4RER+tWZN6Ihu3jxMxPfvf8MHHyhBJCBtTcLM3nH3Q4A7gXFmNhgwQt/E9TmKT0SKzbRpYeTShx/CIYfAY49BkyZJR1W0MjU3NQRw995mNgT4IyFJXOvu3+QgNhEpJitXwr33QufOUKcOPPUUnHGGptRIWKYkUdfMji9l+75mhru/GldQIlJkxo+Hc84JiwD97W+heWmLLZKOSigjSQBHEmoPa3NASUJEKufnn8PyoXffHdaZfuUVOL6076aSlExJ4gvN9Coisfnww9D3MG0anH12aGpq0CDpqGQtmUY3qSFQRKrekiXhorg//zlM7T1gADzxhBJENZUpSZyesyhEpDgMGBAm5OvePSSKTz4JI5ik2kqbJNx9Ui4DEZECtmABnHVWmKG1Tp1w7UO3brDhhklHJmXI9mI6EZGKefnlMKXG88+HdR/GjQtXT0teyGb50g2An1Km5agB1Hb3ZXEHJyJ5bO5cuOQSePXVMM/SO+/A7rsnHZWUUzY1iUFAnZTHdQhLmYqI/J57uBCuWbMwU+tdd8GoUUoQeSqbCf5qu/uPJQ/c/Uczq5PpABEpUrNnQ/v28O67YfTS44/DzjsnHZVUQjY1iaVm9sucvGa2J/BTfCGJSN5ZvTqsMd2iBYwYEUYvDRmiBFEAsqlJdAReMrOvo8dbAieXdZCZ1QaGAutF53nZ3TuZWQPgBaAJMBs4yd0XljtyEakepk4NF8WNGAGHHQaPPgrbbpt0VFJFyqxJuPtHwC7ARUAHoKm7j83itZcDB7r77kBL4FAz25uwDsUgd9+J0N+hdSlE8tHKlXD77WGVuGnT4NlnQx+EEkRByTRV+IHu/l4pk/ztlM0Ef+7uQElfRq3o5sAxQLto+9PAEODa8ocuIokZOzbUHiZMgJNOChPybbZZ0lFJDDI1N7UF3gOOKmVfVhP8mVlNYCywI9Dd3UeZ2ebuPhfA3eeaWamfLDNrD7QH2FbfTESqh59+gltuCfMsbb459O0Lxx6bdFQSo7RJIuo/qAG85e4vVuTF3X010NLM6gF9zaxFOY7tAfQAaN26tVfk/CJShYYOhfPOg+nTQy3i3nuhXr2ko5KYZeyTiC6gu6SyJ3H3HwjNSocC35rZlgDRz+8q+/oiEqPFi6FDB2jbFlatgoEDw9BWJYiikM0Q2HfN7Coz28bMGpTcyjrIzBpGNQjMbH3gL8CnwGvAmdHTzgT6VSx0EYndm2+GYa2PPgodO4YJ+Q46KOmoJIeyGQJbsqbExSnbHNi+jOO2BJ6O+iVqAC+6+xtmNgJ40czOBf4HnFjOmEUkbt9/D//8Jzz3XJh3afhw2HvvpKOSBGSTJJq6+8+pG6JrIDJy94lAq1K2zwf0VUSkOnKHl14Kcy4tXAj/+hfccAOst17SkUlCsmluGp7lNhHJZ19/DccdByefDI0bh2GunTsrQRS5TNdJbAE0AtY3s1b8ulLdxvx2wj8RyWfuYWW4K68MK8X93/+F/od1smlokEKX6VPwV+AsYGvg/pTti4EbYoxJRHJl1iw4/3x4770weunxx2HHHZOOSqqRTNdJPE3oeD7B3V/JYUwiEreSCfluuglq1oTHHgvXQNTQOmTyW9l8Ij40s15m9haAmTWLRiaJSD6aPBn22w+uuAIOOACmTAnTeytBSCmy+VQ8CQwAtooef0aYGVZE8smKFXDrrdCqFcyYEZYTff112HrrpCOTaiybJLFpNC3HGgB3XwWsjjUqEalaH30Ee+4JnTrBCSeE6b1POQXMyj5Wilq2iw5tQriAjmi670WxRiUiVWPZMrjqqnAh3IIF0K8f9O4NDRsmHZnkiWzGuF1BmEpjBzP7EGgI/C3WqESk8oYMCZ3RM2eGPod77oG6dZOOSvJMmUnC3T82s7bAHwjXSkxz95WxRyYiFbNoEVxzDfToAdtvH4a3HnBA0lFJnsp0Md3aiw2V2DmbRYdEJAFvvAEXXghz54aL4269Fero2lepuEw1iZeB8dENfr3iGrJcdEhEcmTePLj88tDf0KIFvPoqtGmTdFRSADIliROAk4HdCNN593b3GTmJSkSy4w59+sBll4Vmps6d4brrYN11k45MCkTa0U3u3tfd/05YxnQmcJ+ZDYv6J0QkaV99BUcfHYaybr89fPxxmLVVCUKqUDZDYH8mDHldDGwAlDlNuIjEaM2aMI1Gs2YwaBDcf39Y76FF1qsDi2QtU8f1AcA/gDbAQKCbu4/JVWAiUooZM8KEfEOGhBFLPXvCDjskHZUUsEx9EoOAicAwYD3gDDM7o2Snu18Wc2wiUmLVKujaFW6+OTQn9ewJ556rK6YldpmSxNk5i0JE0vvkk5AQPvoo9EE8/DA0apR0VFIkypoqXESSsnw53HFHuNWvDy+8ACeeqNqD5JSWnhKpjkaNCrWHyZPhtNPggQdg002TjkqKkCaQF6lOli4N6zzss0+47qF/f3j2WSUISYxqEiLVxXvvhZFLs2bBRRfBXXfBxhsnHZUUuTKThJk1BM4HmqQ+393PiS8skSLyww9w9dVhfemddoL334f99086KhEgu5pEP+ADwrUSWmxIpCr16xdqDd9+GxJF586w/vpJRyXyi2ySRB13vzb2SESKyXffhfmWXngBdtsNXnsNWrdOOiqR38mm4/oNMzs89khEioE7PPccNG0KfftCly4wZowShFRb2dQkLgduMLPlwErClOHu7upREymPL78Maz28+WZYTrRXrzD/kkg1ls3KdBvlIhCRglUyId8114T7XbvCJZdAzZpJRyZSpkwT/O3i7p+a2R6l7Xf3j+MLS6RATJ8e1pkeOhT+8pewpOh22yUdlUjWMtUkrgDaA/eVss+BA2OJSKQQrFoVpvDu1Alq14YnnoCzztKUGpJ3Ms3d1D76qRXURcpjwoQwpcbYsXDccdC9O2y5ZdJRiVSIpuUQqSrLl4epvFu3Dp3UL70U1ppWgpA8pmk5RKrCiBGh9jB1Kpx5ZmhqatAg6ahEKk01CZHK+PFH6NgR9tsvTM739tvw1FNKEFIwsqpJmFkjoDG/nbtpaFxBieSFgQPDhHyzZ4chrXfcARtpxLgUlmwm+LsbOBmYwq9zNzmgJCHFaeFCuOqqMGJp553hgw/gT39KOiqRWGRTkzgW+IO7L485FpHqr29f6NAB5s2D66+Hf/0rDHEVKVDZJIlZQC1ASUKK1zffwKWXwssvQ8uWYWqNVq2SjkokdtkkiWXAeDMbREqicPfLMh1kZtsAzwBbAGuAHu7ezcxuIaxPMS966g3u/mYFYheJn3tYGa5jR1i2LPQ7XHUV1KqVdGQiOZFNkngtupXXKuBKd//YzDYCxprZu9G+B9z93gq8pkjufPEFXHABDBgQRi89/jjsskvSUYnkVDYT/D1tZusCO0ebprn7yiyOmwvMje4vMbOpQKPKBCuSE2vWwMMPw3XXhcf//nfoh6ihEeNSfMr81JtZO2A60B14GPjMzMq1tqKZNQFaAaOiTZeY2UQze8LM6qc5pr2ZjTGzMfPmzSvtKSJV79NPw9Khl14aRixNnhyGtypBSJHK5pN/H3CIu7d19/2BvwIPZHsCM9sQeAXo6O6LgUeAHYCWhJpGaRMI4u493L21u7du2LBhtqcTqZiVK+HOO0On9JQp4YK4t96Cxo2TjkwkUdn0SdRy92klD9z9MzPLqtcuet4rwPPu/mp0/Lcp+3sCb5QvZJEqNm4cnHMOjB8PJ5wADz0EW2yRdFQi1UI2NYkxZtbLzNpFt57A2LIOMjMDegFT3f3+lO2ps50dB0wqb9AiVeLnn8O1Dn/8Yxji+sorYYirEoTIL7KpSVwEXAxcRli6dCihb6Is+wGnA5+Y2fho2w3AP8ysJeGq7dnABeWKWKQqDBsWJuT77DM4+2y47z6oX2r3mEhRy2Z003Lg/uiWNXcfRkgqa9M1EZKcJUvghhvCGg+NG8M778DBBycdlUi1lWn50hfd/SQz+4Twrf833H23WCMTqWoDBkD79mGth8sug9tugw03TDoqkWotU03i8ujnkbkIRCQ2CxbAP/8JzzwDTZvChx/CPvskHZVIXkjbcR1dDAfQwd2/SL0BHXITnkgluIeO6KZN4T//gZtuCiOZlCBEspbN6KbSGmwPq+pARKrU3LlhOOuJJ8I228CYMdClC6y3XtKRieSVTH0SFxFqDDuY2cSUXRsBw+MOTKRC3OHJJ+HKK8MQ17vvhiuugHW0Uq9IRWT6y/kP8BZwJ3BdyvYl7r4g1qhEKuLzz0PH9MCB8Oc/hwn5dt657ONEJK1MfRKL3H020A1YkNIfsdLM9spVgCJlWr0aunWDFi1g5MgwOd+QIUoQIlUgmz6JR4AfUx4vjbaJJG/KlFBr6NgR2rYNE/JddJEm5BOpItn8JZm7/3KdhLuvIbsrtUXis3JluM6hVSuYNi0sDNS/P2y7bdKRiRSUbJLELDO7zMxqRbfLCUuaiiRj7Fho3RpuvhmOPTbUJk47Day0C/xFpDKySRIXAvsCc4CvgL2A9nEGJVKqn36Ca6+FNm1g3jzo2xdeeAE23zzpyEQKVjZzN30H/D0HsYikN3QonHceTJ8efv7f/0G9eklHJVLwykwSZvYkpc/ddE4sEYmkWrw4LCP6yCOw/fZheOtBByUdlUjRyKYDOnVRoNqENSC+jicckRRvvgkXXghz5oS5l7p0gQ02SDoqkaKSTXPTK6mPzaw3MDC2iES+/z4kheeeg2bNYPhw2EuX5ogkoSKDyXcCNM5Qqp576Ihu1gz69IFOneDjj5UgRBKUTZ/EEkKfhEU/vwGujTkuKTZffw0dOkC/fmF466BBsOuuSUclUvSyaW7aKBeBSJFyh1694KqrYPnyMGqpY0dNyCdSTWSaBXaPTAe6+8dVH44UlVmz4Pzz4b33wpQajz8OO+6YdFQikiLT17X7op+1gdbABEKT027AKOBP8YYmBWv1anjwQbjxRqhVCx57LFz7oPmWRKqdTLPAHuDuBwBfAHu4e2t33xNoBczIVYBSYCZPhv32C2s8HHRQeNy+vRKESDWVzV/mLu7+SckDd58EtIwtIilMK1bArbeGCflmzgzLib72Gmy9ddKRiUgG2fQOTjWzx4HnCKObTgOmxhqVFJaPPoJzz4VPPoFTToGuXaFhw6SjEpEsZFOTOBuYDFwOdASmRNtEMlu2DK6+GvbeGxYsCDWH559XghDJI9kMgf3ZzLoTrrJ2YJq7r4w9MslvQ4aEzuiZM+GCC8Ja03XrJh2ViJRTmTUJM2sHTAceAh4GPjOz/eMNS/LWokUhKRxwQHj83nvw6KNKECJ5Kps+ifuAQ9x9GoCZ7Qz0BvaMMzDJQ6+/HpYOnTs3XBzXuTPUqZN0VCJSCdn0SdQqSRAA7v4ZUCu+kCTvzJsXOqSPPhoaNICRI8OV00oQInkvm5rEWDPrBTwbPT4VGBtfSJI33KF3b7jssrDuQ+fOYe2HdddNOjIRqSLZJIkLgYuBywhXXA8l9E1IMfvqq9C09MYbYZbWXr2gefOkoxKRKpYxSZhZDWCsu7cA7s9NSFKtrVkDPXuGoa2rV8MDD8Cll0LNmklHJiIxyNgn4e5rgAlmpvUjJKwvfeCBYbW4Nm3CxXEdOypBiBSwbJqbtgQmm9loYGnJRnc/OraopHpZtSpcJX3zzbDeemG21nPOAbOkIxORmGWTJDrHHoVUXxMnhik1xowJo5ceeQS22irpqEQkRzKtJ1Gb0Gm9I/AJ0MvdV+UqMEnY8uVwxx3hVr9+WFb0xBNVexApMplqEk8DK4EPgMOAZoT5m6TQjRwZag9TpsBpp4Wmpk02SToqEUlApiTRzN13BYiukxidm5AkMUuXwk03Qbdu0KgR9O8Phx+edFQikqBMSeKXSfzcfZWpmaGwDRoUlhL9/HPo0AHuvBM23jjpqEQkYZmGwO5uZouj2xJgt5L7Zra4rBc2s23MbLCZTTWzyWZ2ebS9gZm9a2bTo5/1q6owUgE//BBma/3LX2CddeD996F7dyUIEQEyL19a0903jm4bufs6Kfez+Q+yCrjS3ZsCewMXm1kz4DpgkLvvBAyKHksS/vtfaNYMnnoKrr0WJkyA/TXBr4j8KraFhd19rrt/HN1fQljNrhFwDKFTnOjnsXHFIGl8+y2cdBIcdxxsthmMGgV33QXrr590ZCJSzeRk9XkzawK0AkYBm7v7XAiJBNgszTHtzWyMmY2ZN29eLsIsfO7w7LOh9tCvH9x2W1hadE/N+i4ipYs9SZjZhsArQEd3L7Mvo4S793D31u7euqGWu6y8//0PjjgCzjgD/vAHGDcObrwRamnWdxFJL9YkYWa1CAnieXd/Ndr8rZltGe3fEvguzhiK3po14Srp5s1Dp3S3bvDBB6E2ISJShkxXXC8hrGn9u12Al9V5bWHMbC9gqrunziD7GnAmcFf0s195g5YsffZZGLn0wQdw8MHQowc0aZJ0VCKSR9ImCXffqJKvvR9wOvCJmY2Ptt1ASA4vmtm5wP+AEyt5HlnbqlVw//3QqRPUrg1PPglnnqkpNUSk3LKZ4A8AM9sMqF3y2N3/l+n57j6MUOsozUHZnlfKacKEMEPrxx+H0Uvdu8OWWyYdlYjkqTL7JMzsaDObDnwOvA/MBt6KOS4pr59/DlNqtG4Nc+bAyy/Dq68qQYhIpWTTcd2FcDHcZ+6+HaEW8GGsUUn5DB8OrVrB7bfDqaeGiflOOCHpqESkAGSTJFa6+3yghpnVcPfBQMt4w5Ks/PgjXH45/OlPsGwZvP12uHq6QYOkIxORApFNn8QP0bUOHwDPm9l3hCk3JEnvvAPt24frH0om5NuosmMNRER+K5uaxDHAT0BH4G1gJnBUjDFJJgsXwtlnw1//GkYuDR0KDz2kBCEisSgzSbj7UmBT4FBgPtAnan6SXHv11XAR3LPPwvXXw/jxoalJRCQm2YxuOomw4NCJwEnAKDP7W9yBSYpvvoG//S10Rm+xBYweHZYVrV277GNFRCohmz6JG4E/uvt3AGbWEBgIvBxnYEKYkO/pp+GKK0LH9O23w9VXa74lEcmZbJJEjZIEEZlPjmaPLWqzZ8MFF4QO6v32g8cfh112SToqESky2SSJt81sANA7enwyupguPmvWhKukr78+TKPx0ENw0UVQQ3lZRHKvzCTh7leb2fHAnwjTbPRw976xR1aMPv00TMj34Ydh9NJjj0HjxklHJSJFLO3XUzPb0cz2A3D3V939Cnf/JzDfzHbIWYTFYOXK0BG9++7haumnnoK33lKCEJHEZWrD6AosKWX7smifVIVx46BNm7AA0NFHw9SpmrFVRKqNTEmiibtPXHuju48BmsQWUbH4+efQ7/DHP4Yhrq+8Ai+9BJtvnnRkIiK/yNQnkWkQ/vpVHUhRGTYMzj03LAp09tlw331Qv37SUYmI/E6mmsRHZnb+2hujxYLGxhdSAVuyBC65BP78Z1ixIgxvfeIJJQgRqbYy1SQ6An3N7FR+TQqtgXWB42KOq/AMGBAm5PvyyzBz6223wYYbJh2ViEhGmZYv/RbY18wOAFpEm/u7+3s5iaxQzJ8frph+5hlo2jQMb91nn6SjEhHJSjbXSQwGBucglsLiHjqjL74YFiyAm28OI5jWWy/pyEREspb1GtdSDnPnhuTQty/suWfoe9h996SjEhEpN831UJXcQ0d0s2bhYrh77oGRI5UgRCRvqSZRVT7/PHRMDxwI++8PPXvCzjsnHZWISKWoJlFZq1dDt27QokWoNTz8MAwerAQhIgVBNYnKmDIlXBQ3ciQcdhg8+ihsu23SUYmIVBnVJCpixQro0gVatYLp08Nyov37K0GISMFRTaK8xowJtYeJE+Hkk+HBB2GzzZKOSkQkFqpJZOunn+Caa2CvveD77+G//4U+fZQgRKSgqSaRjfffh/PPD01L558fhrbWq5d0VCIisVNNIpPFi8PSoe3ahVFMAwdCjx5KECJSNJQk0nnzTWjePCSFK64IfRAHHZR0VCIiOaUksbbvv4fTToMjjoC6dWH48LDewwYbJB2ZiEjOKUmUcIcXXghTarz4InTqBB9/HDqqRUSKlDquAebMgQ4d4LXXwnKivXrBrrsmHZWISOKKuybhHuZYatYM3n0X7r0XRoxQghARiRRvTWLmzDCcdfDgMHqpZ0/YccekoxIRqVaKryaxejXcf3+oLYwdC489BoMGKUGIiJSiuGoSkyaFKTVGj4Yjj4RHHoGtt046KhGRaqs4ahIrVkDnzrDHHjBrFvznP6GTWglCRCSjwq9JjB4dag+TJsEpp0DXrtCwYdJRiYjkhdhqEmb2hJl9Z2aTUrbdYmZzzGx8dDs8rvMDcNttsM8+sHAhvP46PP+8EoSISDnE2dz0FHBoKdsfcPeW0e3NGM8PO+wQRjBNnhz6IEREpFxia25y96Fm1iSu18/KP/4RbiIiUiFJdFxfYmYTo+ao+umeZGbtzWyMmY2ZN29eLuMTEZFIrpPEI8AOQEtgLnBfuie6ew93b+3urRuqH0FEJBE5TRLu/q27r3b3NUBPoE0uzy8iIuWT0yRhZlumPDwOmJTuuSIikrzYOq7NrDfQDtjUzL4COgHtzKwl4MBs4IK4zi8iIpUX5+im0oYV9YrrfCIiUvWKY1oOERGpECUJERFJy9w96RjKZGbzgC8qePimwPdVGE6+Kebyq+zFq5jLn1r2xu5eqWsI8iJJVIaZjXH31knHkZRiLr/KXpxlh+Iuf1WXXc1NIiKSlpKEiIikVQxJokfSASSsmMuvshevYi5/lZa94PskRESk4oqhJiEiIhWkJCEiImnlZZJIszRqSzMbGS2LOsbM2kTbT01ZLnW8ma2J5o/CzPY0s0/MbIaZPWhmllCRslbOstcys6ejMk41s+tTjin0sq9rZk9GZZxgZu1Sjsm7skPa8u9uZiOi8rxuZhun7Ls+KuM0M/tryva8K395ym5mm5jZYDP70cweWut18q7sUO7yH2xmY6PtY83swJRjyl9+d8+7G7A/sAcwKWXbO8Bh0f3DgSGlHLcrMCvl8WhgH8CAt0qOr8638pQdOAXoE92vQ5hUsUmRlP1i4Mno/mbAWKBGvpY9Q/k/AtpG988BukT3mwETgPWA7YCZQM18LX85y74B8CfgQuChtV4n78pegfK3AraK7rcA5lSm/HlZk3D3ocCCtTcDJd+i6gJfl3LoP4De8Mu05Ru7+wgP794zwLGxBFyFyll2BzYws3WA9YEVwOIiKXszYFB03HfAD0DrfC07pC3/H4Ch0f13gROi+8cQviAsd/fPgRlAm3wtf3nK7u5L3X0Y8HPqk/O17FDu8o9z95K/g8lAbTNbr6Llj20W2AR0BAaY2b2EZrR9S3nOyYQ/HoBGwFcp+76KtuWjjpRe9pcJ5Z1LqEn8090XmFlrCr/sE4BjzKwPsA2wZ/RzDYVTdghrshwN9ANOJJQRQplGpjyvpJwrKZzypyt7OoX0Nw/Zlf8EYJy7LzezCpU/L2sSaVxE+Ce4DfBP1pqW3Mz2Apa5e0mbXmltcfk6Hjhd2dsAq4GtCE0OV5rZ9hRH2Z8g/BGMAboCw4FVFFbZITQzXGxmY4GNCLVFSF/OQip/urKnU0hlhzLKb2bNgbv5dd2eCpW/kGoSZwKXR/dfAh5fa//fiZqaIl8BW6c83prSm6jyQbqynwK87e4rge/M7EOgNfABBV52d19FSBoAmNlwYDqwkMIpO+7+KXAIgJntDBwR7fqK336zLClnwXzuM5Q9nYIpO2Quv5ltDfQFznD3mdHmCpW/kGoSXwNto/sHEv4hAGBmNQjVsT4l29x9LrDEzPaOevjPIFTb8lG6sv8PONCCDYC9gU+LoexmVicqM2Z2MLDK3acUWNkxs82inzWAm4BHo12vAX+P2qK3A3YCRhdS+TOUvVSFVHZIX34zqwf0B6539w9Lnl/h8ifda1/Bnv7ehHb2kvbVcwmjGcYS2qJHAXumPL8dMLKU12lNaNebCTxEdAV6db6Vp+zAhoRv15OBKcDVRVT2JsA0YCowkDBlct6WPUP5Lwc+i253pZYFuDEq4zRSRrHkY/krUPbZhI7eH6PnN8vXspe3/ISEsRQYn3LbrKLl17QcIiKSViE1N4mISBVTkhARkbSUJEREJC0lCRERSUtJQkRE0lKSEClFdG3JMDM7LGXbSWb2dpJxieSahsCKpGFmLQjXmbQCahLGmx/qv17BWp7Xqunuq6s2QpH4KUmIZGBm9xAuTNog+tmYMOX8OsAt7t7PzJoAz0bPAbjE3YdbWMOiE+EiqJbu3iy30YtUnpKESAbR1B4fEyZPewOY7O7PRVMfjCbUMhxY4+4/m9lOQG93bx0lif5ACw/TdYvknUKa4E+kyrn7UjN7gTC9w0nAUWZ2VbS7NrAtYf6ohyyseLga2DnlJUYrQUg+U5IQKdua6GbACe4+LXWnmd0CfAvsThgMkrrYzdIcxSgSC41uEsneAODSknWBzaxVtL0uMNfd1wCnEzq5RQqCkoRI9roAtYCJ0YL0XaLtDwNnmtlIQlOTag9SMNRxLSIiaakmISIiaSlJiIhIWkoSIiKSlpKEiIikpSQhIiJpKUmIiEhaShIiIpLW/wOTDfE6XD1rjAAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"y_predict = regr.predict(X)\n",
"plt.plot(X, y_predict, color='red')\n",
"plt.title(\"French Coal Production Linear Regression\")\n",
"plt.xlabel(\"Year\")\n",
"plt.ylabel(\"Coal Production in Metric Tons\")"
]
},
{
"cell_type": "code",
"execution_count": 9,
"id": "4d365538",
"metadata": {},
"outputs": [],
"source": [
"X_future = np.array(range(1918, 2022))\n",
"X_future = X_future.reshape(-1, 1)"
]
},
{
"cell_type": "code",
"execution_count": 10,
"id": "4143b866",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"Text(0.5, 0, 'Year')"
]
},
"execution_count": 10,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAckAAAEWCAYAAAANV2yLAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAAA6pElEQVR4nO3deZyW8/7H8denfadNKpKtyJZkX5LooGTPekTTToQQfp2E42RLOJUiTraOhNNMSaUkS5IkJFFpT6V93+bz++O6htuYe5qpmfuae+b9fDzmMfe1f77X9rnW72XujoiIiPxVsagDEBERKaiUJEVEROJQkhQREYlDSVJERCQOJUkREZE4lCRFRETiKPBJ0swmmVm7AhBHXTNzMysRwbTzZR6Y2QNm9lJej3cP0xxjZm0SOc3CyMxuMLNxUceRW+E2dET4+wUz65mAad5sZp9m0/1yM1tsZpvM7MT8jkeSyx6TpJktMLOt4QqU8VcrEcHtDTOrZGb9zGxRGOvcsLlaPk83dj6tMLNXzKxCfk4zN8zsXDNbEtvO3R9z9/xIvnF3Su5+kbsPzetp7o1wh705XGZLzayvmRWPOq6ccPc33L15Xo83XE/Sw3my0czmmNkteT0dAHfv5O6P5CCm/D5Qfgq4zd0ruPuMfJzOX5hZKTMbEe4/3MzOzdS9dHgwscLM1phZmpnVDrvVybRf3hSO4+6Y4a83s4Xhev4/M6uSTSw1zSzVzJaF46mbqXttMxsZxrHEzDpl6j44XF/SzezmTN3MzB4Nt7P14TI9JptYYrfN1WY2wcyu2fMc/X34v+zv9lZOzyQvCVegjL9lmQJK+NlVVsysFDABOAa4EKgEnAGsBk5JQAiXuHsFoBFwMvB/WcRYIOZVUbKHeX5CuMyaANcAbRM8/YJoWThPKgH3AS+aWYPMPSVhueI5BJiVVYcElfFT4Ebg1yy63QGcDhwP1ALWAc8DuPui2P0ycByQDrwDECahQcDfgRrAFmBANnGkAx8AV8bp/jrwSziuFsBjZtY0pvtMoAvwdRbDXk2wbZ0NVAGmAK9lEwv8sW3WB/4D/NvMeu1hmLzn7tn+AQuA87No78CtwM/AL2G7lsA3BAvyc+D4TOPpDnwLrAfeAsrEdL80HHYDMA+4MGw/CXgE+AzYCIwDqsWJtR2wAqiQTXmODse5jmDDaBXTrQUwI4xhMfBQTLe6YZlL5GQ+AU8Co7KZV+2BucAaIBWoFTPsBcCP4Xz6N/Ax0C7s9hDwery4CFbAV4BlwFrgf0B5YCvBRrAp/KuVxbhahfNkXTiPjs7p8ss0L24GPo3TbVJMWW4m2EE8Fcb6C3BRTL/7AUOA5cBS4FGgeNjtcGAiwQHQb8AbwP6Z4r0vjHd7VsstnG9HxDQPB/rHNGe3PjcK15WNwNvh/Hg07HYusCSc/q8EO4NiQA+CdXt1OK0qYf9lCHZAq8NpTQNqxMyj+eF0fgFuyGoeExwMTguXzTTgjEzzPKfb0LnAkkztVgFXhdP8DHiGYL19FCgdLr9FBNveC0DZmGHvCZffMoKd5O/znGDH92h2+wDgn8BuYBvBevvvsN+jgPFhHHOA1jHjqUqwTW0AvgzL/pf1MYx9UxjTZmBevHWHPW8b94T9byZYZ2sAY8L5/SFQOQf72iXAuZnaDQSeyLSPmhNn+F7ARzHNjwFvxjQfDuwAKu4hjhLhPKkb065C2K56TLvBwGtZDP8pcHOmdvcBw2OajwG2ZRPDn7bNsN1V4XpQNWy+BZgdzuP5QMewfbz93SkEyXkdwTr5b6DUHpdLDhbcAuInyfEEO+WyBDuNlcCpQHGgTThs6ZjxfBkGWyUsXKew2ykEG/cFBDuT2sBRMRv4PKBeOJ1JQJ84sf4XGJpNWUoSJKYHgFLAeeEMrh+zgzgujOF4go3+srBbXXKYJIGDCTaoR+LMq/MIduyNCDbU54HJYb/VCDbuq8J47wR2kfMkOZpgh105HL5JNju/38cVzt/N4TIoCdwbzqtSe1p+WcyLm8l5ktxJcMBQHOhMsDO1sPv/CI6EywMHhNPP2BCOCGMtDVQHJgP9Mi2Pb8JlUTZOLLE77KMINpw7w+a46zPBurOQ4Ci/JHAFwc4nNknuAh4P+y8LdAO+AA4K2w0ChoX9dwTSgHLhtE4iOIsrT7AuZKyfNYFjMs/jcHmsJThjKAFcFzZXjZnnOd2GziVcTwi2g8vDZVQ/nOYuoGs4nbJAP4KEVAWoGJbjX+HwFxJsQ8eGZXmTOEmSPe8D2sXEWJ7gIPaWMI5GBNtTxrz5L8FBSPlw2kuJsz5mtUMm07pDzraNLwgSY22C9eZr4MRwWU8EeuVgX5tVkmxMcGBSK1w/3iRmPc/U7zxikhMwErgvUz+bgJP2EEdWSbJi2O6AmHYvAjOyGD6rJHlIOE/qhfPwCeB/OV0mMfvvXYQH0gQHDIcDRnAlaAvQKJv93UnAaWH56hLsw7rtcbnkYMEtCGfsuvDvfzGFOC+mv4GESSGm3Rz+2EkvAG6M6fYE8EL4exDwTDY71f+Lae4CfBCn3/HE2fjD7mcTHNkXi2k3jJgzxkz998uIi5wlyYz5tJDgskbZOPNqCH8+OqxAsCOqC9wEfBHTzQg2nj0mSYKdaDpZHLXGWWl+HxfQkz8f6RUj2Lmcu6fll8W0bibnSXJuTLdyYVkOJNjhbOfPZyXXEXOknGm8lxGzwYbxtt3Duu0ESWhz+HsYfxzUxV2fgXPCeWMx3T7lz0lyB3++UjIbaBbTXDNc5iUIzrD+dKYa9lM+XJ+uJFOi589J8u/Al5m6TyHcUZG7bejccB1aR3CW9g1wbcw0F2VaNzcDh8e0O50/rpa8TMz2SLCDjJck97QPiE2S1wCfZOpnEMGZVPFwvh4V0+0xcp8k28Y052TbuCGm+zvAwJjmrmSTEGL6yypJViJYL50gQcwgvAKRqb+zCfY/FWLaTSDTgWxs3NnE8ZckGbOOP09w5aNRuH785ayWrJNkKeDZmHL8Ahya02US0/7X2Hmdqdv/gDti1uMl8cYf9tMNeG9PyyWn19svc/cPs2i/OOb3IUAbM+sa064UwRFQhthr7ltiuh0MvJ/N9DMPF++BmNUEO594agGL3T09pt1CgqM/zOxUoA/B0WcpgqPAt7MZX2bx5hP8eV7VIua6vbtvMrPVYRy1Yvt1dzez2GGzczCwxt3X5iLm2JgWxkw3PZxu7Zh+4i2/ffH7ON19i5lBsHyrEBw5Lg/bQbBzWgxgZgcAzxHsHCqG3TKXOyfzrRHBEfjVBMu+PEFyzm59dmCph1tanGmtcvdtMc2HAO+ZWey6t5vgYOA1gmX3XzPbn+DS64Puvjl8WKE7MMTMPgPudvcfM03rT8su9Pt6HcrpNgTBPcmD4nSLLWd1ggOb6THLyAgSVUZc0zPFFM+e9gGxDgFONbN1Me1KEMzH6uHv2Dizm248mbfXPW0bK2J+b82ieW8f4htIkJSqEhyQ3EtwGffUTP21Ad5x900x7TYRJNlYlYCNZnZ2OB6Ahe4e9yGaGDcA/QnmzXyCWxx/uVcdRy+C5zQOJlgXbwQmmtkx7r4lJyMws5IEy3dN2HxRON56BNt/OeC7bIavB/QlODsvR7CeTI/Xf4Z9fQUk807in+6+f8xfOXcfloPxLCY4bd5XHwJ/M7PycbovAw42s9hy1yE4uoLgUkYqcLC770dwf8XIG7HzahnBhg5AGG/VMI7lBCtSRjeLbSbYUMrFNB8Y83sxUCXc0WY3/axkjiljukvjDpG/FhMkq2ox61OlmI35XwRlOt7dKxFsdJmX1Z7KHPQUGE5w9vWPmOnHW5+XA7UtJjPw52WU1bQXE1wmih1fGXdf6u473b23uzcguLfYkuCKAu4+1t0vIDj4+5HgEldmf1p2odj1Oi/Flus3ggRwTEyZ9vPgYQvItC6HMcWT3T4gq3n5caZ5WcHdOxPcP92Vi+nGk932msht4wTgP+6+xt23E5zJnRL7tL6ZlSU4yMv81PiscPiM/g4jOPD/yd0/8T8e+slJgsTdF7p7S3ev7u6nEuyzvsxFOd5y9yXuvsvd/0NwSyinSRaCe9a7gC/NrDTBGftTBPfv9yc4yMrYJrPa9gcSbENHhvuMB8jB/j0v35N8EehkZqeGj/uWN7MWZlYxB8MOAW4xs2ZmVix81PiovYjhNYIN6B0zOyocV1UL3ge8GJhKeDRmZiXDx60vIbiHAcEZyRp332ZmpwDX70UMOfEmQXkbhgv7MWCquy8guKd4jJldET5Zdzt/ToTfAOeEj3/vB9yf0cHdlxMcHQ4ws8phGc8JO68AqobDZGU40CJcBiWBuwmS1Od7WUYzszKxf7kZOCzLOOBpC17rKWZmh5tZk7CXioSXty14JP6evYwzVh+gg5kdSPbr8xSCs8DbzKyEmV3Knp+efgH4p5kdAmBm1cPhMLOmZnacBa+fbCC4XLjbzGqYWavwIGp7WN7dWYz7faCeBY/7lwjPPhsAo/ZxfmQrvCLzIvBMeGaf8ZrA38JehgM3m1kDMytHcNQfT3b7gBXAYTH9jiIo79/DdbykmZ1sZke7+27gXeAhMysXPpXbZh+Lmtfbxp9Y8JpHxvZRKtxeMnbe04CbzGy/cNpdCM70f4sZxeUEl8c/yjTqN4BLzOzscB16GHjX3TdmE0sZgkQKEBsXZna0mVW04LWVG4HmBGdmGd1Lhf0bUDIsR0aOmQZcHa7Txczs7/zxjMie5k8VM8s4i33c3Vfzx5W+VcAuC84qY1+Jymp/V5Fg+9oUrlud9zRtyMMk6e5fETyA8W+Cy15zCe5h5GTYLwluwj9DcPP+Y/56ZJyT8WwHzic4WhjPH0+3VSNIQjsInlK7iOAoeABwU8zlqy7Aw2a2keCMYnhuY8hhnBMI7nO8Q3C0fThwbdjtN/649LcaOJLgxn3GsOMJHsz5luBSQeYd4d8JdrI/EjxA0C0c7keCexvzzWydZXrX1d3nEJyNPU8wby4heKVlx14W8wyCs4zf/yz3j9PfRLAx/ECwTo3gj8vpvQkula4nOLB4dy/j/J27f0ew7t2T3foczpMrgBSCndONBMthezajf5bgKsW4cP36gj8umR0Ylm0Dwb3LjwkuuRYj2CEvI7jE1IRgHc0c92qCs8+7CdaZe4GWmXak+eU+gnnzhZltILiaUz+MawzBff2JYT8T441kD/uAZ4GrzGytmT0X7uSbE2wzywgu32U8JAVwG8HlzV8J7nu+si8FzIdtI7M5BNtIbWBs+Duj7N0Jnuj8mSAhXEyQFGO1AV7NdPkfd58FdCJIlisJksRf1p9MthIcjEGwD9ka0+1vBJdZ14bjvdDdV8V0Hxf2fwbBk69bCe7fQ7B8ZvLH0+J3Ale6+7psYplpZpsI1p12BA/V/SMs20aCE4jhYTzXE2xfGWXPan/XPexvI8HB3Vt7mBfAH08Risg+MLOpBA8y7dMOWUQKlgJfLZ1IQWRmTczswPDyZhuCV4Y+iDouEclbhaXGDJFEq09wqacCwdOxV4X3UUWkENHlVhERkTh0uVVERCSOInm5tVq1al63bt2owxARSSrTp0//zd2rRx1HIhXJJFm3bl2++uqrqMMQEUkqZrY3tRclNV1uFRERiUNJUkREJA4lSRERkTiUJEVEROJQkhQREYlDSVJERCQOJUkREZE4lCRFRIqA3em76TulL58u+jTqUJKKkqSISCH3/crvOePlM7h73N28PevtqMNJKkmVJM3sTjObZWbfm9mw8MvXVcxsvJn9HP6vHHWcIiIFwY7dO+g9qTeNBjVi/tr5DLtyGP0u7Bd1WEklaZKkmdUm+BJ1Y3c/FihO8GXyHsAEdz8SmBA2i4gUaV8u/ZKTBp/EQx8/ROtjWjP71tlce+y1mFnUoSWVpEmSoRJAWTMrAZQDlgGXAkPD7kOBy6IJTUQkelt2bqH7uO6cPuR01m5dS9p1abx+xetUK1ct6tCSUtJUcO7uS83sKWARsBUY5+7jzKxGxsdu3X25mR2Q1fBm1gHoAFCnTp1EhS0ikjAf/fIR7dLaMX/tfDqd1InHL3icSqUrRR1WUkuaM8nwXuOlwKFALaC8md2Y0+HdfbC7N3b3xtWrF6kvvYhIIbd+23o6pHXgvFfPo5gVY1KbSQxsOVAJMg8kzZkkcD7wi7uvAjCzd4EzgBVmVjM8i6wJrIwySBGRREqbk0an0Z34ddOvdD+9O72b9qZcyXJRh1VoJFOSXAScZmblCC63NgO+AjYDbYA+4f+RkUUoIpIgqzav4vYPbue/3/+X4w44jpHXjqRxrcZRh1XoJE2SdPepZjYC+BrYBcwABgMVgOFmlkKQSK+OLkoRkfzl7gz7fhi3j7mdDds30Pvc3vQ4qwelipeKOrRCKWmSJIC79wJ6ZWq9neCsUkSkUFu8fjGdR3dm9M+jOe2g0xjSaggNqjeIOqxCLamSpIhIUZTu6QyePph7x9/Lbt/NM397hq6ndKV4seJRh1boKUmKiBRgP6/+mfZp7fl44cecf9j5DG45mEMrHxp1WEWGkqSISAG0K30Xz0x5hn9M+geli5dmSKsh3NLwFtWYk2BKkiIiBczMX2eSkprC9OXTueyoy+h/cX9qVawVdVhFkpKkiEgBsX3Xdh6d/Ch9PutDlbJVeOuqt7i6wdU6e4yQkqSISAEwZfEUUlJTmP3bbG464Sb6Nu9L1XJVow6ryFOSFBGJ0OYdm3lw4oM8N/U5Dqp0EO9f/z4XHXlR1GFJSElSRCQiH87/kPZp7VmwbgG3nnwr/2r2LyqWrhh1WBJDSVJEJMHWbl1L93Hdefmbl6lXtR6Tb57M2YecHXVYkgUlSRGRBHpv9nt0eb8LqzavoseZPeh1bi/KlCgTdVgSh5KkiEgCrNi0gq5juvL2D2/T8MCGjL5+NI1qNoo6LNkDJUkRkXzk7rz27Wt0+6Abm3du5p/n/ZN7zriHksVLRh2a5ICSpIhIPlm0fhEdR3Xkg7kfcMbBZ/DSJS9xdPWjow5LckFJUkQkj6V7OgOnDaTHhB64O89d+BxdTu6iCsmTkJKkiEgemvPbHNqltePTRZ9ywWEXMPiSwdTdv27UYcleUpIUEckDO3fv5OkpT/PQpIcoW7Isr1z6Cm1OaKMq5ZJcpEnSzCoDB7v7t1HGISKyL7759RvajmzLjF9ncMXRV/Dvi/5NzYo1ow5L8kCxRE/QzCaZWSUzqwLMBF4xs76JjkNEZF9t27WNByc8SOPBjVm2cRkjrh7BO63fUYIsRKI4k9zP3TeYWTvgFXfvZWY6kxSRpPLZos9ol9aOH3/7kZsb3szTzZ+mStkqUYcleSzhZ5JACTOrCbQGRkUwfRGRvbZpxyZuH3M7Z79yNlt3bmXsjWN55dJXlCALqSjOJB8GxgKfuvs0MzsM+DmCOEREcmXs3LF0HNWRResXcdspt/FYs8eoUKpC1GFJPkp4knT3t4G3Y5rnA1cmOg4RkZxas3UNd429i6Ezh1K/an0+ueUTzqxzZtRhSQIkPEmaWXWgPVA3dvru3jbRsYiI7Mm7s9+ly+gu/LblNx446wF6NumpCsmLkCgut44EPgE+BHZHMH0RkT36ddOv3Pb+bbwz+x1OPPBEPrjxAxoe2DDqsCTBokiS5dz9vgimKyKyR+7O0JlDuWvsXWzZuYU+zfpw9xl3U6KY6l4piqJY6qPM7GJ3fz+CaYuIxLVg3QI6jurIuHnjOKvOWbx0yUvUr1Y/6rAkQlEkyTuAB8xsB7AzbOfuXimCWERE2J2+m/7T+vPAhAcwM/pf3J9OjTtRzKJ4S04Kkiiebq2Y6GmKiMQze9Vs2qW14/PFn3PhERcyqOUg6uxXJ+qwpICI5CK7mbUCzgkbJ7m7KhUQkYTauXsnT3z2BA9PfpgKpSrw6mWvcuPxN6pCcvmTKF4B6QOcDLwRtrrDzM5y9x6JjkVEiqavl39N25FtmbliJlc3uJrnL3qeGhVqRB2WFEBRnEleDDR093QAMxsKzACUJEUkX23duZXeH/fmqc+fonr56rzb+l0uP/ryqMOSAiyqZ5r3B9aEv/eLKAYRKUI+WfgJ7dLa8dPqn2jbsC1PNX+KymUrRx2WFHAJS5JmNs7dmwP/AmaY2UeAEdybvD9RcYhI0bJh+wbu//B+Bnw1gEP3P5Txfx/P+YedH3VYkiQSeSZZHcDdh5nZJIL7kgbc5+6/JjAOESkixvw8ho6jOrJkwxK6ndqNR897lPKlykcdliSRRCbJ/czsiizan2FmuPu7CYxFRAqx1VtWc+fYO3nt29doUL0Bn6d8zmkHnRZ1WJKEEpokgZYEZ4+ZOZBtkjSz+sBbMa0OA/4BvBq2rwssAFq7+9p9D1dEko27M+KHEdw25jbWbF1Dz3N68uDZD1K6ROmoQ5MklcgkuXBfvvTh7nOAhgBmVhxYCrxH8FTsBHfvY2Y9wmbVDStSxCzbuIxb37+V//34P06qeRLj/z6e42scH3VYkuQSmSTz8g3dZsA8d19oZpcC54bthwKTUJIUKTLcnZdnvMzd4+5m++7tPHH+E9x5+p2qkFzyRCLXor/n4biuBYaFv2u4+3IAd19uZgdkNYCZdQA6ANSpoyqnRAqD+Wvn0yGtAxN+mUCTQ5rw4iUvcmTVI6MOSwqRhNXe6+7f58V4zKwU0Ap4O5fTH+zujd29cfXq1fMiFBGJyO703fT7oh/HDTyOL5d+ycAWA5nYZqISpOS5ZLwecRHwtbuvCJtXmFnN8CyyJrAywthEJJ/NWjmLlNQUpi6dSosjW/BCyxc4qNJBUYclhVTCvwNjZuXN/vj+jJkVM7NyuRjFdfxxqRUgFWgT/m4DjNz3KEWkoNmxewePfPwIJw46kblr5vL65a+Tdl2aEqTkqyjOJCcA5wObwuZywDjgjD0NGCbTC4COMa37AMPNLAVYBFydp9GKSOSmLZ1GSmoK3638jmuPvZbnLnyO6uV120TyXxRJsoy7ZyRI3H1TTs8k3X0LUDVTu9UET7uKSCGzZecWHpr0EE9PeZoDKxzIyGtH0qp+q6jDkiIkiiS52cwaufvXAGZ2ErA1gjhEpACbtGAS7dPaM3fNXNo3as+TFzzJfmX0PQRJrCiSZDfgbTNbFjbXBK6JIA4RKYDWb1vPfR/ex6Dpgzi88uFMvGkiTQ9tGnVYUkQlPEm6+zQzOwqoT1DBwI/uvjPRcYhIwTPqp1F0GtWJ5ZuWc/fpd/Nw04cpVzI3z/WJ5K1EfirrPHefmEUl50eqgnORom3V5lV0G9uNN797k2MPOJZ3r3mXU2qfEnVYIgk9k2wCTAQuyaLbHis4F5HCx915a9ZbdB3TlfXb1vNQk4e4/+z7KVW8VNShiQAJTJLu3it8P3KMuw9P1HRFpGBaumEpXd7vQuqcVE6pfQpDWg3h2AOOjToskT9J6D1Jd083s9sAJUmRIird03np65e4Z/w97Ny9k6ebP80dp95B8WLFow5N5C+ieLp1vJl1J/gG5OaMlu6+JoJYRCSB5q6ZS/u09kxaMImmdZvy4iUvcniVw6MOSySuKJJkxjclb41p5wQfURaRQmhX+i6e/eJZen7Uk5LFS/LiJS+ScmIKZnn5BT2RvBdFkjza3bfFtjCzMhHEISIJ8N2K70hJTWHasmlcUu8SBrYYSO1KtaMOSyRHEl7BOfB5DtuJSBLbvms7vT7qRaPBjViwbgHDrhzGyGtHKkFKUknke5IHArWBsmZ2IkFFAgCVCCo5F5FCYuqSqaSkpjBr1SxuOO4G+l3Yj2rlqkUdlkiuJfJy69+Am4GDgL4x7TcADyQwDhHJJ1t2bqHnxJ70m9qPWhVrMeq6UbSo1yLqsET2WiLfkxwKDDWzK939nURNV0QSY+IvE2mf1p75a+fTuXFn+pzfh0qlK0Udlsg+ieKe5GdmNsTMxgCYWYPwW5AikoTWbVtH+9T2NHu1GcWsGJPaTGJAiwFKkFIoRJEkXwHGArXC5p8IvgwiIkkmdU4qxww4hpe/eZl7z7iXbzt9S5O6TaIOSyTPRJEkq4XV0qUDuPsuYHcEcYjIXlq5eSXXjriWS/97KdXKVWNqu6k8fsHjlC1ZNurQRPJUVB9drkpQgQBmdhqwPoI4RCSX3J03v3uTOz64gw3bN/DwuQ9z31n3qUJyKbSiSJJ3AanA4Wb2GVAduCqCOEQkFxavX0yn0Z14/+f3Oe2g0xjSaggNqjeIOiyRfBXFR5e/NrMm/PHR5Tn66LJIwZXu6QyePph7x9/Lbt/NM397hq6ndFWF5FIkJLIygcwfW85QTx9dFimYfl79M+3S2jF54WTOP+x8BrcczKGVD406LJGESeSZ5Ajgm/AP/qhxB/TRZZECZVf6LvpO6UuvSb0oXbw0Q1oN4ZaGt6hCcilyEpkkrwSuAY4HRgLD3H1uAqcvIjkw89eZpKSmMH35dC476jL6X9yfWhVr7XlAkUIoYa+AuPt77n4t0ASYBzxtZp+G9ydFJGLbd22n58SeNH6xMYs3LGb4VcN5t/W7SpBSpEXxdOs2glc+NgB1AH0mSyRiUxZPISU1hdm/zeamE26ib/O+VC1XNeqwRCKXyAd3mgLXAacAHwLPuvtXiZq+iPzVph2b+L+J/8dzU5/j4P0OZswNY7jwiAujDkukwEjkmeQE4FvgU6A0cJOZ3ZTR0d1vT2AsIkXe+Hnj6TCqAwvWLeDWk2/lX83+RcXSFaMOS6RASWSSvCWB0xKRONZuXUv3cd15+ZuXqVe1HpNvnszZh5wddVgiBVKiP5UlIhF6b/Z7dHm/C6s2r+L+s+7nH03+QZkSeixAJJ4oHtwRkQT7ddOvdB3TlRE/jKDhgQ0Zff1oGtVsFHVYIgWekqRIIebuvPbta3T7oBtbdm7hsfMeo/sZ3SlZvGTUoYkkBSVJkUJq4bqFdBzVkbHzxnLmwWfyUquXOKraUVGHJZJUEp4kzaw60B6oGzt9d2+b6FhECqN0T2fgtIH0mNADd+f5i56ny8ldKGZRfD5WJLlFcSY5EviE4F1JfWxZJA/N+W0O7dLa8emiT2l+eHMGtRxE3f3rRh2WSNKKIkmWc/f7IpiuSKG1c/dOnp7yNA9NeohyJcvxn0v/w00n3KQKyUX2URRJcpSZXezu7+d2QDPbH3gJOJbgyyFtgTnAWwSXbxcArd19bV4FK1LQzVg+g5TUFGb8OoOrGlzF8xc9z4EVDow6LJFCIYqbFHcQJMqtZrbBzDaa2YYcDvss8IG7HwWcAMwGegAT3P1Iglp9euRL1CIFzLZd23hwwoOc/OLJLNu4jHdav8PbV7+tBCmShxJ+Junue1XvlZlVAs4Bbg7HswPYYWaXAueGvQ0FJgG6nCuF2meLPiMlNYU5q+dwS8NbeLr501QuWznqsEQKnURWcH6Uu/9oZlm+wezuX+9hFIcBq4BXzOwEYDrBWWkNd18ejmO5mR0QZ/odgA4AderU2ctSiERr4/aNPDDhAfpP60+d/eow9saxND+8edRhiRRaiTyTvIsgST2dRTcHztvD8CWARkBXd59qZs+Si0ur7j4YGAzQuHFjz+lwIgXFuHnjaJ/WnsXrF9P1lK78s9k/qVCqQtRhiRRqiay7tUP4v+lejmIJsMTdp4bNIwiS5AozqxmeRdYEVu57tCIFx5qta7hr7F0MnTmU+lXr88ktn3BmnTOjDkukSEiat4vd/VdgsZnVD1s1A34AUoE2Ybs2BO9hihQKI34YwdH9j+aN797gwbMf5JtO3yhBiiRQslVL1xV4w8xKAfMJPr9VDBhuZinAIuDqCOMTyRPLNy7ntjG38e7sd2lUsxFjbxxLwwMbRh2WSJGTVEnS3b8BGmfRqVmCQxHJF+7Of775D3eNu4utO7fSp1kf7j7jbkoUS6pNVaTQiGTLM7PawCH8ue7WyVHEIlJQLFi3gA5pHRg/fzxn1TmLIa2GUK9qvajDEinSoqjg/HHgGoL7iRl1tzqgJClF0u703QyYNoD7J9yPmdH/4v50atxJFZKLFABRnEleBtR39+0RTFukQJm9ajbt0trx+eLPufCICxnUchB19tN7vCIFRRRJcj5QElCSlCJr5+6dPPHZEzw8+WEqlKrAa5e/xg3H3aAKyUUKmCiS5BbgGzObQEyidPfbI4hFJOG+Xv41bUe2ZeaKmbQ+pjXPX/Q8B5TPsqIoEYlYFEkyNfwTKVK27txK749789TnT3FA+QN475r3uOyoy6IOS0SyEUUF50PD9xwzHtub4+47Ex2HSCJ9svAT2qW146fVP9HuxHY82fxJ9i+zf9RhicgeRPF067kEX+tYABhwsJm10SsgUhht2L6BHh/2YOBXAzl0/0P58O8f0uwwvdYrkiyiuNz6NNDc3ecAmFk9YBhwUgSxiOSb939+n46jOrJ0w1K6ndqNR897lPKlykcdlojkQhRJsmRGggRw95/MrGQEcYjki9+2/MadY+/k9W9fp0H1Bryd8janHXRa1GGJyF6IIkl+ZWZDgNfC5hsIvg0pktTcneGzhtN1TFfWbltLz3N68uDZD1K6ROmoQxORvRRFkuwM3ArcTnBPcjIwIII4RPLMso3L6Dy6M6lzUmlcqzEftvqQ42scH3VYIrKPoni6dTvQN/wTSWruzpAZQ+g+rjvbd2/nyQuepNtp3VQhuUghkbAt2cyGu3trM/uOoK7WP3F3HXZLUpm/dj7t09oz8ZeJNDmkCS+1eokjqhwRdVgikocSebh7R/i/ZQKnKZLndqfv5rmpz/F/H/0fxa04g1oOol2jdqqQXKQQSliSdPfl4c8u7n5fbLfwyyD3/XUokYJl1spZpKSmMHXpVFoc2YIXWr7AQZUOijosEcknURz6XpBFu4sSHoVILuzYvYNHPn6EEwedyLy183jjijdIuy5NCVKkkEvkPcnOQBfgcDP7NqZTReDzRMUhklvTlk4jJTWF71Z+x7XHXstzFz5H9fLVow5LRBIgkfck3wTGAP8CesS03+juaxIYh0iObNm5hV4f9aLvF305sMKBjLx2JK3qt4o6LBFJoETek1wPrDezZ4E17r4RwMwqmtmp7j41UbGI7MmkBZNon9aeuWvm0qFRB5644An2K7Nf1GGJSIJFcU9yILAppnlz2E4kcuu3rafTqE40HdoUd2fiTRMZdMkgJUiRIiqKN57N3X9/T9Ld081Mb15L5Eb9NIpOozqxfNNyup/end5Ne1OuZLmowxKRCEVxJjnfzG43s5Lh3x3A/AjiEAFg1eZVXP/O9Vwy7BIql63MlJQpPNn8SSVIEYkkSXYCzgCWAkuAU4EOEcQhRZy7M+y7YTQY0IARP4yg97m9md5hOqfUPiXq0ESkgIii7taVwLWJnq5IrCUbltBldBfSfkrjlNqnMKTVEI494NiowxKRAibhSdLMXiHrulvbJjoWKXrSPZ2Xvn6Je8bfw87dO3m6+dPcceodFC9WPOrQRKQAiuKBmVExv8sAlwPLIohDipi5a+bSPq09kxZM4rxDz2Nwy8EcXuXwqMMSkQIsisut78Q2m9kw4MNExyFFx670XfT7oh89P+pJqeKlePGSF0k5MQUzizo0ESngCsKrF0cCdaIOQgqn71Z8R0pqCtOWTaNV/VYMuHgAtSvVjjosEUkSUdyT3EhwT9LC/7+iL4BIHtu+azuPffIYj336GJXLVOa/V/6X1se01tmjiORKFJdbKyZ6mlK0TF0ylZTUFGatmsUNx91Avwv7Ua1ctajDEpEklMivgDTKrru7f52oWKRw2rxjMz0/6km/L/pRu1JtRl03ihb1WkQdlogksUSeST4d/i8DNAZmElxyPR6YCpyVwFikkJn4y0Tap7Vn/tr5dG7cmT7n96FS6UpRhyUiSS5hNe64e1N3bwosBBq5e2N3Pwk4EZibqDikcFm3bR3tU9vT7NVmFLfifHzzxwxoMUAJUkTyRBRPtx7l7t9lNLj792bWMII4JMmN/HEknUd3ZuXmldx35n30atKLsiXLRh2WiBQiUSTJ2Wb2EvA6wdOtNwKzczKgmS0ANgK7gV3u3tjMqgBvAXWBBUBrd1+b92FLQbFy80puH3M7b816i+NrHE/adWmcVOukqMMSkUIoigrObwFmAXcA3YAfwnY51dTdG7p747C5BzDB3Y8EJoTNUgi5O298+wYN+jfgvR/f45Gmj/BV+6+UIEUk30TxCsg2M+tPUMuOA3Pcfec+jPJS4Nzw91BgEnrvstBZvH4xnUZ34v2f3+e0g05jSKshNKjeIOqwRKSQi6IygXMJktkCgqdbDzazNu4+OQeDOzDOzBwY5O6DgRruvhzA3Zeb2QFxptuB8JNcdeqogp9kke7pDPpqEPd+eC/pns6zFz7LrSffqgrJRSQhorgn+TTQ3N3nAJhZPWAYkJNrZme6+7IwEY43sx9zOtEwoQ4GaNy48V++QiIFz8+rf6ZdWjsmL5zMBYddwKCWgzi08qFRhyUiRUgUSbJkRoIEcPefzKxkTgZ092Xh/5Vm9h5wCrDCzGqGZ5E1gZX5ErUkzK70XfSd0pdek3pRpkQZXrn0Fdqc0EZVyolIwkWRJKeb2RDgtbD5BmD6ngYys/JAMXffGP5uDjwMpAJtgD7h/5H5ErUkxMxfZ9I2tS1fL/+ay4+6nP4X96dmxZpRhyUiRVQUSbITcCtwO8E9ycnAgBwMVwN4LzybKAG86e4fmNk0YLiZpQCLgKvzJWrJV9t2bePRyY/y+GePU6VsFYZfNZyrGlyls0cRiVRCk6SZFQOmu/uxQN/cDOvu84ETsmi/GmiWNxFKFKYsnkLb1Lb8+NuP3HTCTfRt3peq5apGHZaISGLfk3T3dGCmmenxUmHTjk10+6AbZ758Jlt2bmHMDWMYetlQJUgRKTCiuNxaE5hlZl8CmzNaunurCGKRiIyfN54OozqwcN1CupzchX81+xcVS+sraiJSsESRJHtHME0pINZuXUv3cd15+ZuXqV+1PpNvmcxZdfQBGBEpmBL5PckyBA/tHAF8Bwxx912Jmr5E773Z79Hl/S6s2ryKHmf2oNe5wSseIiIFVSLPJIcCO4FPgIuABgT1t0oh9+umX+k6pisjfhhBwwMbMvr60TSqme03uEVECoREJskG7n4cQPie5JcJnLZEwN15dear3Dn2Trbs3MJj5z1G9zO6U7J4juqOEBGJXCKT5O+VmLv7Lr3/VrgtXLeQjqM6MnbeWM44+AyGtBrCUdWOijosEZFcSWSSPMHMNoS/DSgbNhvg7q5PyRcC6Z7OgGkD6PFh8MWy5y96ni4nd6GYRfFVNhGRfZOwJOnu+mxDITfntzmkpKbw2eLPaH54cwa3HMwh+x8SdVgiInstildApJDZuXsnT33+FL0/7k25kuVUIbmIFBpKkrJPZiyfQUpqCjN+ncFVDa7i+Yue58AKB0YdlohInlCSlL2ybdc2Hv74YZ747Amql6/OO63f4Yqjr4g6LBGRPKUkKbn26aJPaZfajjmr53Bzw5vp27wvlctWjjosEZE8l8gadzYCnlUn9HRrUti4fSMPTHiA/tP6c8j+hzDuxnFccPgFUYclIpJvEvl0q2qvTmIfzP2AjqM6snj9Yrqe0pV/NvsnFUpViDosEZF8FdnlVjM7APi94k53XxRVLBLfmq1ruHPsnbw681WOrnY0n7X9jNMPPj3qsEREEiLhSdLMWgFPA7WAlcAhwGzgmETHItkb8cMIbn3/VtZsXcODZz9Iz3N6UrpE6ajDEhFJmCjOJB8BTgM+dPcTzawpcF0EcUgcyzcu57Yxt/Hu7HdpVLMR424cxwkHnhB1WCIiCRdFktzp7qvNrJiZFXP3j8zs8QjikEzcnf988x/uGncXW3dupU+zPtx9xt2UKKaHoEWkaIpi77fOzCoQfDLrDTNbCei7khH7Ze0vdBzVkfHzx3POIefw4iUvUq9qvajDEhGJVBRJ8lJgG9ANuAHYD3g4gjgE2J2+m/7T+vPAhAcoZsUYcPEAOjbuqArJRUSIIEm6+2YzqwGcDKwGxrj76kTHITB71WxSUlOYsmQKFx1xEYNaDuLg/Q6OOiwRkQIj4acLZtaa4IPLVwOtgalmdlWi4yjKdu7eyaOTH6XhoIb8tPonXr/8dUZfP1oJUkQkkygutz4InOzuKwHMrDrwITAigliKnOnLptM2tS3frviWa465hucueo4Dyh8QdVgiIgVSFEmyWEaCDK0mgjPaombrzq08NOkhnpryFDXK12DktSNpVb9V1GGJiBRoUSTJD8xsLDAsbL4GGBNBHEXG5IWTaZfajp/X/EzKiSk81fwp9i+zf9RhiYgUeFE8uHOPmV0BnEVQuflgd38v0XEUBRu2b6DHhz0Y+NVADt3/UD78+4c0O6xZ1GGJiCSNRH4F5Aighrt/5u7vAu+G7c8xs8PdfV6iYikK3v/5fTqO6sjSDUvpdmo3Hj3vUcqXKh91WCIiSSWR9wL7ARuzaL8l7CZ54Lctv3HjuzfS4s0WVCpdiSkpU3jmwmeUIEVE9kIiL7fWdfdvM7d096/MrG4C4yiU3J3hs4bTdUxX1m5byz/O+QcPnP2AKiQXEdkHiUySZbLpVjZhURRCyzYuo8voLoycM5LGtRozodUEjqtxXNRhiYgkvURebp1mZu0ztzSzFGB6AuMoNNydl75+iQb9GzB23lievOBJpqRMUYIUEckjiTyT7Aa8Z2Y38EdSbAyUAi5PYByFwvy182mf1p6Jv0ykySFNeKnVSxxR5YiowxIRKVQSliTdfQVwRvj9yGPD1qPdfWKiYigMdqfv5rmpz/HgxAcpUawEL7R4gfYntVeF5CIi+SCK9yQ/Aj5K9HQLg1krZ5GSmsLUpVNpWa8lA1sM5KBKB0UdlohIoZV0px9mVtzMZpjZqLC5ipmNN7Ofw/+Vo44xr+3YvYOHP36YEwedyLy183jzijdJvTZVCVJEJJ8lXZIE7gBmxzT3ACa4+5HAhLC50Ji2dBonDT6JXpN6cfUxV/NDlx+47rjrMLOoQxMRKfSSKkma2UFAC+ClmNaXAkPD30OByxIcVr7YsnML3cd157Qhp7F261pSr03ljSveoHr56lGHJiJSZERRwfm+6AfcC1SMaVfD3ZcDuPtyM0v67z5NWjCJdqntmLd2Hh1P6sjj5z/OfmX2izosEZEiJ2nOJM2sJbDS3ffqnUoz62BmX5nZV6tWrcrj6PLG+m3r6ZjWkaZDmwIw8aaJvNDyBSVIEZGIJNOZ5JlAKzO7mKD2nkpm9jqwwsxqhmeRNYGVWQ3s7oOBwQCNGzf2RAWdU2lz0ug8ujPLNy2n++nd6d20N+VKlos6LBGRIi1pziTd/X53P8jd6wLXAhPd/UYgFWgT9tYGGBlRiHtl1eZVXP/O9bT6bysql63MlJQpPNn8SSVIEZECIJnOJOPpAwwPq7dbBFwdcTw54u4M+34Yt4+5nQ3bN9D73N70OKsHpYqXijo0EREJJWWSdPdJwKTw92ogqb4kvGTDEjqP7syon0Zxau1TGdJqCMcccEzUYYmISCZJmSSTVbqn8+L0F7ln/D3sSt9F3+Z9uf3U2ylerHjUoYmISBaUJBNk7pq5tEttx8cLP+a8Q89jcMvBHF7l8KjDEhGRbChJ5rNd6bvo90U/en7Uk1LFS/HiJS+ScmKKaswREUkCSpL56LsV39E2tS1fLfuKS+tfyoAWA6hVsVbUYYmISA4pSeaD7bu289gnj/HYp49RuUxl3rrqLa5ucLXOHkVEkoySZB6bumQqbVPb8sOqH7jx+Bvp97d+VC1XNeqwRERkLyhJ5pHNOzbT86Oe9PuiH7Ur1Wb09aO5+MiLow5LRET2gZJkHpgwfwLt09rzy7pf6Ny4M33O70Ol0pWiDktERPaRkuQ+WLdtHd3HdWfIjCEcWeVIPr75Y8455JyowxIRkTyiJLmXRv44ks6jO7Ny80ruO/M+ejXpRdmSZaMOS0RE8pCSZC6t2LSC2z+4neGzhnN8jeNJuy6Nk2qdFHVYIiKSD5Qkc2HC/Am0HtGaTTs28WjTR7n3zHspWbxk1GGJiEg+UZLMhSOrHknjWo3p97d+HF396KjDERGRfKYkmQt19qvD2BvHRh2GiIgkSNJ8dFlERCTRlCRFRETiUJIUERGJQ0lSREQkDiVJERGROJQkRURE4lCSFBERiUNJUkREJA5z96hjSDgzWwUszOPRVgN+y+NxFlRFqaxQtMqrshZOeVXWQ9y9eh6MJ2kUySSZH8zsK3dvHHUciVCUygpFq7wqa+FUlMqa13S5VUREJA4lSRERkTiUJPPO4KgDSKCiVFYoWuVVWQunolTWPKV7kiIiInHoTFJERCQOJUkREZE4lCSzYWYvm9lKM/s+pt0JZjbFzL4zszQzqxS2v8DMpoftp5vZeTHDnBS2n2tmz5mZRVGe7OSmrDHd65jZJjPrHtOu0JXVzI4Pu80Ku5cJ2xeqsppZSTMbGrafbWb3xwyTDGU92Mw+CmOfZWZ3hO2rmNl4M/s5/F85Zpj7wzLNMbO/xbQv0OXNbVmTff8UKXfXX5w/4BygEfB9TLtpQJPwd1vgkfD3iUCt8PexwNKYYb4ETgcMGANcFHXZ9qWsMd3fAd4GuhfWsgIlgG+BE8LmqkDxQlrW64H/hr/LAQuAuklU1ppAo/B3ReAnoAHwBNAjbN8DeDz83QCYCZQGDgXmJcuy3YuyJvX+Kco/nUlmw90nA2syta4PTA5/jweuDPud4e7LwvazgDJmVtrMagKV3H2KB2vkq8Bl+R58LuWmrABmdhkwn6CsGe0KY1mbA9+6+8xw2NXuvruQltWB8mZWAigL7AA2JFFZl7v71+HvjcBsoDZwKTA07G0of8R+KcFBwXZ3/wWYC5ySDOXNbVmTff8UJSXJ3PseaBX+vho4OIt+rgRmuPt2ghV3SUy3JWG7ZJBlWc2sPHAf0DtT/4WurEA9wM1srJl9bWb3hu0LY1lHAJuB5cAi4Cl3X0MSltXM6hKcPU0Farj7cgiSC3BA2FttYHHMYBnlSqry5rCssQrL/ikhlCRzry1wq5lNJ7jMsSO2o5kdAzwOdMxolcU4kuW9m3hl7Q084+6bMvVfGMtaAjgLuCH8f7mZNaNwlvUUYDdQi+Dy491mdhhJVlYzq0BwK6Cbu2/Irtcs2nk27QucXJQ1o//CtH9KiBJRB5Bs3P1HgktwmFk9oEVGNzM7CHgPuMnd54WtlwAHxYziIGAZSSCbsp4KXGVmTwD7A+lmto1gYy1sZV0CfOzuv4Xd3ie4x/c6ha+s1wMfuPtOYKWZfQY0Bj4hScpqZiUJ1sM33P3dsPUKM6vp7svDy4srw/ZL+POVoIxyJcU2m8uyFrr9U6LoTDKXzOyA8H8x4P+AF8Lm/YHRwP3u/llG/+Elj41mdlr41NhNwMhEx7034pXV3c9297ruXhfoBzzm7v8ujGUFxgLHm1m58F5dE+CHQlrWRcB5FigPnAb8mCxlDWMbAsx2974xnVKBNuHvNvwReypwbXhv7lDgSODLZChvbstaGPdPCRP1k0MF+Q8YRnB/ZifBEVcKcAfBk2Q/AX34o9ai/yO4n/NNzN8BYbfGBPeB5gH/zhimIP3lpqyZhnuIPz/dWujKCtxI8LDD98AThbWsQAWCp5VnAT8A9yRZWc8iuFT4bcw2eDHBE8kTgJ/D/1VihnkwLNMcYp7qLOjlzW1Zk33/FOWfqqUTERGJQ5dbRURE4lCSFBERiUNJUkREJA4lSRERkTiUJEVEROJQkhRJgPDdw0/N7KKYdq3N7IMo4xKR7OkVEJEEMbNjCd5DPBEoTvCu2oX+R+0nuRlXcXffnbcRikhmSpIiCRRW5bcZKB/+PwQ4jqCKyIfcfWRYYfVrYT8At7n752Z2LtCLoHKAhu7eILHRixQ9SpIiCRRW9/Y1QaXio4BZ7v56WG3YlwRnmQ6ku/s2MzsSGObujcMkORo41oNPO4lIPlMF5yIJ5O6bzewtYBPQGrjEzLqHncsAdQgqmP63mTUk+CpHvZhRfKkEKZI4SpIiiZce/hlwpbvPie1oZg8BK4ATCB6u2xbTeXOCYhQR9HSrSJTGAl3Dry9gZieG7fcDlrt7OvB3god8RCQCSpIi0XkEKAl8a2bfh80AA4A2ZvYFwaVWnT2KREQP7oiIiMShM0kREZE4lCRFRETiUJIUERGJQ0lSREQkDiVJERGROJQkRURE4lCSFBERieP/AcOxoS+/jJBtAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"future_predict = regr.predict(X_future)\n",
"plt.plot(X_future, future_predict, color='green')\n",
"plt.title(\"French Coal Production Linear Regression Predicted from 1870-1918 Data\")\n",
"plt.ylabel(\"Coal Production in Metric Tons\")\n",
"plt.xlabel(\"Year\")"
]
},
{
"cell_type": "code",
"execution_count": 11,
"id": "d6a6b40a",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"[0.43477551] -797.8321632653062\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAckAAAEWCAYAAAANV2yLAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAABCn0lEQVR4nO3debxV4/7A8c+3eZBoVCmZK0PUMWaIhMic6YZcQ6cSRVKmGz9SyE1XxMElJCVFUmmgqCgVukVRmmfN83S+vz+eZ9fudPY5+0xr7XPO9/16ndfZe6291vrutdda3/U861nPElXFGGOMMYcqEnYAxhhjTKKyJGmMMcbEYEnSGGOMicGSpDHGGBODJUljjDEmBkuSxhhjTAwJnyRFZIKI3JcAcdQWERWRYiEsO0/WgYg8ISLv5PZ8M1nmKBFpFeQyCyIRaSkiY8KOI6v8PnSCf/2miDwdwDLvFpFJGYy/QUSWishWETkzr+Mx+UumSVJEFonIDr8BRf6qBxFcdojI4SLyqogs8bHO9+8r5fFyo9fTahF5T0QOy8tlZoWINBaRZdHDVPUFVc2L5BvzoKSqzVS1f24vMzv8AXub/82Wi8i/RaRo2HHFQ1UHqOrluT1fv52k+nWyRUTmicg/c3s5AKraRlWfiyOmvD5R7gW0V9XDVPXnPFzOIUSkhIgM8ccPFZHGacaX9CcTq0VkvYh8KSI1/LhaaY7LW/08OkVN/w8RWey3889FpEIGsVQTkeEissLPp3aa8TVE5AsfxzIRaZNmfIrfXlJF5O4040REnvf72Sb/m56SQSzR++Y6ERkvIrdmvkb3T3/I8S674i1JXuM3oMjfijQBBV66So+IlADGA6cAVwKHA+cD64CzAwjhGlU9DGgAnAU8lU6MCbGuCpNM1nl9/5tdDNwK3BPw8hPRCr9ODge6AG+LSL20H8qH3yuWY4A56Y0I6DtOAu4AVqUzrgNwHnA6UB3YCLwGoKpLoo/LwGlAKvAZgE9CbwF3AlWB7cAbGcSRCowGboox/iNgoZ/X1cALInJJ1PhfgXbAzHSmvRm3b10IVAB+AD7MIBY4sG+eDLwP9BWRbplMk/tUNcM/YBFwWTrDFXgA+BNY6Ic1B37B/ZBTgNPTzOdRYBawCRgElIoaf52fdjOwALjSD58APAdMBrYAY4BKMWK9D1gNHJbB96nr57kRt2NcGzXuauBnH8NS4JmocbX9dy4Wz3oCXgZGZLCu7gfmA+uB4UD1qGmbAnP9euoLTATu8+OeAT6KFRduA3wPWAFsAD4HygI7cDvBVv9XPZ15XevXyUa/jurG+/ulWRd3A5NijJsQ9V3uxh0gevlYFwLNoj5bHngXWAksB54HivpxxwPf4E6A/gYGAEekibeLj3dXer+bX28nRL0fDLwe9T6j7bmB31a2AJ/69fG8H9cYWOaXvwp3MCgCdMVt2+v8sir4z5fCHYDW+WX9BFSNWkd/+eUsBFqmt45xJ4M/+d/mJ+D8NOs83n2oMbAszbC1QAu/zMlAb9x2+zxQ0v9+S3D73ptA6ahpO/vfbwXuILl/neMOfM9ndAwAugP7gJ247bav/2wdYKyPYx5wS9R8KuL2qc3ANP/dD9kefexbfUzbgAWxth0y3zc6+89vw22zVYFRfn2PA46M41i7DGicZlg/4KU0x6h5MabvBnwb9f4F4OOo98cDu4FymcRRzK+T2lHDDvPDKkcNSwE+TGf6ScDdaYZ1AQZHvT8F2JlBDAftm35YC78dVPTv/wn87tfxX0CyHx7reHc2LjlvxG2TfYESmf4ucfxwi4idJMfiDsqlcQeNNcA5QFGglZ+2ZNR8pvlgK/gv18aPOxu3czfFHUxqAHWidvAFwEl+OROAnjFi/QTon8F3KY5LTE8AJYBL/Qo+OeoAcZqP4XTcTn+9H1ebOJMkUBO3Qz0XY11dijuwN8DtqK8B3/nPVsLt3C18vA8De4k/SX6FO2Af6ae/OIOD3/55+fW7zf8GxYHH/Loqkdnvl866uJv4k+Qe3AlDUaAt7mAqfvznuDPhskAVv/zIjnCCj7UkUBn4Dng1ze/xi/8tSseIJfqAXQe34zzs38fcnnHbzmLcWX5x4EbcwSc6Se4FXvSfLw10BH4EjvbD3gIG+s8nA18CZfyyGuJKcWVx20Jk+6wGnJJ2HfvfYwOuxFAMuN2/rxi1zuPdhxrjtxPcfnCD/41O9svcCzzol1MaeBWXkCoA5fz36OGnvxK3D53qv8vHxEiSZH4MuC8qxrK4k9h/+jga4PanyLr5BHcSUtYvezkxtsf0Dsik2XaIb9/4EZcYa+C2m5nAmf63/gboFsexNr0kmYQ7Manut4+PidrO03x2AVHJCfgC6JLmM1uBhpnEkV6SLOeHVYka9jbwczrTp5ckj/Hr5CS/Dl8CPo/3N4k6fu/Fn0jjThiOBwRXE7QdaJDB8a4hcK7/frVxx7COmf4ucfxwi/yK3ej/Po/6EpdGfa4fPilEDZvHgYP0IuCOqHEvAW/6128BvTM4qD4V9b4dMDrGZ8cSY+f34y/EndkXiRo2kKgSY5rPvxqJi/iSZGQ9LcZVa5SOsa7e5eCzw8NwB6LawF3Aj1HjBLfzZJokcQfRVNI5a42x0eyfF/A0B5/pFcEdXBpn9vuls6y7iT9Jzo8aV8Z/l6NwB5xdHFwquZ2oM+U0872eqB3Wx3tPJtu24pLQNv96IAdO6mJuz8BFft1I1LhJHJwkd3NwTcnvQJOo99X8b14MV8I6qKTqP1PWb083kSbRc3CSvBOYlmb8D/gDFVnbhxr7bWgjrpT2C3Bb1DKXpNk2twHHRw07jwO1Jf8lan/EHSBjJcnMjgHRSfJW4Ps0n3kLV5Iq6tdrnahxL5D1JHlP1Pt49o2WUeM/A/pFvX+QDBJC1OfSS5KH47ZLxSWIn/E1EGk+dyHu+HNY1LDxpDmRjY47gzgOSZJR2/hruJqPBn77OKRUS/pJsgTQJ+p7LASOjfc3iRq+Knpdpxn3OdAhajteFmv+/jMdgWGZ/S7x1rdfr6rj0hm+NOr1MUArEXkwalgJ3BlQRHSd+/aocTWBkRksP+10sRrErMMdfGKpDixV1dSoYYtxZ3+IyDlAT9zZZwncWeCnGcwvrVjrCQ5eV9WJqrdX1a0iss7HUT36s6qqIhI9bUZqAutVdUMWYo6OaXHUclP9cmtEfSbW75cT++epqttFBNzvWwF35rjSDwN3cFoKICJVgP/gDg7l/Li03zue9dYAdwZ+M+63L4tLzhltzwosV7+nxVjWWlXdGfX+GGCYiERve/twJwMf4n67T0TkCFzV65Oqus03VngUeFdEJgOdVHVummUd9Nt5+7drL959CNw1yaNjjIv+npVxJzYzon4jwSWqSFwz0sQUS2bHgGjHAOeIyMaoYcVw67Gyfx0dZ0bLjSXt/prZvrE66vWOdN5ntxFfP1xSqog7IXkMV417TprPtQI+U9WtUcO24pJstMOBLSJyoZ8PwGJVjdmIJkpL4HXcuvkLd4njkGvVMXTDtdOoidsW7wC+EZFTVHV7PDMQkeK433e9f9/Mz/ck3P5fBvhfBtOfBPwbVzovg9tOZsT6fERObwFJe5DorqpHRP2VUdWBccxnKa7YnFPjgCtEpGyM8SuAmiIS/b1r4c6uwFVlDAdqqmp53PUVIXdEr6sVuB0dAB9vRR/HStyGFBkn0e9xO0qZqPdHRb1eClTwB9qMlp+etDFFlrs85hR5aykuWVWK2p4Oj9qZe+C+0+mqejhup0v7W2X2nd2HnMG40te/opYfa3teCdSQqMzAwb9Resteiqsmip5fKVVdrqp7VPVZVa2Hu7bYHFejgKp+rapNcSd/c3FVXGkd9Nt50dt1bor+Xn/jEsApUd+pvLrGFpBmW/YxxZLRMSC9dTkxzbo8TFXb4q6f7s3CcmPJaH8Nct+oD7yvqutVdReuJHd2dGt9ESmNO8lL22p8jp8+8rnjcCf+f6jq93qg0U88CRJVXayqzVW1sqqegztmTcvC9xikqstUda+qvo+7JBRvkgV3zXovME1ESuJK7L1w1++PwJ1kRfbJ9Pb9frh96ER/zHiCOI7vuXmf5NtAGxE5xzf3LSsiV4tIuTimfRf4p4g0EZEivqlxnWzE8CFuB/pMROr4eVUUdz/gVcBU/NmYiBT3za2vwV3DAFciWa+qO0XkbOAf2YghHh/jvu8Z/sd+AZiqqotw1xRPEZEbfcu6hzg4Ef4CXOSbf5cHHo+MUNWVuLPDN0TkSP8dL/KjVwMV/TTpGQxc7X+D4kAnXJKaks3vKCJSKvovKxP77zIGeEXcbT1FROR4EbnYf6QcvnpbXJP4ztmMM1pPoLWIHEXG2/MPuFJgexEpJiLXkXnr6TeB7iJyDICIVPbTISKXiMhp4m4/2YyrLtwnIlVF5Fp/ErXLf9996cx7JHCSuOb+xXzpsx4wIofrI0O+RuZtoLcv2UduE7jCf2QwcLeI1BORMriz/lgyOgasBo6L+uwI3Pe902/jxUXkLBGpq6r7gKHAMyJSxrfKbZXDr5rb+8ZBxN3mEdk/Svj9JXLw/gm4S0TK+2W3w5X0/46axQ246vFv08x6AHCNiFzot6H/A4aq6pYMYimFS6QA0XEhInVFpJy421buAC7Hlcwi40v4zwtQ3H+PSI75CbjZb9NFRORODrQRyWz9VBCRSCn2RVVdx4GavrXAXnGlyuhbotI73pXD7V9b/bbVNrNlQy4mSVWdjmuA0RdX7TUfdw0jnmmn4S7C98ZdvJ/IoWfG8cxnF3AZ7mxhLAdat1XCJaHduFZqzXBnwW8Ad0VVX7UD/k9EtuBKFIOzGkOccY7HXef4DHe2fTxwmx/3Nweq/tYBJ+Iu3EemHYtrmDMLV1WQ9kB4J+4gOxfXgKCjn24u7trGXyKyUdLc66qq83Clsddw6+Ya3C0tu7P5Nc/HlTL2/0nWm9PfhdsZfsNtU0M4UJ3+LK6qdBPuxGJoNuPcT1X/h9v2Ome0Pft1ciNwL+7gdAfud9iVwez74Gopxvjt60cOVJkd5b/bZty1y4m4KtciuAPyClwV08W4bTRt3Otwpc9OuG3mMaB5mgNpXumCWzc/ishmXG3OyT6uUbjr+t/4z3wTayaZHAP6AC1EZIOI/Mcf5C/H7TMrcNV3kUZSAO1x1ZurcNc938vJF8yDfSOtebh9pAbwtX8d+e6P4lp0/olLCFfhkmK0VsAHaar/UdU5QBtcslyDSxKHbD9p7MCdjIE7huyIGncFrpp1g5/vlaq6Nmr8GP/583EtX3fgrt+D+31+5UBr8YeBm1R1Ywax/CoiW3Hbzn24RnX/8t9tC64AMdjH8w/c/hX57ukd7x71n9uCO7kblMm6AA60IjTG5ICITMU1ZMrRAdkYk1gSvls6YxKRiFwsIkf56s1WuFuGRocdlzEmdxWUHjOMCdrJuKqew3CtY1v466jGmALEqluNMcaYGKy61RhjjImhUFa3VqpUSWvXrh12GMYYk6/MmDHjb1WtHHYcQSqUSbJ27dpMnz497DCMMSZfEZHs9F6Ur1l1qzHGGBODJUljjDEmBkuSxhhjTAyWJI0xxpgYLEkaY4wxMViSNMYYY2KwJGmMMcbEYEnSGGMKgb2pe/n3D//m+8Xfhx1KvpKvkqSIPCwic0RktogM9A/1rCAiY0XkT///yLDjNMaYRDJjxQzOeeccOo3pxJDfhoQdTr6Sb5Kkf/r8Q0CSqp4KFMU9dLUrMF5VTwTG+/fGGFPobdm1hY6jO3L2O2ezcstKPr35U1698tWww8pX8lu3dMWA0iKyByiDeyr540BjP74/MAH3tHRjjCm0vpj7Be1HtWf55uW0SWpDjyY9KF+qfNhh5Tv5Jkmq6nIR6QUsAXYAY1R1jIhUjTzHT1VXikiV9KYXkdZAa4BatWoFFbYxxgRq+eblPDjqQYbNHcZpVU5jcIvBnFfzvLDDyrfyU3XrkcB1wLFAdaCsiNwR7/SqmqKqSaqaVLlyoerE3hhTCOxL3cdrU1+j7ut1GT1/ND2b9GRG6xmWIHMo35QkgcuAhaq6FkBEhgLnA6tFpJovRVYD1oQZpDHGBO2XVb+QPCKZacuncfnxl9Pv6n4cd+RxYYdVIOSbkiSumvVcESkjIgI0AX4HhgOt/GdaAV+EFJ8xxgRq2+5tdB7TmaSUJBZtXMTHN37M6JajLUHmonxTklTVqSIyBJgJ7AV+BlKAw4DBInIvLpHeHF6UxhgTjJF/jqTdV+1YvGkx9ze4n56X9aRC6Qphh1Xg5JskCaCq3YBuaQbvwpUqjTGmwFu5ZSUdv+7I4DmDqVupLt/d/R0XHnNh2GEVWPkqSRpjTGGVqqmkzEih67iu7Ny7k/9r/H881ugxShYrGXZoBZolSWOMSXCz18wmeUQyU5ZO4dJjL6Xf1f04qeJJYYdVKFiSNMaYBLVjzw6e++45Xp7yMuVLlqf/9f258/Q7cW0XTRAsSRpjTAIau2Asbb9qy4INC2hVvxW9Lu9FpTKVwg6r0LEkaYwxCWTNtjU88vUjDPjfAE6scCLj7xrPpcdeGnZYhZYlSWOMSQCqyn9//i+dx3Zm6+6tPH3R0zxx4ROUKlYq7NAKNUuSxhgTsrl/zyV5RDLfLf6OC2tdyFvN36Ju5bphh2WwJGmMMaHZuXcnPSf1pMekHpQtXpZ3rnmHf575T4pIfuoMrWCzJGmMMSGYsGgCySOS+WPdH7Q8rSWvXP4KVQ+rGnZYJg1LksYYE6B129fReWxn3vvlPY478ji+vuNrLj/+8rDDMjFYkjTGmACoKh/O+pBOYzqxcedGHr/gcZ666CnKFC8TdmgmA5YkjTEmj/257k/aftWW8QvHc97R5/FW87c4reppYYdl4mBJ0hhj8sjufbt5efLLPPfdc5QqVoo3rnqD5KRka5iTj1iSNMaYPDBpySSSRyTz29rfuLnezfS5sg/VylULOyyTRZYkjTEmF23YsYGu47qSMjOFWuVr8dU/vuKqE68KOyyTTZYkjTEmF6gqg+YMouPojqzdvpZO53Xi2cbPUrZE2bBDMzlgSdIYY3Jo4YaFtBvZjtHzR5NUPYlRLUdxZrUzww7L5IJQk6SIHAnUVNVZYcZhjDHZsWffHnr/2JtnJjxD0SJF6XNlHx446wGKFikadmgmlwSeJEVkAnCtX/YvwFoRmaiqjwQdizHGZNfUZVNpPaI1s1bP4rqTr+O1Zq9Rs3zNsMMyuSyMdsjlVXUzcCPwnqo2BC4LIQ5jjMmyTTs30X5ke8579zzWbV/HsFuH8fltn1uCLKDCqG4tJiLVgFuAJ+OdSEROBgZFDToO+BfwgR9eG1gE3KKqG3IrWGOMAdcwZ+jvQ3lo9EOs3LKSB89+kOcvfZ5yJcuFHZrJQ2GUJP8P+BqYr6o/ichxwJ+ZTaSq81T1DFU9A2gIbAeGAV2B8ap6IjDevzfGmFyzZNMSrvvkOlp82oIqZasw9b6p9GnWxxJkIRB4SVJVPwU+jXr/F3BTFmfTBFigqotF5DqgsR/eH5gAdMl5pMaYwm5v6l76TuvLU988haL0atqLDud2oFgRuzGgsAij4U5l4H5c9ej+5avqPVmYzW3AQP+6qqqu9PNYKSJVYiy3NdAaoFatWlkP3BhTqMxYMYPWI1ozc+VMrj7xal6/6nWOOeKYsMMyAQvjdOgL4HtgHLAvqxOLSAlc69jHszKdqqYAKQBJSUma1eUaYwqHLbu28K9v/8V/pv2HKmWrMLjFYFrUa4GIhB2aCUEYSbKMquakOrQZMFNVV/v3q0Wkmi9FVgPW5DxEY0xhNHzecNqPbM+yzctok9SGHk16UL5U+bDDMiEKo+HOCBHJSUeGt3OgqhVgONDKv26FK6kaY0zclm9ezk2Db+K6T66jfKnyTL5nMm9c/YYlSIOoBlvzKCJbgLLAbmCPH6yqengc05YBlgLHqeomP6wiMBioBSwBblbV9RnNJykpSadPn579L2GMKRD2pe6j3/R+PDH+Cfak7qHbxd3odF4nihctHnZoCUlEZqhqUthxBCmM1q3ZbjOtqtuBimmGrcO1djXGmLj9uupXWo9ozbTl02h6XFP6Xd2P4yscH3ZYJsGE0o5ZRK4FLvJvJ6jqiDDiMMYUPtt2b+OZCc/Q+8feVChdgY9u+Ih/nPYPa5hj0hXGLSA9gbOAAX5QBxG5QFWtEwBjTJ4a+edI2n3VjsWbFnPfmffxYtMXqVC6QthhmQQWRknyKuAMVU0FEJH+wM9YTznGmDyyausqOozuwOA5g6lTqQ4T757IRcdclPmEptALq9uII4BI4xprPmaMyROpmkrKjBS6juvKjr07eLbxs3Rp1IWSxUqGHZrJJwJLkiIyRlUvB3oAP4vIt4Dgrk1mqWMAY4zJzOw1s0kekcyUpVO4pPYlvNn8TU6qeFLYYZl8JsiSZGUAVR3onyl5Fi5JdlHVVQHGYYwpwHbs2cFz3z3Hy1NepnzJ8rx33Xu0qt/KGuaYbAkySZYXkRvTGX6+iKCqQwOMxRhTAI37axxtRrRhwYYF3FX/Ll65/BUqlakUdlgmHws0SQLNcaXHtBSwJGmMyZY129bQaUwnPpr1ESdWOJHxd43n0mMvDTssUwAEmSQXZ/FJH8YYkyFV5b8//5fOYzuzdfdWnr7oaZ648AlKFSsVdmimgAgySdoFAWNMrpn791ySRyTz3eLvuKDWBbzV/C3qVa4XdlimgAkySd4Z4LKMMQXUrr276DGpBz0m9aBM8TKkNE/h3gb3UkTCeF6DKegCS5KqOjuoZRljCqYJiyaQPCKZP9b9we2n3k7vK3pT9bCqYYdlCrCwOhMwxpi4rdu+js5jO/PeL+9x7BHHMrrlaK444YqwwzKFQBh9t5YFdkR1S1cEKOWf8GGMMfupKh/N+ohHxjzCxp0b6dqoK09f/DRlipcJOzRTSIRRkhwPXAZs9e/LAGOA80OIxRiToP5c9ydtv2rL+IXjOafGOaRck8LpVU8POyxTyISRJEupaiRBoqpb/cOUjTGG3ft28/Lkl3nuu+coWawkb1z1BslJydYwx4QijCS5TUQaqOpMABFpCOwIIQ5jTIKZtGQSySOS+W3tb7So14I+V/ahernqYYdlCrEwkmRH4FMRWeHfVwNuDSEOY0yC2LBjA13HdSVlZgq1ytdixO0juPqkq8MOy5jgk6Sq/iQidYCTcR0MzFXVPUHHYYwJn6oyaM4gOo7uyNrta+l0XieebfwsZUuUDTs0Y4BgH5V1qap+k04n5ydaB+fGFD4LNyyk3ch2jJ4/mqTqSYxqOYozq50ZdljGHCTIkuTFwDfANemMi6uDcxE5AngHONVPcw8wDxgE1AYWAbeo6obcCNgYk/v27NtD7x9788yEZyhapCivXvEq7c9uT9EiRcMOzZhDBNnjTjd/T+QoVR2czdn0AUaragsRKYG7feQJYLyq9hSRrkBXoEvuRG2MyU1Tl02l9YjWzFo9i+vrXM9rzV7j6MOPDjssY2IKtE2170CgfXamFZHDgYuAd/28dqvqRuA6oL//WH/g+hwHaozJVZt2buLBkQ9y3rvnsW77OobdOoxhtw6zBGkSXhitW8eKyKO4KtJtkYGquj6T6Y4D1gLviUh9YAbQAaiqqiv9PFaKSJX0JhaR1kBrgFq1auX4SxhjMqeqDP19KA+NfoiVW1by4NkP8vylz1OuZLmwQzMmLqKqwS5QZGE6g1VVj8tkuiTgR6CRqk4VkT7AZuBBVT0i6nMbVPXIjOaVlJSk06dPz3rwxpi4Ldm0hPYj2/PlH19yxlFnkNI8hbNqnBV2WCYHRGSGqiaFHUeQwihJ1lXVndEDRCSeJ6QuA5ap6lT/fgju+uNqEanmS5HVgDW5G64xJiv2pu7ltamv8fS3T6MovZr2osO5HShWxJ6nYPKfMPp5mhLnsIOo6ipgqYic7Ac1AX4DhgOt/LBWwBe5EaQxJutmrJjBOe+cwyNjHuHi2hczp90cOp3fyRKkybeCvE/yKKAGUFpEzsR1JABwOK6VajweBAb4lq1/Af/EJfrBInIvsAS4OVcDN8ZkauvurTz9zdP8Z9p/qFK2CoNbDKZFvRaISOYTG5PAgjy9uwK4Gzga+HfU8M242zgypaq/AOnVhzfJYWzGmGwaPm847Ue2Z9nmZbRJasMLTV7giFJHhB2WMbkiyPsk+wP9ReQmVf0sqOUaY/LG8s3LeWj0Qwz9fSinVjmVQS0GcV7N88IOy5hcFcaFgski8i5QXVWbiUg94DxVfTeEWIwxWbQvdR/9pvfjifFPsCd1Dy9c+gKPnv8oxYsWDzs0Y3JdGA133gO+BiLPv/kD92QQY0yC+3XVr5z/3/N5cNSDnHv0ucxuO5vHL3zcEqQpsMJIkpV8t3SpAKq6F9gXQhzGmDht272Nx8Y+RsOUhizauIgBNw7g6zu+5vgKx4cdmjF5KqyHLlfEdVCOiJwLbAohDmNMHEb9OYp2I9uxaOMi7jvzPl5s+iIVSlcIOyxjAhFGknwEd2/j8SIyGagMtAghDmNMBlZtXUXH0R0ZNGcQdSrVYeLdE7nomIvCDsuYQIXx0OWZInIxBx66PM8eumxM4kjVVN6e8TZdxnVhx94dPNv4Wbo06kLJYiXDDs2YwAXZmUDahy1HnGQPXTYmMcxZM4fWI1ozZekULql9CW82f5OTKp4UdljGhCbIkuQQ4Bf/Bwd63IE4H7psjMkbO/bs4PnvnuelKS9RvmR53r/ufe6qf5f1mGMKvSCT5E3ArcDpuP5VB6rq/ACXb4xJx7i/xtH2q7bMXz+fVvVb0evyXlQqUynssIxJCIHdAqKqw1T1NuBiYAHwiohM8tcnjTEBW7ttLXcOu5OmHzZFEMbdOY73r3/fEqQxUcJo3boTd8vHZqAWEM9jsowxuURVef+X93l07KNs2bWFpy96micufIJSxWxXNCatIBvuXALcDpwNjAP6qKo9+diYAM37ex7JI5KZuHgiF9S6gLeav0W9yvXCDsuYhBVkSXI8MAuYBJQE7hKRuyIjVfWhAGMxplDZtXcXPSf15IVJL1CmeBlSmqdwb4N7KSJhdLplTP4RZJL8Z4DLMsZ4ExdNJHlEMvPWzeP2U2+n9xW9qXpY1bDDMiZfCPpRWcaYgKzbvo7OYzvz3i/vcewRxzK65WiuOOGKsMMyJl8Jo+GOMSYPqSoD/jeAh79+mI07N9KlURf+dfG/KFO8TNihGZPvWJI0pgCZv34+bb9qy7i/xnHu0eeS0jyF06qeFnZYxuRbliSNKQB279tNrym9eO675yhRtARvXPUGyUnJ1jDHmBwKPEmKSGXgfqB29PJV9Z6gYzGmIJi8ZDLJI5KZs3YON9e7mVevfJXq5apnPqExJlNhlCS/AL7H3SuZpYcti8giYIufbq+qJolIBWAQLukuAm5R1Q25GK8xCWnDjg10HdeVlJkp1Cpfiy9v/5LmJzUPOyxjCpQwkmQZVe2Sg+kvUdW/o953Bcarak8R6erf52T+xiQ0VWXwnMF0GN2BtdvX8si5j/DsJc9yWInDwg7NmAInjAsWI0Tkqlyc33VA5PaS/sD1uThvYxLKoo2LuPrjq7nts9uoWb4m0++fzitXvGIJ0pg8EkaS7IBLlDtEZLOIbBGRzXFOq8AYEZkhIq39sKqquhLA/6+S3oQi0lpEpovI9LVr1+b4SxgTpD379vDy5Jep93o9vl/yPb2v6M2P9/7ImdXODDs0Ywq0wKtbVbVcDiZvpKorRKQKMFZE5mZhuSlACkBSUpLmIAZjAjV12VRaj2jNrNWzuO7k63it2WvULF8z7LCMKRSC7OC8jqrOFZEG6Y1X1ZmZzUNVV/j/a0RkGK6z9NUiUk1VV4pINWBNrgZuTEg279rMk+Of5PWfXqd6ueoMvWUoN9S9IeywjClUgixJPgK0Bl5JZ5wCl2Y0sYiUBYqo6hb/+nLg/4DhQCugp///RW4GbUzQVJVhc4fx4KgHWbllJQ+c9QDdm3Tn8JKHhx2aMYVOkH23tvb/L8nmLKoCw0QEXNwfq+poEfkJGCwi9wJLgJtzI15jwrBk0xLaj2zPl398Sf2q9Rl26zDOrnF22GEZU2jlmx53VPUvoH46w9cBTYKPyJjcszd1L32n9eWpb55CUXo17UWHcztQrEi+2UWNKZBsDzQmZDNXzuT+L+9n5sqZXHXiVbx+1evUPqJ22GEZY7AkaUxotu7eyr++/Rd9pvahStkqDG4xmBb1WuAvKRhjEkAoSVJEagDHcHDfrd+FEYsxYfhy3pc8MPIBlm5eSpuGbehxWQ+OKHVE2GEZY9IIo4PzF4Fbgd840HerApYkTYG3fPNyHhr9EEN/H8oplU9h8j2TOb/m+WGHZYyJIYyS5PXAyaq6K4RlGxOKfan7eHP6mzw+/nH2pO7hhUtfoNP5nShRtETYoRljMhBGkvwLKA5YkjSFwq+rfiV5RDJTl0+l6XFN6Xd1P46vcHzYYRlj4hBGktwO/CIi44lKlKr6UAixGJNntu/ZzrMTnuWVH16hQukKfHjDh7Q8raU1zDEmHwkjSQ73f8YUWKPnj6bdV+1YuHEh95xxDy9f/jIVSlcIOyxjTBaF0cF5fxEpAZzkB81T1T1Bx2FMXli1dRUPf/0wn8z+hJMrnszEuydy0TEXhR2WMSabwmjd2hj33MdFgAA1RaSV3QJi8rNUTeWdme/QZVwXV83a+Fm6NOpCyWIlww7NGJMDYVS3vgJcrqrzAETkJGAg0DCEWIzJsTlr5pA8IpnJSyfTuHZj3rz6TU6udHLYYRljckEYSbJ4JEECqOofIlI8hDiMyZEde3bQ/fvuvDT5JcqVLMd/r/0vd59xtzXMMaYACSNJTheRd4EP/fuWwIwQ4jAm28b/NZ42X7Vh/vr53FX/Lno17UXlspXDDssYk8vCSJJtgQeAh3DXJL8D3gghDmOybO22tXQa04kPZ33ICRVOYNyd42hynD2ExpiCKozWrbuAf/s/Y/IFVeX9X97n0bGPsmXXFp668CmevOhJShUrFXZoxpg8FFiSFJHBqnqLiPwP11frQVT19KBiMSYr5v09j+QRyUxcPJFGNRuRck0K9SrXCzssY0wAgixJdvD/mwe4TGOybdfeXfSc1JMXJr1AmeJleKv5W9zX4D6KSJGwQzPGBCSwJKmqK/3LdqraJXqcfzJIl0OnMiYcExdNJHlEMvPWzeO2U2+j9xW9Oeqwo8IOyxgTsDBOiZumM6xZ4FEYk471O9Zz7xf30rh/Y3bv282olqMYeNNAS5DGFFKBJUkRaeuvR9YRkVlRfwuB/wUVh0kcA4DauI2wtn8fFlXlo1kfUadvHfr/2p8ujbowu91srjzhyhCjMsaELchrkh8Do4AeQNeo4VtUdX28MxGRosB0YLmqNheRCsAg3HF2EXCLqm7IraBN3hgAtMY9EgZgMXAH7sJ1H9zNs0FZsH4Bbb9qy9i/xnJOjXMYd9c4Tq9q7ciMMQGWJFV1k6ouwh0D16vqYlVdDOwRkXOyMKsOwO9R77sC41X1RGA8Bydgk6Ce5ECCjLYOlzyDKFXu3rebHt/34NR+pzJ1+VT6NuvL5HsmW4I0xuwXxjXJfsDWqPfb/LBMicjRwNXAO1GDr8N1mI7/f33OQzR5bUkG47bjkmhemrJ0Cg1TGvLEN09w9YlX8/sDv/PA2Q9QtEjRPF6yMSY/CSNJiqruv09SVVOJv9r3VeAxIDVqWNVIy1n/v0q6CxVpLSLTRWT62rVrsxW4yT21MhmfURLNiY07N9JmRBsa/bcRm3dtZvhtwxlyyxCql6ueR0s0xuRnYSTJv0TkIREp7v86AH9lNpGINAfWqGq2+nlV1RRVTVLVpMqVrY/NMA3g4KqE9GSWRLNKVRk0exB1+tbh7Zlv88i5jzCn3RyuOfmaXF6SMaYgCaPv1jbAf4CncD3vjMddhspMI+BaEbkKKAUcLiIfAatFpJqqrhSRasCaPIrb5IK0DXbSUwbonovLXLhhIQ+MfIBR80fRsFpDRrYcSYNqDXJxCcaYgiqMvlvXALdlY7rHgcdh/4ObH1XVO0TkZaAV0NP//yLXgjW5LlaDnSK4OvSjcT9kbrRu3bNvD32m9qHbhG4IQu8retP+7PYUKxLGuaExJj8K/GghIu+Rft+t92Rzlj2BwSJyL+5S1s05CM/koQG4Wz3SE7nI/D3uXp6cmrZ8Gq2/bM2vq3/l2pOvpW+zvtQsXzMX5myMKUzCOKUeEfW6FHADsCIrM1DVCcAE/3odYM8qSnCRatZYKgF/A5tzuJzNuzbz5Pgnef2n16lWrhqf3fIZN9S5wR6EbIzJljCqWz+Lfi8iA4FxQcdhgjMAVw++L8b4MsA9wEvAlmwuQ1UZNncYD456kJVbVvLAWQ/QvUl3Di95eDbnaIwx4ZQk0zqR3G/MaBLAAFzPD+sy+VwKcCzZT5JLNy2l/aj2DJ83nPpV6zP0lqGcc3RW+qcwxpj0hXFNcgvumqT4/6uwJ4AUOPG0YgU4BtdIZ7Z/n5Xq1n2p++g7rS9PffsUqZrKS5e9xMPnPWwNc4wxuSaM6tZyQS/TBCuz6tWI6Fs9IhtFvCXJn1f+zP1f3s+MlTNodkIz3rj6DWofUTvrwRpjTAYCS5IikuGNaao6M6hYTN6JlCAzS5BFcdWskVs94k2SW9evotvQh3h1xWdULlOZT276hFtOucUa5hhj8kSQJclX/P9SQBLwK67K9XRgKnBBgLGYXDQAd//jEtz9jvGUIKMTJBxIkhlWt375JQ8PvJ13Tt5G8ql30/Pq3hxR6ojsBW2MMXEI8ikgl6jqJbhb5Rr4LuIaAmcC84OKw+SuSMlxMe4Cc2YJsiKHJkiA4rizp3RLkitWQIsWcO21PP1ndSYn9ePNm96zBGmMyXNhtHCoo6r7H7KsqrNF5IwQ4jC5IFYPOmkVxT2iJaOedMqRJknu2wdvvQWPPw67d8Pzz1Orc2dqlSiR/YCNMSYLwkiSv4vIO8BHuMLHHRz8fEiT4KKrVw/pOikd6VWvpudwoqpbZ82C1q1h6lS47DLo1w9OOCG7IRtjTLaE8RSQfwJzcLfQdQR+88NMPpC2ejUjgrvFI54ECa4kuXv7dujaFRo2hAUL4MMPYcwYS5DGmFCEcQvIThF5HdfLjgLzVHVP0HGY7Im3ejXe0mO0JqNH80i7drBwIdxzD7z0ElSsmK04jTEmNwRekvRP8PgT6Au8AfwhIhcFHYfJnngfhrwdV+IcEM+HV62C22+nV7Nm7CpRAiZMgHfftQRpjAldGNWtrwCXq+rFqnoRcAXQO4Q4TCYG4J7IUcT/H0DW+g/cjit5xpSaCm+/DXXrwtChDHnmGZr/+itcfHE2IzbGmNwVRpIsrqrzIm9U9Q/cHQAmQQzAPZXjDg5ce1yMKxlGnngdr5glz99+c8mwdWs44wyYNYsx3bqxrmTJ7AdujDG5LIwkOUNE3hWRxv7vbWBGCHGYdEQa5qTXKfl2YCTQLgvzO6TkuXMnPP20S4y//QbvvQfffAMnn3xw61ZjjEkAYdwC0gZ4AHgI1wDyO9y1SROyePpcXQLU8K9LAzsy+Gx036wAjB8PbdrA/Plw553wyitQufL+0eX8/PaSGI+nMcaYQEuSIlIEmKGq/1bVG1X1BlXtraq7gozDHCrePleL4EqTRwJv427xiNzq0TbN+/2tW9euhVat3P2OAGPHwgcfHJQg4UDXdFtz/G2MMSZ3BHrCrqqpIvKriNRS1XgbSpoAxHtrxz7gW1xDnpZkcouHKvTvD48+Cps3w5NPur/SpdP9eOTxyJuBI+KM2xhj8lIYtVrVgDkiMg3YFhmoqteGEIvxsnLGkop7CGiG5s1zVasTJkCjRq57uVNOyXCSrD4uyxhj8loYSfLZEJZpMjCA+J7eES1mqXPXLujZE154AcqUgZQUuPdeKJJ5zb4lSWNMognyeZKlcI12TgD+B7yrqnuzOP13QElc3ENUtZuIVAAG4WoAFwG3qOqG3I2+4MroWmQZXOOc9Fq6VkpvZt99B8nJMHcu3HYb9O4NRx0VdyzR1a3GGJMIgmy40x/3HMn/Ac048HzJeO0CLlXV+sAZwJUici7QFRivqicC4/17E6dY1yIjD0Xug0uWaT0U/Wb9erjvPnff465dMGoUDByYpQQJVpI0xiSeIKtb66nqaQAi8i4wLSsTq6pyoOFjcf+nwHVAYz+8PzAB6JLzcAuHWNciUzm4Uc6TuA4FIlKA41RpOXAgPPwwrFsHjz0G3bq5atZssCRpjEk0QZYk93dinpVq1mgiUlREfgHWAGNVdSpQVVVX+vmuBKrEmLa1iEwXkelr167NzuILpFjdzEUPb4m73zG6L5wSCxZw1BVXQMuWULs2zJgBL76Y7QQJVt1qjEk8QSbJ+iKy2f9tAU6PvBaRuI6LqrpPVc8AjgbOFpFT4124qqaoapKqJlVOc39eYTWA9Etth3QCgCtJ7gKK795N1x49mH3qqZz144883bcvTJkC9evnOB4rSRpjEk1g1a2qWjQX57VRRCYAVwKrRaSaqq4UkWq4UqbJRKTBTtrrkRVx1yHT3v+4BDhvyhTeSk7mtNmzGXLTTXTo04eVNWrwXC7FVBJXh25J0hiTKMLouzVbRKSyiBzhX5cGLgPmAsNxvanh/38RSoD5TKwGO4eRTgcBGzfyQdu2TGnUiPKbNnHN8OHcPGQIK2rUyNJTQeJh/bcaYxJJvkmSuE4IvhWRWcBPuGuSI4CeQFMR+RNo6t+bDAzg4EY40Q5qyKMKgwdD3br8IyWF1x5+mHq//caIa64B0q+WzalyHFqSTO+RXcYYE4R804+0qs4Czkxn+DqgSfAR5U+RatZY9pcMFy+Gdu1g5Eho0IAiI0ZQoWFDKuFKoLVwCTLDbumyIW2STFstHHlkF3mwbGOMSSs/lSRNDkWe8hGrt5wywAt790KvXlCvHkyc6DoEmDoVGjakJa63hlT/P7eT1ADgD1x9eW3/Pr1q4Uwf5myMMbkkyB53tuDuazxkFO42yMPTGWdySTxP+fhs2jSubN0afv0VrrkG+vaFWrl91THj+CKPg4mUGGMldOsd3xgThCBbt5bL/FMmr2T0lI9ymzfT56mnuLJvX6hWDYYMgRtvBJFQ49uO6/knvcQeTOo2xhR2oV2TFJEqQKnIe3t0Vt6KtXKv+/xzXm/fnuorVrhrkN27Q/nygcYGsePbh6sGjk6gedFgyBhj0hP4NUkRuda3RF0ITMRd3hoVdBwFUaQVqODOfqL/p/2hj166lGHXX8/nN9xAmYoVkR9+cNWrISRIiF0yjDy8OVKmLULUw5yNMSaPhdFw5zngXOAPVT0W1zJ1cghxFCiRa3qRWzv2xfhfZN8+HurTh9/q1ePyMWOY+dJLHDl9OpxzTpDhHqI7h3akHikx3oi7mF0L12go3EiNMYVJGElyj79to4iIFFHVb3FP9TA5kNE1x4gzZ85k6jnn0KdjR6ZecAFj5syhQefOULx4ECFmqCWuhBh5BFc1DpQYV/thkT4IT8TulzTGBCOMa5IbReQw4HtggIisAbLV4bk5IKMLumW3buXZbt3o+OqrrK1cmdsGDuTHW29lUYANc+LREjgK15XSxxx4tMtK/3981GftfkljTBDCKEleB+wAOgKjgQXANSHEUaDEuqZ31VdfMeeUU+j073/zzn33Uff33xl0220sSbAEGVHV/18dNWyV/78rzWftfkljTF4LvCSpqttEpCpwFu6h96N89avJpgHApjTDqq1YQZ8OHbh5yBDm1KtHo0mTmNKo0f7xiXoLReQ5Z9G91K9K74OeNYk2xuSlMFq33oJ74PLNwC3AVBFpEXQcBUWkwc5G/15SU2n7xhv8Xrcu13z5JU92786ZP/98UIJM5FsoKuI2yugkuTLGZyFxk70xpmAI45rkk8BZqroG3NM9gHHAkBBiyfeiG+yc+r//kdK6Nef9+COTmjThgjffpPsJJ1DPf24Jedfnam4pimu8k7a69XDchWu7X9IYE6QwrkkWiSRIb11IcRQIS4DS27fTo2tXZjZowAnz53PnBx9w0dixcMIJAHne52puq8qh1a3H4lq7HoO7ZzJy/2SifxdjTP4WRklytIh8DQz072/FOhPItju+/ppn2rbluIUL+e8//0nnl19mfcWKHBN2YDlQhUNLktVwCdGSojEmSGE03OksIjcCF+AKBSmqOizoOPK91atZ+MgjfPDxx8w76SQaf/stExs3BvJ/NWRV4Meo9yuBU0KKxRhTuAVWzSkiJ4hIIwBVHaqqj6jqw8A6ETk+qDjyvdRUeOcddtWpQ/UhQ+j2zDOcPmvW/gRZkfxfDRldkkz1r48KLxxjTCEW5LXAVzn0ofPg2mK8GmAc+dfvv8PFF8P99/Nz/frU//VX/q9bN3aXLLn/I4eRvxMkuJLkNv+3AdiDq241xpigBZkka6vqrLQDVXU6rpcxE8vOnfCvf0H9+vDbb/Duu5z/7bfMq1PnkI8WhPsGI/dKruXA7R9WkjTGhCHIJFkqg3GlA4siH4g8zaMIcPs337D59NPhuefg1lth7lwG3HMPRWL0mFMQ7huM7nUn0pGAJUljTBiCTJI/icj9aQeKyL3AjADjSGiRzgG2/v0377VqxcAmTVibmsr4sWPhww8ZULkyrUn/QcT5vcFORKQkORi4zb++HevQ3BgTvCBbt3YEholISw4kxSSgBHBDZhOLSE3gA1yhIhXXKraPiFQABuEKX4uAW1R1Q24HH4QBQCtVWn7wAa906kT5TZvo/sQTPP/UU1QtXZpFxH7aR1Hyf4OdiEhJsi+w279egXVobowJnqhqsAsUuYQDTz2ao6rfxDldNaCaqs4UkXK4RHs9cDewXlV7ikhX4EhV7ZLRvJKSknT69OnZ/Qp5YgDw0h9/0LtNGy799lsmn38+yW+9xZxT3aoS3JlBEdyzFdOKjC8IdhK7/v0Y3JmQMSZ4IjJDVZPCjiNIYdwn+S3wbTamW4lvx6GqW0Tkd6AG7qkijf3H+gMTgAyTZMLZtYvVL77ItO7d2VG6NMlvvsnb99+PFjlQG14EaOf/p1fVWhCuRUZkdPG6IDRMMsbkH2H0uJNjIlIbOBOYClT1CRRVXSkiVWJM0xpfY1erVgKllO+/Z1Pr1jwydy6f3HorHV99ldVHHdpMZR/QL8YsCsq1yIiMrj0m0C9njCkE8l2fqf6BzZ8BHVV1c7zTqWqKqiapalLlypXzLsB4rV/P/Pvvh4suYv3OnTQbOZLbP/kk3QSZkYJ0LRIONFxKT0E7GTDGJL58lSRFpDguQQ5Q1aF+8Gp/vTJy3XJNrOkTgip8/DE76tal9nvv8eJjj3Hq7NmMbtYsW7NLpeAkSCgcDZOMMflHvkmSIiLAu8DvqvrvqFHDgVb+dSvgi6Bji9uCBXDlldCyJfNq16bhjBl0ffFFtpctm+1ZFrTqx1jXHAvayYAxJn/IN0kSaATcCVwqIr/4v6uAnkBTEfkTaOrfJ5Y9e6BHDzj1VPjhB+jbl6QpU5hVv37MSSriqhczUhCrH2Ml/YJ2MmCMyR/yTcMdVZ2Eu9MhPU2CjCVLfvgBWreG2bPhxhvhP/+BGjU4GlgcY5IyQB//ugPugZtpVfSfKWilq+64a5L2cGVjTCLITyXJhBXdjVxt/56NG6FdO2jUCDZtguHD4bPPoEYNwF17S0/0UzxaAn8DH3Hww4Y/8sMLWoIE953s4crGmERhSTKHIq0xF+Nu8l+syohPP2V13brse+st+nToQLk5c6h9zTX7b20YADzuX0d+gIySX0vcDfSp/n9BTxiF7fsaYxJXvqluTVTRrTGPWbSIvu3b0/yrr5jRoAFXjRjBzIYNAdgK3OH/oqVyoDrRkoExxiQWK0nm0BKg6N69dOrVizmnnELjCRPo2Ls350yduj9BZmY7satfjTHGhMdKkjl01U8/8Vzr1pz5yy982bw5D7z+Okuz0aOPdbdmjDGJx5JkNg3avJntTz3F8L59WXXUUdw0ZAhDb7wRYjznMTN2i4MxxiQeS5LZMPHzz7mgfXuqrVjBG+3a8WT37mwuXz7b87NbHIwxJjHZNcmsWLYMbriBi2+4gb8rVuT8KVN4sG/f/QkyO2XI6Fs+jDHGJBZLklkxdy58/TVdXnyRpOnTmXruuQeNVg69n1E5cJ8juD5IoeDf72iMMQWBVbdmxWWXweLFDKpcmb3pjI71QOBIxwDGGGPyFytJZlXlynQHSqcZbNcVjTGm4LEkmQ0tOfDYEes6zRhjCi6rbo3TANwN/4tx1xX34c4w+nNoLzrGGGMKBkuScYj0zxrpfm6f/58KJONKk1aKNMaYgseqW+MQ3T9rWtalnDHGFFyWJOOQWZdx1qWcMcYUTJYk45BZl3HWpZwxxhRMliTj0B13i0d67NYPY4wpuCxJxqEl7haP9HrNsVs/jDGm4LLWrXGyXnOMMabwyVclSRH5r4isEZHZUcMqiMhYEfnT/z8yzBiNMcYUHPkqSQLvA1emGdYVGK+qJwLj/XtjjDEmx/JVklTV74D1aQZfh+v4Bv//+iBjMsYYU3DlqyQZQ1VVXQng/1dJ70Mi0lpEpovI9LVr1wYaoDHGmPypICTJuKhqiqomqWpS5cqVww7HGGNMPlAQWreuFpFqqrpSRKoBazKbYMaMGX+LyOIAYotXJdzzlxOVxZdziR6jxZcziR4f5E6Mx2T+kYKlICTJ4bgnV/X0/7/IbAJVTaiipIhMV9WksOOIxeLLuUSP0eLLmUSPD/JHjIkoX1W3ishA4AfgZBFZJiL34pJjUxH5E2jq3xtjjDE5lq9Kkqp6e4xRTQINxBhjTKGQr0qSBVhK2AFkwuLLuUSP0eLLmUSPD/JHjAlHVDXsGIwxxpiEZCVJY4wxJgZLksYYY0wMliTzQIyO2M8QkR9F5Bff88/ZfnhxEekvIv8Tkd9F5PGoaRr64fNF5D8iInkYX30R+cEv70sROTxq3OM+hnkickVex5fVGEWkqYjM8MNniMileR1jVtehH19LRLaKyKOJFp+InO7HzfHjS+VlfFmNMej9RERqisi3fllzRKSDHx7zgQpB7ydZjTGM/aRAUFX7y+U/4CKgATA7atgYoJl/fRUwwb/+B/CJf10GWATU9u+nAecBAoyKTJ9H8f0EXOxf3wM851/XA34FSgLHAguAonkZXzZiPBOo7l+fCiyPmib0dRg1/jPgU+DRRIoP18p9FlDfv6+YgL9xoPsJUA1o4F+XA/7w+8JLQFc/vCvwYlj7STZiDHw/KQh/VpLMA5p+R+wKRM7cywMrooaXFZFiQGlgN7BZXO9Bh6vqD+q24g/Ipc7bY8R3MvCdfz0WuMm/vg53cNqlqguB+cDZeRlfVmNU1Z9VNbI+5wClRKRkAq1DROR64C8fX2RYosR3OTBLVX/1065T1X2J9BsT8H6iqitVdaZ/vQX4HahB7AcqBL6fZDXGMPaTgsCSZHA6Ai+LyFKgFxCpLhoCbANWAkuAXqq6HrexL4uafpkflldmA9f61zcDNf3rGsDSdOIIOr6MYox2E/Czqu4iQdahiJQFugDPpvl8QsQHnASoiHwtIjNF5LGQ4ssoxtD2ExGpjSuFTSX2AxVC3U/ijDFamPtJvmJJMjhtgYdVtSbwMPCuH342sA+ojqum6SQix+GqPdLKy/t17gEeEJEZuKqb3X54rDiCjg9ixwiAiJwCvAgkRwalM48w1uGzQG9V3Zrm84kSXzHgAqCl/3+DiDQJIb6MYgxlPxGRw3DV5B1VdXNGH40RR56vwyzEGPl82PtJvpKvetzJ51oBHfzrT4F3/Ot/AKNVdQ+wRkQmA0nA98DRUdMfzYEq2lynqnNx1W6IyEnA1X7UMg4usUXiWBZkfJnEiIgcDQwD7lLVBX5woDFmEN85QAsReQk4AkgVkZ24A1sixLcMmKiqf/txI3HXCj8KMr5MYgx8PxGR4rjfaICqDvWDYz1QIZT9JIsxJsR+kt9YSTI4K4CL/etLgT/96yXApeKUBc4F5vpqki0icq5vaXYXcXTenl0iUsX/LwI8BbzpRw0HbvPXLo4FTgSmBR1fRjGKyBHAV8Djqjo58vlEWYeqeqGq1lbV2sCrwAuq2jdR4gO+Bk4XkTL+mt/FwG+J9BsT8H7i5/Uu8Luq/jtqVOSBCnDwAxUC30+yGmOi7Cf5TtgthwriHzAQd+1kD+4s7V5cNdYMXAu4qUBD/9nDcCXLOcBvQOeo+SThrtEsAPrie0jKo/g64FrH/YHrJF6iPv+kj2EeUa3e8iq+rMaIO5huA36J+quSSOswarpnOLh1a0LEB9zht8HZwEsJ+BsHup/g9lfFtfqNbFNX4Vr+jsed5I4HKoS1n2Q1xjD2k4LwZ93SGWOMMTFYdasxxhgTgyVJY4wxJgZLksYYY0wMliSNMcaYGCxJGmOMMTFYkjQmAP7+vkki0ixq2C0iMjrMuIwxGbNbQIwJiIicirvX70ygKO4+tSv1QM8nWZlXUVXdl7sRGmPSsiRpTIB813TbgLL+/zHAabguIp9R1S98Z9Uf+s8AtFfVKSLSGOiGuwH/DFWtF2z0xhQ+liSNCZDvUm0mruPuEcAcVf3Idxk2DVfKVCBVVXeKyInAQFVN8knyK+BUdY9jMsbkMevg3JgAqeo2ERkEbAVuAa4RkUf96FJALVw/v31F5Azcky9OiprFNEuQxgTHkqQxwUv1fwLcpKrzokeKyDPAaqA+rnHdzqjR2wKK0RiDtW41JkxfAw/6Jy8gImf64eWBlaqaCtyJa+RjjAmBJUljwvMcUByYJSKz/XuAN4BWIvIjrqrVSo/GhMQa7hhjjDExWEnSGGOMicGSpDHGGBODJUljjDEmBkuSxhhjTAyWJI0xxpgYLEkaY4wxMViSNMYYY2L4f2ARyMZXlg7MAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.plot(X, y, color='cyan', marker='o')\n",
"plt.title(\"French Coal Production 1870-1919\")\n",
"plt.xlabel(\"Year\")\n",
"plt.ylabel(\"Coal Production in Metric Tons\")\n",
"\n",
"regr = linear_model.LinearRegression()\n",
"regr.fit(X, y)\n",
"print(regr.coef_, regr.intercept_)\n",
"\n",
"y_predict = regr.predict(X)\n",
"plt.plot(X, y_predict, color='red')\n",
"plt.title(\"French Coal Production Linear Regression\")\n",
"plt.xlabel(\"Year\")\n",
"plt.ylabel(\"Coal Production in Metric Tons\")\n",
"\n",
"X_future = np.array(range(1918, 2022))\n",
"X_future = X_future.reshape(-1, 1)\n",
"\n",
"future_predict = regr.predict(X_future)\n",
"plt.plot(X_future, future_predict, color='green')\n",
"plt.title(\"French Coal Production Linear Regression Predicted from 1870-1918 Data\")\n",
"plt.ylabel(\"Coal Production in Metric Tons\")\n",
"plt.xlabel(\"Year\")\n",
"plt.show()"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "257dd31c",
"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.8.3"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment