{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Chapter 17 - Adsorption"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Page 986 Example 17.1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      " Specific surface = 1.747*10**5 m**2/kg \n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAEPCAYAAACKplkeAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd4VGX2wPHvoReRpnQUlI6oFAERJYgFG2BZxRV/oliQ\nLl12EcSCILqiC+66iLoWVlZQcUU6offeERRpEnpvKef3x70xAySTmWRmbiY5n+fJkzvv3Ln35Bo8\nebuoKsYYY0wo5PI6AGOMMdmHJRVjjDEhY0nFGGNMyFhSMcYYEzKWVIwxxoSMJRVjjDEhE7akIiJj\nRSRORNb5lL0tIptEZI2ITBSRoj7vvSwiP4vIZhG5y6e8voisc98b6VOeX0S+dssXi8jV4fpZjDHG\nBCacNZVPgJYXlU0DaqvqDcBW4GUAEakFPAbUcj8zWkTE/cyHQAdVrQpUFZHka3YADrnlfwOGhfFn\nMcYYE4CwJRVVnQccuahsuqomuS+XABXc49bAOFWNV9UdwDagkYiUBYqo6lL3vH8DbdzjVsBn7vEE\noEVYfhBjjDEB87JP5RlgsntcDtjt895uoHwq5XvcctzvuwBUNQE4JiIlwhmwMcYY/zxJKiLyF+C8\nqn7lxf2NMcaER55I31BE2gP3cmFz1R6gos/rCjg1lD2kNJH5lid/5ipgr4jkAYqq6uFU7meLmxlj\nTAaoqqR/1oUiWlNxO9n7AK1V9azPW5OAtiKST0QqA1WBpaq6DzguIo3cjvsnge99PvOUe/wIMDOt\n+6qqfakyaNAgz2PIKl/2LOxZ2LPw/5VRYaupiMg4oBlwhYjsAgbhjPbKB0x3B3ctUtVOqrpRRMYD\nG4EEoJOm/FSdgE+BgsBkVZ3iln8MfC4iPwOHgLbh+lmMMcYEJmxJRVUfT6V4rJ/z3wTeTKV8BVAn\nlfJzwKOZidEYY0xo2Yz6HCQmJsbrELIMexYp7FmksGeReZKZtrNoICKa3X9GY4wJNRFBs3pHvTHG\nmOzNkooxxpiQSbejXkSKATcDlQAFduCM2joW1siMMcZEnTT7VETkVpw5JZWAVcBeQICyQF2c5DJc\nVedHItCMsj4VY4wJXkb7VPzVVB4Eeqnqz2ncsBrQEcjSScUYY0xwMvOHuI3+MsYY84f1+9fT9aeu\nxLaPDXlNBQARqQ3cxoV9KvNUdUOwNzPGGJM1HT17lMGxg/lq3VcMajaIWGIzdJ00R3+JyJMishQY\nAZQBfsFJKGWBESKyTETaZeiuxhhjsoQkTeKTVZ9Qc1RNTsefZkOnDXRu2DnD1/NXUykOtFDVE6m9\nKSKXA+0zfGdjjDGeWr53OV0mdwHgh8d/oEG5Bpm+pvWpGGNMDnPw9EEGzBzAD1t/YGiLofzfDf9H\nLrmw4SpsM+pFpKKIfCsiB9yvCSJSIb3PGWOMyVoSkhIYtXQUtUbVolDeQmzqvIn2N7a/JKFkRiCr\nFH8CfEnKisBPuGV3hiwKY4wxYTV/53y6TO5C8YLFmfXULK4rdV1Y7pNu85eIrFHVG9Iry6qs+csY\nk5PtPbGXvtP7Mue3Obxz1zv8qdafcPez8iucC0oeckeC5RaRPO6Ir4PB3sgYY0zknE88z9sL3ub6\nD6/nqqJXsanzJh6t/WhACSUzAmn+egb4AHjXfb0QeDpsERljjMmUadun0e2nblxT/BoWdVhE1ZJV\nI3Zvf2t/CdAGqAKs89nGN6pY85cxJqfYcXQHPaf2ZE3cGt67+z3ur3Z/hmsm4Wj+Gg30AEoAQ0Tk\nlQxFZowxJqzOxJ/h1dhXafBRA+qVrceGTht4oPoDYW/qSo2/5q/bgOtVNVFECuEsHDkkMmEZY4xJ\nj6oyacskekztQYNyDVj5wkquKnqVpzH5SyrnVTURQFVPixcpzxhjTKq2HNxC9ynd2XlsJ2MeGEOL\na1p4HRLgv0/lDLDNp+haYLt7rKp6fZhjCwnrUzHGZCcnz5/k9bmv8/Gqj3m56ct0bdiVvLnzhvw+\n4dhPpUbytTMWkjHGmFBRVcatH0ff6X1pcU0L1r24jjKXlfE6rEv4SyofAVOAn1R1c4TiMcYYc5G1\ncWvp+lNXTpw7wfg/jadJxSZeh5Qmf81fZYGWwN1AdWAJ8BMwQ1VPRSzCTLLmL2NMtDpy5giDYgfx\nn/X/YUjzITxX7zly58odkXtntPkroFWKRSQ30Ai4B7gdOAtMVdXhwd4w0iypGGOiTfIeJ3+Z9Rce\nrPEgr9/+OiULlYxoDGFLKiJyhaoevKjsSuAuVf0y2BtGmiUVY0w0WbpnKV0mdyFPrjz8/d6/U69s\nPU/iCHlSEZEHgLFAApAIPKaqCzIVpQcsqRhjosH+U/sZMHMAk3+ezFt3vEW769uFdEn6YIVjRv2b\nwK2qWhZ4GBgaZEBjRSRORNb5lJUQkekislVEpolIMZ/3XhaRn0Vks4jc5VNeX0TWue+N9CnPLyJf\nu+WLReTqYOIzxpisICEpgQ+WfEDt0bUpmr8omzpvSnXTrGjhL+qE5FFfqroEKBLktT/B6ej31R+Y\nrqrVgJnua0SkFvAYUMv9zGifyZYfAh1UtSpQVUSSr9kBOOSW/w0YFmR8xhjjqbm/zaXeP+vx/Zbv\nmdN+Du/c/Q5FCxT1OqxM8Tek+EoR6UnKPBXf16qq76b9UVDVeSJS6aLiVkAz9/gzIBYnsbQGxqlq\nPLBDRLYBjUTkN6CIqi51P/NvnEUup7jXGuSWTwD+7i8eY4zJKvYc30Of6X1YsGsB79z1Dg/XfNiT\ndbrCwV9NZQxO7eQy98v3dbC1lmSlVTXOPY4DSrvH5YDdPuftBsqnUr7HLcf9vgtAVROAYyJSIoNx\nGWNM2J1LOMew+cO44R83cG3xa9nYaSOP1Hok2yQU8FNTUdXB4byxqqqIWA+6MSZHmLJtCt2ndKda\nyWosfnYxVUpU8TqksEgzqYjIYOBDn5rFxe+XBTqq6qDU3k9DnIiUUdV97uf3u+V7gIo+51XAqaHs\ncY8vLk/+zFXAXhHJAxRV1cOp3XTw4MF/HMfExBATExNEyMYYk3G/HvmVl6a+xPr96xnZciT3VbvP\n65BSFRsbS2xsbKav429I8f1ALyAfsBL4Hac/pQxQDzgHjFDVyWle3OlT+UFV67ivh+N0rg8Tkf5A\nMVXt73bUfwU0xGnWmgFUcWszS4BuwFLgR+B9VZ0iIp2AOqr6ooi0BdqoattUYrAhxcaYiDsdf5ph\n84cxatkoet7ck54396RAngJehxWwcE5+rAjcglMrAPgNWKCqu9P+FIjIOJxO+Stw+k9eAb4HxrvX\n2gE8qqpH3fMH4GxdnAB0V9Wpbnl94FOgIDBZVbu55fmBz4G6wCGgraruSCUOSyrGmIhRVb7d/C09\np/akUYVGjLhzBBWLVkz/g1lMWJdp8blJWVX9PdibeMmSijEmUjYf3Ey3n7qx98RePrjnA5pXbu51\nSBkWjsmPqfkx2BsYY0x2d+LcCfpM68Otn9zKvVXvZdULq6I6oWRGsEkl+4x7M8aYTFJVvlj7BTVG\n1eDgmYOsf3E9PRr3CMumWdHC3+TH1PwrLFEYY0yUWb1vNV1/6sqZ+DNMeHQCjSs09jqkLCGoPpU/\nPiRymaqeDEM8IWd9KsaYUDp85jADZw3km03f8Frz1+hQt0PE9jiJpEj1qSTbmMHPGWNMVEpMSuSj\nFR9Rc1RNADZ13sTz9Z/PlgklM/xNfuzl53MZXabFGGOizuLdi+kyuQsF8hRgarup3FjmRq9DyrL8\n9am8AYwA4i8qFzJewzHGmKgRdzKO/jP7M237NIbdMYwn6jyRrdbpCgd/SWUV8J2qLr/4DRHpEL6Q\njDHGW/GJ8YxeNprX571O+xvas6nzJi7Pf7nXYUUFf0nlaZyZ6qm5KQyxGGOM52b/OpuuP3WlbJGy\nzG0/l5pX1vQ6pKiSodFf0cRGfxljArHr2C56T+/Nkt1LePfud3mwxoM5uqkr5KO/3O2A06yRiEgj\nEfkk2BsaY0xWcib+DEPmDOHGf95I9ZLV2dh5Iw/VfChHJ5TM8Nf89Tegj4g0BrZw4SrF1YGFOB35\nxhgTdVSVbzZ+Q5/pfWhYviErnl9BpWKVvA4r6gWySnF+nJWArwYUZ5XiNap6NvzhZZ41fxljLrZm\n3xq6T+nOkbNHGNlyJDGVYrwOKcvJaPOX32VaROR6VV0rIqdV9euMh2eMMd47ePogA2cNZOLmibwa\n8yrP1nuWPLmCXa3K+JPefJMOIlIVsCHExpioFZ8Yz8jFI6k5qib5cudjU+dNdGzQ0RJKGPibUT8I\nJ+ksAb4QkUGq+mrEIjPGmBCYvn06Pab2oFyRcsQ+FUvtUrW9Dilb89unIiKtgHuAn1R1UsSiCiHr\nUzEmZ9p2eBu9pvViw/4NvHv3uzxQ7QEb0RWEcC0o2UhVX8QmOxpjosSJcyfoP6M/jcc0pkmFJmzo\ntIFW1VtZQokQm/xojMkWkjSJz9d8zoBZA7jjmjsY2mIo5YqU8zqsqBWW0V/GGBMNluxeQrcp3QCY\n+OhEGlVo5HFEOZclFWNM1Np7Yi/9Z/Rn5q8zGdpiKO2ub0cusUXUveRvmZZHRaRAJIMxxphAnE04\ny9B5Q6nzYR3KFynP5s6b+b8b/s8SShbgr6byZ2CUiEwBxgFTVTUxMmEZY8ylVJXvt3xPr2m9qFOq\nDkufXcq1Ja71OizjI70hxUWBB4G2wI3Ad8A4VZ0TmfAyzzrqjckeNuzfQI+pPdh7Yi/v3f0ed157\np9chZWsZ7agPePSXiFwBPAx0BkqoaoVgb+YFSyrGRLfDZw4zaPYgvt7wNQNvG0jHBh3Jmzuv12Fl\ne+Gap5J88eLAQ8BjQAngv8HeyBhjgpGQlMCHyz6k5qiaJGoiGztvpGujrpZQsjh/y7QUIaXpqx4w\nCXgNiLU//Y0x4TT719l0n9KdEgVLMK3dNG4oc4PXIZkApdn8JSKHgORO+mmqej6SgYWKNX8ZEz12\nHN1B72m9Wb53OSPuGsHDNR+2mfAeCcfkxwqqeiYTMRljTEBOnT/FW/PfYvTy0fRo1IPPH/ycgnkL\neh2WyQB/fSrfiMifRKTQxW+ISGEReUxEJmfkpiLysohsEJF1IvKViOQXkRIiMl1EtorINBEpdtH5\nP4vIZhG5y6e8vnuNn0VkZEZiMcZ4R1X5at1X1BhVg+1HtrP6hdUMbDbQEkoU89f8VQroAjwCJHLh\ndsJ5gK+BUap6IKgbilQCZgE1VfWciHwNTAZqAwdVdbiI9AOKq2p/EakFfIWzqGV5YAZQVVVVRJYC\nXVR1qZvg3lfVKRfdz5q/jMmCVuxdQbcp3TibcJaRLUfS9KqmXodkfIS8+UtV9wOvAK+ISBmc7YQB\nflPVfRkLE4DjQDxQSEQSgULAXuBloJl7zmdALNAfaI0zNyYe2CEi24BGIvIbUERVl7qf+TfQBqcf\nyBiTRcWdjGPAzAFM3jaZ15u/Tvsb25M7V26vwzIh4m+Zlqki8pKI1FDVfaq6xP3KTEJBVQ8D7wA7\ncZLJUVWdDpRW1Tj3tDigtHtcDtjtc4ndODWWi8v3uOXGmCzofOJ5RiwcQe3RtSlesDibO2+mQ70O\nllCyGX8d9e2BlsBgEamOswPkT8AMVT2V0RuKyLVAD6AScAz4r4i08z3HbdoKWZvV4MGD/ziOiYkh\nJiYmVJc2xgTgx60/8tLUl6hasioLnllA9Suqex2SuUhsbCyxsbGZvk5AM+pFJDfQCGcXyNuBszhr\ngQ0P+oYijwF3quqz7usngcbudZur6j4RKQvMVtUaItIfQFXfcs+fAgwCfnPPqemWPw40U9WOF93P\n+lSM8cjmg5vpObUn249s5293/417q97rdUgmQGGdUa+qiaq6UFUHquotOBMi9wR7M9dmoLGIFBRn\nAPodwEbgB+Ap95yncNYZA2fSZVsRyScilYGqwFK3Ge64iDRyr/Okz2eMMR46evYovab24tZPbuWO\na+5g3YvrLKHkEH73UxGRvMBdwG04zVUK7ADm4Yz+CpqqrhGRfwPLgSRgJfARUAQYLyId3Hs86p6/\nUUTG4ySeBKCTT9WjE/ApUBCYfPHIL2NMZCUmJTJ21VgGzh7IA9UeYP2L6yl9Wen0P2iyDX9Digfi\nLCC5CFhKypDiskBDnCarb1T19ciEmjHW/GVMZMz7bR7dp3SnUN5CjGw5kvrl6nsdksmEkK9SLCKt\ngB/S+j+yiOQC7lfVScHeNJIsqRgTXruO7aLvjL7M3zmf4XcMp+11bW1plWwgEkvfXwagqieDvYmX\nLKkYEx5n4s/w9sK3GblkJJ1v6ky/W/pROF9hr8MyIRKOtb+SL1wHZ2JhSff1AeApVV0fdJTGmKin\nqnyz8Rv6TO9Dw/INWfH8CioVq+R1WCaLSDep4HSi91TV2QAiEuOWNQljXMaYLGjNvjV0n9KdI2eP\n8GmbT4mpFON1SCaLCSSpFEpOKACqGisiVsc1Jgc5cOoAA2cP5NvN3/JqzKs8W+9Z8uQK5H8fJqcJ\nZJ7KryIyUEQqiUhlEfkr8Eu4AzPGeC8+MZ6Ri0dSa3Qt8ufOz6bOm+jYoKMlFJOmQH4zngaGABPd\n1/OAZ8IWkTEmS5i2fRo9pvSg/OXliX0qltqlansdkokC/rYTLgh0BKoAa3H6VeIjFZgxxhvbDm+j\n17RerN+/nnfvepdW1VvZEGETMH/NX58B9YF1OGt+jYhIRMYYT5w4d4J+0/vRaEwjbq5wMxs7baR1\njdaWUExQ/DV/1VTVOgAiMgZYFpmQjDGRlKRJ/HvNvxkwcwB3Xnsn615cR7ki5bwOy0Qpf0klIflA\nVRPsrxVjsp9FuxbRfUp3RIRvH/uWRhUaeR2SiXL+lmlJBE77FBUEzrjHqqqXhzm2kLAZ9cZcasfR\nHfSf0Z/5O+fzZos3aXd9O3JJQIuWmxwiHNsJ23ZsxmQzx88dZ+i8oXy08iO6NezGx60+tqVVTEj5\nG/21ApiPs9tjrKqejVhUxpiQSkhK4OOVHzModhAtq7Rkbce1lL/cdt82oeev+Ssv0BRnS+EY4DAw\nBfhJVbdGKsDMsuYvk9NN2z6NXtN6UaJgCd69611bkt4EJBKrFJfHSTB348xdWayqnYK9YaRZUjE5\n1cYDG+k9rTdbD23l7Tvfpk2NNjY82AQsLElFRK4Erga2qepRn/LcQGNVXZCRYCPJkorJaQ6cOsDg\n2MGM3zieAU0H0LlhZ/Llzud1WCbKhHyPehF5FtgAfABsEZHWye+5e9Zn+YRiTE5yLuEcIxaOoNbo\nWuTOlZvNnTfz0s0vWUIxEeVvnspLQG1VPSAi1wBfAd9HJixjTKBUlQmbJtBvRj9qX1mbeU/Po8YV\nNbwOy+RQ/pLKeVU9AKCqv4hI/gjFZIwJ0LI9y+g5rSfHzx3no/s/osU1LbwOyeRw/pJKBRF5H0hu\nUyvv81pVtVvYozPGpGrXsV0MmDWAmb/M5LXmr9H+xvbkzmVTy4z3/CWVPoBvD/cK97VcVG6MiZCT\n508yfMFwRi0bxYsNXmRLly0UyV/E67CM+YO/GfWfRjAOY4wfiUmJfLbmMwbOHkjzSs1Z9cIqrip6\nlddhGXMJfzPqxwIfqmqqqxOLSCOgo6o+Ha7gjDEw69dZ9JrWi0J5CzHx0Ym26KPJ0vzNqK+D0wTW\nGNgC/I7T9FUGqA4sBEao6vrIhJoxNk/FRKuth7bSZ3of1sWtY9gdw3ik1iM2edFETNhm1Lujvuri\nTIJU4DdgTbSsBWZJxUSbw2cOM2TOEL5Y+wV9b+lLt0bdKJCngNdhmRwm5KsUJ1PVc8Bi98sYEybn\nE88zetlo3pz3Jo/UeoSNnTdSqnApr8MyJijpJhVjTHipKpO2TKLP9D5cW+JaZj81m9qlansdljEZ\n4klSEZFiwBigNk6T2tPAz8DXOM1sO4BHk9cbE5GXgWeARKCbqk5zy+sDnwIFgMmq2j2iP4gxmbTq\n91X0mtaL/af28/4979OySkuvQzImU/xu9SYiuUVkRBjuOxInCdQErgc2A/2B6apaDZjpvkZEagGP\nAbVwVkkeLSm9lR8CHVS1KlBVROxfpIkKe0/s5Znvn+GeL+/hsdqPsbrjaksoJlvwm1RUNRFoKiEc\nciIiRYFbVXWse48EVT0GtAI+c0/7DGjjHrcGxqlqvKruALYBjUSkLFBEVZe65/3b5zPGZEmn408z\nZM4Q6nxYh1KFS7GlyxZeaPACeXJZS7TJHgL5TV4NfC8i/yVlz3pV1YkZvGdl4ICIfALcgDNTvwdQ\nWlXj3HPigNLucTkuHCSwGygPxLvHyfa45cZkOUmaxJdrv2TArAE0qdiE5c8tp3Lxyl6HZUzIBZJU\nCuDs+nj7ReUZTSp5gHpAF1VdJiLv4TZ1JVNVFZGQjQMePHjwH8cxMTHExMSE6tLGpGveb/PoOa0n\nuSQXXz/yNU0qNvE6JGMuERsbS2xsbKavE/DOj6EiImWARapa2X3dFHgZuAZorqr73Kat2apaQ0T6\nA6jqW+75U4BBOPNlZrv9MojI40AzVe140f1snorxxPbD2+k3ox/L9i5jaIuhtL2uLbnEb4uzMVlG\nyDfp8rlwRRH5VkQOuF8TRKRCxsIEVd0H7BKRam7RHTibgf0APOWWPQV85x5PAtqKSD4RqQxUBZa6\n1zkuIo3cPp8nfT5jjGeOnj1Kn2l9aDSmEfXK1mNz5838uc6fLaGYHCGQ5q9PgC+BR93XT7hld2bi\nvl2BL0UkH7AdZ0hxbmC8iHTAHVIMoKobRWQ8sBFIADr5VD064QwpLogzmmxKJmIyJlPiE+P5aMVH\nDJk7hFbVWrG+03rKXFbG67CMiahAlmlZo6o3pFeWVVnzlwk3VWXyz5PpPb035YuU55273uGGMlHx\nz8OYNIVtmRbgkIg8ibOdsABtgYPB3siY7Ghd3Dp6TevFzmM7eeeud7i36r226KPJ0QJp5H0Gpylq\nH85KxX/Caa4yJseKOxnHCz+8wB2f30Gr6q1Y9+I67qt2nyUUk+MFsqDkDuCB8IdiTNZ3NuEsf1v0\nN95Z9A5P3fAUmztvpnjB4l6HZUyWYdN4jQmAqvL1hq/pP6M/9cvVZ/Gzi6lSoorXYRmT5VhSMSYd\ni3Ytoue0npxPPM9nbT6jWaVmXodkTJZlScWYNOw4uoP+M/ozf+d83mzxJu2ub2dzTYxJR7pJRURq\nA7cBlXCWqd8BzFPVDWGNzBiPHD93nKHzhvLRyo/o1rAbH7f6mML5CnsdljFRIc2k4g4j7gocApYC\nv+AMKS4LjBCRK4CRqvpFJAI1JtwSkhIYu2osg2IH0bJKS9Z2XEv5y22NUmOC4a+mUhxooaonUntT\nRC4H2ocjKGMibfr26fSc1pOSBUvy459/pF7Zel6HZExUiviCkpFmM+qNPxsPbKTP9D5sPbSVt+98\nm9bVW9tcE2MI04x6EckL3MWFfSq/AXOBqaqaEHyoxnjvwKkDDI4dzPiN4xnQdADfPvYt+XLn8zos\nY6JemkNZRGQgsAy4H2e737E4OzJuwZkMuVxE/hqJII0JlWNnj/HG3DeoNboWeXLlYXPnzbx080uW\nUIwJEX81lTXA62m0HY0VkVw4CceYLO/o2aO8v+R9Plj6AfdUuYcFzyygWslq6X/QGBMU61Mx2drh\nM4d5b/F7jF42mgeqP8CApgOoWrKq12EZk+WFfJMuEekiIle6x1VEZK6IHBWRJSJSJzPBGhNuB08f\n5C8z/0LVD6qy98Relj63lE9af2IJxZgw8zc9uJOqHnCP3wf+hjPMuB/wj3AHZkxG7D+1n37T+1H9\n79U5dOYQK55fwZhWY7im+DVeh2ZMjuCvTyW3z/GVqvqtexwrIkXCGJMxQdt3ch8jFo5g7KqxPH7d\n46x+YTUVi1b0Oixjchx/NZVvRORTEbkG+FZEXhKRq0XkaWBnhOIzxq+9J/bSY0oPao2qRXxiPOte\nXMeo+0ZZQjHGI3476t0E0hG4FsgP7Aa+A95S1WMRiTCTrKM+e9p9fDfD5g/jy3Vf0v7G9vRp0oey\nRcp6HZYx2UZGO+pt9JeJKjuP7WTovKF8veFrOtTtQO8mvSl9WWmvwzIm2wnnHvW+N/lIVZ8P9ibG\nZNavR35l6PyhTNg0gefrPc+WLlu4svCVXodlolRiIhw8CPv2QVxcyveTJ1M/39/fpRl5L6tfLzOC\n3U/lptCHYEzath/ezpvz3uS7Ld/xYoMX2dplKyULlfQ6LJMFJSXBoUOXJork777Hhw5BsWJQpgyU\nLp3yvUgRSGvpN39LwmXkvax+vYwKqvlLRKaoasvQhhBe1vwVnbYe2sqb897kf1v/R5eGXejeqLvt\nBZ8DqcLhw4EligMH4PLLL00UvsfJ36+8EvLm9fqny9oi0vwVbQnFRJ9NBzbxxrw3mLp9Kt0admNb\nt20UK1DM67BMCKnC0aOpJ4aLv+/fD4ULp54Yqle/8HWpUpDPlnDzXJo1FREZC3yoqsvSeL8R0FFV\nnw5jfJlmNZXosGH/Bl6f9zozf5nJS41fonPDzlye/3KvwzIBUoXjxwNLFHFxUKBA6oni4u+lSjnn\nmsgL+egvdymWPkBjnJWJf8fZ+bEMUB1YCIxQ1fUZDToSLKlkbWvj1vLa3NeY99s8et7ckxcbvEiR\n/Da3NitQdTqu/SUI3+958waWKEqXhoIFvf7pTHrCNqRYRPIDdYGrSdlPZY2qns1IoJFmSSVrWvX7\nKl6b+xqLdi+i98296digo+0DH2GHD8OqVbB1a9o1ClUnGQTST1HY/vNlK2FJKiJyvaquTf6eqQg9\nYkkla1m+dzlD5gxhxe8r6NukL8/Vf45CeQt5HVa2t38/rFzpfK1Y4Xw/eBDq1nX6JsqWTT1RXHZZ\n6EcHmegQrqQyEvg70EVVu2civtSunRtYDuxW1QdEpATwNU6NaAfwqKoedc99GXgGSAS6qeo0t7w+\n8ClQAJicWoyWVLKGJbuXMGTuENbGraXfLf14tt6zFMhjjeWhpgq//56SOJK/TpyAevWgfn3ne716\nULUq5PK3UJPJ0cLRpzIIuAJ4AvgCOKSqr2Yqyguv3xOoDxRR1VYiMhw4qKrDRaQfUFxV+4tILeAr\nnDky5YE+pgVdAAAWw0lEQVQZQFVVVRFZipPwlorIZOB9VZ1y0X0sqXho4a6FvDrnVTYf3MzLTV/m\n6RufJn+e/F6HlS2ows6dF9Y+Vq6EhISU5JH8vXJlq3GY4IR8SLGqvioirdxzZqjqpMwE6EtEKgD3\nAm8APd3iVkAz9/gzIBboD7QGxqlqPLBDRLYBjUTkN5yEtNT9zL+BNsAFScV4Y+5vcxkyZwjbj2xn\nQNMBPHXjU7ZlbyaowvbtlzZh5cuXkjief975XrGiJRDjnfTmqTRS1RdF5DUgZEkFZ2+WPoDvmNHS\nqhrnHscByQs6lQMW+5y3G6fGEu8eJ9vjlhuPqCqxO2J5dc6r7D6+m7/c+hfaXd+OvLltllkwEhPh\n558vrH2sWuVM7EuuffTo4RyXtTU0TRaTXlIZ737/b6huKCL3A/tVdZWIxKR2jtu0FbI2q8GDB/9x\nHBMTQ0xMqrc1GaSqzPhlBkPmDiHuZBx/ve2v/LnOn8mTK9hVgHKehATYtOnCGsiaNc78jOS+j/79\nne9X2lJnJoxiY2OJjY3N9HUi3lEvIm8CTwIJOB3slwMTcfpMYlR1n4iUBWarag0R6Q+gqm+5n58C\nDMIZ2jxbVWu65Y8DzVS140X3sz6VMFFVpm6fypA5Qzhy9gh/vfWvtL2uLblz5U7/wznQ+fOwfv2F\nHejr1kGFChd2oNetC8VtRRrjsajrqHfv0Qzo7Y7+Gu7eY5ibSIpd1FHfkJSO+ipubWYJ0A1YCvyI\nddRHhKry488/MmTOEE7Hn2bgbQN5pNYjlkx8nDnjJAzf/o9Nm+Caay7sQL/hBqdZy5isJqo66i++\nlfv9LWC8iHTAHVLsxrFRRMYDG3FqN518skQnnCHFBXGGFFsnfRipKpO2TGLI3CHEJ8bzSrNXeKjm\nQ+SSnD0u9dQpWL36wiasbduc+R/JtY9nnnESSCGbkmOyufSav95Q1b+IyGuqOjCCcYWM1VQyL0mT\n+G7zdwyZMwQR4ZXbXqF1jdY5MpkcO+YkEN9O9B07oHbtC5uwrrvO1qwy0c12fkyDJZWMS9IkJmyc\nwGtzXyN/nvy8ctsr3F/tfiSHjFc9dMgZdeXbhPX773D99Rc2YdWqZcuom+wnHH0qxXDmibTBGd6r\nwH5S9qg/mvFwI8eSSvASkxIZv2E8r897nSL5ivBKs1e4p8o92TqZ7N9/6Sz0Q4ecTvPk2kf9+k6T\nVm7rOjI5QDiSyjRgJs5ExDi3Y7ws8BRwu6relZmAI8WSSuASkhL4z/r/8Prc1ylZqCSDmg3izmvu\nzHbJ5MABWLz4wiRy6tSly5hUqWLLmJicKxxJZauqVgv2vazGkkr64hPj+XLdl7wx7w3KXlaWQc0G\ncXvl27NNMjlyBObMgdmzYdYsZ2mTxo0vXMqkUiWbhW6Mr3Ds/PibiPQFPkue6S4iZXBqKjszFqbJ\nSs4nnufzNZ/z5vw3ubro1Yx5YAzNKjVL/4NZ3PHjMHeuk0Rmz3ZGYt18MzRvDmPGOEkkj83LNCYs\n/NVUSuD0qbQiZcmUOJzlWt5S1cMRiTCTrKZyqXMJ5/h09acMnT+UaiWrMfC2gdx69a1eh5Vhp07B\n/PkpSWTjRmjY0EkizZvDTTfZNrPGBMtGf6XBkkqKswlnGbtqLG/Nf4vapWrzym2vcHPFm70OK2hn\nzsDChSlJZM0apxkrOYk0bmzDeY3JrIgmFRF5WlU/CfqDHrCkAifOnWDMyjG8s+gd6paty8DbBtKw\nfEOvwwrYuXOwZElKElm+3BnWm5xEmjSxSYXGhFqkk8ouVa0Y9Ac9kJOTStzJOD5Y+gH/XPFPbq98\nO32b9KV+ufpeh5Wu+Hgnccya5SSRJUugRo2UJNK0KRSxbeyNCauQd9SLyDo/nysV7I1M5Gw/vJ0R\nC0fw9YavaXtdWxZ1WESVElW8DitNiYnOsN7kmsiCBc4aWbffDt27w623QrFiXkdpjAmEvzEwpYCW\nwJFU3lsYnnBMZqzYu4JhC4Yx69dZdGzQkc1dNlOqcNbL/0lJsHZtyhDfefOclXqbN4fnnoMvvoCS\nJb2O0hiTEf6Syo/AZaq66uI3RGRO+EIywVBVpv8ynWELhrH10FZ6Nu7Jx60+pkj+rNM+pAobNqTU\nRObMgSuucJLIk0/Cxx87+4cYY6Kfjf6KUglJCfx3w38ZvnA48Ynx9L2lL49f93iW2GVRFbZuTUki\ns2c7fSDJfSIxMVDe9ug0JksLS0e9iDyhql+KyOOqOi5TEXokuyWV0/GnGbtqLO8seoeril5F3yZ9\nubfqvZ7OfleFX365MInkyZOSRJo3h6uv9iw8Y0wGhGNGPUA5EXkUqJCxsEyoHDx9kFFLRzFq2Sia\nXtWUcQ+Po3GFxp7Fs3PnhUkkPj4lgQwZ4nS027InxuQ86e38WADoA7wNnA31zo+REO01lR1Hd/Du\nonf5Yu0XPFzzYXo36U31K6pHPI69ey9MIidOOM1YyYmkenVLIsZkJ+Fq/uoN7AHKq+qITMTnmWhN\nKmv2rWH4wuFM2TaFZ+s+S/fG3SlXpFzE7r9/P8TGpiSRAwegWbOUJFK7tiURY7KzcDV//a6q40Tk\n8QzGZYKgqsTuiGXYgmGsjVtLj8Y9GH3vaIoWKBr2ex8+fOFKvrt3O/NDmjeH5593tsK1ZeCNMemx\n0V9ZQGJSIt9u/pbhC4Zz/Nxx+jTpQ7vr25E/T/6w3fPYsQtX8t2+HW65JaUmUreureRrTE5mC0qm\nISsnlbMJZ/ls9WeMWDSCKwpdQb9b+tGqequw7P1++vSFSWTTJmjU6MKVfG1LXGNMMksqaciKSeXI\nmSN8uPxDPlj6AQ3KNaBvk740vappyIcFHz8OP/4IEyfCtGlOE1aLFk4SadQI8oevImSMiXLhWPtr\nNPCVqs7PVGTmD7uO7eK9xe/xyepPaFW9FdOfnM51pa4L6T0OHYJJk2DCBKdmcttt8NBD8OGHzix2\nY4wJJ3+t5luBt0WkHPA1MC61JVtM+jYe2MjwBcOZtGUS7W9sz5qOa6hYNHSLPP/+O3z7rVMjWbYM\n7rwTnngCvvwSioa/j98YY/6QbvOXiFQC2gKPAYWAr3ASzNZwBxcKXjZ/zd85n2ELhrFszzK6NuxK\np5s6Ubxg8ZBc+9dfnUQyYYLTP3LffU6N5O67bW8RY0zmRaRPRUTqAp8AdVQ1d7A380Kkk0qSJvHD\nlh8YvnA4cSfj6N2kN0/d8BQF8xbM9LU3bXJqIxMmOEN+W7eGhx92loi37XKNMaEUtqQiInmAe3Fq\nKy2A2Tg1le8zEmikRSqpnEs4x5frvuTthW9TOG9h+t3Sj4dqPkTuXBnPvaqwerWTRCZOdDreH3rI\n+Wra1Ib8GmPCJ+RJRUTuwkkk9wFLgXHAJFU9mZlAIy3cSeX4ueP8c/k/eW/Je9QpVYe+t/SleaXm\nGR7JlZQEixc7SWTiRGfW+sMPO1833WQTEI0xkRGOpDIbp/9kgqoezmR8nglXUvn9xO+MXDKSMSvH\ncHeVu+nTpA83lrkxQ9dKSHBGak2Y4PSTlCjhJJGHHnL2YrflUIwxkRaOZVpaq+rxdG5aRFVPBHND\nEakI/BtnZ0kFPlLV90WkBM4os6uBHcCjqnrU/czLwDNAItBNVae55fWBT3EWvpysqt2DiSUjth7a\nytsL3uabTd/Qrk47lj23jMrFKwd9nXPnYMYMpzYyaRJUquQkktmzncUZjTEmGvmrqcwAtgDfA8uT\naysiUhJoALQBqqrqHUHdUKQMUEZVV4vIZcAK91pPAwdVdbiI9AOKq2p/EamFU2O6CSgPzHDvqyKy\nFOiiqktFZDLwvqpOueh+IampLNm9hOELhzPvt3l0uqkTXRp24YpCwU38OHUKpkxxaiQ//QTXXeck\nkgcftP1GjDFZS7hWKb4d+DNwC5C8RO5eYD7wparGBh/qJff4Dvi7+9VMVePcxBOrqjXcWkqSqg5z\nz58CDAZ+A2apak23vC0Qo6odL7p+hpOKqvLTtp8YvmA4O47uoNfNvXim7jMUzlc44GscPQr/+59T\nI5k505nJ/vDDzsitMmUyFJYxxoRdOGbU51XVWcCsTEXmhzsHpi6wBCitqnHuW3FAafe4HLDY52O7\ncWos8e5xsj1ueabFJ8bzn/X/YfjC4eSW3PS9pS9/qvWngLfqPXAAvvvOSSQLFjjLojz0EIwZ4/SX\nGGNMduWvT2WRiOwBfgKmqOqOUN7YbfqaAHRX1RO+o6Xcpq2Iz1g8ef4kY1aO4d1F71KlRBVG3DmC\nu669K6CRXLt3p0xGXL3amYT49NMwfryzP7sxxuQEaSYVVW0gIpWBlsB7IlIBp9lrMjBHVc9l9KYi\nkhcnoXyuqt+5xXEiUkZV94lIWWC/W74H8F3TpAJODWUPF25zXMEtu8TgwYP/OI6JiSEmJuaC9/ef\n2s8HSz7gHyv+QUylGCY8OoGbyt+U7s+xfXvKHJKff4YHHoCePZ1lUgpmfq6jMcZETGxsLLGxsZm+\nTsAz6kUkH3ArTpJpBhxQ1fuCvqHzZ/9nwCFVfcmnfLhbNkxE+gPFLuqob0hKR30VtzazBOiGM4/m\nR4LsqP/lyC+MWDiCcevH8Vjtx+jdpDdVSlRJM3ZV2LAhZVZ7XBy0aeP0kcTE2NLxxpjsI+JL34tI\nBVXdnf6Zl3yuKTAXWIszpBjgZZzEMB64ikuHFA/AGVKcgNNcNtUtTx5SXBBnSHG3VO53SVJZ+ftK\nhi8YzoxfZvBC/Rfo2qgrZS5LvddcFVasSKmRnD2bMqu9SRPIHRWL1RhjTHDCNfqrHvA4cBtQCScJ\n/IaTFL6KhlWLk5OKqjLjlxkMXziczQc381Ljl3iu3nMUyX9ph0diIixcmDKrPX/+lFnt9evbZERj\nTPYXjhn1k4EjwCScWsTvgABlcZqiHsBpogq6CSySRETHrRvH8AXDOZd4jr5N+vJ4ncfJl/vCFRjj\n4yE21qmRfPcdlC6dMqu9dm1LJMaYnCUcScV3iG9a55RS1f3+zvGaiGjTsU3pd0s/7q167wVb9Z45\nA9OnO7WRH36AqlVTmraqpN21Yowx2V5Y+1TcyYg34TR/Lc3qicTXxX0qJ07A5MlOIpk6FerWdZLI\ngw9ChQp+LmSMMTlIOJe+fxR4G5jjFt0G9FHV/wYdpQdERA8dUn74wWnaio2FW25xmrZatYJSpbyO\n0Bhjsp5wJpW1wB3JtRMRuRKYqarXZyjSCBMRvfxypUULp0Zy//1QrJjXURljTNYWjlWK/7g2cMDn\n9SG3LGrs3QuFA1+uyxhjTAYFklSmAFNF5CucZPIYztItUcMSijHGREagHfUPAU3dl/NU9duwRhVC\nkd6j3hhjsoNwrFJcDaeDvgrO7Pc+GZlBb4wxJufwt+P5WOB/wMPASuD9iERkjDEmavnrU7lMVf/l\nHm8WkSy/JIsxxhhv+UsqBdy1v8DpoC/ovhacLU9Whj06Y4wxUcXfMi2xpKwiDG4ySX6hqs3DGlmI\nWEe9McYELxxrf5VX1VQ3vYomllSMMSZ44VqluCQwG2euynxVTchUlB6wpGKMMcEL134qBYEY4B6g\nCbCLlD3rd2Ys1MiypGKMMcGLyM6PInINToK5Gyijqg2DvWGkWVIxxpjghaP5S4A2uJMfk7fw9Xk/\nv6qey0iwkWRJxRhjgheOpPIhUAtYCLQA/qeqQzIVpQcsqRhjTPDCkVQ2ANeraqKIFMLpqK+X6slZ\nmCUVY4wJXkaTir9lWs6raiKAqp4mypa7N8YYE3n+aipngG0+RdcC291jjaZNuqymYowxwQnHJl01\nMxGPMcaYHMjv6K/0/sQP5ByvRUGIxhiT5YSjTyVWRPq4+6pcfLPqItIPmBPsDY0xxmRf/moq+YEn\ngMeB64ATOJ31lwHrgS+Br1T1fGRCzRirqRhjTPDCOqNeRHIDV7gvDyaPCosGllSMMSZ4EVmmJRpZ\nUjHGmOCFo08lKohISxHZLCI/u/08xhhjPBLVScVtlvs70BJnSZnHRcSGQqchNjbW6xCyDHsWKexZ\npLBnkXlRnVSAhsA2Vd2hqvHAf4DWHseUZdk/mBT2LFLYs0hhzyLzoj2plMfZ4yXZbrfMGGOMB6I9\nqVgPvDHGZCFRPfpLRBoDg1W1pfv6ZSBJVYf5nBO9P6Axxngoxw0pFpE8wBac/V72AkuBx1V1k6eB\nGWNMDuVvQcksT1UTRKQLMBXIDXxsCcUYY7wT1TUVY4wxWUu0d9T/Ib1JkCJSQ0QWichZEenlRYyR\nEsCzeEJE1ojIWhFZICJRsTdORgTwLFq7z2KViKwQkdu9iDMSAp0oLCI3iUiCiDwUyfgiKYDfixgR\nOeb+XqwSkb96EWckBPJ74T6PVSKyXkRi/V5QVaP+C6fpaxtQCcgLrAZqXnTOlUAD4HWgl9cxe/ws\nbgaKusctgcVex+3hsyjsc1wHZ96T57F78Sx8zpsF/A942Ou4Pfy9iAEmeR1rFnkWxYANQAX39RX+\nrpldairpToJU1QOquhyI9yLACArkWSxS1WPuyyVAhQjHGCmBPItTPi8vAw5GML5ICnSicFfgG+BA\nJIOLsECfRU7YQj2QZ/FnYIKq7gZQVb//RrJLUrFJkCmCfRYdgMlhjcg7AT0LEWkjIpuAn4BuEYot\n0tJ9FiJSHud/KB+6Rdm1wzWQ3wsFmrhNo5NFpFbEoousQJ5FVaCEiMwWkeUi8qS/C0b16C8f2fWX\nPyMCfhYi0hx4BrglfOF4KqBnoarfAd+JyK3A50D1sEbljUCexXtAf1VVERGy71/qgTyLlUBFVT0t\nIvcA3wGXbFiYDQTyLPIC9XCmbhQCFonIYlX9ObWTs0tS2QNU9HldESfj5kQBPQu3c/5fQEtVPRKh\n2CItqN8LVZ0nInlEpKSqHgp7dJEVyLOoD/zHySdcAdwjIvGqOikyIUZMus9CVU/4HP8kIqNFpISq\nHo5QjJESyO/FLpx9tM4AZ0RkLnADkGpSyS7NX8uBqiJSSUTyAY8Baf1DyK5/fSVL91mIyFXARKCd\nqm7zIMZICeRZXOv+VY6I1APIhgkFAngWqnqNqlZW1co4/SovZsOEAoH9XpT2+b1oiDP9IrslFAjs\n/53fA01FJLeIFAIaARvTumC2qKloGpMgReQF9/1/ikgZYBlwOZAkIt2BWqp60rPAwyCQZwG8AhQH\nPnT/3cSrakOvYg6XAJ/Fw8D/iUg8cBJo61nAYRTgs8gRAnwWjwAvikgCcJoc/HuhqptFZAqwFkgC\n/qWqaSYVm/xojDEmZLJL85cxxpgswJKKMcaYkLGkYowxJmQsqRhjjAkZSyrGGGNCxpKKMcaYkLGk\nYkwYiEiiu1T4OhEZLyIFfd77h4g0EZFPReQXn2X3G3sZszGhYEnFmPA4rap1VbUOcB7o6PNeI2Ax\nzrpLvVW1LtAfyDETEE32ZUnFmPCbD1QBEJGawBZVTXLfS142aJ7POT3dGs46d+UHY6JGtlimxZis\nSkTyAPeQsr3APcCUVE59AFjrrj/WHmefi1zAEhGZo6qrIxCuMZlmNRVjwqOgiKzCWW9uB/CxW34X\nKUlFgLfd857F2dvmVmCiqp5xNxCb6JYZExWspmJMeJxx+0r+4K7wWkxV97lFyX0qE33OuYMLV9IW\nbL8gE0WspmJM5DTH2f/d18VbMcwD2ohIQREpDLRxy4yJClZTMSY8Uqtd3AOM93eeqq4SkU+BpW7R\nv1R1TejDMyY8bOl7YyJERFYADVU10etYjAkXSyrGGGNCxvpUjDHGhIwlFWOMMSFjScUYY0zIWFIx\nxhgTMpZUjDHGhIwlFWOMMSFjScUYY0zI/D83MNuIXKP7/AAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f21b64ee750>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      " B = 47.00 \n",
      "\n",
      " Total surface area = 2.040 *10**5 m**2/kg\n",
      "\n",
      " V = 0.0*10**(-6)m**3/kg\n"
     ]
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "from matplotlib.pyplot import plot,title,show,xlabel,ylabel\n",
    "\n",
    "\n",
    "#For 1 m3 of pellet with a voidage ε, then\n",
    "        #Number of particles = (1 − ε)/(π/6)(15 × 10**−9)**3\n",
    "        #Surface area per unit volume = (1 − ε)π(15 × 10−9)2/(π/6)(15 × 10**−9)**3\n",
    "        #= 6(1 − ε)/(15 × 10**(−9)) m**2/m**3\n",
    "        #1 m**3 of pellet contains 2290(1- ε)kg solid and hence:\n",
    "specific_surface = 6/(15*10**(-9)*2290)\n",
    "print\"\\n Specific surface = %.3f*10**5 m**2/kg \"%(specific_surface*10**(-5))\n",
    "#(a) Using the BET isotherm\n",
    "        \n",
    "        #y = (P/Po)/V *10**(-6)\n",
    "y = [1500, 2660, 3576 ,4283 ,4613, 4615]\n",
    "        #y1 = (P/Po)/V *(1-(P/Po)\n",
    "y1 = [1666 ,3333 ,5109, 7138, 9226, 11358]\n",
    "        #x = P/Po\n",
    "x = [0.1 ,0.2 ,0.3 ,0.4 ,0.5 ,0.6]\n",
    "plot(x,y)\n",
    "xlabel(\"P/Po\")\n",
    "ylabel(\"(P/Po)/(V *10**(-6)) or (P/Po)/(V *(1-(P/Po))\")\n",
    "plot(x,y1)\n",
    "show()\n",
    "\n",
    "       #intercept, 1/VB = 300, and slope, (B − 1)/VB = 13,902 \n",
    "B = (13902/300)+1\n",
    "print\"\\n B = %.2f \"%(B)\n",
    "V = 1/(300*47.34)#       \n",
    "       #Total surface area\n",
    "S = ((70.4*10**(-6)*808*6.2*10**(26)*0.162*10**(-18)))/28\n",
    "print\"\\n Total surface area = %.3f *10**5 m**2/kg\"%(S*10**(-5))\n",
    "        \n",
    "#(b) Using the Langmuir form of the isotherm\n",
    "\n",
    "       #Slope 1/V = 13,902\n",
    "V = 1/13902\n",
    "print\"\\n V = %.1f*10**(-6)m**3/kg\"%(V*10**6)\n",
    "       #which agrees with the value of the BET isotherm\n",
    "#It may be noted that areas calculated from the isotherm are some 20 per cent greater than the geometric surface, probably due to the existence of some internal surface within the particles"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Page 1000 Example 17.2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX8AAAEPCAYAAACqZsSmAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XeYFFX28PHvIUkWQZKCEhUBydHEqAgogoq7CPpjVZRl\nSIK4oLii6KKr4gviIoriCriC7oqBJScbQRBEYIAZgiQRQSTnNDPn/aN6ZGChp2e6q6un+3yepx+6\nuqvqHko8c+fUrXtFVTHGGBNf8ngdgDHGmMiz5G+MMXHIkr8xxsQhS/7GGBOHLPkbY0wcsuRvjDFx\nyNXkLyIlROQzEVknIiki0szN9owxxgQnn8vnHwlMV9U/iEg+oIjL7RljjAmCuPWQl4hcCqxU1Squ\nNGCMMSbH3Cz7VAb2iMiHIrJCRN4XkcIutmeMMSZIbib/fEADYLSqNgCOAc+42J4xxpgguVnz3wHs\nUNXv/dufcV7yFxGbWMgYY3JAVSWU413r+avqr8DPInKN/6OWQPIF9ov61wsvvOB5DBanxWlxWowZ\nr3Bwe7RPH+BjESkAbAYedbk9Y4wxQXA1+atqEtDYzTaMMcZknz3hG4SEhASvQwiKxRleFmd45YY4\nc0OM4eLaOP+gGhdRL9s3xpjcSETQaL3ha4wxJnpZ8jfGmDhkyd8YY/ziqQxtyd8YY4ApG6bQZGwT\nTqed9jqUiHB7nL8xxkS1o6eP0n9Wf+ZumctH931EgbwFvA4pIqznb4yJW0t3LKX+mPqkpqeyKnEV\nN151o9chRYz1/I0xcSc1PZVXFr7C29+/zei7RnN/zfu9DiniLPkbY+LKpv2b6PJFF4oVKMbK7iu5\notgVXofkCSv7GGPigqoydsVYmn/QnM61OzPz/2bGbeIH6/kbY+LAnmN76Pbfbmw7uA3fwz5qlanl\ndUies56/MSamzfhxBvXG1OPaUtey9PGllvj9rOdvjIlJx88cZ+Ccgfx343/5uMPHJFRK8DqkqGI9\nf2NMzFmxawUN32vIgZMHSEpMssR/AdbzN8bEjLT0NF7/9nVGfDeCkW1G0vn6zl6HFLUs+RtjYsK2\ng9vo8kUX8ufJzw9//oGKl1b0OqSoZmUfY0yupqpMSJpA4/cbc++19zL3T3Mt8QfBev7GmFxr/4n9\nJE5NJGVPCnO7zKVuubpeh5RrWM/fGJMrzd0yl7rv1uXKYley/M/LLfFnk/X8jTG5ysnUkwyaO4jP\n1n3Gh/d8SMsqLb0OKVey5G+MyTWSfk3ioc8f4rrS15GUmETJQiW9DinXsrKPMSbqpWs6byx+g5Yf\ntWTgjQP59x/+bYk/RNbzN8ZEtZ8P/czDXz7M6bTTLHt8GZUvq+x1SDHBev7GmKj1ydpPaPheQ+6o\ncgcLHllgiT+MrOdvjIk6B08epNf0Xvyw8wdmPDSDhlc09DqkmGM9f2NMVPFt81H33bpcVvAyVnRf\nYYnfJdbzN8ZEhVOpp3j+6+f5aPVHjG0/lruq3+V1SDHN9eQvItuAw0AacEZVm7jdpjEmd0n+LZmH\nPn+ISiUqkZSYROkipb0OKeZFouevQIKq7o9AW8aYXCRd0xm1bBR/++ZvvHr7q3St3xUR8TqsuBCp\nso/91zTGnGPnkZ088uUjHD51mCWPLaFayWpehxRXInHDV4G5IrJcRLpFoD1jTJSbnDKZ+mPqc2PF\nG1nUdZElfg9Eoud/o6ruEpHSwBwRWa+qCyPQrjEmyhw+dZi+M/uyaPsipnSaQtMKTb0OKW65nvxV\ndZf/zz0i8gXQBPg9+Q8ZMuT3fRMSEkhISHA7JGOMB77d/i1dvuhCyyotWdl9JUULFPU6pFzD5/Ph\n8/nCek5R1bCe8JyTixQG8qrqEREpAswGXlTV2f7v1c32jTHeO512mpcWvMTYFWMZc/cY7qlxj9ch\n5XoigqqGdC816J6/iBQEVFVPZeP8ZYEv/Hfv8wEfZyR+Y0zs+/6X7+k6pSuVS1RmVeIqyhUt53VI\nxu+iPX8RyQPcC3QGbsC5OSw44/WXAB8DX4bSdbeevzGx6cSZEwzxDWFc0jjebP0mnWp3siGcYeR2\nz9+HU5t/A1iV0eMXkUuA+kB74EngllACMMbElsU/L6brV12pU7YOqxNXU7ZoWa9DMhcQqOd/SVYl\nnmD2yeJ46/kbEyOOnT7GX+f/lU+TP2XUnaO4v+b9XocUs1zt+WdO6v6btVcB6cAOVT12/j7GmPjl\n2+bjsSmP0bxCc9b2WEupwqW8Dslk4aLJX0SKAd2ATsDlwG6cmn9ZEdmHU/N/X1WPRiJQY0z0OXLq\nCE/PfZopG6bwTtt3aHdtO69DMkEK9ITvl8ARoJ2qVlHV5qraTFUrA3cDx4CvIhGkMSb6zNo0i9rv\n1OZ02mnW9lxriT+XcXWcf5aNW83fmFzn4MmDPDXrKeZtncd77d6jVdVWXocUd8JR8w9qbh8RuUdE\n/p//ZT/ejYlTUzdOpfbo2lyS7xLW9FhjiT8Xy7LnLyKvAo1xavyCcw9guaoOCrlx6/kbkyvsO76P\nvjP7smTHEj5o/wEJlRK8DimuhaPnH0zyXwPUU9U0/3ZenHH/14fSsP9clvyNiXKTUybTZ0YfHqj1\nAENvG0qRAkW8DinuRWp6BwVKAPv82yX8nxljYtjuo7vpPaM3a3av4bOOn3FDxRu8DsmEUTA1/78D\nK0RkvIiMB34AXnE3LGOMV1SViWsmUufdOlS9rCqrEldZ4o9BwZR9LgcK4NT9Ffg+Y5rmkBu3so8x\nUWXnkZ0kTk1k68Gt/LP9P2l8ZWOvQzIX4OpoHxFpJyJ7gDXAMmCvqk4JV+I3xkQPVeXDlR9S7916\n1C9Xn+Xdllvij3GBav6vADer6noRaQq8DrSITFjGmEjZfmg7f/7vn/nt2G/M6TKHuuXqeh2SiYBA\nNf9UVV0PoKpLgeKRCckYEwnpms67y9+l4XsNueXqW1j6+FJL/HEkUM+/tIj0xxnbf/62qupw16Mz\nxrhiy4EtPD7lcY6dOYbvYR+1ytTyOiQTYYF6/mOBYkBR/yvzdjH3QzPGhFu6pjPyu5E0eb8Jbau3\nZXHXxZb441Sgnv9GYLaq7o1UMMYY92zYu4HHpjxGHsnDkseWUL1Uda9DMh4K1PO/Cvi3iCwSkSEi\n0lRsHTZjcp3U9FSGfTuMmz68iU61O+F7xGeJ3wQ1zr840BJoDTQB1gMzgFmqujukxm2cvzGuWvvb\nWrp+1ZVilxRjbLuxVL6sstchmTCIyNw+F2i0FnAn0EpVQ5rSz5K/Me44k3aG1759jZFLR/LybS/T\nrUE3W0A9hria/EWkLPAsUA1YDfxdVQ+H0tgF2rDkb0yYrdy1kq5TulK+aHnG3D2GipdW9DokE2Zu\nz+c/ATgK/ANndM9boTRkjHHXqdRTDJ4/mNb/ak2/pv2Y9uA0S/zmogKN9imnqn/1v58pIisjEZAx\nJvuW/bKMrl91pVrJaiQlJlG+WHmvQzJRLlDyFxEpmfEeyJtpG1Xd72pkxpgsnThzghd8LzAhaQJv\ntnmTB2o9YLV9E5RAyb84zvTNmWVsK1DFlYiMMVlSVaZunEr/2f1pUL4Bq3uspkyRMl6HZXIRW8Dd\nmFwmZU8KT856ku2HtjOi9QjaVGvjdUgmwiK2gHumBoeE0pgxJuf2n9jPEzOeoMW4FrSt3pbViast\n8Zscy1byB+5xJQpjzEWlpqcy+vvRXPf2daSmp7Ku1zqeaPoE+fPm9zo0k4sFs4ZvSPwLvi8Hdqhq\nO7fbMyaWzN86n74z+1K6cGnmdJlDnbJ1vA7JxIjsJv9GOWijL5CCzQRqTNC2HNjCX2b/hVW/ruKN\nVm9wX437bBSPCatslX1UNQ1ARJ4PZn8RqQDchTMdtP3LNSYLR04d4dl5z9Lk/SY0vqIxKb1S6HBd\nB0v8JuyyW/PP0C3I/UYAA4D0HLZjTFxI13QmJE2gxts1+OXIL6zusZpBNw+iYL6CXodmYtRFyz4i\nciTAcYWyOrGI3A38pqorRSThYvsNGTLk9/cJCQkkJFx0V2Ni0nc7vqPvzL4AfN7xc5pWaOpxRCba\n+Hw+fD5fWM8ZaGK37UATVf31At/9rKoBJw0RkVeALkAqUBDnobHJqvqnTPvYOH8Tt345/AuD5g1i\n3tZ5vHr7qzxU5yHySE5/GTfxxO1x/h/hLOhyIZOyOrGqPquqFVW1MtAJmJ858RsTr06mnuSVha9Q\n9926VCxekQ29N9ClbhdL/CaiLlr2yTSp24W+G5iDtqyLb+KaqvL5us/5y5y/0KB8A5Z1W0aVy2yW\nFOONQGWfKqq6JeDBIlVVdXOOG7eyj4kTq3evpu/Mvuw7vo8327zJbZVv8zokk4u5vZjLp0ARYArO\nQ1q7cIZrlscZ798eOKKqnXLcuCV/E+P2Ht/L4PmD+Xz95wxpMYRuDbuRL4/rz1aaGOf6Mo4iUg2n\nXn8jcLX/45+ARcCkrH4zyLJxS/4mRp1JO8Po70czdOFQHqz9IC8kvEDJQiWzPtCYIHiyhm84WfI3\nsWjWpln0m9WPqy69ihGtR1CzdE2vQzIxJhzJP9A4fwGGqupfReRvqjo4lIaMiXUb923kqdlPsX7v\neka0HkHb6m3tyVwTtS46tszfJV8mIm/j1PyNMRdw6OQhBswewA0f3ECLq1uwtsda7r7mbkv8JqoF\n6vm/AJQEOgOpIlJPVV+MWGTGRLm09DQ+XPUhg78eTNvqbUnumUzZomW9DsuYoGR1w3ck8A6QqKr9\nwt641fxNLrXwp4X0ndmXwvkLM7LNSBpe0dDrkEwccXuoZ36glqquEpG6qpoUSkMXacOSv8lVth/a\nzsA5A1n882Jev+N1WzDdeMLt6R1+AXqLyO3A6lAaMSa3O37mOEN8Q6g/pj41Lq/B+t7r6VS7kyV+\nk2sFSv41cW70DgZ2iMhIEWkWmbCMiQ6qyidrP6HGqBqs37ueld1XMiRhCIXzF/Y6NGNCEtQ4fxG5\nAugIPACUAT5V1WdDbtzKPiaK/bDzB/rO7MvxM8cZ2WYkN199s9chGQNE+CEvESkGdAD6A+VVtUwo\nDfvPacnfRJ09x/bw7LxnmfrjVP526994tN6j5M2T1+uwjPmd2zV/RKSQiHQUkc+BTcBtwNPAFaE0\nakw0Sk1P5e1lb1NrdC2KFijK+l7rebzB45b4TUwKNM5/InAHsAD4GHhIVU9EKjBjIunb7d/Sa3ov\nShQswfyH51O7TG2vQzLGVYGmF5yJM77/cKSCMSbSfj36K0/PfZp5W+bxRqs3bOimiRuByj55gOMX\n+1JECojIo+EPyRj3nUk7w4glI6g9ujblipSzoZsGgHi6BRmo518U+F5E1gPfA7/izOdfDmc+/xrA\n+65HaEyY+bb56D29N1cUu4JFXRdR4/IaXodkPHDyJCQnw6pVzispCX76CbZtg3joA2Q1vYPgzOV/\nE2fX882Yz39xqEN1bLSPiaQdh3cwYM4AFv+8mBGtR3Bfjfuspx8n9uxxknvmRL9pE1SvDvXqOa+6\ndZ3X5Zd7HW3W3J7e4UFglqruC6WBgI1b8jcRcDrtNCOWjGDY4mH0aNSDQTcPsoe0YlR6OmzefDbJ\nZyT6I0fOJvmMRF+zJhQs6HXEOePqfP44Pf3/iEgBYC4wA1hm2drkJrM3z6bPjD5UL1md7x7/jmol\nq3kdkgmT48dh7dpzE/2aNVCq1Nkk362bk+grVYqPUk52ZPmQl4gUB1oCbYDGwHqcHwSzVHV3SI1b\nz9+45KeDP9F/dn9W7lrJyDYjaXdtO69DMiH49ddzyzarVjm1+Ro1zu3R16kDl13mdbTu82QZRxGp\nBdwJtFLVViE1bsnfhNnJ1JO8sfgN3vzuTfo27cuAGwdQMF8u/d0+DqWlwcaN/5voT58+N8nXq+ck\n/gIFvI7YGxFL/iJyHc4C7gr8pKrrQ2k003kt+ZuwmbZxGn1n9qVO2ToMbz2cSiUqeR2SCeDoUVi9\n+mxdftUqp4xTvrxTqsmc6CtUsLJNZm7f8K0MPAnchTO9806coZ7lgQrAVGCEqm7LceOW/E0YbN6/\nmX6z+rFh7wb+cec/aF2ttdchmfPs3AkrV56b6HfsgFq1zk30depA8eJeRxv93E7+/8YZx+9T1TPn\nfZcfuBV4XFU75rhxS/4mBMfPHOfVRa8y+vvRDLhhAP2a9eOSfJd4HVbcU4Uff4SFC+Gbb5w/Dx2C\nhg3PTfTXXgv5Ag05MRflSc0/nCz5m5xQVb5c/yVPznqSZhWa8UarN6hQvILXYcWttDRnlE3mZJ8/\nP9xyC9x8s/NnjRqQJ+A0kiY7IpL8/UM9ewC3+D/yAe+e/9tAjhq35G+yaeO+jfSZ0Ycdh3fwjzv/\nwW2Vb/M6pLhz+jQsX3420S9eDGXLOok+I9lffbXV6N0UqeT/Ac7zAONxav5dgFRVfTyUhv3ntuRv\ngnL09FFe/uZl3l/xPs/e/Cx9mvQhf978XocVF44ehe++O5vsv/8errnmbKK/6SYn+ZvIiVTyX62q\ndbL67CLHFsSZEvoSoADwlaoOyvS9JX8TkKryn5T/8NTsp0iolMDrLV+nfLHyXocV0/btg0WLzpZx\nkpOhfv2zZZwbboBLL/U6yvjm9hO+GVJFpJqqbvI3WhVIDebkqnpSRG5V1eMikg9YJCI3qeqiEGI2\ncSL5t2T6zOjDvhP7mNhhoi2j6JIdO5xEn5Hst2+HZs2cZD9sGDRpAoUKeR2lCbdgkv8AYL6IbPVv\nVwKCnspZVTOmhS4A5AX2ZydAE38OnzrMi74XmbB6Ai+0eIHERonky2PDQsLhQiNxDh8+W69/9FGn\nl2+jcGJfwP/EIlIaOAQ0xVm4HWCDqp4MtgERyQOsAKoC76hqSg5jNTFOVfl4zcc8Pfdp2lRtQ3LP\nZMoUCXmp6LiW1Uicp5+2kTjxKtAyjo8DrwCbgSrAn1X1q+w2oKrpQD0RuRSYJSIJqurL+H7IkCG/\n75uQkEBCQkJ2mzAxIOnXJHrP6M3xM8eZ3HEyzSo08zqkXCnQSJz27Z0yjo3EyX18Ph8+ny+s5wz0\nkFcykKCqe0SkCjBRVUP6P1JEBgMnVPUN/7bd8I1zB08eZPD8wfw75d+8lPCSLZieTTYSJz65fcP3\ntKruAVDVLSKS7UcnReRynGGhB0WkEM6C8C/mLFQTS9I1nXGrxvHsvGe5t8a9pPRMoVThUl6HFfWO\nHoVvvwWfz3mtXn12JM7AgTYSxwQvUPKvICJv4YztB7gy07aq6hNBnL88MN5f988DfKSq80KK2OR6\ny3cup/f03gBMe3AaDa9o6HFE0evIEWfY5YIFTrJfu9aZJiEhAV5+2RmVU9jWpTE5EKjs8wjOLJ6Z\nf7XI2FZVHR9y41b2iSsHTx7k2XnP8sX6L3jltld4uN7D5BG705jZ4cPnJvvkZGjUyEn2LVo4yd6G\nXRq3yz5XADNUdWUoDRijqnya/Cn9Z/Wn/bXtSemZwmWF4mDFjSAcOnRusk9JgcaNnWT/+uvQtGnu\nXWrQRLdAPf9OOKt31QOSgOnAbFU9ELbGrecf8zbv30zP6T3ZdWQXY+4eQ/OKzb0OyVOHDjk3ZjOS\n/bp1zkNUCQnOq0kTS/Yma5Ga3kGA+jg/CO7A+W1hDjBTVZeF1Lgl/5h1Ou00byx+g+FLhvP0jU/T\nr1m/uJyL5+BBJ9n7fE7C37DB6c23aHE22V9is1CbbPJqGcdLcX4ItFbVbiE1bsk/Ji38aSGJ0xKp\nXKIyo+4aFVcrah04cG6y37jRqdNnJPvGjS3Zm9C5vZjLNcAwoBqwGviLqv4SSmMXaMOSfwzZd3wf\nT899mpmbZjKyzUg6XNcBifGnifbvP5vsfT7YtAmaNz97g7Zx4/hdZ9a4x+3kvwhnGueFQDuguap2\nCKWxC7RhyT8GqCofrf6IgXMG0rFWR4beNpTil8TmWnz79zsPVGUk+y1bzk32jRpZsjfuczv5r1LV\nepm2V6pq/VAau0AblvxzuQ17N9BjWg8OnTrEmLvH0OiKRl6HFFb79p2b7LdudR6kyrhB27ChM1eO\nMZHk9lDPgiLSIKMtoJB/O2Oc/4pQGja528nUk7y66FVGLRvF4FsG06tJr5iYeVMVVqyATz6BWbPg\np5/OJvv33oMGDSzZm9gQqOfvw3mo6/ePMm+r6q0hN249/1xp/tb5JE5NpHaZ2rx151sxsX7uunUw\naZKT9NPToVMnaNfO6dnb9MYm2tgC7iaifjv2G0/NfopvfvqGUXeOot217bwOKSTbtjnJ/pNPYM8e\neOAB6NzZqdvH+H1qk8uFI/lf9Nl6EWkRRAAh9/5N9EvXdMauGEvt0bUpW6QsyT2Tc23i//VX+Mc/\nnFJO48bOD4CRI53Vq4YPdz6zxG/iQaBfaNuJyDBgLrAc2IXzw6Ic0AhoCXztf5kYlfxbMonTEjmT\ndoY5XeZQt1xdr0PKtgMH4PPPnbLODz845ZznnoM77rD6vYlfAcs+IlIMuAe4Ebja//FPwCKcxdiP\nhtS4lX2i1vEzxxn6zVDeX/E+Lya8SPeG3XPVPPvHjsGUKU7CX7AAWrZ0Sjpt29rEaCb3s5q/ccXM\nTTPpNb0Xja9ozIjWIyhfrLzXIQXl1CmYOdOp4c+Y4Yy/79wZ7r0XisfmYwcmTlnyN2G168gunpz1\nJMt+WcbotqNpU62N1yFlKTXVGX8/aRJ8+SXUru0k/Pvvh9KlvY7OGHdY8jdhka7pjFk+hud9z9Ot\nQTeeu+U5CueP3hVCVGHJEifh/+c/ULGiMzTzgQegQu4fdWpMltx+yAv/ClzNVHVxKI2Y6JX0axLd\np3YnX558+B72UatMLa9DuiBVSEpyEv6nnzqrV3Xu7MyrU72619EZk/sEM6XzOdM8hLVx6/l75ujp\nowzxDWFC0gReuf0VutbvGpWram3c6NTwJ02CkyedHn7nznD99TYk08Qv13v+fnNF5A/AZMvUseG/\nG/5Lnxl9uOXqW1jbcy1lipTxOqRz/Pyz07ufNAl27oSOHeGf/3SmRraEb0x4BNPzPwoUBtKAk/6P\nVVVDHj9hPf/I2nF4B0/MeIK1v63lnbbvcHuV270O6Xd79jj1+08+cdatve8+p4efkAB5c88IU2Mi\nwm74mqCkpqcyatkohn4zlN5NevPMTc9QMJ/3awUeOuSM0Jk0Cb77Du66y0n4rVvbtMjGBBKpsg8i\ncg9wC87EbgtU9b+hNGoiZ/nO5XSf2p3ilxTn267fcu3l13oaz4kTMHWqk/DnzYNbb4VHH4XJk6FI\nEU9DMyauBFP2eRVoDHyMM7NnJ2C5qg4KuXHr+bvm8KnDDJ4/mE+TP+X1O16nS50unq2qpQqLF8P4\n8fDZZ85MmZ07Q4cOUKKEJyEZk6tFquffFqinqmn+RscBq4CQk78JP1Xl83Wf029WP1pVaUVyz2RK\nFS7lSSzbt8OECc4rb154+GFYvdrG4hsTDYJJ/gqUAPb5t0tw7jz/JkpsO7iN3tN7s+XAFiZ2mMjN\nV98c8RiOHXMmURs/HlaudB68+ugjaNLERuoYE02CSf5/B1aIyNc4ZZ8WwDOuRmWy5XTaaUYsGcGw\nxcPo37w/nz/wOQXyRu6OaXo6LFoE48bBF1840yV37+7MnlnQ+/vKxpgLCGq0j4hcgVP3V+B7Vd0V\nlsat5h+y+Vvn02t6L6pcVoW32rxF1ZJVI9b21q1OSWf8eOeJ20cegYcegvK5Yx44Y3ItG+oZx3Ye\n2clTs59iyc9LGNlmJO2vbR+RG7pHjzo3bceNc8bjd+7s1PIbNLCyjjGR4upKXuEgIhVF5GsRSRaR\ntSLyhJvtxYMzaWcYvmQ4dd6pQ5USVUjplcI9Ne5xNfGnp8P8+U6Sr1DBKe088QTs2AFvveWM3rHE\nb0zu4mrPX0TKAeVUdZWIFAV+AO5V1XX+763nnw0Lf1pIr+m9KFu0LKPuHOX6mP1Nm5ySzoQJcNll\nTlnnwQehTHTNBmFM3HF1qKeITAcmAl/mdMUuVf0V+NX//qiIrAOuANbl5HzxavfR3QycO5D5W+cz\nvNVw/lDzD6719A8dcqZZGDcOfvzRSfZffQX1XJnazxjjlUBln/eAu4GtIvJvEblPRHI8hEREKgH1\ngaU5PUe8SUtP4+1lb1P7ndqUKVyGlJ4p/LHWH8Oe+NPSYM4c52bt1VfD9OkwcKBT1hkxwhK/MbHo\noj1/Vf0S+FJEigDtgIeBd/2/EUxS1dnBNuIv+XwG9D3/t4ghQ4b8/j4hIYGEhITsxB+zvtvxHT2n\n9aT4JcVdm2d/w4azZZ1y5Zya/siRcPnlYW/KGBMCn8+Hz+cL6zmzVfMXkbrAeOB6VQ1qrkURyQ9M\nBWao6pvnfWc1//PsPb6XQXMHMe3HaQy7YxgPXv9gWHv6Bw860yWPGwfbtsH//Z+T9GvXDlsTxhiX\nRWR6B/9N2444c/qUBz7F+S0gmAAF+ABIOT/xm3OlazpjV4xl8NeD6VSrE+t6rePSgpeG5dypqU5Z\nZ/x4Z4HzVq3gueec2TPzBTW1nzEm1ly05y8if8ZJ+DWAycAkYEl2uuoichPwDbCas1NCDFLVmf7v\nrecP/LDzB3pO70m+PPkYfddo6parG5bzJic7Cf9f/3LWuX3kEWe6hZIlw3J6Y4xH3O75N8OZ2mGe\nqqbn5OSqugiXnyXIzQ6cOMBz859j8rrJ/P32v/NwvYdDXkpx3z5nQZRx45xVsLp0caZOvu668MRs\njIkNgZJ/XaAUUFVEZqrqtsiEFPvSNZ0JSRMYNG8Q99W4j5ReKZQslPPuuKrzENY778DcuXDnnTB0\nKLRsaatgGWMuLOANXxGpDLQBWgMVgIXADJwFXU6F3Hgcln1W715Nz2k9OZ12mtFtR9PoikY5Ptep\nU04vf/hwp67fp4+zwLnNkW9MbIvo3D7+Mf434/wwaAHsUdW2ITUeR8n/8KnDPP/180xcM5Ghtw3l\nsfqPkTdPzrrle/fCu+/C6NFw/fXQv79zE9emWDAmPkRsGUcAVT0NzPO/EBFbkiMIqsqktZMYMGcA\nd1a7k5ReKVxeOGcD6devhzffdIZqdugAs2fbEE1jTM4ETP4i0gDojLN+byWcETs/4YzgmQjscDm+\nXC1lTwrVMpiSAAAOyElEQVS9pvfi4MmDfPbHz2hesXm2z5FRzx8+HJYvhx49nB8CZcu6ELAxJm4E\nGuo5HTgATAGWAbtwFnMpDzTBeeq3RCiln1gt+xw9fZSXFrzEh6s+5IUWL9CjUY9sl3jOr+c/+aQz\n/UKhQi4FbYzJNVyt+YtIWVXdnUUAZVT1txw3HmPJX1WZvG4y/Wf1J6FSAsPuGEbZotnrols93xiT\nFVdr/pkTv/8p34yVvJZlJPxQEn+s2bhvI31m9GHnkZ38q8O/uOXqW7J1vNXzjTGRlOUTRSLSEWcm\nzj/iTPOwTET+6HZgucXxM8d5bv5z3PDBDbSu2poVf14RdOJXdR7AatsWWrRwJldbvx4++MASvzHG\nXcGM9nkOaJzRyxeR0jgjfv7jZmDRTlWZsmEK/Wb1o+mVTUlKTOLK4lcGdWzmev6ZM05p57PPrJ5v\njImcYJK/AHsybe/zfxa3thzYwhMznmDzgc2MbTeW26vcHtRx59fzX3vNmVzN6vnGmEgLZiKZmcAs\nEXlERB4FpuM85Rt3Tqae5KUFL9Hk/SbcdNVNJCUmBZX416+HxESoXh22boVZs5xXmzaW+I0x3siy\n56+qA0SkA3CT/6MxqvqFu2FFnxk/zqDPjD7ULVeXFd1XcNWlVwXc//zx+YmJNj7fGBM9Aq3hew0w\nDKiGMyXzAFWNu4e6fj70M31n9mX17tWMumsUbaq1Cbi/1fONMblBoLLPP3FW4LofWAG8FZGIokRq\neirDlwyn/pj61C1bl7U91wZM/Hv3OjNpVq4MEyc69fzkZHj8cUv8xpjoE6jsU1RV3/e/Xy8iKyMR\nUDRY9ssyuk/tTslCJVn82GKuKXXNRfc9f3z+rFnOzVxjjIlmgZJ/Qf/cPuCM7ink3xZAVXWF69FF\n2KGTh/jr/L8yed1kht0xjIeuf+iC6+daPd8Yk9sFmt7Bx9mlF8Gf9DM2VPXWkBuPkukdVJXPUj7j\nyVlPclf1u3i15asXXFzlQvV8m2/HGBNpbs/tc6Wq/hLKybNsPAqS/9YDW+k1vRfbD23n3bvf5aar\nbvqffc4fn//kkzY+3xjjnXAk/0A3fN8XkaUi8qqIJIhI0HP/5wZn0s7w2qLXaPx+Y26+6mZWdF/x\nP4l/714YMMDG5xtjYk+gid3uEpFCQALQAXhDRH7GecBrpqpuj0yI4bf458V0n9qdCsUrsKzbMqpc\nVuWc7w8fhhEj4K23oGNHWLsWrgxu5gZjjMkVgl7GEUBEqgB34qzpW05Vm4TUeITLPgdOHOCZuc8w\n9cepjGg9gj/W/OM5N3RPnHAWQX/tNWca5SFDoGrViIVnjDFBcbXsI477RGSAiLQGUNUtqvq2qrbH\nWc83V1BVJq6ZSM3RNcmXJx8pPVPoWKvj74n/zBl47z245hpYuNCZafOjjyzxG2NiV6A6/migJrAY\n+JuINFXVlzK+VNVTbgcXDpv2b6LHtB7sObaHLx/4kqYVmv7+XXq6Mz7/+efh6qudJ3GbNg1wMmOM\niRGBRvskA3VUNU1ECgOLVLXBBXfOaeMuln1OpZ7i9W9fZ+TSkQy6aRB9m/UlXx7nZ50qTJ0Kzz0H\nBQvCK6/A7cFNzGmMMZ5zdSUv4LSqpgGo6nG50NNOUWrBtgUkTkukesnq/zMJm88Hzz7r3NR9+WVo\n395G7hhj4k+gnv8JYFOmj6oCm/3vVVXrhNx4mHv+e4/vZcCcAczdMpe32rzFvTXu/b2u//338Ne/\nwubN8NJL0KkT5M3emurGGBMV3O75XxfKiSNJVRm3ahzPzHuGB2s/SErPFIpdUgyAlBSnvLN0KQwe\nDI89BvnzexywMcZ4LFDy/ymrbrlk0XUXkX8CbYHfVNWV6c7W711P4tREjp4+yoyHZtCgvHNbYutW\nZ6jmjBkwcCB8/LFNw2CMMRkCPeHr8w/z/J8pLUXkWhF5GliQxfk/BAJPgJ9DJ1NP8vzXz3Pzhzdz\n/3X3s/TxpTQo34Bdu6BXL2jUCCpVgh9/hL/8xRK/McZkFqjn3wp4CHhbRGoDR3AmdysKrAU+BloG\nOrmqLhSRSmGJNJO5W+bSY1oP6paty6ruq7iy+JXs3w+vv+6M13/0UWeWzdKlw92yMcbEhkDTO5zC\nWdDlnyKSF7jc/9XejFFAkbb76G76z+7Pt9u/5e273qbtNW05etQZtTNihDOfflISVKzoRXTGGJN7\nBLOAO6qapqq7/a+IJ/50Tee9H97j+neup0KxCiT3TOb2q9oyciRUq+asmLVkidPrt8RvjDFZ83ym\nziFDhvz+PiEhgYSEhHO+X7N7DYnTEknXdOb+aS41S9Vh/Hh48UWoW9eZZbNu3cjGbIwxkeTz+fD5\nfGE9Z7YmdstRA07N/78XGu0TaLDQ8TPHeWnBS3yw8gOG3jqUx+p34/PJeRg8GMqXd57KveEGV0M3\nxpio5PY4/5CJyCSgBVDKPx3086r6YVbHTf9xOr2m96J5heasTlzDyoXlaNwI8uWDUaOgZUt7KtcY\nY0LhavJX1c7Z2X/nkZ30m9mPFbtWMObuMRTa2YqOd8G+fTB0KNx3nyV9Y4wJh6Bu+LotLT2NUctG\nUffdulxT6homNF/DiN6t+NOfoFs3WLPGGcljid8YY8LD8xu+K3etpPvU7hTMV5APb1nAhP9Xk38u\ncubh+eorKFDA6wiNMSb2eN7zb/NxG/5YOZGqC308endNGjZ0nsrt1csSvzHGuMXznn/7HWt59e+l\n6dHDSfolSngdkTHGxD7Pk38RSpOSAmXLeh2JMcbED9fH+QdsPMILuBtjTCxwdQF3Y4wxscuSvzHG\nxCFL/sYYE4cs+RtjTByy5G+MMXHIkr8xxsQhS/7GGBOHLPkbY0wcsuRvjDFxyJK/McbEIUv+xhgT\nhyz5G2NMHLLkb4wxcciSvzHGxCFL/sYYE4cs+RtjTByy5G+MMXHIkr8xxsQhS/7GGBOHLPkbY0wc\nsuRvjDFxyJK/McbEIVeTv4i0EZH1IvKjiDztZlvGGGOC51ryF5G8wCigDVAT6Cwi17nVnpt8Pp/X\nIQTF4gwvizO8ckOcuSHGcHGz598E2KSq21T1DPAJcI+L7bkmt/yDsDjDy+IMr9wQZ26IMVzcTP5X\nAj9n2t7h/8wYY4zH3Ez+6uK5jTHGhEBU3cnRItIMGKKqbfzbg4B0VX0t0z72A8IYY3JAVSWU491M\n/vmADcDtwE5gGdBZVde50qAxxpig5XPrxKqaKiK9gVlAXuADS/zGGBMdXOv5G2OMiV5ujvPP8gEv\nEXnL/32SiNTP9Pk2EVktIitFZJlbMQYTp4jUEJElInJSRJ7KzrFRFGc0Xc+H/P+9V4vItyJSJ9hj\noyTGaLqW9/jjXCkiP4jIbcEeG0VxRs31zLRfYxFJFZH7s3tsFMQZ/PVU1bC/cMo8m4BKQH5gFXDd\nefvcBUz3v28KfJfpu61ASTdiy0GcpYFGwFDgqewcGw1xRuH1bA5c6n/fJuO/e6SuZygxRuG1LJLp\n/fU4z9VE47/NC8YZbdcz037zganA/dF4PS8WZ3avp1s9/2Ae8GoPjAdQ1aVACREpm+n7kO5khytO\nVd2jqsuBM9k9NkrizBAt13OJqh7yby4FKgR7bBTEmCFaruWxTJtFgb3BHhslcWaIiuvp1wf4DNiT\ng2O9jjNDUNfTreQfzANegfZRYK6ILBeRbi7FmFUMbh6bXaG2Fa3X8zFgeg6PzalQYoQou5Yicq+I\nrANmAE9k59goiBOi6HqKyJU4ifadTLEFdWwYhRJnxvugrqdbo32CvYt8sZ9QN6nqThEpDcwRkfWq\nujBMsWUWyt3uSN4pD7WtG1V1VzRdTxG5FegK3JjdY0MUSowQZddSVb8EvhSRm4GPRKSGC7EEDCGo\nnc6LE7jW/1U0Xc83gWdUVUVEOJufou3/9YvFCdm4nm71/H8BKmbarojzEyzQPhX8n6GqO/1/7gG+\nwPlVyKs43Tg2u0JqS1V3+f+Miuvpv4H6PtBeVQ9k51iPY4y6a5kproU4nbmS/v2i8t9mRpwiUsq/\nHU3XsyHwiYhsBe4HRotI+yCPjYY4s3c9XbppkQ/YjHPTogBZ3/Btxtkbf4WBYv73RYBvgVZexZlp\n3yGce8M36GM9jjOqridwFc4NrWY5/Tt6GGO0XcuqnB2u3QDYHI3/NgPEGVXX87z9PwQ6ROP1DBBn\ntq5n2IPPFNSdOE/4bgIG+T/rDnTPtM8o//dJQAP/Z1X8f+FVwNqMY72KEyiHU4M7BBwAtgNFL3Zs\ntMUZhddzLLAPWOl/LQt0bDTFGIXXcqA/jpXAQqBxpK9lKHFG2/U8b9/fk2q0Xc+LxZnd62kPeRlj\nTByyZRyNMSYOWfI3xpg4ZMnfGGPikCV/Y4yJQ5b8jTEmDlnyN8aYOGTJ38QFEblERBaISBUROeGf\n8jZZRN7xPyKfsd90EblSRHz+aXVXicgiEbkmwLnLisj0i31vTDSy5G/ixUM409+m4cyaWB+oA9QE\n7gUQkUJAKVX9BWeOlQdVtR7O7LPDLnZiVd0NHBCRBu7+FYwJH0v+Jl50Br4i0yRYqpoGLAaq+T9K\nAL6+wLELM/YRkWEissa/YEbHTPtM8bdhTK5gyd/EPBHJC9RW1Y3nfV4YuB1Y7f/oTmBm5l38f7YD\nVotIB6Auzm8MLYFhIlLOv88y4BZ3/gbGhJ8lfxMPLgeOZNquKiIrgUXAVFWd5f/8Bv9n4CT+j/37\nNQcGADcBE9XxG7AAaOzffxfOZFzG5ApuzedvTLTJPOf5Zn/N/+yXIlWAn1U11f9RRs1/RaZ9zj9P\nxn4Zn9tEWSbXsJ6/iQd7cWY4DeROnFWmMjs/0S8EHhCRPP7FMm7BKfcAlAd+CjVQYyLFev4m5qlq\nmoisFZFrgVNcuIfeGuh9/qHnnecLEWmOMwW5AgP85R9wFs34JryRG+Mem9LZxAUReQQoq6qvXeC7\nS4CFqprjVaRE5GPgDVVdmfMojYkcK/uYeDERaJv5ga4MqnoqxMRfBihhid/kJtbzN8aYOGQ9f2OM\niUOW/I0xJg5Z8jfGmDhkyd8YY+KQJX9jjIlDlvyNMSYO/X+upVVCFtJM8gAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f21b6435750>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      " S = 323527 m**2/kg\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEZCAYAAAB8culNAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XecVPW9//HXh6UoCio2FBEsqCh6ARUxihkbLqKixiho\nYkuu3NyoMYmCmt8Nm3JvgiWJiZigoMaYiCUWQCyxjIUoRRFFQUElUqJYUFBRgf38/vieZYfZwuzu\nnD1T3s/HYx47p3Dmc47jfvbbzd0RERHJ1CbpAEREpPAoOYiISB1KDiIiUoeSg4iI1KHkICIidSg5\niIhIHUoOUlLM7Fwze6aBY2eZ2SOtHVNrMrOUmS1JOg4pfkoOkhgzW2xmR2fta/CXe0u5+1/d/bg4\nrl3fvcQlekbrzWy1mX1iZnPMbGgzrnOrmf0ijhil+Ck5SJI8euWFmbXN17WaoVn3YpFmfN50d+8E\nbA1MBO4ys62bcR2Reik5SKHZ6BesmV1uZovMbJWZvWpmJ2ccO9fMppvZb8zsA2BMPf/+ajN7xsw6\nZ5dKzKzazEaa2RtmttLMrs84tqeZPWVmH5vZ+2Y2KdcbMLMqM/tLxnbP6LPaRNtpM/ulmU0HPgN2\nN7PzzOy16D7fNLMLNvUxAB6mOLgF2BzYvZ5Yekeft9LM5pnZidH+C4AzgVFRCeSBXO9PykOSf2mJ\nQPRLrpHtRcDh7v6umZ0O3G5me7j7e9HxAcDfgB2A9sBwCH+RAzcCuwDHuvsXDfyBPhQ4CNgKeMHM\nJrv7o8AvgIfd/etm1j46J1e5lCC+BQwBXif8kfYeMNTd3zazI4CHzGyWu89p7CJRaem7wGpgIXBg\nxrF2wBRgAnAMMAh4wMwOcvcbzexQYIm7/7QJ9yZlQiUHSZIB90d/1a40s5XAODJ+ubr7Pe7+bvT+\nLsIvwEMyrrHc3ce5e7W7fxHtawdMIlS5nJixvz6/dvdV7r4EeBLoG+3/CuhpZt3c/St3/2cT76sx\nDtzq7vOjuNe5+zR3fzu6z6eBRwm/zBsyMHpe/wbOAE5x99XZ5wBbuPuvo894EpgKjMiIszlVWlIG\nlBwkSQ4Mc/dtal7Af5PxC8vMzo4aXGuSRx9g24xr1NczZ0/gRODn7r5uEzG8m/H+c6BT9H5UFMfM\nqDrmvCieK6NqmNVmdkMT7jXbRnGb2RAze97MPozu83g2vs9sz0fPbHt3/5q7P1HPOTtnfw7wr2g/\n5LG9R0qPkoMUmszE0INQNfR9oEuUPOax8V+79f2Cmw+cT6ia2as5Qbj7e+5+gbt3A0YCN5jZ7u7+\nf+7eKXr9dwP//FOgY8Z21/o+ouaNmXUA/g5cBewQ3ec0Wv5X/XKge1aDdw9gWXYMItmUHKSQbUH4\nBfYB0Cb6671PLv/Q3ScBVwKPmVmdhtoGZCamb5rZLtHmx1Ec1Tle5yXgCDPrbmZbAVc09lmEtpL2\nhPusNrMhwOAcP6sxMwiloVFm1s7MUsAJhCo3CO0cuT4bKTNKDlJoNnQJdffXgGuB5wjVP32AZ+s7\nt4F/fxvwc+CJqBSSfX6D/5bQAP28ma0GHgAudvfFOcSOuz8G3Am8DMwiNArX91lE568GLgbuAj4i\ntAk01ntoU91ma+L4ilC9NgR4H7ge+La7vxGdNxHYN6qyu3cT9yZlxuJc7MfMKoHfARXABHcfm3V8\nG+Bmwl8vXwDnu/ursQUkEhMz+xA40t1fTjoWkXyIreRgZhWEv1QqgX2BEWbWO+u0K4EX3f0/gLOB\n6+KKRyQuZjaY8P/SwqRjEcmXOKuVBgCL3H2xu68l1HMOyzqnN6H7IO7+OqHr4PYxxiSSV9HguD8B\n/+nua5KORyRf4hwE142Nu9EtZeP+6QBzgVOBZ81sAKEnxS6E+lGRgufuw5OOQSQOcZYccmnM+DWw\ntZnNAS4E5gDrY4xJRERyEGfJYRnQPWO7O6H0sEHUS+P8mm0zext4K/tCZqb+2CIizeDuzRovE2fJ\nYTbQK5p0rD1hiP/kzBPMbKvoGGb2n8BT7v5pfRdzd73cGTNmTOIxFMpLz0LPQs+i8VdLxFZycPd1\nZnYh8AihK+tEd59vZiOj4+MJvZhujUoG84DvxBWPiIjkLtZZWd39IeChrH3jM94/B+wdZwwiItJ0\nGiFdZFKpVNIhFAw9i1p6FrX0LPIj1hHS+WJmXgxxiogUEjPDC7BBuuQ8+yyszp4xX0SkBCk55OiT\nT2DoULj77qQjERGJn5JDjsaPB3d4/vmkIxERiZ+SQw6+/BKuuw5+8xslBxEpD0oOObj9dth/fzj7\nbHjrLVi1KumIRETipeSwCdXVcPXVMGoUtG8P/frBrFlJRyUiEi8lh02YPBk6dYIjjwzbAweqaklE\nSp+SQyPcYezYUGqoWaJdyUFEyoGSQyOefRY++ABOPbV236GHhuSgMXkiUsqUHBpx1VVw6aVQUVG7\nb+edYfPN4c03k4tLRCRuSg4NmDcPZs8OPZSyqWpJREqdkkMDrrkGLroolBKyHXooPPdc68ckItJa\nYp2yu1gtWRJ6KTVUdTRwYBj7ICJSqjQraz1+/OPw89pr6z/+xRew7bbw/vvQsWOrhSUi0iQFOyur\nmVWa2QIzW2hmo+s5vp2ZPWxmL5nZPDM7N854crFyJdxyC1xyScPnbLZZGDE9e3brxSUi0ppiSw5m\nVgFcD1QSlgMdYWa9s067EJjj7n2BFHCtmSVa1fXHP8JJJ0H37o2fp0ZpESllcZYcBgCL3H2xu68F\nJgHDss75N9A5et8Z+NDd18UYU6PWrIHf/x4uu2zT5yo5iEgpizM5dAOWZGwvjfZlugnYz8yWA3OB\nH8QYzybddhscfDDst9+mz63psVQETTYiIk0WZxVOLr82rwRecveUme0B/MPM/sPd66y3VlVVteF9\nKpWKZZ3YqVNhdJ2Wkfrtumv4+c470KNH3kMREWmydDpNOp3Oy7Vi661kZgOBKnevjLavAKrdfWzG\nOdOA/3X36dH248Bod5+dda1W6a1UXR3mULIc2/ZPOQWGD4czzog3LhGR5ijU3kqzgV5m1tPM2gNn\nAJOzzlkAHANgZjsCewNvxRhTo9q0yT0xQGh30GA4ESlFsSWHqGH5QuAR4DXgTnefb2YjzWxkdNr/\nAQeZ2VzgMWCUu38UV0z5VjMJn4hIqdEguBb47DPYYQf46CPo0CHpaERENlao1Uolb4stYO+94cUX\nk45ERCS/lBxaSFVLIlKKlBxaSIPhRKQUKTm0kJKDiJQiJYcW2nPP0DC9fHnSkYiI5I+SQwuZqfQg\nIqVHySEPlBxEpNQoOeSBlg0VkVKjQXB5sGoV7LxzWCioXbukoxERCTQILmGdO8Nuu8HLLycdiYhI\nfig55Ikm4RORUqLkkCcaKS0ipUTJIU/UY0lESomSQ57ssw988AGsWJF0JCIiLafkkCdt2sAhh8CM\nGUlHIiLSckoOeaSqJREpFbEmBzOrNLMFZrbQzEbXc/xSM5sTvV4xs3VmtnWcMcVJPZZEpFTENgjO\nzCqA1wlrRC8DZgEj3H1+A+efAFzi7sfUc6ygB8HV+Ogj6NkzDIarqEg6GhEpd4U6CG4AsMjdF7v7\nWmASMKyR888E7ogxnth16RJGSr/6atKRiIi0TJzJoRuwJGN7abSvDjPrCBwH/D3GeFqFqpZEpBS0\njfHaTakHOhF41t0/buiEqqqqDe9TqRSpVKrZgcWpZjDcyJFJRyIi5SadTpNOp/NyrTjbHAYCVe5e\nGW1fAVS7+9h6zr0PuNPdJzVwraJocwCYOxeGD4f59basiIi0npa0OcSZHNoSGqSPBpYDM6mnQdrM\ntgLeAnZx9zUNXKtoksO6daHtYfHi8FNEJCkF2SDt7uuAC4FHgNcIJYP5ZjbSzDIrXU4GHmkoMRSb\ntm3hoINg5sykIxERaT6t5xCDK6+E9u0ho5lERKTVFWTJoZxppLSIFDuVHGKwYgXsvTd8+GGYc0lE\nJAkqORSYHXYIjdGvv550JCIizaPkEBNVLYlIMVNyiIlGSotIMVNyiImWDRWRYqYG6Zh89VVod/j3\nv6FTp6SjEZFypAbpAtS+PfTtq8FwIlKclBxipKolESlWSg4xUo8lESlWanOI0bJloWppxQqwZtX6\niYg0n9ocClS3brD55vDmm0lHIiLSNEoOMVPVkogUIyWHmGkwnIgUIyWHmKnkICLFSA3SMfviC9h2\nW3j/fejYMeloRKScFGyDtJlVmtkCM1toZqMbOCdlZnPMbJ6ZpeOMJwmbbQZ9+sDs2UlHIiKSu9iS\ng5lVANcDlcC+wAgz6511ztbAOOBEd+8DnBZXPElS1ZKIFJs4Sw4DgEXuvtjd1wKTgGFZ55wJ/N3d\nlwK4+wcxxpMYjZQWkWITZ3LoBizJ2F4a7cvUC+hiZk+a2Wwz+3aM8SRm4ED45z9h3bqkIxERyU3b\nGK+dSwtyO6A/cDTQEXjOzJ5394XZJ1ZVVW14n0qlSKVS+YmyFfToAf37w0UXwQ03aLS0iMQjnU6T\nTqfzcq3YeiuZ2UCgyt0ro+0rgGp3H5txzmhgc3evirYnAA+7+z1Z1yra3ko1Vq2Cww+Hb38bLrss\n6WhEpBwUam+l2UAvM+tpZu2BM4DJWec8ABxuZhVm1hE4BHgtxpgS07kzPPgg/P73cM89mz5fRCRJ\nsVUrufs6M7sQeASoACa6+3wzGxkdH+/uC8zsYeBloBq4yd1LMjkAdO8OkyfD4MGwyy6hLUJEpBBp\nEFwCpk2D73wHpk+H3XdPOhoRKVWFWq0kDTj+ePif/wk/P/oo6WhEROpSySFBl14Ks2bBo49Chw5J\nRyMipaYlJQclhwRVV8M3vxnmXLrtNnVxFZH8UrVSkWrTBv7yF3jjDcgYxiEikrg4B8FJDjp2DD2Y\nDj00NE6fc07SEYmIKDkUhB13DGMgUqnQ3fWoo5KOSETKnaqVCkTv3nDnnTBiBLxWsiM9RKRYKDkU\nkFQKrr4ahg6F995LOhoRKWdKDgXm7LNDu8OJJ8LnnycdjYiUK3VlLUDuIUF8+incfTdUVCQdkYgU\nI3VlLTFmMGECrFwJo0YlHY2IlCMlhwLVvj3ce2+Yh2ncuKSjEZFyo66sBWybbUIX18MPDwsGnXBC\n0hGJSLlQm0MRmDEjJIZHHgkryomI5EJtDiXukEPgT3+CYcNgyZJNny8i0lKxJgczqzSzBWa2MFoS\nNPt4ysw+MbM50ev/xRlPMfvGN+CSS8IYiFWrko5GREpdnGtIVwCvA8cAy4BZwAh3n59xTgr4kbuf\ntIlrlXW1Ug13+P734c03YepUaNcu6YhEpJAVarXSAGCRuy9297XAJGBYPedpouocmYU1qNu2DUlC\n+VJE4hJncugGZNaQL432ZXLga2Y218ymmdm+McZTEtq2hUmTwiJBV12VdDQiUqri7Mqay9+1LwLd\n3f1zMxsC3A/sFWNMJaFTp1CtdOih0LWrpvkWkfyLMzksA7pnbHcnlB42cPfVGe8fMrMbzKyLu9dZ\nWbkqYzWcVCpFKpXKd7xFpVu30LX1uOPgk0/g4ouTjkhEkpZOp0mn03m5VpwN0m0JDdJHA8uBmdRt\nkN4RWOHubmYDgLvcvWc911KDdAP+9S849lg480wYM0ZLjYpIrZY0SMdWcnD3dWZ2IfAIUAFMdPf5\nZjYyOj4eOA34npmtAz4HhscVT6nq0QOeeQYqK+HDD+G668LyoyIiLaER0iXik0/CNN/du8Ott6qb\nq4i0rOSQU3Iws77AIEIj8zPuPrc5H9ZcSg65WbMGTj89dHG9666wPrWIlK9YxzmY2Q+A24HtgR2B\n281MzZ8FaPPNw0yu22wTGqo//jjpiESkWG2y5GBmrwAD3f2zaHsL4Hl3378V4quJQSWHJqiuDlNt\nPP106NG0445JRyQiSWiNEdLVDbyXAtSmTWiYPvXUMN334sVJRyQixSaX3kq3ADPM7F7CVBcnAzfH\nGpW0mBn89KfQpQsMGgQPPwz77Zd0VCJSLHJtkO7Pxg3Sc+IOLOvzVa3UAn/9K/z4x/DAA2H6bxEp\nD7H0VjKzgcB4YE/gZeA77v5as6NsASWHlnvwQTjvvJAojj026WhEpDXE1eYwDrgU2Bb4DfDb5nyA\nFIahQ+Hvf4ezzoJ77kk6GhEpdI0lhzbu/g93/8Ld7wZ2aK2gJB6DBsGjj4Z5mG66KeloRKSQNdYg\nvZWZnUrteguZ2+7u98YeneRd377w1FMweDB89BGMrrM+n4hI420Ot7LxtNuWue3u58Ua2caxqM0h\nz5YtCwli6FAYO1YT9omUorgapL9GGOyW+LgGJYd4fPghHH887L8/jB8PFRVJRyQi+RRXg/TZwAtm\nNsnMzjWzrs0LTwrVttvC44+HQXKnnw5ffpl0RCJSKHKZPqM3MAQYDGwNPAE8DEx39/WxR4hKDnH7\n8svQi+njj+G++8JKcyJS/GKflTXjgzoCRxKSxaHufmBzPrSplBzit349jBwJr7wC06aFUoWIFLe4\n2hw2B/6L2kFwE919XbOjbAElh9bhDpdfDlOmhC6vu+ySdEQi0hJxtTn8GTgQeAU4Hri2GYFVmtkC\nM1toZg12mjSzg81sXdRVVhJiFnounXtuGBOxcGHSEYlIUhob59C7ZlpuM5sAzGrKhc2sArgeOAZY\nBswys8mZa0hnnDeW0I6hDpUFYNSoMGHf178ept3o1y/piESktTVWcthQhdTM6qQBwCJ3X+zua4FJ\nwLB6zrsIuAd4vxmfITH57nfhD38IiwY980zS0YhIa2ssORxgZqtrXsD+Gdurcrh2N2BJxvbSaN8G\nZtaNkDD+GO1Sw0IB+cY34G9/Cz+nTk06GhFpTQ0mB3evcPdOGa+2Ge8753DtXH7R/w64PGptNlSt\nVHCOOSY0UH/3uzBhQtLRiEhryWWxn+ZaBnTP2O5OKD1kOhCYZGHuhu2AIWa21t0nZ1+sqqpqw/tU\nKkUqlcpzuNKQQw4JS46ecAIsWBAarTWaWqTwpNNp0ul0Xq7VpHEOTbqwWVvgdeBoYDkwExiR3SCd\ncf4twJT6JvRTV9bC8NFHcNppsMUWobpJg+VECltrrCHdZFEj9oXAI8BrwJ3uPt/MRprZyLg+V+LT\npUtYbrRr17A29TvvJB2RiMQltpJDPqnkUFjc4be/hWuugXvvhYEDk45IROrTatNnJEXJoTBNmQLn\nnx+6vA4fnnQ0IpJNyUESM3cunHRSWJ96zBitCyFSSJQcJFHvvgsnnwy77QY33wybb550RCICBdog\nLeWja1d48snw/sgjQ7IQkeKm5CB5sfnmoXvrkCFhXMTcuUlHJCItoWolybtJk+Cii0IV04knJh2N\nSPlqSbVSnCOkpUwNHx7aH049Fd54A370IzVUixQblRwkNu+8E0oOBx8M48ZBhw5JRyRSXtQgLQVp\n113h2Wdh5cowonrx4qQjEpFcKTlIrDp1gnvuCVVNhxwCDz2UdEQikgtVK0mrefppGDEiTP/9059q\nZleRuGkQnBSNd9+FM86AzTaDv/4Vttsu6YhESpfaHKRodO0Kjz8OffvCgQfCjBlJRyQi9VFykFbX\ntm1YMOi660JvpnHjwkyvIlI4VK0kiVq0KKxR3acP3HhjWEhIRPJD1UpStPbcE557Dtq1gwEDwjKk\nIpK8WJODmVWa2QIzW2hmo+s5PszM5prZHDN7wcyOijMeKUwdO8Itt8All8CgQXD33UlHJCJxriFd\nQVhD+hhgGTCLrDWkzWwLd/8ser8/cJ+771nPtVStVCZeeCGsU33yyXDVVaFEISLNU6jVSgOARe6+\n2N3XApOAYZkn1CSGyJbABzHGI0XgwANDgli4MEz/vWxZ0hGJlKc4k0M3YEnG9tJo30bM7GQzmw88\nBFwcYzxSJLp0gcmTw/TfBx0ETzyRdEQi5SfOWVlzqgdy9/uB+81sEPAXYO/6zquqqtrwPpVKkUql\nWh6hFKw2beAnPwlTbnzrW3D00fDLX0KPHklHJlK40uk06XQ6L9eKs81hIFDl7pXR9hVAtbuPbeTf\nvAkMcPcPs/arzaGMrVoF11wTxkOcdx5ceWUoXYhI4wq1zWE20MvMeppZe+AMYHLmCWa2h1mY6d/M\n+gNkJwaRzp3h5z+HefNg9WrYe2+4+mr44oukIxMpXbElB3dfB1wIPAK8Btzp7vPNbKSZjYxO+wbw\nipnNAa4DhscVjxS/nXaC8ePDBH7Tp4ckcdttsH590pGJlB6NkJai9eyzcNllsGZNmI5j8GCtOCeS\nSbOyStlyh/vugyuugO7dw9iI/v2TjkqkMBRqm4NI7MzCWtXz5oU5moYOhbPOgrffTjoykeKm5CAl\noV07+N73wuC5vfYK4yN+9CP4UN0bRJpFyUFKypZbwpgx8OqroTfTPvuE9og1a5KOTKS4KDlISera\nFW64ITRaz5wZejbdeqt6NonkSg3SUhb++U8YNSoMqBs7Fior1bNJSp96K4nkwD3M2TR6NOy8c+jZ\ndNBBSUclEh/1VhLJgRkMGxZ6Ng0fDiedBGeeCUuWbPrfipQbJQcpO23bwgUXhJ5NvXpBv35heo7P\nP086MpHCoeQgZWuLLeBnPwvrR7z6KvTuDXfdFaqfRMqd2hxEIk8/DT/4AXTqBNddF0oUIsVMbQ4i\neXDEETB7dlg/YsiQUPX0/vtJRyWSDCUHkQwVFSEpLFgQBtTtuy/89rfw1VdJRybSupQcROqx9dbw\nm9/AM8/Ao4/CAQfAQw8lHZVI61Gbg0gOpk2DH/4Q9twzlCT22ivpiEQ2TW0OIjE7/nh45RU46ig4\n7DC49FL45JOkoxKJT+zJwcwqzWyBmS00s9H1HD/LzOaa2ctmNt3MDog7JpHmaN8efvzjMIjuk0/C\npH4TJmi+JilNsVYrmVkF8DpwDLAMmAWMcPf5GeccCrzm7p+YWSVQ5e4Ds66jaiUpOC++GLq+fvZZ\n6Po6aFDSEYlsrJCrlQYAi9x9sbuvBSYBwzJPcPfn3L2mgD4D2CXmmETyon//MDZi1KiwwNDw4fDO\nO0lHJZIfcSeHbkDmzDVLo30N+Q4wLdaIRPLILCSFBQtCNVO/flBVBR9/nHRkIi3TNubr51wXZGZH\nAucDh9V3vKqqasP7VCpFKpVqYWgi+dOxY0gK558PP/kJ7LYbnHIKjBwJAwZoenBpHel0mnQ6nZdr\nxd3mMJDQhlAZbV8BVLv72KzzDgDuBSrdfVE911GbgxSVFSvC4kI33hgG040cGaqeOndOOjIpJwW7\nnoOZtSU0SB8NLAdmUrdBelfgCeBb7v58A9dRcpCiVF0NTzwB48fDY4/BaaeFRKF1JKQ1FGxyADCz\nIcDvgApgorv/ysxGArj7eDObAJwC1DTlrXX3AVnXUHKQovfuu3DzzXDTTdClS0gSI0aEif5E4lDQ\nySEflByklFRXhyk5xo+Hp56C008PiUKzwEq+KTmIFKnly2HixFCa6NoV/uu/4IwzwloTIi2l5CBS\n5Navh4cfDqWJ6dNDddPIkbD//klHJsWskAfBiUgOKipg6FCYPBleegm22y6sKfG1r8Gf/wxr1iQd\noZQblRxECtS6dWE22PHjYcaM0BV25MiwxoRILlRyEClBbdvCSSfBgw+Gda47d4ZjjglzON1+O3zx\nRdIRSilTyUGkiKxdC1OnhtLECy/At78dVq7bZ5+kI5NCpJKDSJlo1y5My/HwwzBzJmy2GaRS4XXH\nHfDll0lHKKVCJQeRIvfVV/DAA6E0MXduaNg+4QQYPFjTdZQ7dWUVESBMGT5lSqh6mj4dDjkETjwx\nJIvdd086OmltSg4iUsenn4b5nKZMCY3a224bksSJJ8LAgaHBW0qbkoOINKq6GmbPDiWKqVNDCaOy\nMiSK446DrbdOOkKJg5KDiDTJ0qWhNDFlSljN7qCDaksVvXolHZ3ki5KDiDTb55/D44/Xliq23LI2\nURx2WOghJcVJyUFE8sId5sypbdR+881Q7XTCCWE6jy5dko5QmkLJQURisXx5mMJjyhR48kno27e2\n99M++2j500Kn5CAisVuzJiSIqVNDsujQobb6adAgaN8+6QglW0EnBzOrpHYluAn1rB+9D3AL0A/4\nibtfW881lBxECog7vPxybfXT66/DsceGZHH88WFWWUlewSYHM6sgrCF9DLAMmEXdNaS3B3oAJwMr\nlRxEis9779VWPz3+OPTpU1uq2G8/VT8lpZCTw6HAGHevjLYvB3D3X9dz7hjgUyUHkeL25ZeQTtdW\nP5mFRHHCCWEOqA4dko6wfBTyxHvdgCUZ20ujfSJSojp0CD2c/vAHePvtkCC6dYOf/xx22AFOPRXG\njYPnnw/daKUwxT2APm9/7ldVVW14n0qlSKVS+bq0iMTELFQx9ekDl18O778PDz0EzzwDN98M8+fD\nHntA//5w4IHhZ9++YayFNF06nSadTuflWnFXKw0EqjKqla4AqrMbpaNjqlYSKTNffgmvvhrWpnjx\nxfCaNw923XXjhNGvH2y1VdLRFp9CbnNoS2iQPhpYDswkq0E649wqYLWSg0h5W7s2lChqEsYLL4Se\nUTvttHHC6N9fg/I2pWCTA4CZDaG2K+tEd/+VmY0EcPfxZtaV0IupM1ANrAb2dfdPM66h5CBSxtav\nD91lMxPGSy+FmWYzk8WBB8L22ycdbeEo6OSQD0oOIpKtuhoWLqytjqpJHJ061SaMmp877ZR0tMlQ\nchARIQzOe+utjZPFCy+E0dvZCWOXXUp//IWSg4hIA9zD+hXZCcN94+qo/v2hZ8/SShhKDiIiTeAe\nJhXMThhr1tRt9N5jD2gT94iwmCg5iIjkwbvv1rZh1CSMjz8OXWkzq6V69YKKiqSj3TQlBxGRmHzw\nQd1G7xUrwmC9zGqpffYpvHW5lRxERFrRypVhUaTMhLFsGey//8bVUvvum+xU5koOIiIJW7UqjL3I\nTBhvvx1mpc1MGPvv33qTDyo5iIgUoM8+g7lzN54e5Nxz4Yc/bJ3PV3IQEZE6CnnKbhERKUJKDiIi\nUoeSg4iI1KHkICIidSg5iIhIHUoOIiJSR6zJwcwqzWyBmS00s9ENnPP76PhcM+sXZzwiIpKb2JKD\nmVUA1wOrDpy/AAAIUElEQVSVwL7ACDPrnXXO8cCe7t4LuAD4Y1zxlIp8LR5eCvQsaulZ1NKzyI84\nSw4DgEXuvtjd1wKTgGFZ55wE/BnA3WcAW5vZjjHGVPT0xa+lZ1FLz6KWnkV+xJkcugFLMraXRvs2\ndc4uMcYkIiI5iDM55DrfRfbQbs2TISKSsNjmVjKzgUCVu1dG21cA1e4+NuOcPwFpd58UbS8Avu7u\n72VdSwlDRKQZmju3UpxLU8wGeplZT2A5cAYwIuucycCFwKQomXycnRig+TcnIiLNE1tycPd1ZnYh\n8AhQAUx09/lmNjI6Pt7dp5nZ8Wa2CPgMOC+ueEREJHdFMWW3iIi0roIeIZ3LILpSZmaLzexlM5tj\nZjOjfV3M7B9m9oaZPWpmWycdZxzM7GYze8/MXsnY1+C9m9kV0fdkgZkNTibqeDTwLKrMbGn03Zhj\nZkMyjpXys+huZk+a2atmNs/MLo72l913o5FnkZ/vhrsX5ItQFbUI6Am0A14CeicdVys/g7eBLln7\nrgJGRe9HA79OOs6Y7n0Q0A94ZVP3Thhk+VL0PekZfW/aJH0PMT+LMcCP6jm31J9FV6Bv9H5L4HWg\ndzl+Nxp5Fnn5bhRyySGXQXTlILsxfsPAwejnya0bTutw92eAlVm7G7r3YcAd7r7W3RcTvvQDWiPO\n1tDAs4C63w0o/Wfxrru/FL3/FJhPGC9Vdt+NRp4F5OG7UcjJIZdBdKXOgcfMbLaZ/We0b0ev7dH1\nHlBOI8obuvedCd+PGuXyXbkompNsYkY1Stk8i6gnZD9gBmX+3ch4Fs9Hu1r83Sjk5KCWcjjM3fsB\nQ4Dvm9mgzIMeyopl+ZxyuPdSfy5/BHYD+gL/Bq5t5NySexZmtiXwd+AH7r4681i5fTeiZ3EP4Vl8\nSp6+G4WcHJYB3TO2u7Nx1it57v7v6Of7wH2EIuB7ZtYVwMx2AlYkF2Gra+jes78ru0T7Spa7r/AI\nMIHa6oGSfxZm1o6QGP7i7vdHu8vyu5HxLG6veRb5+m4UcnLYMIjOzNoTBtFNTjimVmNmHc2sU/R+\nC2Aw8ArhGZwTnXYOcH/9VyhJDd37ZGC4mbU3s92AXsDMBOJrNdEvwBqnEL4bUOLPwswMmAi85u6/\nyzhUdt+Nhp5F3r4bSbe4b6I1fgihBX4RcEXS8bTyve9G6FnwEjCv5v6BLsBjwBvAo8DWScca0/3f\nQRhZ/xWh7em8xu4duDL6niwAjks6/pifxfnAbcDLwFzCL8Idy+RZHA5UR/9fzIleleX43WjgWQzJ\n13dDg+BERKSOQq5WEhGRhCg5iIhIHUoOIiJSh5KDiIjUoeQgIiJ1KDmIiEgdSg6SqPqmo472DzSz\nm8zsg5rBgBnH7jezb0bve5jZORnHrJ59Z0XzzLxsZtPN7ICMYx3M7Ckza2NmKTOb0ox7GNPYvoam\nVs661xub+JnnmFmPHPZdZGbzo88dG+37DzOb2JTPk/Kj5CBJu4UwiClbJfAQYSXBU2p2mtlWwGHA\nVAtrkB8O9DCzCWa2M2Femex9bwFHuPsBwC+AzF/EZwFT3b26qYFHSedSYDMzuyzarrMPWAv80N33\nAwYS5snqnXGpIdG95vKZO5vZTYRpEAaZ2Z/q2xedeyRhttID3L0PcA2Au88F9jCzHZp6z1JGkh7l\np5dehLnlX8na9wzQCTgBmJax/1zCnDoQ1vx4kLDuxfbRvjYZ+7ar57O2AZZmbP8D2Ct6nwKmRO8P\nBl4kjFTfPjpvHnATsJhonQ3CuujrgNMzrllnX1YM9wNH13Ov50bHHo3ivxC4NIrjOWCb6Pwdo+OT\nqV3Nsb59dwFHNRDDaOD7Sf+316twXyo5SMExs+2AtR5m23wU6G9m20SHhxOmkwD4A/BX4Gbgf6M5\nZa7P2Pd/WfPMAHwHmBZ9TgXQx93fyPr8rxFKICe5+9uExVMe8/DX9z3AruE0G0GY8vhqQkllRH37\nsq7dk9ppprPvFWA/QknpYOB/gVXu3p+QHM6O7ucXhDl17gLG1bcvulYv4Agze97M0mZ2UEYoM4Ej\n6jx8kRpJZye99CKr5ACcCYzO2L4RuADYDngXqMg41gM4J+t6dfZF+48EXmPjv8DfyDieIsw78wrQ\nNWP/HKBHxvaHZKzQB4yp57Pq27clYULJk+u7V0LJ4caMY/8Cdorenwf8NuPYOZkx1bcvuo/rovcH\nA29lHNsHmJH0f3u9CvelkoMUokrg4YztOwglhm8A97v7+poD7v4vd/9z5j+ub1/UCH0ToTSQuapa\n5opZTpj/fg3QPyum+lbWqvm8n21qX31TK0cy79WBLzOOVWdsVwNtM67/Z3f/V9ZnZu9bCtwbHZsF\nVJvZthn3o4nVpEFKDlJQommID/DQaFojDewFfJ/aKqWmXHNXwi/Jb7n7ooxDHxD+mt9wKvAxoZ3j\nV2b29Wj/dOD06FqDCe0WTfn8hqZWzr7XBhPQJo415H7gqOiz9gLau/uH0bGdCCUTkXopOUiizOwO\n4J/A3ma2BLiM0AC7gbs7cDehKuepZnzM/xB+of/RzOaY2czouuuBeWa2d81HRR+3gpAgxpnZwcDP\ngMFRd9vTCFVbq7M/pBGHAd8Cjow+f46ZDQEOJFRZbbhVNv5rPvt9U//SvxnYPYr7DuDsjGMDgKeb\neD0pI5qyWwqKmf0EWOjud7XS551LmO9+bCPntAfWu/t6MzsUGOehkbiln92q95r12WlCb6pyWklQ\nmkDJQcpa9Iv/MeDr3sD/DGa2J6EXUBvCgjvfc/cXWi/K/IraXy529+8mHYsULiUHERGpQ20OIiJS\nh5KDiIjUoeQgIiJ1KDmIiEgdSg4iIlKHkoOIiNTx/wEoXiEQ59McdAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f21b6346e50>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "from matplotlib.pyplot import plot,title,show,xlabel,ylabel\n",
    "\n",
    "\n",
    "#(a) \n",
    "  #For n = 1,\n",
    "  #(P/P0)/V= (P/P0)V1+ 1/B2V1\n",
    "  #where V is the gas phase volume equivalent to the amount adsorbed\n",
    "  #x = (P/Po)\n",
    "x = [0.05, 0.10 ,0.15 ,0.20 ,0.25 ,0.30, 0.35, 0.40]\n",
    "  #y = (P/Po)/V in kg/m**3\n",
    "y = [0.76 ,1.35, 1.85, 2.27, 2.66, 2.94 ,3.21, 3.42]\n",
    "  #y1 = (P/Po)/(V*(1-P/Po)) in kg/m**3\n",
    "y1 = [0.80, 1.50, 2.18 ,2.88, 3.55, 4.20 ,4.94, 5.73]\n",
    "plot(x,y)\n",
    "plot(x,y1)\n",
    "xlabel(\"(P/Po)\")\n",
    "ylabel(\"(P/Po)/V or (P/Po)/(V*(1-P/Po))\")\n",
    "show()\n",
    "  \n",
    "  #The data, which are plotted as (P /P 0)/V against P/P0  may be seen    to conform to a straight line only at low values of P/P0, suggestin    g that more than one layer of molecules isadsorbed.\n",
    "Slope = 12.56\n",
    "V1 = 1/12.56\n",
    "  #The surface area occupied by this absorbed volume \n",
    "S = V1*6.02*10**(26)*0.162*10**(-18)/24\n",
    "print\"\\n S = %d m**2/kg\"%(S)\n",
    "  \n",
    "  \n",
    "#(b) \n",
    "  #P/P0/V*(1 − P/P0)= 1/V1*B2+(B2−1)/(V1*B2)*(P /P 0)\n",
    "  #y2 = (P/Po)\n",
    "  #x2 = 1V**2\n",
    "y2 = [0.05, 0.10, 0.15, 0.20, 0.25, 0.30 ,0.35, 0.40, 0.50 ,0.60 ,0.87, 0.80]\n",
    "x2 = [230 ,183 ,152 ,129, 113, 96, 84 ,73 ,53 ,37 ,26 ,20]\n",
    "\n",
    "plot(x2,y2)\n",
    "title(\"Harkins-Jura Plot\")\n",
    "xlabel(\"1/V**2(kg**2/m**6)\")\n",
    "ylabel(\"P/Po\")\n",
    "show()\n",
    "  \n",
    "  \n",
    "  \n",
    "  "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Page 1015 Example 17.3"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      " m=0.666667\n",
      "\n",
      " z = 3.18 m\n",
      "\n",
      " z = 0.97 m\n",
      "\n",
      " For case (c)\n",
      "\n",
      " z = 0.56 m\n",
      "\n",
      " saturation = 20.5 per cent\n"
     ]
    }
   ],
   "source": [
    "from math import pi\n",
    "\n",
    "#For case (a):\n",
    "  #Cs = 10C which represents a linear isotherm.\n",
    "  #All concentrations move at the same velocity.If z0=0 at t = 0 for a    ll concentrations, the adsorption wave propagates as a step change     from the inlet to the outlet concentration.\n",
    "  #f (Cs ) = 10\n",
    "  #u = 1×10**(−4)/[(π/4)*(0.15**(2)*ε)] m/s\n",
    "  #where ε is the intergranular voidage = 0.4\n",
    "e = 0.4\n",
    "m = e/(1-e)\n",
    "print\"\\n m=%f\"%(m)\n",
    "t = 3600#              #time is in secs\n",
    "z = ((4*10**(-4))/(pi*(0.15**2)*0.4))*(3600/(1+10*(0.6/0.4)))\n",
    "print\"\\n z = %.2f m\"%(z)\n",
    "  \n",
    "  \n",
    "#It may be noted that, when the adsorption wave begins to emerge from the bed, the bed is saturated in equilibrium with the inlet concentration.\n",
    "#Hence: uAεtC0 = zA[εC0 + (1 − ε)Cs ]\n",
    "\n",
    "#For case (b):\n",
    "#Cs = 3C0.3 which represents a favourable isotherm.\n",
    "#As C increases, f(C) decreases and points of higher concentrations are predicted to move a greater distance in a given time than lower concentrations. It is not possible for points of higher concentrationsto overtake lower concentrations, and if z0=0 for all concentrations,the adsorption wave will propagate as a step change similar to case a.\n",
    "\n",
    "#Hence: z = ut/[1+(1/m)(Cs/Co)]\n",
    "Co = 0.003#              #in kmol/m**3\n",
    "z = ((4*10**(-4)*5/((pi)*0.4*0.15**2)))*(3600/(1+(0.6/0.4)*(3/Co**0.7)))\n",
    "print\"\\n z = %.2f m\"%(z)\n",
    "\n",
    "#For case (c):\n",
    "#Cs = (10**(4))*(C**2) which represents an unfavourable isotherm.\n",
    "#f(C)=2* 10**4* C.\n",
    "#As C increases, f (C) increases such that, in a given time, z for lower concentrations is greater than for higher concentrations. Following the progress of the breakpoint concentration,C = 0.003 kmol/m**3,then:\n",
    "#f(0.003) = 60\n",
    "z1 = (4*10**(-4)/(pi*0.15**(2)*0.4))*(3600/(1+(0.60/0.40)*60))\n",
    "print\"\\n For case (c)\"\n",
    "print\"\\n z = %.2f m\"%(z1)\n",
    "\n",
    "\n",
    "#At breakpoint,the bed is far from saturated and:\n",
    "saturation = 100*((4*10**(-4)*3600)/(pi*(0.15**(2)*0.4)))/(0.55*(1+(0.6/0.4)*(9/0.03)))\n",
    "print\"\\n saturation = %.1f per cent\"%(saturation)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Page 1031 Example 17.5"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      " m = 0.89\n",
      "\n",
      " Gs = 4.16*10**(-3) kg/m**2.sec \n",
      "\n",
      " Gs1 = 4.10*10**(-3) kg/m**2.sec\n",
      "\n",
      " ta = 8.3 hrs\n",
      "\n",
      "I= 0.0\n",
      "\n",
      " I=0.000106\n",
      "\n",
      " I=0.000745\n",
      "\n",
      " I=0.001570\n",
      "\n",
      " I=0.002556\n",
      "\n",
      " I=0.003675\n",
      "\n",
      " I=0.004906\n",
      "\n",
      " I=0.005568\n",
      "\n",
      "time = 0.0 h\n",
      "\n",
      "time = 1.1 h\n",
      "\n",
      "time = 3.0 h\n",
      "\n",
      "time = 3.9 h\n",
      "\n",
      "time = 4.6 h\n",
      "\n",
      "time = 5.6 h\n",
      "\n",
      "time = 7.1 h\n",
      "\n",
      "time = 8.3 h\n"
     ]
    }
   ],
   "source": [
    "from sympy import symbols,solve\n",
    "\n",
    "#let x = [f(C)]mean\n",
    "x = (8.4*10**(-2)*1266)/(2.67*10**(-3)*1.186)\n",
    " \n",
    " #e is porosity\n",
    "e = 0.47\n",
    "m = e/(1-e)\n",
    "print\"\\n m = %.2f\"%(m)\n",
    "   \n",
    "#The velocity uc with which the adsorption wave moves through the column may be obtained from equation 17.79\n",
    "uc = 6.2*10**(-6)\n",
    "density = 1266#               #density is in kg/m**3\n",
    "Gs = uc*(1-e)*density\n",
    "print\"\\n Gs = %.2f*10**(-3) kg/m**2.sec \"%(Gs*10**3)\n",
    " \n",
    " \n",
    "#From overall balance:\n",
    "  #We have given eq: Gs1(0.084-0)=0.129(0.00267)\n",
    "  #On solving above eq :\n",
    "Gs1 = symbols('Gs1')\n",
    "Gs2 = solve(Gs1*(0.084-0)-0.129*(0.00267))[0]\n",
    "print\"\\n Gs1 = %.2f*10**(-3) kg/m**2.sec\"%(Gs2*1000)\n",
    "  \n",
    "#Part(b)\n",
    " #time ta for the adsorption zone to move its own length za is given by:\n",
    "za = 0.185\n",
    "ta =(za/uc)/3600\n",
    "print\"\\n ta = %.1f hrs\"%(ta)# \n",
    " \n",
    " #The time taken for a point at a distance z into the zone to emerge is given by:\n",
    " #t = (z1/za)ta\n",
    "yr = [0.0001, 0.0002, 0.0006, 0.0010 ,0.0014 ,0.0018, 0.0022, 0.0024]\n",
    "yre = [0.00005, 0.0001, 0.00032, 0.00062, 0.00100, 0.00133, 0.00204, 0.00230]\n",
    "I=0\n",
    "I1 = 0\n",
    "i=0\n",
    "print\"\\nI= %.1f\"%(I)\n",
    "Y=[]\n",
    "while i <=7:\n",
    "    y.append(1/(yr[(i)]-yre[(i)]))\n",
    "    if i>0:\n",
    "        I = I + (yr[(i)]-yr[(i-1)])*(y[(i)]+y[(i-1)])/2\n",
    "        print\"\\n I=%f\"%(I)\n",
    "    \n",
    "    i=i+1#   \n",
    "\n",
    "z_za = [0, 0.137, 0.362, 0.473, 0.560 ,0.674, 0.852, 1.000]\n",
    "t=[]\n",
    "for zza in z_za:t.append((zza)*ta)\n",
    "j=0\n",
    "while j<=7:\n",
    "    print\"\\ntime = %.1f h\"%(t[(j)])\n",
    "    j=j+1\n"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.9"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}