{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<!--NOTEBOOK_HEADER-->\n",
    "*This notebook contains course material from [CBE20255](https://jckantor.github.io/CBE20255)\n",
    "by Jeffrey Kantor (jeff at nd.edu); the content is available [on Github](https://github.com/jckantor/CBE20255.git).\n",
    "The text is released under the [CC-BY-NC-ND-4.0 license](https://creativecommons.org/licenses/by-nc-nd/4.0/legalcode),\n",
    "and code is released under the [MIT license](https://opensource.org/licenses/MIT).*"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<!--NAVIGATION-->\n",
    "< [Gases with One Condensable Component](http://nbviewer.jupyter.org/github/jckantor/CBE20255/blob/master/notebooks/07.01-Gases-with-One-Condensable-Component.ipynb) | [Contents](toc.ipynb) | [Operating Limits for a Methanol Lighter](http://nbviewer.jupyter.org/github/jckantor/CBE20255/blob/master/notebooks/07.03-Operating-Limits-for-a-Methanol-Lighter.ipynb) ><p><a href=\"https://colab.research.google.com/github/jckantor/CBE20255/blob/master/notebooks/07.02-Vapor-Liquid-Equilibrium-for-Pure-Components.ipynb\"><img align=\"left\" src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open in Colab\" title=\"Open in Google Colaboratory\"></a>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "yP9LBLF6p4qQ"
   },
   "source": [
    "# Vapor-Liquid Equilibrium for Pure Components"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "yP9LBLF6p4qQ"
   },
   "source": [
    "## Summary\n",
    "\n",
    "This [Jupyter notebook](http://jupyter.org/notebook.html) describes the modeling of vapor-liquid equilibrium with Antoine's equation, including the calculation of saturation pressure, saturation temperature, and normal boiling points."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "k8bnw65Zp4qT"
   },
   "source": [
    "## Gibb's Phase Rule"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "fLhbzTN-p4qT"
   },
   "source": [
    "The Gibb's phase rule shows how many independent intensive thermodynamic variables (e.g. $T$, $P$, $\\hat{V}$, or $x_i$) are required to completely specify the state of a substance.\n",
    "\n",
    "$$ F = C + 2 - \\Pi - r$$\n",
    "\n",
    "where\n",
    "\n",
    "$$\n",
    "\\begin{align*}\n",
    "F & = \\mbox{Thermodynamic Degrees of Freedom} \\\\\n",
    "C & = \\mbox{Number of Components} \\\\\n",
    "\\Pi & = \\mbox{Number of Phases} \\\\\n",
    "r & = \\mbox{number of independent reactions at equilibrium}\n",
    "\\end{align*}\n",
    "$$\n",
    "\n",
    "This simple rule that has profound implications for engineering analysis."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "pbzZr-W4p4qV"
   },
   "source": [
    "## Phase Diagram for a Pure Component"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "OIC0QEjVp4qV"
   },
   "source": [
    "For a pure component (i.e., $C = 1$) and no reactions (i.e., $r=0$), the Gibb's phase rule reads\n",
    "\n",
    "$$ F = 3-\\Pi$$\n",
    "\n",
    "which shows:\n",
    "\n",
    "* Two independent thermodynamic variables, such as $T$ and $P$, are sufficient to specify the state of a single phase.\n",
    "* If two phases are in coexistence, then there must be a relationship between $T$ and $P$.\n",
    "* The coexistence of three phases completely specifies the thermodynamic state.\n",
    "\n",
    "These observations are demonstrated in a 2-dimensional phase diagram for a pure substance.\n",
    "\n",
    "![](https://upload.wikimedia.org/wikipedia/commons/thumb/3/34/Phase-diag2.svg/500px-Phase-diag2.svg.png)  \n",
    "[By Matthieumarechal, CC BY-SA 3.0](https://commons.wikimedia.org/w/index.php?curid=4623701)\n",
    "\n",
    "The green line shows the solid/liquid coexistence (the dashed green line showing the anomolous special case of water)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "wElXeCCHp4qV"
   },
   "source": [
    "### Triple Point"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/",
     "height": 336
    },
    "colab_type": "code",
    "executionInfo": {
     "elapsed": 741,
     "status": "ok",
     "timestamp": 1539185455925,
     "user": {
      "displayName": "Jeffrey Kantor",
      "photoUrl": "",
      "userId": "09038942003589296665"
     },
     "user_tz": 240
    },
    "id": "SSK498Rdp4qX",
    "outputId": "41483d2e-2f09-491d-f872-e47218cc502b"
   },
   "outputs": [
    {
     "data": {
      "image/jpeg": "/9j/4AAQSkZJRgABAQAAAQABAAD/2wCEABALDA4MChAODQ4SERATGCgaGBYWGDEjJR0oOjM9PDkz\nODdASFxOQERXRTc4UG1RV19iZ2hnPk1xeXBkeFxlZ2MBERISGBUYLxoaL2NCOEJjY2NjY2NjY2Nj\nY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY//AABEIAWgB4AMBIgACEQED\nEQH/xAAbAAACAwEBAQAAAAAAAAAAAAAAAgEDBAUGB//EADsQAAICAQIDBQYEBQQDAAMAAAABAgMR\nBCEFEjETIkFRcQYUMmGR0UJSgZIVI0NisRYzcqGCweEkNFP/xAAXAQEBAQEAAAAAAAAAAAAAAAAA\nAQID/8QAHREBAQACAwADAAAAAAAAAAAAAAERMQIDMiEiQf/aAAwDAQACEQMRAD8A+fgAAAAAAAAA\nAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA\nAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA\nAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAaVorX4x+pPuFvnD6gZQNkeG3SeFKterf2LHwf\nUL8dX1f2A54HRXBtS/x1fV/YZ8D1SXxVfV/YDmAdL+C6nOHKr6v7DvgGrUOfnpx/yf2A5QHT/gWr\nzjmq+r+wPgeqX4qvq/sBzAOt/p7WYypU/uf2J/07rMZ56V/5P7AcgDrv2c1iWe0o/c/sQvZ/VvpO\nn9z+wHJA669ndW/6lH7n9hl7Na1td+jf+5/YDjAdWz2f1dc3BzpbXlJ/Ylez2sazz0/uf2JkckDt\nQ9mNdPpOj9z+wL2X1z/qUfuf2KOKB2X7Na1f1KP3P7Ef6b1ucc9H7n9gOOB2H7Na3DfPQ/8Ayf2D\n/Tetf4qf3P7AccDsP2b1qx36N/7n9hl7Ma5rPaUfuf2A4oHY/wBN6389P7n9gXs1rWs89H7n9gOO\nB2P9N63fv07f3P7D1ey2vtmoxnRl/wBz+wHEA793shxCmfK7dNJ/2zf2Eh7KcQssUFOjL/uf2A4Y\nHdfspxBPHaaf9z+xEfZXXylyqzT5/wCT+wHDA9CvY3iT/q6b97+xVV7K6+3UWUxs0/NX1zJ4/wAF\nktHDA9E/YziSWXZpv3y+xRH2W1884so2/uf2IOIB237La9f1NP8Auf2Ifsxr1+Oj9z+wHFA7kfZX\nXymoRnQ2+nef2H1PsjxDTXumdumcljpJ/YDgAelr9huLThOTengoLfnlJf8Aoq/0fxHv/wA3Tdzr\n33v/ANAefA7UvZfXRim50Yf9z+xbX7IcRsr54z0/LlreT+wHAA7l3srr6VFys075vKT+xV/pzWL8\ndH7n9gOQB2n7Ma1R5u10/wC5/Yley+uf9Sj9z+wHEA7f+ltf/wD0o/c/sVz9ndZDrOl+kn9gOQB1\n4+zmslHmVlH7n9hV7P6xvHPTn/k/sBygOzH2Z1sullH7n9iP9N63OOej9z+wHHA6n8B1W/fp2/uf\n2D+Aat9J0/uf2A5YHVXANY38VX1f2Cz2f1lcHOUqsJN7Sf2Augu6vQYILur0GwZBHZmyMPMxnWrp\nzSpfIqKoRXgXZ7mMGZ5jLYthZthgQ8S6mzT6WeqrcIPosmLGW/I1aOyVdndlygK24S5ZLdbFrqUk\nproTfFOTbabJpknHGQgykhnicFkWcVJdS+utKC8QKo0ub36BOlRl0NMeXzCTjKQFChFx3Ww9VfLP\nuvYujFNY5diqf8qqbT36CrGDPPZZJ+LNNUe5l9EUwg4x36tNl0ZNVYwRqtEcQqz8ipWPcm2f8gp6\ntorKyb3Xigi1NPYmOHEtr5eXwAyqTjZ6FkJcz+bIsS7RroW01bJ4TAJJYTitxoQmovPQneM90XOS\n5MIDLPvZS8BYSxjJY08MWMHzbgTPlb+ZNeYvMXgWytxfMtx6oNrDAV2TlLvybZdpbZQvhlZWROz5\nUm+ucEqXLJY6gPzNtsKortFLO6J2x8yK31WN8gbYGLTd3jWoXnFM2xMcXy8cf91Z16/1K2z6MwLE\nbJYOhYjlaiaha/U5VV8llZ8RXLGzQsLM158SVJS6og6HDa4R1kbJy2UcojimrpsTVdKjKTz2j67F\nEXLT2Yr3bScS+XD3dp4OVsIzk3yrzKMkuIauzEZ2ZSLdOrNVZOEvxLLb2x8wo00NTU4wytRW94/m\nRcpxfF5VdE12f/QEaah1aiVdkYzlGLcPJ+RVS7q9bOvULld0W2vIs1M5aTV0RWZRhFLn89yL4Sru\nU733m5cjzu1gDn6mck+yksShs35mWcZNJrdIutdkpKU87d3cs06XM01kgSqtuO6Zc3nGDRFKEcbF\nNmWnjAQvOs4M1u72Lt0Z7JcrbwUNWsRb8iiFj7X7juaUHkpb5ms7IC+u3vvxK7pynPZYRCjyPMd0\n0Dg3FvOAKo5jnbJYoZ3WclSk4zwy97LKe4BCPmyL0/drd/wP/BOe6/MqlOUtNavBQf8AgDkQ+Feh\nJMV3F6Aw0g7Gm1UFp1GSOOaNPdCCams+RBv7CM1Ka6GeEe+OrW4d2WwsHhZZUWNbM1cP0vvU+VSS\na8zHFqfiWQU6nmDaYGvU1dlN1tp48ShLlLOWc6uee8mUvmWEVFqS8xq5WYaTWxVflxTisBRJ836A\nO7cPHiL2s8ldkn2q2LYYSyBdXqJRkl1yV6u7MOXze4tSzc35FLg7dTCL6NmWo1uKVVcV8SisodbR\nw/ApslmxuPnhF0toSQWqW3JOL6FceaMmslsN4pizxCSeMtdCsCEnnfwLozeMolYlUpv4ilSzLZeo\nGlKMt2tx4Zj06FNUlh56l8MYeWFNJqTH5MRyJHDHk3FJNAI/hy0K7FFrmQzxKPkK4KUlvlAO5KUP\nkPp052RhFbsqworBq0EnDWVzxsgEvh2cuWS3TKq1GVq5vh5tzbxFqeo58bGWEIu1JPGWAPEXjwYk\nJKLa8y6UErMeTKJ4U/mBvh0MNy5OM0P80GjbB91P5GLXbcR0cvng7dW05OhN7HH1UM3t/M7Mjm6i\nC7aT8zjVJVssAyFLlJjKMuoDRte2XiUdkzuXqmvSae66TxXhpLxZxNOufUOLSeYtL6HQ1Vqs4DV0\n5m0Bz67ZrXrUweHz5foW3/yOK2SfhPmRmpTcY8u7ZZ2c53RVmcy2TfiBp1FvbynWo80Obmi/GInu\n9lz5rZdPFl8Iwos5ZTwnHr/6K9ZbO9ctMXjoBVrZVOrsl+HD2M+mi+dMvfCdYqe07NYSzjO5nhd2\nSaxvy4IL3LLwyG8J4EhZzR5n1YspYQB2kSic1LyFsbfQzyUolRe4xsTiiicJKWMFlEmpZb6lsoqT\nYVni3GxJPbBcq3PfOxVKuMbM5wXJ80Ek8YIjE1i2SbyWwTxuJZBxm/E0VpSjiRRVurceAupsjCix\nR27j/wAF1jUemNkZrq1Zp7pSfSD/AMBXNh8C9CQj8EfQkKhoqswluWmfUvYCmzVTr2rk8Ce/3dG2\nI45FcAjbRxJRW+UzVDi8JYUng5HZoOzRR3ruM1KMVCecEPjFcvFHBlX5COsD0ceKVy2bWPUvr19O\ndsM8q4vwI766SYR696umb8B1qKmscyPHqy1dJseOov8AzMK9hHu0TtT6vCKdNY/eIt9ENGLr4dp6\n5PvYzIrUHCrmeySy2ZaiO2fbqEfM183clvucvTz5bnJ+JbfquVZT6lK6Na/lp+Q0q1NKXyOE+JT6\ndEPDjEoLG4ZdlxcY4T2JphKKz1OUuLRltJM008Qi1tIDUk+d5WC2KkkY1rFLx3L6dXFLcIuVkovo\nW33Pkisb4K1dXOSwWy5W1kKWvEoZF3c1joh4OKyl0Eb3eGBbzR5dyyEsMyNPlLFJpAaZvufqJW3G\n5NPdPIjsxXLK6MTtVzeWVsCNUnmTfjkWUeZvbOxXCzpkJ3OOWgrZV8CMXE3y36WX9+DZp3mpMx8Z\nWKqZ+ViOvV6jPLToz6HP1TxP1R0ZPMTDqYrKOdnyrL2b/QmPde6LI9EMkn1IGrxKa5NpPoaqK6rq\nlXNYUMPrsnkxRbpthLwTyXaiC7XUxzhNqSAa23S6Sc5Pl5u07qj1SKdRre3lVyQa5ZZT8Q02klr7\n+Tk38X4I9FpuGafTRWIKUl+JgcijhmounmSaT/FI7Gm0Venjld6XmzUKQK+h5zjelVOoVsViNnX1\nPSM53Gau10U34wfMgjzEm4x2IU9t2NjKkihpx8Qq1+DFnHuuT3+QtknyJJ4ZW5SxhsoamTbeYrBZ\nzPPTYqg+TbqWt/y8rqAs4d4pbnzYWdi1T55bDqUE8ERncZSecsups5XvHPqTzwWcIRybRRE2pyb2\nQl6S01qf5H/gVxlzZ6CaiMo0WN53g/8AAVgj8EfQkiPwR9CQqGUXl7M9oGdojA7FCIwTgkEAY2F5\nUOQAnIg5CwMAV8hdpKO01dUMdZIg6PB4qN0rpL4Vt6gdTU4nY8PaKK9ZLFEYLzJpasssf4UkTqYc\n2+fhWSNRz1hc3yRnuTlJLyRr7Npc2OpRbHv7FSqewT6iPTl7ckhoN4YRllW0uhKi1E0bMJR2wkBR\nDm82jR2k1WsSb3K+TCGl8Cz4AWU6i1TjuupfLiVkJtYyY6Vm2JNkczl6hWyHFMPdYyMuJRU13jmu\nHyBUZe4R3q9dCSW5f7yn0weejCUZJJmrT8yby3kDsy1C5WvMaMoza3XwnNcZYfe8B65SU01noB1o\nxi4xfiUW5WSuN0lCOzInZLl5mnj5og6eklmmJRxv/wDRz5STLdDJSpyhOLxzw635YZ16/US6bod6\nmL80jHqunoatJLm0tT/tRk1sXLptuZ5bpFEZDTuSxhblEcrO40vmjKroTViY+ot7W1NLflSZTDKi\n2ky7h0O24hXXJdXlgei4bpVptMljvy3ka2SDCFIJFYVDKroKcJRfRrBcyuREeMti67XB+DF5cmvj\nFXZa+W20tzLGai1lNlFV9clhorr5uV8yLrZuU8YwkSkpJp7AIkkl4izexdhRW26JnGM4JpbhWVcq\naa8Q8GNyONq2zgvjCM+qwEUOPcy11K45ztuabY4Si+iGSrgk47+YGd87eWivUZlp7OZ9IP8AwaJW\nJ7mXUzzp7f8Ai/8AAHOj8EfQkIruR9CQ0SRRYXyKLAKZEYGkLkICUQCAkCQwAInBAwEYOvpK1Vo4\nN9Z7nLjFykorq9jszrahCK3SSQGjTQ5dPJ43k8DXQzl+aDGK60vBMR2bPbclbVbKOPCJRZUuVST6\n+BZHvJpiyX0LGaz8mW0LFddyyaeX5ERW+AiIx26DOKW5YojRpnNYUdgrM0LZHc3w0Xiy/wB3j5Ac\nzTVvtVsO6ZOT2OrTTCMnjyFahF7rIRz1pLNnjYdaOb8jot8yXM1FIK4p5cZZA5stK4Nb7+RfChpZ\n6GqcMSUtmx60l3mAUaaHJ/Nl18DTDT0vdR6IhTrxnlzjwGdqnHEEksAaq4VKCSiuhXOMJVyjJLDF\nobdS3/UWyahtkirNHX2UHFPMc5RPEFnQXf8AEjSuTT5ljyLNYs6K7/izfD1EujcP30NP/FGXi1vu\n9XP4ZNHC3nh1L+RTxeCsoxJZWcmuyfapNOFHiacntsXPXdosPYzvRxy+VrAnYyXQ5q3V6nfHO8HV\n4HNPilTz1TODVTNzWEdThtfYa6ixtrEwr2zIJICFziW5HMs48SLU9nn1DZbIgkWSHEYRwvaGnMK7\nkundZw5y5KnLy3Z6ziVSu0dkceGTzNdatcq5dHFphWNayqccZQkLovL5jzV1dlV86+Zrlk0J2t8e\nk2Uewrsi1u1geM49Dx8dZqY/jLIcT1MPJgesebHiKSwWQWFv1PLw43bH4oFq47l5lFoI9DbW7Fhs\nRUpV4ycePHYeOS2PGqH+IDfZXKL2jsY7V/Jt/wCL/wADri9U4YU0VW6uDptxKLzF/wCArLHHZR88\nEMiL7kfQlhSMpmXSKZhFTIJZAAAAgAkgkoB4py6ComOckGjRL/8AIjJ9I7nVql2mH1edzmaVYrsf\n6HQ0O2GwrXN4jgztb/oaJrniyi3Z4+QVU5+ONxY2pJp758DZRonLvWbLwL/cqs7LDDLlPMukGW0a\nSc92mjpdnCHRJiuXkRSxorr67g7I5aB97dlMo4fdA0Rll9ULOayymvOSa0pTeU2vkEaNPZDv822E\nUu5Sn0wjQowipPG+Ohm5YuWcFFkYQmsjwXLlKLwJG1Q8DXVm3HLs/JijK231Q8VGdb8GhdVCdc8S\nWCdJXOx92Lkl1wgJT5cfI0wqXOrV0cXlGO9OTwnhDad3ReFnGPEDTCxyhHk280Z7W22pPoaqquVL\nLwF1VfNnqRTcPsc4NPflNlq5tNYvOLMOixCxxSwmjZKa7KXozXHaVVwhv+HVryyLxDLqeds7CcIt\nT0WG+kmiOI2p14W+507fVSacnl7ObWcltbj44K7ZQ508BHEpJLKOSr1JZXKi5SlHDxjG5XVBRa33\nLrZJY3yFezpn2lMJ/mimMzHwe3teG0y8lg2sIVrKwVwjhfNFgvRkRDIZJDArsWYteZ5S5S0+psit\nsNr9D1jPO8cq5b1NfjW4HjeLVcmvsfhLcwuJ1uMR70JePQ5mCqr5A7MsQwFHZoh1ryL8EYQGd1B2\nKLwAzOrHQWUJKL7zxg1OIs49yXoB0I/BH0DIR+CPoDClkUzLWVSArYo7EYQBnAEASSQAE5wTzbCj\nQXNNRXi8Abqu5TFee7N+m7taZknFKxQTykbVF92MVlvwCtlC5o+bZbXpoyasnu14GjS6bsNMpS3m\n/DyIlJVrDIFlsVSsipYfUWTcnzKWzfQluOcvdgVxqnN53x5lj7KKxJ5ZVO+a2Tx6CRjKc1nIGlVx\nlDYzXYhtjY0PurBktzZPHgEELKnnKa2IjKC+CbTHdPLVst2Ea8RWxVXVRjOMuWX1EdeJdU/1Hqrk\n6nyoiOnfVvcBVCeN1leDHrdnMnF4wXQhKMWsbfMtr00pYxlkFkboXVSq1MVjG0vFFWnnKmMqabGo\nSfxYNdXDueXeaXqzXVw2v8LlLHVRQyjjX0SqyuXL8zOrXzbvB3bdPyywk015mezh6uTylnzSwFYV\nPuLxItk87vBsfCpwWO1i15BZw29tdmk0COS75US7RZ6mHUcVusfZrEYt7nY13DrqKJSshiK32Z5m\n6D53lbN7G+NwVqr192nn3HmL6pln8RlbdGPm9zHbL+Xy8vTxE00efVQWd3JYOnL7fNR2JqP4Yt/M\nVTk5LY1djbXFppxfoU2KyLTUmv0OIsbnCKbwvUiNi88t+Ja9bZKEY2whJLbdFirourxXHkkFd32W\nsb0VlcusJ5+p2zzfsyp06y6qfSccr9D0bCIFe5LFZEBAAAhx+PR/lVy8ng7LOXxlRlpkpPGZbMK8\nZxWOYnGO5xJYWH54ONOPLNxKFRIoAAAAAAABAtnwS9GOJZ8EvRgb4/BH0BjQS5I7+AsgEZXNYLGV\ny3Cq2Ix2K0EKBOASAjA2CeUJAKaNFHm1MfqZjfwiPNqpfKDA0VRc9Sj0Ol00aVGySzP/AAcnhdal\nrG34Hf7OWFLwIpZWdzDMdsnLZs0WsyWbsBI7Npbj8uxFaaL1XlZAzzqaQ1fcb8yxflZbXS3PddAM\n81Jw+YtdO+WdGVKfyKZYWyAqcO588i8ixuWz+FCqK6gWwX8tKKwWVQSe+BHLlgkivtANklBbpZY0\nbcLbBjVqwQ558SDbLUYW8/8Ast0vE1p3JqKln5nGlycz5mOlX2Lks5yWDbqeIucpT6Z3wV069R5s\nyl8PgjDLL25SyqiT6rARts4rFPuwfTxGXFprl7qRnWkg45eSx0VpL5BU67VW3aKxY7jR5LVvE5Ye\nEeqtaemtj4cp5PWR7z3N8Jkqu2xdnj/sr083DVVSXhJGiOl5q4yzsNVVCNkG93zI6I9W7rH3nTKW\nxktlKzLencUdlWQkoNLGIrYW6tfFFvPXY4q43u7clHkmtvIdaSfK8ZfodF6fVybsi3NeTY9F7rtU\nbK8PyRBn4O3TxGlyzvmL+h6tmBqmfLLsowmnmLN2cxTKhWKS2KyIM5II6EN7hUSOPx5N6aMV55Op\nbdCGzl3vLxOZqX20sz6MDyevzNVvxzhnO4hX2d6/uimd7itMampR6NnI4w8zoeMdzAHOAAKAAAAA\nCAAWz4JejGFs+CXowN0fgXoD6BD4Y+gMKRiMdiSArkKMxWEGQTFABnIXOQAAR0+BLm1so+dbOajf\nwWzs+JV/PYDscOXLrLI+ayd7Ga474wcedXu/EYT/AAzOqt4rciq1FSzlCy0ufhaZdlJ4bQvPB+oF\nUdO49UxuzwtmO54xhg23uVFNlWd11Loyxv1yVzsaxsSpPClhehFNNtv7FXZPPQt94S6xFd6lssoB\nZ1S5VsipVyzv/kuk5cuX0K8ptJgNbFrbK6Gdxb8S65pyKtwhHHbdjcryDx45G5kBU4Lm3eToafTT\ns02YQ5o5MKw2dLR65aWl1uLafkBVyxhLOFsy6xQnNSg8Ra3+RgtnPvSTGrstlXbsmoRTfzCtsOWS\nxzRCxRr5W3FmONlaypwZNs6njDaSAfU/7U5KLUGmk/NnkdXPnk8Hq7tZ22mjTKXdj026HmtRRCMn\nJTzv0wduvGEGmlN8kGnjA+FGUPk8/MWiTdqa6NYLfdrVZnrnZNCj0elufYJt5eNiyOtuhZySqRnr\nVem0EarFm7rnyMktS44aba8jlR3o8TvjjFCS+bwZ9RrJ3y71defDGMo5M7ouzCjLpnr0M83ZzZh1\nIrtfxR1dybT9UdbhXEFqIOE2k18O/U8Y+0c4yk2zVDVyhLKfK10wB7iTFZ5yj2k7OKV8ed+a2I1H\ntPlONFOH5yYR6CyyFScrJxil4tnG1vHVHMNJHnfjP7HC1Gr1GpsjK+cpKXReBdoI519UWs43x5gd\nLh9N1l8brpSfMnht5ybNXjlW2GngnTYhCxOO8Nsp/QzXXSnzc2cb4QVx+JzU+VeUji8a2spi+vId\n+WlWo7KX5W+Y8/x6alxGSXSMUhEc4AAoAIyAE5IAgCci2f7cvRkkT+CXowN8fgXoDCPwL0IYUshG\nMxGAssCMZisIUCSGAAAABbRZ2N0LF+GWSogD33LXqaIt7prKfkTW5KCTefmvE53s/q1foezk8zht\nj5HUrqjFNL4Xul5EUOMZJYe/kJJ1wfTcuwl12+aFlWpPK3Az9rByNFU4Nddyq6mCaykPGMUly4As\ns5W+hWiWxc7gHIpbMbsYrwFzuMp5eACUMehmlHNixt6m1LBOIyaykEc62OLe/lehE5RSxFs23Ux5\n28FE9LneJRkU3ncsVqxuskTonHZora5XuQXJ1Po2n80M/VMzdehdCahH4VJvxCpsmvFD1W1pWx5s\nOccGebbmttglhT3XgEXQSbeJIa2puHRMxy7rym0NOcuTuyYUt8GlFcvXyC7hyhFczbbRV21y6S5s\neZojxRWJR1FeMLaUTUGSOlUH8HQ00p1ruRSbJlqKmv5clv8APARnyxy19NxaHlC+fhn9Spuda5ZL\n9GjVF2yg5xjiC8X4lFtdtj5mtiCqFn8zmcV6ETvw3hYHVE4v4Sq2qSW8WQI7JPoJbzRWX4k8uFvk\nSe6wBW5po38JdT19asipwkmmn6GFxwmiynmjOMoZUl0aKjvdppKIaeVd0UsSjGUllx36NGCUfetU\n56eLjmT/AFDS8KutvSvUq0+9mS6nc0+np0lcUknPmxzeLIIpUqdJyvpJ7sexQxnPd3Rdp9JLUJx+\nGKl1DUaePbdnDourIrBJRqolZ0guiPCaq13amyx/iZ6z2l1saNL2Nb3l3UeOLEKSDFKJyBGCQIAk\nAAWfwS9GMRPHZS33wBtj8C9AIi+4vQGwpWIx30EYQrEGZACrqTKONs5IAogAYYIIAlkAbOF6uWk1\nSecRltI9rp5q2C5euM+p89PRcC4ltGqcu9Hp80Qenim0RiMtpLfzRq0XJOptzXeXQqnypgYbNPPm\n2eR4VuuOWXyjlPHgVOWwUjeRZFnNGxZxghwXgwK09xpPG66k9mwdbaCK1qHnDHjd3lkqnTvvlEKH\n6gWW3NSfkV+8rLwRNPHmJ7u5PK2KHnqHjrn1Mdtj5u8jZHTpbvcd01tbxTA5q3ezwWKUo7Pf5mmW\nmh4LAiplHdAVStSnuiXJTeU/Apsli7lkhe0rjb1a28ALZLKytsGedkknncvralF8s0/kxsuvL7JP\n54A5vb9d8ExsblusrzR1FpXdBOdMcM0UaCiFkZ9nh+SCuVGKfWDS+aLUkl3VsdiylYzKKQnu8ZQw\ntgOdPU2ckYJ4jHoiYa2xLDSfqTPRWu7Eej8WZLIThKSfgBtetUvDAvvEZdZfoc9czeOZBKMkuhBq\nnPK2a3F5orCwmzPVTKabXgbNNoJTeZPARr4do4aib7SOVy5SXVs6Nelqrrq/lRhyWd543fkVabs6\nUsLn5fmPqe0nRZNvCe4FrujpkoSnzTcm8LcrhKU5RT88mLT04xOfU6+joUouctkiK3UydVax06sx\n6q5U0znJ4b3ZrlOMK+XxPH+0/FMZ09cszfxY8AOJxXWPV6uUs91bIwkEo0iWthSQAgGSQwAgAAlE\nT/25ejJFs/25ejA2x+BegMiPwr0JYCsRjiMBWKPtncWWM7dAFAAKIIJYEEEDEYAgsrk67FOLw10F\nSGSA9Twniatgot8sl1WTpK1uWct/qeMpcq5qcXho72h4hGxqNm0v8kV3Fe4vKWU+okr4SWML9CuE\nk0Sq6+bPLuA1cqZS5XLlkWdk85TyiidKlJyjs2RF2Q2UmgL+zy8KW5DrmujyVdpLm67jq9RliXXz\nAnf8SFaiWuyLxl5T8SXVCSyn18mEZ3EFFrO5eqljZ5EcflgCvp1Fc14FvJkmdLjFScdmFZufvdBu\nbubD7J7xGbra+FblRznCTlzSjnIr0sJTzjc3JRWcrYRJKS8yZVOn0dT+KCZqjRWvwlLlKKzEftZJ\nZYyNXd5PhwirMYPbBmd1jXVlMnJ53A6F90Jd2Hh4iVtZyzBTGbn8joVw7u4ESks5wc/W6eE5c8W4\nvH1OjOvbKM06lKSy8JAcf3eUm+73vNGyjR8yzN/obuzikngWTUU30SCCNNVcsxisGW7V12SdS544\nfWJRqNU2/wCW3uRp6JXXRrgt57gdTh1DWmct3GUvEtss58rPcXgbYVLT6TkeFhGWupJbolUaehz7\n8+7BeLOjUm15RXQyxk1CPaNKEP0RweN+0vKpafRy36OYGvj/AByrSVypokpXvy/CeJsnK2cpzblK\nTy2yJTlZNylJuT6tkpbFQuCScAUQAAAEEkMAAAABbP8Abl6DC2fBL0YGyL7i9AZC+FegBUCslisC\nBWSyHnGWEQBDBFEgAAQCQEogaKHiiEOgGWxLk1uhGyMgdTRcWlU1C7ePn5Hd0+qruipQmmjx2R6r\n7KZc1c3F/Iivbxlknqef0fHN1G9Y/uR16NXXdHMJxfowi2cN+aOz/wAkOtS69RlInIGbknW2luhL\nZ2xWYJ/obHuQ+XGHFAc6Oture6aNMOKr+pHJdiD2cER2NL+KOPRAUvi0ObCiWW8W7WmEOXeKwUWa\nStzzFELT56SxHySCj3mT6lsG5JMWFEY/P1LVttgZD4I5VkOYMkEvYVh4gEQSkvICGUPFqLyVajVS\nra5VlMZkYT6oAp1E54ysIt6lXNFLoJZfiO3/AERVs7VDGeniZrbeeSx8OSmTnZ8WfQsjHdLH6Ioo\n9255vlOrpY16RxlHv2YafkhIVwis2SUF5ZFt4todHHOeeS+eAjo1wu1L5p5x/wBGbiHFNJw9PvKy\na8EzznEvae/Up11Pkh5R2ODZbKyWZybfzGFdXinHdRrm455KvCMTk5yxckoodLAyFiOugQCjEZ2A\nUAAAAMkMAAAAgif+3L0Yws/9uXowNcfhXoDIj8K9AYCsVjMgBSG8ksgCAACiAJIABkQSiCxDttpZ\n8BIjgQxcDAApGRmKwAau2dUs1ycX8hAA6dHGr69rMTX/AGdOjjWns2k3B/M8yTkD2leprsWYTT/U\nbnPFRnKDzGTXozXVxHUw/qN+pB6rmDmOBXxm1bThF+hfDjNf44NAdfK8yub5f/hz3xPSz6za/Qsr\n1mml/Wjj1IroR70EwZXp3TYttTX9R7Yxh/Vg/SSAjO7Y0Xl4Mrab+NP0ZZHEV1/7INOEurQu3mjI\n7It/EvqHb1rrZH6lGvbzIbRmeq00Vl3RKZcT0sOtiYRsnLyI76XXJzJcb06e2X+hVZx2Ge5XJ+pR\n1eVvqxW4x6vLOFbxi+b7qjEx267UWfFY/wBBhXortXVX8UkjFdxuFe1SbfyODKTfVtisI3ajimov\n/FhGOUnJ5k2/UUOpQATgnAC4GSDBIEodPYRDICRRhWBAAAEADBAQSBKy9gIFn/ty9GM0LP8A25ej\nA1R+FegMI/CvQGArIJZAEPoKS2QBAAQ3kAJIACSUKSgLYjlSHQDBh46D01yusjXBNyk8JI9DrNLR\nwbRvTShXbdYu/Lq//gHmmKbFTXGvmt5u9nlwZZAIGCQAXAYGDAC4JJwGABSYOQYIwQLJiNjtCuJR\nHO10bDtZrpJ/UjlI5QH7axdJyX6k+8Wv+pL6lWGSkwH7WfjKX1I55ebFwwwwJ5n5gmRyk8oBknIc\noYAlPYVlk5Oc3JpZfkLgBMBgfBGAFwSiUicAKSTgMARgCWGAJQ3gQiQIIZLFACMksUCQIJAAAgAI\nn8EvRkiz+CXowNUfhXoDCPwL0BgKKx2KwFACGAMXJLFAnIEEgAyFQyAZDp4EQ8FzPC3YHo+E1w4V\no1xW9KVkk1TB/wCWcrUameosd1zc7JPL8jRrdRPUQhCycYcsfhb6HOnPwX1CpsslOWZP9PIqZJDC\nIAkAIGIwAEkAAEADIAGQwACMEEhgBSUGAAgCQwBBJKQYAgCcMlICAwNgMALgMDEAQlhgyQAjAYJJ\nAXAYGABUSAADEGZAEEEsgAIyAAAAAEET+CXoySJ/BL0A1R+FehDOetdaljEPoHv1vlD6BW8Vsw+/\nW+UPoR77Z5R+gRtAxe+WeUfoHvlnlH6AbGQY/e7PKP0D3uzyiFbAMfvdnlEPep+UfoBsJyYvep+U\nfoHvU/KP0CNyLabOzlzJJv5nN97s8o/QPfLPKP0A6lt0rZZlj9BcnO99t8o/QPfrfKH0A6WQOd7/\nAG+UPoHv9v5YfRgdIDm+/wBv5YfRh7/b+WH0YHSwBzf4hb+WH0ZP8Qt/LD6Mo6AHO/iFv5YfRh7/\nAG/lh9CDoYbQYMH8RuxjEPoyPf7fyw+jA6GAwc/3+38sPow9/t/LD6MDfgDB7/b+WH0Ye/W/lh9G\nBvFwYvfrfyw+hHv1v5YfQDdglIwe/W/lh9A9+t/LD6MDoYDBg9/t/LD6B7/b+WH0A6GCcHO9/t/L\nD6E/xC38sPowOhgDn/xC38sPoyP4hb+WH0YHQIMHv9v5YfQPf7fyw+gG8Dn+/W+UPoHv1vlD6AdA\nDn+/W+UPoHv1vlD6AdADn+/W/lh9A9+t/LD6AdADn+/W+UPoHv1vlD6Ab2RsYffbPKP0I98s8o/Q\nDcyGYvfLPKP0D3yzyj9ANpBj98s8o/Qj3uzyj9CjaRkx+9WeUQ96n5RA1iz+CXoZvep+USHqZtNY\njuBSAAFAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA\nAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA\nAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA\nAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAB\n/9k=\n",
      "text/html": [
       "\n",
       "        <iframe\n",
       "            width=\"560\"\n",
       "            height=\"315\"\n",
       "            src=\"https://www.youtube.com/embed/BLRqpJN9zeA?rel=0\"\n",
       "            frameborder=\"0\"\n",
       "            allowfullscreen\n",
       "        ></iframe>\n",
       "        "
      ],
      "text/plain": [
       "<IPython.lib.display.YouTubeVideo at 0x7f495c0eef60>"
      ]
     },
     "execution_count": 1,
     "metadata": {
      "tags": []
     },
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from IPython.display import YouTubeVideo\n",
    "YouTubeVideo(\"BLRqpJN9zeA\",560,315,rel=0)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "t_SvvZWjp4qc"
   },
   "source": [
    "### Critical Point"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/",
     "height": 336
    },
    "colab_type": "code",
    "executionInfo": {
     "elapsed": 369,
     "status": "ok",
     "timestamp": 1539185457464,
     "user": {
      "displayName": "Jeffrey Kantor",
      "photoUrl": "",
      "userId": "09038942003589296665"
     },
     "user_tz": 240
    },
    "id": "sP_C_vKFp4qe",
    "outputId": "88fabe2b-4170-4e41-bcc8-9b3844543b6f"
   },
   "outputs": [
    {
     "data": {
      "image/jpeg": "/9j/4AAQSkZJRgABAQAAAQABAAD/2wCEABALDA4MChAODQ4SERATGCgaGBYWGDEjJR0oOjM9PDkz\nODdASFxOQERXRTc4UG1RV19iZ2hnPk1xeXBkeFxlZ2MBERISGBUYLxoaL2NCOEJjY2NjY2NjY2Nj\nY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY//AABEIAWgB4AMBIgACEQED\nEQH/xAAbAAABBQEBAAAAAAAAAAAAAAAAAQIDBAUGB//EAEgQAAEDAgQDBAUJBgUDAwUAAAEAAgME\nEQUSITETQVEGImFxFDKBkZIXQlJUobHB0dIVIzNicuE0Q1OC8AcW8USisiQ1Y3PC/8QAGQEBAQEB\nAQEAAAAAAAAAAAAAAAECAwQF/8QAJREBAQEBAAICAgICAwEAAAAAAAERAhIhAzFBUQQiMmEUcYET\n/9oADAMBAAIRAxEAPwDz9CEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEI\nBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEI\nBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEI\nBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCEIBCfwneCOGeoQMQn8M+COGfBAxCfwz4JMh8EDUJ2Q+CcI\nifnNHmgjQp202b/Oi9pP5J4oZD6skZ8nIKqFa/Z8/wDKkNFKNy37UFZCseiSdW+9J6JJ1aggQrBo\n5Bzak9EkHNqCBCuxYc+UDLPDnJtkJIP3WWlB2PxSf1eCPNx/JBgIXVt/6f4u7/OpB5vd+lL8nuLf\nWKL43/pVyjk0LrPk9xb6xRfG/wDSl+T3FvrFF8b/ANKZRySF1p/6eYva/pFF8b/0rPxDsnieHtzT\niIt6tcSPuUXKwkKd1LI020SejP6tRfGoUKf0STq1SR4fLI8NDmAnqT+SmnjVRC0xgdSTYSQ+8/kp\nW9m6x20kHxH8k2L4dfpjoW8OyVef82m+J35J3/Z+If61L8Tv0psTw6/Tn0Lom9jMRcbCel+N36U5\n3YnEmtJ41Lp/O79KbF8Ov05tC1T2frePwWmN7/Am3vIVr/tDEc2USUhf9ATAkexVnMYCFvnshiIO\nUSU5kH+XmId7iNfYq1L2eq6mokg4sEUrBfJI4gu8tEMZKFsP7N1jN5IPiP5KF2C1Ld3xe8/kiYzU\nK+cJnHz4vefyR+yJ/pxe8/ki4oIV/wDZM/04vefyR+yZ/pxe8/khlUEK/wDsif6cXvP5KWLAqiU2\nbNAD4uP5IZWWhbY7L1pNuLTjzcfyVWowarppAybIwfS1I+5ExnIXSwdicQqIxJFVUTmnUHO79Kp4\nj2YxHDSOOIy07PYSR9yauVjIV1uGTO2fH7z+SkbgtQ7aSL3n8lNjU+Pq/hnIWuzs7Vv2lg+I/kr9\nJ2GxKrjL4qijsOr3fpTS8dc/ccyhdb8nmLfWKL43/pR8nmLfWKL43/pVZxySF1vyeYv9Yovjf+lH\nyeYv9Yovjf8ApQyuSQut+TvF/rFF8b/0pfk6xf6zQ/G/9KajDMgJ1ji+1ITGf8gex66ipwvDqKk4\n1U1wbsA3clZVsEkv36mL+poP3IMzhxH5srfZcJvChO04Hg5tlqigwyTWLEmtP87S1H7He/8AgVsE\nvlKCqMr0W/qyxu8nJpo5/oX8iFpSYHXtF+AHjqAD9yrSUFVCO9TvaPC4QUnQyt3jd7kwgjkrn70c\n5R7bo40o3ff+pigpaoVviA7shd9iLRuGtOf9rroKzZHt9VxHkVIKmYbSFScODmJGebUnAhd6tQ0f\n1aIG+mS8w0+bU4VgPrRD2FJ6I4+pIx3kU00kw+bf2oJxVQHdsjfcU+7X/wAN4f4HQqk6GVu8brdb\nJmo8Cgtu0Oy3ez/aiowt7YZyZqXbKfWZ5H8FzrKgkASa+KcRzGoQezUtRDV07KinkEkbxcEKRebd\nksddhlYIZnXpZjZwJ9Q/SXpO4C3KUIQhVC3TJYo5ozHK3M124KchFcfjPZYtzTUfebzbzC5SamdE\n4hwII6r1vVcv2j7PumcailZcn1mhY65dee/24kR25p7CGPDhyT5YHQvLXjKRyKZ3fpLm6tRh1BHN\naNMbrGppA9gb9FaEL7KOs9xrcQNRnkeO6NUyji4+t7rXigZE3YXVkZ3GSyWaGVvFYcp5rSMjREXE\n6AJ0zBK3KRos6pzxwyRHmCAlhawq2oiqpTFw3tc192kHQ+YVqbDmVMTXTSZXn517KKCqEtNAZYST\nGS0ybkm+xUdXUZpy3VoHI6LvzzMeHq3yWqOslgkdhWKyG7RmpqjmPC/RNragVTWTOs2qhOrh84Kh\nUTGWrps50jbuVA6oGeSxvdZ6mVqVuNc2ohzEm5WbOLPIVikktA3Xkq1QbvuFzrSFyZcpz9lESopx\ndqlzKO6LqqkzJzS6+l1GOq1sLwyprXhzIrMA3NgqI6Wskjs2Rpc1a8bo5o9hIw7gqGowaSNnezRk\nbE7Kux09C4CdvdPzgP8Al01ldp6aSgvLQOLoybvgd94WvS1MGIwFjmhwJIex41HmFRp5Q5ocw+5T\nSQCa09O4Q1LRqRzH4pRjYz2ckp71FAC+I6mPm1Yccp5rvqCuM5dHK3hzsAGUnQ+I8Fk452c4zzU0\nEYD/AJ8Q5+IWcd/j7s+2JTyEkLr+zodwHOOxKw6Ds7WvkHGZwmjcldfS0zaaFsbNgEka+XuWYlsj\nRKEclt5iITU4bJecJS20KaTYpyZzV5iWuJx6MVMlNSnYkuKh/wC3KRwAu8HqCrLzx8cfzETQ327r\nTACyza56TstHb93UOH9TbqtJ2ZnbqyWM+dwurWdLjNJFJlc5xbctzBtxcbqxNc+cFxOHWIn/AGSW\nSXx2n+fU28e8uppaynqxeCTNvoQR96nJaNyBfqU/7NcYcXxFmlQyKT/9kIuj9rRv0nw2nd4tu1dm\nWBw1FwoZKCll9enjP+0Ia5L0nCZPXopoz/I8O+9JwcGl9Woli8Hx3+5dJLgVA8fwMt/okhVX9mKR\n2rZJW+0H8EVjDCqWT+BikHgC4tP2pTgVaR+6fHMPBwcr8nZX/TqR/uYqsnZqtYbxmN3k6yGxRlwa\nuj9elv8A7fyVd1NNF60D2eTiFqeiY3S6sNQAPovv+KQ4pjMP8UucB/qRf2QZJdIG6ukA/m1CvUz9\nnOY0ttle0i4IXSYfFHiuFtnqKeNshJBytsDYrm5LRySNZ6mYgeSQZ+IUno1S9rdWbt8jsoI5LaHZ\nXK2YyvF/mtDfcqJFignabFeldjcW9Pw4U0rrz04tr85vL3bLzKN3IrUwfEZMMrY6qLXI7vD6TTuE\nV62UKKnqI6unjqIXZo5GhzT4KRdIlCLoQqgQhIorz7tbBwsTebb6rAaLFdj21htJHLbcLkRrouV+\n3p49xPAS1wIWrF347jcBZMQ1AWpTuMDm8xzWcddxqYbKWkLbY7MAsCnGWS7CCDtZbEEpDe8FZWVx\njQqeIQh0TjspDVtjveyzK/Ec4IBAFktX7YUGIvoKuVrImTMk0lgf87xBW23GcOqaThy0M5la3QSs\nBt7VylawTTF99bqq+eqjGX0h2Xa2danVjz/JzNXJZmtmke43J0b4BVA7vF1/NQ+MkjfIG5Ub5c12\nsFgpusLvp5Jysl4YA0uL3KjlxN72FuUA9QVTbC957rSSei0oOzuIzNDmwd08y4D8VcTVQV0jdibH\nkdUNrDoTfxC0D2YxEf5Q+IKGXs9iMYuadxH8uqYbTI6yN+92lWGODtje6y5qWaB1pI3MPQiyKdzm\nTs1PrBTFnTYjLA8cQkNJtp962cPfKSWcVzmX08gsihnjbVtiqGh0T7gm2rTbQ+9aTZWw0Qa24cND\n4Hotcz9p1WnLVVFE01DM8kI/iMLsw9xSVFRBLC2SMg00+zf9MqhRV3Cc++sbmkOHVVIXCOhlANrO\n0CvUhFuhqnU9QYybtW3FMHgOaVyrZL1GYbWVyKvMMwN7t5hc9V0j2tqWtkacszNWu8Vboqt8oyym\n0rPW6f3CyYprhr2G4OquOdnaJWeu3dt7X8FbPysrbY7MOh5jonXWfS1QfZ1xrqddx+av8lrm7MOj\nkJL6pbKYBCEKKQlFgUu6Nld9J+XEYYOJPNPe+d5K1As7CyxtMxoIuBrr03WgEcqZUPeynkdG0ueG\nnKBzKz8Lw0Q0WSqjY97nFxDgHWUXaB7nshpmOyl7szj0A6qevr2QYc99NMx7xZrSHB3tWoitgjDJ\nWVdSW5WgmONlrZRfa3uUcjDiuM1EDnO4EceQ25G4/H7loU8j6XCBPO4ueGZ3Fx1vvZJg0r5qITPj\nYxz3EksbbNruhrMdTHDsTp6Wlme6SQguLvmtv/YrozqNNFlUUsGIYnLVMhIMQyB9/W35KSqxJ4ll\np6OLi1Edi4O2t4dVKKr6rGKQF1SITA06yEA6eQIWtSVIqKNlQ5uQObmIPILIm/aGLU5gkp/R8rgc\nztA7wWjiEb2YU+Knb83J5N5/YlE9PUR1UDJoj3HbJwex0jo2vBewAub0uszDq6jp6NkDZQZWjbKR\nmcl7PRPEE1TLcvmfud7D/hTGpWoQeSTLfcJyUbrLVMqXtpqGRwFgxpP4rgybBdh2jm4WFPbzkIaF\nxkjrNPkrGVOQ3cSonqQ6lMk0CKRujgrTTaM+YUT2gBT08ZmlihG8jw0e0oPVez8Ho2BUcR34YcfM\n6/ir505qBrwxrWt2aLBLxVqemrEwuUE6qLjHqo3S+Kup4rOnVISFV41kx0ynlWsjP7VQ8fDw4C+U\nrgLWK9HqyJqZ7DzC8+nhcKnINybBZ6dOLixQxF15CNG7LQiZnuXbW1QyHhRsiHLfzT5zwoQwbu1K\njHfW1QM8kTrNJLRspmYnIBbMR7VXl0F1VebFZx0ny/tpOr3HUuKqz1ecHVUi8qrUT5BYHUpOS/L+\nj6iry6N36qiXFxudSmnUq3S0plNzstOFtqOGB8ps0LUpcNY3WU3PRTQRBgysC0KenDtZHADxKLhK\naEMNoox7AtaGKqAHdLVNRSQRxC1gOpGqt+lwdT7lZyWqktQKXKJ3lpOxQyup3bVDb+OizcaqBNPZ\nvqjZZmXVRp1ErI547PZHKzxAKyKvs5R1FzA/0eTpu3+yqwPkjN2PLfIrVp6l8tmvaCdrhJSxzddR\nS0c4jqAGOOocDcewrp8KxGklaxlfDkkd3ZXWux+mjvPxVbFYGzQ8KZpuPVPMLAdJVUQLBI4x8uis\nuMWOjxymomZf2eBc3zdAubmmyt4bToDqRzVeetmm9eVxHRNjbnGa+iW2rItQvIanZrlRXsLID1lW\nzhNZleIXnuu2W/TyZJLP22K4lsha4OabELqqOo9Jpo5eYFnLUGqX+jTi9iwna261ad925SQeYPgs\nYOMsWUnUDQqxQzgC19WHW+5U+rq/bWRqENsRdLdb1nMASpNylWeiBNO6ckI1Tn7K4ptAxsZHGymS\n7dbc9SPaQpTSSiRzg4HNfmRYm4v7svuUlRTcSanL2v8A3bi4chfxVlHJBU0kNUwMmZmb52VRuC0j\nCbB+U2u0u0Ot1pJA5pcWhwJG4uqIKymFVTOhJsDb7DdZsdHikL2tbUsMN7Fg0sPDRbSapLgp4TRm\ngoxE4gvJzOttdZUVTNhtVNJNSyyySm7nAEAa7DTVdFZCuhlJOKmmjmDS0PF7Hkpr6JnPRKdVFxkY\n42GlgY+CBjZpH2a5rRcLWgYWU8bHOLi1oBJNySo5qeKd0Zkbm4bszfNTDRNQtk5nrJvJOiFyVGvd\njB7WTaU8PiXH7vxK5ac6WWz2km4mJubfSNob+P4rClcb7qz6IiG6a4XkaPFOCSPWbyCKkdq5anZ2\nETY5SAi4Y7OfYL/ksrmuh7Ix3xCab6Edh7T/AGRY7YyI4pVbOm8RVpa4vimulVUvKaZFBZMvionT\neKrmTxUZkQWHTFYcdNxMTfIRoz71pGTRVKmqjoojIQC950b1UNPDLvVWp7zyVVvi9YOJHHMyI82M\nIHvViHDK6UWjxGKST/RdLr7ijOqsoVSRgV6WOWKQxVEeSQbhRSQncagoMyf900u9yzHOLjcq7iT/\nAN5wxpbdUmi5VRNTQmV4AW1DEGtDW7BQ0cHCjBt3itFrBGwdeajUIA2MeKmppP3rS7UX2Kpufmdo\npYnWI1UdHaQCExNLWCxCmDWHkFk4TUZ4shOoWm03C3rLl8XaBXSAC2qptatrG6UiQTDZ26y2sWas\nIxuq1cMicZc4ZmyqnBCXuAaLkrpaCjdTxC1sx1KkWoqqOnrYeHKMj+RPJcvVUrqed0UzQRzHIhds\n6NrhaRg81RrsKjqYHhh7wF236rrZM9OevOcXoDRyNfHd0EmrHH7R7FUppsj8p2K6V7WyQzUlSDlI\n0P0Hcj+C5WVhjkLTuCsH00TukumU7+LF/MNCnFRTw5bOAVGWcwOOj9vNYgVinkMM7JAdWm6DtYrt\nPkmlxp60OGgKcx4ka2RtrOF0lW3NA19gcpVI36VwLS297c1MQszDKgyNbc3IGVagU53CkCVIQg7L\nWJpbISBLopftY5CPHaEi/GLP6mkKzHiFNL6k8b/aFwOV7b2zBALj4+Yurjnj0QOjeLhrSPBI2KEO\nLw0hx31Xn7Jpoz3HuZ/SSPuVuPFa6OwbUyeRN/vUwyO2LGHZx9oTXRmxylt+VzZcqztBXsIzOY8e\nLVai7TPuBJTsPk4j8ExcjejhmaDns7+lI4EbghZbO0lMdHxSt8rH8VZj7QUTjYyPYP5mkISRZe8R\nszHZJHIJGhzb2KWPFKKX1KmJ3mfzUrHRPHcyEfylE8TEqlyNPI+9Ry0zZMv7xzS250HhbVTUsBJB\nTmusw6boEOVoGe9ha5UNdKKeilf9FhI9yEjiMQm41ZPJf1nm3kqEpFwAp5OirSm7zbZaaINkkXzn\nexBNmlLHpEPE3UD27rqeyjAykmlI1e+3sA/uuWbsV12DN4WGQi24Lveixr59U0v1UHEKXOmtJg/R\nRufqmOconOUEpco3OUbnphf1QPc/RYpqnzVuZgaZL90uFwweS03uu0jwWThzadj6hz5CKhrhw25b\ngjnqrJqdVv0tLi9R+9hxKbO3Zj5C1p9gSEuqpf2djT4mVzv4FQ0b+BI0IUVZVtMDGxPIadwsutL2\nU8dVZjsj8jWu8Qul5kc57XjUiSOWir7+kw6RyHfTl5JlO8OizaDqOhVfE3xMxNro3Z7tBc7qVC+Y\nxYdUOGhJsPaudbYNVIZal7zzKsYbTmaoaLaDUqnzW/gcWWN7+ugQXWQ2dmNrNUdS+zfEq4fUCzqq\n+fwCjURtKmY5V27qViNxp0NQYnghbsVYHAX0XMRGxWnBJmbZCxtl7JGZXAOaeSrfs2mc/MC5vhdV\nWOI0BspmTuGm6MtClp4Kf+GNepVxsiy2VA5qdk7TsVTWgXg7ppaN2lVg+6lbIrKma5ztHScGdtQ0\naSb+a4/ForTCQDR4v7V6PjMQqMOkB3Z3guFxOLNSHTWM3USsmiflmDTs7RXXixWbG7JI13QrXeLt\nuNiFCIAnhIAntCK6vBJOLhjL+tGS1aJbnge3wusTsy/u1EZ8CuhhbcHxBWk3FbCpcsmS+5Btyuui\nBsuZoxlq3M10HJbz6uGBrOM9rS4aXU5/yxb9LN0DdMY9kmrHAjqnhpG616jP2ckISoWZcax5NQUb\n62VzA/JlaSCRfXorr8KqA/JHOx5sOfImw++6ZRU8UlOHvbmcSdcxFlLUDgU7nxzTMLQLDPcLWubK\nLraFrTbwUjYZHsDxTvLSbBwBsVEdVs07qyCnEbY4ntLbCxAc3c8/NBkFo2vbzF0mQfNcPatWtrCK\ncxPpOEXjKLG4sLfd+KzISwTM41+HmGa3TmgYYzuLOHgUZXtNg13sC1J4MPkgc6mcBJlzBpfbppqf\nMqnPTmBjCCSToel9Nvf9iYKhkINj7jqntqHN2Av/AC6fcrtNQVNVDxGEBoNhmvr5JJcOqow7NAxw\nabEgDxP4FMEMWJ1MR7s8tvB/53VpnaCsZ/mkj+cA/kqUtOYTaWBzNSL6jUaFRFjDsXD7VMG2ztPP\n89kZ94TMRxr02idG0Nbe1+/+Cx+E2+jxbxFk2SEtbmBBHgUwQvdqoDckkqV50UPNRTZPUUlrNA6B\nMfq9rfFSHUoHtaXWaNyV2kUfCibGNmtAWF2dpYZ5pXTBrjGAWgn7V0gdHexe33plv0sRIup+GDsk\nMVwplaQXTCVK6E30UTmkaIpl7nVNsEOsml1iEEkUd5WgrGqY5MNxJ7mjM1wPd6tK26Z152q7WUcV\nZDkddp5PG4TcZvtSdWYI7D299xe1v8N0Lg4+GYaLCqDxqhspY6KmYe41+5VqrgxiheW05bLHyIY1\nxHvF1mTRYnVS5p45nu/mbYLV6tZzA9+eRz3czskxBxbQtbf1nXKswYe9hzTnb5oUOLgcBttgVj8t\nMhvrBdJhP+EH9S5xp7wXR4Qf/ox4ErSRp5NB5KlUMuHLUZ3omHwVCe3FcPFZajNA1UrAkkbaQp7B\nayrcSsBVqB2UqBoUrNNVGlxrz1UjXlVc2gKVr/FNYxcD/YpGv9qqtdontIv0VRcZO5p0Jsp46zXv\nAHyWeL8jdLmtuERpvnjkhe0m12kWK5SeMPjnZ1YfuWzmBae9y5rIcLuk/od9yQck8WcVs0rTJRxn\nwsoYKFhmDqgnKfmt39q3aR0EEPDio2vbfncqsxkcNOaxbDI6SUubNAYnHYs5KKow59OQ4EPjd6rx\nzRU3Z3u1kjfpM/FdNT+sPaudwZhZXXP0SuggJvfoiqcH/wBxdqRodgug4EUrAHxtNhsQueppLYo8\nh1u6faujZ6g8lJ/kv4Un4WGuzUz3Qn+Q6e5NE1dSazM40Y3c0aj2LTGqHAEZTzXXzv1WEVPUsnjz\nsIIUw2WfPTCmc6oidkv6wGxV2N5fE13ULPWb6PbxuWOUODmOLNOWl03iVZGWSVzmdCbrXnxKCeF7\nfRWNcQA1xscu9zt0T3PwmQ3IezRtrNsNCL7HmPuVyMswOLWhWRjpa60tP5lrlO+PDzE17Xuvkc4t\nzga30FtUgwps087aaZsjYhuRbMd7D2KWCrUYg2sfHw2ua1u91JQyRx1WeZwaANCeqnOBThx4T4ZB\nmygtdvpdV5sMqg0tdA462OXW3u8wmehqtfSzG9oHeIAWTIYzUyiMWaDoAqkmGyxXLmSMt9JpCdEM\ngtfVPY2qQzVFOM1TKxre6wN0AUzvS4WF7azNZuz2Dl46rBNVVwD9xIcvTcIdilbJGYpALO3OWxU9\ni1VVktY1hla0ZSSMotublW8NJjiJdRmdjnA5t9tNvesljiG3vr4q9DjMEEbYpWPBaLXGoKu4Lr30\nXEbJPC5vcs4PhtmNxrpzIv5FYkp08FoVmLU89IY4ZCXOI0IWbI67R5KaISQDqotzdK/UpAFFNbrN\n5BSNFymM3efYpYh3kGz2ZpG1eISGQXZGzUeK6w4fTFtuCz3LK7F09qSonI1e8NHsH910gYVZbPpr\nIyzhcQ/hvki/pcmOp66LWORkw6PFitcs1RkWvO/kxiCrYJMlSx0D/wCbY+1TugDhduo6rQmpo52F\nkrA9p6hUhhskVxBUvY3k094BXOej3FKSAt5FVntN9VpOp61vzopPMEKvJDVfOp2exynh/s1Up7tq\nWea2GnksswzBwORrba7q0+TKLhZvOCSZ+V5HVVZZLgpKl5JBv4KuXErKVDM7e6y8UGamcbbEFaUg\nKqzx8SJzTzCDnRutzCJP3L2cwbrEIyuIO4WlhMgbUBp2eLKo6qmOamb4aKnVC0xVigJ70Z802sju\nLjko3FCZuYZuiSMKZnQphbw3W5KNynhSNTARZPCjRzrkWTQHgbJzTdSBWM0xshG6mbJfxSta0jUJ\nOADsbIylaWk7kKRrnW0II3sqwikbsbhFy3dpCqJ5XDIbtsVkyzcLMQdSr0sx4Tm3vdYtWf3oaUKn\noIo3zNNQXZCdm7ldI51Ex0YjpHRMy75iCVhYNTvmmDyCGuvlJ523stmarEDjCQDl0uSuk52MXr2k\n4Uc0T3xEysb67T6zVQbUikmMLjnhfqLqP080eJse3uxzNLSqmJsdCXtk0IsRfxWbMV0FNGwzhzOi\n0SBHAT1Wbg3fp4n8+E0nz/4FcqpLRG+wUaU8MHErZHEuudBYbrWqcSFNLwg0OI6lZmBsObiHNYm5\n6eH4qpWS8Wre87Erna1HRRYk19rxn2EFJLi8Dc3ddcciNVi0sJklia5ps86W6LRmkYwvLXHhNGYn\ny5KeVMmoq7FWVFMI23a69z5J1LXPmpmsYdNieqygTO4uLiXu1seSv4bIyF3CkG50IVvf9cn2ePt5\n5dbGG4I6qpTV1U7KWm2D37nyWO3RdNi9PNVYFhklK10kEcVnhmtnWAufcV6HBiz0QdiJpqB5qhpl\nc0etopcQwaswuGOaoyNDzlGV1zdbHZFsMVLXVL82ZgAJYLua219P+clg4rJDJWvNNNLLDoQ6U3JN\ntUVDHPLGQWSOBGgsdrq5FWYhI0lhc8EZDZgN9t/cFnhdDhGCy1UUVS+WOKGMXvINyfd4c1BRnxOo\nLZWSxCN0jszj3mn5vU/yhOkxOObJxoL/AL0PdqDoCTYewgb7ALpMIwSKkleH1Laru3jYHEDx0v5L\nO7SYe5kjpBSMhhyd1zbesNTfzH3K4jNmnwyZriIS06AXbbTrod/O6R1Lh0pbwagtLjpd213ADfwJ\n9yyiluoJpYm+luhhdxBmytPVOlwioO8Ad/S4FR0z3xzCWOMyFguW+Cttx6K9pIZGnw1UoznUnBfl\nfG5jt7EWTZDZWKusbVVWeO+QCwuqkh3U/CoidUpNwU3mh57pKiiMWjHjqpo/VJUY0ACnjbdgAGrj\nYIPQezNPwcCp9NXgvPtP5LVsm00Ip6eKEfMaG+4KU2WpG/ozKkAT7aabpcoT6S+0ZCaWqU32SZVc\nVCWaKJzLq0QmELNgz5YlQnaWuc3RbTm3WbiMdu8OahWe4ZoyOdrqAbaq1CHZ7W0B+xIaeziDfQoy\nqPaoTGtJ0OmoUToVFcpiVOYKkkDuv1ChgeWvBBsQbhdHidCailOUd9mrVzQ7p8QqjsKCQSxxTt9v\nnzV+aMOFxsdlzOB1wil4UhtG87/RK7GFoezI72HoUxYwpYzG+yTK14s5adTTHVrhYrPcwxus5RrT\nBA4Huu0UjYb7uTgVIExdqPLYWCAFIQCmEaoyc11ipGuUCcHWQWmuT+6d9lVa/wAUSy2bZBBVObmN\ntAsR8wlqXEHRp+xWMVqxGwxtPeO6xoZHccFupOlkZdZhNQ50NNSCwdHK63Uhyt12E1XpTnBwId10\nWTSu9GkY5zM4abjqFoTY8ZKbI6V97+0Lpz3+Kzef0oV0bppoYGjvNKbjlSairJvm2F+tgh2IWjtC\n3vEWLjuoWxj1nm8h5dFnrra1I6nA7toG35C3uUmIE8MNbqXmwCbQt4dNFHtZuqTP6RXZcuZrB7ui\nz9Np4g2kw2WUg6jK26z4a3itEbo2PHiPxXQPo4JGMhkuWjUt6lVa7CaOKkkfDDle0XuCsYptMWRD\njajLo25vfyUc8bnRnhs4jHu25Af+UlNTS1FDG2MtdG3vAHdW4qWqY0NDGjrcrGUUyyx4b4WB7tA5\nvzeStUdC3jWbfufOTXtdHM9kmYZm3NtQVp0LckIsBY6pJ7XXkSt0mJ1lHG+Onnc1j9C3cKotgDCp\noLeo5jRrfKX2aNtLakn3L1vOqYZic+GVHFgsbizmnZwUeI1ba2rdO2BkAdbuM2UdLJHFVRSSszxt\ncC5vUK5V1FJUUr3Oc+SquAx7mZLjnoDZBnArucAxGllw2nidMyOSA6h3Pl+K4iJrXXzG1tkr2ZAC\nHA3WpLmmu9rMVpI6pzJIaaQAizpHAcvIrGxvHG1eHmCKAwxk2BJuXa628PzXNtqJmjuyPaPBxCa6\nRz3Xe4uPUm6zoChIi6CxSVclI5zow05hYhw3HRWDiEUlzUUMLybat7pus8HwRdQTVkscs+eJhjZY\nANJvZUpDoVM891V3rNU0JH3uB1KcAkvmlA6KKfu5bfZ6gdXYrTsI/dxfvHnwBWKwXcF6R2aw30HD\nxI9tpZ7Od4DkP+dUWNgi5SFKghdS0jQkG6X5yANVDSHVIU8hNtrZSXWr6NTXDwUhFk0lTDURHXRV\nquPPCQBqNVZee8VE43Cgy2ss4E781K+O7Q7mNCpTDaW/IqzHECCDz0V/CM/hXTXQq7Pw6aJ8kpys\nYNSsOTEamodeD9xGNiRdx/JS4LZh0XN4/hRgf6TEO471h0K6KmjxKQZo6x5PR0wH2Ep80k8YMOJ0\n4ex2mcNAP2aFQcCwlh8F1uAYwwhlLUuPRkhP2FZOM4OaO1RA4SUsh7rh83wKy2uLD4Ij1J7Y5mBk\nuhHquG4WbWUTmeu0FvJ42KwMI7QPpmiKovLF1+c38109HXx1EZdTytkYRqP7JjWsd0WU6JwC13U8\nMpvlyHw2URw5xPdcw/YpgzbKMrUOHS9G/EE39lyu3dG3zcgzEoBJ0C0RQRMbmkkLj9FuyjlDI2G1\nmNCCpYNb4rOrqxtOw3IMhGgTa/F448zIe876S5+WV8zySSSUxKSaV0rySbklaWG0eW0sg15BNw6g\nNxJKLdAVqlummiVJDHFV3tDt1ZeNFC4WUaRt0VvD4jNVsFu63UqrbVbuDUhjizu0L9fYqNBz+DC5\n53OydhEQANRJbXvadOiqyk1dSGMvw2GxNtCeit1kwhbwW211JCzVi6yoLn5irbZQ9uVwuDobrAjn\nLTursVVstaamFHPROMlGc8e/DPJXqKtbUXDmmORu7Sq8dVpa6JRHN3wcko2cFm/6VoBrTmuAS7dO\naA1oAFgOSpUtaQ4RVAyv5O5OVy4KssR5NiEc5kbUywCJtQM7LHQhOw51IG1LauwJhdwibmz+Wyl4\npFLDBiD7U8V3MbbXXlfdZNTWte60MYY0bWXTXJLyV7FKqOuqI3wtLGsiayzrX08isMzPPNIJnjmp\npjSa0Bj8w73zQVHYk2VVlU8aHZW4ahpcHCwIW51uSmJuI0Wa9trC1rKOXIQCwW6phBLjfdStLYwW\nvbr4rrvl6ZzEISqweC62wUDwGusDcLHXOfldIhCFgMkPJQuOqkkOqi5rKl2RGblx9iL2SR+p5lRV\nzDoTUVsMQF872t95XrBAAA6LznsjBxscgJ2ju8+wfmvRTqdFrlfwXkhIN0hK2ycgd1MJ8UF1xYlZ\n6ah5f0TSSU241SB19lZhdPNuabuglIPci4a5mqiLOqmPimkhSxUbo7t0UkH2pb3To7NddSJ057tR\nUB9VHS3tHG0PkA5k7BV6Wp9HcIo4IJJd++M5HTwCkx6JzsZky2zPLA2/Lu2H2plJGKM1Mc7g2SNx\nBvzKSbWbcXXVET5BFiNIyIu2lhblt+aJ3SYVKIKj/wCoophca3Fuo6FZ0lS6WNzZHl7SNj+CqCol\nlwtrXPc8MJ0PJXqJzWhRSxtqZKZ1paaXTK4aOCzMc7OPpS6ejBfFuWblv5p1HLmliI3aQF1hcXCy\n5ytvL75TzCmgrJIXh7HFrh85psV21f2dpK67svCk+k3n7Fz1b2Wrackw2nb4aFaTE1J2nnjs2XJK\nP5hY+9akPaemc274nt8iCuSmw6shP7ymlb/tKr2c02II81DXeDtJh/05PgUUvaahbewkd7LLhru6\nlWaejfMMz3ZWnbqUq63artQXaQwho6uN1jVmKVFVpJISOgWhFhkBb/Ac7xJN0+HCaT0lvGa9rb3I\nupplYsFJUVbv3TC4fS5LbocGZTgPl7z/ALl1cVFAIAYA0N8FSmhyuOiGM50Vhpoo3MV17NFEWIKZ\namOZvdXOGk4Be4Ac0EFFSGeYC3dG62amXhtEEX8R2nkoiRRRBkbc0ruit0NJkJln7zjY3Klaiejg\n9GiFyeTrHnpumSxCRxcRqrGrt04MTBnupw3ZK0ZSr/CukMF+SYIo7dVaYo2RBu4U8QaNwUDjCyRu\nV4uPuSNqH0bxHU3dCdGy9PAqzG0HZSmNr2Fj25mncFMV5RjcstayGoBL2Nblf4OudT5iyyREXDZX\n6SpMbrHVpFiDsQrLqSw49IC5o1czm3y6hbcWIYyEhaei1QyGqGVv7ubkDs78lVex0bi14IcNwUwU\n0NJabhWS0HcJj4huFFWKae9iNHDmnuJcbu3VFpLHgq/FLYh4aDpbVb5v4t9JT44846IdC4AuuLBM\nzkE2NvJLxX/SK6bzntn2QtIFyEl0Zid0hNmlc60gdzTW7pXIbtdYDXnulPAsAOgTH94gdSn80Vr4\nHWPoJ+PHvaxB5hd5h+JQ1sQdGRmtq08l5dJM6KRltgFp4biEjZozE4h1wLBFlekZxfRAfc2VBlWC\nAU/0gLW1ci25wTMwVbjA7lMM6aq5nCUPAGioCcpeMpp9LvEATg/RUBIU8SdU0Wc2uqQlQh5JTwbo\nJQU7N0UVkahJ6Ptldo6YvMdU3XKMj/LkVn01LHPGTxmmbd8cz8pd4tO3vXSyAPaWuAIIsQVk1GHv\nBvFEyoYPmvNnDyPNWXGbGPVU0cDWiCRz5nf5RsSzxNioX3ipRE3XNoPLqtKRk7RljoBH1JY5x/JJ\nS4PV1Dy97C3+Z/5LHXW+o1OUGEURMzegOZxXTRtUEFK2mblYNeZPNW4gnMs+xIGaJCxSBODbrSq5\njuFA+igef3kTHebVfLEhYp9jmsZpMPo6Z0gpYzI45WXb9qzKGnjJa+eQsadrC5K0u1DSaqnZ80NL\nln0dOyorAWkljRfVZz2z9N6nZQm0cFRMwnQFwFvsSVsD6YgVjWyxO9WVvL2qs6tpadxjip2Ajd3O\n6G4hLUwvpDo06i/PTRapKjgq/RZ+G15dETotOZgcMw5rmSHB+Uk5muXSQPz0zL9FnRVcxRcK7lcc\nzdDIi435DmtIpiHVI9/C7sYzSnTT5qmle55LIenrBPhijhH0nrNrRKOk0Ek+rtdVoFxeRflooGFz\nlMxpSRUjQpWhIxuo0U7GG60HMbongDmE5oslyqqThNI2TTAL91TBpsE4BBXa1zFOx1wnWBTQ3KdF\nB4jG9aVBVuheLGyyBdpU8b1XFt1dCyrbxqUAS2u5oPreI8VUZI2rAiqDllGjJDz8CnUVTlIBd5FW\n66kFW0zwttM31mj53j5qoyJYnxPLXixCYr0TxVsEMpAlAsx55+BVSRjo3ljgQ4aEFRVd7dVLAe7b\nokcEQ6EoL88MEccbo6kSOdYuaGkZVZIw6YiziywAGuXmN977n3LMGyLEb81WbCnfRNf6qVNkOwRp\nC5OGjU06lPI0CyqMayeACkbq4JjPnH2KemjMkzWjcmyB3BMr7AXW7g+ECJ7ZpRruArtHQMiN8q02\nNRZCZNEpFraqQC6QtRpESUlynuGlky2miIUHVP5qMKTxQPCXdNGyAdUVMw8lM0iyrNJupgdEE7NQ\nnWumMdZLmQIRpdRO0UznaKJ1rIFjma13f26rSgawsDmWI8FivCbHUS07s0btOnIqzIjdkgZINRY9\nVXMBjOm3VQU+Mwus2YcN3XktBr2ytux7XN8Fv1U9q4FzopWs5p5jaeVilyWHVS86aba26S11JtyT\nSU8DXP8AaSmN4ahg1Ycp9qyTK88PhtY2SMBji0AZr89F2E8bJoXMkbcOXL1tC6lmzHPlabte3cea\nz9XRnTUszKg3c03O4Kko6rjYiwjSOPS/kisxCaeDhS1odHzGXU/Ysx07cnDpwRmFi78lnrqLIsOn\n41XIY/nONl01K1wjYzU2ACw8MoOERNMMrRsDuVrHERGMsTdepWeb72rV58TI23ld/tG6pzTGQZb5\nWcmt/FVnSyynM5x1TmRkkK9XUSN9XKwWB6KWOMmydHGrUcd7JIGxxqyxit09MzIC4XupxBGNmrTW\nKjGKZjeanEbRyRlA2CoYBfknNFkqVXE8oEIQpi6EqELKvDXtTGmx1Uzgo3tVcUsb7LUoqstsCTfk\nbrFa7krEchBVGtiFKJYzVwaEayN//oKC/p1Pr/iIxr/O381NRVeVwubnoVFVw+h1bJoTaN+rdb26\nhVFGybGO8VflhaSXAaO6dVUa0GQNJygm1+iip3Vcr6VlO4jhs1At/wA6q83FIXsyywWFg0huxH/L\n+9VvQQ6r4MEzJW2zF4GgCSTDqiNmbKHeDTc72t71XO+NMq307yw08eQ/O8f+aqlJ6ynLHMdZ7S09\nCLKu895RuTDRunu0TWesiT1So0RvqDxWrgMPFxGPT1bu9yy7WAC6PstD+8llPJoHvQdC1ilAKdG1\nStatKhDT7E6ymDUGO+wUrSq5t0wjkrZiPRMdEeiIrZbbJWnqFMY0mQXUU21hdIBqpMuiaRYoFCkY\ndQowlYS0kHVBYB0Sh+iYE4aiyBx1CaQQE4IKCBw6qJ7dD0Vgi6rvBCVKpStsSq4nlgdeKRzD4FXZ\nNd1TkZqVBcg7Q1UekgbIPEWKvw9pIHfxInNPMg3XOObqmlqvlR17MboX/wCaR5gqUYnRkXE7fauK\nKBor5mO1dilH86dvsBVebGKHLl1f7FyepUjdlPMxcrDh9Q4uFI0k8yqsbI4zeKJjPIIAupWtXNSW\nc/dSsj8FJGzu9VPGzwVDI4bq1FClYyymaCLEIhzGCylaLJrRpa6laLrUWLdOe5qp1Xg0U11c1TkJ\nGne6VLMSXTQlS6IK15e2fEiEWKUBW2E0IQi65tvG6qmdDIWkKo9q3IpY6uPhy2D+TiqdTh8sRJy3\nbyIWnFkuGqcxylfFysoSwhRW12eoP2ricdLxuECCS61zYdFqdocHlwyLgvfxYz34pALbbgjyK5am\nqZaSdk0LyyRhu1w3BXT4j2vjxXBvRqqny1TSC17fVPI+WhKujHbMBSFumYmwVcNOjrHKdj1SE5ms\nt9FW4ql7o4qaSQCBr82rbgfikLVdj3Mddji02tcKyMRqdQXh1+rQrbo8PlF87GXNyQSDpvp42HLn\n4KlPTsE0rIpAcrjYE8vNVz2X8GVVS6oeXvte1tPO/wCKouU79GqByjpDo+ZSPF8vilZo3zSP/ieQ\nUUo1K7Ls5Hkw+/03krj4m5nhd5hsXBoYW88t/eqsXoxqp2NUTNLFTNvbwSatO0unBInbDRbs1kwi\nxSFvVSgX1TsizeWog4aYYxdWg1NdHqpiqZYW6phbfVWy3koizkoICCAlaNNVI4aJmVA7dPb4Jo0U\njQgdpbxTTcpwGiXRURHRQv3VlwCheLhQVyy99FA+HoraYfFBnSQkX0ULmWWoWAqGSEclkZrmpOG5\nXjAk4RQVGsJuE4R23Vnh2S8O4UEcbL+Sma3VKxllMxlygWFpcbNC0Iaa/rXso4YTmGVwBbuArcbc\nw11Xy/5H8u748OvPG+6ifGG6a6DWyRrbWII15KZ9wNvO6h7kZMjm5jfey4fH/L+Tm/bV+OJmttrZ\nSN3TGyNsCGPcDskFQ1psGONj0X1uf5Hx38uPjVqE95WsthqqFNWROflsW26rQuC0arpz8k6v9S6b\nYJUc0LpaYEqS/JCgVCRB2SIE0jVKEqv+NPuPGmOa51wcvgT+KvwVU0bbWzM94WfkjO92nw1CVjHt\nuY5QPAGxVcmhJVQPHep2EqjOWu9WEMCHSVZFjM8gf/k/uo3McTdzwT5oqAxNJufsS8Jh2ClDQNd1\no1oqsRk9IbR5Wtjbm4TNLAbm3kmGoo5IpaKKjbTxMla8njk2JHQqOpop6VjXShovpYOBINr2I5aF\nV1YrKx9WWukZGHj1nNbYuPUrSK6EiVAyTQBQu3Usp72ii5rKpWDRo57ply4uJ8lKDk7x5BRD1fPV\nRVrD4zNUsYOZAXfsaG7bLjuzsRfXxkfNu4+xdmzzVipWG5UrVC299FM1az0aeByCcAUNITgE9oUJ\nyAAkVWlCLXQlHRCIyzVMLQpybphbYeKzdX0rvaLX3ScMKUg3sgjVTFRZE8BOIJSWPiof9AdLX809\nrRfokAN/xUmQkaK5U+0Tm6qu9vgrhabbKF8bhyTBTcCE07KdzT0Ubme5FRFJZPc1M1CgMqTJdPGo\nT8uqiIeElEWmysNaOafkUVWEYUjWd4KXhpzWa8wOadeoIKqWDDqWSrl7ot0Gp5Kxh9ZFWUQqGABp\n36C3K65urpn4h6bPXukjijBfEz5xYNBpyTMKrZq7CZ8PpGMjeG903sbdPFfH7/jS8W7711nbqfSo\nZD+7lY7rZwTIXiafV3dXJ4RJJh5mNXLw3SAWjzfaQtulxEMiqKl0jXRxizBcAnquPf8AGvF/r7jU\n79N4WAs3l0TXNG7lmYPizcSa54aWtbvqLK66ojkLmRtEz28gRouN+PvnrLPbUp3AaAdCR0VqB7xY\nOGnJZLsTbHUNgqmOicRe4uQFPTU9eyt4ktQ10HJvUL1fDevj68uvTPVlnps77JdeigMuXZROrcrd\nNTey9v8AzuI5+FXDe2yFBT1BmLm7WNlYI0Xp+L5p8nOxnr+v2RCACUe1dU2BCEIrxhBSIO604luh\nIhUKrLK2RlE6laG5Xuu4kAna1tdue3VVUKAQhCoEBCTkUETtTdNbq9Kdikj3usqmcQ2GS+5Oiits\nFJK0cOK+51TGi7gPFFdL2Xis6SToAPf/AOF0zdNbLJ7PQhlBm+m4n8FsNHVFPaRa6kCYLW0Txf3r\ncRIBqpQNNd0yM25KQi/NMqyh2mgStGmpRa9glAVSkA5p4aDulATwLBVDAwJC3kpLBOski1Bw7lBi\nvorFrIshqvwSnCOwU5SAaJgjEYvdPDQPalASoaSw6JCxp5JwueSENRGBh5JhpYzyVhCGqjqCM7KF\n+G29U3WkEqlka8mK6ikZyuEghI3C20jmtPILPimscRkpwabbLSNOwpjqfpqpjSiGlUcTrZKM0zIW\n3knlDNr2HMraMJsuZ7SYlFS1EYis6pgdlv8AQzDWwXP5fU+j8s/BxLJi1Qys78k92ASGxIP9gEk2\nNR4ZVupKGna2IEAjLq/2qR9JUftd9fNMeGzKGvvblofK6ZjstGY4qqma0VGziw39q8H9e+5LNmf+\nN5kUcWhqf2jx5wZKd9nh56HYXVjDy6vpquGQlsB1DWt2t0uq/Hq67D2yyyk8F9mMcbZr/wDhX3uk\nHCeIWZ5GkBlrAaWLit9bnj+U/wBpcOp5KaF1PG1ga+xdICAQDr9y1cNhqHxcKNjYoNRnFi53issM\nZVubS0LGPMQu5+vfXU0kPo1MyMANIAuB1Xn+b5v/AJz191rmaSCihgZktxLc36qxoAo3G976AKtU\n10dKQJNCfVvoD7V4vK/LcdcxNLKQTawbzcVnYtiQw+njla1pzHS6oYhVvxCnfFCC0g31da1uSpVV\nXwsBgp61rXyvcSM2trHRenj+Neep5MXv8R0FJjdM6ip6mQ5OMbW6Hmt9hD2hwdcEaLgZHUtZgETn\nObTuY8hoaN/7LpcIM/EBEzTSCFoYDudNSvV8HM469Mde2wSboFkgN9krTqvpes9MYcEzO29hulvY\n35LDr8S9HnLQbe1eb5+7xNi8+3mnNCqekv6NR6S/o1enXLGjSUz6upZBHYFx3OwHMq+/Ap2zGMSR\n3zWF7j5xaOXMgrn/AEl/RqsNxiuaGgTus12YAkkA9U1MXpsOqYmtdkD2uIDSxwN77ab6206qGWnm\nhLhJE9hbvmbayYzHatnAGWE8FzXDuesW+rfrZTR9psQYWk8J+UWsWmx1B2BA+aB5IZVdCllx+aVp\nBpKMXsNIugtzPl7lVmxAyyZxTwRi3qsDrfaU09pCkcbNVb0l/RqHVD3CxDU1cPdsnR6C/M6KuZXH\nkE5tQ5trBuhuoqzM3LPlvewRA28zfNVnVD3PLiBcp0VU+J+YNaT4oPRsPi4dDAzowH8VcYAuGb2w\nxBosIaXT+V36k4ds8RH+RSfA79SK71jU7UFcEO2uJA34NJ8Dv1Jf+98S/wBCk+B36kn+zXfs0Uzd\nQvOh23xMG/ApPgd+pPHbvFB/kUfwO/Ut+UR6JbTRDQea88/79xS1vR6P4HfqSjt9igP+Hovgd+pP\nIeitTxtqbLzj5QMV+r0XwO/Uj5QMVv8A4ei+B36lfKI9I0Tl5t8oOLfV6L4H/qR8oOLfV6L4H/qT\nyg9JQvNvlBxb6vRfA/8AUj5QcW+r0XwP/UnlB6SdUcl5t8oWLfV6L4H/AKkfKFi31ei+B/6k8or0\npC81+ULFvq9F8D/1I+ULFvq9F8D/ANSeUR6UDbZIvNvlCxb6vRfA/wDUj5QcW+r0XwP/AFJ5Qekp\nQvNflBxb6vRfA/8AUj5QcW+r0XwP/UnlB6TdKvNflBxb6vRfA/8AUj5QsW+r0XwP/UnlFelIXmvy\nhYt9Xovgf+pL8oWLfV6L4H/qU8oj0lLfS1l5k7/qBi7iCIqRtuQY7X/3Jr+32MOcCGUrLcmsOvvK\neUV6bnjEgY57cx1Db6rz/EmwDEMRqOHne15cCTtfb7Vz9R2gramqdUvyCU/OGbTy1TarHKuqka+R\nsVwLGzSM3nquPyb19NSteaulxLCOESeI12h87aKbDW8Cl4U7A6R79C75o2WDS41U0smaJkVvo2Nj\n9qfPj9VUVPHfFAHXvZoIA8N1wvxde+Z9NeTbqDIcSIaHPggOZoAuBporJeKl4icZJKp+kjoyLAHk\nuem7R1cty2GniJblJjaRce9Oou0tTQvD4qWkLwLZnNdf/wCSl+Lrx+vZ5R6JhWGw4dARHcl2t3bh\nXSeq89Pb3FD/AOno/gd+pQz9tcTmjLeHTMJ+cxrgR/7l8/r+F8/fW9Ok+TmOwxDEWcaSizlkzhZj\nuV1Rjpp/QagYlmLIgCw3ub+C5R3aaskY0Sw08rmuzB72uzfepHdrcRdPxXMgJtbKWut7rr0c/wAX\nviZyzfkjpHTMmhDobFtTdoHzgbDdZ9bMaTDm0dVCyYtcSx79wN7LBjx6piDgyKABxJIseftRU47U\n1NuJFBoAPVPL2r0z4uvUrPlHRV0dHWYVQinlbBLb1bb+fRauGQ1NFiVPSAl8bafvu5X8FwlVi89S\n5p4UMQa0DLG0gac91pQds8SgFmxUzvFzXE//ACTn4u+bP0eUeoxOuLp99F5mO32Kgf4ej+B36k4f\n9QcWH/p6L4H/AKl6ucjEr0x3eYuQxZxbiTcoDja1zyWH8oeLWt6PRfA/9SzKvtPXVcvEfHA0/wAr\nT+a598+Uwl9sVCELogQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQh\nCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQh\nCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQh\nCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQh\nCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCAQhCD//2Q==\n",
      "text/html": [
       "\n",
       "        <iframe\n",
       "            width=\"560\"\n",
       "            height=\"315\"\n",
       "            src=\"https://www.youtube.com/embed/-gCTKteN5Y4?rel=0\"\n",
       "            frameborder=\"0\"\n",
       "            allowfullscreen\n",
       "        ></iframe>\n",
       "        "
      ],
      "text/plain": [
       "<IPython.lib.display.YouTubeVideo at 0x7f495bc2eac8>"
      ]
     },
     "execution_count": 2,
     "metadata": {
      "tags": []
     },
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from IPython.display import YouTubeVideo\n",
    "YouTubeVideo(\"-gCTKteN5Y4\",560,315,rel=0)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "ffvdVsm5wigL"
   },
   "source": [
    "### Saturation Vapor Pressure\n",
    "\n",
    "At a given temperature, the **saturation vapor pressure** is the pressure of a vapor when in equilibrium with the liquid phase. \n",
    "\n",
    "By itself, the term 'vapor pressure' may refer to a component of a gaseous mixture that is not in equilibrium with a liquid. Adding the extra adjective 'saturation' makes it clear that we're referring to a vapor-liquid equilibrium."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "GoJ_G3QfxjLB"
   },
   "source": [
    "## Clausius-Clapeyron Equation\n",
    "\n",
    "The Clapeyron equation is a relationship between the saturation vapor pressure of a pure substance and absolute temperature $T$.\n",
    "\n",
    "\\begin{equation}\n",
    "\\frac{dP^{sat}}{dT} = \\frac{\\Delta\\hat{H}_v}{T(\\hat{V}_g - \\hat{V}_l)}\n",
    "\\end{equation}\n",
    "\n",
    "where $\\Delta\\hat{H}_v$ is the latent heat of vaporization. This is the amount of heat needed to vaporize one mole of the substance in a constant pressure environment.\n",
    "\n",
    "This equation is derived by comparing the change in free energy of the liquid phase compared to the free energy of the vapor phase. The changes must be equal for the two phases to remain in equilibrium which leads to this equation.\n",
    "\n",
    "**Suggested Study Exercise:** *Compare the units on the left and right hand sides of this equation. Convince yourself they are the same.*\n",
    "\n",
    "Under typical process conditions, the molar volume of vapor is much larger than the molar volume of liquid, i.e., \n",
    "\n",
    "\\begin{equation}\n",
    "\\hat{V}_g \\gg \\hat{V}_l\n",
    "\\end{equation}\n",
    "\n",
    "In turn, the molar volume of the saturated vapor can be approximated by the ideal gas law\n",
    "\n",
    "\\begin{equation}\n",
    "\\hat{V}_g \\approx \\frac{RT}{P^{sat}}\n",
    "\\end{equation}\n",
    "\n",
    "Putting these together,\n",
    "\n",
    "\\begin{equation}\n",
    "\\frac{dP^{sat}}{dT} = \\frac{P^{sat}\\Delta\\hat{H}_v}{RT^2}\n",
    "\\end{equation}\n",
    "\n",
    "Integrating,\n",
    "\n",
    "\\begin{equation}\n",
    "\\int{\\frac{1}{P^{sat}}dP^{sat}}\n",
    " = \\int{\\frac{\\Delta\\hat{H}_v}{RT^2}dT}\n",
    "\\end{equation}\n",
    "\n",
    "Assuming $\\Delta\\hat{H}_v$ is independent of $T$ (which is only an approximation) gives the **Clausius-Clapeyron** equation.\n",
    "\n",
    "Doing an indefinite integration gives\n",
    "\n",
    "\\begin{equation}\n",
    "\\ln{P^{sat}} = K -\\frac{\\Delta\\hat{H}_v}{RT} \n",
    "\\end{equation}\n",
    "\n",
    "where $K$ is a constant of integration.\n",
    "\n",
    "When integrating between specific bounds on $P$ and $T$,\n",
    "\n",
    "\\begin{equation}\n",
    "\\int_{P_1}^{P_2}\\frac{1}{P^{sat}}\\,dP^{sat}\n",
    " = \\int_{T_1}^{T_2} \\frac{\\Delta\\hat{H}_v}{RT^2}\\,dT\n",
    "\\end{equation}\n",
    "\n",
    "one obtains\n",
    "\n",
    "\\begin{equation}\n",
    "\\ln{\\frac{P_2^{sat}}{P_1^{sat}}} = -\\frac{\\Delta\\hat{H}_v}{RT}\\left(\\frac{1}{T_2}-\\frac{1}{T_1}\\right)\n",
    "\\end{equation}\n",
    "\n",
    "This equation provides a useful means of estimating heats of vaporization from measurements of vapor pressure.\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "dssZoqitp4qh"
   },
   "source": [
    "## Antoine's Equation"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "a0GqDFk0G2C2"
   },
   "source": [
    "### Vapor Pressure"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "iZJDCsBKp4qi",
    "slideshow": {
     "slide_type": "-"
    }
   },
   "source": [
    "Antoine's equation is used to estimate the saturation pressure (also called vapor pressure) of pure substances between the triple and critical points. Louis Charles Antoine, an engineer working in the French Navy, in 1886 published the equation as method for representing the vapor pressure of water. \n",
    "\n",
    "A common form of the equation is\n",
    "\n",
    "$$\\log_{10}P^{sat} [mmHg]= A - \\frac{B}{T[^{\\circ}C] + C}$$\n",
    "\n",
    "where pressure is units of millimeters of mercury (mmHg, also called torr), and temperature in degrees Celcius.\n",
    "\n",
    "Compared to the Clausius-Clapeyron equation, this is an empirical expression with three parameters that can be adjusted to fit data for a variety of compounds. The algebraic structure is similar to the Clausius-Clapeyron equation, which gives considerable confidence that Antoine's equation should work well over reasonably wide temperature regimes.\n",
    "\n",
    "** Example:** In the following cell, create a Python function to compute the vapor pressure of water in the range of 60$^\\circ$C and higher. Verify by computing the vapor pressure of water at 100$^\\circ$C."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/",
     "height": 34
    },
    "colab_type": "code",
    "executionInfo": {
     "elapsed": 507,
     "status": "ok",
     "timestamp": 1539263118367,
     "user": {
      "displayName": "Jeffrey Kantor",
      "photoUrl": "",
      "userId": "09038942003589296665"
     },
     "user_tz": 240
    },
    "id": "wHgGVSt-7TiZ",
    "outputId": "899650a8-70e4-40a1-ea18-85afc08e238e"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Vapor pressure of water at 100C = 759.98304330841 mmHg\n"
     ]
    }
   ],
   "source": [
    "# water (data from Appendix B.4, Felder et al.)\n",
    "A = 7.96681\n",
    "B = 1668.21\n",
    "C = 228.0\n",
    "\n",
    "def Psat(T):\n",
    "  return 10**(A - B/(T + C))\n",
    "\n",
    "print('Vapor pressure of water at 100C =', Psat(100.0), 'mmHg')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "dPVDqmiGHEsD"
   },
   "source": [
    "### Saturation Temperature"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "TP7X4kCu7MKB"
   },
   "source": [
    "An alternative form of the equation is to calculate the saturation temperature as a function of pressure\n",
    "\n",
    "$$T^{sat}[^{\\circ}C] =\\frac{B}{A-\\log_{10}P[mmHg]} - C$$\n",
    "\n",
    "Values for the constants $A$, $B$, and $C$ are tabulated in various references, including the NIST Chemistry Webbook. The values of the constants depend on the units used for pressure and temperature, and whether the logarithm is computed for base $e$ or base 10. \n",
    "\n",
    "**Example:** On a given day, the air pressure at the Denver airport is 623 mmHg. What is the boiling point of water at that pressure?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/",
     "height": 34
    },
    "colab_type": "code",
    "executionInfo": {
     "elapsed": 401,
     "status": "ok",
     "timestamp": 1539185463542,
     "user": {
      "displayName": "Jeffrey Kantor",
      "photoUrl": "",
      "userId": "09038942003589296665"
     },
     "user_tz": 240
    },
    "id": "5fE6RK1jp4ql",
    "outputId": "c94a2f78-099e-48d9-b1da-908ea39ce358"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Boiling point of water at 623mmHg = 94.52632667664744 deg C\n"
     ]
    }
   ],
   "source": [
    "from math import log10\n",
    "\n",
    "A = 7.96681\n",
    "B = 1668.21\n",
    "C = 228.0\n",
    "\n",
    "def Tsat(P):\n",
    "    return B/(A - log10(P)) - C\n",
    "  \n",
    "print('Boiling point of water at 623mmHg =', Tsat(623.0), 'deg C')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "tbmJZgG3kEE5"
   },
   "source": [
    "### Working with Wide Temperature Ranges\n",
    "\n",
    "Standard practice is to specify a range of temperatures over which a particular set of constants is known to offer an accurate representation. Multiple ranges may be pieced together to obtain saturation pressure over wider ranges.\n",
    "\n",
    "**Example:** Create a python function `Psat` to compute the vapor pressure of water from the triple point to the critical point, and use it to create a corresponding plot."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "colab": {},
    "colab_type": "code",
    "id": "7Y76qMSap4qv"
   },
   "outputs": [],
   "source": [
    "# Antoine's equation for water from 1 to 374 degrees C\n",
    "\n",
    "def Psat(T):\n",
    "  if (1 <= T < 100):\n",
    "    return 10**(8.07131 - 1730.63/(T + 233.426))\n",
    "  elif (100 <= T <= 374):\n",
    "    return 10**(8.14019 - 1810.94/(T + 244.485))\n",
    "  else:\n",
    "    return float('nan')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "l3oBDG0KnZfE"
   },
   "source": [
    "The following code uses the `numpy` and `matplotlib.pyplot` libraries to create the requested plot."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/",
     "height": 407
    },
    "colab_type": "code",
    "executionInfo": {
     "elapsed": 981,
     "status": "ok",
     "timestamp": 1539185833734,
     "user": {
      "displayName": "Jeffrey Kantor",
      "photoUrl": "",
      "userId": "09038942003589296665"
     },
     "user_tz": 240
    },
    "id": "3ZL2UBfbm-Xe",
    "outputId": "e9336bec-a991-412c-90ab-a23eee0f5892"
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmoAAAGGCAYAAADLvK0bAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJzs3Xd4VFX+x/H3TCa9JyQhgSRAgBBC\n7xBAqoCFoqgRUVlZFdHdtaCwqCBFsexaQFkVy09EJQu6igqCgqAivXcIJUASIL33zO+PwOAYIJRM\nGp/X8+wjc+beO985G+DDPfecYzCbzWZEREREpMYxVncBIiIiInJhCmoiIiIiNZSCmoiIiEgNpaAm\nIiIiUkMpqImIiIjUUApqIiIiIjWUgpqIcPfdd/PZZ5+Va4+JieHuu++u8nrmzJlDp06dGDx4MIMH\nD2bQoEFMnTqVvLy8Kq/Flk6cOMHAgQMZNmyYVfv//d//MXbsWKu2yZMnl2ubOHEib7zxxiU/48iR\nI2zatKlyChaRKqegJiLcdtttfPvtt+Xav/nmG2677bZqqAgGDRrEDz/8wA8//MC3335LUlIS77zz\nTrXUYitbtmzBz8+Pb775xqq9R48ebNmyhcLCQqtjY2NjKSoqsrRt2LCBHj16XPIzfvrpJwU1kVpM\nQU1EGDJkCPv37+fEiROWtpMnT7Jv3z6GDBkCwKJFixgyZAg33ngj99xzD/Hx8UDZ3a9Jkybx8MMP\n07dvX6Kjo0lJSQEgISGBsWPHMmjQIG655Ra+/vpry7V79uzJSy+9xOjRoyusz8HBgbvuuou1a9cC\nMGnSJGbNmsWtt97KsmXLKCwsZObMmQwaNIh+/frx7rvvWs5dsGABQ4YMYfDgwYwcOZJDhw5dsj08\nPJxTp05Zzj/3esOGDURHR/OPf/yDp556CigLQbfeeiv9+/fngQceIDU19YL1L1u2jFtuuYXBgwdz\n3333cfz4cbZt28a//vUv9u7dy9ChQ62Ob968OW5ubuzYsQMou/Pm5OREs2bN2LVrFwDHjh0jPT2d\n9u3bU1payrRp0yzf/+mnn6aoqIhVq1bx3nvvMX/+fF5++WWg7C7p4MGD6devH08++ST5+fkX7FMR\nqRkU1EQENzc3BgwYYHVn59tvv6V///64ubmRkpLC9OnT+fjjj1mxYgUhISHMnTvXcuyKFSt47rnn\n+PnnnwkODua9994D4Pnnn6dLly4sX76c9957j5kzZ3Ly5EkA0tPTiYiIYMGCBZdVY1FREQ4ODpbX\n69atY/HixQwZMoR58+YRGxvLt99+y3fffcfy5cv5+eefyc7O5q233mLRokX88MMPjB07ltWrV1+0\nvSJ79+4lOjqaf//735w4cYJnnnmGf//736xcuZKuXbvywgsvlDsnISGB559/nnfeeYcffviBPn36\nMGXKFNq3b8+TTz5Ju3btWLJkSbnzoqKiWL9+PQDr16+nffv2dOzYkQ0bNljaOnfujIODAz/++COb\nN2/mu+++Y9myZezZs4elS5fSr18/Bg4cyH333cekSZPYvHkzb731Fp988gmrVq3Czc2Nt95664J9\nKiI1g4KaiADlhz+XLFliGfb09fVly5Yt1K9fH4BOnTpZ3X3r2rUrwcHBANx4441s27aNoqIifv/9\nd0aNGgVAgwYN6Nq1qyV8FBUVMXDgwMuqLTs7m88//9zq+O7du+Po6AjAzz//zKhRo3BwcMDFxYVh\nw4axYsUKHB0dMRgMLF68mOTkZIYMGcKDDz540faKODk50b17dwB++eUXunTpQvPmzQGIjo5m1apV\nlJSUWJ2zdu1aunbtSmhoKAB33HEHGzZsoLi4+JKf1aNHD9atWweUDXF26dKFzp07WwW1qKgooGyY\n+Msvv8Te3h5HR0dat25t9f/POatWreKmm24iICAAKHs2ccWKFRfsUxGpGUzVXYCI1AzdunWjoKCA\nHTt2YDQaycvLo1u3bgCUlJQwe/ZsSxDJycmhcePGlnO9vLwsv/bw8CAzM5P09HTMZjPu7u5W750b\nHrSzs8PNze2i9SxfvpwtW7YAYG9vz8CBAxkzZozlfU9PT8uvs7KymDVrFq+//joAhYWFtGnTBnt7\ne/7v//6Pd999lzlz5hAeHs7UqVMJDw+/aPul/PkzN2/ezODBgy1tbm5upKen4+vra2lLS0vDw8PD\n8trd3R2z2UxaWtolP6tHjx48++yz5ObmsnHjRv75z3/i7u7Ovn37KCwsZOPGjYwfPx6A1NRUZsyY\nwd69ezEYDCQnJ3P//feXu2ZWVhY//vgjv/32GwBms9nqmbc/fj8RqRkU1EQEAKPRyLBhw/juu++w\ns7Nj2LBhGI1lN92XLl3KqlWrWLBgAT4+Pvz3v/+1uvv2x9CRkZGBp6cn3t7eGI1Gy2ugXIi5lEGD\nBvHiiy9e1rH+/v488MAD9O3bt9x7LVu2ZPbs2RQWFvLBBx8wdepUFi5ceNF2o9FouSuWkZFxyc/s\n0aMHs2fPvmRtvr6+bNu2zfI6IyMDo9GIt7f3Jc/z8/OjSZMmfP/997i7u1v6LSIigiVLlmA0Gi13\n89544w1MJhPffvstDg4OlmfoLlTziBEjmDhx4iU/W0RqDg19iojFbbfdxqpVq1i5cqXVbM+UlBQa\nNGiAj48PaWlpLFu2jJycHMv7W7ZsITExESi7E9axY0dMJhM9e/YkJiYGgOPHj7N58+YKZylejf79\n+7No0SJKSkowm83MnTuXX375hQMHDvD3v/+dwsJCHBwcaNWqFQaD4aLtUBaQ9u/fD8CXX35pCat/\n1rNnTzZv3mwZYty5cyczZ84sd1xUVJTVcQsXLiQqKgqTqeJ/J0dFRTF//ny6du1qaevSpQuffvqp\nVT+mpKTQvHlzHBwc2L9/P9u2bSM3NxcAk8lEVlYWAP369WPFihWWu5o//fQT77//foV1iEj10R01\nEbEIDQ3F39/f8utzbrnlFr7//nsGDhxIcHAwjz/+OI888ggvv/wyrq6u9OjRg2nTprFv3z6CgoJ4\n9tlnAZg2bRrPPfccX331Ffb29sycOZPAwEDLhILKMmrUKE6ePMnNN9+M2WymVatW3H///bi4uNCw\nYUNuueUW7O3tcXV1ZcqUKTRv3vyC7QBPPPEEL7zwArNnzyY6Ovqiw7P+/v7MmDGDRx99lKKiIlxd\nXZk8eXK54+rXr8/MmTMZP348RUVFNGzYkBkzZlzW9+rRowcfffQRjz76qKWtS5cuvPXWW/zlL3+x\ntD3wwANMnDiRr776ik6dOjFx4kSeffZZ2rRpQ9++fZkwYQLx8fHMnj2bcePGce+991JaWoqvry/T\npk27kq4WkSpmMJvN5uouQkRqrzlz5nDq1KnLHqYUEZHLp6FPERERkRqqRgx9JiUl8d5771FcXEx0\ndDQtWrSo7pJEREREqp1N76gdPHiQAQMGWC1o+dJLL3HXXXcRHR3Nzp07AVi8eDENGjTAyckJPz8/\nW5YkIpXsb3/7m4Y9RURsxGZBLTc3lxkzZlgWhwTYuHEjcXFxxMTE8OKLL1r+cE9MTOTGG2/krrvu\nYv78+bYqSURERKRWsVlQc3BwYN68eZYZZFC2PcmAAQMACAsLIyMjg+zsbHx9fTGbzbi4uJCXl2er\nkkRERERqFZs9o2YymcqtE5ScnExkZKTltY+PD0lJSYwcOZLZs2dTUlLCww8/XOG1i4tLMJnsKr1m\nERERkZqkWicTnFsZpEGDBrzyyiuXfV5aWu5VfZ6fnztJSVlXdW5dpP6wpv6wpv6wpv6wpv6wpv6w\npv6w5ufnXvFBF1Gly3P4+/uTnJxseX3mzBlNHhARERG5iCoNalFRUSxfvhyAPXv24O/vf8lNmUVE\nRESuZzYb+ty9ezevvPIK8fHxmEwmli9fzpw5c4iMjCQ6OhqDwcDUqVNt9fEiIiIitZ7NglqrVq34\n9NNPy7VPmDDBVh8pIiIiUqdoCykRERGRGkpBTURERKSGUlATERERqaEU1ERERERqKAU1ERERkRpK\nQU1ERESkhlJQExEREamhFNREREREaigFNREREZEaSkFNREREpIZSUBMRERGpoRTURERERGooBTUR\nERGRGkpBTURERKSGUlATERERqaEU1ERERERqKAU1ERERkRpKQU1ERESkhlJQExEREamhFNRERERE\naigFNREREZEaSkFNREREpIZSUBMRERGpoUzVXYCIiEhdYDabiYn5jO+/X0JxcTElJSV06dKdceMe\nw83NrdzxX34ZQ2pqKg8++Ah79uzG0dGRpk2bWbVfqcTEBKKjR7BmzYZy740ceStmsxlHR0dKS804\nOjpw771/YcCAQZe85uXWs3LlCrp164Gra/nvKldPQU1ERKQS/Oc/c9i2bQuvv/42fn7+5OXl8dZb\n/+KZZx7nnXfmYTAYrI6//fa7LL9eunQJbdq0o2nTZlbtlW3KlJm0bdsOgOPH4xg37gGaNGlKkyZh\nFz3ncuv58MP3aN26rYJaJVNQExERuUbp6eksXhzDxx9/hp+fPwDOzs488cQzbNq0AbPZzEcfvU9y\nchKxsQcZOHAwWVlZJCWdoUWLCH744Xt+++0X0tJSycnJISnpDJMmPU98/EleemkayclJuLt78PTT\nkwkPb8Hx48eYNWsGmZkZFBcX89e/jmPgwMFXVHNISCgdO3Zmy5ZNNGkSxtatm3n77TfIz8/H1dWN\np56aSIsWLfnww/cs9Tz22EP07NmbNWt+JjExgbZt2/PCCy8ya9Z0jh+P429/e5jJk19gwIBetujm\n65KeURMREblGO3bswN/fn9DQRlbtjo6O9OzZG6Ox7K/bdevW8tprs7nzzlGWY4YPH0lERCTjx/+d\n6OjRVue/+upLDBgwiJiYr7nvvgeYMWMKAG+//RY9evTis88W889/TuHll2dQXFx8xXUXFxdjb29P\nbm4uU6ZM4vHHn+bzz7/knnvu44UXnqO0tLTcOWvX/sqbb77DF198xdatm9m1aweTJ08FYM6c9yx3\n7KRyKKiJiIhco/T0dLy9fSo8rmXLVnh5eV3WNQsKCti2bbPlGbJevW7g/ff/D4CXX/43o0bdC0Cb\nNu0oLCwkOTn5imo+cGA/O3duo3v3KPbu3Y2fnz9t2pSFrD59+pORkU5iYkK58/r06Y+joxPOzs4E\nB4dw+vSpK/pcuTIa+hQREblG3t7eJCcnVXich4fHZV8zKyuT0tJSy0QEg8GAi4sLABs2rGP+/A9J\nS0vHaDRgNpsxm8vf/fqz6dOfs0wm8PHxYfr0lwkIqM+uXTtwd7euzc3NnbS0tHLX+OPECKPReMG7\nblJ5FNRERESuUbt27UhNTeHAgf2Eh7ewtBcXF/Phh+9x//1jr/iaHh6eGAwGMjIy8PLywmw2Ex9/\nkvr1A5kyZRLTp8+ie/eeFBYW0r9/1GVd84+TCf7Ix8eXjIwMy2uz2UxWVgY+PhXfJRTb0tCniIjI\nNfLw8GDUqPuYOXMKJ0+eACA/P59XX32RQ4cO4OTkdMnzTSYT2dlZVm0ODg507tyNpUu/Bcruok2Y\n8A/y8vLIy8ujRYuWACxa9AX29vbk5eVedf0REZGkpqawe/dOAH76aTl+fv4EBgZd9jXs7OzKfQe5\ndrqjJiIiUgnGjn0YDw9PJk16kpKSEoxGIz173sBTT02q8Nzevfsyd+5bJCTE4+LiammfNOk5pk9/\nnv/9bzEeHh688MJM3N3dGTXqPv7yl3vw9vbm/vvH0qvXDTzzzBO8+uqbV1W7s7Mz06fP4vXXXyU/\nPw8vL2+mTXup3JIil9Kv30DGjRvLxInPER1921XVURkK8os5ciCJ3OwCXNwcaRLuh6NT7Y07BrPZ\nbK7uIq5UUtLVJXY/P/erPrcuUn9YU39YU39YU39YU39YU39Yq67+2PJ7HFvXxVFcdP65OZO9kQ7d\nQ+nYI7TK6znHz8/9qs+tvRFTRERE5Kwtv8ex8Zej5dqLi0ot7dUZ1q6WnlETERGRWq0gv5it6+Iu\neczWdXEUFlz5WnPVTUFNREREarUjB5KshjsvpLiolMP7K15CpabR0KeIiIjUKjn5RRw/nc3x01kc\nP51F6tE0LmeFutycQpvXVtkU1KRGeOyxh7CzM/HWW3Mv+H7Pnp3461/HMWbMXyvtM1988QV27txO\nTMzXV12XiIjYVkZOIXGnMok7lUXc2XCWnJFvdUw9wOMyBgldXB1sVKXtKKhJrfDNNz9YTVmvDP/4\nx4Sr2htPREQqn9lsJi2r4Gwgy7L8Nz27/F0we5ORhn5uhAa4ERLgTqC3C798ueuSw58meyNhLfxs\n+RVsQkFNagVf33qVfs0/boMiIiJVx2w2k5KZT9ypLI6dOh/KsnKLyh3r5GBHaIA7ofXdCQ1wJyTA\njfq+LtgZre+gZXcPveCsz3M6dA/FwbH2xZ7aV7Fcl/489Llw4QIWLvyM7OwsIiIiefjhxxg37i88\n//x0Bg266YLDmj/9tJwXXniWRYuWEBgYVO6Y/fv3M2nSZA4fPoSvrx9jxlz5li8iImLt3J2yo4lZ\nxJ3O5FhiWTjLzisfylydTJZAdu6/ft7OGC9j4d1zS2/UxHXUroWCmtQ669b9xttvv8m99/6FIUNu\n4ejRI7z++svXdM2ioiIeeeQRPDy8+M9/PsJoNPLBB//h2LGjhIU1q6TKRUTqvvTsAo6czmbHgTMc\nO5XFsVOZF7xT5upkolGgB43OBrJG9d3x9XS6ot0Q/qxjj1Bad2zA4f1J5OYU4uLqQFgLv1p5J+2c\n2lu5XLd++GEpgYFBPPzwowCEhISSnp7Ga6+9dNXX3LZtCwkJCUyZ8qJlQ+Xnn5/B0KE3VkrNIiJ1\nUXZeEcdOZXI0MYtjiZkcO5VFWlZBuePO3SlrVL8smFVGKLsYB0cTEW0DK/261UVBTWqduLhjRERE\nWrW1bNnqmq557NgRAJo1a25pc3NzIzS00TVdV0SkrigoKiHuVBZHEzM5mlg2hHkmPa/ccU4OdjQP\n8SbI14XGgR6E1nfHz0ah7HpQI4LanDlzOHXqFB4eHgwdOpSIiIjqLklqsNzcHJycnKzaXFxcrvGa\nuRgMBhwdHa3anZ2v7boiIrVRaamZhOQcjpwNZUcTMjmZlEPpn7YHtzcZCfF3o1GgB40D3Wkc6EGA\njwsB/h7a+7SS2DSoHTx4kPHjxzNmzBhGjx4NwEsvvcSOHTswGAxMnjyZNm3aAODk5ERRURH+/v62\nLEnqAGdnZ3Jzc63asrIyrV4bDAbMf/oDJS+v/L/8znFycsJsNlNQUGAV1rKzs/D29q2EqkVEaqZz\nD/sfScgsC2YJZUOYBUUlVscZDBDs70bjQPeyYFbfgwZ+rpjstMmRLdksqOXm5jJjxgy6d+9uadu4\ncSNxcXHExMRw+PBhJk+eTExMDHfeeSdeXl4kJSXxySef8OSTT9qqLKkDgoNDiY09aNW2c+d2q9cu\nLq5kZVn/a27v3t0XvWZISCMA9u/fS9u27QFISjpDXNwxBTURqVMKCks4diqzLJglZHI4IeOCa5XV\n83SicaAHjQM9aBLkQWiAO44OdtVQ8fXNZkHNwcGBefPmMW/ePEvbunXrGDBgAABhYWFkZGSQnZ3N\nkSNH6NKlCx4eHhQW1r7tHaRyFBcXkZKSXK7d3d16Y5CBAwezevVK5s37D4MH30xs7EGWLv3O6pjw\n8BYsXryQRYsW0qNHTzZt2sCePbsu+tkdO3aiXr16vPPOW0yYMAmAd999Gy8v70r4ZiIi1aPUbOZ0\nau7ZQJbJkfiMCw5hujiaaBzkQZNAD8t/PWrhKv51kc2CmslkwmSyvnxycjKRkecfAvfx8SEpKYn8\n/HwmTZqEyWTioYceqvDa3t4umExXl+r9/Nyv6ry6qqb0h4ODiY0bNzJs2OBy77377rsAuLo64ufn\nzsiRQ0lKimfBggUsXLiAtm3b8txzkxkzZgzu7k74+bkzatQdxMbu45NPPuCDD/5Dv379eOaZpxk3\nbhw+Pq74+bnj5GSPnZ3xbB+4M3fuXKZNm8bDD/8Ff39/xo8fz6pVq8jNza0x/VTVrtfvfTHqD2vq\nD2s1oT9y84s4EJfG/rg09selciAujZw/rVdmNBpoEuRJeKg34SHehId6E1TPDaOxch/2rwn9URcY\nzH9+kKeSzZkzB29vb0aPHs3zzz/PDTfcYLmrdvfdd/PSSy/RuHHjK7rm1T6g6Ofnrocb/6Au9Udi\nYgJ33DHUsuDt1ahL/VEZ1B/W1B/W1B/WqqM/zGYzp1JzORxfNnx5OD6D+KQc/vyXupebA2FBnjRp\n4EFYkGeVDGHq58PatYTWKp316e/vT3Ly+aGtM2fO4OdX+/bdEhERqWoFRSUcTcgkNj6D2PiyYJaT\nb71fsZ3RQEiAO00beBLWwIOmDTzx8XC6yBWlNqjSoBYVFcWcOXOIjo5mz549+Pv7a79FERGRC0jL\nKigLZScziI1P5/jpbEpKre+Xebo50DTIk7CzwaxRfXfsr/LRIKmZbBbUdu/ezSuvvEJ8fDwmk4nl\ny5czZ84cIiMjiY6OxmAwMHXqVFt9vFxnAgOD+O23zdVdhojIVSktNROfnMOhk+mWcJackW91jMEA\nIQFuNGvgRVjDsrtlvh5aSLaus1lQa9WqFZ9++mm59gkTJtjqI0VERGqFouISjiZmcehkOgdPlA1l\n5hVYD2M6OdgRFuRB04ZeNG3oSZNAD5xr8Z6VcnX0/7iIiIiN5eQXEXsyg4Mn0zl0MoNjiZkUl1gP\nY/p6ONGsoSdNG3rStIEnDf0qfyam1D4KaiIiIpUsPbuAgyfSOXAinUMn0svNxjQADf1caRbsRbOG\nnjRv6KWH/uWCFNRERESugfnsorLrdiVy8EQ6B0+kczrNess6k52BxoEeNGvoRfPgsof/XZ3sq6li\nqU0U1ERERK7AufXLDpwNZQdPpJOaWWB1jKODHc0aeNI82IvmwV40DtRsTLk6CmoiIiKXcC6Y7T+e\nzoHjaew/nk5mjvV2h27O9jQ9G8zCQ7wICXDDzqjNyuXaKaiJiIj8gdlsJjEl1xLKDpwoH8w8XOxp\nHuJNeLAX4cFetI2oT0pKdjVVLHWZgpqIiFzXzGYzp9Py2B+Xxr64NA4cTyMz13p/TA9XB1qElIWy\n8BBvAn1drNYv0+xMsRUFNRERue4kZ+Sx79zm5cfTScuyfsbMEsxCvGkR4kV9HxctLCvVQkFNRETq\nvPTsAssds/3H00hKt171383Znhah3kSEeNEi1FvBTGoMBTUREalzcvOLOXAijb3HysJZQnKO1fvO\njibCg72ICPUmItSbID9XjApmUgMpqImISK1XVFzK4fgM9salse9YKkcTsyg1n19i1sHeSPNgLyJC\nvGkR6k1ogLueK5NaQUFNRERqnVKzmROns9kbl8q+Y2kcPJFOYXGp5X2jwUDThp60DPWmZSMfmgR5\nYLLTchlS+yioiYhIrZCWVcCeo6nsPZbKnmOpZP1pZmYDP1dahvoQ0ahs2QxtYC51gX6KRUSkRioo\nLOHAiXRLOIv/03Nm3u6OtGzkTWQjHyJCvfF0c6ymSkVsR0FNRERqhHPDmbuPprDnaCqx8RkUl5x/\nzszR3o7wEC8iG/vQqrGPZmbKdUFBTUREqk1WbiF7jqWy+0gqu4+mWu0AYAAa1Xe3BLOwBp56zkyu\nOwpqIiJSZUpKSzmakMWuIynsPprCscQszH9439vd0RLMWjbywc3ZvtpqFakJFNRERMSm0rML2HUk\nhV1HUtl7NJXcgmLLeyY7A82DvWjV2JdWTXxoUM9Vw5kif6CgJiIilaq01MyRhEx2Hklm5+EUjp+2\n3qw8wNuZVk18ad3Eh/Bgbxwd7KqpUpGaT0FNRESuWWZuIXuOpLLzSAq7j6SQk3/+rpmDyUhEqDet\nw3xp1cQXfy/naqxUpHZRUBMRkStmNpuJO5XFjthkdh5J4WhCptWzZv5ezrQJ86VNmC/hIV7Ym3TX\nTORqKKiJiMhlKSgqYd+xNLbHJrP7aCqpmec3NjfZGQgP8aZNk7JwFuDjUo2VitQdCmoiInJRqZn5\n7Dicwo7YZPbFpVH0h22avNwcaNu0Hm3CfGkZ6qNnzURsQEFNREQsSs8OaW47lMzO2GSOn7GeCNA4\n0J22YfXo0zkEdwejZmiK2JiCmojIda6ouIR9cWlsP5TM9thk0rPPLzrraG9HZGMf2p593uzcNk1+\nfu4kJWVVV8ki1w0FNRGR61B2XhE7YsuC2e4jqRQUlVje83Z3pF2zerRvWk8TAUSqmYKaiMh14kx6\nHtsOJrHtUDKHTqZj/sM0zRB/t7Jw1syPkAA3DWmK1BAKaiIidZTZbObEmWy2Hkxi68FkTiadf97M\nzmggPNSL9s38aNe0Hr6eTtVYqYhcjIKaiEgdUlpqJjY+g60Hk9h2KImk9PNLaDg52NEmzJf2zfxo\n3cQXFyf9FSBS0+l3qYhILVdcUsq+uLSz4SyZzJzzkwE8XOxp18yPDs39iAj1xt5krMZKReRKKaiJ\niNRChUUl7DmayuYDZ9gem0LeHzY6r+fpRIfmZeGsaQNPjEY9byZSWymoiYjUEgVFJew6nMLmA2fY\ncTiFgsLzMzUb+LnS8Ww4C/bXZACRukJBTUSkBssrKGbXkRQ27z/DziMpFBad3xkgNMCdTi386Bju\nT31t2SRSJymoiYjUMHkFxeyITWbT/jPsPppqtW1TkyAPOoX70zHcDz8v52qsUkSqgoKaiEgNUFBY\nwo7DyWzaV3bn7Fw4MwDNGnpawpmPh5bRELmeKKiJiFSTwqISdh1JYeO+M+w4nGw1rNmsoSddIgLo\n0NwPb3fHaqxSRKqTgpqISBUqKi5l99EUNu07w7bYZKsJAWENPOjcIoDOLfwVzkQEUFATEbG50lIz\n+46nsWHvabYcSLJaSqNxoDudWwTQqYUf9Tz1zJmIWFNQExGxAbPZzJHETDbsOc2m/WfI+MMitCH+\nbnSO8KdzRAD+mhAgIpegoCYiUonik7LZsO80G/aettq+yd/bmW4tA+jaMoBAX9dqrFBEahMFNRGR\na5Samc/6vadZv+e01cbnXm4OdIkoC2eN6rtrEVoRuWIKaiIiVyE3v5gtB86wbs8pDhxPx3y23dXJ\nRMdwf7q2DCA82EvbN4nINVHkznQ6AAAgAElEQVRQExG5TMUlpew+ksq6PafYHptsWevMZGekfbN6\ndIsMoHUTX0x22vhcRCqHgpqIyCWcmxTw5a9HWbP1JNl5RZb3WoR40T2yPh3D/XFx0h+nIlL59CeL\niMgFpGTk8/vuRH7ffYrTaXmW9qB6rnSPDKBby/r4emqXABGxrRoT1JKSkhgxYgSrV6/GZKoxZYnI\ndaSgqIStB5L4bVci++PSLM+debg60LdjMO2a+BAS4KZJASJSZWyaiA4ePMj48eMZM2YMo0ePBuCl\nl15ix44dGAwGJk+eTJs2bQD4+OOP6dy5sy3LEREpx2w2c+hkBmt3JbJp/xnyz+4UYLIz0qF5PXq0\nCiSysTf1AzxJSsqq5mpF5Hpjs6CWm5vLjBkz6N69u6Vt48aNxMXFERMTw+HDh5k8eTIxMTF88803\n3HjjjSxcuNBW5YiIWEnOyOP33af4fdcpzqSfH9oMC/IgqnUgnSP8cXWyr8YKRURsGNQcHByYN28e\n8+bNs7StW7eOAQMGABAWFkZGRgbZ2dns2LGD48ePs2/fPr7//nuGDRt2yWt7e7tgMtldVV1+fu5X\ndV5dpf6wpv6wVtf6o7CohHW7EvlxYxw7DiVb2n08nOjXKZh+nYIJDrj4d65r/XGt1B/W1B/W1B+V\nw2ZBzWQylXvWLDk5mcjISMtrHx8fkpKSmDJlCgDx8fHcfPPNFV47LS33qmry83PX0MUfqD+sqT+s\n1aX+OH46i193JLJ+7yly8sv22bQ3GenQ3I+o1vVpGepjWe/sYt+5LvVHZVB/WFN/WFN/WLuW0Fqt\nT+2bzWar1y+//HI1VSIidU1ufhEb9p7mlx2JxJ0+/xdGaH13ercJpGvLAFw0tCkiNVyVBjV/f3+S\nk88PN5w5cwY/P7+qLEFE6jCz2cyB4+n8ujOBzQeSLAvSujqZ6NayPr3aBhJyiaFNEZGapkqDWlRU\nFHPmzCE6Opo9e/bg7++Pm5tbVZYgInVQZk4ha3clsmZHAmf+sOZZRKg3vdoG0rG5H/ZX+VyriEh1\nsllQ2717N6+88grx8fGYTCaWL1/OnDlziIyMJDo6GoPBwNSpU2318SJSx5nNZvbHpbF6ewJbDyZR\nUlr2KIW3uyNRrQPp2SYQfy/naq5SROTa2CyotWrVik8//bRc+4QJE2z1kSJyHcjKLWTtrlOs2R5v\n2THAYIB2TevRu10QbZr4aiN0EakztAWAiNR4ZrOZgyfSWb09gS0HzlBccv7uWa82gfRuG4SPh7Zz\nEpG6R0FNRGqs3Pwi1u46xert8SSmlC3LYwDahPlyQ7sg2oT5Ymc0Vm+RIiI2pKAmIjXO8dNZrNoa\nz/q9pygsKpu56eXmQK82QfRuG6TN0EXkuqGgJiI1QlFxKZsPnOHnrfHExmdY2ls28qZv+4a0a6a7\nZyJy/VFQE5FqlZKRz+rt8fyyI4Gs3CIAnB3tiGoVSN8ODQj0da3mCkVEqo+CmohUObPZzN5jaaza\nepLtscmc26SkoZ8b/To2oHvL+jg6aN0zEZGLBrWEhITLukBQUFClFSMidVtBYQm/7znFT5tPWCYH\n2BkNdIrwp1+HBjRt4InBoKU1RETOuWhQGz58OC1atCi3H+cfHThwgI0bN9qkMBGpO5Iz8li1pWx4\nM7egbFN0LzcH+rZvQO+2QXi6OVZzhSIiNdNFg1p4eDjz58+/5Mn33ntvpRckInXDubXPftp8kq2H\nkizDm2FBHgzoFEzHcD9MdpocICJyKRcNaud2FThx4kS59+zs7AgICLjgzgMicn0rKi5h/d7T/LT5\nJCfOZANlw5udW/ozoGMwTYI8qrlCEZHao8LJBA899BBxcXE4OztjNBrJzc0lICCAnJwcpk+fzqBB\ng6qiThGp4TJzClm19SQ/b4u3zN70cLGnT/sG9GnfAC8Nb4qIXLEKg9oNN9xAVFQUvXr1AmDt2rVs\n3LiRe++9l0ceeURBTeQ6l5Ccw4pNx/l992mKS8oWpw0NcGdAp4Z0iQjA3qThTRGRq1VhUNu1axeT\nJk2yvI6KiuL999/niSeewGTS6h4i1yOz2cy+uDSWbzzBriMpQNnWTu2a1mNQl2CaB3tp9qaISCWo\nMGmVlpayYMECunTpgtFoZNu2baSnp7N169aqqE9EapDiklI27jvNio0nOH72+TMHk5Go1oEM7BxM\nfR+Xaq5QRKRuqTCovfrqq8yePZuYmBhKS0sJCwvjtddeo7CwkBdffLEqahSRapabX8TP2+JZueUk\n6dmFAHi4OtC/Q9nzZ+4uDtVcoYhI3VRhUAsODua1116rilpEpIZJyyrgx00n+Hl7PAWFJQA0qOfK\njZ2D6RYZgL1JuweIiNjSRYPaqFGjLvmMyWeffWaTgkSk+iWm5PDFqlhWbT5BSWnZAmgRod4M7hpC\nq8Y+ev5MRKSKXDSoPf7445Zfm81mZs2axeTJk6ukKBGpHocTMli2/jjbDiZhpmyCQKcW/gzpGkLj\nQK1/JiJS1S4a1Lp06WL12t3dvVybiNR+ZrOZ3UdTWbY+jv3H0wEw2Rno3zmEPm0CCdAEARGRanPZ\n62toqEOkbiktNbNp/xmWro+z7CDg7GhHn/YNGNgpmGaN65GUlFXNVYqIXN+0EJrIdaa4pJR1e06x\ndF0cp9PyAPB0dWBg52D6tGuAi5P+WBARqSku+ify008/bXUX7fDhwzzzzDOW16+++qptKxORSlVU\nXMJvu8oCWkpmPgD1PJ24qXsoUa3qawaniEgNdNGg1qNHD6vX3bt3t3kxIlL5CopKWLM9gR82xFnW\nQAv0deHm7qF0bRmAnVFbPImI1FQXDWojRoyoyjpEpJLlFRSzautJVmw6YdkkvaGfG7dGNaJjcz+M\nRj13KiJS01X4MMp7773HBx98QHZ22cPGZrMZg8HAvn37bF6ciFy53Pwiftx8kp82nyAnvxiAxoHu\n3NqjMW2b+mpikIhILVJhUPv666/5+uuvqV+/flXUIyJXKTe/mB83n2DFphPkFZQFtOYNPbk1qjEt\nG3kroImI1EIVBrVmzZpRv3597Oz0oLFITZRXcDagbTxB7tmAFhHqzdCoRoSHeFdzdSIici0qDGrD\nhw9n6NChREZGWoW1WbNm2bQwEbm0vIJiftpykhUbj1uGOFuEeDGsZ2MFNBGROqLCoDZr1iyGDRtG\nQEBAVdQjIhU4N0nghw3nA1rzYC+G92xMi1AFNBGRuqTCoBYSEsJjjz1WFbWIyCUUFJaw8mxAy84r\nm8XZtKEnw3s2JiJUz6CJiNRFFQa1tm3bMnv2bDp06GA19Kl11USqRlFxKau3x/P978fIPLvMRlgD\nD4b3bKJJAiIidVyFQW3Tpk1W/4WyfT8V1ERsq6S0lHW7T/PNb0ctOwk0DvRgRK/GRDb2UUATEbkO\nVBjUPv3006qoQ0TOMpvNbDmQxP9+PUJiSi4ADeq5clvvJrRrVk8BTUTkOlJhUPv999/5/PPPycrK\nwmw2W9rnz59v08JErjdms5k9x1L5cs0R4k5lAWV7cY7o1YSuLQO0k4CIyHWowqD2wgsv8Mgjj2jB\nWxEbij2ZwZdrDnPgRDoAnm4ODO3RiF5tgzDZaS9OEZHrVYVBrVGjRtr3U8RG4pOy+XLNEbbHJgPg\n6mTipm6h9OvYEEd7LTItInK9qzCo3XnnnTz77LO0b98ek+n84cOHD7dpYSJ1WVpWAd/8doRfdyZi\nNoOjvR0DOwczuEswLk721V2eiIjUEBUGtXfffRdnZ2cKCwstbQaDQUFN5CrkFRSzbMNxVmw8TmFx\nKXZGAze0D+LWqMZ4ujpUd3kiIlLDVBjU7O3tNfNT5BoVl5Tyy44EvvntKFln10Lr2NyP2/uEUd/H\npZqrExGRmqrCoNavXz/Wr19Phw4drIY+jUY94CxSEbPZzNaDSSxefZjTaXkANG3gyZ19m9K0oWc1\nVyciIjVdhUFt7ty55OWV/QVjMBgwm80YDAb27dtn8+JEarPYkxnE/HyIw/GZAAR4OzOyT1M6NNda\naCIicnkuGtRSUlLw9fVl27ZtFz353DEicl5Seh6LVh9m8/4zALi72DOsZ2N6a6kNERG5Qhf9W+OJ\nJ56o8OTLOUbkepFXUMyXaw7z7LwNbN5/BgeTkVt7NOLlh7vTr0NDhTQREbliF72jtm3bNvr372+1\nG8E554ZAU1JSbFqcSG1QajazdlciX605QkZO2ezobpEBjLwhDB8Pp2quTkREarOLBrVdu3ZVWRFb\ntmxh4cKFFBUVMXbsWFq3bl1lny1yLQ4cT2PhyljiTpdt+dQkyIO7+zcjrIEmCoiIyLWrcDLBtTh4\n8CDjx49nzJgxjB49GoCXXnqJHTt2YDAYmDx5Mm3atMHNzY2ZM2dy4MABNm7cqKAmNV5Seh6Lfo5l\n84EkALzdHbmjTxhdWgZg1EQBERGpJDYLarm5ucyYMYPu3btb2jZu3EhcXBwxMTEcPnyYyZMnExMT\nQ3h4OGvWrOHDDz9k5syZtipJ5JrlFRSzdH0cyzeeoLikFAeTkSHdQhncNURbPomISKWzWVBzcHBg\n3rx5zJs3z9K2bt06BgwYAEBYWBgZGRlkZ2dz+PBhevfuTevWrXn77beZMmWKrcoSuSpms5lN+8/w\nxcpDZGSXPYfWPTKA2/UcmoiI2FCFQS0jI4N3332XpKQk/vWvf7Fq1SratWuHj4/PpS9sMlktkAuQ\nnJxMZGSk5bWPjw9JSUlkZGQwZcoUcnNzGTp0aIVFe3u7YDJd3d0LPz/3qzqvrlJ/WLtQf8QnZfPu\nV7vYfqhsmLN5iBcPDm9Ni9BL/x6oC/TzYU39YU39YU39YU39UTkqDGrPPfccnTt3tqynVlhYyMSJ\nE63ulF2tczNKe/fuTe/evS/7vLS03Kv6PD8/d5KSsq7q3LpI/VGmJDeX7C2bcCzOp8DkhFvHzti5\nuFBQVML36+L4YUMcxSVmXJ1MjOwTRq+2QRgNhjrfd/r5sKb+sKb+sKb+sKb+sHYtobXCoJaamsp9\n993Hjz/+CMDgwYP57LPPrurD/P39SU5Otrw+c+YMfn5+V3UtkcqQ8t0SUpd9j7mgwNJ2ZuHnFHTt\ny/ycEJIz8gHo2SaQkX3C8HDRxukiIlJ1LmsFzqKiIsuWN8nJyeTmXt0draioKJYvXw7Anj178Pf3\nx83N7aquJXKtUr5bQsrXX1mFNABzQQEOv/xAs6Mbaejnxj9Hd+CBmyIU0kREpMpVeEftnnvuYeTI\nkSQlJTFu3Dh27drFs88+W+GFd+/ezSuvvEJ8fDwmk4nly5czZ84cIiMjiY6OxmAwMHXq1Er5EiJX\nqiQ3l9Rl31/ymF6ZewmLfhB7V9cqqkpERMSawXyhrQf+5NSpU2zbtg0HBwdat26Nv79/VdR2UVc7\n7q0xc2vXc39k/LqG0598XOFxAWMewLPn5T8/WZdczz8fF6L+sKb+sKb+sKb+sGbTZ9Qef/xx3nzz\nTYYMGXLVHyJS0xRnZFzecenpNq5ERETk4ioMag0bNmTx4sW0b98eB4fzz+gEBwfbtDARW4rPM3I5\nq5+ZvLxsXouIiMjFVBjUli5dWq7NYDCwcuVKmxQkYktpWQUsWHGAvQdNPGow4WAuvuixBkdH3Dp2\nrvCaPXt24uabh/LPf55fqHnr1s189NH7vP32+5VS95Xo2bMTX331Pf7+AVbtjz32ECdOHMf17DN3\nzs4ujBv3KJ07d7vk9das+Zm1a39h8uSpPPbYQ9x663A6dOjEk08+xqef/rdSav7ww/dYtOgLfHx8\nASgtNdOpU2cee+wJnJwuHqmTks5cVh3Hjx8jNTWVdu06VEq9IiJVpcKgtmrVqqqoQ8SmSs1mVm+L\nZ/Hqw+QXluDk7Exu5z44bPzpouf4DLkZO2fny7r+9u1bOXhwP82bt6iskm1i/Pi/M2jQTQBs2LCO\n55+fxNKlqzAaLz4B/IYb+nLDDX2t2vz8/CstpJ3Tp09/Jk16HiibaT5lyiQ+/ngejzzyt4uec7l1\nrFmzmpKSYgU1Eal1KgxqzzzzzAXbX3311UovRsQW4pNz+GTZfmLjy55La9+sHvcMbI6Pxw2kBHmU\nW0fN4OiIz5Cb8b2l4l0yznn44ceYPfv1C95BKy0tZd68/7BmTdk/eiIjW/PkkxNxdnbmscceonXr\ntvzyy89MmvQ8S5b8j3r1/Ni9eydHjx7m1ltHEBTUgEWLvji7f+7LREREkpqawsyZL3DqVAJFRUXc\nfvudREePvqJ+6dChE9nZ2WRkZODt7U1s7CH+/e9ZZGRk4ODgyKRJz9CiRTuWLv2W5cuX8dZbcy3n\nJiYmEB09gjVrNrB06bf8/vtvuLq6smPHdkwmO6ZPf5kmTcJITExg8uQJZGdn06VLN5KSztCnT39u\nuunWS9Zmb2/P0KEjmDfvPzzyyN/IzMzgtddmERt7EDs7OwYPvpnRo8dcVh0JCfEsWPAx9vb2ZGVl\n8be/PXFF/SQiUp0qXEete/fulv916tSJkpISAgMDq6I2kWtSVFzK178e4YWPNhIbn4GnqwPjh7fi\nsdtaW/bn9L1lKE3+9SYBYx4g5J67CRjzAE3+9eYVhTSAfv0GYDab+fnn8nfoVq36kQ0bfufDDxfw\n6af/JTs7i5iY84tGHziwn08//S+tW7cFYMOG33n11TeZPftdPv98PunpacyfH0OfPv1ZvHghAJ98\n8iFBQUF8/vmXvPnmXN577x1Onz512fWazWa+/34JrVq1wdvbm9LSUl54YTK33XYnn3/+JZMmPcdT\nTz1Fbm7OZV1v/fq1jBhxBwsXfkX79p1YtOgLAN555006d+7GokVL6Nq1B5s3b7zsGouLi7G3L3su\n9r333sHd3Z0vvviKuXM/4H//W8yOHdsvq46ePXvTu3dfRo6MVkgTkVqnwjtqI0aMsHp955138vDD\nD9usIJHKcDQxk4++30d8clnQuKFdEHf0CcPFyb7csXbOznj27H3N08n//veneP75ifTo0cuqfd26\n3xg8+Baczw6j3nTTrSxc+BljxvwVgO7do6yGHjt16oqzszONG4dRWlpKVFTZ9cLCmvLDDwcAePzx\npykpKQGgQYOG+Pj4kpiYQEBA/UvWOHfubD755ENycnIoKCjgueemAWV3yFJSUhgwYBAALVq0JCgo\niH379l7Wd2/UqAktWkQAEB4ezqpVZYF1x47t3H//WAB69+6Dr+/l7USSm5vDV18tsgy5rlu3llde\neQMADw9PbrihL5s2refmm60D9cXqEBGprSoMaqWlpVavExMTOXbsmK3qEbkmRcWlLFl7lGXrj1Nq\nNhPg7cyYIS0ID/G2+WeHh7egXbsOxMR8RqtWbSztaWnpuLufX0PH3d2DtLRUy2sPDw+r67i4uABl\nk3aMRiPOzmWvjUajJZzt27fHchfNaDSSkpJc7vfqhfzxGbVTp07x9NN/x2Aw4OnphZubu2UHknN1\n/bHOS3F1Pb/DiNFoR0lJWS1ZWZm4u3ta3rvUlnGrV69k586yu2T29vb07t2XO+8cBUB6ehru7uf7\nyd3dg+TkpMuuQ0SktqowqLVs2dLyh7fZbMbd3Z0HH3zQ5oWJXKk/3kUzADd2DmZE7yY42ttVWQ0P\nPfQoY8eOJjAwyNLm4+NDZub5ddsyMjIssxuv1vTpU7jrrlEMH347BoOB4cOvfJ3D+vXr06NHL9av\n/5277x5NVlYGZrPZ8vs9PT0dHx9fTp1KvOo6XV1dycs7v+VcSkryRY/942SCP/P29iEzM5369cvu\nGFZGH4qI1AYVPqO2f/9+9u3bx759+9i/fz+bNm3ioYceqoraRC5LUXEpX645zIvztxCfnEOAtzOT\nRncgun+zKg1pAPXq1eO22+7go4/OTyro0aMXy5cvIz8/n+LiYr7//hu6d4+6ps9JT08lPLwFBoOB\nZcu+Iz8/j7y8vCu6Rm5uDps3b6Rx4yYEBgbh5+fPypUrANi1awfJyclEREReU50REZGsWvUjAGvX\n/nrBu2CXo0ePXixZ8j+gLED+8svPV9SHJpOJ7Ozsq/psEZHqVOEdtd27d5OUlETfvn1544032L59\nO3/729/o1KlTVdQncklHEzP5aOk+4pOq7y7an0VHj+bbb7+2vO7btz+HDx9i7NjRmM1mOnToxMiR\n0df0GX/96zgmT34aT09Phg69jaFDb+PVV2cyd+6Hlzzv3DNqULZW2YABNzJixEgMBgPTpr3Ea6/N\n4uOP5+Hk5Mxbb71lea7uao0f/3emTXuOlStX0K1bD1q1amM1vHq5HnzwEf7971mMGnU7RqORe+65\nn5YtW5GYmHBZ50dF9WLatOc4dSqBmTM1Y11Eao8K9/qMjo7m5ZdfJjk5mblz5zJ58mSmT5/O/Pnz\nq6rGcrTXZ+Wozf1xoWfRHrg5gmYNr34ngdrcH7ZQWf3xx+HUv/71Pu6//wF69epzzdetavr5sKb+\nsKb+sKb+sGbTvT4dHR1p1KgRMTEx3HnnnTRt2vSSi2OK2NqxU5l8+H3NuosmF/bOO2+Rl5fHhAmT\niIs7RlzcUcLDI6q7LBGRWqPCoJaXl8eyZcv46aefePTRR0lPTyczM7MqahOxUlJaytJ1cSxZe4yS\n0rK7aH+5KYLmwdqPs6a66657mDFjCnfdNRyj0ciTT04st7WViIhcXIVB7cknn2T+/Pk88cQTuLm5\nMWfOHMaMGVMFpYmcdyYtl3nf7eVwfNk/EgZ0asjtN4TpLloNV69ePasdDURE5MpUGNS6detGq1at\ncHNzIzk5me7du9Ohg/bLk6phNpv5dWciX6w8REFhCd7ujjxwcwSRjXyquzQRERGbqzCozZgxgxYt\nWjBw4ECio6Np1aoVS5YsYfr06VVRn1zHMnML+WTZfrYdKlt7q3MLf+4dFI6bc/ndBUREROqiCmcF\n7N27lzvuuINly5YxYsQI3nzzTeLi4qqiNrmO7TyczJQPN7LtUDLOjnY8eEtLxg2LVEgTEZHrSoV3\n1M6t3rF69Woef/xxAAoLC21blVy3CgpLiPk5ltXb4gEID/Zi7C0R1PO8tvW8REREaqMKg1rjxo25\n6aab8PHxISIigq+//hpPT8+KThO5YkcTM3n/272cTs3FzmjgthuaMKhzCEbjlS+QKiIiUhdUGNRm\nzpzJwYMHCQsLA6Bp06a8+qpW9pbKU2o288OG4/zvlyOUlJppUM+VB29tSUjA1S8QKCK2NWbMGM6c\nSeKTTxZe8P3jx+MYNep2nnjiaW6//a4qrk6k7qjwGbXs7GyWLFnCs88+C8CZM2coLi62eWFyfcjI\nKeSNmO0sXn2YklIzAzsFM2VMJ4U0kRpu+PDhHD4cS2zsoQu+v2LFMuzt7Rk4cHAVVyZSt1QY1J57\n7jkCAwM5ceIEUPZ82sSJE21emNR9e46mMvWjjew5loabsz2P39GGuwc0w96ktdFEarpBgwbh4uLK\nihVLL/j+ihXLiIrqhYeHHpURuRYVBrXU1FTuu+8+7O3LZtsNHjyY/Px8mxcmdVdxSSmLVx/m9Zjt\nZOYU0iLEi2kPdKFNWL3qLk1ELpOzszN9+/bnxx+XU1paavXerl07SEiIZ8iQWwFYvXolY8feS79+\nPRg8uC9PPPEoR47EWo5///25DBs2iK1bNzN69B307dudu+++jbVrf7W67pdfxnD33bfRp083brll\nADNmTCEtLc3y/iOPjOXFF1/glVdepH//KLZu3WzDHhCpGpe1aWdRUZFlU+Xk5GRyc3NtWpTUXcnp\nebzy2VaWro8DA4zo1ZgJ0e3xdnes7tJE5AoNGXILSUlnygWi5cuX4evrS9eu3Tl69AhTpvyTLl26\nsWDBIubO/QB7ewcmTnzS6jGarKwsPvnkQyZNep4PP/yUkJBQpkyZRHJy2TqKX34Zw+zZr3PXXfew\nYMEipk2bxZ49u5g48Qmrz96+fRt2dnZ89tliWrZsZftOELGxCoPaPffcw8iRI4mNjWXcuHEMGzaM\nsWPHVkVtUsds3n+GqR9v4nBCJt7ujkwc1YFboxprVqdILdW2bXuCghqwfPn54c/i4mJ+/vlHBg4c\ngslkIigoiE8+WcgDDzxEUFADmjQJ4447oklMTODEieOW8woLC7n//rG0atWGJk2a8sQTz1BQUMBv\nv60BYOHCzxk4cDDDh99Ow4bBdOzYmX/84yn27t3N/v17LdfJzEznH/94ivr1A3Fycqq6zhCxkQpn\nfd5000106NCBbdu24eDgwPTp0/H396+K2qSOKCwqYeGq82ujtW9Wj7/cFKHFa0VqOYPBwJAht/DF\nFwuYMGESjo5OrF+/loyMDG666RYAHB2diI09yGuvvcSJE8fJz8+zDJVmZmZaXS8y8vwdsPr1A/H0\n9OT06VNkZmaQmBjPPffca3V8RETZ8QcO7KdFi5YANGkSZnlUR6QuqPCO2uOPP079+vUZMmQI/fv3\nV0iTK5KQnMPM+ZtZvS0ek52BUQOa8dhtrRXSROqIwYNvJi8vl19+WQ2UDXuGh0fQpElTAH76aTnT\npj1HaGhjXnnldT7++HMmT55a7jomkwlHR+s7YM7OLmRlZZGTkwOAi4ub1fsuLi4AlvfLznGttO8m\nUhNUeEetYcOGLF68mPbt2+Pg4GBpDw4OtmlhUvtt3Heaj5fup6CohABvZ8YNa0VofS27IVKXBAYG\n0a5dB378cTlRUb1Yu/ZXHn3075b3f/ppOY0aNWHixGctbQcPHih3neLiYoqKiqzuhuXm5uLu7o6r\na1lAy8nJtjrn3Gs3N+sAJ1KXVBjUli4tP/XaYDCwcuVKmxQktV9xSSn//TmWnzafBKBrywDuGxSO\ns2OFP24iUgvddNOtvPbaLFat+gmzudRq7bSiouJyu9n8+OOys78yW7Xv3Lmdjh07A5CYmEBmZgYh\nIaF4eHgQFNSAXbt2MGLESMvxu3fvBCAioqUNvpVIzVDh35yrVq2qijqkjkjNzOc/3+zmcHwmdkYD\n0f2b0a9DA8usYRGpe850TiEAAB+aSURBVG64oR+vv/4K8+bNLbd2WsuWkSxY8H+sW7eWhg2D+e9/\nv8DT0wsoC1pNmzYDyoY+P/54Hg4ODri6uvLuu2/j7OxC7959ABg16j7eeONVWraMpHv3npw8eYLZ\ns1+nY8cuNGsWXuXfWaSqXDSoZWdnM3fuXI4cOULnzp25//77MZl0R0Qubu+xVN5bsoes3CK83R0Z\nP7wVYQ202KVIXefi4kKfPv1Ztuw7y9pp59x992ji4o4xdepknJycuOWWYfz1r+PIyEjngw/exdm5\n7Dkzo9HI2LEP869/vczx48cIDAzixRdftQx7Dh9+OyUlxSxaFMPbb7+Ju7sHvXv3Yfz4v5erR6Qu\nMZjNZvOF3pgwYQL+/v507tyZFStWEBAQwOOPP17V9V1QUlLWVZ3n5+d+1efWRZXVH6VmM0vXxfG/\nX49gNkPLRt48NDQSDxeHik+uQfTzYU39YU39Ye3/27vT+KjKu//jn0kmCyGBsCRhE6XsJhCgaBtQ\nFhtBFAWjkKhIuSsCfxprRSqUWrFFoCACCuVWor7aEjQhYBGLEhFEUSkIBLKApqCyxEAWQvZ1cv4P\nuJkyEJZAJjOZ+b4fMSdnrvPLz+sFX892NWQ/Vq9eRUJCPNu3f9Ug4zmC5oct9cNWUND135992VNk\nWVlZLFmyBIAhQ4YwadKk6z6IuK7Simre/OAQB4/mAzB60C2MvUPvRhMREWkIlw1qF17m9PTU2oty\nqWOnivnrP9PIK6ygua+ZyaNvJbybloESERFpKJcNahff/K2bweVCOw/+yJqPM6mx1HJzSADTHwwj\nKLCZo8sSkSZoypTpTJky3dFliDilywa1lJQUhg0bZv2cn5/PsGHDMAwDk8nEjh07GqE8cTY1lloS\ntv2H7fvPrTIwJLwDj93dHS+zzrqKiIg0tMsGtS1btjRmHdIEFJVV8frGdL45fhazp4kJI3oyJLyD\no8sSERFxWZcNah07dmzMOsTJHT9dzIoNaeQXVdCyuTe/jupDN716Q0RExK70YjS5qr3f5PDm5kNU\nVdfSpX0AsVF9aRXg4+iyREREXJ6CmlxWrWGwcef3/OurHwCICG3HpFE9dT+aiIhII1FQkzqVV9YQ\n98EhDhzJw2SC8cO7MeK2m/T0r4iISCNSUJNLnC4oY8WGNH7MK8XPx8y0saGEdWnj6LJERETcjoKa\n2Mj4/gyvv59OaUUN7dv48ZuH+xLSys/RZYmIiLglBTWx2rbvJO98kolhQL9ubXny/ltp5qMpIiIi\n4ihO8a9wSkoKSUlJWCwWHn/8ccLCwhxdkluprTVI2P4fPtl7EoDRg25m7J0/wUP3o4mIiDiUhz0H\nz8zMJDIykvj4eOu2BQsWEB0dTUxMDKmpqQA0a9aMuXPnMmnSJPbu3WvPkuQiFVU1rHwvjU/2nsTT\nw8Tk0b2JGtJVIU1ERMQJ2C2olZWVMW/ePCIiIqzb9uzZw7Fjx0hMTGT+/PnMnz8fgF69elFdXc07\n77zD2LFj7VWSXCS/sJy/rN3PgSN5NPc1MzOmH4PC2ju6LBEREfk/drv06e3tTVxcHHFxcdZtu3bt\nIjIyEoCuXbtSWFhISUkJhmHw8ssvM2PGDAIDA686dqtWfpiv811eQUEB1/U9V/P9j4X8+X8/J6+w\ngvZtmzN38s/pGOTv6LIcTvPDlvphS/2wpX7YUj9sqR8Nw25BzWw2YzbbDp+Xl0doaKj1c+vWrcnN\nzeWf//wnpaWlrFq1ioEDBzJy5Mgrjl1QUHZdNQUFBZCbW3xd33UlqUfz+N/3M6isstCtU0ueiuqD\nN4bb90bzw5b6YUv9sKV+2FI/bKkftm4ktDr0YQLDMACYMWOGI8twK9v3n2Tt1nNPdg7t34lHf9FV\nKw2IiIg4qUYNasHBweTl5Vk/5+TkEBQU1JgluK3aWoPE7UfYuvcEAA8MvoXJD/YlL6/EwZWJiIjI\n5dj1qc+LDR48mOTkZAAyMjIIDg7G31/3RdlbZZWFle+lsXXvCTw9TDxxX2/G3vkTLQclIiLi5Ox2\nRi09PZ1FixaRlZWF2WwmOTmZFStWEBoaSkxMDCaTiblz59rr8PJ/isqqeDXpIN9nF9Pc10xsVB96\ndm7l6LJERETkGtgtqIWFhbFmzZpLts+cOdNeh5SL5J4tZ2niAU4XlNO2pS/PjA+nfZvmji5LRERE\nrpFTrEwgDe/46WKWrTtIYWkVNwX7M2N8OC39fRxdloiIiNSDgpoL+vZ4Aa9tSKW80kKvzoHERvXF\nz1f/qUVERJoa/evtYvZ9m8Mbmw5RY6llYM8gnrw/FC9zoz4zIiIiIg1EQc2FfJqSRXzytxjAXQM6\n8mhkDzw89GSniIhIU6Wg5gIMw+D9L75n05c/APDgnV0YPegWvX5DRESkiVNQa+Jqaw3iP/6WHQd+\nxGSCiSN7MrRfR0eXJSIiIg1AQa0Jq66xsHrTIfZl5uJl9mDqA6EM6KGVHkRERFyFgloTVV5Zw2vr\nU/n2xFma+Zh5+uG+9Lgp0NFliYiISANSUGuCSsqrWbbuIN9nFxHo782M8f3oFKyluERERFyNgloT\nU1RaxZKEA5zMLaFtS19+90h/ggKbObosERERsQMFtSakoLiSJQkpZOeXEdLaj9/F9KN1C19HlyUi\nIiJ2oqDWROSdLeflhBRyz1bQKag5z8b0p2Vzb0eXJSIiInakoNYEnDpTxpKEFM4UVXJLuwBmRPfD\nv5mXo8sSERERO1NQc3Inc0tYknCAotIqunVqyW8fDte6nSIiIm5C/+I7sWOninkl8QAl5dX0vrkV\nv3moLz7eno4uS0RERBqJgpqTOnKykGVJByivtNC3axumjw3D20shTURExJ0oqDmhwz+c4bUNaVRW\nWxjYM4gpD4Ri9vRwdFkiIiLSyBTUnEz69/ms2JBGdU0tEaHt+NV9vfD0UEgTERFxRwpqTiTjhzPW\nkDa0XwceH9kTD5PJ0WWJiIiIgyioOYnDP5zhtfWpVNfUMqxfByYopImIiLg9XVNzAt8eL+DVDedC\n2pDw9gppIiIiAiioOVzmibMsT0qlqrqWwX3aMfGeXgppIiIiAiioOdS5V3AcpLLaQkRoO/5nVG+F\nNBEREbFSUHOQoz8WsnTdASqrLPz81hCeuK83Hh4KaSIiIvJfCmoO8H12EUsTD1BRZeH23sE8MVoh\nTURERC6loNbIjp0q5pWEcysODOwZxJP336r3pImIiEidlBAa0fHTxSxJSKGssoYBPc6tOKCQJiIi\nIpejlNBITuSUsCThAKUVNfTr1pZpY7QslIiIiFyZkkIjOHWmjCUJKZSUV9O3axv+39gwhTQRERG5\nKqUFOztTVMErCSkUl1UTeksrfv1gH7zMaruIiIhcnRKDHRWXVfFK4gHyiyrp2qEFv45SSBMREZFr\np9RgJ+WVNSxbd5Ds/DI6BjXn6XHh+HpraVURERG5dgpqdlBdY2HFhlR+OFVMUKAvz0b3w7+Zl6PL\nEhERkSZGQa2BWWpref39DL45fpaW/t48G9OfQH8fR5clIiIiTZCCWgOqNQz+9uE3pPwnj+a+Zp6N\n7kdwYDNHlyUiIiJNlIJaAzEMg4Rt/+HL9FN4e3nw9LhwOgX5O7osERERacIU1BrIB1/9wCd7T+Lp\nYSI2qg/dOrZ0dEkiIiLSxCmoNYBt+06ycef3mEww9YFQwrq0cXRJIiIi4gIU1G7QroxTrN2aCcAv\n7+nFwF7BDq5IREREXIWC2g1I/y6ft/51GIBxw7syJLyDgysSERERV6Kgdp2Ony7mrxvTqTUM7vlZ\nZ0b97GZHlyQiIiIuRkHtOpwpqmB50kEqqyzc3juYh4d1dXRJIiIi4oIU1OqprKKGZUkHOVtSRY+b\nAnnivlvxMJkcXZaIiIi4IAW1eqix1LJqYxpZuaW0b+NHrBZZFxERETtyipSRk5PD008/TVJSkqNL\nuSzDMPj7R99w6IcCWvh58dtx4Vq/U0REROzKrkEtMzOTyMhI4uPjrdsWLFhAdHQ0MTExpKamnivC\nw4Po6Gh7lnLD3v/ie5tVB4K0NJSIiIjYmd2CWllZGfPmzSMiIsK6bc+ePRw7dozExETmz5/P/Pnz\nAWjbti2enp72KuWG7Uz9kU1f/oDJBNMeCKNL+xaOLklERETcgNleA3t7exMXF0dcXJx1265du4iM\njASga9euFBYWUlJSgr9//dbEbNXKD7P5+oJdUFBAvfZP+TaHf2z5FoCpD/bl7kFdruu4zqq+/XB1\n6oct9cOW+mFL/bClfthSPxqG3YKa2WzGbLYdPi8vj9DQUOvn1q1bk5ubS1paGu+++y7FxcUEBgZy\n9913X3HsgoKy66opKCiA3Nzia97/RE4JC+P3Yak9966023u0rdf3nV19++Hq1A9b6oct9cOW+mFL\n/bClfti6kdBqt6B2LQzDACAiIsLmEqkzOP+utIoqC7f10rvSREREpPE16lOfwcHB5OXlWT/n5OQQ\nFBTUmCVck/LKGpYnpVJQXEn3Ti2ZPLq33pUmIiIija5Rg9rgwYNJTk4GICMjg+Dg4Hrfn2ZvtbUG\nb2zK4GRuCe1a+/HUQ33xus774URERERuhN0ufaanp7No0SKysrIwm80kJyezYsUKQkNDiYmJwWQy\nMXfuXHsd/rqt/+woqUfz8W/mxW/H611pIiIi4jh2C2phYWGsWbPmku0zZ8601yFv2Ffp2WzZfRxP\nDxPTx4YRrHeliYiIiAM5xcoEzuDoj4X87aNzr+F4NLI7vW5u5eCKRERExN0pqAEFxZWsfC+NGkst\nw/p3ZPiATo4uSURERERBrarawsr3UiksqaLnTYE8Gtnd0SWJiIiIAG4e1AzD4G9bvuH77GLatvRl\n+oNhmD3duiUiIiLiRNw6lXy0+zj/zjiNj5cnv3moLwF+3o4uSURERMTKbYPagSN5bNhxFIAn77+V\nTsHO9T43EREREbcMall5pazelIEBPHhnFwb0cL7VEUREREQcutZnYymvKSclJ42anEpqyswkf1xN\nRZWFgb2CGT3oFkeXJyIiIlInlw9qW37YRvKxT6myVFm3GV08adOmF0/cNxST1vAUERERJ+XSlz63\n/LCND75LtglpACZPC2WtMvg0a4djChMRERG5Bi4b1Mprykk+9ukV90k+9inlNRWNVJGIiIhI/bhs\nUEvJSbvkTNrFqixVpOSkNVJFIiIiIvXjskGtsLL4mvYrqiqycyUiIiIi18dlHyZo6RNwTfu18G5h\nl+MvWbKQ/fv3ApCVdZK2bYPw8fEB4M03/4GfX3Ob/efNe4HhwyO5444hlx3zrbfeIDc3h9mz/3jN\ndcyf/yK7dn1BixYtMQwDk8nE0KF3MXnytCt+79ChdN5883WWLl15xf0yMtLx8fGhWzctvSUiItLQ\nXDao9Q/uQ9J/Nl3x8qe3pzf9g/vY5fgzZ/7e+ueHH76fP/5xHuHh/S67/x//+Ge71HHu+DFMmjQZ\ngNLSEn77218TEhLCk0/+z2W/c+utYVcNaQAffriJvn37KaiJiIjYgcte+mxmbsbIm4dfcZ+RNw+n\nmdm3kSqyFRs7hTfe+CuPPfYwaWkHiY2dQnLyh2Rn/8iIEUN5551/8Pjj4xkz5h527txxyfdzck7z\n3HPPEBMTRUxMFLt2fXlNx23e3J9Ro0azZ89uAE6dOsWMGbE88kgUjz8+no8++hcA+/fvJTp6LHDu\nTN7SpYv4/e9nMm7cGJ58ciJ5eXls3LieLVs2s2rVayQkxDdMY0RERMTKZc+oAdxzyy8ALnmPmren\nNyNvHm79uaN8++03rFmzDg8P27xcVlYKmFizZh1ff72bP/3pD0RE3GGzz/z5LxIW1pfFi5dx8uQJ\npkyZxLvvbqBly8CrHrempgYvLy8AFi+eT//+P2Xp0pWcOpXNpEmP0K/fgEu+8+mn24iL+zshIe2Y\nNesZNm9+n1/+8gk++eRj7r9/LCNH3nv9jRAREZE6uXRQg3NhbWinwaTkpGHxqsSz2of+wX0cdibt\nQhERgy8JaeeNHj0GgNtu+xk1NTWcPHnC+rPy8nL279/LvHl/AaBTp5sID+/HV199wahRo694zIKC\nM2zevIlf/vIJqqur2bt3N3/+80IA2rVrT//+A9m372s6dOho873w8P60a9cegO7de3L69Knr+6VF\nRETkmrl8UANoZvZlUIfbCAoKIDf32p4GbQwtWtT9IIPJZLL5mb9/AMXF/306tbS0BMMwmDbtV9Zt\n5eXlDBhwW53jrV+fwMcffwSAr68v998/hrvuiuTs2bMYhoG//38XpA8ICKCgoOCSoHbhPh4eHtTW\n1tbjNxUREZHr4RZBrakxDIPCwrPWy5jFxUW0aNHS+vPAwFZ4enry5ptr8PPzu+p4Fz5McKFWrVrh\n4eFBUVGRNRgWFRXSunXrBvpNRERE5Ea47MMETd3WrVsA2LPn3/j4+HLTTZ2tPzObzUREDGbjxg0A\nVFRUsGDBn+p9OdJsNnP77T/n/fffA869RuTAgRQGDry9XmOUlDjPWUoRERFXojNqTsjT05Pq6hom\nTBhPcXEhs2Y9f8m9bDNn/p7Fixfwr39tBGDEiFGEhLSr97Fmzvw9ixbN56OPPsDLy4vZs58nJKQd\nWVknr+n7Q4YMZ9WqV/nxxyyeempGvY8vIiIil2cyDMNwdBH1db33mTnbPWp1yc7+kZiYB/nss912\nP1ZT6EdjUj9sqR+21A9b6oct9cOW+mErKOjaXsJfF136FBEREXFSCmoiIiIiTkpBzcm0b9+hUS57\nioiIiPNTUBMRERFxUgpqIiIiIk5KQU1ERETESSmoiYiIiDgpBTURERERJ6WgJiIiIuKkFNRERERE\nnJSCmoiIiIiTUlATERERcVIKaiIiIiJOSkFNRERExEkpqImIiIg4KQU1ERERESeloCYiIiLipBTU\nRERERJyUgpqIiIiIk1JQExEREXFSCmoiIiIiTsrs6AIAUlNTSUhIwDAMYmNj6dixo6NLEhEREXE4\nu55Ry8zMJDIykvj4eOu2BQsWEB0dTUxMDKmpqQC8++67vPjii0yfPp2kpCR7liQiIiLSZNjtjFpZ\nWRnz5s0jIiLCum3Pnj0cO3aMxMREjh49ypw5c0hMTKSmpgZvb2+CgoLIz8+3V0kiIiIiTYrdzqh5\ne3sTFxdHcHCwdduuXbuIjIwEoGvXrhQWFlJSUkKzZs2orKzk1KlTtG/f3l4liYiIiDQpdjujZjab\nMZtth8/LyyM0NNT6uXXr1uTm5hIdHc2LL76IxWJhxowZVx07KCjguuu6ke+6IvXDlvphS/2wpX7Y\nUj9sqR+21I+G4dCHCQzDACA0NJSFCxc6shQRERERp9Oor+cIDg4mLy/P+jknJ4egoKDGLEFERESk\nyWjUoDZ48GCSk5MByMjIIDg4GH9//8YsQURERKTJsNulz/T0dBYtWkRWVhZms5nk5GRWrFhBaGgo\nMTExmEwm5s6da6/Di4iIiDR5JuP8jWIiIiIi4lS0hJSIiIiIk3KKJaTsbcGCBRw8eBCTycScOXPo\n27evo0tqdLt37+bpp5+me/fuAPTo0YPJkyfz3HPPYbFYCAoK4uWXX8bb29vBldpXZmYm06dPZ9Kk\nSUyYMIHs7Ow6e7Bp0yb+/ve/4+Hhwfjx4xk3bpyjS7eLi/sxe/ZsMjIyCAwMBOCJJ55g2LBhbtOP\nxYsXs2/fPmpqapg6dSp9+vRx6/lxcT+2b9/utvOjvLyc2bNnk5+fT2VlJdOnT6dXr15uOz/q6kdy\ncrLbzo/zKioqGD16NNOnTyciIqJh5ofh4nbv3m1MmTLFMAzDOHLkiDF+/HgHV+QY//73v42nnnrK\nZtvs2bONDz/80DAMw3jllVeMtWvXOqK0RlNaWmpMmDDBeP755401a9YYhlF3D0pLS40RI0YYRUVF\nRnl5uXHfffcZBQUFjizdLurqx6xZs4zt27dfsp879GPXrl3G5MmTDcMwjDNnzhhDhw516/lRVz/c\neX5s3rzZWL16tWEYhnHy5EljxIgRbj0/6uqHO8+P85YuXWpERUUZGzZsaLD54fKXPi+3GoKcO8v2\ni1/8AoDhw4eza9cuB1dkX3WtllFXDw4ePEifPn0ICAjA19eXAQMGsH//fkeVbTd19aMu7tKP2267\njVdffRWAFi1aUF5e7tbzo65+WCyWS/Zzl37ce++9PPnkkwBkZ2cTEhLi1vOjrn7UxV36AXD06FGO\nHDnCsGHDgIb798Xlg1peXh6tWrWyfj6/GoI7OnLkCNOmTeORRx7hyy+/pLy83Hqps02bNi7fF7PZ\njK+vr822unqQl5dH69atrfu46pypqx8A8fHxTJw4kWeeeYYzZ864TT88PT3x8/MDYP369QwZMsSt\n50dd/fD09HTb+XFeTEwMM2fOZM6cOW49P867sB/gvn9/ACxatIjZs2dbPzfU/HCLe9QuZLjpQ663\n3HILsbGxjBo1ihMnTjBx4kSb/zt2175c6HI9cKfejBkzhsDAQHr37s3q1atZuXIl/fv3t9nH1fvx\nySefsH79et5++21GjBhh3e6u8+PCfqSnp7v9/EhISODw4cP87ne/s/ld3XV+XNiPOXPmuO382Lhx\nI/369eOmm26q8+c3Mj9c/oyaVkM4JyQkhHvvvReTyUTnzp1p27YthYWFVFRUAHD69OmrXgJzRX5+\nfpf0oK454y69iYiIoHfv3gDcddddZGZmulU/du7cyeuvv05cXBwBAQFuPz8u7oc7z4/09HSys7MB\n6N27NxaLhebNm7vt/KirHz169HDb+bFjxw62bdvG+PHjSUpKYtWqVQ3294fLBzWthnDOpk2beOut\ntwDIzc0lPz+fqKgoa28+/vhj7rzzTkeW6BCDBg26pAfh4eGkpaVRVFREaWkp+/fvZ+DAgQ6utHE8\n9dRTnDhxAjh3f0X37t3dph/FxcUsXryYN954w/rUmjvPj7r64c7zY+/evbz99tvAuVtqysrK3Hp+\n1NWPF154wW3nx/Lly9mwYQPr1q1j3LhxTJ8+vcHmh1u88HbJkiXs3bvXuhpCr169HF1SoyspKWHm\nzJkUFRVRXV1NbGwsvXv3ZtasWVRWVtKhQwcWLlyIl5eXo0u1m4tXywgJCWHJkiXMnj37kh5s2bKF\nt956C5PJxIQJE3jggQccXX6Dq6sfEyZMYPXq1TRr1gw/Pz8WLlxImzZt3KIfiYmJrFixgi5duli3\n/eUvf+H55593y/lRVz+ioqKIj493y/lRUVHBH/7wB7Kzs6moqCA2NpawsLA6/w511374+fnx8ssv\nu+X8uNCKFSvo2LEjd9xxR4PMD7cIaiIiIiJNkctf+hQRERFpqhTURERERJyUgpqIiIiIk1JQExER\nEXFSCmoiIiIiTkpBTUTEQbZt28bUqVNJTU11dCki4qTcbgkpEXGsxYsXk5aWRmVlJYcOHbIuMfPQ\nQw8xduxYB1dXt/fff58xY8Y0+Lh79uxh5cqVLF++nL59+1q3Hzp0iNdee42srCzrmpuxsbEMGjSo\nwWsQEeem96iJiEOcPHmSRx99lM8//9zRpVyRxWLh3nvvtb5hvCHt2LGDhIQE68tT4dwb35999lkW\nLFjA4MGDAcjMzOTo0aOMGjWqwWsQEeemoCYiDlFXUFuzZg0fffQRFouFn/zkJ8ydO5eDBw/y+uuv\n065dO9LS0ggPD6dnz55s3bqVs2fPEhcXR7t27di9ezfLly+nQ4cOZGVlERAQwLJly/D397/suKtW\nrcLHx4e7776b1NRUvvvuO6qqqggPD+f5558HYNasWWzevJnbb7+dqVOnsnz5ct59910AZs+ezU9/\n+lM6d+5sM1ZFRcUlx/P19b1qT6qrqxk5ciTPPfcc99xzj30aLyJNiu5RExGnkJqaytatW1m7di2J\niYkEBASQlJRk/dmsWbPYsGEDH3zwAS1atGDNmjWEhoayZcsW6xgZGRk899xzJCQkEBgYyHvvvXfF\ncdPT01m8eDGRkZH07NmTtWvXkpSUxBdffEFmZiZwbn3L1q1bW9c1vJzzY50PkXUd72p27twJwMiR\nI+vdPxFxTbpHTUScwu7duzl+/DgTJ04EoKysDLPZTI8ePejatat1YfDAwEDrfW0hISGUlJRYx+jW\nrRshISEADBgwgMOHD1NZWXnZcbt06UJgYCAWi4Xs7Gyio6Px9vYmNzeXgoKCetV/fqykpKQ6j3ct\nDh8+TGhoKCaTqV7HFhHXpaAmIk7B29ubu+66ixdeeMFm++7du/H09LTZduHnC+/euPjPJpPpiuN6\neXkBsHnzZtLS0li7di1ms5moqKg6a7w4QFVXV1v/fH6syx3vWjRv3pza2tp6f09EXJcufYqIUxgw\nYACff/45paWlAKxdu5aUlJR6jfHdd9+Rk5MDwL59++jZs+c1jZufn0+XLl0wm82kp6dz/Phxqqqq\nAPDw8KCmpgYAf39/Tp8+jWEYlJeXc/DgwQb9PYYOHcrXX39NWlqadds333zDZ599Vq8+iIjr0Bk1\nEXEKffr04bHHHuPxxx/Hx8eH4OBgoqKi6vWOsW7durF06VKOHTtGy5YtGTt2LH5+flcd95577mHa\ntGlMmDCBAQMG8Ktf/YqXXnqJdevWERwcTNu2bYmKiiI+Pp6ePXvy4IMP0rlzZ+sl2Gv5Pa5Fly5d\nWLZsGS+99BJlZWVUVVXRoUMHnn322WvugYi4Fj31KSIu4fxTn+efyBQRcQW69CkiIiLipHRGTURE\nRMRJ6YyaiIiIiJNSUBMRERFxUgpqIiIiIk5KQU1ERETESSmoiYiIiDgpBTURERERJ6WgJiIiIuKk\nFNREREREnJSCmoiIiIiT+v/90JufCa9jWQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f4958b41940>"
      ]
     },
     "metadata": {
      "tags": []
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "# Use Psat(T) to construct to show the vapor-liquid equilibrium diagram\n",
    "T = np.linspace(1,374)\n",
    "plt.figure(figsize = (10,6))\n",
    "plt.semilogy(T,[Psat(T) for T in T],linewidth=2)\n",
    "plt.xlabel('Temperature $^{\\circ}C$')\n",
    "plt.ylabel('Pressure [mmHg]')\n",
    "plt.title('Vapor Pressure of Water')\n",
    "\n",
    "# Additional annotations\n",
    "plt.semilogy(0.01,4.58,'o',markersize=10)\n",
    "plt.annotate('Triple Point', xy=(10,4.58), xytext=(20,3))\n",
    "plt.semilogy(100,760,'o',markersize=10)\n",
    "plt.annotate('Normal Boiling Point', xy=(100,760), xytext=(110,500))\n",
    "plt.semilogy(374,1.67e5,'o',markersize=10)\n",
    "plt.annotate('Critical Point', xy=(374,1.67e5), xytext=(300,2e5))\n",
    "plt.text(250,200,'Vapor',fontsize=16)\n",
    "plt.text(75,100000,'Liquid',fontsize=16)\n",
    "plt.ylim([1,1000000])\n",
    "plt.xlim([-10,400]);"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "-WcGP529p4q5"
   },
   "source": [
    "**Example:** The catapults on aircraft carriers require steam at 520 psig. What is the minimum operating temperature?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/",
     "height": 51
    },
    "colab_type": "code",
    "executionInfo": {
     "elapsed": 424,
     "status": "ok",
     "timestamp": 1539189411367,
     "user": {
      "displayName": "Jeffrey Kantor",
      "photoUrl": "",
      "userId": "09038942003589296665"
     },
     "user_tz": 240
    },
    "id": "yPN6phi8p4q5",
    "outputId": "cdb4c155-7537-4bbe-b75a-74607d58a709"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Operating Pressure = 27651.671203048452 mmHg\n",
      "Minimum Operating Temperature = [245.16090162] deg C\n"
     ]
    }
   ],
   "source": [
    "# convert pressure to absolute mmHg\n",
    "P = 520             # psig\n",
    "P = 520 + 14.696    # convert psig -> psia\n",
    "P = P*760/14.696    # convert psig -> mmHg\n",
    "\n",
    "# function to solve\n",
    "def f(T):\n",
    "  return Psat(T) - P\n",
    "\n",
    "# import a root-finding algorithm, provide initial guess\n",
    "from scipy.optimize import fsolve\n",
    "T = fsolve(f,200)\n",
    "\n",
    "print(\"Operating Pressure =\", P, \"mmHg\")\n",
    "print(\"Minimum Operating Temperature =\", T, \"deg C\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "97UEv7qSqyts"
   },
   "source": [
    "## Normal Boiling Points\n",
    "\n",
    "The **normal boiling point** of a pure component is the temperature at which the saturation vapor pressure is equal to one atmosphere.  In other words, it is the saturation temperature corresponding to a pressure of one atmosphere.\n",
    "\n",
    "Not every species will have a normal boiling point. For example, if the pressure at the triple point is above one atmosphere, then at one atmosphere  the species will sublime directly from the solid phase to vapor phase without boiling. Carbon dioxide is one example of such a compound.\n",
    "\n",
    "**Example:** Propylene glycol is a commonly used heat transfer fluid for industrial applications. What is the normal boiling point of propylene glycol?\n",
    "\n",
    "**Solution:** The Antoine coefficients are available on the [NIST Webbook entry for propylene glycol](https://webbook.nist.gov/cgi/cbook.cgi?ID=C57556&Mask=4&Type=ANTOINE&Plot=on). NIST, however, reports constants for the case where temperature is in Kelvin (K) and pressure in bar.\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/",
     "height": 34
    },
    "colab_type": "code",
    "executionInfo": {
     "elapsed": 395,
     "status": "ok",
     "timestamp": 1539195894280,
     "user": {
      "displayName": "Jeffrey Kantor",
      "photoUrl": "",
      "userId": "09038942003589296665"
     },
     "user_tz": 240
    },
    "id": "M6nIFp7t2LAy",
    "outputId": "0b94b481-95ea-4355-b120-0d728e49e034"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Normal boiling point of propylene glycol = 188.04733800941904 deg C\n"
     ]
    }
   ],
   "source": [
    "# propylene glycol (pressure in bar, temperature in K)\n",
    "\n",
    "A = 6.07936\n",
    "B = 2692.187\n",
    "C = -17.94\n",
    "\n",
    "from math import log10\n",
    "\n",
    "def Tsat(P):\n",
    "  P = 1.01325*P/760         # convert to pressure from mmHg\n",
    "  T = B/(A - log10(P)) - C  # compute temperature in K\n",
    "  T = T - 273.15            # convert temperature from K to C\n",
    "  return T\n",
    "\n",
    "print('Normal boiling point of propylene glycol =', Tsat(760), 'deg C')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "colab_type": "text",
    "id": "Ub9aY0GVMt2u"
   },
   "source": [
    "**Example:** What is the saturation pressure of water at the normal boiling point of propylene glycol. Why do you think it might be preferred to water?\n",
    "\n",
    "**Solution:** The advantage of propylene glycol as a heat transfer fluid is that it can be used for transfering heat in the range from 100C to 188C without the need for handling fluids under pressure.\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {
    "colab": {
     "base_uri": "https://localhost:8080/",
     "height": 34
    },
    "colab_type": "code",
    "executionInfo": {
     "elapsed": 250,
     "status": "ok",
     "timestamp": 1539195895469,
     "user": {
      "displayName": "Jeffrey Kantor",
      "photoUrl": "",
      "userId": "09038942003589296665"
     },
     "user_tz": 240
    },
    "id": "uRjufAbQMSqP",
    "outputId": "aeab44b2-7f27-4927-9191-00b320b5a3e6"
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Saturation pressure of water at 188.05C =  158.99201785868928 psig\n"
     ]
    }
   ],
   "source": [
    "# Antoine's equation for water from 1 to 374 degrees C\n",
    "\n",
    "def Psat(T):\n",
    "  if (1 <= T < 100):\n",
    "    return 10**(8.07131 - 1730.63/(T + 233.426))\n",
    "  elif (100 <= T <= 374):\n",
    "    return 10**(8.14019 - 1810.94/(T + 244.485))\n",
    "  else:\n",
    "    return float('nan')\n",
    "  \n",
    "P = 14.696*(Psat(188.05)/760 - 1)\n",
    "\n",
    "print('Saturation pressure of water at 188.05C = ', P, 'psig')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<!--NAVIGATION-->\n",
    "< [Gases with One Condensable Component](http://nbviewer.jupyter.org/github/jckantor/CBE20255/blob/master/notebooks/07.01-Gases-with-One-Condensable-Component.ipynb) | [Contents](toc.ipynb) | [Operating Limits for a Methanol Lighter](http://nbviewer.jupyter.org/github/jckantor/CBE20255/blob/master/notebooks/07.03-Operating-Limits-for-a-Methanol-Lighter.ipynb) ><p><a href=\"https://colab.research.google.com/github/jckantor/CBE20255/blob/master/notebooks/07.02-Vapor-Liquid-Equilibrium-for-Pure-Components.ipynb\"><img align=\"left\" src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open in Colab\" title=\"Open in Google Colaboratory\"></a>"
   ]
  }
 ],
 "metadata": {
  "colab": {
   "collapsed_sections": [],
   "name": "Vapor-Liquid_Equilibrium_for_Pure_Components.ipynb",
   "provenance": [],
   "toc_visible": true,
   "version": "0.3.2"
  },
  "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.7.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}