{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Heat Capacity of Cementite ($Fe_3C$)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Bengt Hallstedt, Dejan Djurovic, Jörg von Appen, Richard Dronskowski, Alexey Dick, Fritz Körmann, Tilmann Hickel, Jörg Neugebauer, Thermodynamic properties of cementite, Calphad, Volume 34, Issue 1, March 2010, Pages 129-133, ISSN 0364-5916, http://dx.doi.org/10.1016/j.calphad.2010.01.004. (http://www.sciencedirect.com/science/article/pii/S0364591610000052)\n",
    "\n",
    "The TDB file used here differs slightly from the published TDB to ensure compatibility with pycalphad's TDB parser. All phases except cementite are omitted. The numerical results should be the same."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2022-02-19T19:16:53.923299Z",
     "iopub.status.busy": "2022-02-19T19:16:53.923299Z",
     "iopub.status.idle": "2022-02-19T19:16:53.933979Z",
     "shell.execute_reply": "2022-02-19T19:16:53.932979Z"
    }
   },
   "outputs": [],
   "source": [
    "TDB = \"\"\"\n",
    " ELEMENT C    GRAPHITE                   12.011     1054.0      5.7423 ! \n",
    " ELEMENT FE   BCC_A2                     55.847     4489.0     27.2797 ! \n",
    " TYPE_DEFINITION % SEQ * !\n",
    " TYPE_DEFINITION A GES AMEND_PHASE_DESCRIPTION @ MAGNETIC -3 0.28 !\n",
    " PHASE CEMENTITE_D011 %A 2 3 1 !\n",
    " CONSTITUENT CEMENTITE_D011 : FE : C : !\n",
    " PARAMETER  G(CEMENTITE_D011,FE:C;0) 0.01 +GFECEM; 6000 N !\n",
    " PARAMETER  TC(CEMENTITE_D011,FE:C;0) 0.01 485.00; 6000 N !\n",
    " PARAMETER  BMAGN(CEMENTITE_D011,FE:C;0) 0.01 1.008; 6000 N !\n",
    " FUNCTION GFECEM      0.01  +11369.937746-5.641259263*T-8.333E-6*T**4;\n",
    "        43.00  Y  +11622.647246-59.537709263*T+15.74232*T*LN(T)\n",
    "       -0.27565*T**2;\n",
    "       163.00  Y  -10195.860754+690.949887637*T-118.47637*T*LN(T)\n",
    "                  -0.0007*T**2+590527*T**(-1);\n",
    "      6000.00  N !\n",
    "\"\"\""
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Do some initial setup, including reading the database."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2022-02-19T19:16:53.938980Z",
     "iopub.status.busy": "2022-02-19T19:16:53.937980Z",
     "iopub.status.idle": "2022-02-19T19:16:54.901699Z",
     "shell.execute_reply": "2022-02-19T19:16:54.901699Z"
    }
   },
   "outputs": [],
   "source": [
    "# Only needed in a Jupyter Notebook\n",
    "%matplotlib inline\n",
    "# Optional plot styling\n",
    "import matplotlib\n",
    "matplotlib.style.use('fivethirtyeight')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2022-02-19T19:16:54.901699Z",
     "iopub.status.busy": "2022-02-19T19:16:54.901699Z",
     "iopub.status.idle": "2022-02-19T19:16:56.661508Z",
     "shell.execute_reply": "2022-02-19T19:16:56.661508Z"
    }
   },
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "from pycalphad import Database, calculate\n",
    "\n",
    "db = Database(TDB)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Compute the molar heat capacity at all temperatures from 1K to 2000K with a step size of 0.1K.\n",
    "\n",
    "We do this with the `calculate` routine instead of `equilibrium` because the cementite phase has zero internal degrees of freedom. Since there's nothing to minimize, we can do the computation faster with `calculate`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2022-02-19T19:16:56.721543Z",
     "iopub.status.busy": "2022-02-19T19:16:56.661508Z",
     "iopub.status.idle": "2022-02-19T19:16:56.741641Z",
     "shell.execute_reply": "2022-02-19T19:16:56.741641Z"
    }
   },
   "outputs": [],
   "source": [
    "result = calculate(db, ['FE', 'C'], 'CEMENTITE_D011', T=(1, 2000, 0.1), P=101325, N=1, output='heat_capacity')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2022-02-19T19:16:56.761125Z",
     "iopub.status.busy": "2022-02-19T19:16:56.761125Z",
     "iopub.status.idle": "2022-02-19T19:16:56.981139Z",
     "shell.execute_reply": "2022-02-19T19:16:56.981139Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnYAAAGTCAYAAAClNV9VAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABj5ElEQVR4nO3dd1gU59oG8HtYioggiICNooKKiphgiyW2RMWuaNSY2EvUL1FzbESjxsRgi8aCaIItJyaxR2NNEWPsJkcBG2JUxIYILkiH3fn+QBaGOsAuC8v9uy4vmXdmZ593gd2HtwpKpVIEEREREVV4RvoOgIiIiIi0g4kdERERkYFgYkdERERkIJjYERERERkIJnZEREREBoKJHREREZGBYGJHREREZCCY2OlAeHi4vkPQK9af9a/MWH/WvzJj/fVffyZ2RERERAaCiR0RERGRgWBiR0RERGQgmNgRERERGQgmdkREREQGgokdERERkYFgYkdERERkIIzlXiiKIq5fv46wsDDExMRAEATY2tqiUaNGaNasGQRB0GWcRERERFSEIhO706dPY+fOnTh27BgSEhIgiqLkvCAIqFatGnr16oWRI0eic+fOOguWiIiIiApWYGL3+++/Y+nSpbh69Src3d3x/vvvo2XLlnBxcYG1tTVEUYRSqURERASuXr2KoKAgDBw4EJ6envj000/RvXv3sqwHERERUaVXYGI3atQovP/++9i0aRMaN25c4A3atGmDoUOHAgDCwsKwZcsWjBo1Co8ePdJ+tERERERUoAITu9DQUNja2hbrZo0bN8aKFSswd+7cUgdGVN5FvMzAxusJqGYiYLqHJaxMOReJiIj0q8DEztbWFklJSahatarsm/3666/o0aNHsRNCooomJUPEoBPPcfelCgDwNFkN/442eo6KiIgqu0KbGIYMGYKkpCRZN9q3bx/ee+89rQRFVN5tvpmgSeoAIOhRih6jISIiylRoYnf58mUMHTq0yORu+/btmDRpEho1aqTV4IjKI5VaxKrgl5Kyx0lqPUVDRESUrdDELjAwEJcuXSo0uVu7di1mzpwJLy8vHD58WCdBEpUnhx+k4GW6WPSFREREZazQxG7AgAHYsmVLgcndkiVLsHjxYnTp0gU///wzrK2tdRkrkd6JoojVIS/zlNtV4cQJIiLSvyI/jfr374+tW7dKumVFUcSsWbOwZs0a9O3bF7t27SrWJAuiiurk41QEx6TnKWf7HRERlQeythTr168ftm7dinHjxmHIkCGoW7cu9u7dixEjRmDDhg0wMmJrBRk+URSx8mre1rrMc2UcDBERUT5kZ2R9+/bF9u3b8c8//2Dfvn344IMPsHHjRiZ1VGmceJiCC8/S8j2nZpsdERGVA4W22GXtKJFTjRo1EBcXhzt37uQ5LwgCdu/erd0IicoBlVrEkr/jCzzPFjsiIioPCk3sbt26BUEQJGUmJiaoWbMmwsLC8lyf+1oiQ7H3XjJuKDMKPM+8joiIyoNCE7vQ0NCyioOo3ErOELH0f9LWup71zHDiYarmmIkdERGVBxwgR1SEddde4kFC9i4TxgLwyetWkmvYFUtEROVBiRO7+Ph4TJs2Dbdv39ZmPETlSsTLDKzJtW7d+CYWqG8pbexmYkdEROVBiRO7lJQU/Pjjj3jy5Ik24yEqV+ZfikNKdmMdalYxgu9rVsg9nJR5HRERlQel6ooV2UxBBuxEZAoOP0iRlC1uZQVrM6M8vzj8TSAiovKAY+yI8hGXpsbMcy8kZa3sTPCua+YOK7lb7NT8I4eIiMqBEid2pqam6NChA/eHJYO04FIcHiepNccKAVjVzhpGrzI6AdLMjmkdERGVB7K2FMuPtbU1Dh8+rM1YiMqFk49S8N/wJEnZR82roWVNU81x7hUb2WBHRETlAbtiiXKITVHh/85Iu2AbVzfG3JbS5U2MOHmCiIjKoWK12N24cQObN2/G1atXER8fD7VaLTkvCAKuXr2qzfiIyowoiph6RinpgjUSgA0dbVDFWJrJ5R1jVxYREhERFU52i9358+fRrVs3HDt2DLVq1cL9+/fh4uKC2rVrIzIyEhYWFmjfvr0uYyXSqc03E3E8UjoL9qPm1dDa3jTPtXm6YnUYFxERkVyyE7ulS5fC0dERly9fxsaNGwEAH3/8MY4fP45jx47h0aNHGDJkiM4CJdKl4Jg0LLwcJylrZWeC+bl2mMjCMXZERFQeyU7srl69ivfffx/Vq1eHkVHmw7K6Ytu2bYvRo0dj6dKluomSSIdiUlR472Qs0nKMLLAyERDYuQZMcg+me4Vj7IiIqDySndgJgoDq1asDAKpWzVzLKzY2VnPe1dUVN2/e1HJ4RLqVrhYxJigWkTn2ggWAdR1s4GJZ8BBUIfcgO3DBbiIi0j/ZiZ2TkxPu378PADAzM4OzszOCgoI058+dO4caNWpoPUAiXVpwKQ5/PU2TlE1sYoGB9c2LfS+mdUREpG+yE7tu3brhwIEDmuPRo0dj586d6N+/P/r164ddu3Zh6NChOgmSSBe+u52IzTcTJWXtHUzxZdvqsh7PcXZERFTeyF7u5D//+Q98fHyQnp4OExMTzJgxAyqVCgcPHoRCocC8efPw8ccf6zJWIq35/WEKZp5TSsrqWSiwo2vB4+pyMxIAVY5kTg1Aob0QiYiIik12YmdtbY2WLVtqjgVBwKxZszBr1ixdxEWkM1efp2F0UKwkKauiAL7vVgN25vJTs9yJ3VfBLzHvtfxn0RIRkfaJogiVCGSogQxRlP6vFpEhZv+frgZUOcrS1YDq1bXpr8pVOb7OUIuF3leV65xKDXSpIsBNz69JibcU04azZ89i/fr1CA4OxpMnT+Dv74+RI0cCANLT0/HFF1/gt99+w/3792FpaYlOnTph0aJFcHR01NwjNTUVCxYswL59+5CSkoI333wTX331FerWrauvalE5dv9lBt75PQaJGdkZmQBg85s1JFuGyeFRwwT/e56uOV529SWSM0S8K68nl4hIJ0QxM2nJSjpU4qtkJZ/kJF0tahKUdDVeXZ872RGR/irRkd7r1bkc930WYwIrZVyepChdLU3AciZZ2eey7iU9lxVTujr/GMoT1yb639CrwMRu+fLlxb6ZIAiYM2eO7OsTExPRtGlTjBgxAh988IHkXFJSEoKDgzFr1ix4eHggPj4eCxYswJAhQ3D27FkYG2eG7uvri6NHj2LLli2wsbHB/PnzMWzYMPz5559QKNgxRtmik1UY8msMniVLd0zxa1sdA1yKP1liQ0cbeB+NRlxa9jvL2msJ+F8NU2xzUqFmFf78EVUEWclJeo4kJD0rscmRbGR+nZ1gpOe4NuNV8vMwSgFbdWK+12nuLRZyLse98ivPyK/8VWtRViKm351wTICHCfoMQK9U5SDRFJRKZb5h2NjYFP9mgiBZAqU46tatixUrVmha7PJz69YttGvXDmfPnkWzZs0QFxcHV1dX+Pv745133gEAPHz4EB4eHti7dy+6d+9eolhKKzw8HG5u+m6M1Z/yWP8XqWr0ORaNGy8yJOUfNq+Gz1uXvIntYlQqBv8qbQEEAHtzIyz0ssKIhlWhkDlmz1CUx+9/Waqs9Ve/SlZu3r4Dp/oNkK4G0l4lHmlqEWkqUVKW/qosLcfXmvL8zmc97lWZJLEScydJ+SVTBSdS5eCzmAzE541S8WGHBnqNocAWuxcvXhR0Sm9evnwJIHO8H5C5aHJ6ejq6deumuaZevXpo3LgxLl68qLfEjsqXuDQ1Bv/6PE9S51PfHJ+1Kt2YuLYOZjjsXRPv/BaD6JTslsBnyWr83xklNl5LwNTm1TCkftU8+80SFYcoZrbGpKpEpKpEpKiANJWIlFfHmjK1tCxVlfmYdHXOxEpEmirH11lJTq6ytDyPyy5Le1WW8aos+2+bqsCFp/p8qaiSMRIAYwEwMRKgMAKMBQHGRoCJkHUsPWdiBBgbCVAImf+bvLpGofk681oTIwHGr67J+l8hACZGmddm3TfrscYCUD9d/z/7eh1jVxxpaWlYsGABevXqpRk/9+zZMygUCtja2kqutbOzw7Nnzwq8V3h4uE5jLavnKM/KS/2TVcCH180QHC/tFm1nrcLHtWPw752YUj9HNQCBzQXMDzPFtZfS57mhzMD/nVHikwsv0NVWhe41M/CalRqG3ktbXr7/2qIWkZnIZP2vFpCmBlLVQLooZP6vBlLVAtLUChx9dvdVWea5zMcKrx6b/fjc98s8znHdq+Os+6vzLLJDlD+FIGYmIwKg0PwTJcfGmv/zL1cIYq7rsr4Wc12X9z7GucoVee5ReDx57yF93pznylWniInu3/+K6hGoEIldRkYGJk2ahLi4OPz444+lvp+uu0kqa1dMlvJS/+QMEcN+j0FwfKqkvGMtU+x+2xZVjbU3yNUNwKlmIr4KeYmvrsYjTZS+08RnCDgYZYyDUcYwNQJa2ZniDQdTeNQwRVMbYzSwMoZxuXp3Kjltfv+zBoHnboFKUYmalqm0Vy1V2V9nnUOOx2Rdm+OxGa/OqQu+LvN8Zgxk2EyyWmhetfRktepIynMeC9nHmddlfp2U8BK21lYFXJP3XprWpfyuE7K/LugemlajHPdSCPnvjlMWysv7v76Uh/rLTuxatGhR5A+KIAi4evVqaWOSyMjIwPjx43Hjxg0cPnxYsruFvb09VCoVYmJiULNmTU15dHQ03njjDa3GQRVLYroaI/6Ixekn0qSujZ0pfnxLu0ldFmMjAXNbWuENxTN8H1sDe+8l5zuIOU0NnItKw7mo7B0vzBSZ6+g5VjNGPQsF6lkoYFvFCDXMMv/ZmBnB2swIVY0FVFEIMDcWZK+3V1xqMTOxSX01xin11ZimVM3XmclVujr/cw+fGqN68kvN4yVdhuq83YjZ3YnIcT77sVRxmBgBxhBhZmwEU4UAUyNpYpKzzNRIgIlCgImQWZ5VlvPrrMdKzudOjIzySX6EfBKoHMmXiaZLLfNYm4lQeHgM3NyKP0adSFtkJ3YdOnTI84OvUqkQGRmJixcvwt3dHS1atNBqcOnp6Rg3bhxu3ryJw4cPw8HBQXK+ZcuWMDExQVBQkGbXi0ePHiEsLAxt27bVaixUccSnqTHs9xicj5JuFeZpa4Ldb9vC0kS309FrVxHxTeca+OT1DHx7MxE/3EnEi9TCh2enqoB/41X4N15+JqMQAHOFgCrGmX+1GwmA0avfUaOsY2SWqV+t9ZT5L3sJBNWr2XQq8dXMOq3MqDMF7sWX9iaUDyMBqKIQYKbI/N/UKDPRN1UIqKIAzBRC9j+j7OtMXh2bKjITnMwEKStpevX/qzLjV2WmrxIhU6OsxEuacEmSNEXmz6AgCOWixYKoMpOd2AUEBBR4LjQ0FD4+PpqZqXIlJCTg7t27AAC1Wo2HDx8iJCQENjY2qF27NkaPHo0rV67gxx9/hCAIiIqKAgBYWVnB3Nwc1atXx/vvv49FixbBzs5Os9xJs2bN0KVLl2LFQoZBmaqGz6/P8U+O9eUAoKmNMQ70sIW1WdmtMeRiaYylbapjcSsr/PUkFT/fT8afj1MRkaCdZiiVCCRkiEgobws5GagqCrxKoKSJU1YilZVgZSQlwLa6paRM7mOrvEqoqrxqmc1M3LKvM5TueiLSHa2MsfPw8MCYMWOwaNEi/Pnnn7Ifd+XKFfTr109z7OfnBz8/P4wYMQLz5s3D0aNHASBPkpZzIWM/Pz8oFAqMHTtWs0Dxpk2buIZdJRSTosLAEzEIjZUmdS1qmOBAT1vU0NOMBRMjAd3qVkG3ulUAAA8SMnDuaRpCY9Nx/UU6brxIz7O2HmUzMcrdOoXsBEmTKGUmP1WMs68zy92ClTNx0iRS0qSroMeaGsnvqgsPj4WbW42iLyQi0gGtTZ6wt7dHWFhYsR7TqVMnKJXKAs8Xdi6LmZkZVq5ciZUrVxbrucmwRCWpMPDEc9xUSpc0aWVngr1v1yzTlrqiOFUzhpOrMYbnKFOmqvEwUYWHiRmITFDhSZIKL1LViE1VIzYl8/+4tMyxZ8kZIpJVok4XITV9lSiZvmpZMtV05WV3zWUlVKaa5Cfz2qT4ODjY2mS3UBkhV3JVRHImKc/uXiYioqJpJbGLjY3Ff//7X9SpU0cbtyMqlkeJKgw4/hx34qVJ3RsOpmUypk4brF9Njmhew0TW9VmzRZNVmbM71cgcG6cWMxM+EYD4arycShQz13kyEmAkAAoha92lzFl0RlnLCbwaUF7ageTh4c/h5sZ91YiI9EF2YpezyzSnuLg4hIeHIy0tDZs3b9ZaYERyRLzMQP/jz/OMW+tSxww7u9WARQVI6kpCEF4NcFcIqF68LW6JiMiAyU7s1Gp1nr/iBUGAs7MzunTpgvfeew+NGjXSeoBEBfk3LjOpe5QkTep61DPDd11tudMDERFVOrITuyNHjugyDqJiuaVMx4DjzxGVa9JBX6cq2NqlBkwVTOqIiKjyqRA7TxDlFBKThkEnYhCTKk3qfOqbY9ObNjpbuJeIiKi8K1ZiJ4oiTp8+jYiICCiVSoiidFqeIAj46KOPtBogUU7/RKdh8K/PEZcm/dl717Uq1newhoJJHRERVWKyE7urV69i7NixiIiIyJPQZWFiR7p0PioV7/wWg5fp0p+/cY0tsOqN6lwWg4iIKj3Zid306dMRGxuLNWvWwMvLC1ZWVrqMi0jiz8cpGPFHLJJy7bIwtZkFlraurrcNr4mIiMoT2YldWFgY5s+fj9GjR+syHqI8fo1MwftBMUjNtRPXrBaWmP+6JZM6IiKiV2Qndg0bNiywC5ZIVw7dT8b4P2ORnmvHrQWvW2GWp6V+giIiIiqnZK/e6uvri2+//RaRkZG6jIdI48C9JIw9lTepW9qmOpM6IiKifMhusevbty9SUlLQpk0bdOrUCXXq1IFCId1UXRAErFq1SutBUuVz8H4yJvz5AqpcjcSr37DGuCYW+gmKiIionJOd2J0+fRozZ85ESkoKfvvtt3yvYWJH2nA4IhnjT8VKkjojAdjQwRrvujGpIyIiKojsxG7evHmwsrLCd999x1mxpDPHHiRj7KlYZORK6jZ1ssE7DavqLzAiIqIKQPYYu3v37uGjjz5C165dmdSRTvwamYLRQdIxdQIA/45M6oiIiOSQndg1adIEL1680GUsVImdeZqK94NikJZrosT6jtYY4cqkjoiISA7Zid3nn3+OHTt24OLFi7qMhyqh4Jg0jPg97zp1a9tb4z2OqSMiIpJN9hi7tWvXwsLCAt7e3nB1dUW9evXynRW7e/durQdJhutufAaG/Jp3m7Cv3qiO0Y2Z1BERERWH7MTu1q1bEAQB9erVQ0pKCu7cuZPnGu4AQMXxNEmFQSeeIzpF2v+62MsK45tU01NUREREFZfsxC40NFSXcVAl8zJdDZ9fnyMiQdr/Oq1ZNUz3YFJHRERUErLG2CUlJaFfv374/vvvdR0PVQIZahHjT8Xi+osMSfnwhub4vLUVW36JiIhKSFZiV7VqVQQHB0OlUhV9MVER5l+Kw68PUyVlPR2rYH1HGxgxqSMiIiox2bNi27dvj3PnzukyFqoEvr2ZgM03EyVlLW1NsK2LDUyMmNQRERGVhuzEbsWKFfjnn3/w6aef4v79+1Cr1UU/iCiHoEcpmHsxTlJWt6oCP75li6rGsn8UiYiIqACyJ0+0adMGoijC398f/v7+MDIygomJieQaQRDw+PFjrQdJFd+DhAyM//MF1DlWNbEwFvDT27aoXVVR8AOJiIhINtmJ3aBBgzionUokJUPE6KBYxKZmt/IKAAI728CjhknBDyQiIqJikZ3YBQQE6DIOMmDzLipx5Xm6pGzB61bwdjLXU0RERESGiQObSKd+CE/E9ttJkjJvxyqY2YJr1REREWlbsRK7Fy9eYOHChWjXrh3q1KmDOnXqoF27dli8eDFevHihqxipgrobn4HZF6STJepbKhDQicuaEBER6YLsxO7hw4fo1KkT1q9fD3Nzc/Tr1w/9+vVD1apVsXbtWnTq1AkPHz7UZaxUgaSrRUz8MxaJGdmzJcwVAv7bzRbWZmwoJiIi0gXZY+wWL16MuLg4/PLLL+jYsaPk3Llz5zB8+HB89tln+Pbbb7UeJFU8y6++xD+5xtUta1sdzTlZgoiISGdkN52cPHkSkydPzpPUAZmLF0+aNAl//PGHVoOjiulqnBFWh7yUlPVxqoJRjarqKSIiIqLKQXZil5ycjJo1axZ4vmbNmkhOTtZKUFRxJWeIWBJuKlmvrpa5EdZ1sOZyOURERDomO7Fr0qQJ9uzZg9TU1Dzn0tLSsHv3bri7u2s1OKp4ll+NR2SK9MdqYycb2FbhIsRERES6JnuM3YwZMzBmzBh07doV48aNg6urKwAgPDwc27Ztw61bt7Bjxw6dBUrl39XnaVh/LUFSNr6JBbrVraKniIiIiCoX2YndgAEDsGnTJixcuBCzZ8/WdKuJogh7e3sEBASgX79+OguUyrd0tYgPzyqhytEFW7eqAou8rPQXFBERUSVTYGIXHx8PCwsLKBTZXWjDhg2Dj48Prly5gsjISACAo6MjXnvtNRgby84RyQBtvpGA0FjpLNiv2leHlSmXNiEiIiorBX7quri4YP/+/ZrjadOm4e+//4axsTFat26NwYMHY/DgwWjdunWJk7qzZ89i+PDhcHd3h7W1NXbu3Ck5L4oi/Pz80KRJE9SqVQt9+vTBzZs3JdcolUpMmjQJTk5OcHJywqRJk6BUKksUD5VMVJIKy69KZ8EOaWCOXo7cMoyIiKgsFZjYmZmZSSZK/PDDD7h3755WnzwxMRFNmzbFsmXLYG6eNwlYu3Yt/P39sXz5cpw8eRJ2dnYYNGgQXr7MTiImTJiAkJAQ7N27F3v37kVISAgmT56s1TipcEv+F4+X6dl9sBYKEV+2qa7HiIiIiCqnApvaGjVqhK1bt6JmzZqwtLQEAISFheHs2bOF3rBDhw6yn7xHjx7o0aMHAGDq1KmSc6IoIiAgADNmzMCAAQMAAAEBAXBzc8PevXsxduxYhIWF4ffff8fx48fRpk0bAMCaNWvg7e2N8PBwuLm5yY6FSuaf6DTsDJfuBTvJKR325pwFS0REVNYKTOyWLFmCMWPGYMSIERAEAYIgYPXq1Vi9enW+14uiCEEQEBsbq5XAIiIiEBUVhW7dumnKzM3N0b59e1y8eBFjx47FpUuXUK1aNbRt21ZzTbt27WBhYYGLFy8ysdMxURQx76JSUta4ujHeqZ2U/wOIiIhIpwpM7Dp37ozQ0FBcu3YNz549w+jRozF58mS88cYbZRJYVFQUAMDOzk5SbmdnhydPngAAnj17BltbW8nCt4IgoGbNmnj27FmB9w4PD9dBxGX/HPp28rkCl6PNJGUf1kuAsVHlqH9hWH/WvzJj/Vn/ykzX9S+q0arQWQ/VqlVDu3btAGR2sXp7e6Nz587ai05PdN2SVxm6gVVqEe+FPgOQoSnr7VQF77WpWynqXxjWn/Vn/Vn/yor113/9ZU9nPXz4sC7jyMPBwQEAEB0dDUdHR015dHQ07O3tAQD29vaIiYnRdAMDmd2Dz58/11xDuvHTv0kIi8tO6owEYCHXrCMiItKrEi8y9uzZM9SoUQN//vmnNuPRcHZ2hoODA4KCgjRlKSkpOH/+vGZMXZs2bZCQkIBLly5prrl06RISExMl4+5Iu1JVIvyuSJc3GdawKppYm+gpIiIiIgKK0WKXH1EUi76oEAkJCbh79y4AQK1W4+HDhwgJCYGNjQ0cHR0xZcoUrF69Gm5ubnB1dcWqVatgYWGBIUOGAAAaN26Mt956CzNnzsTXX38NAJg5cyZ69uyp96ZQQ/bd7UQ8TFRpjk2NgHktLfUYEREREQGlTOxK68qVK5JtyPz8/ODn54cRI0YgICAA06dPR3JyMmbPng2lUgkvLy/s379fs/wKAAQGBmLOnDnw8fEBAHh7e2PFihVlXpfKIl0tYm2odD/YsY0t4GzJnUeIiIj0Ta+fxp06dSp0lwhBEODr6wtfX98Cr7G2tsY333yjg+goP3vvJkta68wUwMct2FpHRERUHpR4jJ2VlRX8/f3h7u6uzXioHFOLItaESMfWvedmAYeqXIyYiIioPChxi12VKlXw7rvvajMWKucOR6Tgdo6ZsAoB+LB5NT1GRERERDkVmNhFRkaW6IY5lyYhw7I2VNpa59PAHC4cW0dERFRuFPip3KJFC8mODnJpa0sxKl/+jk7DP8/TJWUzPTi2joiIqDwpMLHbsGFDiRI7Mkybb0hnwr5d1wzuNly3joiIqDwpMLEbOXJkWcZB5VhUkgo/30+WlE1qyrF1RERE5U2JBkiJooiYmBgAgK2tLVv2DNy2sESkq7OPG1op0L2umf4CIiIionwVa7mTu3fvYsyYMXByckKjRo3QqFEjODk5Ydy4cZodJMiwpKlEbAtLlJRNdK8GIybzRERE5Y7sFrubN2+iZ8+eSElJgbe3Nxo1agQAuH37No4cOYKTJ0/i2LFjXNfOwJx4mIKo5OzmumrGAt51rarHiIiIiKggshO7xYsXo2rVqjh16hQaNGggOXfv3j14e3vjs88+w08//aT1IEl/doYnSY6HNjSHlWmJ17UmIiIiHZL9CX3+/HlMmDAhT1IHAPXr18f48eNx7tw5rQZH+vU0SYXfHqZIyt53s9BTNERERFQU2YmdSqWCmVnBA+arVKkClUpV4HmqeHb9mwSVmH3c1NoYr9XkEidERETllezEztPTE9999x2USmWec0qlEt999x1atmypxdBIn0RRzNMNO7KRBWdAExERlWOyx9h98sknGDRoEFq1aoV3330Xrq6uAIDw8HD89NNPiIuLw9dff62rOKmM/R2dLtkX1lgAhjU012NEREREVBTZiV3Hjh2xb98+LFiwAOvXr5ec8/T0xNatW9GhQwetB0j6se+etLWup2MV1Kyi0FM0REREJEexFih+8803cfr0aURFRSEyMhIA4OjoCAcHB50ER/qhUov4+Z50p4l3GnKJEyIiovKuRDtPODg4MJkzYOefpeFprrXretSroseIiIiISI5iJXaiKOL06dOIiIiAUqmEKIqS84Ig4KOPPtJqgFT2DuRqrevtVAXmxpw0QUREVN7JTuyuXr2KsWPHIiIiIk9Cl4WJXcWXoRZx8L40sRtUn5MmiIiIKgLZid306dMRGxuLNWvWwMvLC1ZWVrqMi/TkzNNUPE/J7oatbiqgW112wxIREVUEshO7sLAwzJ8/H6NHj9ZlPKRnRyKkO030dTaHmYLdsERERBWB7AWKGzZsWGAXLBkGURRxLFKa2PVzZmsdERFRRSE7sfP19cW3336rWeaEDE9obDoeJmZvC2euENC5NhM7IiKiikJ2V2zfvn2RkpKCNm3aoFOnTqhTpw4UCumCtYIgYNWqVVoPkspG7ta6LnXMOBuWiIioApGd2J0+fRozZ85ESkoKfvvtt3yvYWJXsR17IE3svJ3YWkdERFSRyE7s5s2bBysrK3z33XecFWuAHieqcDUmXXMsAOjlyMSOiIioIpE9xu7evXv46KOP0LVrVyZ1BujXh9LWOi87E9ibc29YIiKiikR2YtekSRO8ePFCl7GQHp18JE3senILMSIiogpHdmL3+eefY8eOHbh48aIu4yE9UKlFnH6SKinrzkWJiYiIKhzZY+zWrl0LCwsLeHt7w9XVFfXq1ct3Vuzu3bu1HiTp1tWYdCjTstcotDET4GlroseIiIiIqCRkJ3a3bt2CIAioV68eUlJScOfOnTzXCAKXxqiIgh5LW+s6164ChRG/l0RERBWN7MQuNDRUl3GQHgU9lo6v61rHTE+REBERUWkUOsZu5cqVTOgMXEK6GpeepUnKujCxIyIiqpAKTey2bduGzp07o2nTppg+fTqOHj2KpKSksoqNysC5p2lIV2cfN7BUwNlSdkMuERERlSOFfoLfuHEDwcHB+PXXX/Hbb7/h+++/h4mJCdq3b4+ePXuiR48eqF+/flnFSjpw5ql0fF1XzoYlIiKqsIpc7sTT0xOzZ8/Gr7/+ivDwcHz99dewsbGBn58fvLy80Lp1a8yfPx+nT59GRkZGWcRMWnQ+SprYdaxlqqdIiIiIqLRkr2MHADVq1MDw4cOxZcsW/Pvvv/jll1/Qp08fBAUFYcCAAWjQoAFGjRqFS5cu6Spe0qLEdDWuPE+XlL3hwPF1REREFVWxErucFAoFOnTogMWLF+PcuXMICQnBwoULkZqaqrVFjFUqFb744gu0aNECDg4OaNGiBb744gtJy6AoivDz80OTJk1Qq1Yt9OnTBzdv3tTK8xu6v6PTkZG9fB3qWypQqyq3ESMiIqqotDZK3tHRERMmTMCECRO0dUt8/fXXCAwMREBAAJo2bYrr169jypQpMDU1xZw5cwBkLpzs7+8Pf39/uLm5YcWKFRg0aBAuX74MS0tLrcViiHJ3w7avxdY6IiKiiqzQxO7AgQPFupmRkRGsrKzQuHFj1KlTp1SBAcClS5fQq1cveHt7AwCcnZ3h7e2Nf/75B0Bma11AQABmzJiBAQMGAAACAgLg5uaGvXv3YuzYsaWOwZCdj5Iuc/KGA8fXERERVWSFJnbjxo2DIAgQRbGwy/IQBAFvvvkmAgMDUbNmzRIH165dO2zZsgW3b99Go0aNcOvWLfz111+YOXMmACAiIgJRUVHo1q2b5jHm5uZo3749Ll68yMSuEOlqEZejpYlde46vIyIiqtAKTex++eWXYt1MFEUkJCTgn3/+wYYNGzBnzhxs3bq1xMHNmDEDCQkJaNu2LRQKBTIyMjBr1ixNd29UVBQAwM7OTvI4Ozs7PHnypMD7hoeHlzgmucriOUrj2ksjJGVkL21iayIi4+k9hEdp5/7lvf66xvqz/pUZ68/6V2a6rr+bm1uh5wtN7Dp27FiiJ/X29oYoiggMDCzR47Ps378fP/30EwIDA9GkSROEhoZi3rx5cHJywqhRo0p836JelNIKDw/X+XOU1skbCQDiNMed6lZFo0b1tHLvilB/XWL9WX/Wn/WvrFh//ddfZ1sMdOzYsdTLnixcuBD/93//Bx8fHwBAs2bNEBkZiTVr1mDUqFFwcHAAAERHR8PR0VHzuOjoaNjb25fquQ1dcIx0mROvmiZ6ioSIiIi0pcDlTr788ksolcpi31CpVOLLL79E165di92Vm1tSUhIUCunyGwqFAmp15h5Yzs7OcHBwQFBQkOZ8SkoKzp8/j7Zt25bquQ1dcIx0fF0LWyZ2REREFV2Bid2RI0fQvHlzTJs2Db///jtSU1MLuhSpqan47bffMHXqVHh4eODYsWNaCa5Xr174+uuvceLECUREROCXX36Bv78/+vbtCyBzksaUKVOwdu1aHDp0CDdu3MDUqVNhYWGBIUOGaCUGQ5SSIeKWUrpLSAtbzoglIiKq6Arsij179iz27NmD9evX44cffoCxsTEaN24MFxcXWFtbQxRFKJVKREREICwsDBkZGfD09MTXX3+t6TotrRUrVmDp0qX4z3/+g+fPn8PBwQGjR4/WrGEHANOnT0dycjJmz54NpVIJLy8v7N+/n2vYFeLGi3Sockx0dqqmgI1ZideqJiIionKi0DF2Q4cOxdChQxEcHIwjR47g8uXLuHr1KmJjYwFkbjHWuHFj9O/fH71790bz5s21GpylpSWWLVuGZcuWFXiNIAjw9fWFr6+vVp/bkIXESsfXebIbloiIyCDImjzh6ekJT09PXcdCZST3+DpPdsMSEREZBPa/VUK5Z8SyxY6IiMgwMLGrZNLVIq6/kCZ2LWowsSMiIjIETOwqmTBlBlJV2ce1zI3gUFVR8AOIiIiowmBiV8nkHV/H1joiIiJDwcSukgnJNb6O69cREREZDiZ2lQyXOiEiIjJcshO7zz77DLdv39ZlLKRjalHM02LHxI6IiMhwyE7s/P390a5dO3Tp0gWbNm1CdHS0LuMiHfg3PgOJGdlbTtiYCahnwYkTREREhkJ2YhcWFoYVK1bAzMwMvr6+aNq0Kd555x3s37+/0H1kqfzIu36dKQRB0FM0REREpG2yEzsbGxtMmDABJ06cwJUrVzBr1izcu3cP48ePh5ubGz788EP89ddfuoyVSilPYsf164iIiAxKiSZPuLi4YO7cubh8+TL++OMPdO3aFTt37sSAAQPQvHlzfPHFF3jy5Im2Y6VS4o4TREREhq3Es2ITExPx008/4fPPP8fhw4dhbGyMXr16oVWrVli/fj28vLxw8OBBbcZKpSCKIkK4RywREZFBMy7OxWq1GidPnsSuXbtw7NgxJCYmwtPTE19++SWGDBkCW1tbAEB0dDTGjRuHBQsWYMCAAToJnIrnQYIKyrTsiROWJgLqW3HiBBERkSGRndjNmzcPBw4cQHR0NBwcHDB+/HgMHz4c7u7uea61s7PDe++9hw8++ECrwVLJ5e6GbV7DBEacOEFERGRQZCd23333Hfr06YPhw4eja9euMDIqvBe3Xbt28Pf3L3WApB1cv46IiMjwyU7swsLCYGlpKfvGzs7OcHZ2LlFQpH1594jl+DoiIiJDI3vyRMeOHXH06NECzx8/fhyenp5aCYq0L5hbiRERERk82YndgwcPkJiYWOD5xMREREZGaiUo0q6nSSo8S1ZrjqsogEbVizVvhoiIiCqAYi13UtguBXfu3ClWVy2VndwTJ5rZmMDYiBMniIiIDE2hzTY//PADfvzxR83xqlWrsGPHjjzXKZVK3LhxA7169dJ+hFRqHF9HRERUORSa2CUnJyMmJkZznJCQkO9sWAsLC4wbNw5z587VfoRUatxxgoiIqHIoNLEbP348xo8fDwBo0aIFli1bht69e5dJYKQ9TOyIiIgqB9kj6ENCQnQZB+lITIoKDxNVmmNjAXC3YWJHRERkiEq8VyxVDLkXJna3MYGZghMniIiIDFGBLXY2NjYwMjLCkydPYGpqChsbm0JnxQKZs2Zzjskj/QvJtX5dC3bDEhERGawCE7s5c+ZAEAQYGxtLjqliyTO+rgYTOyIiIkNVYGLn6+tb6DFVDHmXOmFiR0REZKg4xs6Axaep8W989sQJAUBzttgREREZLNmJ3Zo1a9CzZ88Cz3t7e2P9+vVaCYq0IzTX+LpG1Y1hYcJcnoiIyFDJ/pTfs2cPWrduXeD51q1b46efftJKUKQdXL+OiIiocpGd2N2/fx9ubm4Fnm/YsCEiIiK0EhRpR0iu8XUeTOyIiIgMmuzEzszMDE+fPi3w/JMnT/Ldboz0J/cadtwjloiIyLDJzsTatGmD77//HnFxcXnOKZVKfP/992jbtq1Wg6OSS84QERaXISlrwYkTREREBk32lmLz5s2Dt7c3OnTogClTpsDd3R0AcOPGDWzatAnR0dHYvn27ruKkYrrxIh0qMfvYqZoC1mZsUSUiIjJkshO71157Dbt27cL06dOxYMECzWLFoijCxcUFu3btQqtWrXQWKBVP3m5YttYREREZOtmJHQB07twZV65cQXBwMO7duwcAqF+/Pjw9PbkrRTkTEiudOMFuWCIiIsNX7L45QRDQsmVLDBo0CIMGDULLli11mtQ9ffoUH3zwARo2bAgHBwe0bdsWZ86c0ZwXRRF+fn5o0qQJatWqhT59+uDmzZs6i6eiyLvUCSdOEBERGbpitdgBQHp6Om7fvo34+Hio1eo85zt06KCVwIDMSRk9e/ZEu3btsHv3btja2iIiIgJ2dnaaa9auXQt/f3/4+/vDzc0NK1aswKBBg3D58mVYWlpqLZaKJF0t4voLaWLXgl2xREREBk92YieKIj7//HN8++23SExMLPC62NhYrQQGAOvWrUOtWrWwefNmTZmLi4skpoCAAMyYMQMDBgwAAAQEBMDNzQ179+7F2LFjtRZLRXJbmYHU7J3EYG9uhFpVFfoLiIiIiMqE7K7Yr7/+GmvWrIGPjw82bdoEURSxePFirFmzBu7u7vDw8MCBAwe0GtyRI0fg5eWFsWPHwtXVFR07dsQ333wDUcyc7hkREYGoqCh069ZN8xhzc3O0b98eFy9e1GosFUlIrq3EPDm+joiIqFKQ3WL3/fffo3///vj66681rXKenp7o3Lkzhg8fju7du+PMmTPo3Lmz1oK7f/8+tmzZgqlTp2LGjBkIDQ3F3LlzAQCTJk1CVFQUAEi6ZrOOnzx5UuB9w8PDtRajPp+jIKfvmgDITubqCQkID39RpjHos/7lAevP+ldmrD/rX5npuv6F7QIGFCOxe/jwIaZNmwYAmh0mUlNTAWTuSjFs2DBs3rwZ8+fPL2mseajVarz22mtYtGgRgMxE8u7duwgMDMSkSZNKfN+iXpTSCg8P1/lzFCYiPBpA9qzYLm4OcHMxL7Pn13f99Y31Z/1Zf9a/smL99V9/2V2x1tbWSElJAQBYWVnB1NQUjx490pw3MzPT6vg6AHBwcEDjxo0lZY0aNcLDhw815wEgOjpack10dDTs7e21GktFoRZFXMvdFcuJE0RERJWC7MTO3d0doaGhmQ8yMsLrr7+OLVu24NGjR4iMjMT27du1nqW2a9cOd+7ckZTduXMHjo6OAABnZ2c4ODggKChIcz4lJQXnz5+vtNub3X+pQnx69pYTVqYCnKtx4gQREVFlIDuxGzp0KMLCwjStdgsXLsSdO3fg4eEBT09P/Pvvv1i4cKFWg5s6dSouX76MVatW4e7du/j555/xzTffYMKECQAy19SbMmUK1q5di0OHDuHGjRuYOnUqLCwsMGTIEK3GUlHk3nGiRQ0TLh5NRERUScgeYzdy5EiMHDlSc/zGG2/gwoULOHr0KIyNjdG9e3c0bNhQq8G9/vrr2LlzJ5YsWYKVK1eiXr16+OSTTzSJHQBMnz4dycnJmD17NpRKJby8vLB///5Ku4Zdnh0n2A1LRERUaRR7geKcXFxcMHXqVG3Fkq+ePXuiZ8+eBZ4XBAG+vr7w9fXVaRwVBXecICIiqryKndjdvXsXv/76Kx48eAAgc5zb22+/jQYNGmg9OCoeURTz7YolIiKiykF2YqdSqTB37lxs27Ytz1Zin3zyCcaMGYMVK1ZAoeBAfX15mqxGdEr296aKAnCrXqpGWSIiIqpAZE+e+OKLL7BlyxYMHToUQUFBePDgAR48eICgoCAMHToU27ZtwxdffKHLWKkIwTHS8XXNa5jA2IgTJ4iIiCoL2c05O3fuxMCBA7Fp0yZJecuWLbFp0yYkJydj586dmsWEqezl7Ybl+DoiIqLKRHaLXVJSEjp27Fjg+TfffBNJSUlaCYpKJu/ECY6vIyIiqkxkJ3bt27fHhQsXCjx/4cIFtG/fXitBUcmE5NpxgkudEBERVS6yE7uvvvoKISEh+M9//oOwsDCkp6cjPT0dYWFh+PjjjxEaGorVq1frMlYqxItUNSITVJpjhQC4WzOxIyIiqkxkj7Fr3bo1RFHE7du3sW3bNs1uBqKYuX2VsbExWrduLXmMIAh4/PixFsOlgoTkmjjRxNoYVYw5cYKIiKgykZ3YDRo0iFtTlWN5Jk5wYWIiIqJKR3ZiFxAQoMs4qJSCc4+v48LERERElY7sMXZUvuVuseOMWCIiosqn2NsSPH78GMHBwYiPj8+zAwUAjBgxQiuBkXyJ6WqEx2VIypqzxY6IiKjSkZ3YpaamYtq0aThw4ADUajUEQdBMnMg59o6JXdm7FpsOMcdxQysFrEzZGEtERFTZyP70X7p0KX7++WfMnz8fhw8fhiiKCAgIwIEDB9CtWzd4eHjg7NmzuoyVCpBn/TruOEFERFQpyU7sDhw4gOHDh+Pjjz+Gu7s7AKB27dro0qUL9uzZg6pVq2Lr1q06C5QKlnvHCS5MTEREVDnJTuyePXumWafO2DizBzclJQVAZlfsgAEDcOjQIR2ESEXhxAkiIiICipHY1axZE/Hx8QAAS0tLmJub4969e5rz6enpSExM1H6EVKg0lYibSrbYERERUTEmT3h4eOCff/4BkNlC16FDBwQEBMDT0xNqtRrffPMNPDw8dBYo5e+mMh3pOSYn162qQM0qCv0FRERERHoju8Vu9OjRyMjI0HS/LlmyBAkJCejTpw/69u2LpKQkLF26VGeBUv5yd8N6sLWOiIio0pLdYuft7Q1vb2/Nsbu7O65cuYK//voLCoUC7dq1g7W1tS5ipELk3UqMiR0REVFlVewFinOysrJCnz59tBULlUDupU48uTAxERFRpVVoV2xcXBx8fHzw1VdfFXqTVatWYciQIUhISNBqcFQ4lVrEtdxr2LHFjoiIqNIqNLELDAzEpUuXMGrUqEJvMmrUKFy8eJHr2JWxf+MzkJiRvedEDTMj1LPgxAkiIqLKqtDE7siRIxg4cCDs7OwKvYm9vT0GDx6MgwcPajU4KlyeHSdsTSTbuxEREVHlUmhid/v2bbz22muybuTp6Ynbt29rJSiSJ8/ECY6vIyIiqtQKTexUKpVml4miGBsbIyMjQytBkTy5txLjjhNERESVW6GJXZ06dXDt2jVZN7p27Rpq166tlaCoaKIoIiQ2TVLGiRNERESVW6GJXbdu3bBr1y5ERUUVepOnT59i165d6N69u1aDo4JFJqrwIjV74kQ1YwENrUq1eg0RERFVcIUmdh999BFUKhX69++PS5cu5XvNpUuXMGDAAKhUKnz44Yc6CZLyyj2+rnkNExhx4gQREVGlVmgTj6OjI3bs2IGxY8eiV69ecHZ2RrNmzVCtWjUkJCTgxo0buH//PiwsLLBt2zY4OTmVVdyVXu4ZsdxKjIiIiIrsu+vevTvOnDmDtWvX4vjx4zhy5IjmXO3atTF27Fh8+OGHcHFx0WWclAsnThAREVFusgZlOTk54auvvsJXX32Fly9f4uXLl7C0tISlpaWu46MChMbkmjjBpU6IiIgqvWKPtmdCp3/RySo8TlJrjk2NgCbWTOyIiIgqu0InT1D5lHt8nbuNCUwVnDhBRERU2TGxq4C44wQRERHlh4ldBcSJE0RERJQfJnYVUEjuiRNM7IiIiAjFSOwiIyORnJxc4Pnk5GRERkZqJaiCrF69GtbW1pg9e7amTBRF+Pn5oUmTJqhVqxb69OmDmzdv6jQOfYpLU+PuS5Xm2EgAmtkwsSMiIqJiJHaenp44fPhwgeePHTsGT09PrQSVn8uXL2P79u1o1qyZpHzt2rXw9/fH8uXLcfLkSdjZ2WHQoEF4+fKlzmLRp2u5Jk64WRnDwoQNr0RERFSMxE4UxULPZ2RkQNDRllZxcXGYOHEiNmzYAGtra0lMAQEBmDFjBgYMGICmTZsiICAACQkJ2Lt3r05i0bc8EyfYDUtERESvFKupp6DELS4uDr///jvs7Oy0ElRuWYnbm2++KSmPiIhAVFQUunXrpikzNzdH+/btcfHiRZ3Eom/BHF9HREREBSh0geJly5ZhxYoVADKTukmTJmHSpEn5XiuKIqZOnar1AHfs2IG7d+/im2++yXMuKioKAPIklHZ2dnjy5EmB9wwPD9dukGX4HH8/qYKc+bht8jOEhz/VyXOVRlm8xuUZ68/6V2asP+tfmem6/m5uboWeLzSx8/Lywvjx4wEAgYGB6Nq1Kxo2bCi5RhAEWFhYoGXLlujfv38pw5UKDw/HkiVLcPz4cZiYaK9lqqgXpbTCw8N18hzJGSLunX0sKevdoj6szcrXGDtd1b+iYP1Zf9af9a+sWH/917/QxO7tt9/G22+/DQBITEzEuHHj0KpVqzIJDAAuXbqEmJgYtGvXTlOmUqlw7tw5bN26FRcuXAAAREdHw9HRUXNNdHQ07O3tyyzOsnLzRTpUOYY6OlVTlLukjoiIiPRH9l6xGzdu1GUc+erTpw9ee+01Sdm0adPQsGFDfPzxx3B1dYWDgwOCgoLw+uuvAwBSUlJw/vx5LFmypMzj1bXcW4lxxwkiIiLKSXZil+Xx48cIDg5GfHw81Gp1nvMjRozQSmAAYG1tLZkFCwBVq1aFjY0NmjZtCgCYMmUKVq9eDTc3N7i6umLVqlWwsLDAkCFDtBZHeZF74gR3nCAiIqKcZCd2qampmDZtGg4cOAC1Wg1BEDRLoOScLavNxE6O6dOnIzk5GbNnz4ZSqYSXlxf2798PS0vLMo2jLORd6sRUT5EQERFReSQ7sVu6dCl+/vlnzJ8/H23btkXfvn0REBCAWrVqYcOGDYiOjsamTZt0GSsA4MiRI5JjQRDg6+sLX19fnT+3PmWoRVx/wT1iiYiIqGCyR94fOHAAw4cPx8cffwx3d3cAQO3atdGlSxfs2bMHVatWxdatW3UWaGV3Oy4DKdk7icHe3Ai1qir0FxARERGVO7ITu2fPnqF169YAAGPjzIa+lJQUAJmtZgMGDMChQ4d0ECIB+XTDcuIEERER5SI7satZsybi4+MBAJaWljA3N8e9e/c059PT05GYmKj9CAkAJ04QERFR0WSPsfPw8MA///wDILOFrkOHDggICICnpyfUajW++eYbeHh46CzQyi7PUiecOEFERES5yG6xGz16NDIyMjTdr0uWLEFCQgL69OmDvn37IikpCUuXLtVZoJWZWhQRGsOJE0RERFQ42S123t7e8Pb21hy7u7vjypUr+Ouvv6BQKNCuXbs8a86RdkS8VCE+PXvLCStTAc7VOHGCiIiIpIq9QHFOVlZW6NOnj7ZioQLk7ob1qGEiWTuQiIiICChGVywApKWl4bvvvsPEiRMxcOBABAcHAwCUSiV+/PFHPHr0SCdBVnacOEFERERyyG6xi42NRb9+/XDjxg3Y29sjOjoaSqUSQGbL3dKlS3Hr1i189tlnuoq10sq71AknThAREVFeslvsFi1ahMjISBw/fhznzp3TbCcGAEZGRujfvz9+++03nQRZmYmiiGBOnCAiIiIZZCd2x48fx+TJk9G2bdt8x3c1bNgQDx8+1GpwBDxNViM6Ra05rqIA3KqXamgkERERGSjZid3Lly9Rr169As+npqZCpVIVeJ5KJnc3bDMbExgbceIEERER5SU7sWvQoAGuXLlS4PmTJ09q9pAl7ck7cYLj64iIiCh/xVqg+IcffsDu3buhVmd2DQqCgKSkJCxevBgnT57E2LFjdRZoZZVn4gTH1xEREVEBZA/Wmjx5Mm7duoXJkyfD0tISADBu3DgolUqoVCpMmDABI0eO1FmglVXuNew4cYKIiIgKUqxR+GvWrMHw4cNx4MAB3L17F2q1GvXr18egQYPQvn17XcVYab1IVeNBQva4RYUAuFszsSMiIqL8FXt6Zdu2bdG2bVtdxEK55O6GbWxtjCrGnDhBRERE+SvWzhNUtkI4cYKIiIiKodAWu379+hXrZoIg4NChQ6UKiLLlHl/Xoga7YYmIiKhghSZ2z58/lyxGrFKpEB4eDkdHR1hYWOg8uMoud1csJ04QERFRYQpN7M6fPy85jomJgaurK9atW4fOnTvrNLDKLjFdjdtxGZKy5myxIyIiokIUa4xdfluJkW5cf5EOMcdxA0sFrEw5JJKIiIgKxkyhnArO0w3LiRNERERUOCZ25RR3nCAiIqLiYmJXTuWZEcvEjoiIiIpQ6OSJ6OhoyXFsbCwAIC4uLs+5LHZ2dloKrfJKU4m48YJLnRAREVHxFJrYNWrUKN8JE2PGjCnwMVnJH5XcLWU60tXZx3WqGsHOXKG/gIiIiKhCKDSxmzNnDmfC6kHeblhOnCAiIqKiFZrY+fr6llUclEPuGbEcX0dERERycPJEORSaO7Hj+DoiIiKSgYldOaNSiwiN5VZiREREVHxM7MqZuy8zkJiRveeEjZmAehacOEFERERFY2JXzuRemNjT1pQTWIiIiEgWJnblTJ6JExxfR0RERDIxsStnuOMEERERlZTsxO7ChQtYvXp1gefXrFmDS5cuaSWoykoURQTHpEnKOHGCiIiI5Cp0Hbucli9fDmtr6wLPX7t2DWfOnMG+ffu0EVel9DBRhRep2RMnLIwFNLSS/S0iIiKiSk52i11ISAjatGlT4PnWrVsjODhYK0FlWb16Nbp27QpHR0c0bNgQw4YNw40bNyTXiKIIPz8/NGnSBLVq1UKfPn1w8+ZNrcZRVnJPnPCoYQIjTpwgIiIimWQndklJSUXOzkxISCh1QDmdOXMG48ePx4kTJ3Do0CEYGxtj4MCBePHiheaatWvXwt/fH8uXL8fJkydhZ2eHQYMG4eXLl1qNpSwE5xpf58FuWCIiIioG2Ymdq6srTp48WeD533//HQ0aNNBKUFn279+P9957D02bNkWzZs2wefNmPH/+HBcuXACQ2VoXEBCAGTNmYMCAAWjatCkCAgKQkJCAvXv3ajWWspC7xY4zYomIiKg4ZA/gGjVqFObMmYM5c+bA19cXNjY2AIDY2Fj4+fnh5MmTWLp0qc4CBTJbBNVqtWasX0REBKKiotCtWzfNNebm5mjfvj0uXryIsWPH5nuf8PBwncZZ0uf4X1QV5My1bRKfIDxcLPgB5VhZvMblGevP+ldmrD/rX5npuv5ubm6Fnped2E2cOBGhoaH49ttvERgYCHt7ewDAs2fPIIoi3n33XUyZMqV00RZh3rx58PDw0Iz1i4qKAgDY2dlJrrOzs8OTJ08KvE9RL0pphYeHF/s5nqeo8OzMU82xiRHQo0VDmCoq3hi7ktTfkLD+rD/rz/pXVqy//utfrCmX69atw9ChQ3Ho0CHcv38fAODi4oIBAwagY8eOuohP45NPPsGFCxdw/PhxKBSGt8VW7m7YpjYmFTKpIyIiIv0p9loanTp1QqdOnXQRS4F8fX2xf/9+/PLLL3BxcdGUOzg4AACio6Ph6OioKY+Ojta0KFYU3HGCiIiISqvc7zwxd+5c7Nu3D4cOHUKjRo0k55ydneHg4ICgoCBNWUpKCs6fP4+2bduWdailkmfiBGfEEhERUTEV2GLXt29fGBkZYf/+/TA2Nka/fv2KvJkgCDh06JDWgps1axZ27dqF77//HtbW1poxdRYWFqhWrRoEQcCUKVOwevVquLm5wdXVFatWrYKFhQWGDBmitTjKAnecICIiotIqMLETRRFqtVpzrFari1zHThS1O4MzMDAQADBgwABJ+dy5c+Hr6wsAmD59OpKTkzF79mwolUp4eXlh//79sLS01GosuhSfpsbdlyrNsQCgmQ0TOyIiIiqeAhO7I0eOFHpcFpRKZZHXCIIAX19fTaJXEV3LtTBxo+rGsDAp973kREREVM7Iyh6Sk5M1OzuQ9uWZOMFuWCIiIioBWYmdubk51qxZg4cPH+o6nkopJJYzYomIiKj0ZPf3NW/eHHfv3tVlLJVWSK6JEy1sTfUUCREREVVkshO7Tz/9FDt27MCJEyd0GU+lk5Ih4pYyQ1LGrlgiIiIqCdkLFG/YsAE2NjYYMWIE6tSpAxcXF5ibm0uuEQQBu3fv1nqQhuymMh2qHJOJnaopYGPGiRNERERUfLITu1u3bkEQBNSrVw8A8ODBgzzXFLUcCuXFHSeIiIhIW2QndqGhobqMo9LijhNERESkLezz07OQ2Nw7TnDiBBEREZWM7Ba7nF6+fIn4+HjJzhRZHB0dSx1UZZGhFvMsTswWOyIiIiqpYiV227dvx/r163Hv3r0Cr4mNjS11UJVFeFwGUrJ3EoO9uRFqmbMRlYiIiEpGdhbx3XffYebMmXB2dsaCBQsgiiKmTJmCmTNnwt7eHh4eHli/fr0uYzU4+U2c4AQUIiIiKinZid2mTZvQpUsX7N+/H2PGjAEA9OjRA59++ikuXLgApVKJ+Ph4XcVpkHKPr2M3LBEREZWG7MTu7t276N27d+aDjDIflp6e2eJkbW2NUaNGITAwUAchGq7cM2I5cYKIiIhKQ3ZiZ2FhAVHMXEm3WrVqUCgUePLkieZ8jRo18PjxY+1HaKBEUeQesURERKRVshO7Ro0aISwsDABgbGwMDw8P7Nq1C+np6UhJScGuXbvg7Oyss0ANTUSCCvFp2VtOWJkKcLFU6DEiIiIiquhkJ3a9e/fG8ePHkZKSAgCYNWsWzp07BxcXF7i6uuLixYuYOXOmzgI1NLknTnhw4gQRERGVkuzlTj788EN8+OGHmuM+ffrgyJEjOHjwIIyNjdGrVy907NhRJ0EaopCYXBMn2A1LREREpVSiBYqzvPHGG3jjjTe0FUulwokTREREpG3FTuxevnyJv/76Cw8ePAAAODs7o2PHjrC0tNR6cIZKFEUEc8cJIiIi0rJiJXbr1q3DihUrkJSUpJkhCwBVq1bFnDlzMH36dK0HaIieJKnxLDl7OzZzhYBG1UvVeEpEREQkP7Fbv349Fi1ahI4dO2LChAlwdXUFANy5cwfffvstPvvsMxgZGUnG4VH+ruYaX9e8hjGMjThxgoiIiEpHdmK3efNmdO3aFfv27ZPM3mzWrBn69++PQYMGYfPmzUzsZLiaa3xdS46vIyIiIi2QvdxJbGwsevfune+SHIIgoG/fvoiNjdVqcIYq+Lm0xc6zJsfXERERUenJTuw8PT1x69atAs/fvHkTnp6eWgnK0LHFjoiIiHRBdlfsypUr4ePjA0dHR4wfPx7VqlUDACQkJCAwMBBHjhzBvn37dBaooYhMyEBUjokTZgqgsTUnThAREVHpFZhRtG3bNk+ZIAj47LPP8Pnnn8Pe3h4A8OzZM6jVatjb22P8+PG4cOGC7qI1AOejpN2wr9mawoQTJ4iIiEgLCkzsatasmWc8nZ2dnWY2bJYGDRroJjIDde5pquS4Qy12wxIREZF2FJjYHTlypCzjqDTO5Wqxa1/LTE+REBERkaGRPXmCSu9Jkgq34zI0x0YC0MaeLXZERESkHcUatZ+eno4dO3bg119/1Wwp5uTkhF69euH999+HiQmX7SjMbw9TJMev1zSBpQlzayIiItIO2YmdUqlE//79ERoaCnt7e83YuuDgYPz222/YsWMHDh48CGtra13FWuEdeyBN7HrWq6KnSIiIiMgQyW4u+uyzz3Dz5k34+/vj5s2bOHbsGI4dO4Zbt24hICAAN2/exJIlS3QZa4WWmK7Gn0+kEyd6OjKxIyIiIu2RndgdPXoUEydOxLvvvgsjo+yHCYKA4cOHY8KECZxwUYhDESlIyhA1x3WrKuBRg13XREREpD2yE7u4uDjUr1+/wPP169dHXFycVoIyRDvDEyXHQxua57s9GxEREVFJyU7sGjRogKNHj0IUxTznRFHEkSNHuKZdAa4+T8OZp9JlTka4VtVTNERERGSoZCd2EyZMwKlTp+Dj44Nff/0Vd+/exd27d3HixAn4+Pjg9OnTmDRpki5jrbD8rr6UHL/hYIrG1uyGJSIiIu2SPSt23LhxiImJwapVq3Dq1ClNuSiKMDU1xSeffIIxY8boIMSK7eD9ZJyIlM6G/biFpZ6iISIiIkNWrHXsZs+ejXHjxuHUqVOIjIwEADg6OqJr166oUaOGTgKUKzAwEOvWrUNUVBSaNGkCPz8/tG/fXq8x/fk4BdP+eiEpa21ngrfqcrcJIiIi0r5iJXYAYGtrCx8fH13EUmL79+/HvHnz8NVXX6Fdu3YIDAzE0KFDceHCBTg6OpZpLNdj07Holile3HqGf56nS84ZC8Ca9jacNEFEREQ6IXuMXVhYWJ7lTM6ePYvBgweje/fu2Lhxo9aDk8vf3x/vvvsuRo8ejcaNG2PlypVwcHDA1q1byzyW+HQ1fn9unCepA4Av21RHcy5xQkRERDoiKJXKvNNc8zF06FAIgoDdu3cDAB49eoS2bdvCzMwMdnZ2uH37NjZs2IB3331XpwHnlpaWhtq1a2PLli0YOHCgpnzWrFm4ceMGjh49mucx4eHhOovnTqKAEVfM85RPc07DGMeMfB5BREREJI+bm1uh52V3xQYHB2PatGma4127dkGtVuPMmTOoXbs2RowYgcDAwDJP7GJiYqBSqWBnZycpt7Ozw7Nnz/J9TFEvSmlYJKqAK081x23tTfGplxU61qo84+rCw8N1+hqXd6w/68/6s/6VFeuv//rLTuzi4uJga2urOf7tt9/QqVMn1K5dGwDQs2dPLFy4UPsRVjA1qxhhSaNUuDvXQfMaJqhdVaHvkIiIiKiSkD3Gzs7ODg8ePAAAKJVK/P333+jatavmfGpqakEP1SlbW1soFApER0dLyqOjo2Fvb1/m8ZgqBHjbq/B2vSpM6oiIiKhMyW6x69q1K7755htYWVnhzJkzAIDevXtrzt+6dQt169bVfoRFMDU1RcuWLREUFCQZYxcUFIT+/fuXeTxERERE+iI7sVu4cCHu3LmDTz/9FKampliyZAmcnJwAACkpKfj555/xzjvv6CzQwkybNg2TJ0+Gl5cX2rZti61bt+Lp06cYO3asXuIhIiIi0gfZiZ2dnR2OHTuGuLg4mJubw9TUVHNOFEUcOnQI9erV00mQRRk8eDBiY2OxcuVKREVFwd3dHbt379YknkRERESVQbEXKK5evXqeMnNzc3h4eGgloJKaMGECJkyYoNcYiIiIiPSp0MTun3/+KfYNvby8ShwMEREREZVcoYndW2+9JXv7K1EUIQgCYmNjtRIYERERERVPoYmdv79/WcVBRERERKVUaGJX1rtIEBEREVHJyV6gmIiIiIjKNyZ2RERERAaCiR0RERGRgRCUSqWo7yCIiIiIqPTYYkdERERkIJjYERERERkIJnZEREREBoKJHREREZGBYGJHREREZCCY2BEREREZCCZ2WhQYGIgWLVrAwcEBnTt3xrlz5/QdUqmtXr0aXbt2haOjIxo2bIhhw4bhxo0bkmumTJkCa2tryb+33npLck1qaipmz56NBg0aoE6dOhg+fDgePXpUllUpMT8/vzz1a9Sokea8KIrw8/NDkyZNUKtWLfTp0wc3b96U3EOpVGLSpElwcnKCk5MTJk2aBKVSWcY1KRkPD4889be2tsY777wDoOjXB5D3GpUXZ8+exfDhw+Hu7g5ra2vs3LlTcl5b3+/r16+jd+/eqFWrFtzd3bF8+XKIov5Xnyqs/unp6Vi0aBHat2+POnXqoHHjxpgwYQIiIyMl9+jTp0+en4lx48ZJrimvvxNFff+19X4XGRmJYcOGoU6dOmjQoAHmzJmDtLQ0ndevKEXVP7/3Amtra8yaNUtzTUX+TJDzmVfe3wOY2GnJ/v37MW/ePPznP//B6dOn0aZNGwwdOjTPG15Fc+bMGYwfPx4nTpzAoUOHYGxsjIEDB+LFixeS67p06YKwsDDNvz179kjO+/r64pdffsGWLVtw9OhRvHz5EsOGDYNKpSrL6pSYm5ubpH45k/a1a9fC398fy5cvx8mTJ2FnZ4dBgwbh5cuXmmsmTJiAkJAQ7N27F3v37kVISAgmT56sj6oUW1BQkKTuf/75JwRBwMCBAzXXFPb6APJeo/IiMTERTZs2xbJly2Bubp7nvDa+3/Hx8Rg0aBDs7e1x8uRJLFu2DOvXr8eGDRvKpI6FKaz+SUlJCA4OxqxZs/Dnn3/ihx9+wKNHjzBkyBBkZGRIrh05cqTkZ2LNmjWS8+X1d6Ko7z9Q+vc7lUqFYcOGISEhAUePHsWWLVtw6NAhzJ8/X+f1K0pR9c9Z77CwMPz0008AIHk/ACruZ4Kcz7zy/h7ABYq1pHv37mjWrBnWrVunKXv99dcxYMAALFq0SI+RaVdCQgKcnJywc+dOeHt7A8j86yw2Nha7du3K9zFxcXFwdXWFv7+/ppXn4cOH8PDwwN69e9G9e/cyi78k/Pz8cOjQIZw/fz7POVEU0aRJE0ycOFHzF2tycjLc3Nzw+eefY+zYsQgLC0Pbtm1x/PhxtGvXDgBw/vx5eHt74/Lly3BzcyvT+pTWqlWrsG7dOoSFhcHc3LzQ1weQ9xqVV3Xr1sWKFSswcuRIANr7fm/ZsgWLFy/G7du3NR+eK1euxNatW3Hjxg0IgqCfCueSu/75uXXrFtq1a4ezZ8+iWbNmADJb7Jo2bYqVK1fm+5iK8juRX/218X7322+/4Z133kFoaCjq1asHANi1axc++ugjhIeHw8rKSveVk0HO9/+jjz7CuXPn8Pfff2vKDOkzIfdnXkV4D2CLnRakpaXh6tWr6Natm6S8W7duuHjxop6i0o2EhASo1WpYW1tLys+fPw9XV1d4eXnho48+QnR0tObc1atXkZ6eLnl96tWrh8aNG1eY1+f+/fto0qQJWrRogXHjxuH+/fsAgIiICERFRUnqZm5ujvbt22vqdunSJVSrVg1t27bVXNOuXTtYWFhUmPpnEUUR//3vfzFs2DDJX/MFvT6AvNeootDW9/vSpUt44403JK9h9+7d8eTJE0RERJRRbbQjq5Ui93vCvn370KBBA7Rr1w4LFiyQtGZU9N+J0r7fXbp0CY0bN9YkdUDm9z81NRVXr14ts3qUVkJCAvbv34/Ro0fnOWconwm5P/MqwnuAcakeTQCAmJgYqFQq2NnZScrt7Ozw7NkzPUWlG/PmzYOHhwfatGmjKXvrrbfQr18/ODs748GDB/jiiy/Qv39/nDp1CmZmZnj27BkUCgVsbW0l96oor0+rVq2wceNGuLm54fnz51i5ciV69OiBCxcuICoqCgDy/d4/efIEAPDs2TPY2tpK/gITBAE1a9asEPXPKSgoCBERERg1apSmrLDXp0aNGrJeo4pCW9/vZ8+eoU6dOnnukXXOxcVFV1XQqrS0NCxYsAC9evVC3bp1NeVDhw6Fo6MjatWqhVu3buGzzz7D9evXceDAAQAV+3dCG+93z549y/MzZGtrC4VCUe7rn9PevXuRlpaGESNGSMoN6TMh92deRXgPYGJHsn3yySe4cOECjh8/DoVCoSn38fHRfN2sWTO0bNkSHh4eOHHiBPr376+PULXq7bfflhy3atUKLVu2xA8//IDWrVvrKSr92LFjB15//XV4eHhoygp7ff7v//6vrEOkMpKRkYFJkyYhLi4OP/74o+TcmDFjNF83a9YMLi4u6N69O65evYqWLVuWbaBaZujvd8WxY8cO9O7dGzVr1pSUG8prVNBnXnnHrlgtyPpLK2dTMwBER0fD3t5eT1Fpl6+vL/bt24dDhw4V+ZdE7dq1UadOHdy9excAYG9vD5VKhZiYGMl1FfX1qVatGpo0aYK7d+/CwcEBAAr93tvb2yMmJkYy20kURTx//rxC1T86OhpHjx7Nt9slp5yvDwBZr1FFoa3vt729fb73yDpX3mVkZGD8+PG4fv06Dh48iBo1ahR6/WuvvQaFQiF5TzCE3wmgZO93+X3/s3p+Kkr9Q0JCcOXKlSLfD4CK+ZlQ0GdeRXgPYGKnBaampmjZsiWCgoIk5UFBQZI+9opq7ty5mh/w3MtY5CcmJgZPnjzR/AK0bNkSJiYmktfn0aNHmgGmFU1KSgrCw8Ph4OAAZ2dnODg4SOqWkpKC8+fPa+rWpk0bJCQk4NKlS5prLl26hMTExApV/x9++AFmZmaSv8bzk/P1ASDrNaootPX9btOmDc6fP4+UlBTNNUFBQahduzacnZ3LqDYlk56ejrFjx+L69ev45ZdfNN/nwly/fh0qlUpzraH8TgAle79r06YNwsLCJMt7BAUFwczMrMK0aO7YsQPOzs7o0qVLkddWtM+Ewj7zKsJ7ALtitWTatGmYPHkyvLy80LZtW2zduhVPnz4t1zP+5Jg1axZ27dqF77//HtbW1prxBRYWFqhWrRoSEhKwbNky9O/fHw4ODnjw4AGWLFkCOzs79O3bFwBQvXp1vP/++1i0aBHs7OxgY2OD+fPno1mzZrLeFPQtawxRvXr1NGPIkpKSMGLECAiCgClTpmD16tVwc3ODq6srVq1aBQsLCwwZMgQA0LhxY7z11luYOXMmvv76awDAzJkz0bNnz3Iz+68ooijiu+++w+DBg1GtWjXJucJeHwCyXqPyJCEhQdOyoFar8fDhQ4SEhMDGxgaOjo5a+X4PGTIEy5cvx9SpUzFr1izcuXMHX3/9NebMmaP3GbGF1b927doYPXo0rly5gh9//BGCIGjeE6ysrGBubo579+5h9+7d6NGjB2rUqIGwsDAsWLAALVq00MwQLM+/E4XV38bGRivvd926dYO7uzs++OADfPHFF3jx4gUWLlyIUaNG6X1GbFE//0Dmsjd79uzBRx99lOfntaJ/JhT1maet93xdvgdwuRMtCgwMxNq1axEVFQV3d3d8+eWX6NChg77DKpXcM92yzJ07F76+vkhOTsbIkSMREhKCuLg4ODg4oFOnTpg/f75kxldqaioWLFiAvXv3IiUlBW+++Sa++uoryTXl1bhx43Du3DnExMSgZs2aaNWqFebPn48mTZoAyEx6li1bhu3bt0OpVMLLywurVq1C06ZNNfdQKpWYM2cOjh07BgDw9vbGihUrCnx9y5vTp0+jf//++OOPP+Dl5SU5V9TrA8h7jcqLv/76C/369ctTPmLECAQEBGjt+339+nXMmjUL//vf/2BtbY2xY8di7ty5ek/sCqv/vHnz4Onpme/j/P39MXLkSDx8+BCTJk3CzZs3kZiYiLp166JHjx6YN28ebGxsNNeX19+Jwuq/evVqrb3fRUZGYtasWTh9+jSqVKmCoUOH4vPPP4eZmVmZ1LMgRf38A8D333+P6dOn49q1a6hdu7bkuor+mVDUZx6gvfd8Xb0HMLEjIiIiMhAcY0dERERkIJjYERERERkIJnZEREREBoKJHREREZGBYGJHREREZCCY2BEREREZCCZ2REQEURTRoUMHfPHFFyV6/LfffovmzZsjNTVVy5ERUXEwsSOiUrG2tpb1b+fOnfoOVW9u3boFPz8/RERE6DuUAu3duxf37t3DlClTNGV//fUXrK2tsW/fPsm1GRkZGD16NKytrbFx40YAwHvvvYfU1FRs27atTOMmIiluKUZEpbJ582bJ8fbt2/H3339jw4YNkvLysAekvoSFhWH58uXo2LFjud0Ldt26dejfvz9sbW0LvS4jIwMTJkzAwYMH8eWXX2Lq1KkAAHNzc4wYMQIbNmzApEmTYGTEdgMifWBiR0SlMmzYMMnxqVOn8L///S9PuSFJTEyEhYWFvsNAWloajIyMYGxcurfy4OBghIaGYsGCBYVep1KpMHHiRPz888+SpC7L4MGDsW7dOpw+fVrve34SVVb8k4qIysSePXvQtWtX1KpVC87Ozhg9ejTu378vuaZPnz5o3bo1rl+/jt69e6N27drw9PTE/v37AQDnz5/HW2+9hVq1aqFVq1b4448/JI/38/ODtbU1bt26hQkTJsDJyQnOzs6YMWMGEhIS8sR08uRJ9O7dG3Xr1kXdunXh4+ODkJAQyTVTpkyBg4MDIiIiMHz4cDg6OuKdd94BAFy7dg1Tp05Fy5Yt4eDggAYNGmDcuHGIjIzUPH7nzp0YPXo0AKBfv355uqY9PDwk3Z85X4s+ffpojrO6RXfv3g0/Pz80b94ctWrVwqNHjwAAd+7cwZgxY1C/fn3N/pwHDx6U9b05cuQIFAoF3nzzzQKvUalUmDRpEg4cOJBvUgcALVu2hI2NDQ4fPizreYlI+9hiR0Q6t2bNGixZsgQDBgzAyJEjoVQq8e2336JXr144c+YMatasqbk2Pj4ew4YNw8CBAzFw4EBs27YNEydOhCiK8PX1xdixY+Hj44MNGzZgzJgxuHbtGqpXry55vnHjxqFOnTr49NNPERoaiu3bt+PRo0fYs2eP5po9e/Zg0qRJ6Nq1KxYuXIi0tDRs374dvXv3xsmTJ9GoUSPNtWq1GoMHD4aXlxeWLFkChUIBAAgKCsKdO3cwfPhw1K5dG/fu3cPWrVvxzz//4Pz586hatSo6dOiAyZMnY/PmzfjPf/6juW9Ju6ZXr14NIyMjTJ48GQBQrVo1hIWFoUePHnBwcMD06dNhYWGBw4cPY/To0di8eXORraeXLl1C48aNYW5unu95lUqFyZMnY9++fQUmdVk8PT1x4cKFEtWNiEqPiR0R6VRkZCSWLl2KefPmYe7cuZpyHx8ftGvXDhs3bsTChQs15VFRUdi0aROGDx8OAOjatStat26NCRMm4NixY2jXrh0AoHHjxhg8eDAOHjyIUaNGSZ6zTp062LNnDwRBAAA4ODhg5cqVOHXqFLp06YLExETMnj0b7777Lvz9/TWPe//999GqVSusWLECgYGBmvL09HT07NkTX375peR5xo8fjw8//FBS5u3tjZ49e+KXX37BsGHD4OLigvbt22Pz5s3o0qULOnXqVJqXEwkJCbh48aKkK3jChAmoXbs2goKCNMnZxIkTMWjQIHz22Wd45513NK9FfsLDw+Hp6Vng+c8//xyRkZFFJnUA4OLignPnzhWzVkSkLeyKJSKd+uWXX5CRkYHBgwcjJiZG88/KygpNmzbFX3/9Jbne3Nxc09UJAG5ubqhevTpcXV01SR0AeHl5AUCe7lwgM6nJmch88MEHAIDjx48DyGxpUyqVGDp0qCQmlUqFN954I09MQGbylFvVqlU1XyckJCA2Nhaurq6oXr06rl69KuPVKb7hw4dLkroXL17g1KlTGDhwIJKSkiT16d69Ox4/fow7d+4Ues/Y2FhYW1sXeD46OhoKhQINGjQoMj5ra2ukpaUhPj5edp2ISHvYYkdEOvXvv/8CAFq3bp3veRcXF8lx7dq188yotLKyQt26dSVlWd2vSqUyzz0bNmwoOba1tYW1tTUePHggiWngwIH5xpT7+Y2MjODk5JTnOqVSicWLF+PgwYN48eKF5JyuEpvcr9fdu3chiiKWLVuGZcuW5fuY6OhouLm5FXpfURQLPLdo0SIEBgZizJgx2LNnT6Gtjln3KayFkIh0h4kdEemUWq0GkLlOWn6zN6tUqSI5zhq/lltB5YUlJEXFtHHjRtSpU6fI601MTPKNfcyYMbh48SKmTZuGFi1awNLSEoIgYNy4cZrnKEpBCZBarc53yZDc4+Cynmfq1Kno0aNHvvdq2rRpoTHY2trmmyBnsbe3x4EDB+Dt7Y13330XBw8exOuvv57vtUqlEqamprC0tCz0OYlIN5jYEZFO1a9fHwBQr149NGnSpEye899//4Wrq6vmOCYmBkqlUtPqlhVTzZo1S7wsh1KpxKlTpzBv3jzMmzdPU56SklJokpSbtbU14uLi8pRHRkbKWvMuqwXP2Ni4xHVp1KhRkYsnOzk54eeff4a3tzeGDBmCI0eOwN3dPc91ERERkoknRFS2OMaOiHSqf//+UCgUWLFiRb6tazExMVp/zm+//VbyXJs2bQIA9OzZEwDQrVs3VK9eHatXr0ZaWlqexz9//rzI58hqTctdp40bN+ZprcsaE5dfwle/fn38/fffkjiOHz+Ohw8fFhkDANjZ2aFTp07YsWMHHj9+nOe8nLq0bdsWYWFhSE5OLvQ6Nzc37N27VzNmMr/xjcHBwZV6MWoifWOLHRHplIuLCxYvXoxPP/0UkZGR6NOnD6pXr46IiAgcPXoUgwYNgq+vr1af8/Hjxxg6dCh69uyJa9euYceOHejWrRu6du0KIHPM3po1azBx4kS8+eab8PHxgb29PSIjI/HHH3+gSZMmCAgIKPQ5rKys0LFjR6xbtw7p6elwdHTE+fPnce7cOdSoUUNybYsWLaBQKLBmzRrExcXB3NwcXl5ecHFxwahRo3Dw4EH4+Phg0KBBuHfvHnbv3q1pVZRj9erV6NmzJzp06IDRo0ejfv36iI6Oxt9//42wsDBcuXKl0Mf37t0bfn5+OH36tCb5LUjLli3x008/wcfHBwMHDsSxY8dQu3ZtAMDVq1fx4sULyfp7RFS22GJHRDr34YcfYufOnTAxMcGqVaswf/58HD58GB07dixwAkNpbNmyBdbW1vj8889x4MABjBo1Ctu3b5dcM3jwYPzyyy+oV68eNmzYgLlz52Lv3r1o1KgRxo4dK+t5AgMD0bNnT2zbtg0LFy5EXFwcDh06lGdXCnt7e6xduxZKpRLTp0/H+PHjcfbsWQBA9+7d8cUXX+Dff/+Fr68vLl++jF27dska+5fFzc0NQUFB8Pb2xk8//YRZs2Zhy5YtUKlUmD9/fpGP9/DwgKenJw4cOCDr+dq3b4/vvvsOjx49wuDBgxEbGwsAOHDgAOrVq8ddJ4j0SFAqlcUfeUxEVA75+flh+fLlCAsLg4ODg77DqVD27t2L6dOnIzQ0NE+LoxwpKSlo0aIFZsyYUeRad0SkO2yxIyIi+Pj4wMXFBRs3bizR4//73//C1NQU48eP13JkRFQcHGNHREQQBEHTPVwSEydOxMSJE7UYERGVBFvsiIiIiAwEx9gRERERGQi22BEREREZCCZ2RERERAaCiR0RERGRgWBiR0RERGQgmNgRERERGYj/ByxCPIpc+g+3AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 648x432 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Note: 4 moles of atoms per formula unit (Fe3C1). That's why we multiply times 4\n",
    "fig = plt.figure(figsize=(9,6))\n",
    "fig.gca().set_xlabel('Temperature (K)')\n",
    "fig.gca().set_ylabel('Isobaric Heat Capacity (J/mol-formula-K)')\n",
    "fig.gca().plot(result['T'], np.squeeze(4.0 * result['heat_capacity']))\n",
    "plt.show()"
   ]
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.0"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}