{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "tags": [
     "remove_cell"
    ]
   },
   "outputs": [],
   "source": [
    "from datascience import *\n",
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "from utils import *\n",
    "plt.style.use('seaborn-muted')\n",
    "from matplotlib import patches\n",
    "import csaps\n",
    "import warnings\n",
    "warnings.filterwarnings(\"ignore\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# An Empirical Example from EEP 147\n",
    "\n",
    "Let's take a look at an empirical example of production. The dataset for this section comes from EEP 147: Regulation of Energy and the Environment."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table border=\"1\" class=\"dataframe\">\n",
       "    <thead>\n",
       "        <tr>\n",
       "            <th>Group</th> <th>Group_num</th> <th>UNIT NAME</th> <th>Capacity_MW</th> <th>Average Variable Cost</th>\n",
       "        </tr>\n",
       "    </thead>\n",
       "    <tbody>\n",
       "        <tr>\n",
       "            <td>Old Timers  </td> <td>7        </td> <td>BIG CREEK      </td> <td>1000       </td> <td>0                    </td>\n",
       "        </tr>\n",
       "        <tr>\n",
       "            <td>Fossil Light</td> <td>8        </td> <td>HELMS          </td> <td>800        </td> <td>0.5                  </td>\n",
       "        </tr>\n",
       "        <tr>\n",
       "            <td>Fossil Light</td> <td>8        </td> <td>DIABLO CANYON 1</td> <td>1000       </td> <td>11.5                 </td>\n",
       "        </tr>\n",
       "        <tr>\n",
       "            <td>Bay Views   </td> <td>4        </td> <td>MOSS LANDING 6 </td> <td>750        </td> <td>32.56                </td>\n",
       "        </tr>\n",
       "        <tr>\n",
       "            <td>Bay Views   </td> <td>4        </td> <td>MOSS LANDING 7 </td> <td>750        </td> <td>32.56                </td>\n",
       "        </tr>\n",
       "        <tr>\n",
       "            <td>Old Timers  </td> <td>7        </td> <td>MOHAVE 1       </td> <td>750        </td> <td>34.5                 </td>\n",
       "        </tr>\n",
       "        <tr>\n",
       "            <td>Old Timers  </td> <td>7        </td> <td>MOHAVE 2       </td> <td>750        </td> <td>34.5                 </td>\n",
       "        </tr>\n",
       "        <tr>\n",
       "            <td>Big Coal    </td> <td>1        </td> <td>FOUR CORNERS   </td> <td>1900       </td> <td>36.5                 </td>\n",
       "        </tr>\n",
       "        <tr>\n",
       "            <td>Bay Views   </td> <td>4        </td> <td>MORRO BAY 3&4  </td> <td>665        </td> <td>36.61                </td>\n",
       "        </tr>\n",
       "        <tr>\n",
       "            <td>East Bay    </td> <td>6        </td> <td>PITTSBURGH 5&6 </td> <td>650        </td> <td>36.61                </td>\n",
       "        </tr>\n",
       "    </tbody>\n",
       "</table>\n",
       "<p>... (32 rows omitted)</p>"
      ],
      "text/plain": [
       "Group        | Group_num | UNIT NAME       | Capacity_MW | Average Variable Cost\n",
       "Old Timers   | 7         | BIG CREEK       | 1000        | 0\n",
       "Fossil Light | 8         | HELMS           | 800         | 0.5\n",
       "Fossil Light | 8         | DIABLO CANYON 1 | 1000        | 11.5\n",
       "Bay Views    | 4         | MOSS LANDING 6  | 750         | 32.56\n",
       "Bay Views    | 4         | MOSS LANDING 7  | 750         | 32.56\n",
       "Old Timers   | 7         | MOHAVE 1        | 750         | 34.5\n",
       "Old Timers   | 7         | MOHAVE 2        | 750         | 34.5\n",
       "Big Coal     | 1         | FOUR CORNERS    | 1900        | 36.5\n",
       "Bay Views    | 4         | MORRO BAY 3&4   | 665         | 36.61\n",
       "East Bay     | 6         | PITTSBURGH 5&6  | 650         | 36.61\n",
       "... (32 rows omitted)"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ESG_table = Table.read_table('ESGPorfolios_forcsv.csv').select(\n",
    "    \"Group\", \"Group_num\", \"UNIT NAME\", \"Capacity_MW\", \"Total_Var_Cost_USDperMWH\"\n",
    ").sort(\"Total_Var_Cost_USDperMWH\", descending = False).relabel(4, \"Average Variable Cost\")\n",
    "ESG_table"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This table shows some electricity generation plants in California and their costs. The `Capacity` is the output the firm is capable of producing. The `Average Variable Cost` shows the minimum variable cost per megawatt (MW) produced. At a price below AVC, the firm supplies nothing. At a price above the AVC, the firm can supply up to its capacity. Being a profit-maximising firm, it will try to supply its full capacity. \n",
    "\n",
    "First, lets look at just the Big Coal producers and understand this firm's particular behavior."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table border=\"1\" class=\"dataframe\">\n",
       "    <thead>\n",
       "        <tr>\n",
       "            <th>Group</th> <th>Group_num</th> <th>UNIT NAME</th> <th>Capacity_MW</th> <th>Average Variable Cost</th>\n",
       "        </tr>\n",
       "    </thead>\n",
       "    <tbody>\n",
       "        <tr>\n",
       "            <td>Big Coal</td> <td>1        </td> <td>FOUR CORNERS        </td> <td>1900       </td> <td>36.5                 </td>\n",
       "        </tr>\n",
       "        <tr>\n",
       "            <td>Big Coal</td> <td>1        </td> <td>HUNTINGTON BEACH 1&2</td> <td>300        </td> <td>40.5                 </td>\n",
       "        </tr>\n",
       "        <tr>\n",
       "            <td>Big Coal</td> <td>1        </td> <td>REDONDO 5&6         </td> <td>350        </td> <td>41.94                </td>\n",
       "        </tr>\n",
       "        <tr>\n",
       "            <td>Big Coal</td> <td>1        </td> <td>REDONDO 7&8         </td> <td>950        </td> <td>41.94                </td>\n",
       "        </tr>\n",
       "        <tr>\n",
       "            <td>Big Coal</td> <td>1        </td> <td>HUNTINGTON BEACH 5  </td> <td>150        </td> <td>66.5                 </td>\n",
       "        </tr>\n",
       "        <tr>\n",
       "            <td>Big Coal</td> <td>1        </td> <td>ALAMITOS 7          </td> <td>250        </td> <td>73.72                </td>\n",
       "        </tr>\n",
       "    </tbody>\n",
       "</table>"
      ],
      "text/plain": [
       "Group    | Group_num | UNIT NAME            | Capacity_MW | Average Variable Cost\n",
       "Big Coal | 1         | FOUR CORNERS         | 1900        | 36.5\n",
       "Big Coal | 1         | HUNTINGTON BEACH 1&2 | 300         | 40.5\n",
       "Big Coal | 1         | REDONDO 5&6          | 350         | 41.94\n",
       "Big Coal | 1         | REDONDO 7&8          | 950         | 41.94\n",
       "Big Coal | 1         | HUNTINGTON BEACH 5   | 150         | 66.5\n",
       "Big Coal | 1         | ALAMITOS 7           | 250         | 73.72"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "selection = 'Big Coal'\n",
    "Group = ESG_table.where(\"Group\", are.equal_to(selection))\n",
    "Group"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "tags": [
     "remove_input"
    ]
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAicAAAGECAYAAADtMIAAAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nO3de7gkdX3n8fcHEFBBbk4mhFFARJQYAR0NKqvhoiJRwYRVlDUTY5ZNolmNyWNgTVwSNMFcNJoYDVFkTFBQooGooASJrEbR4Q4CAgrIcJlRQdB4g3z3j6oDnfGcMz3DdJ9fn36/nqefrvpVdfX3N9XMfKj6VVWqCkmSpFZsttAFSJIkDTKcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFE0iaT5D1J/nCh61hXklOSvHmh65A0HMOJpKEluTHJ95N8N8mdST6R5FEzy6vqN6rqhAex/ecluSDJPUnWJvlskhdtmuolTQrDiaQN9cKq2gbYGbgD+OtNsdEkRwIfAT4ALAOWAm8CXrgpti9pchhOJG2UqvoBcAaw90zbuqdPkrwhyW1Jbk3y60kqyWPX3VaSAG8DTqiq91bVd6rqP6vqs1X1P/t1NkvyB0luSrImyQeSbDewjY8kuT3Jd/qjLz87yv5LGh3DiaSNkuRhwEuBL86x/FDg9cAhwGOBX5hnc3sBj6ILO3P51f51IPAYYBvgbwaWnw3sCfwUcDFw6no7IalJWyx0AZImzj8nuRd4OLAWeN4c670EeH9VXQWQ5Hjg6DnW3al/v22e7z0aeFtVfa3f3nHAlUleWVX3VtXJMyv233Vnku2q6jvDdUtSKzxyImlDHVFV2wNbA68BPpvkp2dZ72eAbwzMf2OWdWZ8q3/feZ51fga4aWD+Jrr/wVqaZPMkJya5IcndwI39Oo+cZ3uSGmU4kbRRquq+qvoocB9wwCyr3EY3sHXGo2ZZZ8a1dOHll+dZ51Zg14H5RwP30g3KfTlwON0ppO2A3fp1Ms/2JDXKcCJpo6RzOLADcPUsq3wYeGWSJ/TjU+a8/0lVFd34lD9M8sokj+gHwB6Q5KR+tQ8Bv5Nk9yTbAH8CnF5V9wLbAj+kOwLzsH6ZpAllOJG0of4lyXeBu4G3ACtmxpUMqqqzgXcC5wPX88DA2R/OttGqOoNugO2v0R0luQN4M3Bmv8rJwD8AFwBfB34A/Ha/7AN0p3lWA19hjkG6kiZDuv9hkaTRSvIE4Epgq/5ohyTNyiMnkkYmyYuTbJVkB+CtwL8YTCStj+FE0ij9L2ANcAPdwNnfXNhyJE0CT+tIkqSmeOREkiQ1xXAiSZKaMhG3r3/kIx9Zu+2220KXIUmSNoGLLrrom1W1ZK7lExFOdtttN1atWrXQZUiSpE0gyU3zLfe0jiRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSWrcLst2JclYX7ss23XB+rvFgn2zJEkayq2rb+bQYy8Z63eec+J+Y/2+QR45kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmjCycJNkryaUDr7uTvC7JjknOTXJd/77DqGqQJEmTZ2ThpKqurap9q2pf4CnAfwAfA44FzquqPYHz+nlJkiRgfKd1DgZuqKqbgMOBlX37SuCIMdUgSZImwLjCyVHAh/rppVV1Wz99O7B0TDVIkqQJMPJwkmRL4EXAR9ZdVlUF1ByfOybJqiSr1q5dO+IqJUlSK8Zx5OT5wMVVdUc/f0eSnQH69zWzfaiqTqqq5VW1fMmSJWMoU5IktWAc4eRlPHBKB+AsYEU/vQI4cww1SJKkCTHScJLk4cBzgI8ONJ8IPCfJdcAh/bwkSRIAW4xy41X1PWCnddq+RXf1jiRJ0k/wDrGSJKkphhNJktQUw4kkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJI7LLsl1J8qBf02akz9aRJGma3br6Zg499pIHvZ1zTtxvE1QzOTxyIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJktQUw4kkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1JSRhpMk2yc5I8k1Sa5O8vQkOyY5N8l1/fsOo6xBkiRNllEfOXkHcE5VPR7YB7gaOBY4r6r2BM7r5yVJkoARhpMk2wHPAt4HUFU/qqq7gMOBlf1qK4EjRlWDJEmaPKM8crI7sBZ4f5JLkrw3ycOBpVV1W7/O7cDS2T6c5Jgkq5KsWrt27QjLlCRJLRllONkCeDLw7qraD/ge65zCqaoCarYPV9VJVbW8qpYvWbJkhGVKkqSWjDKc3ALcUlUX9vNn0IWVO5LsDNC/rxlhDZIkacKMLJxU1e3AN5Ls1TcdDHwFOAtY0betAM4cVQ2SJGnybDHi7f82cGqSLYGvAa+kC0QfTvIq4CbgJSOuQZIkTZCRhpOquhRYPsuig0f5vZIkaXJ5h1hJktQUw4kkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJktQUw4kkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1JQtRrnxJDcC9wD3AfdW1fIkOwKnA7sBNwIvqao7R1mHJEmaHOM4cnJgVe1bVcv7+WOB86pqT+C8fl6SJAlYmNM6hwMr++mVwBELUIMkSWrUqMNJAZ9OclGSY/q2pVV1Wz99O7B0tg8mOSbJqiSr1q5dO+IyJUlSK0Y65gQ4oKpWJ/kp4Nwk1wwurKpKUrN9sKpOAk4CWL58+azrSJKkxWekR06qanX/vgb4GPA04I4kOwP072tGWYMkSZos6w0nSR6X5LwkV/bzT0ryB0N87uFJtp2ZBp4LXAmcBazoV1sBnLmxxUuSpMVnmCMnfw8cB/wYoKouB44a4nNLgc8luQz4EvCJqjoHOBF4TpLrgEP6eUmSJGC4MScPq6ovJRlsu3d9H6qqrwH7zNL+LeDgoSuUJElTZZgjJ99MsgfdlTckORK4bf6PSJIkbZxhjpy8mu6qmccnWQ18HfgfI61KkiRNrfWGk/70zCH9oNbNquqe0ZclSZKm1TBX6/xJku2r6ntVdU+SHZK8eRzFSZKk6TPMmJPnV9VdMzP9Q/oOG11JkiRpmg0TTjZPstXMTJKHAlvNs74kSdJGG2ZA7KnAeUne38+/kgce3CdJkrRJDTMg9q1JLueBe5OcUFWfGm1ZkiRpWg314L+qOhs4e8S1SJIkzR1Oknyuqg5Icg/9DdhmFtE9UPgRI69OkiRNnTnDSVUd0L9vO75yJEnStJv3ap0kmye5ZlzFSJIkzRtOquo+4Nokjx5TPZIkacoNMyB2B+CqJF8CvjfTWFUvGllVkiRpag0TTv5w5FVIkiT15g0nSY4AHgtc4b1NJEnSOMw55iTJ3wK/A+wEnJDEIyiSJGnk5jty8ixgn6q6L8nDgP8HnDCesiRJ0rSa72qdH/VX61BV/0F38zVJkqSRmu/IyeP7Z+pAF0z26Odn7hD7pJFXJ0mSps584eQJY6tCkiSpN184OQk4Bzi7qrxLrCRJGov5xpysAO4Ejk9ycZJ3Jzk8ycPHVJskSZpC8z3473bgFOCUJJsBPw88H3hDku8Dn66qPxtLlZIkaWqs9w6xSZ5ZVZ8HvtC/3pTkAGDXURcnSZKmz7wP/uv99Sxt76yqUzd1MZIkSXMeOUnydOAZwJIkrx9Y9Ahg81EXJkmSptN8p3W2BLbp19l2oP1u4MhRFiVJkqbXfANiPwt8NskpVXUTQD8wdpuquntcBUrSKO2ybFduXX3zQpcxry22fCj3/uj7C12GNDbrHRAL/GmS3wDuA74MPCLJO6rqz0dbmiSN3q2rb+bQYy9Z6DLmdc6J+zVfo2Z3zon7LXQJE2mYAbF790dKjgDOBnYHXjHSqiRJ0tQaJpw8JMlD6MLJWVX1Y6BGW5YkSZpWw4STvwNuBB4OXJBkV7pBsZIkSZvcesNJVb2zqnapqsOqcxNw4LBfkGTzJJck+Xg/v3uSC5Ncn+T0JFs+iPolSdIis95wkmS7JG9Lsqp//SXdUZRhvRa4emD+rcDbq+qxdM/uedUGVSxp4u2ybFeSNPGS1J5hrtY5GbgSeEk//wrg/cAvre+DSZYBvwi8BXh9ur8JDgJe3q+yEjgeePcGVS1porV0hYxXU0jtGSac7FFVvzww/0dJLh1y+38FvIEHbuK2E3BXVd3bz98C7DLbB5McAxwD8OhHP3rIr5MkSZNumAGx3+8f9Ad0DwIE1ns3oCQvANZU1UUbU1hVnVRVy6tq+ZIlSzZmE5IkaQINc+TkN4APJNmun78T+NUhPvdM4EVJDgO2pnsmzzuA7ZNs0R89WQas3uCqJUnSojXM1TqXVdU+wJOAJ1XVflV12RCfO66qllXVbsBRwGeq6mjgfB54Ns8K4MyNrl6SJC06c4aTJK9Pcv+VNFV1d1XdneRVSV73IL7z9+kGx15PNwblfQ9iW5IkaZGZ77TO0cD+s7T/A7CKbrDrUKrq34B/66e/Bjxt6AolSdJUme+0zhb9rer/i6r6EeDNASRJ0kjMF042S7J03cbZ2iRJkjaV+cLJnwOfSPLsJNv2r18APg78xViqkyRJU2fOMSdV9YEka4E/Bp5I9yTiq4A3VdXZY6pPkiRNmXnvc9KHEIOIJEkam2HuECtJkjQ2hhNJktQUw4kkSWrKesNJkqVJ3pfk7H5+78E7x0qSJG1Kwzz47xTg/cAb+/mvAqezSG47v8uyXbl19c0LXYYkSeoNE04eWVUfTnIcQFXdm+S+Edc1NreuvplDj71kocuQpso5J+630CVIatgwY06+l2QnuvuckGR/4DsjrUqSJE2tYY6cvB44C9gjyeeBJcCRI61KkiRNrfWGk6q6OMmzgb3oHvh37WwPBJQkSdoU5gwnSX5pjkWPS0JVfXRENUmSpCk235GTF86zrADDiSRJ2uTme/DfK8dZiCRJEgx3E7adkrwzycVJLkryjv7qHUmSpE1umEuJTwPWAr9Md5XOWrqbsEmSJG1yw1xKvHNVnTAw/+YkLx1VQZIkaboNc+Tk00mOSrJZ/3oJ8KlRFyZJkqbTfJcS30N3VU6A1wH/2C/aDPgu8Hsjr06SJE2d+a7W2XachUiSJMFwY05IsgOwJ7D1TFtVXTCqoiRJ0vRabzhJ8uvAa4FlwKXA/sAXgINGW5okSZpGwwyIfS3wVOCmqjoQ2A+4a6RVSZKkqTVMOPlBVf0AIMlWVXUN3UMAJUmSNrlhxpzckmR74J+Bc5PcCdw02rIkSdK0Wm84qaoX95PHJzkf2A44Z6RVSZKkqTXffU4eUVV3J9lxoPmK/n0b4NsjrUySJE2l+Y6cfBB4AXARD9yMbfD9MSOvTpIkTZ35bsL2giQBnl1VN4+xJkmSNMXmvVqnqgr4xJhqkSRJGupS4ouTPHVDN5xk6yRfSnJZkquS/FHfvnuSC5Ncn+T0JFtucNWSJGnRGiac/DzwhSQ3JLk8yRVJLh/icz8EDqqqfYB9gUOT7A+8FXh7VT0WuBN41cYWL0mSFp9h7nPyvI3ZcH9K6Lv97EP6V9Hd9v7lfftK4Hjg3RvzHZIkafFZ75GTqrqpqm4Cvk8XLmZe65Vk8ySXAmuAc4EbgLuq6t5+lVuAXeb47DFJViVZtXbt2mG+TpIkLQLrDSdJXpTkOuDrwGeBG4Gzh9l4Vd1XVfvSPTTwacDjhy2sqk6qquVVtXzJkiXDfkySJE24YcacnED3JOKvVtXuwMHAFzfkS6rqLuB84OnA9klmTictA1ZvyLYkSdLiNkw4+XFVfQvYLMlmVXU+sHx9H0qypH8mD0keCjwHuJoupBzZr7YCOHOjKpckSYvSMANi70qyDXABcGqSNcD3hvjczsDKJJvThaAPV9XHk3wFOC3Jm4FLgPdtZO2SJGkRGiacHA78APgd4Gi6B//98fo+VFWXA/vN0v41uvEnkiRJP2G+B/+9C/hgVX1+oHnl6EuSJEnTbL4xJ18F/iLJjUn+LMlPHAWRJEna1OYMJ1X1jqp6OvBs4FvAyUmuSfJ/kzxubBVKkqSpMuxN2N5aVfsBLwOOoLvqRpIkaZMb5iZsWyR5YZJT6W6+di3wSyOvTJIkTaX5BsQ+h+5IyWHAl4DTgGOqapjLiCVJkjbKfJcSHwd8EPjdqrpzTPVIkqQpN2c4qaqDxlmIJEkSDHf7ekmSpLExnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJktQUw4kkSWqK4USSJDXFcCJJkppiOJEkSU0ZWThJ8qgk5yf5SpKrkry2b98xyblJruvfdxhVDZIkafKM8sjJvcDvVtXewP7Aq5PsDRwLnFdVewLn9fOSJEnACMNJVd1WVRf30/cAVwO7AIcDK/vVVgJHjKoGSZI0ecYy5iTJbsB+wIXA0qq6rV90O7B0HDVIkqTJMPJwkmQb4J+A11XV3YPLqqqAmuNzxyRZlWTV2rVrR12mJElqxEjDSZKH0AWTU6vqo33zHUl27pfvDKyZ7bNVdVJVLa+q5UuWLBllmZIkqSGjvFonwPuAq6vqbQOLzgJW9NMrgDNHVYMkSZo8W4xw288EXgFckeTSvu3/ACcCH07yKuAm4CUjrEGSJE2YkYWTqvockDkWHzyq75UkSZPNO8RKkqSmGE4kSVJTDCeSJKkphhNJktQUw4kkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJktQUw4kkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiSpKaMLJwkOTnJmiRXDrTtmOTcJNf17zuM6vslSdJkGuWRk1OAQ9dpOxY4r6r2BM7r5yVJku43snBSVRcA316n+XBgZT+9EjhiVN8vSZIm07jHnCytqtv66duBpWP+fkmS1LgFGxBbVQXUXMuTHJNkVZJVa9euHWNlkiRpIY07nNyRZGeA/n3NXCtW1UlVtbyqli9ZsmRsBUqSpIU17nByFrCin14BnDnm75ckSY0b5aXEHwK+AOyV5JYkrwJOBJ6T5DrgkH5ekiTpfluMasNV9bI5Fh08qu+UJEmTzzvESpKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJktQUw4kkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmLEg4SXJokmuTXJ/k2IWoQZIktWns4STJ5sC7gOcDewMvS7L3uOuQJEltWogjJ08Drq+qr1XVj4DTgMMXoA5JktSghQgnuwDfGJi/pW+TJEkiVTXeL0yOBA6tql/v518B/HxVvWad9Y4Bjuln9wKuHVFJjwS+OaJtt2Ra+gnT09dp6SdMT1+npZ8wPX2dln7ChvV116paMtfCLTZNPRtkNfCogfllfdt/UVUnASeNupgkq6pq+ai/Z6FNSz9hevo6Lf2E6enrtPQTpqev09JP2LR9XYjTOl8G9kyye5ItgaOAsxagDkmS1KCxHzmpqnuTvAb4FLA5cHJVXTXuOiRJUpsW4rQOVfVJ4JML8d2zGPmpo0ZMSz9hevo6Lf2E6enrtPQTpqev09JP2IR9HfuAWEmSpPl4+3pJktSUqQ0ni/EW+kluTHJFkkuTrOrbdkxybpLr+vcd+vYkeWff/8uTPHlhq59bkpOTrEly5UDbBvcryYp+/euSrFiIvqzPHH09Psnqfr9emuSwgWXH9X29NsnzBtqb/n0neVSS85N8JclVSV7bty+q/TpPPxfjPt06yZeSXNb39Y/69t2TXNjXfXp/IQRJturnr++X7zawrVn/DFoxT19PSfL1gf26b98+kb/fGUk2T3JJko/386Pfp1U1dS+6gbg3AI8BtgQuA/Ze6Lo2Qb9uBB65TtufAcf208cCb+2nDwPOBgLsD1y40PXP069nAU8GrtzYfgE7Al/r33fop3dY6L4N2dfjgd+bZd29+9/uVsDu/W9680n4fQM7A0/up7cFvtr3Z1Ht13n6uRj3aYBt+umHABf2++rDwFF9+3uA3+ynfwt4Tz99FHD6fH8GC92/Ift6CnDkLOtP5O93oP7XAx8EPt7Pj3yfTuuRk2m6hf7hwMp+eiVwxED7B6rzRWD7JDsvRIHrU1UXAN9ep3lD+/U84Nyq+nZV3QmcCxw6+uo3zBx9ncvhwGlV9cOq+jpwPd1vu/nfd1XdVlUX99P3AFfT3Sl6Ue3Xefo5l0nep1VV3+1nH9K/CjgIOKNvX3efzuzrM4CDk4S5/wyaMU9f5zKRv1+AJMuAXwTe28+HMezTaQ0ni/UW+gV8OslF6e6wC7C0qm7rp28HlvbTk/5nsKH9mvT+vqY/HHzyzKkOFklf+0O/+9H93+ei3a/r9BMW4T7tD/9fCqyh+4f2BuCuqrq3X2Ww7vv71C//DrATE9rXqprZr2/p9+vbk2zVt03yfv0r4A3Af/bzOzGGfTqt4WSxOqCqnkz3xOdXJ3nW4MLqjq8tusuzFmu/Brwb2APYF7gN+MuFLWfTSbIN8E/A66rq7sFli2m/ztLPRblPq+q+qtqX7s7fTwMev8Aljcy6fU3yROA4uj4/le5Uze8vYIkPWpIXAGuq6qJxf/e0hpOhbqE/aapqdf++BvgY3V8Od8ycrunf1/SrT/qfwYb2a2L7W1V39H8R/ifw9zxwOHSi+5rkIXT/YJ9aVR/tmxfdfp2tn4t1n86oqruA84Gn053CmLmn1mDd9/epX74d8C0mt6+H9qfxqqp+CLyfyd+vzwRelORGulOJBwHvYAz7dFrDyaK7hX6ShyfZdmYaeC5wJV2/ZkaArwDO7KfPAn6lH0W+P/CdgcPpk2BD+/Up4LlJdugPoT+3b2veOmOBXky3X6Hr61H9CPndgT2BLzEBv+/+PPT7gKur6m0DixbVfp2rn4t0ny5Jsn0//VDgOXRjbM4HjuxXW3efzuzrI4HP9EfL5vozaMYcfb1mIFiHbhzG4H6duN9vVR1XVcuqaje639xnqupoxrFP5xstu5hfdKOnv0p3TvSNC13PJujPY+hGQ18GXDXTJ7rzfecB1wH/CuzYtwd4V9//K4DlC92Hefr2IbpD3z+mO1f5qo3pF/BrdAOxrgdeudD92oC+/kPfl8v7/8h3Hlj/jX1frwWeP9De9O8bOIDulM3lwKX967DFtl/n6edi3KdPAi7p+3Ql8Ka+/TF0/xBdD3wE2Kpv37qfv75f/pj1/Rm08pqnr5/p9+uVwD/ywBU9E/n7XafPv8ADV+uMfJ96h1hJktSUaT2tI0mSGmU4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJpPsl+ekkpyW5oX9G0yeTPG7E3/kzSc7op/dNcthGbuf4JJXksQNtr+vblid5bZK/Glj2d0n+dWD+t5O888H0RdKmYTiRBNx/V8uPAf9WVXtU1VPonhWydP5PPjhVdWtVzdxtcl+6m41trCvo7mQ547/T3ZQQ4PPAMwaW7QNsl2Tzfv4ZwL8/iO+WtIkYTiTNOBD4cVW9Z6ahqi4DLklyXpKLk1yR5HDonrKb5Jokpya5OskZSR7WL3tTki8nuTLJSX3wIcljk/xrksv67e3Rb+fK/rbsfwy8NMmlSV6a5LokS/rPbpbk+pn5Ofwz3ePZSbIH3VNRv9kvuxR4XJKHJtkO+H7f9nP98mfQBRhJC8xwImnGE4HZnj76A+DF1T3x+kDgL2fCBrAX8LdV9QTgbuC3+va/qaqnVtUTgYcCL+jbTwXeVVX70IWB+5/nVFU/At4EnF5V+1bV6XS3AD+6X+UQ4LKqWjtPH+4GvpHuCbFHAacPbP9euluOPxXYH7gQ+CLwjCS7AKmqb/zkJiWNm+FE0voE+JMkl9M972YXHjjV842qmjna8I90z5IBODDJhUmuoHuS6c/2D6bcpao+BlBVP6iq/1jPd58M/Eo//Wt0T3pdn9PogskRdKepBv07XSh6BvCF/jUz7ykdqRGGE0kzrgKeMkv70cAS4ClVtS9wB90DvqB7qN2gSrI18LfAkVX1c8DfD6y/QfojGXckOYju8fNnD/GxjwOvAG6uqrvXWTYz7uTpdMHkamBvDCdSUwwnkmZ8BtgqyTEzDUmeBOwKrKmqHyc5sJ+f8egkT++nXw58jgeCyDeTbEP/aPWquge4JckR/ba3mhmjMuAeYNt12t5Ld1TmI1V13/o60R+N+X3gLbMs/gLdKZ0lVbWmuiefrqUbp+J4E6kRhhNJAPT/UL8YOKS/lPgq4E+BTwLL+1M0vwJcM/Cxa4FXJ7ka2AF4d1XdRXe05ErgU8CXB9Z/BfC/+1NE/w789DplnA/sPTMgtm87C9iG4U7pzPTltFE1doIAAACESURBVKq6eJb2O+nCyFUDzV8Afgq4bNjtSxqtdH8fSdKGSbIb8PF+0Osov2c58Paq+m+j/B5J7dhioQuQpLkkORb4TR64YkfSFPDIiaSJkuSNdDdXG/SRqpptjImkCWQ4kSRJTXFArCRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkpvx/oo8PSE2t03gAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 648x432 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Make the plot\n",
    "plt.figure(figsize=(9,6))\n",
    "plt.bar(new_x_group, height_group, width=width_group, edgecolor = \"black\")\n",
    "# Add title and axis names\n",
    "plt.title(selection)\n",
    "plt.xlabel('Capacity_MW')\n",
    "plt.ylabel('Variable Cost/Price')\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We have created the Big Coal supply curve. It shows the price of electricity, and the quantity supplied at those prices, which depends on variable cost. For example, at any variable cost equal to or above 36.5, the producer `FOUR CORNERS` (the one with the lowest production costs) will supply, and so on. Notably, we observe that the supply curve is also upward sloping since we need higher prices to entice producers with higher variasble costs to produce. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "tags": [
     "remove_input"
    ]
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Price: $30\n",
      "No production\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAicAAAGECAYAAADtMIAAAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAfJ0lEQVR4nO3df7xldV3v8dcbEFRAfo7jyBigIkmWoKOBejPEH2QqWFzEeNikdLlZdjVvt8bbzSi1sB+almmTokOhgCRBKBghyaNUcvglg4D8kAEGmDkqOGj+GvzcP/Y6sJvOOXNmmHXOd5/9ej4e+7HX+q611/58Z21m3nzXr1QVkiRJrdhhvguQJEkaZjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJktQUw4mk7SbJB5L87nzXsbkkH0ny9vmuQ9LsGE4kzVqS25J8J8m3ktyb5JNJnjC5vKp+pare9jC2/5IklyW5P8lEks8mecX2qV7SqDCcSNpaL6+q3YAlwHrgL7bHRpMcB3wcOB1YCiwG3gq8fHtsX9LoMJxI2iZV9V3gHOCQybbND58k+a0kdye5K8kvJ6kkT958W0kCvAt4W1V9sKq+WVU/rKrPVtX/6NbZIcn/S7I2yYYkpyfZY2gbH09yT5JvdqMvP9Zn/yX1x3AiaZskeTTwKuAL0yw/Gngz8ELgycBPz7C5g4EnMAg70/ml7nUk8ERgN+Avh5ZfCBwEPBa4Ejhji52Q1KSd5rsASSPnH5JsAnYFJoCXTLPe8cCHq+o6gCSnACdOs+4+3fvdM3zvicC7qurWbntvAdYkeW1Vbaqq0yZX7L7r3iR7VNU3Z9ctSa1w5ETS1jq2qvYEHgm8AfhsksdNsd7jgTuG5u+YYp1JX+/el8ywzuOBtUPzaxn8D9biJDsmOTXJLUk2Ard16+w7w/YkNcpwImmbVNUDVfUJ4AHgeVOscjeDE1snPWGKdSbdyCC8/PwM69wF7D80/yPAJgYn5f4CcAyDQ0h7AAd062SG7UlqlOFE0jbJwDHAXsD1U6xyNvDaJE/tzk+Z9v4nVVUMzk/53SSvTfKY7gTY5yVZ2a32MeA3khyYZDfgD4GzqmoTsDvwPQYjMI/ulkkaUYYTSVvrH5N8C9gIvANYPnleybCquhB4L3ApcDMPnTj7vak2WlXnMDjB9nUMRknWA28HzutWOQ34W+Ay4KvAd4Ff75adzuAwzzrgy0xzkq6k0ZDB/7BIUr+SPBVYA+zSjXZI0pQcOZHUmySvTLJLkr2AdwL/aDCRtCWGE0l9+p/ABuAWBifOvn5+y5E0CjysI0mSmuLIiSRJaorhRJIkNWUkbl+/77771gEHHDDfZUiSpO3giiuu+FpVLZpu+UiEkwMOOIDVq1fPdxmSJGk7SLJ2puUe1pEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJEmN22/p/iSZ09d+S/eft/7uNG/fLEmSZuWudbdz9Iqr5vQ7Lzr1sDn9vmGOnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSU3oLJ0kOTnL10Gtjkjcl2TvJxUlu6t736qsGSZI0enoLJ1V1Y1UdWlWHAs8E/gM4F1gBXFJVBwGXdPOSJEnA3B3WOQq4parWAscAq7r2VcCxc1SDJEkaAXMVTk4APtZNL66qu7vpe4DFc1SDJEkaAb2HkyQ7A68APr75sqoqoKb53MlJVidZPTEx0XOVkiSpFXMxcvIzwJVVtb6bX59kCUD3vmGqD1XVyqpaVlXLFi1aNAdlSpKkFsxFOHk1Dx3SATgfWN5NLwfOm4MaJEnSiOg1nCTZFXgR8Imh5lOBFyW5CXhhNy9JkgTATn1uvKq+DeyzWdvXGVy9I0mS9F94h1hJktQUw4kkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJKknuy3dH+SPOzXuOn12TqSJI2zu9bdztErrnrY27no1MO2QzWjw5ETSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJktQUw4kkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkpvQaTpLsmeScJDckuT7JEUn2TnJxkpu69736rEGSJI2WvkdO3gNcVFU/CjwduB5YAVxSVQcBl3TzkiRJQI/hJMkewE8BHwKoqu9X1X3AMcCqbrVVwLF91SBJkkZPnyMnBwITwIeTXJXkg0l2BRZX1d3dOvcAi6f6cJKTk6xOsnpiYqLHMiVJUkv6DCc7Ac8A3l9VhwHfZrNDOFVVQE314apaWVXLqmrZokWLeixTkiS1pM9wcidwZ1Vd3s2fwyCsrE+yBKB739BjDZIkacT0Fk6q6h7gjiQHd01HAV8GzgeWd23LgfP6qkGSJI2enXre/q8DZyTZGbgVeC2DQHR2kpOAtcDxPdcgSZJGSK/hpKquBpZNseioPr9XkiSNLu8QK0mSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJktQUw4kkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaslOfG09yG3A/8ACwqaqWJdkbOAs4ALgNOL6q7u2zDkmSNDrmYuTkyKo6tKqWdfMrgEuq6iDgkm5ekiQJmJ/DOscAq7rpVcCx81CDJElqVN/hpIB/SnJFkpO7tsVVdXc3fQ+weKoPJjk5yeokqycmJnouU5IktaLXc06A51XVuiSPBS5OcsPwwqqqJDXVB6tqJbASYNmyZVOuI0mSFp5eR06qal33vgE4F3g2sD7JEoDufUOfNUiSpNHSWzhJsmuS3SengRcDa4DzgeXdasuB8/qqQZIkjZ4+D+ssBs5NMvk9H62qi5J8ETg7yUnAWuD4HmuQJEkjprdwUlW3Ak+fov3rwFF9fa8kSRpt3iFWkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJktQUw4kkSWqK4USSJDXFcCJJkpqyxXCS5LmzaZMkSdoeZjNy8hezbJMkSXrYpn3wX5IjgOcAi5K8eWjRY4Ad+y5MkiSNp5meSrwzsFu3zu5D7RuB4/osSpIkja9pw0lVfRb4bJKPVNVagCQ7ALtV1ca5KlCSJI2X2Zxz8kdJHpNkV2AN8OUk/6fnuiRJ0piaTTg5pBspORa4EDgQeE2vVUmSpLE1m3DyiCSPYBBOzq+qHwDVb1mSJGlczSac/DVwG7ArcFmS/RmcFCtJkrTdzXS1DgBV9V7gvUNNa5Mc2V9JkiRpnM3mDrF7JHlXktXd688YjKJIkiRtd7M5rHMacD9wfPfaCHy4z6IkSdL42uJhHeBJVfXzQ/O/n+TqvgqSJEnjbTYjJ99J8rzJme6hf9/pryRJkjTOZjNy8ivA6Un26ObvBX6pt4okSdJYm83VOtcAT0/ymG7ey4glSVJvpj2sk+TNSU6anK+qjVW1MclJSd40N+VJkqRxM9M5JycCp0/R/rfA6/opR5IkjbuZDuvs1N2q/j+pqu8nSY81SdKc2W/p/ty17vb5LmNGO+38KDZ93+sQND5mCic7JFlcVeuHG5Ms7rkmSZozd627naNXXDXfZczoolMPa75GTe2iUw+b7xJG0kyHdf4E+GSS5yfZvXv9NHAB8KdzUp0kSRo7046cVNXpSSaAPwCexuBJxNcBb62qC+eoPkmSNGZmvJS4CyEGEUmSNGdmc4fYhyXJjkmuSnJBN39gksuT3JzkrCQ7912DJEkaHb2HE+CNwPVD8+8E3l1VT2Zwt9mTpvyUpAVrv6X7k6SJl6T2zOb29dssyVLgZ4F3AG/uLkF+AfAL3SqrgFOA9/dZh6S2tHSFjFdTSO3Z4shJksVJPpTkwm7+kOE7x27BnwO/Bfywm98HuK+qNnXzdwL7TfO9JydZnWT1xMTELL9OkiSNutkc1vkI8Gng8d38V4At3r4+ycuADVV1xbYUVlUrq2pZVS1btGjRtmxCkiSNoNmEk32r6my60Y9u1OOBWXzuucArktwGnMngcM57gD2TTB5OWgqs29qiJUnSwjWbcPLtJPswuM8JSQ4HvrmlD1XVW6pqaVUdAJwAfKaqTgQuBY7rVlsOnLcthUuSpIVpNifEvhk4H3hSkn8DFvFQuNgWvw2cmeTtwFXAhx7GtiRJ0gKzxXBSVVcmeT5wMBDgxqkeCLiFbfwL8C/d9K3As7e6UkmSNBamDSdJfm6aRU9JQlV9oqeaJEnSGJtp5OTlMywrwHAiSZK2u5ke/PfauSxEkiQJZncTtn2SvDfJlUmuSPKe7uodSZKk7W42lxKfCUwAP8/gKp0J4Kw+i5IkSeNrNpcSL6mqtw3Nvz3Jq/oqSJIkjbfZjJz8U5ITkuzQvY5ncDt7SZKk7W6mS4nvZ3BVThg8S+fvukU7AN8CfrP36iRJ0tiZ6Wqd3eeyEEmSJJjdOSck2Qs4CHjkZFtVXdZXUXNpv6X7c9e62+e7DEmS1NliOEnyy8AbGTxB+GrgcODzDJ4yPPLuWnc7R6+4ar7LkMbKRaceNt8lSGrYbE6IfSPwLGBtVR0JHAbc12tVkiRpbM0mnHy3qr4LkGSXqrqBwUMAJUmStrvZnHNyZ5I9gX8ALk5yL7C237IkSdK42mI4qapXdpOnJLkU2AO4qNeqJEnS2JrpPiePqaqNSfYear62e98N+EavlUmSpLE008jJR4GXAVfw0M3Yht+f2Ht1kiRp7Mx0E7aXJQnw/KryRiCSJGlOzHi1TlUV8Mk5qkWSJGlWlxJfmeRZvVciSZLE7C4l/kngxCRrgW/TnXNSVT/Ra2WSJGkszSacvKT3KiRJkjqzuc/JWoAkj2XowX+SJEl92OI5J0lekeQm4KvAZ4HbgAt7rkuSJI2p2ZwQ+zYGTyL+SlUdCBwFfKHXqiRJ0tiaTTj5QVV9HdghyQ5VdSmwrOe6JEnSmJrNCbH3JdkNuAw4I8kGBlftSJIkbXezGTk5BvgO8BsMHvh3C/DyPouSJEnja6YH/70P+GhV/dtQ86r+S5IkSeNsppGTrwB/muS2JH+c5LC5KkqSJI2vacNJVb2nqo4Ang98HTgtyQ1Jfi/JU+asQkmSNFa2eM5JVa2tqndW1WHAq4Fjget7r0ySJI2l2dyEbackL09yBoObr90I/FzvlUmSpLE00wmxL2IwUvJS4N+BM4GTq2pWlxEneSSDy4936b7nnKr6vSQHdtvaB7gCeE1Vff9h9UKSJC0YM42cvAX4HPDUqnpFVX10tsGk8z3gBVX1dOBQ4OgkhwPvBN5dVU8G7gVO2sbaJUnSApSq6v9LkkcD/wq8Hvgk8Liq2pTkCOCUqprxycfLklrde5WSJGkuBK6oqmnvNj+bm7Bt+5cnOya5GtgAXMzgBm73VdWmbpU7gf2m+ezJSVYnMZdIkjRGZnP7+m1WVQ8AhybZEzgX+NGt+OxKYCXAsmXLitX9ZJQkHL3iql62LWlqF516WDP/3bVUy3RGoUZNbXvtu/n4DVx06mH0dnQlmXFxryMnk6rqPuBS4AhgzySToWgpsG4uapAkSaOht3CSZFE3YkKSRwEvYnB/lEuB47rVlgPn9VWDJEkaPX0e1lkCrEqyI4MQdHZVXZDky8CZSd4OXAV8qMcaJEnSiOktnFTVl4D/8jyeqroVeHZf3ytJkkbbnJxzIkmSNFuGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJktQUw4kkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKb2FkyRPSHJpki8nuS7JG7v2vZNcnOSm7n2vvmqQJEmjp8+Rk03A/66qQ4DDgV9LcgiwArikqg4CLunmJUmSgB7DSVXdXVVXdtP3A9cD+wHHAKu61VYBx/ZVgyRJGj1zcs5JkgOAw4DLgcVVdXe36B5g8VzUIEmSRkPv4STJbsDfA2+qqo3Dy6qqgJrmcycnWZ1k9cTERN9lSpKkRvQaTpI8gkEwOaOqPtE1r0+ypFu+BNgw1WeramVVLauqZYsWLeqzTEmS1JA+r9YJ8CHg+qp619Ci84Hl3fRy4Ly+apAkSaNnpx63/VzgNcC1Sa7u2v4vcCpwdpKTgLXA8T3WIEmSRkxv4aSq/hXINIuP6ut7JUnSaPMOsZIkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJktQUw4kkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKb2FkySnJdmQZM1Q295JLk5yU/e+V1/fL0mSRlOfIycfAY7erG0FcElVHQRc0s1LkiQ9qLdwUlWXAd/YrPkYYFU3vQo4tq/vlyRJo2muzzlZXFV3d9P3AIvn+PslSVLj5u2E2KoqoKZbnuTkJKuTrJ6YmJjDyiRJ0nya63CyPskSgO59w3QrVtXKqlpWVcsWLVo0ZwVKkqT5Ndfh5HxgeTe9HDhvjr9fkiQ1rs9LiT8GfB44OMmdSU4CTgVelOQm4IXdvCRJ0oN26mvDVfXqaRYd1dd3SpKk0ecdYiVJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJktQUw4kkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSU+YlnCQ5OsmNSW5OsmI+apAkSW2a83CSZEfgfcDPAIcAr05yyFzXIUmS2jQfIyfPBm6uqlur6vvAmcAx81CHJElq0HyEk/2AO4bm7+zaJEmSSFXN7RcmxwFHV9Uvd/OvAX6yqt6w2XonAyd3swcDN/ZU0r7A13radkvGpZ8wPn0dl37C+PR1XPoJ49PXceknbF1f96+qRdMt3Gn71LNV1gFPGJpf2rX9J1W1EljZdzFJVlfVsr6/Z76NSz9hfPo6Lv2E8enruPQTxqev49JP2L59nY/DOl8EDkpyYJKdgROA8+ehDkmS1KA5Hzmpqk1J3gB8GtgROK2qrpvrOiRJUpvm47AOVfUp4FPz8d1T6P3QUSPGpZ8wPn0dl37C+PR1XPoJ49PXceknbMe+zvkJsZIkSTPx9vWSJKkpYxtOFuIt9JPcluTaJFcnWd217Z3k4iQ3de97de1J8t6u/19K8oz5rX56SU5LsiHJmqG2re5XkuXd+jclWT4ffdmSafp6SpJ13X69OslLh5a9pevrjUleMtTe9O87yROSXJrky0muS/LGrn1B7dcZ+rkQ9+kjk/x7kmu6vv5+135gksu7us/qLoQgyS7d/M3d8gOGtjXln0ErZujrR5J8dWi/Htq1j+Tvd1KSHZNcleSCbr7/fVpVY/dicCLuLcATgZ2Ba4BD5ruu7dCv24B9N2v7Y2BFN70CeGc3/VLgQiDA4cDl813/DP36KeAZwJpt7RewN3Br975XN73XfPdtln09BfjNKdY9pPvt7gIc2P2mdxyF3zewBHhGN7078JWuPwtqv87Qz4W4TwPs1k0/Ari821dnAyd07R8AXt9N/yrwgW76BOCsmf4M5rt/s+zrR4Djplh/JH+/Q/W/GfgocEE33/s+HdeRk3G6hf4xwKpuehVw7FD76TXwBWDPJEvmo8AtqarLgG9s1ry1/XoJcHFVfaOq7gUuBo7uv/qtM01fp3MMcGZVfa+qvgrczOC33fzvu6rurqoru+n7gesZ3Cl6Qe3XGfo5nVHep1VV3+pmH9G9CngBcE7Xvvk+ndzX5wBHJQnT/xk0Y4a+Tmckf78ASZYCPwt8sJsPc7BPxzWcLNRb6BfwT0muyOAOuwCLq+rubvoeYHE3Pep/Blvbr1Hv7xu64eDTJg91sED62g39Hsbg/z4X7H7drJ+wAPdpN/x/NbCBwT+0twD3VdWmbpXhuh/sU7f8m8A+jGhfq2pyv76j26/vTrJL1zbK+/XPgd8CftjN78Mc7NNxDScL1fOq6hkMnvj8a0l+anhhDcbXFtzlWQu1X0PeDzwJOBS4G/iz+S1n+0myG/D3wJuqauPwsoW0X6fo54Lcp1X1QFUdyuDO388GfnSeS+rN5n1N8jTgLQz6/CwGh2p+ex5LfNiSvAzYUFVXzPV3j2s4mdUt9EdNVa3r3jcA5zL4y2H95OGa7n1Dt/qo/xlsbb9Gtr9Vtb77i/CHwN/w0HDoSPc1ySMY/IN9RlV9omtecPt1qn4u1H06qaruAy4FjmBwCGPynlrDdT/Yp275HsDXGd2+Ht0dxquq+h7wYUZ/vz4XeEWS2xgcSnwB8B7mYJ+OazhZcLfQT7Jrkt0np4EXA2sY9GvyDPDlwHnd9PnAL3ZnkR8OfHNoOH0UbG2/Pg28OMle3RD6i7u25m12LtArGexXGPT1hO4M+QOBg4B/ZwR+391x6A8B11fVu4YWLaj9Ol0/F+g+XZRkz276UcCLGJxjcylwXLfa5vt0cl8fB3ymGy2b7s+gGdP09YahYB0G52EM79eR+/1W1VuqamlVHcDgN/eZqjqRudinM50tu5BfDM6e/gqDY6K/M9/1bIf+PJHB2dDXANdN9onB8b5LgJuAfwb27toDvK/r/7XAsvnuwwx9+xiDoe8fMDhWedK29At4HYMTsW4GXjvf/dqKvv5t15cvdf+RLxla/3e6vt4I/MxQe9O/b+B5DA7ZfAm4unu9dKHt1xn6uRD36U8AV3V9WgO8tWt/IoN/iG4GPg7s0rU/spu/uVv+xC39GbTymqGvn+n26xrg73joip6R/P1u1uef5qGrdXrfp94hVpIkNWVcD+tIkqRGGU4kSVJTDCeSJKkphhNJktQUw4kkSWqK4USSJDXFcCLpQUkel+TMJLd0z2j6VJKn9Pydj09yTjd9aJKXbuN2TklSSZ481Pamrm1Zkjcm+fOhZX+d5J+H5n89yXsfTl8kbR+GE0nAg3e1PBf4l6p6UlU9k8GzQhbP/MmHp6ruqqrJu00eyuBmY9vqWgZ3spz03xnclBDg34DnDC17OrBHkh27+ecAn3sY3y1pOzGcSJp0JPCDqvrAZENVXQNcleSSJFcmuTbJMTB4ym6SG5KckeT6JOckeXS37K1JvphkTZKVXfAhyZOT/HOSa7rtPanbzprutux/ALwqydVJXpXkpiSLus/ukOTmyflp/AODx7OT5EkMnor6tW7Z1cBTkjwqyR7Ad7q2H++WP4dBgJE0zwwnkiY9DZjq6aPfBV5ZgydeHwn82WTYAA4G/qqqngpsBH61a//LqnpWVT0NeBTwsq79DOB9VfV0BmHgwec5VdX3gbcCZ1XVoVV1FoNbgJ/YrfJC4JqqmpihDxuBOzJ4QuwJwFlD29/E4JbjzwIOBy4HvgA8J8l+QKrqjv+6SUlzzXAiaUsC/GGSLzF43s1+PHSo546qmhxt+DsGz5IBODLJ5UmuZfAk0x/rHky5X1WdC1BV362q/9jCd58G/GI3/ToGT3rdkjMZBJNjGRymGvY5BqHoOcDnu9fkvId0pEYYTiRNug545hTtJwKLgGdW1aHAegYP+ILBQ+2GVZJHAn8FHFdVPw78zdD6W6UbyVif5AUMHj9/4Sw+dgHwGuD2qtq42bLJ806OYBBMrgcOwXAiNcVwImnSZ4Bdkpw82ZDkJ4D9gQ1V9YMkR3bzk34kyRHd9C8A/8pDQeRrSXaje7R6Vd0P3Jnk2G7bu0yeozLkfmD3zdo+yGBU5uNV9cCWOtGNxvw28I4pFn+ewSGdRVW1oQZPPp1gcJ6K55tIjTCcSAKg+4f6lcALu0uJrwP+CPgUsKw7RPOLwA1DH7sR+LUk1wN7Ae+vqvsYjJasAT4NfHFo/dcA/6s7RPQ54HGblXEpcMjkCbFd2/nAbszukM5kX86sqiunaL+XQRi5bqj588BjgWtmu31J/crg7yNJ2jpJDgAu6E567fN7lgHvrqr/1uf3SGrHTvNdgCRNJ8kK4PU8dMWOpDHgyImkkZLkdxjcXG3Yx6tqqnNMJI0gw4kkSWqKJ8RKkqSmGE4kSVJTDCeSJKkphhNJktQUw4kkSWrK/wcx8YPQWPV5GwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 648x432 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "group_plot(30)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "tags": [
     "remove_input"
    ]
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Price: $37\n",
      "Total Production/Market Supply:  1900\n",
      "\n",
      "Suppliers:  ['FOUR CORNERS']\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAicAAAGECAYAAADtMIAAAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAfQElEQVR4nO3de7gddX3v8fcHAmgBuW7TSCxBRSq1lWi0oJ5awAu1KthyFMtjU6WH04s9Wk9PG09Pra3aYi9aba2WKhpaFJBqoVSwFKk8rYqGm4CAXCRAuGSrYNB6A7/njzUbVuPeOzshs/dv7fV+Pc961sxvZs36/jIL8snMb2ZSVUiSJLVih4UuQJIkaZjhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkrabJO9N8nsLXcfmknwwyVsWug5Jc2M4kTRnSW5J8q0k30hyT5J/TvLYqeVV9StV9eaHsf0XJLk4yX1JJpN8KslLtk/1kkaF4UTS1npxVe0GLAPuBv5ye2w0ybHAR4BTgeXAUuCNwIu3x/YljQ7DiaRtUlXfBs4CDp5q2/z0SZLfTnJnkjuS/HKSSvKEzbeVJMDbgTdX1fuq6utV9f2q+lRV/Y9unR2S/L8k65NsTHJqkj2GtvGRJHcl+Xp39OXH+uy/pP4YTiRtkyQ/BLwc+OwMy48CXg88F3gC8NOzbO4g4LEMws5Mfql7HQ48DtgN+Kuh5ecBBwKPBi4DTttiJyQ1aclCFyBp5PxjkvuBXYFJ4AUzrPcy4ANVdQ1AkjcBx8+w7j7d+52zfO/xwNur6uZue28Ark7yqqq6v6pOmVqx+657kuxRVV+fW7cktcIjJ5K21jFVtSfwCOA1wKeS/PA06z0GuG1o/rZp1pny1e592SzrPAZYPzS/nsE/sJYm2THJSUluSrIJuKVbZ99ZtiepUYYTSdukqh6oqo8CDwDPnmaVOxkMbJ3y2GnWmXI9g/Dy87Oscwew/9D8jwD3MxiU+wvA0QxOIe0BrOjWySzbk9Qow4mkbZKBo4G9gGunWeVM4FVJntSNT5nx/idVVQzGp/xeklcleVQ3APbZSU7uVvsw8JtJDkiyG/BHwBlVdT+wO/AdBkdgfqhbJmlEGU4kba1/SvINYBPwVmD11LiSYVV1HvAu4CLgRh4aOPud6TZaVWcxGGD7agZHSe4G3gKc3a1yCvB3wMXAl4FvA7/RLTuVwWmeDcAXmWGQrqTRkME/WCSpX0meBFwN7NId7ZCkaXnkRFJvkrw0yS5J9gLeBvyTwUTSlhhOJPXpfwIbgZsYDJz91YUtR9Io8LSOJElqikdOJElSUwwnkiSpKSNx+/p99923VqxYsdBlSJKk7eDSSy/9SlVNzLR8JMLJihUrWLdu3UKXIUmStoMk62db7mkdSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJktQUw4kkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiS1Lj9lu9Pknl97bd8/wXr75IF+2ZJkjQnd2y4laPWXD6v33n+SSvn9fuGeeREkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJktQUw4kkSWqK4USSJDXFcCJJkprSWzhJclCSK4Zem5K8LsneSS5IckP3vldfNUiSpNHTWzipquur6pCqOgR4GvCfwMeANcCFVXUgcGE3L0mSBMzfaZ0jgZuqaj1wNLC2a18LHDNPNUiSpBEwX+HkOODD3fTSqrqzm74LWDpPNUiSpBHQezhJsjPwEuAjmy+rqgJqhs+dmGRdknWTk5M9VylJkloxH0dOfga4rKru7ubvTrIMoHvfON2HqurkqlpVVasmJibmoUxJktSC+Qgnr+ChUzoA5wCru+nVwNnzUIMkSRoRvYaTJLsCzwM+OtR8EvC8JDcAz+3mJUmSAFjS58ar6pvAPpu1fZXB1TuSJEk/wDvESpKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJPVkv+X7k+Rhv8ZNr8/WkSRpnN2x4VaOWnP5w97O+Set3A7VjA6PnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJktQUw4kkSWqK4USSJDWl13CSZM8kZyW5Lsm1SQ5LsneSC5Lc0L3v1WcNkiRptPR95OSdwPlV9aPAU4BrgTXAhVV1IHBhNy9JkgT0GE6S7AH8FPB+gKr6blXdCxwNrO1WWwsc01cNkiRp9PR55OQAYBL4QJLLk7wvya7A0qq6s1vnLmDpdB9OcmKSdUnWTU5O9limJElqSZ/hZAnwVOA9VbUS+CabncKpqgJqug9X1clVtaqqVk1MTPRYpiRJakmf4eR24PaquqSbP4tBWLk7yTKA7n1jjzVIkqQR01s4qaq7gNuSHNQ1HQl8ETgHWN21rQbO7qsGSZI0epb0vP3fAE5LsjNwM/AqBoHozCQnAOuBl/VcgyRJGiG9hpOqugJYNc2iI/v8XkmSNLq8Q6wkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJktQUw4kkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJasqSPjee5BbgPuAB4P6qWpVkb+AMYAVwC/CyqrqnzzokSdLomI8jJ4dX1SFVtaqbXwNcWFUHAhd285IkScDCnNY5GljbTa8FjlmAGiRJUqP6DicF/EuSS5Oc2LUtrao7u+m7gKXTfTDJiUnWJVk3OTnZc5mSJKkVvY45AZ5dVRuSPBq4IMl1wwurqpLUdB+sqpOBkwFWrVo17TqSJGnx6fXISVVt6N43Ah8DngHcnWQZQPe+sc8aJEnSaOktnCTZNcnuU9PA84GrgXOA1d1qq4Gz+6pBkiSNnj5P6ywFPpZk6ns+VFXnJ/k8cGaSE4D1wMt6rEGSJI2Y3sJJVd0MPGWa9q8CR/b1vZIkabR5h1hJktQUw4kkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJasoWw0mSZ82lTZIkaXuYy5GTv5xjmyRJ0sM244P/khwGPBOYSPL6oUWPAnbsuzBJkjSeZnsq8c7Abt06uw+1bwKO7bMoSZI0vmYMJ1X1KeBTST5YVesBkuwA7FZVm+arQEmSNF7mMubkj5M8KsmuwNXAF5P8n57rkiRJY2ou4eTg7kjJMcB5wAHAK3utSpIkja25hJOdkuzEIJycU1XfA6rfsiRJ0riaSzj5G+AWYFfg4iT7MxgUK0mStN3NdrUOAFX1LuBdQ03rkxzeX0mSJGmczeUOsXskeXuSdd3rzxkcRZEkSdru5nJa5xTgPuBl3WsT8IE+i5IkSeNri6d1gMdX1c8Pzf9Bkiv6KkiSJI23uRw5+VaSZ0/NdA/9+1Z/JUmSpHE2lyMnvwKcmmSPbv4e4Jd6q0iSJI21uVytcyXwlCSP6ua9jFiSJPVmxtM6SV6f5ISp+araVFWbkpyQ5HXzU54kSRo3s405OR44dZr2vwNe3U85kiRp3M12WmdJd6v6/6KqvpskPdYkSfNmv+X7c8eGWxe6jFkt2fmR3P9dr0PQ+JgtnOyQZGlV3T3cmGRpzzVJ0ry5Y8OtHLXm8oUuY1bnn7Sy+Ro1vfNPWrnQJYyk2U7r/Cnwz0mek2T37vXTwLnAn81LdZIkaezMeOSkqk5NMgn8IfBkBk8ivgZ4Y1WdN0/1SZKkMTPrpcRdCDGISJKkeTOXO8Q+LEl2THJ5knO7+QOSXJLkxiRnJNm57xokSdLo6D2cAK8Frh2afxvwjqp6AoO7zZ4w7ackLVr7Ld+fJE28JLVnLrev32ZJlgM/C7wVeH13CfIRwC90q6wF3gS8p886JLWlpStkvJpCas8Wj5wkWZrk/UnO6+YPHr5z7Bb8BfDbwPe7+X2Ae6vq/m7+dmC/Gb73xCTrkqybnJyc49dJkqRRN5fTOh8EPgE8ppv/ErDF29cneRGwsaou3ZbCqurkqlpVVasmJia2ZROSJGkEzSWc7FtVZ9Id/eiOejwwh889C3hJkluA0xmcznknsGeSqdNJy4ENW1u0JElavOYSTr6ZZB8G9zkhyaHA17f0oap6Q1Utr6oVwHHAJ6vqeOAi4NhutdXA2dtSuCRJWpzmMiD29cA5wOOT/AcwwUPhYlv8DnB6krcAlwPvfxjbkiRJi8wWw0lVXZbkOcBBQIDrp3sg4Ba28W/Av3XTNwPP2OpKJUnSWJgxnCT5uRkWPTEJVfXRnmqSJEljbLYjJy+eZVkBhhNJkrTdzfbgv1fNZyGSJEkwt5uw7ZPkXUkuS3Jpknd2V+9IkiRtd3O5lPh0YBL4eQZX6UwCZ/RZlCRJGl9zuZR4WVW9eWj+LUle3ldBkiRpvKWqZl8heTvwOeDMrulY4BlV9Vs91/agVUmtm68vkyRJvQpcWlWrZlw+UzhJch+Dq3IC7MpDD+/bAfhGVT1qO9c6I8OJJEmLx5bCyWxX6+zeT0nb4GlPg3XGE2mxSMJRay5f6DIAOP+klc3UMpNRqFHT2177biF+A+eftJItnV3ZZsmsi+cy5oQkewEHAo+Yaquqix9WYY3Yb/n+3LHh1oUuQ5IkdbYYTpL8MvBaBk8QvgI4FPgMg6cMj7w7Ntzqv0ikeXb+SSsXugRJDZvLpcSvBZ4OrK+qw4GVwL29ViVJksbWXMLJt6vq2wBJdqmq6xg8BFCSJGm7m8uYk9uT7An8I3BBknuA9f2WJUmSxtUWw0lVvbSbfFOSi4A9gPN7rUqSJI2tGcNJkkdV1aYkew81X9W97wZ8rdfKJEnSWJrtyMmHgBcBl/LQzdiG3x/Xe3WSJGnszHYTthclCfCcqvJGIJIkaV7MerVODW4N98/zVIskSdKcLiW+LMnTe69EkiSJuV1K/JPA8UnWA9+kG3NSVT/Ra2WSJGkszSWcvKD3KiRJkjpzuc/JeoAkj2bowX+SJEl92OKYkyQvSXID8GXgU8AtwHk91yVJksbUXAbEvpnBk4i/VFUHAEcCn+21KkmSNLbmEk6+V1VfBXZIskNVXQSs6rkuSZI0puYyIPbeJLsBFwOnJdnI4KodSZKk7W4uR06OBr4F/CaDB/7dBLy4z6IkSdL4mu3Bf+8GPlRV/zHUvLb/kiRJ0jib7cjJl4A/S3JLkj9JsnK+ipIkSeNrxnBSVe+sqsOA5wBfBU5Jcl2S30/yxHmrUJIkjZUtjjmpqvVV9baqWgm8AjgGuLb3yiRJ0liay03YliR5cZLTGNx87Xrg53qvTJIkjaXZBsQ+j8GRkhcCnwNOB06sqjldRpzkEQwuP96l+56zqur3kxzQbWsf4FLglVX13YfVC0mStGjMduTkDcCngSdV1Uuq6kNzDSad7wBHVNVTgEOAo5IcCrwNeEdVPQG4BzhhG2uXJEmL0GwDYo+oqvdV1T3bsuEa+EY3u1P3KuAI4KyufS2DMSySJEnA3G7Cts2S7JjkCmAjcAGDG7jdW1X3d6vcDuw3w2dPTLIuybrJyck+y5QkSQ3pNZxU1QNVdQiwHHgG8KNb8dmTq2pVVa2amJjorUZJktSWXsPJlKq6F7gIOAzYM8nUQNzlwIb5qEGSJI2G3sJJkokke3bTjwSex+D+KBcBx3arrQbO7qsGSZI0eubyVOJttQxYm2RHBiHozKo6N8kXgdOTvAW4HHh/jzVIkqQR01s4qaovAD/wPJ6qupnB+BNJkqQfMC9jTiRJkubKcCJJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJktQUw4kkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNaW3cJLksUkuSvLFJNckeW3XvneSC5Lc0L3v1VcNkiRp9PR55OR+4H9X1cHAocCvJzkYWANcWFUHAhd285IkSUCP4aSq7qyqy7rp+4Brgf2Ao4G13WprgWP6qkGSJI2eeRlzkmQFsBK4BFhaVXd2i+4Cls5HDZIkaTT0Hk6S7Ab8A/C6qto0vKyqCqgZPndiknVJ1k1OTvZdpiRJakSv4STJTgyCyWlV9dGu+e4ky7rly4CN0322qk6uqlVVtWpiYqLPMiVJUkP6vFonwPuBa6vq7UOLzgFWd9OrgbP7qkGSJI2eJT1u+1nAK4GrklzRtf1f4CTgzCQnAOuBl/VYgyRJGjG9hZOq+ncgMyw+sq/vlSRJo807xEqSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJktQUw4kkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkpvQWTpKckmRjkquH2vZOckGSG7r3vfr6fkmSNJr6PHLyQeCozdrWABdW1YHAhd28JEnSg3oLJ1V1MfC1zZqPBtZ202uBY/r6fkmSNJrme8zJ0qq6s5u+C1g6z98vSZIat2ADYquqgJppeZITk6xLsm5ycnIeK5MkSQtpvsPJ3UmWAXTvG2dasapOrqpVVbVqYmJi3gqUJEkLa77DyTnA6m56NXD2PH+/JElqXJ+XEn8Y+AxwUJLbk5wAnAQ8L8kNwHO7eUmSpAct6WvDVfWKGRYd2dd3SpKk0ecdYiVJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJktQUw4kkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUxYknCQ5Ksn1SW5MsmYhapAkSW2a93CSZEfg3cDPAAcDr0hy8HzXIUmS2rQQR06eAdxYVTdX1XeB04GjF6AOSZLUoIUIJ/sBtw3N3961SZIkkaqa3y9MjgWOqqpf7uZfCfxkVb1ms/VOBE7sZg8Cru+ppH2Br/S07ZaMSz9hfPo6Lv2E8enruPQTxqev49JP2Lq+7l9VEzMtXLJ96tkqG4DHDs0v79r+i6o6GTi572KSrKuqVX1/z0Ibl37C+PR1XPoJ49PXceknjE9fx6WfsH37uhCndT4PHJjkgCQ7A8cB5yxAHZIkqUHzfuSkqu5P8hrgE8COwClVdc181yFJktq0EKd1qKqPAx9fiO+eRu+njhoxLv2E8enruPQTxqev49JPGJ++jks/YTv2dd4HxEqSJM3G29dLkqSmjG04WYy30E9yS5KrklyRZF3XtneSC5Lc0L3v1bUnybu6/n8hyVMXtvqZJTklycYkVw+1bXW/kqzu1r8hyeqF6MuWzNDXNyXZ0O3XK5K8cGjZG7q+Xp/kBUPtTf++kzw2yUVJvpjkmiSv7doX1X6dpZ+LcZ8+IsnnklzZ9fUPuvYDklzS1X1GdyEESXbp5m/slq8Y2ta0fwatmKWvH0zy5aH9ekjXPpK/3ylJdkxyeZJzu/n+92lVjd2LwUDcm4DHATsDVwIHL3Rd26FftwD7btb2J8CabnoN8LZu+oXAeUCAQ4FLFrr+Wfr1U8BTgau3tV/A3sDN3fte3fReC923Ofb1TcBvTbPuwd1vdxfggO43veMo/L6BZcBTu+ndgS91/VlU+3WWfi7GfRpgt256J+CSbl+dCRzXtb8X+NVu+teA93bTxwFnzPZnsND9m2NfPwgcO836I/n7Har/9cCHgHO7+d736bgeORmnW+gfDaztptcCxwy1n1oDnwX2TLJsIQrckqq6GPjaZs1b268XABdU1deq6h7gAuCo/qvfOjP0dSZHA6dX1Xeq6svAjQx+283/vqvqzqq6rJu+D7iWwZ2iF9V+naWfMxnlfVpV9Y1udqfuVcARwFld++b7dGpfnwUcmSTM/GfQjFn6OpOR/P0CJFkO/Czwvm4+zMM+HddwslhvoV/AvyS5NIM77AIsrao7u+m7gKXd9Kj/GWxtv0a9v6/pDgefMnWqg0XS1+7Q70oG//pctPt1s37CItyn3eH/K4CNDP6ivQm4t6ru71YZrvvBPnXLvw7sw4j2taqm9utbu/36jiS7dG2jvF//Avht4Pvd/D7Mwz4d13CyWD27qp7K4InPv57kp4YX1uD42qK7PGux9mvIe4DHA4cAdwJ/vrDlbD9JdgP+AXhdVW0aXraY9us0/VyU+7SqHqiqQxjc+fsZwI8ucEm92byvSZ4MvIFBn5/O4FTN7yxgiQ9bkhcBG6vq0vn+7nENJ3O6hf6oqaoN3ftG4GMM/udw99Tpmu59Y7f6qP8ZbG2/Rra/VXV39z/C7wN/y0OHQ0e6r0l2YvAX9mlV9dGuedHt1+n6uVj36ZSquhe4CDiMwSmMqXtqDdf9YJ+65XsAX2V0+3pUdxqvquo7wAcY/f36LOAlSW5hcCrxCOCdzMM+HddwsuhuoZ9k1yS7T00DzweuZtCvqRHgq4Gzu+lzgF/sRpEfCnx96HD6KNjafn0CeH6SvbpD6M/v2pq32ViglzLYrzDo63HdCPkDgAOBzzECv+/uPPT7gWur6u1DixbVfp2pn4t0n04k2bObfiTwPAZjbC4Cju1W23yfTu3rY4FPdkfLZvozaMYMfb1uKFiHwTiM4f06cr/fqnpDVS2vqhUMfnOfrKrjmY99Otto2cX8YjB6+ksMzon+7kLXsx368zgGo6GvBK6Z6hOD830XAjcA/wrs3bUHeHfX/6uAVQvdh1n69mEGh76/x+Bc5Qnb0i/g1QwGYt0IvGqh+7UVff27ri9f6P4jXza0/u92fb0e+Jmh9qZ/38CzGZyy+QJwRfd64WLbr7P0czHu058ALu/6dDXwxq79cQz+IroR+AiwS9f+iG7+xm7547b0Z9DKa5a+frLbr1cDf89DV/SM5O93sz7/NA9drdP7PvUOsZIkqSnjelpHkiQ1ynAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJD0ryw0lOT3JT94ymjyd5Ys/f+ZgkZ3XThyR54TZu501JKskThtpe17WtSvLaJH8xtOxvkvzr0PxvJHnXw+mLpO3DcCIJePCulh8D/q2qHl9VT2PwrJCls3/y4amqO6pq6m6ThzC42di2uorBnSyn/HcGNyUE+A/gmUPLngLskWTHbv6ZwKcfxndL2k4MJ5KmHA58r6reO9VQVVcClye5MMllSa5KcjQMnrKb5LokpyW5NslZSX6oW/bGJJ9PcnWSk7vgQ5InJPnXJFd223t8t52ru9uy/yHw8iRXJHl5khuSTHSf3SHJjVPzM/hHBo9nJ8njGTwV9SvdsiuAJyZ5ZJI9gG91bT/eLX8mgwAjaYEZTiRNeTIw3dNHvw28tAZPvD4c+POpsAEcBPx1VT0J2AT8Wtf+V1X19Kp6MvBI4EVd+2nAu6vqKQzCwIPPc6qq7wJvBM6oqkOq6gwGtwA/vlvlucCVVTU5Sx82Abdl8ITY44AzhrZ/P4Nbjj8dOBS4BPgs8Mwk+wGpqtt+cJOS5pvhRNKWBPijJF9g8Lyb/XjoVM9tVTV1tOHvGTxLBuDwJJckuYrBk0x/rHsw5X5V9TGAqvp2Vf3nFr77FOAXu+lXM3jS65acziCYHMPgNNWwTzMIRc8EPtO9puY9pSM1wnAiaco1wNOmaT8emACeVlWHAHczeMAXDB5qN6ySPAL4a+DYqvpx4G+H1t8q3ZGMu5McweDx8+fN4WPnAq8Ebq2qTZstmxp3chiDYHItcDCGE6kphhNJUz4J7JLkxKmGJD8B7A9srKrvJTm8m5/yI0kO66Z/Afh3HgoiX0myG92j1avqPuD2JMd0295laozKkPuA3Tdrex+DozIfqaoHttSJ7mjM7wBvnWbxZxic0pmoqo01ePLpJINxKo43kRphOJEEQPcX9UuB53aXEl8D/DHwcWBVd4rmF4Hrhj52PfDrSa4F9gLeU1X3MjhacjXwCeDzQ+u/Evhf3SmiTwM/vFkZFwEHTw2I7drOAXZjbqd0pvpyelVdNk37PQzCyDVDzZ8BHg1cOdftS+pXBv8/kqStk2QFcG436LXP71kFvKOq/luf3yOpHUsWugBJmkmSNcCv8tAVO5LGgEdOJI2UJL/L4OZqwz5SVdONMZE0ggwnkiSpKQ6IlSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUlP8P322l3ET0YiAAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 648x432 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "group_plot(37)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "tags": [
     "remove_input"
    ]
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Price: $50\n",
      "Total Production/Market Supply:  3500\n",
      "\n",
      "Suppliers:  ['FOUR CORNERS' 'HUNTINGTON BEACH 1&2' 'REDONDO 5&6' 'REDONDO 7&8']\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAicAAAGECAYAAADtMIAAAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAfKUlEQVR4nO3dfbhmdV3v8fcHBtAAeRwnZDyAiiRZgo4G6kkBH8hUsDiIcdmkdDg92NE8nRpPJ7PUwk5pWqaRokOhgKRBJBghylUqOTwJCMiDDM/MqOCg+YR+zx9rbbgb996zZ5i19+/e9/t1Xfd1r/Vb61739zfrZubDbz2lqpAkSWrFNgtdgCRJ0ijDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJG01Sd6b5PcWuo6NJflgkrcsdB2S5sZwImnOktyS5FtJvpHk3iT/lOSxU8ur6leq6s0PY/svTHJxkvuTrE/y6SQv3TrVSxoXhhNJm+slVbUTsBdwD/AXW2OjSY4BPgKcCiwHlgFvBF6yNbYvaXwYTiRtkar6NnAWcOBU28aHT5L8dpK7ktyZ5JeTVJInbLytJAHeDry5qt5XVV+vqh9U1aer6r/362yT5P8mWZtkXZJTk+wyso2PJLk7ydf70ZcfH7L/koZjOJG0RZL8CPBy4HMzLD8SeD3wPOAJwHNn2dwBwGPpws5Mfql/HQY8DtgJ+MuR5ecB+wOPBi4DTttkJyQ1aclCFyBp7PxDkgeAHYH1wAtnWO9Y4ANVdQ1AkjcBx8+w7h79+12zfO/xwNur6uZ+e28Ark7yqqp6oKpOmVqx/657k+xSVV+fW7cktcKRE0mb6+iq2hV4BPAa4NNJfnSa9R4D3DYyf9s060z5av++1yzrPAZYOzK/lu5/sJYl2TbJSUluSrIBuKVfZ89ZtiepUYYTSVukqr5fVR8Fvg88e5pV7qI7sXXKY6dZZ8r1dOHl52dZ505gn5H5/wI8QHdS7i8AR9EdQtoF2LdfJ7NsT1KjDCeStkg6RwG7AddOs8qZwKuSPKk/P2XG+59UVdGdn/J7SV6V5FH9CbDPTnJyv9qHgd9Msl+SnYA/As6oqgeAnYHv0I3A/Ei/TNKYMpxI2lz/mOQbwAbgrcDKqfNKRlXVecC7gIuAG3noxNnvTLfRqjqL7gTbV9ONktwDvAU4u1/lFOBvgYuBLwPfBn6jX3Yq3WGeO4AvMsNJupLGQ7r/YZGkYSV5EnA1sEM/2iFJ03LkRNJgkrwsyQ5JdgPeBvyjwUTSphhOJA3pfwDrgJvoTpz91YUtR9I48LCOJElqiiMnkiSpKYYTSZLUlLG4ff2ee+5Z++6770KXIUmStoJLL730K1W1dKblYxFO9t13X9asWbPQZUiSpK0gydrZlntYR5IkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJDVu7+X7kGReX3sv32fB+rtkwb5ZkiTNyZ133MqRqy6f1+88/6SD5/X7RjlyIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNGSycJDkgyRUjrw1JXpdk9yQXJLmhf99tqBokSdL4GSycVNX1VXVQVR0EPA34D+BjwCrgwqraH7iwn5ckSQLm77DOEcBNVbUWOApY3bevBo6epxokSdIYmK9wchzw4X56WVXd1U/fDSybpxokSdIYGDycJNkeeCnwkY2XVVUBNcPnTkyyJsma9evXD1ylJElqxXyMnPwMcFlV3dPP35NkL4D+fd10H6qqk6tqRVWtWLp06TyUKUmSWjAf4eQVPHRIB+AcYGU/vRI4ex5qkCRJY2LQcJJkR+D5wEdHmk8Cnp/kBuB5/bwkSRIAS4bceFV9E9hjo7av0l29I0mS9EO8Q6wkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIknSQPZevg9JHvZr0gz6bB1JkibZnXfcypGrLn/Y2zn/pIO3QjXjw5ETSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJktQUw4kkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkpgwaTpLsmuSsJNcluTbJoUl2T3JBkhv6992GrEGSJI2XoUdO3gmcX1U/BjwFuBZYBVxYVfsDF/bzkiRJwIDhJMkuwE8D7weoqu9W1X3AUcDqfrXVwNFD1SBJksbPkCMn+wHrgQ8kuTzJ+5LsCCyrqrv6de4Glk334SQnJlmTZM369esHLFOSJLVkyHCyBHgq8J6qOhj4JhsdwqmqAmq6D1fVyVW1oqpWLF26dMAyJUlSS4YMJ7cDt1fVJf38WXRh5Z4kewH07+sGrEGSJI2ZwcJJVd0N3JbkgL7pCOCLwDnAyr5tJXD2UDVIkqTxs2Tg7f8GcFqS7YGbgVfRBaIzk5wArAWOHbgGSZI0RgYNJ1V1BbBimkVHDPm9kiRpfHmHWEmS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJktQUw4kkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUlCVDbjzJLcD9wPeBB6pqRZLdgTOAfYFbgGOr6t4h65AkSeNjPkZODquqg6pqRT+/CriwqvYHLuznJUmSAEhVDbfxbuRkRVV9ZaTteuC5VXVXkr2AT1XVAbNtZ0VSawarUpIkzafApSODFj9k6JGTAv45yaVJTuzbllXVXf303cCy6T6Y5MQka5KYSyRJmiCDnnMCPLuq7kjyaOCCJNeNLqyqSjLt0E1VnQycDLBixYpijRlFkjReknDkqssf9nbOP+ngrbKdzf3OwY6uJLMuHnTkpKru6N/XAR8DngHc0x/OoX9fN2QNkiRpvAwWTpLsmGTnqWngBcDVwDnAyn61lcDZQ9UgSZLGz5CHdZYBH0s3dLME+FBVnZ/k88CZSU4A1gLHDliDJEkaM4OFk6q6GXjKNO1fBY4Y6nslSdJ48w6xkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJktSUTYaTJM+aS5skSdLWMJeRk7+YY5skSdLDNuOD/5IcCjwTWJrk9SOLHgVsO3RhkiRpMs32VOLtgZ36dXYead8AHDNkUZIkaXLNGE6q6tPAp5N8sKrWAiTZBtipqjbMV4GSJGmyzOWckz9O8qgkOwJXA19M8r8HrkuSJE2ouYSTA/uRkqOB84D9gFcOWpUkSZpYcwkn2yXZji6cnFNV3wNq2LIkSdKkmks4+WvgFmBH4OIk+9CdFCtJkrTVzXa1DgBV9S7gXSNNa5McNlxJkiRpks3lDrG7JHl7kjX968/oRlEkSZK2urkc1jkFuB84tn9tAD4wZFGSJGlybfKwDvD4qvr5kfk/SHLFUAVJkqTJNpeRk28lefbUTP/Qv28NV5IkSZpkcxk5+RXg1CS79PP3Ar80WEWSJGmizeVqnSuBpyR5VD/vZcSSJGkwMx7WSfL6JCdMzVfVhqrakOSEJK+bn/IkSdKkme2ck+OBU6dp/1vg1cOUI0mSJt1sh3WW9Leq/0+q6rtJMmBNkjRv9l6+D3fecetClzGrJds/kge+63UImhyzhZNtkiyrqntGG5MsG7gmSZo3d95xK0euunyhy5jV+Scd3HyNmt75Jx280CWMpdkO6/w/4J+SPCfJzv3rucC5wJ/OS3WSJGnizDhyUlWnJlkP/CHwZLonEV8DvLGqzpun+iRJ0oSZ9VLiPoQYRCRJ0ryZyx1iH5Yk2ya5PMm5/fx+SS5JcmOSM5JsP3QNkiRpfAweToDXAteOzL8NeEdVPYHubrMnTPspSYvW3sv3IUkTL0ntmcvt67dYkuXAzwJvBV7fX4J8OPAL/SqrgTcB7xmyDkltaekKGa+mkNqzyZGTJMuSvD/Jef38gaN3jt2EPwd+G/hBP78HcF9VPdDP3w7sPcP3nphkTZI169evn+PXSZKkcTeXwzofBD4BPKaf/xKwydvXJ3kxsK6qLt2Swqrq5KpaUVUrli5duiWbkCRJY2gu4WTPqjqTfvSjH/X4/hw+9yzgpUluAU6nO5zzTmDXJFOHk5YDd2xu0ZIkafGaSzj5ZpI96O5zQpJDgK9v6kNV9YaqWl5V+wLHAZ+squOBi4Bj+tVWAmdvSeGSJGlxmssJsa8HzgEen+TfgKU8FC62xO8Apyd5C3A58P6HsS1JkrTIbDKcVNVlSZ4DHAAEuH66BwJuYhufAj7VT98MPGOzK5UkSRNhxnCS5OdmWPTEJFTVRweqSZIkTbDZRk5eMsuyAgwnkiRpq5vtwX+vms9CJEmSYG43YdsjybuSXJbk0iTv7K/ekSRJ2urmcinx6cB64OfprtJZD5wxZFGSJGlyzeVS4r2q6s0j829J8vKhCpIkSZNtLiMn/5zkuCTb9K9j6W5nL0mStNXNdinx/XRX5YTuWTp/1y/aBvgG8FuDVydJkibObFfr7DyfhUiSJMHczjkhyW7A/sAjptqq6uKhippPey/fhzvvuHWhy5AkSb1NhpMkvwy8lu4JwlcAhwCfpXvK8Ni7845bOXLV5QtdhjRRzj/p4IUuQVLD5nJC7GuBpwNrq+ow4GDgvkGrkiRJE2su4eTbVfVtgCQ7VNV1dA8BlCRJ2urmcs7J7Ul2Bf4BuCDJvcDaYcuSJEmTapPhpKpe1k++KclFwC7A+YNWJUmSJtZs9zl5VFVtSLL7SPNV/ftOwNcGrUySJE2k2UZOPgS8GLiUh27GNvr+uMGrkyRJE2e2m7C9OEmA51SVNwKRJEnzYtardaqqgH+ap1okSZLmdCnxZUmePnglkiRJzO1S4p8Cjk+yFvgm/TknVfWTg1YmSZIm0lzCyQsHr0KSJKk3l/ucrAVI8mhGHvwnSZI0hE2ec5LkpUluAL4MfBq4BThv4LokSdKEmssJsW+mexLxl6pqP+AI4HODViVJkibWXMLJ96rqq8A2SbapqouAFQPXJUmSJtRcToi9L8lOwMXAaUnW0V21I0mStNXNZeTkKOBbwG/SPfDvJuAlQxYlSZIm12wP/ns38KGq+reR5tXDlyRJkibZbCMnXwL+NMktSf4kycHzVZQkSZpcM4aTqnpnVR0KPAf4KnBKkuuS/H6SJ85bhZIkaaJs8pyTqlpbVW+rqoOBVwBHA9cOXpkkSZpIc7kJ25IkL0lyGt3N164Hfm7wyiRJ0kSa7YTY59ONlLwI+HfgdODEqprTZcRJHkF3+fEO/fecVVW/n2S/flt7AJcCr6yq7z6sXkiSpEVjtpGTNwCfAZ5UVS+tqg/NNZj0vgMcXlVPAQ4CjkxyCPA24B1V9QTgXuCELaxdkiQtQrOdEHt4Vb2vqu7dkg1X5xv97Hb9q4DDgbP69tV057BIkiQBc7sJ2xZLsm2SK4B1wAV0N3C7r6oe6Fe5Hdh7hs+emGRNkjXr168fskxJktSQQcNJVX2/qg4ClgPPAH5sMz57clWtqKoVS5cuHaxGSZLUlkHDyZSqug+4CDgU2DXJ1Im4y4E75qMGSZI0HgYLJ0mWJtm1n34k8Hy6+6NcBBzTr7YSOHuoGiRJ0viZy1OJt9RewOok29KFoDOr6twkXwROT/IW4HLg/QPWIEmSxsxg4aSqvgD80PN4qupmuvNPJEmSfsi8nHMiSZI0V4YTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJktQUw4kkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkpg4WTJI9NclGSLya5Jslr+/bdk1yQ5Ib+fbehapAkSeNnyJGTB4D/VVUHAocAv57kQGAVcGFV7Q9c2M9LkiQBA4aTqrqrqi7rp+8HrgX2Bo4CVverrQaOHqoGSZI0fublnJMk+wIHA5cAy6rqrn7R3cCy+ahBkiSNh8HDSZKdgL8HXldVG0aXVVUBNcPnTkyyJsma9evXD12mJElqxKDhJMl2dMHktKr6aN98T5K9+uV7Aeum+2xVnVxVK6pqxdKlS4csU5IkNWTIq3UCvB+4tqrePrLoHGBlP70SOHuoGiRJ0vhZMuC2nwW8ErgqyRV92/8BTgLOTHICsBY4dsAaJEnSmBksnFTVvwKZYfERQ32vJEkab94hVpIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJktQUw4kkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1ZbBwkuSUJOuSXD3StnuSC5Lc0L/vNtT3S5Kk8TTkyMkHgSM3alsFXFhV+wMX9vOSJEkPGiycVNXFwNc2aj4KWN1PrwaOHur7JUnSeJrvc06WVdVd/fTdwLJ5/n5JktS4BTshtqoKqJmWJzkxyZoka9avXz+PlUmSpIU03+HkniR7AfTv62ZasapOrqoVVbVi6dKl81agJElaWPMdTs4BVvbTK4Gz5/n7JUlS44a8lPjDwGeBA5LcnuQE4CTg+UluAJ7Xz0uSJD1oyVAbrqpXzLDoiKG+U5IkjT/vECtJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSZLUFMOJJElqiuFEkiQ1xXAiSZKaYjiRJElNMZxIkqSmGE4kSVJTDCeSJKkphhNJktQUw4kkSWqK4USSJDXFcCJJkppiOJEkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1BTDiSRJaorhRJIkNcVwIkmSmrIg4STJkUmuT3JjklULUYMkSWrTvIeTJNsC7wZ+BjgQeEWSA+e7DkmS1KaFGDl5BnBjVd1cVd8FTgeOWoA6JElSgxYinOwN3DYyf3vfJkmSRKpqfr8wOQY4sqp+uZ9/JfBTVfWajdY7ETixnz0AuH6gkvYEvjLQtlsyKf2EyenrpPQTJqevk9JPmJy+Tko/YfP6uk9VLZ1p4ZKtU89muQN47Mj88r7tP6mqk4GThy4myZqqWjH09yy0SeknTE5fJ6WfMDl9nZR+wuT0dVL6CVu3rwtxWOfzwP5J9kuyPXAccM4C1CFJkho07yMnVfVAktcAnwC2BU6pqmvmuw5JktSmhTisQ1V9HPj4Qnz3NAY/dNSISeknTE5fJ6WfMDl9nZR+wuT0dVL6CVuxr/N+QqwkSdJsvH29JElqysSGk8V4C/0ktyS5KskVSdb0bbsnuSDJDf37bn17kryr7/8Xkjx1YaufWZJTkqxLcvVI22b3K8nKfv0bkqxciL5sygx9fVOSO/r9ekWSF40se0Pf1+uTvHCkvenfd5LHJrkoyReTXJPktX37otqvs/RzMe7TRyT59yRX9n39g759vySX9HWf0V8IQZId+vkb++X7jmxr2j+DVszS1w8m+fLIfj2obx/L3++UJNsmuTzJuf388Pu0qibuRXci7k3A44DtgSuBAxe6rq3Qr1uAPTdq+xNgVT+9CnhbP/0i4DwgwCHAJQtd/yz9+mngqcDVW9ovYHfg5v59t356t4Xu2xz7+ibgt6ZZ98D+t7sDsF//m952HH7fwF7AU/vpnYEv9f1ZVPt1ln4uxn0aYKd+ejvgkn5fnQkc17e/F/jVfvrXgPf208cBZ8z2Z7DQ/ZtjXz8IHDPN+mP5+x2p//XAh4Bz+/nB9+mkjpxM0i30jwJW99OrgaNH2k+tzueAXZPstRAFbkpVXQx8baPmze3XC4ELquprVXUvcAFw5PDVb54Z+jqTo4DTq+o7VfVl4Ea633bzv++ququqLuun7weupbtT9KLar7P0cybjvE+rqr7Rz27Xvwo4HDirb994n07t67OAI5KEmf8MmjFLX2cylr9fgCTLgZ8F3tfPh3nYp5MaThbrLfQL+Ockl6a7wy7Asqq6q5++G1jWT4/7n8Hm9mvc+/uafjj4lKlDHSySvvZDvwfT/d/not2vG/UTFuE+7Yf/rwDW0f1DexNwX1U90K8yWveDfeqXfx3YgzHta1VN7de39vv1HUl26NvGeb/+OfDbwA/6+T2Yh306qeFksXp2VT2V7onPv57kp0cXVje+tuguz1qs/RrxHuDxwEHAXcCfLWw5W0+SnYC/B15XVRtGly2m/TpNPxflPq2q71fVQXR3/n4G8GMLXNJgNu5rkicDb6Dr89PpDtX8zgKW+LAleTGwrqoune/vntRwMqdb6I+bqrqjf18HfIzuL4d7pg7X9O/r+tXH/c9gc/s1tv2tqnv6vwh/APwNDw2HjnVfk2xH9w/2aVX10b550e3X6fq5WPfplKq6D7gIOJTuEMbUPbVG636wT/3yXYCvMr59PbI/jFdV9R3gA4z/fn0W8NIkt9AdSjwceCfzsE8nNZwsulvoJ9kxyc5T08ALgKvp+jV1BvhK4Ox++hzgF/uzyA8Bvj4ynD4ONrdfnwBekGS3fgj9BX1b8zY6F+hldPsVur4e158hvx+wP/DvjMHvuz8O/X7g2qp6+8iiRbVfZ+rnIt2nS5Ps2k8/Eng+3Tk2FwHH9KttvE+n9vUxwCf70bKZ/gyaMUNfrxsJ1qE7D2N0v47d77eq3lBVy6tqX7rf3Cer6njmY5/OdrbsYn7RnT39Jbpjor+70PVshf48ju5s6CuBa6b6RHe870LgBuBfgN379gDv7vt/FbBiofswS98+TDf0/T26Y5UnbEm/gFfTnYh1I/Cqhe7XZvT1b/u+fKH/j3yvkfV/t+/r9cDPjLQ3/fsGnk13yOYLwBX960WLbb/O0s/FuE9/Eri879PVwBv79sfR/UN0I/ARYIe+/RH9/I398sdt6s+gldcsff1kv1+vBv6Oh67oGcvf70Z9fi4PXa0z+D71DrGSJKkpk3pYR5IkNcpwIkmSmmI4kSRJTTGcSJKkphhOJElSUwwnkiSpKYYTSQ9K8qNJTk9yU/+Mpo8neeLA3/mYJGf10wcledEWbudNSSrJE0baXte3rUjy2iR/PrLsr5P8y8j8byR518Ppi6Stw3AiCXjwrpYfAz5VVY+vqqfRPStk2eyffHiq6s6qmrrb5EF0NxvbUlfR3clyyn+juykhwL8BzxxZ9hRglyTb9vPPBD7zML5b0lZiOJE05TDge1X13qmGqroSuDzJhUkuS3JVkqOge8pukuuSnJbk2iRnJfmRftkbk3w+ydVJTu6DD0mekORfklzZb+/x/Xau7m/L/ofAy5NckeTlSW5IsrT/7DZJbpyan8E/0D2enSSPp3sq6lf6ZVcAT0zyyCS7AN/q236iX/5MugAjaYEZTiRNeTIw3dNHvw28rLonXh8G/NlU2AAOAP6qqp4EbAB+rW//y6p6elU9GXgk8OK+/TTg3VX1FLow8ODznKrqu8AbgTOq6qCqOoPuFuDH96s8D7iyqtbP0ocNwG3pnhB7HHDGyPYfoLvl+NOBQ4BLgM8Bz0yyN5Cquu2HNylpvhlOJG1KgD9K8gW6593szUOHem6rqqnRhr+je5YMwGFJLklyFd2TTH+8fzDl3lX1MYCq+nZV/ccmvvsU4Bf76VfTPel1U06nCyZH0x2mGvUZulD0TOCz/Wtq3kM6UiMMJ5KmXAM8bZr244GlwNOq6iDgHroHfEH3ULtRleQRwF8Bx1TVTwB/M7L+ZulHMu5Jcjjd4+fPm8PHzgVeCdxaVRs2WjZ13smhdMHkWuBADCdSUwwnkqZ8EtghyYlTDUl+EtgHWFdV30tyWD8/5b8kObSf/gXgX3koiHwlyU70j1avqvuB25Mc3W97h6lzVEbcD+y8Udv76EZlPlJV399UJ/rRmN8B3jrN4s/SHdJZWlXrqnvy6Xq681Q830RqhOFEEgD9P9QvA57XX0p8DfDHwMeBFf0hml8Erhv52PXArye5FtgNeE9V3Uc3WnI18Ang8yPrvxL4n/0hos8AP7pRGRcBB06dENu3nQPsxNwO6Uz15fSqumya9nvpwsg1I82fBR4NXDnX7UsaVrq/jyRp8yTZFzi3P+l1yO9ZAbyjqv7rkN8jqR1LFroASZpJklXAr/LQFTuSJoAjJ5LGSpLfpbu52qiPVNV055hIGkOGE0mS1BRPiJUkSU0xnEiSpKYYTiRJUlMMJ5IkqSmGE0mS1JT/D60oln8R0K8XAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 648x432 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "group_plot(50)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now we will look at all the energy sources. They have been colored according to source for reference."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "tags": [
     "remove_input"
    ]
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Price: $30\n",
      "Total Production/Market Supply:  2800\n",
      "\n",
      "Suppliers:  ['BIG CREEK' 'HELMS' 'DIABLO CANYON 1']\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiYAAAGECAYAAAAC8us+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nO3df5xmdV338dcbll8CCygbwrA/yNQizF+rgliGmhqJYBFRhGvRg82iJK3U6lHe9tP7TkNv1HYLdU0UV/wFpika6h0CuiAIiCj+WHYZhCXFJRUF/dx/XGfwcpyZnZ2Za67v7Lyej8d5zHXOuc45n2vOzOx7v+d7vidVhSRJUgt2G3YBkiRJYwwmkiSpGQYTSZLUDIOJJElqhsFEkiQ1w2AiSZKaYTCRdjFJ3pTkb7rXP59k67BrkqTpMphIC1SSjyb5epK9ZrGPSvLNJP/TN/3pXNY5CEmelOQTSb6R5GtJLkvyuGHXJWn2lgy7AEk7L8kq4GeBbwDPBt4xi909sqpunoOyJpVkSVXdN0f7Wgq8D3g+sBHYk9734jtzsf9xx5qzuiVNjy0m0sL0XOAK4E3AmkEcIMnLkmxM8uYkdye5IcnqvvWHJXlnkm1JvpzkD8dte2GStyTZDjwvyT5JNnStPDcm+dOxy0xJ/iTJO8cd/zVJXj1BaQ8DqKq3VdX3qurbVfWhqvpMt91uSf4iyeYkd3T1H9Ct+5FLW0m+kuRpU9S9e5I/S/LF7vtwVZLl3ft/MsklXavNTUlO6dvv8Uk+221za5I/ntUJkRYJg4m0MD0XOL+bnpHkkAEd59nABcCBwEXAudD7xx+4GLgWGAGeCpyd5Bl9254IXNhtez7wV8Aq4MeBXwB+s++9bwGemeTAbv9LgFOBN09Q0+eB73Uh5xeTHDRu/fO66bjuWPuN1T1N4+t+IfDrwPHAUuC3gW8l2Re4BHgr8GNdva9LcmS3n/OAtVW1P3AU8J87UYO0aBlMpAUmyZOAlcDGqroK+CLwG7PY5dVJ7uqb+sPFf1XV+6vqe8C/AY/slj8OWFZVL6+q71bVl4B/ofeP85jLq+o9VfX9qvo2cArwd1X19araCrxm7I1VdRvwceBXu0XPBO7sPt8PqartwJOA6o65LclFfeHsNOBVVfWlqvof4KXAqV3YmY7xdf8O8BdVdVP1XFtV/w08C/hKVb2xqu6rqk8D7+z7DPcCRyZZ2n3mq6d5fGlRM5hIC88a4ENVdWc3/1ZmdznnMVV1YN/0wb51X+17/S1g7+4f+JXAYf2BBvgzoL/lZsu44xw2btn49Rv4QSvKb9ILQhOqqhur6nlVdTi91ojDgHP6jrO57+2b6fWnm26r0vi6ltMLf+OtBJ4w7ntwGvDgbv2v0Gtl2ZzkY0mOmebxpUXNzq/SApJkH3otD7snGQsNewEHJnlkVV07T6VsAb5cVQ+d4j3jH11+G3A48Nlufvm49e8BXp/kKHqtEdO6O6iqPpfkTcDabtEovdAwZgVwH3A7vdDygLEVSXYHlu2g7i3AQ4DrJ1j+sar6hUnq+hRwYpI9gLPoddQd/5kljWOLibSwnAR8DzgSeFQ3/RTw/+j1O5kvnwTuTvLirlPr7kmO2sEtuxuBlyY5KMkIvX+s71dV99Dr2/FW4JNVdctEO+k6nL4oyeHd/HJ6fUCu6N7yNuCPkhyRZD/g74C3d3fXfJ5eq88vdYHhL+gFu6n8K/DXSR6anp9J8iB6dwY9LMnpSfbopscl+akkeyY5LckBVXUvsB34/g6OIwmDibTQrAHeWFW3VNVXxyZ6nTtP24l+FP2uzQ+PY3LOjjbo+pw8i14w+jJwJ71/wA+YYrOXA1u793+YXggZf4vvBuARTHEZB7gbeAJwZZJv0gsk1wMv6ta/odv+492x7gH+oKv7G8DvdbXeCnyzq2kqr6IXqj5EL2CcB+xTVXcDT6fXr2aU3mWvV/CDoHM68JXu7p7fpXeZR9IOpGp8q6UkDV6S5wOnVtWT+5atAD4HPLjr5CppkbHFRNK8SHJokmO7cUYeTq+F491963ejd2vuBYYSafGy86uk+bInsA44AriL3vgorwPoxgS5nd4dNM8cVoGShs9LOZIkqRleypEkSc0wmEiSpGYsiD4mBx98cK1atWrYZUiSpDlw1VVX3VlV4wc3BBZIMFm1ahWbNm0adhmSJGkOJNk82Tov5UiSpGYYTCRJUjMMJpIkqRkGE0mS1AyDiSRJaobBRJIkNcNgIkmSmmEwkSRJzTCYSJKkZhhMJElSMwwmkiSpGQYTSZLUDIOJJElqhsFEkiQxsnwlSe6fRpavHEodS4ZyVEmS1JTRrbdw5rpN98+vX7t6KHXYYiJJkpphMJEkSc0wmEiSpGYYTCRJUjMMJpIkqRkGE0mS1AyDiSRJaobBRJIkNcNgIkmSmmEwkSRJzTCYSJKkZhhMJElSMwwmkiSpGQYTSZLUDIOJJElqhsFEkiQ1w2AiSZKaYTCRJEnNMJhIkqRmGEwkSVIzDCaSJKkZBhNJktQMg4kkSWqGwUSSJDXDYCJJkpphMJEkSc0wmEiSpGYYTCRJUjMMJpIkqRkDDSZJ/ijJDUmuT/K2JHsnOSLJlUluTvL2JHsOsgZJkrRwDCyYJBkB/hBYXVVHAbsDpwKvAP6pqn4C+DpwxqBqkCRJC8ugL+UsAfZJsgR4AHAb8BTgwm79BuCkAdcgSZIWiIEFk6q6FfhH4BZ6geQbwFXAXVV1X/e2rcDIoGqQJEkLyyAv5RwEnAgcARwG7As8cye2PzPJpiSbtm3bNqAqJUlSSwZ5KedpwJeraltV3Qu8CzgWOLC7tANwOHDrRBtX1fqqWl1Vq5ctWzbAMiVJUisGGUxuAY5O8oAkAZ4KfBa4FDi5e88a4L0DrEGSJC0gg+xjciW9Tq5XA9d1x1oPvBh4YZKbgQcB5w2qBkmStLAs2fFbZq6q/gr4q3GLvwQ8fpDHlSRJC5Mjv0qSpGYYTCRJUjMMJpIkqRkGE0mS1AyDiSRJaobBRJIkNcNgIkmSmmEwkSRJzTCYSJKkZhhMJElaZEZWrCDJD02tGOiQ9JIkaXhWjYyweXR0wnVrLr74h+Y3nHDCfJS0QwYTSZJ2UZtHR7lzzdofWX7whnVDqGZ6vJQjSZKaYTCRJEnNMJhIkqRmGEwkSVIzDCaSJKkZBhNJkha4VctHfmRckpbGJtkZ3i4sSdICt3nrKNvPOeVHli89e+MQqpkdW0wkSVIzDCaSJKkZBhNJktQMg4kkSWqGwUSSJDXDYCJJkpphMJEkSc0wmEiSpGYYTCRJUjMMJpIkqRkGE0mS1AyDiSRJC9CK5SsW9MP6JmMwkSRpAdqydQuXnXsZl5172bBLmVMGE0mS1AyDiSRJaobBRJIkNcNgIkmSmmEwkSRJzTCYSJKkZhhMJElSMwwmkiSpGQYTSZLUDIOJJElqhsFEkqQFYtXykV3y+Tj9DCaSJC0Qm7eOsv2cU9h+zinDLmVgDCaSJKkZBhNJkho2snzlLn/5pp/BRJKkho1uvYUz123izHWbhl3KvDCYSJKkZhhMJElSMwwmkiSpGQYTSZLUDIOJJElqhsFEkiQ1w2AiSZKaYTCRJKkxq0Z2/WfiTMZgIklSYzaPjnLnmrXcuWbtsEuZdwYTSZLUDIOJJElqhsFEkiQ1w2AiSZKaMdBgkuTAJBcm+VySG5Mck+SBSS5J8oXu60GDrEGSpIVg1crFeydOv0G3mLwa+I+q+kngkcCNwEuAj1TVQ4GPdPOSJC1qm28Z5TtXnMF3rjhj2KUM1cCCSZIDgJ8DzgOoqu9W1V3AicCG7m0bgJMGVYMkSVpYBtlicgSwDXhjkk8n+dck+wKHVNVt3Xu+ChwywBokSdICMshgsgR4DPD6qno08E3GXbapqgJqoo2TnJlkU5JN27ZtG2CZkiSpFYMMJluBrVV1ZTd/Ib2gcnuSQwG6r3dMtHFVra+q1VW1etmyZQMsU5KkubFi+Yr7O7Du7KSeJYPacVV9NcmWJA+vqpuApwKf7aY1wD90X987qBokSZpPW7Zu4bJzL5vRtseedewcV7MwDSyYdP4AOD/JnsCXgN+i10qzMckZwGbglAHXIEmSFoiBBpOqugZYPcGqpw7yuJIkaWFy5FdJktQMg4kkSbPQ3+FVs2cwkSRpFsY6vM6006t+mMFEkiQ1w2AiSZKaYTCRJEnNMJhIkqRm7DCYJPmRoegmWiZJkjRb02kx+b/TXCZJkjQrk478muQY4InAsiQv7Fu1FNh90IVJkqTFZ6oh6fcE9uves3/f8u3AyYMsSpIkLU6TBpOq+hjwsSRvqqrNAEl2A/arqu3zVaAkSVo8ptPH5O+TLE2yL3A98NkkfzLguiRJ0iI0nWByZNdCchLwAeAI4PSBViVJkhal6QSTPZLsQS+YXFRV9wI12LIkSdJiNJ1gsg74CrAv8PEkK+l1gJUkSZpTOwwmVfWaqhqpquOrZzNw3DzUJklSs5avWE6SYZexy5nOyK8HJHlVkk3d9Ep6rSeSJC1aW7dsZd3l5w67jF3OdC7lvAG4Gzilm7YDbxxkUZIkaXGaaoC1MQ+pql/pm/9fSa4ZVEGSJGnxmk6LybeTPGlspnuA37cHV5IkSVqsptNi8rvAm5Mc0M1/HXjewCqSJGlIRpavZHTrLcMuY1HbYTCpqmuBRyZZ2s17q7AkaZexamSEzaOj98+fuW7TtLZbv3b1oEpa1Ca9lJPkhUnOGJuvqu1VtT3JGUnOnp/yJEkarM2jo9y5Zi13rlk77FLE1H1MTgPePMHyfwN+ezDlSJKkxWyqYLKkG37+h1TVdwFHlJEkSXNuqmCyW5JDxi+caJkkSdJcmCqY/B/g35M8Ocn+3fTzwPuAf5yX6iRJ0qIy6V05VfXmJNuAlwNH0Xui8A3AX1bVB+apPkmStIhMebtwF0AMIZIkaV5MZ+RXSZKkeWEwkSQtOstXLCcJiTeZtmY6Q9JLkrTgjaxYweiWLffPr7v8XADWHnPWsErSBHbYYpLkkCTnJflAN39k/4iwkiQtBKNbtrDm4otZc/HFwy5FU5jOpZw3AR8EDuvmPw84JL0kSZpz0wkmB1fVRuD7AFV1H/C9gVYlSZIWpekEk28meRC9cUxIcjTwjYFWJUmSFqXpdH59IXAR8JAklwHLgJMHWpUkSVqUdhhMqurqJE8GHk7v4X03TfRwP0mSWjL+LhwtDJMGkyS/PMmqhyWhqt41oJokSZq1sbtwxmw44YQhVqPpmqrFZKozWIDBRJIkzampHuL3W/NZiCRJ0nQGWHtQktckuTrJVUle3d2lI0mSNKemc7vwBcA24Ffo3Y2zDXj7IIuSJEmL03RuFz60qv66b/5vkvzaoAqSJEmL13RaTD6U5NQku3XTKfSGqJckSZpTU90ufDe9u29C79k4b+lW7Qb8D/DHA69OkiQtKlPdlbP/fBYiSdLOGFm+ktGttwy7DM2x6fQxIclBwEOBvceWVdXHB1WUJEk7Mrr1Fs5ct2nS9evXrp7HajRXdhhMkvwO8ALgcOAa4GjgcuApgy1NkiQtNtPp/PoC4HHA5qo6Dng0cNdAq5IkSYvSdILJPVV1D0CSvarqc/Qe6CdJ0qyNLF9Jkp2etGuaTh+TrUkOBN4DXJLk68DmwZYlSVosdtRXZDL2Idk17TCYVNVzupcvS3IpcADwHwOtSpIkLUpTjWOytKq2J3lg3+Lruq/7AV8baGWSJGnRmarF5K3As4Cr+MFAa/1ff3zg1UmSmrVqZITNo6PDLkO7mKkGWHtWer2LnlxVjmAjSbuwmYaMO9esnfWxD96wbtb70K5jyj4mVVVJ/h14xDzVI0nqrFo+wuatE4eFB+y1B9/6zr1zerydDRkGCg3CdO7KuTrJ46rqUwOvRpIWsEEMkb79nFMmXL707I2TrpuJpWdvnLN9SbMxnWDyBOC0JJuBb9L1MamqnxloZZK0wMz0ttfJeDusFqPpBJNnzOYASXYHNgG3dv1WjgAuAB5Er2Pt6VX13dkcQ5Lmmg+Ik4ZjOuOYbAZI8mP0PcRvJ7wAuBFY2s2/AvinqrogyT8DZwCvn8F+JWmHVixfwZatW2a07c62ftjCIc3edB7i92zglcBhwB3ASnpB46ense3hwC8Bfwu8sLvL5ynAb3Rv2QC8DIOJpFmaqqPoZedettP7O/asY2dbkqQZmM6lnL+m90ThD1fVo5McB/zmNPd/DvCnwP7d/IOAu6rqvm5+KzCyE/VK0oQ2bx2dsDOonTqlhWU6D/G7t6r+G9gtyW5VdSmww/bKJM8C7qiqq2ZSWJIzk2xKsmnbtm0z2YUkSVpgptNicleS/YCPA+cnuYPe3Tk7cizw7CTH0+ubshR4NXBgkiVdq8nhwK0TbVxV64H1AKtXr65pHE+SJC1w02kxORH4NvBH9B7e90XghB1tVFUvrarDq2oVcCrwn1V1GnApcHL3tjXAe2dQt6Rd0IrlK3bqsff9k6Rdw1QP8Xst8Naq6u81tmEOjvli4IIkfwN8GjhvDvYpaQGZavjzmXRUBTurSruKqS7lfB74xySHAhuBt1XVp2dykKr6KPDR7vWXgMfPZD+SFo5VK0fYfMvkz16ZaPhzhziXNNVD/F4NvDrJSnqXYt6QZB/gbfRCyufnqUZJC9DmW0b5zhVnTLhur6NtKJU0sR32MamqzVX1iqp6NPDrwEn0xjGRJEmaUzsMJkmWJDkhyfnAB4CbgF8eeGWS5syqkZEZdyrdd8897Ywqad5M1fn1F+i1kBwPfJLe823OrKrp3CosqSGbR0d3+pH2Yw7esG5G29pfRNJMTNX59aXAW4EXVdXX56keSZK0iE3V+fUp81mIJEnSdAZYkyRJmhfTGZJeWrCmeuLssD1g7z341j337vR2e++1N/d8554BVCRJw2cw0S5tsifOtmDp2RsnHedjKnsdfd5Oj47qqKiSFgov5UiSpGYYTCRJUjMMJpIkqRkGE0mS1AyDiSRJaobBRAO3fMXyGT+nZe999p7xtj6vRZIWHm8X1sBt3bKVdZefO6Nt1x5z1oy3HdtekrRw2GIiSZKaYTCRJEnNMJhIkqRmGEwkSVIzDCaSJKkZBpNFYtXKkRnfcrvvPnt6y64kaV54u/AisfmW0Rk9yRZ6T7Od6bZj20uSNB22mEiSpGYYTCRJUjMMJpIkqRkGE0mS1AyDiSRJakaqatg17NDqpDYNuwhJkjQnAldV1eqJ1tliIkmSmrEwxjF57GNhk20ms5FkqOOYrLv83Bltu/aYs2a87dj22885ZcbbD9LSszfO6Pu619Hncdm5l+3UNseedSx3rlm708cCOHjDuhlte/CGdZN+vr2OPm/CfR68Yd1Of7Yxx5517ITneunZG2e0z2PPOpYz1+3c3531a1fv9DY72t9kP79Lz944pz/bS8/euNPneaY/GxPtZybftx19v9evXc2aiy++f37DCSfcP7/hhBPu/9uy9piz7v8cO1PL+rWrWXf5uaw95qz7f8aOPevYWf0Mj/3O9P+O9NfU/zPR/7M91c//ZL9r/d8b6H1P+j/7+rWrGdhVlSkG37TFRJIkNcNgIkmSmmEwkSRJzTCYSJKkZhhMJElSMwwmkiSpGQYTSZLUDIOJJElqhsFEkiQ1w2AiSZKaYTCRJEnNMJhIkqRmGEwkSVIzDCaSJKkZBhNJktQMg4kkSWqGwUSSJDXDYCJJkpphMJEkSc0wmEiSpGYYTCRJUjMMJpIkqRkGE0mS1AyDiSRJaobBRJIkNcNgIkmSmmEwkSRJzTCYSJKkZhhMJElSMwYWTJIsT3Jpks8muSHJC7rlD0xySZIvdF8PGlQNkiRpYRlki8l9wIuq6kjgaOD3kxwJvAT4SFU9FPhINy9JkjS4YFJVt1XV1d3ru4EbgRHgRGBD97YNwEmDqkGSJC0s89LHJMkq4NHAlcAhVXVbt+qrwCHzUYMkSWrfwINJkv2AdwJnV9X2/nVVVUBNst2ZSTYl2bRt27ZBlylJkhow0GCSZA96oeT8qnpXt/j2JId26w8F7pho26paX1Wrq2r1smXLBlmmJElqxCDvyglwHnBjVb2qb9VFwJru9RrgvYOqQZIkLSxLBrjvY4HTgeuSXNMt+zPgH4CNSc4ANgOnDLAGSZK0gAwsmFTVfwGZZPVTB3VcSZK0cDnyqyRJaobBRJIkNcNgIkmSmmEwkSRJzTCYSJKkZhhMJElSMwwmkiSpGQYTSZLUDIOJJElqhsFEkiQ1w2AiSZKaYTCRJEnNMJhIkqRmGEwkSVIzDCaSJKkZBhNJktQMg4kkSWqGwUSSJDXDYCJJkpphMJEkSc0wmEiSpGYYTCRJUjMMJpIkqRkGE0mS1AyDSYNGVqwgyZxOkiQtBEuGXYB+1OiWLay5+OI53eeGE06Y0/1JkjQItphIkqRmGEwkSVIzDCaSJKkZBhNJktQMg4kkSWqGwUSSJDXDYCJJkpphMJEkSc0wmEiSpGYYTCRJUjMMJpIkqRkGE0mS1AyDiSRJaobBRJIkNcNgIkmSmmEwkSRJzTCYSJKkZhhMJElSMwwmkiSpGQYTSZLUDIOJJElqhsFEkiQ1w2AiSZKaYTCRJEnNMJhIkqRmGEwkSVIzDCaSJKkZBhNJktQMg4kkSWqGwUSSJDXDYCJJkpphMJEkSc0wmEiSpGYYTCRJUjOGEkySPDPJTUluTvKSYdQgSZLaM+/BJMnuwGuBXwSOBH49yZHzXYckSWrPMFpMHg/cXFVfqqrvAhcAJw6hjoEYWbGCJLOaJElarJYM4ZgjwJa++a3AE4ZQx0CMbtnCmosvntU+NpxwwhxVI0nSwpKqmt8DJicDz6yq3+nmTweeUFVnjXvfmcCZ3ezDgZsGVNLBwJ0D2rdmx3PTNs9Puzw3bfP8wMqqWjbRimG0mNwKLO+bP7xb9kOqaj2wftDFJNlUVasHfRztPM9N2zw/7fLctM3zM7Vh9DH5FPDQJEck2RM4FbhoCHVIkqTGzHuLSVXdl+Qs4IPA7sAbquqG+a5DkiS1ZxiXcqiq9wPvH8axJzDwy0WaMc9N2zw/7fLctM3zM4V57/wqSZI0GYeklyRJzVi0wcRh8YcjyVeSXJfkmiSbumUPTHJJki90Xw/qlifJa7pz9Jkkj+nbz5ru/V9IsmZYn2ehS/KGJHckub5v2ZydjySP7c73zd22jiC4EyY5Py9Lcmv3O3RNkuP71r20+17flOQZfcsn/HvX3YRwZbf87d0NCZqGJMuTXJrks0luSPKCbrm/P7NVVYtuotfp9ovAjwN7AtcCRw67rsUwAV8BDh637H8DL+levwR4Rff6eOADQICjgSu75Q8EvtR9Pah7fdCwP9tCnICfAx4DXD+I8wF8sntvum1/cdifeSFNk5yflwF/PMF7j+z+lu0FHNH9jdt9qr93wEbg1O71PwPPH/ZnXigTcCjwmO71/sDnu3Pg788sp8XaYrJLD4u/AJ0IbOhebwBO6lv+5uq5AjgwyaHAM4BLquprVfV14BLgmfNd9K6gqj4OfG3c4jk5H926pVV1RfX+yr65b1+ahknOz2ROBC6oqu9U1ZeBm+n9rZvw7133v++nABd22/efa+1AVd1WVVd3r+8GbqQ3srm/P7O0WIPJRMPijwyplsWmgA8luaob3RfgkKq6rXv9VeCQ7vVk58nzN1hzdT5Gutfjl2v2zuouB7xh7FIBO39+HgTcVVX3jVuunZRkFfBo4Er8/Zm1xRpMNDxPqqrH0Hu69O8n+bn+ld3/DLxVrBGejya9HngI8CjgNuCVwy1ncUuyH/BO4Oyq2t6/zt+fmVmswWRaw+Jr7lXVrd3XO4B302tmvr1rtqT7ekf39snOk+dvsObqfNzavR6/XLNQVbdX1feq6vvAv9D7HYKdPz//Te9ywpJxyzVNSfagF0rOr6p3dYv9/ZmlxRpMHBZ/CJLsm2T/sdfA04Hr6X3vx3qirwHe272+CHhu15v9aOAbXRPpB4GnJzmoa8Z+erdMc2NOzke3bnuSo7v+DM/t25dmaOwfvc5z6P0OQe/8nJpkryRHAA+l13lywr933f/mLwVO7rbvP9fage5n+jzgxqp6Vd8qf39ma9i9b4c10esh/Xl6vdX/fNj1LIaJ3l0B13bTDWPfd3rXuj8CfAH4MPDAbnmA13bn6Dpgdd++fpte576bgd8a9mdbqBPwNnqXA+6ldw37jLk8H8Bqev9wfhE4l25QR6dZnZ9/677/n6H3j92hfe//8+57fRN9d3BM9veu+538ZHfe3gHsNezPvFAm4En0LtN8Brimm47392f2kyO/SpKkZizWSzmSJKlBBhNJktQMg4kkSWqGwUSSJDXDYCJJkpphMJEkSc0wmEi6X5IHJ7kgyRe75xm9P8nDBnzMw5Jc2L1+VJLjZ7iflyWpJD/Rt+zsbtnqJC9Ick7funVJPtw3/wdJXjObzyJp9gwmkoD7R7J8N/DRqnpIVT0WeCk/eAjZQFTVaFWNjT76KHqDVM3UdfRGNh3zq/QG8wO4DHhi37pHAgck2b2bfyLwiVkcW9IcMJhIGnMccG9V/fPYgqq6Fvh0ko8kuTrJdUlOhN4TVZN8Lsn5SW5McmGSB3Tr/jLJp5Jcn2R9F3pI8hNJPpzk2m5/D+n2c303XPrLgV9Lck2SX0vyhSTLum13S3Lz2Pwk3kPv8fIkeQjwDeDObt01wMOS7JPkAODb3bJHdOufSC+8SBoig4mkMUcBV02w/B7gOdV7KvRxwCvHggbwcOB1VfVTwHbg97rl51bV46rqKGAf4Fnd8vOB11bVI+kFgbHHw1NV3wX+Enh7VT2qqt4OvAU4rXvL04Brq2rbFJ9hO7AlyVH0Wk7e3rf/+4BPA48Djqb3iPorgCcmGaE33PeWH92lpPlkMJG0IwH+Lsln6D37Y4QfXN7ZUlVjrQxvoff8EIDjklyZ5DrgKcBPdw9wHKmqdwNU1T1V9a0dHPsN9B5eBr3nibxxGvVeQC+UnETv0lS/T9ALRE8ELu+msXkv40gNMJhIGnMD8NgJlp8GLAMeW1WPAm4H9u7WjX/YViXZG3gdcHJVPQL4l77375SuBeP2JE8BHg98YBqbvX77MToAAAFdSURBVA84HbilqraPWzfWz+QYeqHkRuBIDCZSMwwmksb8J7BXkjPHFiT5GWAlcEdV3ZvkuG5+zIokx3SvfwP4L34QQu5Msh9wMkBV3Q1sTXJSt++9xvqk9Lkb2H/csn+l1xrzjqr63o4+RNcK82LgbydYfTm9yzjLquqO6j3FdBu9fin2L5EaYDCRBED3j/RzgKd1twvfAPw98H5gdXdZ5rnA5/o2uwn4/SQ3AgcBr6+qu+i1klwPfBD4VN/7Twf+sLss9AngwePKuBQ4cqzza7fsImA/pncZZ+yzXFBVV0+w/Ov0gsgNfYsvB34MuHa6+5c0OOn9LZKknZNkFfC+roPrII+zGvinqvrZQR5HUhuWDLsASZpMkpcAz+cHd+ZI2sXZYiJpQUny5/QGTuv3jqqaqE+JpAXGYCJJkpph51dJktQMg4kkSWqGwUSSJDXDYCJJkpphMJEkSc34/7/d7KiOBIZXAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 648x432 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAATIAAACRCAYAAABAOTQxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAeWElEQVR4nO2dd7xV1ZXHvz9AmthRQfIAC2pUrGCPJmINUYktaogaUFCDDrElRhML9uho7GUsGCO2MdaomSQ6GnVQLDO2GI0hYiFqjAVULKz5Y+0jhxd4XIRzzj3P9f187ufde+697/zuKWvvvfZaa8vMCIIgqDMdqhYQBEGwoIQhC4Kg9oQhC4Kg9oQhC4Kg9oQhC4Kg9oQhC4Kg9oQhC9oNkkzSKlXrCMonDFnQMJImS9q6ah1B0JowZEEQ1J4wZMECI+lbkp6U9I6khyStnXtvfUlPSHpf0o2Srpd0UoPfnSzpCEn/J+nd9N2uufePlPS6pNckjSjvFwfNRhiyYIGQtB5wBTAaWAa4BLhNUhdJnYFfA1cBSwMTgG838t3cLvYAtgdWBNYG9kvf3R44AtgGGADEkPdLTBiyYEEZBVxiZhPN7DMzGw/MADZOj07AuWb2iZndDDzS4HczzjWz18zsbeB2YN20fQ/gSjN72symA8cX+SOD5iYMWbCg9AMOT0PDdyS9A7QAK6THqzZ7ZYIpDX43Y2ru+QdAj/R8hVb/628L5+cEdSQMWbCgTAFONrMlc4/uZjYBeB3oI0m5z7c0+N158Xqr/9V3gX9JUFvCkAXzyyKSumYP4DLgQEkbyVlU0lBJiwEPA58BYyR1krQzsGHuf7X13XlxA7CfpDUkdQeOW8i/M6gRYciC+eU3wIe5xzDgAOB84J/AiySHvJl9DOwCjATeAYYDd+B+MMxs0ty+Oy/M7C7gHOAP6Xt/WPCfFtQVRWHFoEwkTQQuNrMrq9YStB+iRxYUiqQtJfVKQ8t98RCKu6vWFbQvOlUtIGj3rIb7sxYFXgJ2M7PXq5UUtDdiaBkEQe2JoWUQBLVnnoZM0hWS3pD0dBmCgiAI5pd5Di0lbQFMA642s7Ua+ac9e/a0/v37L7i6IAiCHI899thbZrZs6+3zdPab2f2S+s/Pzvr378+kSZPm5ytBEATzRNIcU9EWmo9M0ihJkyRNevPNN+fru/379EFS6Y/+ffq0raulfF39W+ahqV9Fx6rf3HW19G2pRFNL35a5aurTt28lmvr0bTtTqk9Lv2p0tfSbq6a+LdUcq74tCy+rrKFZS3mP7I5Gh5aDBg2y+emRSeKtfUc3/PmFRc/xl9DW75fEe+fsUaIiWHzsDfPUNON/RpaoyOmy8eVz1SWJSx4+v2RFMHqTMW1q2vf220tWBON33HGe52/UJeWPVi4dPajNY/Xg+Q+WrAg2G7NZm8dqTkh6zMwGtd4es5ZBENSeMGRBENSeRsIvJuBVDFaT9Iqk8sc1QRAEbdDIrOVeZQgJgiD4osTQMgiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2tOQIZO0vaTnJb0o6cdFiwqCIJgf5mnIJHUELgB2ANYA9pK0RtHCgiAIGqWRHtmGwItm9pKZfQxcB+xcrKwgCILGkZm1/QFpN2B7M9s/vf4esJGZjWn1uVHAqPRyNeD5hS93jvQE3ippX43SjJqgOXWFpsZpRl1la+pnZsu23thpYf13M7sUuHRh/b9GkTTJzAaVvd+2aEZN0Jy6QlPjNKOuZtHUyNDyVaAl9/oraVsQBEFT0IghexQYIGlFSZ2BPYHbipUVBEHQOPMcWprZp5LGAPcAHYErzOyZwpU1TunD2QZoRk3QnLpCU+M0o66m0DRPZ38QBEGzE5H9QRDUnjBkQRDUnjBkJSGpt6Slq9YRLBwkxb3TRHxpToakhRYz9wU5Czhb0saSulWsZTaym1JSL0ndq9aTkdLjkLSdpNUlqWpNAJIWN7OZkpaRdImkRSrUsmRV+24LScMlfbes/bV7QyZptfR0lKRDK5RyMPBX4FTgKEmrNoFxRZLSTdkFOAePE6RqoyGpg5l9JqkFuAiYbmYmaQVJvarUBuws6SrgKuAvZvZJmecyZ+C3Bw5OYVHNxjvA/pK2glmai6JdGzJJiwKjJb0EHG5m56btXUvW0cHM3jGz44ExeIDxecB+qRdUpdHI9j0OeMXM/pyKAhwtaXhVosxsZnp6EnCumU2RtD8wHjhO0hIVavsl8BGwFdBX0hJm9imApA1K2P9n6emZwH+b2ceStkhpgpWSXctmdgdwBbCnpC45zYXQrg2ZmU03s8OA54Cekm6Q1GJmHwFI2rtoX4ekjqnH01PScmb2TMpbPQ2vKPJLYM0iNbRF0rY4sDlwi6RvAz8DegGbSlqsbE2tzslEoKukO4F+eCWWDsAmZetK2jLDfwawE7AIMEnSIZIGADdLWqEEHbsBz5nZgynO88fArpIOLHrfc9GTnbOlJS0qqW8y+DOAazP/cFH3W7uNI8u1DCbpq8Cf8aHTPunvR8A6ZrZnkRosHWBJE4DuwKfArcCvzex9SfsCN5nZ9KJ0NKBzEfymPAgPkj4E+AvwILCXmf2pIl3nA/8FDAQ6mtkJSetEYBczm1yynkXSMHJ1/Di9ZmZvSxoEnA5MAx4wszNTL3xmm/9wwbSskva5AXAN8B/AAGC4me1b1H4b0HUusBbwGrAEcCdwPHBKNiIqZL/t1ZBlSNoJ6AzckwzHSsC/A+8Cx5jZKwXuu0Pq8RwDLI33vn4L/A7vVUwAbkvG9nOjVyaSVgUOxG+Gt4CZZvaKpBOBr5jZiDK15Y7ZLsAeZpYNTWak1vxc4CMzO6IMPTldmRHrCDyGN4R/BR4CbjazVyUtaWbvpM8XfsyS/3CdNIxD0r3AOWZ2a9GGtJWO2X6rpGWAD4BvAUvi1/4BwNlmdkEhGtqjIcvdDAcB+wGGtxC/ASaY2XRJXbMhZsFaeuAGbBTwc+A+4G7gfuAGMzu2aA1tkYZBPwJeN7PT0rZV8KHKsWY2tcybIu2/I3AZsBTe88p6tcvjQ8rbytST9n0h8MekaWkzGydpd2BLYCbwf6Rrq6D9d0yTH6sA3wQWw3v2LwCfAcOA3c3sO0Xsvw1d+VHHPknX62Z2c6vPrQycDHzfzD5c2DrapY8sGbFOwHbADma2Md772Q74D0nblGTEOpjZNOAIvAVfHLjbzKbi9druzD5XtJbWutLf9YGuuE9sI0l3SVrZzF4EjqzCiCWWBZ4BlgF+IekbAGb2dzO7pQIj1h14BNgW2Bs3HJjZjfiwaTKwfJHugZyz/Gp8RnAffOLjKGAV3A0wIuktdIawFUr7PBr3+fYCxkm6XdKQ3OdWBHoUYcQAMLN2+cC7tc8BI3Lbspt2SMH77pD+dk5/l05/DwNeBG7CexVVHp8uwL3Ay8CRwKZ4T3E80LsCPdkxU+64DcIbgYvwyZGVKjxeXYH+wLHpHJ4NDMi93y3/OxbyvrOR057Aden5o8CuwP8CjwPrVnhsegJP4UUlrgfGAqcALwEnpc8snd0HRTza69BS+I25I7Ay3u2/xcyeKlnHqbjj8xXgDfwktwCrAzea2WvZkKFMXTl9O+HxbQLuwmPIDsMN/b0VaToZP0YDgNG4S2AwsA1wvpm9VLKebEi3OvCquZ91XWAv4KvAk3iIyKdWcE9R0ja4H3MwMNDMDpG0NXAosKuZfVLk/tvQtSawKvAs8Esz2zCF8IwFzivlvqvKihfcQgh3MnbBe2Yn4LM6hwOdCt531rPYH59xWxefEfwJcFxrnRUdn+WArfFwhhZgD7x8OcCSFejJjtl2+ERIN+AfwOq589mt4mvqDnzG+bDctqHpulqiwP3uCXyj1bad8N70ENxP9r20vWOJx2NFvLEZkF4vhvdYr8J7X0cA15Slp930yHIt5zbA9/BgxZuAa/GhwDDgYzO7piQ9JwOPm9l/Jp/UQHw4cqKZ3VeGhja07QVsBvTA/StT8N7YCWb2u4r8YplD/df4mg/rmNkBkrYFhgMjreQeR+6aGoEb/Q648VoCONjM7sl9dqEfM0mj8Yb4+vx1m/y/hwHr4w3zbgtzvw3oagH+E3dLvGJmY3MTbOfhCxa9C4w1s2dLuZ6qaN0Kbikewlv2nsBPgSeA3fGLsOje2LrAaun53viM0ra5928HvlXRccl6PZ1wX8ay+A25Gx6B/QRwRsXnbis80HQSsFjaNh4Pk6lKU3fc17NWbts4PLzgVxTkT8R9co8Da+S2DccLnH4Nd5wDdE9/y+yNjQd+iPfCfo3Pxp+B+587443ikvnrruhHu+mRASQ/xsXAMJsVz7MdPoU/uuB9H4BPxV9uyb8kaSTeavbH/QdbmtmGReqYi7b8FPnFQG/gE+CnZvZc2t4P9wF9WlaPbA7xR0vjwcpr4VP1m+FD3s2K1jIHbVkPowNwPu77eTi91x2frewDPGQFxEZJGosbz2z1MuEzk8/ix+UyvIePlXgTS9oCH+n0Nu+tTsI7D5Nwv/RU81S8Uqk8aXlByW4Gef7ky3jU976SrjGzf+BO9kJXeUn7Hg0cYGZPZLrwE/4UPlT6FA/m/HzIUqSmueg8BpiOJ65fC7yU4tymm9nfss+VYcQySYCl4Nc+wLtmto+kvfGb9QU83KBU5MnYnSTdYx4EOwm4SdJpeODwIXhoyHnAsZIuM1/zdWHyN1LqWgqnaAHONLObJQ3EY/+6WlHhDHOnN264xkrqjRuuQ9P1/jweevEVKzDQfE7U3pDhw6RP8XSNv+DDt6OAFeXpLAPxqfsiOQh40syeyFry1Eq+K2kmnrYyOftwWUasVU+sKz5begI+m3S2ebR85v/5cRma8qQez9p43ffj8FzBkXhP8dqy9eR4Cw9rOEDScviwaSLeEO2ApyIdjE/gTCzAiIH3vI6QtJmZPQhMljQ1vXcg8IKZfVhGBkErbsdnkjcBvgE8kRmudI11KtuIQTuJ7JfUB7jazIak171wJ+k/gffM7L8K3v+ueODt/skROxM+v1HXxScaTrJUIaFM5LWz3kvP98JXie9lZl9P2/4InGZmd1RwUyBpc2BFM/tlatX3xns87wM7VdDjyHR1wuO01sP9rbea2e2alS61EnCJmW1ToIYD8WPxW9wvZ3hy/+n4sHNmFecsaeuJj3Q2AZbHw1BG4bO695U96qitIUsHcoSZnZFa8UPw1uJCM3u9ZC2r4Q7zH5nZH9O2LDfvctz39LMyNSUN2+LD2vFm9l4y8OcCH+I9joG4U73sWa9sNnBDPLygH+5zei756JYCNjGz35SpK2mbzT8or2ixCT4T1wm4Mxm0JYEuZvb3gvVsgBuIXfDsgrfw1LY7q4xBzOlbBT82u+I5sKUVU5xNR40N2cV45P4E/CR3BtbA/WSPAH+0EtKQcnqy1vMePEByJu7nORNYs+zWM91oLWb2lKRT8NijA4H/wVNZlsaHCHeZ2ZtVOPglPZ70rIi36A/jISulD01akz8e8qKTX8Xjth41s/sr0NMLj1d7vux9t9LRrfWQNvVeVwcmm9m0KgxsLQ1Z8n0dg09RbwjsZ2Yvpx7IpnjA57N476zMZOd86zkReBtvPe8ovavtZYY3MbMx6fVI3JA9hQfmTilLSytdWW9sbzy04NjkGhiON0Rv4hUcSjVm6do5DzjUUnxYmrHsZF64cH1gUzM7vwQtrWdzZ7t2yh5Oymv4TZHXqpPlEsKrGtr+C1ZRfM6CPnAj9izuS/lObvtSwPeBLSrU1gsf0lUSuZ80rIdXa7gEWCRtWw44EQ8QPrlCbZ2B3wOvAxvktg/Gh+dV6ToYT1a/gRQPmHvvWWD79Hyhn1dSHFi6bralVcxjVdcSHiu2Gx4OM5mUAdJsj8oFfMGDmwV3/hCvAjAReAD4egVa1Op1x7beL1nbcummPKjV9nWzY1WWPjzU4vjcuVsWjxl7Bl+YZdEKj5NyzxfB3QEvAaenbWOA35ak5Xq8usZleDxdKQGl89C0Hu6ueR5Pvdso994GWUNZ5aN2Q8vc0GRRPOXok7R9LN6iPof30Arzj+U0rIY7qv9guRnJpuluA/JSKscCF5nZDRXq6AlsjEerH4IPbz+WtBaeAzsET/u5o0KNnw/h5PWzzsRjufoAg83TbRa6i0Czgm//DS9ycCs+e/tVPF/3AjN7o4IhZf54rIkfi4F4L+1R3C82wAqsstwwVVvSL9hCbIxXa5iKd3mXStsXB8aUqKOpWk+8N7En7m/ql9u+Ge5U/3ZFuvrjPcNu6flEPBJ8r9xndsQr0lZ17LJGvUP+POLDvFHZewXuvwvuCtg8t20XPKD7WrzaRVXH5lJgu/S8N/BdPLbubmDjoo9NI4/a9Mjkq+bsiFd5vRLPc/tffIawBU8NuqwEHU3ZeiZtm+PBwC/jfo278NnTKXg1kI3x2Kw3StZ1FfCUmZ2V27YLHgT7CnCqpbCVknVtiEfIH2YpsyE7b3PqeRVxTiV1thRQK+mneIjHueZZKUi6Bg+X6YAb1LKCqbPrfFd8aL0tHsc2HI9rewN3o8xoihFIlVZ0PluFHfBYrXHAjcByufd2wsfvo0rS0sytZ1aUcBng2+m4jcerFbyAh2SUqWdLvKLpTnN5/3i8QSo0oX8u++6Wjs3zuFHtOofPFOZDxAOlxzIr8Xtl4OakZQRwIZ6itSjuo+pewTG6DNgIr5Ly83SuniOXzN4Mj9r0yAAkZWkRW+ELP1yFJ+1m8SyFhjg0a+uZ05e1onPqTfTAb8r3ywwFkfR74FXcDSB8oY6HW32mi5nNKENPbp8DmRWAuyZ+k/bBHfylpUelUKLj8Yot38d706PxbALwKP5xeB38cSXqEn6+foBXEJ6M+1knSLoar1RbesDy3KiFIctFyWc5XVvjzuEueGs60cyeLLKLK2kY7t+51Mw+SM7grIWagqdr9MDzLu/FZwU/KELL/FBWoOtc9v0DPGRhR3nd/U3wY/gXPL7urxXpEu7IPxlPj3osbR+G94bewGuzPVSwhh6pYemKxx+OwBelOd7M3k6fWx6Ph9y1KC2tdbW+hyRtilfAfUSeUH+ymRW+EPH80PSGLNfL6Icv1rGNmb2eIp23wx3Zz5nZ2SVoadrW08wsBXU+m4x9ZQYsp2trvPDen9LrFfAA5g3xoe+TeL5iVYZ2MXxVq8dwH+vEtP0kvL78wQXtdzB+/ayB91Czqih98N7PUDx3eFzaXmYPOruWRuDnaBruRsl6r7filV9vrCKCf240vSHLkHQDMMk8t7KTzVqivg+Amb1a4L6bsvVspbEb3kP8sxW4EOoXobVRla+lOQSYYWZXVKfs80mkI/BihfcDV5jZ5NwNXUTl12xt02wBkfF4aMPVqdHeHM+BvatkI5aFFQ3C/XOXA0fjwcvPALcBU8p2AzRE1U66Rh74EPIKoE96nUWqH4mnjRS578H4Sb0PT1/JtvfBAzpfwMvOZNtLq9Q5F60P4w7kThVr+RcneX5bOqelB1IyK4L+e8APc9tXBn6Bx7kNLnD/I4Fr0/PFcJ/qmXjv9EVylWgrPHdX437oofgoaG/cHXA5Fazp0JDmqgXMx8E9He/iDkqvF8f9U30L3u9v8ZCGb+CVQgfji/5mEeqb4yV8SjdiOQ09SDNueKT1aUD/Cs/V+rnnlWU2zEFXNgLphs+8ZWXJtyPFtOGrNRU5U/kM8HR6PhrvAWbvnQgcVdGxGZQal8XT9b0EHr7TL71/DfDdqs/h3B5NW1gx17WXOT/KovflpUOm4iVVXi5Qw0jgLfPh7GJ4DNtH+ApEx0oaZrn4JyvZX2CzhjyHAF+T9DB+o+wAbC9plLmDtsyqG32Be9Lw6aeWlm9rhlij3P73w3teL0s6HK9d9zdJMywlRBeodyBwnqRpeL5wS+69D/HA6lJJ0QCHmdnuwAzgquROeRk4Wl6KajA+kdUU57I1TekjyxmxzvhFtyoeXzMND/BcD6/IWuh6eZKewa//teQr2mxkZtlqzicC08zsjCI1NKCxK17xYzoe1f8ifjOsipcyOq4iXbfgs5QT8MVDpqftld8EaeLoKNyAXYSH8eyEn9+RJWq4EFgb74VNkPQAXmPvhTInayTdja+zermkr+Eji/vkqyVdBXyMh1uMbyYH/2xU3SWcSzc3GzJdhE+RX49HgP8cbxlK8a3g8WAX4Ab0U3Ir5uBO0KsrPj6d8BVr+ufey3xAq+JO2lKChFvtO7sBjsEXzHgJX8+gGa6tzfEsh3VIFVLScXyUCtJtcD/U83gDfVYF+x8JTMj2i7trNmr1mcp8rQ3/jqoFtHGAV8Hjw8AX8TgIj5r/B76qcplastCPKczypTzArMVJS80zY1ZP+tJk7N/DJ0PWaPW5cfjq3GXpygzstcw+AbInPmT5HdVEp2cGdufsmkqvu+K+oOHZcaIinx4ejtG5bA24K+Kp9PxgPE4y//5ZeGn00o/J/Dya1keGt5KnyPPhljKziyRdiRuUp8sUYp6HN1TSUODfJf0KX7yj1CEAzDbsXg+vPLtZikxfHrhfntd4El6n7TZ8qFmGrvyQ4xFg8eRn6Whm1yWNM6yCIOGcrqPx0k9IOhoPebgXT9/KKoMIzyksW+OVSVfZ8X95n103vIdP0nIEsIKZTZ3bl5uFDlULyJMufADM7E9mdivw9/TePvgN+rRVVO7XzO40s9Xw7vjR2eaSNWT7G4pfgLvgAbhD8XUOd8dLIn9mZo+a2T+L1iTpK8Bt8gqi4GsnbI0nZK+SosGH4sP0SpC0LF4csZc8lawnbnC/jrsqPoZSl8KbI2Xv33zFrx/gRv0u4BH5IjXguaA/g8+r5TYtTePszzuBJe2Hd/v/gYc/DMXjWlbCa429WZXOPFX0xkiBuel1Z+Aw3Hf3b5J+ArxtZheX7CzujSfNZwnip+JDyVPxir3vAw9bCWWi20LSGDw85Vkz+7l8KbqLzWzTKnU1E9moAxiAjzoOb4YskXnRTIYsS0U6Bg9OXARflGJL3I8xDR+mVLI8WNWktJb98JbzFjM7J21fFV+rcgbeu1jbfMWkKsoIDcDLB30NL2l0Id6JLGLdx0b0ZMPw3rhz/0ncj/eafKWmm3FDdn3TzsZVhKTvA78yL35Z+UzzvGgaQwYgX4r+N2b2dUln48uonSlpO7yn8WjFEisjl9byCJ7WcjXu35iAz1AuBXxgHjdWZlrLRrhxXQ4PmvwEb833wmcvf29mvyj7ZmiVbnMevhDMcviw9yx8aPlNM7uoLE11pA69MWg+Q9YLOBQPddjKzDZP2x8FTjSz26vUVxUpMHeIme2dAnPfwLMMtsEX8hhmZn+uQNexwBb4zNdKuKP4JjxMRrjB7Wdmp5etLafxeuB2M7sm9RhPwmfhtsyNApq+xxG0TeWzlvI1A3ub2WQzmyrpEbyUyg2S+uNVYd/6shqxxGHMmlTYG4/7ORI+r9QwDC89XBqp0RkBrJPz2W2QdAxNmm6mghnA3JCyGx4ovEjqob0AfEfSBEmrmNmLMNsESlBTKjVk6Wa4COgk6SM8kfdOvNu/Dj58egZPgv4y01Zay3QqSGvBnfu/N68I0gOYbl7Xa4ik6/Bl3u4pW5RmFd88xXwh2VvxSrlPSnob6I6XESp0hfCgXKqeUj0Rz+caiQ+XhqRt95rZIXid8LFVhVs0C62myO9h9inyb+KBr2VPkT8I9JSXVJqWekCLpfceAbbKh9OUyAm4PyxrKG/DK5QcD5yCV7g4ORngjhXoCwqgsh5ZyukaaGabpNfbp7c6A49LOi6bmQucZgnMTQbqBbwhvFDSuWb2dDbExP1m91Ywa5ol+V+cjPr1wHAzO1nSinhlh6lm9ncoP8k/KI7KnP2ShuNDxy3xxRUONLNh6b3BeK3wg76s4RaNUPUUuTzx+UC8jNDbeL7iQGDHKmKz5En+M81soLzM9rpmdkCEVrR/Kp21lLQMXh3zm/h6lBem7TsDPzGzjSoTVyOqnCJP8Vjb4Mn8mwPXAfeb2RMVaOmAh1rsS/IlmpdF74jHj31StqagHJoi/EK+2vQt+KKtI/Gk4wurcBYH9Sf1FC/AS+QcZymPMWi/NIUhg8/9LnviPbQHzGxIxZKCmpN8iWfhrotN8QDrpg/uDOafpjFkGSl/sIelBT2CYEGRdCieilRJqlRQPE1nyIKgKOqSbhPMP2HIgiCoPVUHxAZBECwwYciCIKg9YciCIKg9YciCIKg9YciCIKg9YciCIKg9/w/vv35ABZt+lgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 360x72 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "ESG_plot(30)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "tags": [
     "remove_input"
    ]
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Price: $50\n",
      "Total Production/Market Supply:  18650\n",
      "\n",
      "Suppliers:  ['BIG CREEK' 'HELMS' 'DIABLO CANYON 1' 'MOSS LANDING 6' 'MOSS LANDING 7'\n",
      " 'MOHAVE 1' 'MOHAVE 2' 'FOUR CORNERS' 'MORRO BAY 3&4' 'PITTSBURGH 5&6'\n",
      " 'ORMOND BEACH 1' 'ORMOND BEACH 2' 'MORRO BAY 1&2' 'MANDALAY 1&2'\n",
      " 'CONTRA COSTA 6&7' 'HUNTINGTON BEACH 1&2' 'PITTSBURGH 1-4'\n",
      " 'EL SEGUNDO 3&4' 'ENCINA' 'REDONDO 5&6' 'REDONDO 7&8' 'COOLWATER'\n",
      " 'ETIWANDA 1-4' 'SOUTH BAY' 'EL SEGUNDO 1&2' 'HUMBOLDT'\n",
      " 'HUNTERS POINT 1&2' 'HIGHGROVE']\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiYAAAGECAYAAAAC8us+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nO3de5QlZX3u8e8Dw01gAGWC0MyFGDUhGG+jgpgY1HghIpgQQkJwTMhiYkIi0SRqkpV4cvecaNCDmiFBHSOKiDcwGkWDeoKADggCIoqXYYZGGCI4RAVBf+ePXY2bprunp7t377env5+1avWuql1Vv9013f3MW2+9lapCkiSpBTsNuwBJkqQxBhNJktQMg4kkSWqGwUSSJDXDYCJJkpphMJEkSc0wmEg7mCRvS/I33eufT7J52DVJ0nQZTKQFKsknk9yRZLdZ7KOSfCfJ//RNfzKXdQ5Ckqcl+UySbyf5VpJLkjxp2HVJmr0lwy5A0vZLsgr4WeDbwAuA98xid4+tqhvnoKxJJVlSVffN0b6WAh8CXgKcB+xK73txz1zsf9yx5qxuSdNji4m0ML0IuAx4G7BmEAdI8uok5yV5e5K7klyXZHXf+oOSvDfJliRfT/IH47Y9P8k7kmwFXpxkjyTru1ae65P8ydhlpiR/nOS9447/hiSvn6C0RwFU1buq6gdV9b2q+lhVfaHbbqckf55kY5Lbuvr36dY96NJWkm8kedYUde+c5E+TfLX7PlyRZHn3/p9MclHXanNDkhP69nt0ki9229yc5I9mdUKkRcJgIi1MLwLO6abnJDlgQMd5AXAusC9wAXAm9P74AxcCVwMjwDOB05M8p2/bY4Hzu23PAf4SWAX8OPALwG/0vfcdwHOT7NvtfwlwIvD2CWr6MvCDLuQ8L8l+49a/uJuO6o6111jd0zS+7pcBvwYcDSwFfgv4bpI9gYuAdwI/1tX7piSHdvs5G1hbVXsDhwH/uR01SIuWwURaYJI8DVgJnFdVVwBfBX59Fru8MsmdfVN/uPivqvpwVf0A+Dfgsd3yJwHLquqvqur7VfU14F/o/XEec2lVfaCqflhV3wNOAP6uqu6oqs3AG8beWFW3AJ8GfqVb9Fzg9u7zPUBVbQWeBlR3zC1JLugLZycBr6uqr1XV/wCvAk7sws50jK/7t4E/r6obqufqqvpv4PnAN6rqrVV1X1V9Hnhv32e4Fzg0ydLuM185zeNLi5rBRFp41gAfq6rbu/l3MrvLOU+oqn37po/2rftm3+vvArt3f+BXAgf1BxrgT4H+lptN445z0Lhl49ev50etKL9BLwhNqKqur6oXV9XB9FojDgLO6DvOxr63b6TXn266rUrj61pOL/yNtxJ4yrjvwUnAw7v1v0yvlWVjkk8lOWKax5cWNTu/SgtIkj3otTzsnGQsNOwG7JvksVV19TyVsgn4elU9cor3jH90+S3AwcAXu/nl49Z/AHhzksPotUZM6+6gqvpSkrcBa7tFo/RCw5gVwH3ArfRCy0PGViTZGVi2jbo3AY8Arp1g+aeq6hcmqetzwLFJdgFOo9dRd/xnljSOLSbSwnIc8APgUOBx3fRTwP+j1+9kvnwWuCvJK7pOrTsnOWwbt+yeB7wqyX5JRuj9sb5fVd1Nr2/HO4HPVtVNE+2k63D68iQHd/PL6fUBuax7y7uAP0xySJK9gL8D3t3dXfNleq0+v9gFhj+nF+ym8q/AXyd5ZHp+JsnD6N0Z9KgkJyfZpZuelOSnkuya5KQk+1TVvcBW4IfbOI4kDCbSQrMGeGtV3VRV3xyb6HXuPGk7+lH0uzoPHMfkjG1t0PU5eT69YPR14HZ6f8D3mWKzvwI2d+//OL0QMv4W3/XAY5jiMg5wF/AU4PIk36EXSK4FXt6tf0u3/ae7Y90N/H5X97eB3+1qvRn4TlfTVF5HL1R9jF7AOBvYo6ruAp5Nr1/NKL3LXq/hR0HnZOAb3d09v0PvMo+kbUjV+FZLSRq8JC8BTqyqp/ctWwF8CXh418lV0iJji4mkeZHkwCRHduOMPJpeC8f7+9bvRO/W3HMNJdLiZedXSfNlV2AdcAhwJ73xUd4E0I0Jciu9O2ieO6wCJQ2fl3IkSVIzvJQjSZKaYTCRJEnNWBB9TPbff/9atWrVsMuQJElz4Iorrri9qsYPbggskGCyatUqNmzYMOwyJEnSHEiycbJ1XsqRJEnNMJhIkqRmGEwkSVIzDCaSJKkZBhNJktQMg4kkSWqGwUSSJDXDYCJJkpphMJEkSc0wmEiSpGYYTCRJUjMMJpIkqRkGE0mS1AyDiSRJYmT5SpLcP40sXzmUOpYM5aiSJKkpo5tv4tR1G+6fP2vt6qHUYYuJJElqhsFEkiQ1w2AiSZKaYTCRJEnNMJhIkqRmGEwkSVIzDCaSJKkZBhNJktQMg4kkSWqGwUSSJDXDYCJJkpphMJEkSc0wmEiSpGYYTCRJUjMMJpIkqRkGE0mS1AyDiSRJaobBRJIkNcNgIkmSmmEwkSRJzTCYSJKkZhhMJElSMwwmkiSpGQYTSZLUDIOJJElqhsFEkiQ1w2AiSZKaYTCRJEnNMJhIkqRmDDSYJPnDJNcluTbJu5LsnuSQJJcnuTHJu5PsOsgaJEnSwjGwYJJkBPgDYHVVHQbsDJwIvAb4p6r6CeAO4JRB1SBJkhaWQV/KWQLskWQJ8BDgFuAZwPnd+vXAcQOuQZIkLRADCyZVdTPwj8BN9ALJt4ErgDur6r7ubZuBkUHVIEmSFpZBXsrZDzgWOAQ4CNgTeO52bH9qkg1JNmzZsmVAVUqSpJYM8lLOs4CvV9WWqroXeB9wJLBvd2kH4GDg5ok2rqqzqmp1Va1etmzZAMuUJEmtGGQwuQk4PMlDkgR4JvBF4GLg+O49a4APDrAGSZK0gAyyj8nl9Dq5Xglc0x3rLOAVwMuS3Ag8DDh7UDVIkqSFZcm23zJzVfWXwF+OW/w14MmDPK4kSVqYHPlVkiQ1w2AiSZKaYTCRJEnNMJhIkqRmGEwkSVIzDCaSJKkZBhNJktQMg4kkSWqGwUSSJDXDYCJJ0iIzsmIFSR4wtWKgQ9JLkqThWTUywsbR0QnXrbnwwgfMrz/mmPkoaZsMJpIk7aA2jo5y+5q1D1q+//p1Q6hmeryUI0mSmmEwkSRJzTCYSJKkZhhMJElSMwwmkiSpGQYTSZIWuFXLRx40LklLY5NsD28XliRpgdu4eZStZ5zwoOVLTz9vCNXMji0mkiSpGQYTSZLUDIOJJElqhsFEkiQ1w2AiSZKaYTCRJEnNMJhIkqRmGEwkSVIzDCaSJKkZBhNJktQMg4kkSWqGwUSSpAVoxfIVC/phfZMxmEiStABt2ryJS868hEvOvGTYpcwpg4kkSWqGwUSSJDXDYCJJkpphMJEkSc0wmEiSpGYYTCRJUjMMJpIkqRkGE0mS1AyDiSRJaobBRJIkNcNgIknSArFq+cgO+XycfgYTSZIWiI2bR9l6xglsPeOEYZcyMAYTSZLUDIOJJEkNG1m+coe/fNPPYCJJUsNGN9/Eqes2cOq6DcMuZV4YTCRJUjMMJpIkqRkGE0mS1AyDiSRJaobBRJIkNcNgIkmSmmEwkSRJzTCYSJLUmFUjO/4zcSZjMJEkqTEbR0e5fc1abl+zdtilzDuDiSRJaobBRJIkNcNgIkmSmmEwkSRJzRhoMEmyb5Lzk3wpyfVJjkjy0CQXJflK93W/QdYgSdJCsGrl4r0Tp9+gW0xeD/xHVf0k8FjgeuCVwCeq6pHAJ7p5SZIWtY03jXLPZadwz2WnDLuUoRpYMEmyD/BzwNkAVfX9qroTOBZY371tPXDcoGqQJEkLyyBbTA4BtgBvTfL5JP+aZE/ggKq6pXvPN4EDBliDJElaQAYZTJYATwDeXFWPB77DuMs2VVVATbRxklOTbEiyYcuWLQMsU5IktWKQwWQzsLmqLu/mz6cXVG5NciBA9/W2iTauqrOqanVVrV62bNkAy5QkaW6sWL7i/g6s2zupZ8mgdlxV30yyKcmjq+oG4JnAF7tpDfAP3dcPDqoGSZLm06bNm7jkzEtmtO2Rpx05x9UsTAMLJp3fB85JsivwNeA36bXSnJfkFGAjcMKAa5AkSQvEQINJVV0FrJ5g1TMHeVxJkrQwOfKrJElqhsFEkqRZ6O/wqtkzmEiSNAtjHV5n2ulVD2QwkSRJzTCYSJKkZhhMJElSMwwmkiSpGdsMJkkeNBTdRMskSZJmazotJv93msskSZJmZdKRX5McATwVWJbkZX2rlgI7D7owSZK0+Ew1JP2uwF7de/buW74VOH6QRUmSpMVp0mBSVZ8CPpXkbVW1ESDJTsBeVbV1vgqUJEmLx3T6mPx9kqVJ9gSuBb6Y5I8HXJckSVqEphNMDu1aSI4DPgIcApw80KokSdKiNJ1gskuSXegFkwuq6l6gBluWJElajKYTTNYB3wD2BD6dZCW9DrCSJElzapvBpKreUFUjVXV09WwEjpqH2iRJatbyFctJMuwydjjTGfl1nySvS7Khm15Lr/VEkqRFa/Omzay79Mxhl7HDmc6lnLcAdwEndNNW4K2DLEqSJC1OUw2wNuYRVfXLffP/K8lVgypIkiQtXtNpMflekqeNzXQP8Pve4EqSJEmL1XRaTH4HeHuSfbr5O4AXD6wiSZKGZGT5SkY33zTsMha1bQaTqroaeGySpd28twpLknYYq0ZG2Dg6ev/8qes2TGu7s9auHlRJi9qkl3KSvCzJKWPzVbW1qrYmOSXJ6fNTniRJg7VxdJTb16zl9jVrh12KmLqPyUnA2ydY/m/Abw2mHEmStJhNFUyWdMPPP0BVfR9wRBlJkjTnpgomOyU5YPzCiZZJkiTNhamCyf8B/j3J05Ps3U0/D3wI+Md5qU6SJC0qqZr8QcFJnge8EjiM3hOFrwP+oao+Mj/l9axOanp9pCVJUusCV1TVhLc1TXm7cBdA5jWESJKkxWs6A6wN3xOfCBtsM5Ekzb0k998qvP/6dds1jsm6S89k7RGnccmZlwBw5GlH3v96ex152pHcc1lvlI7dDj97wprOWruarWecAMDS0897wHHHlvdbevp5E94Gvf/6day58MIHLFt/zDEP+OxnrV3NVFdVZmWKpzJPZ0h6SZJ2KMtXLCcJmeIPpIZjYbSYSJI0SyMrVjC6adP98+suPROAtUecNqySNIFttpgkOSDJ2Uk+0s0f2j8irCRJC8Hopk2sufDCB13CUFumcynnbcBHgYO6+S8DDkkvSZLm3HSCyf5VdR7wQ4Cqug/4wUCrkiRJi9J0gsl3kjyM3jgmJDkc+PZAq5IkSYvSdDq/vgy4AHhEkkuAZcDxA61KkiQtStsMJlV1ZZKnA4+m9/C+GyZ6uJ8kSS0ZfxeOFoZJg0mSX5pk1aOSUFXvG1BNkiTN2thdOGPWH3PMEKvRdE3VYjLVGSzAYCJJkubUpMGkqn5zPguRJEmazgBrD0vyhiRXJrkiyeu7u3QkSZLm1HRuFz4X2AL8Mr27cbYA7x5kUZIkaXGazu3CB1bVX/fN/02SXx1UQZIkafGaTovJx5KcmGSnbjqB3hD1kiRJc2qq24Xvonf3Teg9G+cd3aqdgP8B/mjg1UmSpEVlqrty9p7PQiRJ2h4jy1cyuvmmYZehOTadPiYk2Q94JLD72LKq+vSgipIkaVtGN9/Eqes2TLr+rLWr57EazZVtBpMkvw28FDgYuAo4HLgUeMZgS5MkSYvNdDq/vhR4ErCxqo4CHg/cOdCqJEnSojSdYHJ3Vd0NkGS3qvoSvQf6SZI0ayPLV5JkuyftmKbTx2Rzkn2BDwAXJbkD2DjYsiRJi8W2+opMxj4kO6ZtBpOqemH38tVJLgb2Af5joFVJkqRFaapxTJZW1dYkD+1bfE33dS/gWwOtTJIkLTpTtZi8E3g+cAU/Gmit/+uPD7w6SVKzVo2MsHF0dNhlaAcz1QBrz0+vd9HTq8oRbCRpBzbTkHH7mrWzPvb+69fNeh/acUzZx6SqKsm/A4+Zp3okSZ1Vy0fYuHnisPCQ3Xbhu/fcO6fH296QYaDQIEznrpwrkzypqj438GokaQEbxBDpW884YcLlS08/b9J1M7H09PPmbF/SbEwnmDwFOCnJRuA7dH1MqupnBlqZJC0wM73tdTLeDqvFaDrB5DmzOUCSnYENwM1dv5VDgHOBh9HrWHtyVX1/NseQpLnmA+Kk4ZjOOCYbAZL8GH0P8dsOLwWuB5Z2868B/qmqzk3yz8ApwJtnsF9J2qYVy1ewafOmGW27va0ftnBIszedh/i9AHgtcBBwG7CSXtD46WlsezDwi8DfAi/r7vJ5BvDr3VvWA6/GYCJplqbqKHrJmZds9/6OPO3I2ZYkaQamcynnr+k9UfjjVfX4JEcBvzHN/Z8B/Amwdzf/MODOqrqvm98MjGxHvZI0oY2bRyfsDGqnTmlhmc5D/O6tqv8GdkqyU1VdDGyzvTLJ84HbquqKmRSW5NQkG5Js2LJly0x2IUmSFpjptJjcmWQv4NPAOUluo3d3zrYcCbwgydH0+qYsBV4P7JtkSddqcjBw80QbV9VZwFkAq1evrmkcT5IkLXDTaTE5Fvge8If0Ht73VeCYbW1UVa+qqoOrahVwIvCfVXUScDFwfPe2NcAHZ1C3pB3QiuUrtuux9/2TpB3DVA/xeyPwzqrq7zW2fg6O+Qrg3CR/A3weOHsO9ilpAZlq+POZdFQFO6tKO4qpLuV8GfjHJAcC5wHvqqrPz+QgVfVJ4JPd668BT57JfiQtHKtWjrDxpsmfvTLR8OcOcS5pqof4vR54fZKV9C7FvCXJHsC76IWUL89TjZIWoI03jXLPZadMuG63w20olTSxbfYxqaqNVfWaqno88GvAcfTGMZEkSZpT2wwmSZYkOSbJOcBHgBuAXxp4ZZLmzKqRkRl3Kt1z113tjCpp3kzV+fUX6LWQHA18lt7zbU6tquncKiypIRtHR7f7kfZj9l+/bkbb2l9E0kxM1fn1VcA7gZdX1R3zVI8kSVrEpur8+oz5LESSJGk6A6xJkiTNi+kMSS8tWFM9cXbYHrL7Lnz37nu3e7vdd9udu++5ewAVSdLwGUy0Q5vsibMtWHr6eZOO8zGV3Q4/e7tHR3VUVEkLhZdyJElSMwwmkiSpGQYTSZLUDIOJJElqhsFEkiQ1w2CigVu+YvmMn9Oy+x67z3hbn9ciSQuPtwtr4DZv2sy6S8+c0bZrjzhtxtuObS9JWjhsMZEkSc0wmEiSpGYYTCRJUjMMJpIkqRkGE0mS1AyDySKxauXIjG+53XOPXb1lV5I0L7xdeJHYeNPojJ5kC72n2c5027HtJUmaDltMJElSMwwmkiSpGQYTSZLUDIOJJElqhsFEkiQ1w2AiSZKaYTCRJEnNMJhIkqRmGEwkSVIzDCaSJKkZBhNJktQMg4kkSWqGwUSSJDXDYCJJkpphMJEkSc0wmEiSpGYYTCRJUjMMJpIkqRkGE0mS1AyDiSRJaobBRJIkNcNgIkmSmmEwkSRJzTCYSJKkZhhMJElSMwwmkiSpGQYTSZLUDIOJJElqhsFEkiQ1w2AiSZKaYTCRJEnNMJhIkqRmGEwkSVIzDCaSJKkZBhNJktQMg4kkSWqGwUSSJDVjYMEkyfIkFyf5YpLrkry0W/7QJBcl+Ur3db9B1SBJkhaWQbaY3Ae8vKoOBQ4Hfi/JocArgU9U1SOBT3TzkiRJgwsmVXVLVV3Zvb4LuB4YAY4F1ndvWw8cN6gaJEnSwjIvfUySrAIeD1wOHFBVt3SrvgkcMB81SJKk9g08mCTZC3gvcHpVbe1fV1UF1CTbnZpkQ5INW7ZsGXSZkiSpAQMNJkl2oRdKzqmq93WLb01yYLf+QOC2ibatqrOqanVVrV62bNkgy5QkSY0Y5F05Ac4Grq+q1/WtugBY071eA3xwUDVIkqSFZckA930kcDJwTZKrumV/CvwDcF6SU4CNwAkDrEGSJC0gAwsmVfVfQCZZ/cxBHVeSJC1cjvwqSZKaYTCRJEnNMJhIkqRmGEwkSVIzDCaSJKkZBhNJktQMg4kkSWqGwUSSJDXDYCJJkpphMJEkSc0wmEiSpGYYTCRJUjMMJpIkqRkGE0mS1AyDiSRJaobBRJIkNcNgIkmSmmEwkSRJzTCYSJKkZhhMJElSMwwmkiSpGQYTSZLUDIOJJElqhsFEkiQ1w2DSoJEVK0gyp5MkSQvBkmEXoAcb3bSJNRdeOKf7XH/MMXO6P0mSBsEWE0mS1AyDiSRJaobBRJIkNcNgIkmSmmEwkSRJzTCYSJKkZhhMJElSMwwmkiSpGQYTSZLUDIOJJElqhsFEkiQ1w2AiSZKaYTCRJEnNMJhIkqRmGEwkSVIzDCaSJKkZBhNJktQMg4kkSWqGwUSSJDXDYCJJkpphMJEkSc0wmEiSpGYYTCRJUjMMJpIkqRkGE0mS1AyDiSRJaobBRJIkNcNgIkmSmmEwkSRJzTCYSJKkZhhMJElSMwwmkiSpGQYTSZLUDIOJJElqxlCCSZLnJrkhyY1JXjmMGiRJUnvmPZgk2Rl4I/A84FDg15IcOt91SJKk9gyjxeTJwI1V9bWq+j5wLnDsEOoYiJEVK0gyq0mSpMVqyRCOOQJs6pvfDDxlCHUMxOimTay58MJZ7WP9McfMUTWSJC0sqar5PWByPPDcqvrtbv5k4ClVddq4950KnNrNPhq4YUAl7Q/cPqB9a3Y8N23z/LTLc9M2zw+srKplE60YRovJzcDyvvmDu2UPUFVnAWcNupgkG6pq9aCPo+3nuWmb56ddnpu2eX6mNow+Jp8DHpnkkCS7AicCFwyhDkmS1Jh5bzGpqvuSnAZ8FNgZeEtVXTffdUiSpPYM41IOVfVh4MPDOPYEBn65SDPmuWmb56ddnpu2eX6mMO+dXyVJkibjkPSSJKkZizaYOCz+cCT5RpJrklyVZEO37KFJLkryle7rft3yJHlDd46+kOQJfftZ073/K0nWDOvzLHRJ3pLktiTX9i2bs/OR5Ind+b6x29YRBLfDJOfn1Ulu7n6GrkpydN+6V3Xf6xuSPKdv+YS/77qbEC7vlr+7uyFB05BkeZKLk3wxyXVJXtot9+dntqpq0U30Ot1+FfhxYFfgauDQYde1GCbgG8D+45b9b+CV3etXAq/pXh8NfAQIcDhwebf8ocDXuq/7da/3G/ZnW4gT8HPAE4BrB3E+gM9270237fOG/ZkX0jTJ+Xk18EcTvPfQ7nfZbsAh3e+4naf6fQecB5zYvf5n4CXD/swLZQIOBJ7Qvd4b+HJ3Dvz5meW0WFtMduhh8RegY4H13ev1wHF9y99ePZcB+yY5EHgOcFFVfauq7gAuAp4730XvCKrq08C3xi2ek/PRrVtaVZdV77fs2/v2pWmY5PxM5ljg3Kq6p6q+DtxI73fdhL/vuv99PwM4v9u+/1xrG6rqlqq6snt9F3A9vZHN/fmZpcUaTCYaFn9kSLUsNgV8LMkV3ei+AAdU1S3d628CB3SvJztPnr/BmqvzMdK9Hr9cs3dadzngLWOXCtj+8/Mw4M6qum/ccm2nJKuAxwOX48/PrC3WYKLheVpVPYHe06V/L8nP9a/s/mfgrWKN8Hw06c3AI4DHAbcArx1uOYtbkr2A9wKnV9XW/nX+/MzMYg0m0xoWX3Ovqm7uvt4GvJ9eM/OtXbMl3dfburdPdp48f4M1V+fj5u71+OWahaq6tap+UFU/BP6F3s8QbP/5+W96lxOWjFuuaUqyC71Qck5Vva9b7M/PLC3WYOKw+EOQZM8ke4+9Bp4NXEvvez/WE30N8MHu9QXAi7re7IcD3+6aSD8KPDvJfl0z9rO7ZZobc3I+unVbkxze9Wd4Ud++NENjf/Q6L6T3MwS983Nikt2SHAI8kl7nyQl/33X/m78YOL7bvv9caxu6f9NnA9dX1ev6VvnzM1vD7n07rIleD+kv0+ut/mfDrmcxTPTuCri6m64b+77Tu9b9CeArwMeBh3bLA7yxO0fXAKv79vVb9Dr33Qj85rA/20KdgHfRuxxwL71r2KfM5fkAVtP7w/lV4Ey6QR2dZnV+/q37/n+B3h+7A/ve/2fd9/oG+u7gmOz3Xfcz+dnuvL0H2G3Yn3mhTMDT6F2m+QJwVTcd7c/P7CdHfpUkSc1YrJdyJElSgwwmkiSpGQYTSZLUDIOJJElqhsFEkiQ1w2AiSZKaYTCRdL8kD09ybpKvds8z+nCSRw34mAclOb97/bgkR89wP69OUkl+om/Z6d2y1UlemuSMvnXrkny8b/73k7xhNp9F0uwZTCQB949k+X7gk1X1iKp6IvAqfvQQsoGoqtGqGht99HH0BqmaqWvojWw65lfoDeYHcAnw1L51jwX2SbJzN/9U4DOzOLakOWAwkTTmKODeqvrnsQVVdTXw+SSfSHJlkmuSHAu9J6om+VKSc5Jcn+T8JA/p1v1Fks8luTbJWV3oIclPJPl4kqu7/T2i28+13XDpfwX8apKrkvxqkq8kWdZtu1OSG8fmJ/EBeo+XJ8kjgG8Dt3frrgIelWSPJPsA3+uWPaZb/1R64UXSEBlMJI05DLhiguV3Ay+s3lOhjwJeOxY0gEcDb6qqnwK2Ar/bLT+zqp5UVYcBewDP75afA7yxqh5LLwiMPR6eqvo+8BfAu6vqcVX1buAdwEndW54FXF1VW6b4DFuBTUkOo9dy8u6+/d8HfB54EnA4vUfUXwY8NckIveG+Nz14l5Lmk8FE0rYE+LskX6D37I8RfnR5Z1NVjbUyvIPe80MAjkpyeZJrgGcAP909wHGkqt4PUFV3V9V3t3Hst9B7eBn0nify1mnUey69UHIcvUtT/T5DLxA9Fbi0m8bmvYwjNcBgImnMdcATJ1h+ErAMeGJVPQ64Fdi9Wzf+YVuVZHfgTcDxVfUY4F/63r9duhaMW5M8A3gy8JFpbPYh4GTgpqraOm7dWD+TI+iFkuuBQzGYSM0wmEga85/AbklOHVuQ5GeAlcBtVXVvkqO6+TErkhzRvf514L/4UQi5PclewPEAVXUXsDnJcd2+dxvrk9LnLmDvccv+lV5rzHuq6gfb+hBdK8wrgL+dYPWl9C7jLKuq26r3FNMt9Pql2L9EaoDBRBIA3R/pF5sN/BgAAADRSURBVALP6m4Xvg74e+DDwOrussyLgC/1bXYD8HtJrgf2A95cVXfSayW5Fvgo8Lm+958M/EF3WegzwMPHlXExcOhY59du2QXAXkzvMs7YZzm3qq6cYPkd9ILIdX2LLwV+DLh6uvuXNDjp/S6SpO2TZBXwoa6D6yCPsxr4p6r62UEeR1Iblgy7AEmaTJJXAi/hR3fmSNrB2WIiaUFJ8mf0Bk7r956qmqhPiaQFxmAiSZKaYedXSZLUDIOJJElqhsFEkiQ1w2AiSZKaYTCRJEnN+P90W51luaH1+AAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 648x432 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAATIAAACRCAYAAABAOTQxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAeWElEQVR4nO2dd7xV1ZXHvz9AmthRQfIAC2pUrGCPJmINUYktaogaUFCDDrElRhML9uho7GUsGCO2MdaomSQ6GnVQLDO2GI0hYiFqjAVULKz5Y+0jhxd4XIRzzj3P9f187ufde+697/zuKWvvvfZaa8vMCIIgqDMdqhYQBEGwoIQhC4Kg9oQhC4Kg9oQhC4Kg9oQhC4Kg9oQhC4Kg9oQhC9oNkkzSKlXrCMonDFnQMJImS9q6ah1B0JowZEEQ1J4wZMECI+lbkp6U9I6khyStnXtvfUlPSHpf0o2Srpd0UoPfnSzpCEn/J+nd9N2uufePlPS6pNckjSjvFwfNRhiyYIGQtB5wBTAaWAa4BLhNUhdJnYFfA1cBSwMTgG838t3cLvYAtgdWBNYG9kvf3R44AtgGGADEkPdLTBiyYEEZBVxiZhPN7DMzGw/MADZOj07AuWb2iZndDDzS4HczzjWz18zsbeB2YN20fQ/gSjN72symA8cX+SOD5iYMWbCg9AMOT0PDdyS9A7QAK6THqzZ7ZYIpDX43Y2ru+QdAj/R8hVb/628L5+cEdSQMWbCgTAFONrMlc4/uZjYBeB3oI0m5z7c0+N158Xqr/9V3gX9JUFvCkAXzyyKSumYP4DLgQEkbyVlU0lBJiwEPA58BYyR1krQzsGHuf7X13XlxA7CfpDUkdQeOW8i/M6gRYciC+eU3wIe5xzDgAOB84J/AiySHvJl9DOwCjATeAYYDd+B+MMxs0ty+Oy/M7C7gHOAP6Xt/WPCfFtQVRWHFoEwkTQQuNrMrq9YStB+iRxYUiqQtJfVKQ8t98RCKu6vWFbQvOlUtIGj3rIb7sxYFXgJ2M7PXq5UUtDdiaBkEQe2JoWUQBLVnnoZM0hWS3pD0dBmCgiAI5pd5Di0lbQFMA642s7Ua+ac9e/a0/v37L7i6IAiCHI899thbZrZs6+3zdPab2f2S+s/Pzvr378+kSZPm5ytBEATzRNIcU9EWmo9M0ihJkyRNevPNN+fru/379EFS6Y/+ffq0raulfF39W+ahqV9Fx6rf3HW19G2pRFNL35a5aurTt28lmvr0bTtTqk9Lv2p0tfSbq6a+LdUcq74tCy+rrKFZS3mP7I5Gh5aDBg2y+emRSeKtfUc3/PmFRc/xl9DW75fEe+fsUaIiWHzsDfPUNON/RpaoyOmy8eVz1SWJSx4+v2RFMHqTMW1q2vf220tWBON33HGe52/UJeWPVi4dPajNY/Xg+Q+WrAg2G7NZm8dqTkh6zMwGtd4es5ZBENSeMGRBENSeRsIvJuBVDFaT9Iqk8sc1QRAEbdDIrOVeZQgJgiD4osTQMgiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2hOGLAiC2tOQIZO0vaTnJb0o6cdFiwqCIJgf5mnIJHUELgB2ANYA9pK0RtHCgiAIGqWRHtmGwItm9pKZfQxcB+xcrKwgCILGkZm1/QFpN2B7M9s/vf4esJGZjWn1uVHAqPRyNeD5hS93jvQE3ippX43SjJqgOXWFpsZpRl1la+pnZsu23thpYf13M7sUuHRh/b9GkTTJzAaVvd+2aEZN0Jy6QlPjNKOuZtHUyNDyVaAl9/oraVsQBEFT0IghexQYIGlFSZ2BPYHbipUVBEHQOPMcWprZp5LGAPcAHYErzOyZwpU1TunD2QZoRk3QnLpCU+M0o66m0DRPZ38QBEGzE5H9QRDUnjBkQRDUnjBkJSGpt6Slq9YRLBwkxb3TRHxpToakhRYz9wU5Czhb0saSulWsZTaym1JSL0ndq9aTkdLjkLSdpNUlqWpNAJIWN7OZkpaRdImkRSrUsmRV+24LScMlfbes/bV7QyZptfR0lKRDK5RyMPBX4FTgKEmrNoFxRZLSTdkFOAePE6RqoyGpg5l9JqkFuAiYbmYmaQVJvarUBuws6SrgKuAvZvZJmecyZ+C3Bw5OYVHNxjvA/pK2glmai6JdGzJJiwKjJb0EHG5m56btXUvW0cHM3jGz44ExeIDxecB+qRdUpdHI9j0OeMXM/pyKAhwtaXhVosxsZnp6EnCumU2RtD8wHjhO0hIVavsl8BGwFdBX0hJm9imApA1K2P9n6emZwH+b2ceStkhpgpWSXctmdgdwBbCnpC45zYXQrg2ZmU03s8OA54Cekm6Q1GJmHwFI2rtoX4ekjqnH01PScmb2TMpbPQ2vKPJLYM0iNbRF0rY4sDlwi6RvAz8DegGbSlqsbE2tzslEoKukO4F+eCWWDsAmZetK2jLDfwawE7AIMEnSIZIGADdLWqEEHbsBz5nZgynO88fArpIOLHrfc9GTnbOlJS0qqW8y+DOAazP/cFH3W7uNI8u1DCbpq8Cf8aHTPunvR8A6ZrZnkRosHWBJE4DuwKfArcCvzex9SfsCN5nZ9KJ0NKBzEfymPAgPkj4E+AvwILCXmf2pIl3nA/8FDAQ6mtkJSetEYBczm1yynkXSMHJ1/Di9ZmZvSxoEnA5MAx4wszNTL3xmm/9wwbSskva5AXAN8B/AAGC4me1b1H4b0HUusBbwGrAEcCdwPHBKNiIqZL/t1ZBlSNoJ6AzckwzHSsC/A+8Cx5jZKwXuu0Pq8RwDLI33vn4L/A7vVUwAbkvG9nOjVyaSVgUOxG+Gt4CZZvaKpBOBr5jZiDK15Y7ZLsAeZpYNTWak1vxc4CMzO6IMPTldmRHrCDyGN4R/BR4CbjazVyUtaWbvpM8XfsyS/3CdNIxD0r3AOWZ2a9GGtJWO2X6rpGWAD4BvAUvi1/4BwNlmdkEhGtqjIcvdDAcB+wGGtxC/ASaY2XRJXbMhZsFaeuAGbBTwc+A+4G7gfuAGMzu2aA1tkYZBPwJeN7PT0rZV8KHKsWY2tcybIu2/I3AZsBTe88p6tcvjQ8rbytST9n0h8MekaWkzGydpd2BLYCbwf6Rrq6D9d0yTH6sA3wQWw3v2LwCfAcOA3c3sO0Xsvw1d+VHHPknX62Z2c6vPrQycDHzfzD5c2DrapY8sGbFOwHbADma2Md772Q74D0nblGTEOpjZNOAIvAVfHLjbzKbi9druzD5XtJbWutLf9YGuuE9sI0l3SVrZzF4EjqzCiCWWBZ4BlgF+IekbAGb2dzO7pQIj1h14BNgW2Bs3HJjZjfiwaTKwfJHugZyz/Gp8RnAffOLjKGAV3A0wIuktdIawFUr7PBr3+fYCxkm6XdKQ3OdWBHoUYcQAMLN2+cC7tc8BI3Lbspt2SMH77pD+dk5/l05/DwNeBG7CexVVHp8uwL3Ay8CRwKZ4T3E80LsCPdkxU+64DcIbgYvwyZGVKjxeXYH+wLHpHJ4NDMi93y3/OxbyvrOR057Aden5o8CuwP8CjwPrVnhsegJP4UUlrgfGAqcALwEnpc8snd0HRTza69BS+I25I7Ay3u2/xcyeKlnHqbjj8xXgDfwktwCrAzea2WvZkKFMXTl9O+HxbQLuwmPIDsMN/b0VaToZP0YDgNG4S2AwsA1wvpm9VLKebEi3OvCquZ91XWAv4KvAk3iIyKdWcE9R0ja4H3MwMNDMDpG0NXAosKuZfVLk/tvQtSawKvAs8Esz2zCF8IwFzivlvqvKihfcQgh3MnbBe2Yn4LM6hwOdCt531rPYH59xWxefEfwJcFxrnRUdn+WArfFwhhZgD7x8OcCSFejJjtl2+ERIN+AfwOq589mt4mvqDnzG+bDctqHpulqiwP3uCXyj1bad8N70ENxP9r20vWOJx2NFvLEZkF4vhvdYr8J7X0cA15Slp930yHIt5zbA9/BgxZuAa/GhwDDgYzO7piQ9JwOPm9l/Jp/UQHw4cqKZ3VeGhja07QVsBvTA/StT8N7YCWb2u4r8YplD/df4mg/rmNkBkrYFhgMjreQeR+6aGoEb/Q648VoCONjM7sl9dqEfM0mj8Yb4+vx1m/y/hwHr4w3zbgtzvw3oagH+E3dLvGJmY3MTbOfhCxa9C4w1s2dLuZ6qaN0Kbikewlv2nsBPgSeA3fGLsOje2LrAaun53viM0ra5928HvlXRccl6PZ1wX8ay+A25Gx6B/QRwRsXnbis80HQSsFjaNh4Pk6lKU3fc17NWbts4PLzgVxTkT8R9co8Da+S2DccLnH4Nd5wDdE9/y+yNjQd+iPfCfo3Pxp+B+587443ikvnrruhHu+mRASQ/xsXAMJsVz7MdPoU/uuB9H4BPxV9uyb8kaSTeavbH/QdbmtmGReqYi7b8FPnFQG/gE+CnZvZc2t4P9wF9WlaPbA7xR0vjwcpr4VP1m+FD3s2K1jIHbVkPowNwPu77eTi91x2frewDPGQFxEZJGosbz2z1MuEzk8/ix+UyvIePlXgTS9oCH+n0Nu+tTsI7D5Nwv/RU81S8Uqk8aXlByW4Gef7ky3jU976SrjGzf+BO9kJXeUn7Hg0cYGZPZLrwE/4UPlT6FA/m/HzIUqSmueg8BpiOJ65fC7yU4tymm9nfss+VYcQySYCl4Nc+wLtmto+kvfGb9QU83KBU5MnYnSTdYx4EOwm4SdJpeODwIXhoyHnAsZIuM1/zdWHyN1LqWgqnaAHONLObJQ3EY/+6WlHhDHOnN264xkrqjRuuQ9P1/jweevEVKzDQfE7U3pDhw6RP8XSNv+DDt6OAFeXpLAPxqfsiOQh40syeyFry1Eq+K2kmnrYyOftwWUasVU+sKz5begI+m3S2ebR85v/5cRma8qQez9p43ffj8FzBkXhP8dqy9eR4Cw9rOEDScviwaSLeEO2ApyIdjE/gTCzAiIH3vI6QtJmZPQhMljQ1vXcg8IKZfVhGBkErbsdnkjcBvgE8kRmudI11KtuIQTuJ7JfUB7jazIak171wJ+k/gffM7L8K3v+ueODt/skROxM+v1HXxScaTrJUIaFM5LWz3kvP98JXie9lZl9P2/4InGZmd1RwUyBpc2BFM/tlatX3xns87wM7VdDjyHR1wuO01sP9rbea2e2alS61EnCJmW1ToIYD8WPxW9wvZ3hy/+n4sHNmFecsaeuJj3Q2AZbHw1BG4bO695U96qitIUsHcoSZnZFa8UPw1uJCM3u9ZC2r4Q7zH5nZH9O2LDfvctz39LMyNSUN2+LD2vFm9l4y8OcCH+I9joG4U73sWa9sNnBDPLygH+5zei756JYCNjGz35SpK2mbzT8or2ixCT4T1wm4Mxm0JYEuZvb3gvVsgBuIXfDsgrfw1LY7q4xBzOlbBT82u+I5sKUVU5xNR40N2cV45P4E/CR3BtbA/WSPAH+0EtKQcnqy1vMePEByJu7nORNYs+zWM91oLWb2lKRT8NijA4H/wVNZlsaHCHeZ2ZtVOPglPZ70rIi36A/jISulD01akz8e8qKTX8Xjth41s/sr0NMLj1d7vux9t9LRrfWQNvVeVwcmm9m0KgxsLQ1Z8n0dg09RbwjsZ2Yvpx7IpnjA57N476zMZOd86zkReBtvPe8ovavtZYY3MbMx6fVI3JA9hQfmTilLSytdWW9sbzy04NjkGhiON0Rv4hUcSjVm6do5DzjUUnxYmrHsZF64cH1gUzM7vwQtrWdzZ7t2yh5Oymv4TZHXqpPlEsKrGtr+C1ZRfM6CPnAj9izuS/lObvtSwPeBLSrU1gsf0lUSuZ80rIdXa7gEWCRtWw44EQ8QPrlCbZ2B3wOvAxvktg/Gh+dV6ToYT1a/gRQPmHvvWWD79Hyhn1dSHFi6bralVcxjVdcSHiu2Gx4OM5mUAdJsj8oFfMGDmwV3/hCvAjAReAD4egVa1Op1x7beL1nbcummPKjV9nWzY1WWPjzU4vjcuVsWjxl7Bl+YZdEKj5NyzxfB3QEvAaenbWOA35ak5Xq8usZleDxdKQGl89C0Hu6ueR5Pvdso994GWUNZ5aN2Q8vc0GRRPOXok7R9LN6iPof30Arzj+U0rIY7qv9guRnJpuluA/JSKscCF5nZDRXq6AlsjEerH4IPbz+WtBaeAzsET/u5o0KNnw/h5PWzzsRjufoAg83TbRa6i0Czgm//DS9ycCs+e/tVPF/3AjN7o4IhZf54rIkfi4F4L+1R3C82wAqsstwwVVvSL9hCbIxXa5iKd3mXStsXB8aUqKOpWk+8N7En7m/ql9u+Ge5U/3ZFuvrjPcNu6flEPBJ8r9xndsQr0lZ17LJGvUP+POLDvFHZewXuvwvuCtg8t20XPKD7WrzaRVXH5lJgu/S8N/BdPLbubmDjoo9NI4/a9Mjkq+bsiFd5vRLPc/tffIawBU8NuqwEHU3ZeiZtm+PBwC/jfo278NnTKXg1kI3x2Kw3StZ1FfCUmZ2V27YLHgT7CnCqpbCVknVtiEfIH2YpsyE7b3PqeRVxTiV1thRQK+mneIjHueZZKUi6Bg+X6YAb1LKCqbPrfFd8aL0tHsc2HI9rewN3o8xoihFIlVZ0PluFHfBYrXHAjcByufd2wsfvo0rS0sytZ1aUcBng2+m4jcerFbyAh2SUqWdLvKLpTnN5/3i8QSo0oX8u++6Wjs3zuFHtOofPFOZDxAOlxzIr8Xtl4OakZQRwIZ6itSjuo+pewTG6DNgIr5Ly83SuniOXzN4Mj9r0yAAkZWkRW+ELP1yFJ+1m8SyFhjg0a+uZ05e1onPqTfTAb8r3ywwFkfR74FXcDSB8oY6HW32mi5nNKENPbp8DmRWAuyZ+k/bBHfylpUelUKLj8Yot38d706PxbALwKP5xeB38cSXqEn6+foBXEJ6M+1knSLoar1RbesDy3KiFIctFyWc5XVvjzuEueGs60cyeLLKLK2kY7t+51Mw+SM7grIWagqdr9MDzLu/FZwU/KELL/FBWoOtc9v0DPGRhR3nd/U3wY/gXPL7urxXpEu7IPxlPj3osbR+G94bewGuzPVSwhh6pYemKxx+OwBelOd7M3k6fWx6Ph9y1KC2tdbW+hyRtilfAfUSeUH+ymRW+EPH80PSGLNfL6Icv1rGNmb2eIp23wx3Zz5nZ2SVoadrW08wsBXU+m4x9ZQYsp2trvPDen9LrFfAA5g3xoe+TeL5iVYZ2MXxVq8dwH+vEtP0kvL78wQXtdzB+/ayB91Czqih98N7PUDx3eFzaXmYPOruWRuDnaBruRsl6r7filV9vrCKCf240vSHLkHQDMMk8t7KTzVqivg+Amb1a4L6bsvVspbEb3kP8sxW4EOoXobVRla+lOQSYYWZXVKfs80mkI/BihfcDV5jZ5NwNXUTl12xt02wBkfF4aMPVqdHeHM+BvatkI5aFFQ3C/XOXA0fjwcvPALcBU8p2AzRE1U66Rh74EPIKoE96nUWqH4mnjRS578H4Sb0PT1/JtvfBAzpfwMvOZNtLq9Q5F60P4w7kThVr+RcneX5bOqelB1IyK4L+e8APc9tXBn6Bx7kNLnD/I4Fr0/PFcJ/qmXjv9EVylWgrPHdX437oofgoaG/cHXA5Fazp0JDmqgXMx8E9He/iDkqvF8f9U30L3u9v8ZCGb+CVQgfji/5mEeqb4yV8SjdiOQ09SDNueKT1aUD/Cs/V+rnnlWU2zEFXNgLphs+8ZWXJtyPFtOGrNRU5U/kM8HR6PhrvAWbvnQgcVdGxGZQal8XT9b0EHr7TL71/DfDdqs/h3B5NW1gx17WXOT/KovflpUOm4iVVXi5Qw0jgLfPh7GJ4DNtH+ApEx0oaZrn4JyvZX2CzhjyHAF+T9DB+o+wAbC9plLmDtsyqG32Be9Lw6aeWlm9rhlij3P73w3teL0s6HK9d9zdJMywlRBeodyBwnqRpeL5wS+69D/HA6lJJ0QCHmdnuwAzgquROeRk4Wl6KajA+kdUU57I1TekjyxmxzvhFtyoeXzMND/BcD6/IWuh6eZKewa//teQr2mxkZtlqzicC08zsjCI1NKCxK17xYzoe1f8ifjOsipcyOq4iXbfgs5QT8MVDpqftld8EaeLoKNyAXYSH8eyEn9+RJWq4EFgb74VNkPQAXmPvhTInayTdja+zermkr+Eji/vkqyVdBXyMh1uMbyYH/2xU3SWcSzc3GzJdhE+RX49HgP8cbxlK8a3g8WAX4Ab0U3Ir5uBO0KsrPj6d8BVr+ufey3xAq+JO2lKChFvtO7sBjsEXzHgJX8+gGa6tzfEsh3VIFVLScXyUCtJtcD/U83gDfVYF+x8JTMj2i7trNmr1mcp8rQ3/jqoFtHGAV8Hjw8AX8TgIj5r/B76qcplastCPKczypTzArMVJS80zY1ZP+tJk7N/DJ0PWaPW5cfjq3GXpygzstcw+AbInPmT5HdVEp2cGdufsmkqvu+K+oOHZcaIinx4ejtG5bA24K+Kp9PxgPE4y//5ZeGn00o/J/Dya1keGt5KnyPPhljKziyRdiRuUp8sUYp6HN1TSUODfJf0KX7yj1CEAzDbsXg+vPLtZikxfHrhfntd4El6n7TZ8qFmGrvyQ4xFg8eRn6Whm1yWNM6yCIOGcrqPx0k9IOhoPebgXT9/KKoMIzyksW+OVSVfZ8X95n103vIdP0nIEsIKZTZ3bl5uFDlULyJMufADM7E9mdivw9/TePvgN+rRVVO7XzO40s9Xw7vjR2eaSNWT7G4pfgLvgAbhD8XUOd8dLIn9mZo+a2T+L1iTpK8Bt8gqi4GsnbI0nZK+SosGH4sP0SpC0LF4csZc8lawnbnC/jrsqPoZSl8KbI2Xv33zFrx/gRv0u4BH5IjXguaA/g8+r5TYtTePszzuBJe2Hd/v/gYc/DMXjWlbCa429WZXOPFX0xkiBuel1Z+Aw3Hf3b5J+ArxtZheX7CzujSfNZwnip+JDyVPxir3vAw9bCWWi20LSGDw85Vkz+7l8KbqLzWzTKnU1E9moAxiAjzoOb4YskXnRTIYsS0U6Bg9OXARflGJL3I8xDR+mVLI8WNWktJb98JbzFjM7J21fFV+rcgbeu1jbfMWkKsoIDcDLB30NL2l0Id6JLGLdx0b0ZMPw3rhz/0ncj/eafKWmm3FDdn3TzsZVhKTvA78yL35Z+UzzvGgaQwYgX4r+N2b2dUln48uonSlpO7yn8WjFEisjl9byCJ7WcjXu35iAz1AuBXxgHjdWZlrLRrhxXQ4PmvwEb833wmcvf29mvyj7ZmiVbnMevhDMcviw9yx8aPlNM7uoLE11pA69MWg+Q9YLOBQPddjKzDZP2x8FTjSz26vUVxUpMHeIme2dAnPfwLMMtsEX8hhmZn+uQNexwBb4zNdKuKP4JjxMRrjB7Wdmp5etLafxeuB2M7sm9RhPwmfhtsyNApq+xxG0TeWzlvI1A3ub2WQzmyrpEbyUyg2S+uNVYd/6shqxxGHMmlTYG4/7ORI+r9QwDC89XBqp0RkBrJPz2W2QdAxNmm6mghnA3JCyGx4ovEjqob0AfEfSBEmrmNmLMNsESlBTKjVk6Wa4COgk6SM8kfdOvNu/Dj58egZPgv4y01Zay3QqSGvBnfu/N68I0gOYbl7Xa4ik6/Bl3u4pW5RmFd88xXwh2VvxSrlPSnob6I6XESp0hfCgXKqeUj0Rz+caiQ+XhqRt95rZIXid8LFVhVs0C62myO9h9inyb+KBr2VPkT8I9JSXVJqWekCLpfceAbbKh9OUyAm4PyxrKG/DK5QcD5yCV7g4ORngjhXoCwqgsh5ZyukaaGabpNfbp7c6A49LOi6bmQucZgnMTQbqBbwhvFDSuWb2dDbExP1m91Ywa5ol+V+cjPr1wHAzO1nSinhlh6lm9ncoP8k/KI7KnP2ShuNDxy3xxRUONLNh6b3BeK3wg76s4RaNUPUUuTzx+UC8jNDbeL7iQGDHKmKz5En+M81soLzM9rpmdkCEVrR/Kp21lLQMXh3zm/h6lBem7TsDPzGzjSoTVyOqnCJP8Vjb4Mn8mwPXAfeb2RMVaOmAh1rsS/IlmpdF74jHj31StqagHJoi/EK+2vQt+KKtI/Gk4wurcBYH9Sf1FC/AS+QcZymPMWi/NIUhg8/9LnviPbQHzGxIxZKCmpN8iWfhrotN8QDrpg/uDOafpjFkGSl/sIelBT2CYEGRdCieilRJqlRQPE1nyIKgKOqSbhPMP2HIgiCoPVUHxAZBECwwYciCIKg9YciCIKg9YciCIKg9YciCIKg9YciCIKg9/w/vv35ABZt+lgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 360x72 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "ESG_plot(50)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Look at the thin bars concentrated on the right end of the plot. These are plants with small capacities and high variable costs. Conversely, plants with larger capacities tend to have lower variable costs. Why might this be the case? Electricity production typically benefits from economies of scale: it is cheaper per unit when producing more units. Perhaps the high fixed cost required for electricity production, such as for equipment and land, is the reason behind this phenomenon. "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    " "
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}