{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "###### The latest version of this IPython notebook is available at [http://github.com/jckantor/CBE20255](http://github.com/jckantor/CBE20255) for noncommercial use under terms of the [Creative Commons Attribution Noncommericial ShareAlike License](http://creativecommons.org/licenses/by-nc-sa/4.0/).\n", "\n", "J.C. Kantor (Kantor.1@nd.edu)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Antoine's Equation" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Initialize the IPython workspace with with default settings for plots." ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": true }, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "%matplotlib inline\n", "\n", "plt.rcParams.update({'figure.figsize':(10,6)})\n", "plt.rcParams.update({'font.size':12})\n", "plt.rcParams.update({'lines.linewidth':1.6})" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Introduction" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/jpeg": "/9j/4AAQSkZJRgABAQAAAQABAAD/2wCEABALDA4MChAODQ4SERATGCgaGBYWGDEjJR0oOjM9PDkz\nODdASFxOQERXRTc4UG1RV19iZ2hnPk1xeXBkeFxlZ2MBERISGBUYLxoaL2NCOEJjY2NjY2NjY2Nj\nY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY//AABEIAWgB4AMBIgACEQED\nEQH/xAAaAAEAAwEBAQAAAAAAAAAAAAAAAQQFAgMG/8QARxAAAQMCAwQECwUGBQUAAwAAAQACAwQR\nBRIhEzFBURQiYZEGFTJSU3FygZKhsiM0NVSxQnN0wcLRJDNi8PElQ4Ki4RZjZP/EABgBAQEBAQEA\nAAAAAAAAAAAAAAABAgME/8QAKBEBAQACAgICAgEDBQAAAAAAAAECERIxAyEiQVFhE4GxwQQyodHx\n/9oADAMBAAIRAxEAPwD7xVan75R+076SrN1VqfvlJ7TvpK4RpaRLhLqaGPVxmjxqKrZI7LUHJIwN\nvo1pOnFeVHWRYnjkVVQtfs2QuZM97C2+osPXoVcxIkV2GkBxAmINhe3VK8pqWooawVGHMDo5XDbw\nk2HtjtXql3jN96ZWa3FqSgkEczyZCL5GNzED1BeAx6kkgZLC2aUyOysY2M5ibX3HsVN+C1zn1gbW\ntY2ocSXiP7Q33DNwC9WUNeIKSxp2S0pswDMQW5cuum9SYePXZutCDEqeakfUOJhZGS2QSjKWHtU0\n9cyorZ6eNp+xa0l99DfVUjhb3YPUUrps1ROc75C3QuJ/TReuHUlTHWz1VY6PaSNbGGxg2Abf+6zc\ncNWxWVMaoYBiHR3RNjElQJM4JNsx3LVpcNOyo31E5kmpzma4NA0Itl9X9lUOD1Lo6+M1VoZtoYo2\njS7testLDop4KUMqpdrKXEkgGwB3ALeefx+NSR5OqZh4QMpQ77E0+dwy7jmtvWgqNPSytxOoq55G\nuztEcbWjyWi5171eBBXHP60oszDMor8VlNgdsGk9gY1aRIH/AAsKuw2td05sTmOiqZWvDLkE2tmB\nPIgLXim9y3WyrkGO0NTWNpYXve998rgw5TbkVpXGuo03r5+PC8R2wq46iGnky7MRNZmZGzs7V03w\ndy7aLxhPsJjeRlus51rXzX3cbLWXjw+qbrRxSqnooG1EUbZImG8w/aDeYVuORssTJGG7HgOaeYKx\negYpNSspKirhEGXI8sjOZw/ktmFjIYI42XyMaGi/ILGWMmP7HooVeurI6GmdPKHEN0s0XJPJZs+N\nVUDWTy4c+Knc5ou94z6m25THx5ZdG2nXMqJKZwpZdnMNWki4NuB7Cs3Dq2XFK95c0xxU8WzliPGQ\nnUeoAfNadXUdHpJpgwvMbC4NAN3cgsvBsPnpZ+kPkN6hpdOw+eTcfLRbw1wu/wChe2xFFHDGI4mN\nYwbmtFgFgYdi1LQNqKWYS5mVEnkxki2Yq/glU6oppto8vdHPIy55A/8A1c4PFJFPiO1je0OqnOZm\nG8HiE1rlMhTna7FJoX0VA6ECQONVI3I63Gw3m6166gp8QhENUwvjBvYOISetZBWU9M5ri6fMGkbg\nQL6qxnB57r7lMs8vWvQqSYVQSxhj6SJwDQ0Et1sN2q94KaGmblhjawWA0423LyxGubQUb6lzC8MI\nuBpvNlZDgQCs252e+hk4nhpfI+djampLyBsW1Gza3tXWBYfJRCd0tPDCZCLCNxcfeSfUtW6XTneP\nE0zPCCCWaiifEwvME7JS1u9wG+ypT4y6sraGGjZI2F8wD5HtLb2F8oW+XAb79ywW1gxOvwyUwvij\nvI9mb9qw36Lr4vc9zraVvrExyoZS4thkzmvcAZLtY0uJGXkrseIh2LzUL25S1jXRmx62+/douaul\nkmxahqWluzhD81xrqLaLGE45fL8f4WvGnbUVmLsrXUpp4o4iwF9s77nlyWsl0usZ5cr0qUUXS6xo\nSuJGNkjcx18rhY2Nl1dLqjCNY/wfa+CqEktPvp5N59g/37VewalfTUhfPbpFQ8yy23Bx1t7tyqbP\nx1XSGpid0KnJYxjhpI/cXW5BeWMCLDcNpYqaaSN8cwETnONhrc35i116bOesZ3e2W+sTGRLUYvh9\nPSyRsnZnkJcL5Ra17K1BizamqbFBBM+MAl85YWsHqvvVrosHTemZPt8mzzdl7rnjvx5e4vbypMNi\npaaWISSPdNcySk9ZxPHsXhTYFSQ080MpfUCZwc8ym9yNykYvHDNVR1uWHYO0INw5p1Hv7FoB4Lcw\nvYi+5S3yY3Z6U4aOlwoVM8LXtjcMzo26gW5BeHj+mIGSnrH33WgctQOF/ks6XBMNnnfPJBme8ku6\nx1PqUlm95exoRSbWJrw1zcwvZwsQs6srnUmLwNnkEVG+M9YgWL77ieGi0ImRwxtjiaGsaLAAbkmM\nQYdtlyDfmGgTGyX3PQysIxd2IYlWxW+xjIMTuY3X9+9bKxKJ8f8A+RVjmWybBhGUaWXjT4p0Osro\n5HVVQ0PAiYWlziba203Lrn4+VtwnUibWPCR9JBhk7psjZZG5WEt1cQdAqscfTq2knw6kdAyN4Ms5\nGTOOLbbyr2MRTVsdNTxxHJLIDK62rANf5LUuFJnxwkn7X7FKi6XXnVKKLpdNCVVqfvlH7TvpKsqt\nU/fKP2nfSVYLKlQigKVClBCIiCVCKUEIiIClFCCVCIglFCICIlxe19UAgEWIuFgSVTcW8IY6OMnY\nUX2jzcWe/cB6h+qtYlUS1k3i2idZzh9vLv2TeXrKuUtBTUjYhFE0OjZka7jbfv8AWu2Pwm739J2s\novKpqY6WNskpIaXtYLDiTYKaioipYHzzvDI2C5JXLVVMUMcIcImNYHOLjYbyeK7XMUjZomSMvleA\nRcWKxHVeIx4hVMgmZUmF2c074y05Duyu571vHG5b9o0cRpG1PR5HS7Lo8okzdliLfNe9NUw1cImp\n5BJGdzgs6vqIMR8Hppmh5YW3LWi7gQd1lmQjFq/D6YspujxjNmjicIi48HerXcF0njuWPu9G1zwj\npJpIhJ0yRsTpYm7INFtXAXWtSU76eIskqJJyTfM+1/Vos6OgqK/BWUeIh8MsZb12vDsxbx/5Wvo0\nAErOfqcfwR0i8zNEN8jB/wCQXa5WWdqLJo8JlpcREu2DqePOY2W1aXWuPVp81rL5qsqcaifIXStg\nG3axn2YLS0nQ3v3rt4ZctyXSV9Jkbnz5Rmta9tbKVktra+hmYzEo4nwPcGNnhuLE+cCtZcssbFFK\nhSsgoREEooRAUOa1ws5oPrClFQAAFhoEQkNFyQBzKAhwuCCOxPYpOwiifXOrJIs8xIN3EkAjdpuV\n3cFl4/JI2CCNsz4YpZmxyyM3tB7eGtl5V2Kijhlp42OqzBCTLI5waBppfmTyC68M/Jr7Tcjzhxam\nbj1ZtqljWNZHHHrcHef5rbjkZK0Oje17TxabhfJ+Ds1LW4j1aBkcUlOA1rrHVpNz8/kvp6WkgpA8\nU8bY2vdmIbuut/6jxzDLje0l2sKni0ElThlRDELyOZZo59itoSALleeXV20yMGY+WqqaySnfBmay\nNrHi1gBr81r2F721XAmicbCRhPtBd3W/Jbct60kEUqFzVKKEQSiKECw5BVagAVlJ7TvpKtqrU/fK\nP2nfSVYLOUcglhyRSoIyjkEyjkOSlEEWHJMo5BSiCLDklhyUogiw5JYclKIIsOSWHJSiDnKOQU2H\nJUsUqZaeOFlPl208ojaXagbyT3Aq6tWWTYWHJLDkilZHlM+KCN00uVrWAkuPBfM0dVU4rVy1FHHr\nKMjZHg5IY7/Nxtey+nmZHJC9krQ6MjrNPEL52irm13hFTwUTXR0lLAS5mXLqdBcL1eCbxyuuma26\nGghooBGwZjfM57tS53MqzlbyClSvNcrbutK9ZSR1lJJTyCzXttcbxyI9S+eo5X4hicVBXyMtSEm3\nCoeP1tfUc1u4s+pZhs5o4zJUFtmAb7nS/uVKPCYqrBIaeSOSCVgJa4mz2O11uOe9d/HlMcbv/wAS\ntcNaNzQLLCqqumw/wgqHzAky08eVgFy5wLrAD3LUZUNp5aejleXTPjJDtwNrX/VeYpX+PHVRYDHs\nAwE7w4OJ/ms4ax3y/BXGC08kdGXzxhj5ZHS5PMvwV/Zs8wbrbuHJUhiTW427D36F0QkjPPmFfWc9\n73fsiMjfNCzcUjwuaSKHEHtY6znszPyjhfXctNcSQxTW2sbH23Zmg2WZdKymYFg9SYZ4Io3MjNwY\n3XDj2nitiw5BcxxsjbljY1jeTRYLpLlaIyjkFh+EjqljaeOCjZNEZGftWub+St5ZOMueavDIWjqv\nqMzvcLrp4brOJeni+mq6toqMTMdPTR2mdBH1iS3XrHjuWvBJFUQMli1je27Ta2hXjirsuFVZIJ+x\nfoBc7iuME/BqO3oW/omV5Ycv2LthyCWHJSi4qiw5JlHIKUQRYckyi1rKUQUqfEKaenmmByCFxbJn\n0ykc11hk76uhjqJIhGZLua3/AEnd8ljeElNIKmI0TXulqurNE0i0jG6n38Fqx4nSMo4JZCadkjtm\n1sgsQ7db5Lvlh8ZcftNuPCCLaYFWtYBfZHuGqoMhGCinqYHO6HNlbNESSGk7nNWjitTB4smG1adq\nxzGWN8xynRUJxU1dHSU7zBBSPbGXSPf1nWF7Ae5b8f8At99bStDGMP8AGOGzUzC1j3DquIvYqtR4\nBSNihfWwRzVTWBr3m5DiOPatKaphg2QkkDdq4MZ/qJ3Beq488pjxXSvDh9HA5joqaJjmCzSG7gvf\nK0aWCOcGNLjuAuV8zWVNNWNlNBV4hNM+7mthvZp5aq4YXyX3S+n0+UcgvKppoqqB8MzczHixCiia\n9lHC2TNnDBfMbm/aV7Ln1VYlLhmFTmak8XtY6mIYXOtmcODrjmkYOG4xT0UUr5aeoY47J7sxitax\nB5diq0cVVPj+MU8lQ+IkNLHs3ht9PlotenwqClZIafM2oe0tMzjmd816c/hbLe5/eMxesOSWHJZ2\nD1EssdTHPKJTTymPagWzCwPfqrtPPFUwtlgeHxu3OC8+WNxulelhyCWHIKUWVRYcksOSlEBVan75\nR+076SrSq1P3yj9p30lWCypUKVAUIpQQpRQglQilBClFCAilEGPXPDsdp7kkUsEkxF+OgH81Tqca\nqKbBKSVzvt6nUvy5sg37hxtZaxw2N2KurnuLnGHZZLaAX3rmpp6TD8KqdnA1kQY5xa0di9OOWHxx\n1vpljYNj8ceHSPdFVzBj3Oe8i9m8yd3uC2m4rA/FG0EfWfszITfRu7T5qjhVEyt8EaamDzE2WIXc\n3nvKsxYDRsbCXB75Y37Qyl3Wc7t7N2nYteS+Ld9a7Jt5sxprMPmqahnWjndC1rP2iDoq3grCZjWY\no9mV9XKS0Zr2aCrz8CoZC/asfI10m0DHOJa1x329at0NHFQUjKaG+zZe1+9ZueEwsx7v9jV291Kh\nF5mhFKIKE2FQT4rDiEhcZIWZWN4Dfr81eRFq5W9jyNNE6rbVFv2zWFgdfgTf+S9UUqbEKVClQQil\nQglQpRBBAIIIuDwXMcbImBkbQ1jRYNAsAukQEREBSiICKFRw7EDVUD6qZojDHPBtyB3rUxtmx0aJ\nzsYFa6Tqth2bGcje5P6KcTbeje5tK2qkb5EbgLE+9WY5GyxtkYbtcAQeYTax7XZZ27S2bJfW3Oyv\nK7lv0j53CsFqnhvjFrIY2Zy2KMi+Z+834WG5bc9FDJRinfE2VrG2a14vwsF6moiFSKfP9qWl4b2D\nS/zXnDWxzVk9MARJBYm/EHiFvLyZ5ZcjUZeD4dXDorsRLWtpGZYomkHW1rk/ILcWbjtVUUVCKmmc\n0bN4zNcNHAm1uxe+HV7K+DOGmORpyyRu8ph7VM+Wc5/RPwtqtiNa3D6KSpcwvDLdUG17myzq7Eai\nn8IIKUTRsp3Rh7g4au3314HQL57F8Sp8Qpp5Gmt6zmlmdxEflaiw03fqu3h/01yym+qlr7prg5oI\nKlZNPhNHh2SoiqJYGNIc4GY5HacQdFfZW0kjwyOqgc47mtkBJXnykl9e2nbYImzunDAJXDKXcSF6\nIvm6aoxSCmmqnVAm6PK4TwPtcAcj6lrHC5ze03poN8H6IPlOaoyyOLyzakNud+gWjBDHTwtihYGM\naLBo3BTE8SRMeNzgCu1i5ZX1aqERFkEUog563MdyrVF+mUlz+276Sraq1P3yj9p30lUWLO5i9uXF\nOtzHcukQcjNxI7ks7mN3LiukQc9bmN/JAH8SOPDuXShBHX5jdyTrX3i1+S6RBzZ1tSL+pLO5juXS\nIObOvvFr8ks628buSlSg5s6+8dyy/CZ0jMAqy0jybHTgTZay5exsjS17Q5p3gjQrWGfHKZfhFDA4\n3R4JRt0H2LdLdiv9bgRv5KRopUyy5ZW/lXNnX3juTrcx3KUU2I6194tfkos/mOPBdog5s62hHcnW\nvvG/kukQc9fmO5LP5jhw710ibHPW5jfySzuY3cuK6UJsR1uY7kGbmO5dIg5s+28Xty4oc3AjfyXS\nIOetzHclnW3i9uXFdIg563Mb+XBOtzHcpUoOetzHcnW5jfy4LpEHNn23i9uSdfmOHDvXSIOHB2U6\njjwWVhA2uBxxMkjuY3tLTvzX1/VbCoU2FU1NWy1Mea8otkv1W332HbZbxykxsRSwl1fJhNKKaWCz\nGljjI0k6Gy6xCimNO+smqYoqmnaXNmjYRZoG466hasEEdPEIoWBjBuAXUkbJY3RyNDmOFi07it3z\nfL1DT5SgxSqiLq+ppn2naXvly3a1otly69puFqVodDidFiDC3JKRC8kfsnUe+60qukirKY08t9mS\nCQ3Tcb2XhitA6to2U8TmxgSMdqNwBut3y4ZZS61/0mlXwnbMcBqD1XZbOdYcAQf5KvBPU0tSypqY\n4oen5WNLTmawgaX3XJW7NEyeF8Ugux7S1w5gqniWHuqqKOCnkbCYntcwubmtlWfHnjrhVs+1aeip\nMQxl20BdNTBjja4FtbDfvvqtCWkjnp3QTMjdE8EOaG2vdc4fQsoYSwOdI97s0kjt7zzKtLGfkt1J\nfUNPCalZPBsZGsdHp1SOSyH4BK6dskc1NBkcSx0MFnN991vIpj5LiaVaOklpo3h1VJO92pfJz7By\nXztFh75cerqOtqJJfJmy7mPGm8DXivrFwIo2yulDAHuABdxIC1h5bjv9mmVjcclNG2tgqXxSsLWN\njaOq8X0aR71qHatiJc9gIbqSNAe/cvOvooq+DZSue0BwcHMdZwI4gqi3wepdRNPVztP7MsxI+VlO\nWNxko1G5jc5mkG1rBTZ194t6kYxsbGsY0Na0WAHBSuaos628buSdbmN/JdIghVan75R+076Sraq1\nP3yj9p30lBZRFKghFKIIRSBdTlWuNNuUXWVMqcam3KLrKmVONNuUXWVMqcabcousqZU4025RdZUy\npxptyi6yplTjTblF1lTKnGm3KLrKmVONNuUXWVMqcabcousqZU4025RdZUypxptyi6yplTjTblF1\nlTKnGm3KLrKmVONNuUXWVMqcabcousqZU4025RdZUypxptyi6yplTjTblF1lTKnGm3KLrKoIsmrF\nQilFkQilEEIiEhouSAO1UERrg4AtIIPEIgIvnKPwnY0v8YlrATeMsad1yDf1ELfpqiKqgZPA8Pje\nLtcOK3n48sO4m3oilFzVF1VqfvlH7TvpKtqrU/fKP2nfSVYLN0uiKBdLoiDoblKgblK7zpkREVBE\nRAREQERZkGMxzOgLqeaOKckMldbKT36bkGmi8OnUZjL+lQZAbF20FgUfWUzWOcJ4nW4B4uSRcD3h\nB7osunx+hqJI2NfYu3lxADOrm1VueugiopKtr2yxMFyYyCgsovKWqp4A0zzxxZt2d4bfvUdLpvzE\nPPyxzt+qD2REQERZ9RjVHT1JgeZCWmz3NYS1nrPBBoIudrHmc3O3M0XIvqAo2sWTPtGZQbXzaX5I\nO0Xj0unzStMrAYbZ7m2W67MsQcGmRmYi4GbUoO0VXxjS9CbV7UCBxADiO225e+1jzObtGZmgFwvq\nAdyDtFw2aJxAbKw3vazhw3oJYi1zhIwtb5RzaD1oO0XJljGW8jet5Ou/1LydWUzHRAzM+1cWssbg\nkb/0Qe6Ll0kbAS57WgC5JPBA9jnFrXtLgLkA6gIOkREBERAREQERQS0byB60Erl5sull4ziclAY2\nQU5nkeC/L/pFr+/UJxufqDRus2Crllx6qpyfsYomkNtxK8qjFqapcKOCrfDO9wZmbHmAPFt911lY\nDFJWY1XzxVk2SGRrMxdm2gF9D3e5b8fh+GWWXr1/ktfWXS65kLmxuLBmcASBzKy2Yweg1b6iEw1N\nKOvFmvw0se1eeYXLpWtfTj3LOxagirY2yTOfkia47P8AZdccfUvfC45osNp21Ej5JsgLy83NzrZe\nGPPqY8JnkpHMa5rSXZhfq2N7dq1hvHPWJ9MXBGGjfhTKaRzpJ4c08RddobvDrcCvp5xI6nkELskh\nacriNx4FZ+A4dBRUUbmXfNIxpe9xud271di1Frz5y57hOmHSeDsPRyMRd0idzy5z23aNd404LRwy\ngZhtMaeOR748xLQ79kch2K4izn5c8+6aRdLoi5KlVan75R+076SrNhyVWpH+MpPad9JVFpFGUcgp\nsOSAiWHJMo5BB0NylQNyldp0yIix341JYtipHOlMrmMY45SQG5ie7gqNhFk0ON9NrGRCAtikDtm8\nu1JaATp71blxSggldHLVxMe3e0u1CC2io+OcN/OwfEnjnDfzsHxILy+bjwCaCnp3QtYZwxzJWveS\n3UHrDkVreOcN/OwfEnjnDfzsPxIMalwWuhmM7oYZCC0iN8gIdoRfQdvJdjBq5j2tDY3MLo5C7Nax\nawttb3rW8c4b+dg+JPHOG/nYfiQZdPgdQzD6unfHDmlZHkN+IaARu01HzXvDhlY3AKqkfkMspdkF\nxpfdcgfyV3xzhv52H4k8c4b+dh+JBlYrh+KYk1odBDGGxlgAlvqba7uxebfB2p6cx5cwU4cAWg65\nPKI+JbPjnDfzsHxJ45w387B8SC8io+OcN/Ow/EnjnDfzsHxILywjHiFHUVNPT0cU8dVK6Rsrj1W3\n35hZX/HOG/nYPiTxzhv52H4kGHLhdaRK0wvMw2rnzB2krXAgAAfpwsorMJqC2VrIJMge1zQw/wCg\nA6cdVu+OcN/Ow/EnjnDfzsHxIMKegrGAhuHtdIQwuIuWnq2t26815Ybg1WKaR9XA8ziaMMJd+wDY\n+6y+i8c4b+dh+JPHOG/nYfiQZvQJD4NdFNK8SRSDMy3lWcCSOeizzg9dnqssclyx99dHX8kDmvov\nHOG/nYfiTxzhv52D4kGI7CZoJpnxwlrGztLSDuZkN/duVSDDaqpoYpaeGSKARMEjMxzSkX6wHvC+\nmOMYY4EGtgIO/rKG4vhbGhrayBrQLAB25Bjw4RMaR+0gL5GwO2IebWJcSBv32sq1PhdWZTNPRHLt\nbtYBuBYRcAnTWy+i8c4b+dh+JPHOG/nYfiQZNNRVFSWMqqR7I2thjIJOuXNcphOF1dNjG0eHhjC/\nM4nRwPk25rW8c4b+dh+JPHOG/nYPiQXkWfJjmGRxl5rYiBwBuSvKl8I8LqmktqRHbhKMpPerq9jV\nRUfHOG/nYPjQ4zhg310A/wDNQXl41FVDTGITPDTK8MZ2uPBUHeEmFCoZCKpri4XzN1aPWeC8sXrK\nCrwyWSKriL6ctlY5vWs4Hq6dp0W5jd6o1hUQ9J6NtBtg3Pk42va6w/CeihDWVzw979pGzIXkNtfX\nvuvTwcqGTOqHTm2IvdmnYRbLyA7B/Nc+FkxNNBRws2lRNKHRs4Oy6ldMMbj5JE+no+DxTX0baN8h\njqJNm+Bz8wtbyhfdZMZp6ybFsNdSsaGROLpJTwG63vXvRQCKsEldOJa+VhLRawY0bw0e9aLgCdQF\njPPXvskYTvBtjt1dO3I8viAtaMnf61fjhocHpXPaGwxADO48baaq7Yclw+Jj22cxpHaF575MsvWV\na0xMOxN9TjdSynZLNTOy3kdo2PQ7r81V8K6R1XWU8FF97kYS9odbMxuov71s4ZhzqJ1Q6aczumcD\nctA0AsFc2EO22uyZtBpny6967fy44eTlh9JrcfLST1UeIdHgfUPqG0gjLL3aXloF78AN5K+kNMZM\nP6LPI5xdHke/idLEr3DWgWAFvUpyjkFzz8vLXrSyPnqulrIcYjOGU4DTDkfM9/VvuBI7AFt0cJp6\nOGEvMhYwNLz+0ea9creQ5oGgCwAAWcvJcpJTSUSw5JYclhREsOSWHJBKq1P3yj9p30lWVWqfvlH7\nTvpKQWVKhFBKhEQdDcpUDcpXedMiz6rCYqkvcZZGPMheHMNiLtDSPVYLQRUZ1Jg1PR1QnidJZoIY\nwuu1t7XI9dlFSxpx6iu0axS309laSzqn8eov3Uv9KC9so/Rt7lDmwsAzNYLmwuOK9FXrQ4xxloJy\nyNcbDgCpldTaybeuzjvbI2/qU7KPzG9yryufIQ5gcGaZuBtdVyZxTy5DI5xa8Mvvvc2XP+WNcF4t\nhDS4tZYbzYaKRHGRcMbb1LOqmVToqhoL7OjeGZd+bgtKMWjaDyC1jnyuks0jZR+jb3Jso/Rt7l2i\n2y42Ufo29ybKP0be5dog42Ufo29ybKP0be5dog42Ufo29ybKP0be5dog42Ufo29ybKP0be5dog42\nUfo29ybKP0be5dog42Ufo29ybKP0be5dog42Ufo29ybKP0be5dog42UfmN7lDWwuvlbGbGxsBoVF\nTTx1UDoZgSx2+ziPmF8tFDQU9LO8SzR1PSJGwiKR2ZxBsBbj710wwmUR9Xso/Rt7lwDTuldENkZG\ni7mC1wPUvNj6qLC2vdGJapsYuwG2Z1lh01PiNRjlcXVMdNNs48+zZm01tYlJhuW76Gh4Qtjhwt9U\nGMz0zhKy7dCR/wAqmK+mlxaOsEZNJHEYxIyIkOedTuG4W+a9vC2KolwqNkGZwMzBIAN7f+bK5idP\nVnDuj4UYoXmwudA1vG3at4zHjN/exZpn01VC2aAMex245Vm42yg2sAnoJKuUasijjJBvz4d6t4PS\nVFFS7Gd0RA8kRg6c7k71eXO6xyuhnYXSGPavloaamz2s2PU+/RU6yiNT4T02VloIotpINwcbnL67\nFbqpMH/W5D//ADt+opz97HhjVA+pNNLTRgVEczftAcpazj6/Ur76eKSoinc0GSK+Q8r716opyutK\nrT0bJ62mqnEh1PmyjnmA/srDt6lQ7esZ22EcqVCLg0IiICIiAiIgIiICIiCVVqfvlH7TvpKsWPnf\nJVqi/TKS5v13fSVYLSlcWdbytbclNjz+Sg6Rc2d53yUWdbyuHLig9BuUqBuUrvOmREVTEsQjw6Fs\nkjHPzOygAgcL8fUqLazqn8eov3Uv9KUmMQVdQ2JkcgD75HuFg4i1xz0uomlifjtK1j2OeyOQOANy\n3yd6DSREQFAaGjQW1upRAREQEREBERAREQEREBERAREQFhYlSf8AUIdvidUwzOdswxzWtZbXXTVb\nq+dxulgxLwhw+iqWSGMRvfo6wO7+y6eKbqVcoKupArY8xrhTkCN4s0vNtQTu0XmMWraeutiMMFNS\nlhIdnu4Hlv17lpdChbR9FhBgjtYbI5SPUvClweipZDK2LaTE3Mspzu7yryw7sFiiqTWUrZxG6MPu\nWh28jgffvWTgeFupaOoqZWWrZnvfmLQS3U2st1FmZ2Syfaq2Gy1E1BDJVs2c5HXba1iuaTD201XU\nVJlfLJPa5dbqgbgLK2inLv8AYIi8pamCFzWyysY53ktLgCfUOKxboeqKmK18sjmQU8hA/wC5IMrf\nnqVl1Ne5zJGuqhMd52J2ccY03vU5b6i6/LZqq2mpIy+eVrQPeSqGH1jK3FpZY2PazYNALxbNqdyz\nMLwR1XKyebqUzSCGhttof7dp1K22aY1IBoBTtsP/ACK1qztF5ERAUOUriQOJ0dbTTRZy6EoubOv5\nXHkgDuLufBcWnSLmzvO+SWdffx5IOkXIDvO4cks7zvkg6Rc2Nxr8ks63la25IOkXNnXPW+SWdcdb\n5IOkXFnedw5KbOv5XHkg6VWp++UftO+kqyqtT98o/ad9JVgtKVCKCUUKUEjcpUDcpXedMioYph5r\n2xZZGtdGSQHNzNNxbcr6KjJoMF6HUxybfPHFmyNLdbuAvc+5ek7GNx6kysa0uilLiGi58nitJZ1T\n+PUX7qX+lBooiICIiAiIgIiICIiAiIgIiICIiAhIaCSbAbyUWRUl2L1bqONxbRwm072/9x3mDs59\nysmwOI1tZO3xbTtNKHgPnkPlC+uUdnNaxa0vDi0ZhoDbUIxjY2BjGhrW6ADgpVt30CIuJJY4ml0j\n2tA1NysW67HaKmyvNRG51JA+S24vGRp954KG09bPBlqqhsTideji2nK5upy/C6WKiohpoy+eVsbR\nxcVXfXSSNb0OmfNm3OccrR36/JesNFBCBlaXEbnPJce8qKuvgo8okcS92jWNF3H3J8qenElNVVBb\ntKkxNG9kI1J9Z4KtNLQ0k5EEHSaw7mtGZ3vPBcTOq6lhkrJBQUo1yh3XcO08FzTtnlYIsNgFHSO8\nqZ4+0d2gfzKuPj+0uVeVXJLK/Y1bjJM7VlFA63xO5K1S4TctfWZC1tiynjFo2Hn2lXKKhgoYskLS\nXEdaR2rn9pKsre9dJoVFn43J/Dt+oq8qLPxuT+Hb9RWVXkREBQ5SodvWcuiIRQi4tJRQiCUUIglE\nRARQiCUUIglVan75R+076SrKq1P3yk9p30lWQWlKi/Yl+xNCUUXS/Ymh0NylQNyzsSxboM8cTIHT\nOdbMGndc2C7TplpIvOnmZUQRzR3yvaHC68cQxCDDohJOT1jZrW7yqLSzqn8eov3Uv9K9KbFaSqqN\njC8uJF2ut1XcwDzC5qB/1uiP/wCuQfSgvoiICIiAiIgIiICIiAiLPxHETRysaNl1hez3WJ7BpqVn\nLKYzdWTbQReElXHFEySS7Q7fp5PrXsx2djXgEZhexGqsyl9Q1Uoi4mmjgjMkrwxrRcklVFHF6qRg\nio6Vw6VUnK0+Y3i73BW6OlioqVkEIsxg7zxKpPnjmnZVUdM+omDSxr75WgGx1v8A2Xu6CrqA3aT7\nBu9zYdSezMf7JfJLJMV1+VmWeKAAyysYDoMzgLqs6vLphFT000vN+XK0e87/AHL0bQ0zZNoYg9/n\nP6x+asblj5Xv0elNsVbLJmmnZHH6ONtz8RXcGH00MhkbHmlO+R5zO7yrK5kkZEwvkeGNG8k2VmMN\nul5VNVBSR555AwcBxPqHFZdXjfUJpw2OEGxqJtG+4byq9MyeecSwUz5pCPvdUbNG7yW/8LcxtZ2t\nz1lTPGXsc2hpRvnm0J9QVejdZzhhlI6Rx8qrqSQHdo4lXYcJZnEtbK+rlBuM56rfU3d3rPx7D+l4\npSsDi3aRuaw5iAHNsRoOBFwumExt0NGDCmbRs9a/pdQNQ546rfZbwWgs6lwkUsjHx1VS0DV0ZkzN\nPZrqtFYy79VRERZBUWfjcn8O36iryos/G5P4dv1FBeREQFDlK5ebcCs5dEQpUJfsXHTQiX7Cl1dC\nUUX7Ev2JoSii6X7E0JRRfsS/YpoES/YiuhKq1P3yj9p30lWbjmqtSf8AGUntO+kpBaUqLjmEuOag\nlFFxzS45oOhuXznhMY21sJmpjOwMJtmLcpvzC+jG5fP+Es5ZNFGyeSNzmgWEpY03NtefH1LvOmW3\nRva+khc1rWtLAQ1u4KpjOGvxFkWzcwGMnR97G4tw5KxhpYcOpzE0sZsxZp3jRWVRjYfg01JVxPfM\nx8cWct0OYlwF79y9JYI4/COmlY2z5YpC887ZbLVWdU/j1F+6l/pQaKIiAiIgIiICIiAiKHvZGwve\n4NaNSSbAIJXmYGmo2x1IbltYW371XlxGIZRBHJUudu2Tbj4tw71Eja+ctySMpmbz1c7vVyWOUq6s\nd1dHBKOu7Zgm7sptm9aqU0zRWSGndLVfZNaXfslwvx3A+pWhhtOZxPKHTSDyTI4kN9Q3BWwABYaB\nZmH301y/qpNjrp2O20rKe+4Q6uHvOnyXpDh9PCB1XSEftSOLj81ZRb4z7Z3RFDnBrSXEADUkqjJj\nNAzRk4mPKEbQ/Jbkt6RfS4A1WBPj9U4hlJQO2p3NlPW323Dd77L1gwuurY74xVuc0/8AZh6rffbe\nrx/NTbrE8eZTvNPQxmqqt2VguG+uyqMwjEcQftMRmEbT+xo4j1cB3XW9T0sFJGI6eJkTBwaLL1Tl\nJ0aUaXCKSmDTkMr2/tynM75q8iKW29qKnV0kk9bRzska1tO5znAjV1xZXESXQIiKAiIgKiz8bk/h\n2/UVeVFn43J/Dt+ooLyIiAocpXLiAdSs5dAii4S45hcWkoouOYS45oJRRcc0uOYQSii45pcc0Eoo\nuOaXHMIJRRcc0uECyq1P3yj9p30lWlWqfvlH7TvpKsFlLIigIiIOhuXzfhTK6Opha42ie23VAvv1\n3j1fNfSDcvm/CAibE2QsL87WBxyxGS2vYu86Zb9EJBRwiYgyZBmLdy9l4ULGx0ULGZ8rWADOLH3h\nessjIY3SSuDGNFy47gqOlnVP49Rfupf6VbZV08k5hZMx0oFywHVVKn8eov3Uv9KDRREQERCQ0XJA\nQEVR+IQNmbFHnlkPCNpNvWeCgGvll3RQQ8jdzz/ILHOfS6W3Oa0EuIAHEqoMTp35xTl1Q9lwWxjj\n6zopiw+Nri6aSSd53mR1x3bvkrTWtY0NaAGjcBwT5X9HpTacQqIjmaykJ5HO636A96luGwuYG1Dn\n1JGt5TdWZJ4om5pJGNHMlUZcdoGOa1kpnceELS8/Ja/j37NtBrWsFmtAHIBSsw4jVzNJo8OlN9zp\niGDuOqiNmNy2201JAOOzY5zvmbLfDXaNRVqnEaOlF56mNmtt9yqviYSvz1dbVTnzc+Rvc2ytQYdR\n01jDTRMI4hov3p8RVOMCUWoaSep5ODcrO8oBi9SzrOgoweA+0cB69Ar1TUw0sZknkaxvad6+fqfC\nPpEwipHOZGdMzWZ5CeQHD3qz9RNrVVS0VM1vjGqnq5D5MbneUeQaLfNdwUVRUNADBh9MN0cVs7h2\nngvGiNRCNpTYTLI92+aeZocVaOI18ILqjC3hg3mKQPI9y1cc76t/5RdpaSCkZkgjDRxO8n1le68K\nSrhrYRLA/M3cRuIPIjgV7rl00IiICIiAiIgIiICIiAqLPxuT+Hb9RV5UWfjcn8O36igvIiIC5cul\nDt6zl0OUspULi0WRFKCEspUICIpQRZLKVCAilQgKrU/fKP2nfSVbVWp++UftO+kqwWURSoIRSiDo\nbl890avp64yR0YkLal8u0DwC5rrgN/TuX0A3KV3nTLmNznRtc9mRxGrb3sfWs/HMPlxKj2UU2ztc\n5bXDzbQLSRUYWHYZVwYkx8wYWRukftAfLL7cOyysyQbPwgppM8rtpHIbOfdrfJ3DgtRZGMurOnUj\ncP2XSCyQja+TbS/8lYNZzg1pc4gAakk7lVkxCIMBgDqkncIhmB9+4L5p9Pish/6sKwsBvamDXD9V\ncg8UPLWTV1S598obI9zCD6hayknvVq69NSplrHRZhLDRMtcvks4t/kq0suEve01NYyqeNzc+YX9k\naKxDg+FnrshZJrvLy8fMq5HTQRf5UMbPZaAkxxl97pus441EC1tJQ1cwOnUgLQO+y9Ol4nL/AJWH\nsj/fS/2utJFqWTqIyzDjMr7uq6enbbcyPP8AqujhTpQek19VJ2Mdsx3BaSJMrOh8zhkuDSVUUDMM\nMbn5mxySNDsxbv1X0DZaZjsrXxtIdksLeVa9vWs3DcBipacl4Dap2YGaMm4BN9Lrwgw2UY6S3aGC\nFmbPLrnltYO7bCy6WY3clRvoq1BthRM6S4mQXzE8dTqvGbFGGTY0cZqpuIbo1vrK49KuySMiYXyP\naxg3ucbALONfUVt24bF1L2M8mjR6hxUxYdJUPEuJS7Zw1ETdGNP81pABoAAAA3AIj53FsMbFBFNU\nyvnLpmNnkcbWYSbgDgLkK3jNGxmGRwUsMbftWBrQ3Tf+i06iCOpgfDM0OjeLOBWcKDE4WbKnxFmz\nHkmWLM4Dle4W97x0aeVEyowpmwMW1kmL5A1hsxpAvlC8enVpxSFsglijdYyRgXDXWNm5u3QqyaXG\nRKHmqpJSNWl0RGXT1r0EGMuHWq6Vh7ISf5q8PXcRVdmMRxejiMUzSRNCDfaNBsf/ACWzBMyogjni\nN2SNDmnsKyosEqLPbNiUxieSTHEA0a7+ZWtDEyCFkMTcsbAGtHIKWzX7V2iIsKIiICIiAiIgIiIC\nos/G5P4dv1FXlRZ+Nyfw7fqKC8iIgKHb1Khyzl0OUUouLSEUoghFKIIRSiCEUoghFKIObHziq1QP\n8ZSa367vpKtqrU/fKP2nfSVYPfKbeWd1uCnKfOKlSmxzY+cVGU28s7rcO9domxLdylQNyldp0yIi\n8p6iGmYHzysjaTYFxtqqPVZ1T+PUX7qX+lW2VVPJM6Fk0bpGi5aHa2VSp/HqL91L/Sg0V5zU8NQ2\n00TJBa3WaCvREGa/AqK5MIkpi7fsXlo7tygUOJU4HR8R2oGmSoZcW9YsVpoi7rLNVi9Pm21BHUDe\nHQPy/Jxuhx6mjaDUw1NPzzxGw9/FaihwDhZwBHIoK8GI0VS4iCqikI4Bysgg7iqs+GUNQby0kTjz\ny6rMqaGhw5toq6ppXHVscchd/wCvFD03VTq8Shp3iJt5p3eTFHqffyWSIvCGoh6s8bY3bg8bN9vd\ney9qWSbC4i1+EPuDq+Bwfm7ddUTSz0SrxAXr37GK/wDkRneP9R/stCGGKnjEcMbWMHBosqMeOUD3\nBr5TC4i9pWlnzOiuw1EM4vDNHIBxY4FDWnoiIgIiICIiAiIgIiICIiAiIgIiICIiAqLPxuT+Hb9R\nV5UWfjcn8O36igvIiIC5kBO5xGnBdKHb1nLojixv5R3oGm3lnjyXSLltpzlPnFMpv5R3rpE2OQ0+\ncUynzjwXSJsc5TcdY70ym3lHcukTY5sbnrFMpuOsV0ibHGU28s7uxTY38o710ibEKrU/fKP2nfSV\nbVWp++UftO+koLKIpUEIpRBI3KVA3KV3nTIsvHYpZYI2w0u2cSRmtfIDvsOa1EVHz+F4bPBiETzB\nsoozI65cCSHAZR6xZWpIns8I6d7pnPa+OQtYRozyVrLOqfx6i/dS/wBKDRREQERQ97Y2lz3BrRvJ\nKCV5VFTDSx7SeRrG8yqJxGerkMeGw5m8aiQWYPVzXrTYYyOUT1L3VNR579zfUOCI8ekVuIfdWdGg\nP/dlHWI7ArNJh0FK4vAMkx3yyG7j/ZW0Q0IiIrmSOOUWkY145OF1RlwTD5NRBs3cDE4s/RaCIbZb\n8KqGi1PilSxo4Os5L4zTm7hTVbbcCYz/ADWoiLtkjGnxNHTMPqoTeznBuZg7bq1T4tQVLssVUwkC\n+un6q4vCooqWqINRTxSkeewFD09mua9oc0gg7iFkeE9TUU1DA6me9r3TtadnvIIOi9X4FS2+wfPT\n63GykIHduXlPhVa4xlmIbVsZzNZURhwuNxuNUlWa2jDK6SlpAzFJHCdxL2McLvycL246FeYx10tf\nTR0+z2FTYsc/qloB61+08FzUx4z06GeaFs0cOoZA4AHtIdrddS12GGXaVtHLDM5wJMkR0Ld2qu4v\n7bFPUw1Ue0geHsuRcc16rPwuShbG+OjqY5QXueQHDS5v3arQ3qMUREQFnYlNM6dlJBNsM0bpHy5b\nkAW0HbqtFVq6hirmNEjnsc3yXxuyub71ZqX2M6LFnQ0IcXNqixxDpP8AL0HO/wC1qrEGMMqMRFMy\nF2UtBzEgHUX3ctd/NQcCpDGxuaa7b3fnu597XufcN3JerMKgbURTZ5XbIDKwu6oIFr252WrxReRE\nWFFRZ+Nyfw7fqKvKiz8bk/h2/UUF5ERAUO3qVDt6zl0RyilFxaQilEEIpRBCKUQQilEEIpRBHuVW\np++UftO+kq2qtT98o/ad9JVgs37Ev2IpUEe5L9ilEEjcpUDcpXedMiIioLOqfx6i/dS/0rQDmlxa\nHDMN4vqs+p/HqL91L/Sg0UVSrxKClcIheWc+TEzVxVcUtZX2dXP2MJ12ETtT7Tv7Im3dRijA/Y0b\nDVz8oz1W+s8FyzDpah4lxKUS21bC3Rjf7q/DBFTxhkMbWNHABdoIaA1oDQABwClERRERAREQEREB\nERAREQEREBCA4WIBHIoiCpUYZQ1JvNSxOO6+WxVYYQ+Cxoq+ohtua87RgHqK1ERd1l3xqnuLU1W2\n973MbrcrWsjsXlhy9Jw2pjHFzRnaO5aiIbUIcaw+YkCpYwjeJOofmrrJGSsD43te07nNNwVxNS08\n/wDmwxv4dZoKovwGivmg2tO+9w6OQ6eoHRD000WX0HEYXk0+JF7SLZZ2B1u6ygVGMwR/a0cFQ4b3\nRSZb+4hDTVRZfj2JjgKmlqqfS5dJH1R7wrMGK0NQSIqqMkb7m36oaq2qLPxuT+Hb9RV5UW/jkn8O\n36iiLyIiAuXnsXShyzl0OUv2KUXFpF+xFKIIv2JfsUogi/Yl+xSiCPcl+xSiCL9iX7FKIIzDmqtQ\nQaykt5zvpKtKtU/fKP2nfSVRZzDmlwiIFxzTMOaIg6G5SoG5Su06ZERFRg0dNM3Hi8wvaRJM58hb\no5ptk148dOC9MUpMRqa2PYSxRts4NcAbgEC91tEgAkmwCrQYhS1EUksUzSyPyydLdvqQUaWhxCkb\naEUIcfKeQ8ud6yve2L+dQ9z1apaqGsi2sD87b2OliDyIXsgz7Yv51F3PS2L+dRdz1oIgz7Yv51F3\nPUPOLtY516I2F7Br9Voqu+upo6ttM6S0rtwsbdgvuQVmHF3sa69ELi9i19wpti/nUXc9WG19M+rN\nK2UbYcLceIvzVhBn2xfzqLuelsX86i7nrQRBn2xfzqLuelsX86i7nrQXlU1MVJFtJnZW3A0FyTyA\nQUWuxd0j2f4MZbalr7G67ti/nUXc9WXVtM2j6W6Vogtmz8F7AhwBG46oKFsX86i7npbF/Oou560E\nQZ9sX86i7npbF/Oou560EQZkjsYZk+5uzODdGv07V3bF/Ooe56sU9fTVG0MUoIj8onQW5+rTeppa\nyCsY50Di4NNjcEfqgrWxfzqLuelsX86i7nrQRBn2xfzqLuelsX86i7nrQRBnkYvbyqLuevOCTGJo\nWyEUjMw8lzXgj5q5NW00FRHBLKGySeSP97kdXUzKttK6QCV3C38+aCvbF/Oou56WxfzqLuetBEGf\nbF/Oou56rz4fWVJvNBhjzzMbrrYRB854hr2yB0FVHTtG5sb3lo7bFXsNpa2DEZDWVBqLxC0mQNG/\ndp/vVaj3tYwve4NaBcknQLwfX0zKNtU6UbF1srgL3vyCLtYRcxSMmibJG4OY8XBHELpEFy4gHU8F\n0odvWcuiObjmmYc0Rcmi45pcc0RAuOaZhzREC45pcc0RAzDmlxzREC45pcc0RAVap++UftO+kq0q\ntT98o/ad9JSCyiKVBCIiDoblKgblK7zpkREVHjVwmoo5oA7KZY3MvyuLLI8XVz6KohdFCx0zWsuH\n3tlGh3a3K3UQU8NhniZM+paxkksheWsNwNAP5K4iICIiAseow6pdipnjDDE6VkxJdYgtaW5bdvNb\nCIMmCnrXYuJ6imhZC2+TJJq0kauOmpO5ayIgIiICpYnTzTxwOpw10kMzZQ1xsHWvpf3q6iD5+owa\nslwaGiEjWmFvO4ef/i3YWuZDG15Bc1oBI5rtEBF5sqIX5MkrHZ75bHfbeu3ODGlziA1ouSeCCUXm\nJoiYwJGkyC7NfKHYvRBhUeFVcIe1+S2w6O0h+8ZnHNu037lZwOgnoWSCYNaHBoDWuzagWLr9vLgt\nREBERAREQZOK4fUVVQTCI8krGMc5zrFmV+bTndeMmGVc+Lx1jmRs6zHO698uW9xa2t+fBbiICIiA\niIg5lY2SJzHNDg4WsVk9CrG4dSU7YYi6m2b9X6OIOo3aetbCIK+HwOpqKKJ/ltGut9d6sIiAocpU\nO3rOXRHKIi4tClQiAiIgIiICIiAiIgKtU/fKP2nfSVYsfOKrVA/xlJrfru+kqwWlK4ym3lndbgps\nfOKg6UKLHzimU28t263BB2NylQNyld50yIiKgiIgIiICIiAiIgIiICIiAiIgISA0k7giIPmaN76S\numq+hSGKdr3U7fMsddOF9634KhlQCz9sNBe23ML3cA5paRcHQhcsYxnktA3bkGFTQuw819Q1j5ej\nEsgY46NZ5Rt2a/JbzTdoPMI9jZGOY8Xa4WIPEKdAgIiICIiAiIgIiICIiAiIgIiICIiAocpXD2kn\nRxGnBZy6WClc5TfyjvugaQPLJ38lxVKKMp888OSZTfyjvQdIucp84plPnFBKKMpv5R3plNrZju3o\nJRRlNz1imU3HWKCVK5ym3lHcmU38o70HSq1P3yj9p30lWVVqfvlH7TvpKsFpSoRQSihEEooRBKKE\nQSihEEooRBKKEQSihEEooRBKKEQSihEEooRBKKEQSihEEooRBKKEQSihEEooRBKKEQSihEEooRBK\nKFRqsWpaSo2M7nNda98pI+SuhfRZYx+gLb7RwPIsKjx/RtcWyF7Dcfs30te+ncrxqNVFmePaEh2S\nRznAE2yEXsLqG4/QWGeVzHHeHMOh7U40aiLLPhBh4IBlNzyaTZS/HKJrY3iTM199QN1rf3TjRpos\n2LGqSaXIwkjM1gdbQk3/ALLyq8djilfDBHtXht2nMA0pxo11C+arJsUlrG0rKrN1tdmMh7d2thuv\nzvyVWrOJ7duyNTHHISIxE5+7nY+virwNvqaipjp2guu5xNg1gu4+5cRYjSS6NqGB3muNiOyxWbDD\n4sfTymAuldA5r8tyXSXabX71UpqWelrJ8zhTjILNbE6UWHHlwV4wfTXFr3FlGZvMd6+Z2tLe01C5\nzDfWDOLntBstSgo8OnjZUx0TY3A6Zm2IIWbjpV/pEHpo/jCrVE8Jq6Q7aOwc79seaURXSLPSIPTx\n/GE6RB6aP4wiJxDpEHpo/jCdIg9NH8YRE4h0iD00fxhOkQemj+MIicQ6RB6aP4wnSIPTR/GEROId\nIg9NH8YTpEHpo/jCInEOkQemj+MJ0iD00fxhETiHSIPTR/GE6RB6aP4wiJxDpEHpo/jCdIg9NH8Y\nRE4h0iD00fxhOkQemj+MIicQ6RB6aP4wnSIPTR/GEROIdIg9NH8YTpEHpo/jCInEOkQemj+MJ0iD\n00fxhETiHSIPTR/GE6RB6aP4wiJxDpEHpo/jCdIg9NH8YRE4h0iD00fxhOkQemj+MIicQ6RB6aP4\nwnSIPTR/GEROIdIg9NH8YTpEHpo/jCInEOkQemj+MJ0iD00fxhETiHSIPTR/GE6RB6aP4wiJxDpE\nHpo/jCdIg9NH8YRE4j5yqoJDUSyQmmexzyQHzm9r9hH+wuPF1QbG9HrvG3dYf+yItbHTcOk6xeaQ\n9U2Andv71zLhTs0ZY+meQ27s0/VJ5c/miJscHC5g4HLRG28CY2PzXRwyTKyzaMuscxdMd9z28rIi\nbEeKJA4lslLbKAA6XjbU6HndR4tmzWLKN1hcnbHX5oiuxJwqQZHDojjlOZplNr3PbysojwybZguF\nHck6GY6dzkRNjuHDXNlZtW0jozIM2WY3DezVbda2GSlc2lkp2TADZuLh1feiKUZkXjOGSWVj8PEk\nrruOYf74fNescmKmMsfV0rdQMwcL2O/uUoghkmLMcR0yke0EgF7hc679Fznxd7XXrKRhztcLOG7i\nPV/dSiCNrjG1z9Jo9B5OfQ/7v+i2mVEWRueeIutqcw3oil9j/9k=\n", "text/html": [ "\n", " \n", " " ], "text/plain": [ "" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from IPython.display import YouTubeVideo\n", "YouTubeVideo('Adr9_2LnQdw') " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Vapor Pressure of Pure Components" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The calculations in this notebook are for a representative system of two components, acetone and ethanol. We start by creating two functions to estimate vapor pressure for the individual species using Antoine's equation." ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [], "source": [ "# Antoine's equations\n", "A = 'acetone'\n", "B = 'ethanol'\n", "\n", "def PsatA(T):\n", " return 10**(7.02447 - 1161.0/(T + 224))\n", "\n", "def PsatB(T):\n", " return 10**(8.04494 - 1554.3/(T + 222.65))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Create plots for the vapor pressure of each component as a function of temperature." ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnYAAAGKCAYAAACB7AwGAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xd8VFX6x/HPQ28RFVgURIrCAhYUUFEUccVeEVFsiILY\nf66yFmygIIJiQ9YCCohIkaIoCqwo2BusIFJEVPqqNCEBQknO749zB4dhkkzITCaZfN+v17ySueeW\nZ860Z84951xzziEiIiIixV+pZAcgIiIiIvGhxE5EREQkRSixExEREUkRSuxEREREUoQSOxEREZEU\nocROREREJEUosRMpIsxsfzObbGabzcyZ2SH52HaUmU1LZHwiIWbWycx+MbMsM3sl2fEUJWb2mZm9\nlOw4iivVX8EpsSvCgi/5b3Ioq2BmG8ysb2HHFavgDeqC2w4zW2pmj5lZxWTHVkTdChwHnAQcDKyJ\nXMHM2uY36SsoM1sV9jxuM7NFZvYvM0vq54eZ1TWzl81smZltN7PVZjbNzC5MZlxFlZl1MbNdcdhP\nWeBV4A2gDnBXQfcZb2Y2NEg6b9zH7QvyPrsQuGdfjptHTKMKkkSbWbew93Hk7Zh4xhpjPL3NbGmU\nooTUX0mixK5oGwIcZ2bNopR1AKoCSf21bF7ZXFYZiU9SGgJ9gB5A/1z2l9u+ksrMyiX4EA2B+c65\nH5xzvznnshN8vPx4DP88HgEMBZ4E7ijIDgvyXJtZc2Au0AL4J3AkcAYwBRhkZlUKEpvkqjZQCXjf\nObfGObc5coUYPhcSxszSgE5AP6B7YR/fObchWp0UEdvx7+PI2w/JDCpcEa+/4sE5p1sRveET7+XA\n4Chls4CpYfevBr4BNgPr8F9wh4eVHw444EpgJrAN+BnoGLHfJsBUIANIB94BGoSVdwMygXb4L9ad\nwBk5xP8Z8FLEsuHAiuD/dkFM5wCfB/u9ISg7DvggiOMPYAJQJ2w/dYBJwWMNPZa7wsovCeLbCvwJ\nfA0cHXHcgyJic8DVEfV1BTAt2M9jQVkj4K1gvxuB6cAReTyX5YAn8K1wO/AfpJeHla8Kjhe6zYiy\nj8Mj1tm9HjAqiPOm4DWzGXgbqBGxj7OAL4M6W41veTkwj9hXAfdFLJsJfBL+mogorxfEd3JBn+so\n8RiwIHh+y0Qp3w8oHfxfFZ+Irg2O+Q1wepQ67RTEsBVYBJwcvMamAVuC5+uksO1Cj+c8YHaw7/lA\n24hYTgI+Dep7Y/A8VQ8r7wssxr9efwzqYCZwWMR+8no/5LqfsHjDb68EZW2AL/Dv981BvbbLoe67\nRdnPyeTyuQBcF8S2A1gJPBp6fsI+J4bgE7G1+PfVo/jPv0eA34PH/GiMn5s34d/vFYFNQIsojyET\nOIW/PiO+Da1H7u8zA+4Ffg0ez8/A7bl97gX3XwZ6BY9lA/5zsHLEdlcB84LYfgUGApXCykeFnrP8\nPm85vU9zWO+f+Pf8Vvx3QRfCPi+j7Ye93++l8Y0OP/PX53NfoFwur6MHc6i/vD47ywTb34RvRc4I\nXmf3RMSY43dCqt2SHoBueTxB8DD+C6Fi2LKGwQu5fdiyrvgvmcOA5sB7+A/TskF56MNqNT5Z+Tvw\nOJAFHBWsUyl4Q/8n2EdL4BP8F0VoP92Cbb4B2gbHq55D7NESuxeA34L/Q182C4Hzgfr41oCj8F+m\nDwONgaOBicHjCX0wvB/EeUzwofIPoFNQVhv/xXJXsM8m+A/NIyKOG0titzKor/rBcQ7Gf8kMDuJs\nDLyI/0Kqlsvz+Aw+Cb0Unxg+BGQTJAJADfyX9UfAQcABUfZROvhwcsCx4evhP/T/DP4eAbTGJ3jD\nw7Y/E/+hdmvw+I4Pnt+P8ngNRkvs3ge+CntNxJrY5fu5jhJPy2BfnWJ4/7wF/BI89ibB87YdaBjx\nPC/FnwJqBLwbPOYPgYv4K5FfRpBIhj2eH4Fzgab4L+sMoGawTq3g/uv4FsU2+IR0Zlh8fYN1puLf\nc8cA30WsE8v7Idf94L8c7wB2Ba+bg/AJcFl88vNkUBcN8a+x1jnUZ0WgVfDYzw32U5YcPheC+svC\nn1prhH8v/Qn0ivic2IRP7BoBNwT7n4r/jGoEXB8si/ojMiLG/wI3B/8PBYZElIdi/Rj/PmmC/yxZ\nin+P5fY+uwP/HuoW1NUt+NfTtTl97gX3/8Qnan8Hzo5SB92A9fgf6A2AU/EJTPj7d3dil9/nLaf3\naZR1OuA/O+8Iey7+IP+JXVn8GZoTgrKL8UntQ2Gvo4H4BDb0eqycQ/3l9dkZSuz+h/8ePDyI3wGn\nxvKdkGq3pAegWx5PkH9B7gI6hy0bgP/1sldrRdg6NYIX9gnB/dAXWK+I9b4JfXgAN+K/HA4MKz84\n+OC6Mrgf+qV1Ygyx736D4n/ptg4+0EYFy0JfjldEbDcqtE7YsopBHOcH9xcQ/MKLctzjgv0ekkN5\nfhK7nhHr9AU+i1hm+C/923I4Xhr+l2b3iOXvAv+JeNzT8qjTttEeW7Dtb4QlQ8ADwMqI56NvxHYN\ngv0dmcsxdyd2+C+9C4PH0zfsNRFrYpfv5zpKPFcG+8r11zb+S9QBZ0Y8V/MIvuzDnufbwtY5MVh2\nR5TXVOOIx3Nt2Dplg7rqFdx/HJ9clw1bp0Ww3Ulhr6edhP0owH/hZPHXj6lY3g+x7KcLsCtiP6HP\niZPzej+HbROqs1Zhy6J+LuBbh0dHLOuBT1RDSfJnwOyIdX4EvotYtgDon0dsx+FbvPYP7p+Mb9FK\nixLr0WHLWgfLQi2cbYn+Pvsf0C9i2fPAkoj3WWRi99+IbYYCn0a8x7pFrPMPfAKTFuVx7svzFnrc\nGRG3dWHrfAW8FrHds+Qzscvh+HcDi8Lu9waWRlkv/Hsjz89O/krsno5Y5yegT8T7N+p3Qqrd1Meu\niHPOrca3vt0Au/sldQGGOed2d4Q2s+Zm9nbQkTwd/0sIoG7ELr+MuP85voWH4O8PzrkNYcf/H/4N\nckTYNtn400+x6GpmGfgP24/xpy0j+2ZFDhA5DuhoZhmhG75FrCz+lyn4X3EPm9lXZtbfzE4O2/6/\nwAxgkZlNMrP/K8Bgg2ixnRARWzr+tF3Dvbb2GgaxfxKx/GP2rNeCWuic2xF2fw1QM+x+S+BfEbF/\nHxZjbh4J1t8GvAkMw/8iz699ea4jWYzHCtXtp6EFzn/Kf8re9T4v7P/fgr/fR1n2t4jtdr+fnHM7\n8af0wt9PXwbLQ+vMwSc14cdf6ZxbH3Z/Df40ZI3gfqx1lNd+9uKcWwuMAGaY2ftmdq+Z5fVayEm0\nz4WmRH/dV8K3nITMi1jnN/as/9CyyPqP1B2Y7Jz7E8A59xm+pejKiPV2sWe/stBApZrkwMwOxLcs\nRXs8h5lZ+Vzimhtxf/d708wOxv+AHxTxHL+Lf60fHrmzAjxvO/CtueG348PKm+JP74b7LIb97sXM\nbjKzb83sj+Dx9GHv76O85OezM8c6Jr7fCUVemWQHIDEZAkwxsyb4N151wgZNBJ2F/4PvU9MF/wFY\nGv/BlYgO/zvDv6zyMB5/CmkHsCY8GQ2zJeJ+KfyH1pNR1l0H4Jx7xcym4k9rtAWmm9l451wX51yW\nmZ2J/8BqB1wGDDCzS5xzU/FfQBCWIOTS0TtabNPx/VAibcphH4VlR8R9x54DpErhP1zHRNn2tyjL\nwg3C9xPaBvzP7TmwI5u9k6381OcIcnmuo/gx+NuUvb/891X469nlsiwRP4ajPW/hx4q1jvLaT1TO\nuevM7Gn86eozgD5mdrNz7tU84o6Un8+FvbaNDCuHZTk+luBz8AqgYsTo31L4hO/lsGVZEa/hRD6/\nkPt7M/T3NvZOYMB3B9nLPj5vzjkXbSRqfuT5fjezK4Dn8P0RP8X3AeyE72eYKDnWcQzfCSlFiV3x\nMBVYgW+1a4Jvgl4WVt4UqAbc75z7CcDM2uSwr1b4JDDkJHy/J/CnOq43swNDrXbBr8mG7Puoqc37\n8EEyG3+aJNftgtbMV4FXzWw68LqZ3eKc2xq0zHwd3B4zsxn4pHcqvs8I+D5Q/wv+b56P2K7At45s\nj3Gbn/BfUm3w/aJCTiX/9Rr68Cqdz+0A5uD7lOzLB/v6XLb7AyhnZtXCWozyU595PtcR5uAHONxn\nZhMifywEX/Bb8a9n8J3k/xOUWXD/q3wcLzetgCXBvsviW9dCP7oWAFeaWdlQwmNmLYDK5O9535c6\nimYHUMrMLHh/7Oacm48f/PFUMKVGd/x7q6AW4l/34fOSnYp/fn6NusW+uwr/w6NVxPIawEdm1iJo\nMY3FXu8z59wGM/sN/3jC54w8FX9KMdbPg0hr8J9DjZxzw/KzYQKet4X474TwJLh1xDqxvN/b4E+v\nPxtaYGb1I9bZQd6fY3H77MzjOyGl6FRsMRD8snwV34H4THwLXrhl+DfJ/5lZAzM7A3gqh911N7PL\nzayRmT2GPz33TFD2Or4P3FgzO9bMWgJjg/1PiONDystjwFFmNtLMjgse0z/M7HkzqwtgZi+Y2dlm\ndpiZHQG0B5Y557aa2Slm9oCZHW9mhwb1cSR/JbA/4vu0PGJmfzezU/AdeWMxCKgAvG1mJ5tZveBv\nPzM7IdoGzrl0fKf9fmbWIaj7h/CDXfrls26WE4zGNLO/mdl++dj2IaCDmQ00s2PM7HAzO8fMhlvB\npnL5Ct8SNyC0T+DBGLfN87mOFHxAd8b36/nKzC4ys4Zm1sTMbsa34lV0zv2IH/TwkpmdEbR4P48f\ngBDr852X+4M6bIJPXg7gryTmefwPrmFmdkTwOnsNmOWci+wSkZt811EOfsW3tJxvZjXMrErw+n/c\nzFqbnxfwJPwX+cLcdxWzx4HLzeye4DnqhH8dPpFD631BdAcmOT9dUPhtJv4UeX7mtMvpffY48E8z\n6xo8npuD4+b3fbxb8Hp+ALjTzHoGr5W/m1l7M3sx2jYFed7M7KAot9Bp5KfwP0ZuDx5fV/Y+jR3L\n+/1H4BgzuyBY5y78QJpwvwK1gs/p6hZlftN4fXbG8J2QWpLdyU+32G78NYgi6qAJfNPyUnxftv/i\nOw079h4McBW+f0ImfrTg5RH7CU13sgXfsTbqdCcxxrzXqNiI8qiDGIKyZsGx/8T/Cl+K/xUZ6hT9\nEv7X3Db8aLIpQJOg7KjgMfyO72C+DD/gJLwT+4lBPW3D980IdZ6OrK9WUWKrB4zG93MK7f91oG4u\njzVyyP4CIkZ1EsPgiWC9nsF+soiY7iRivS7s3Vn+VPzI24zgOV6IT+xL53K8vUbFRlnnQvwv6m3B\n834O0QdP5Pu5zuWY9fGd0JcHdboG3zIXPhVCaLqTdcFzldN0J+EDAeqFxx4sOyRY1jbi8VwQvI62\n41sQ/hER40lBfeQ63UnENm2J6OidVx3lYz/PB69bh29ZrI1PflcHj2F1sN/9cqn3nAZPRP1cwP8g\nDU13soro051Ejp6fRdjUHsGyGcCIHI4R6hx/eg7lPQgGUUSLNYfnPNr7LNp0J/8Xsa9ogyciH19v\nIgYO4Ee1fhU8v5vxo5pzGiC2L89btClGQreLw9a7K3jc2/Dvp+uIeO+S9/u9XPD62oDvojKKYFR2\n2D7K4xsONpC/6U72+Ozkr8ETkZ+nu19DxPCdkEo3Cx60pDgzOxyfCJ3onIvXaSiREsnM2uHnlTvY\nOZdX/0SRYkuv9eJHp2JFREREUoQSOxEREZEUoVOxIiIiIilCLXYiIiIiKUKJnYiIiEiKKLETFFev\nXt3Vq1cvocfYsmULlStXTugxShrVaXypPuNPdRpfqs/4U53GV2HV55w5c9Y553K8RGBIiU3s6tWr\nx+zZsV7udN/MmjWLtm3bJvQYJY3qNL5Un/GnOo0v1Wf8qU7jq7Dq08yWx7KeTsWKiIiIpAgldiIi\nIiIpQomdiIiISIpQYiciIiKSIkrs4Imc7Ny5k1WrVpGZmVngfVWtWpVFixbFIarUU6FCBQ455BDK\nli2b7FBERERShhK7CKtWrSItLY169ephZgXaV3p6OmlpaXGKLHU451i/fj2rVq2ifv36yQ5HREQk\nZehUbITMzEyqVatW4KROcmZmVKtWLS6toiIiIvIXJXZRKKlLPNWxiIhI/CmxExEREUkRSuxERERE\nUkShJXZmdpuZzTaz7WY2IqKskpm9YGbrzGyTmX0SVmZmNsDM1ge3ARZ2Hs/M6pnZTDPbamaLzaxd\nYT2mVNC7d2+uvvrqZIchIiIicVCYLXZrgL7AsChlQ4ADgSbB3zvDyroDFwPNgKOBC4Abw8rHAN8B\n1YAHgAlmludFckVEREQKYtMmeO+9g3Eu2ZH8pdASO+fcJOfc28D68OVm1hi4EOjunFvrnMtyzs0J\nW+Va4Cnn3Crn3GpgINAl2LYR0Bzo5Zzb5pybCHwPdEj8I0qO/v37c9hhh5GWlkbTpk156623dpcN\nHTqUJk2a7C7773//C8CaNWvo0KEDNWrUoH79+gwaNAiAadOm0a9fP8aNG0eVKlVo1qzZ7vUvvPBC\nDjzwQA4//HCGDh26+xi9e/fmsssuo3PnzqSlpXHEEUcwe/bs3eU5HUtERCSVpKfDuefCwIF/5+uv\nkx3NX4rCPHbHA8uBR8zsGuB/QO8gSQM4ApgXtv68YFmo7BfnXHoO5Xsws+74FkBq1qzJrFmz9lqn\natWqpKf/tbuuXSvw00/7lv86VxGzrDzXa9gwm1dfjW3qj1q1ajF16lRq1qzJW2+9xdVXX83cuXP5\n8ssv6dWrF6NHj6Z58+b88ssvlC1blk2bNnHeeedx7rnnMmTIEFavXs1FF11EnTp1aNeuHT169OCX\nX37hlVdeAfzcex07dqRJkyb8+OOPLFmyhIsuuoiDDz6YU089le3bt/POO+8watQoBg0aRJ8+fbj5\n5pv56KOPyM7OzvVYkTIzM6M+B7nJyMjI9zaSM9Vn/KlO40v1GX+q04LLzCzFffcdzbx5+9O16wIy\nM9dSVKq0KCR2hwBHAhOBWsCJwHtmttA5twioAmwKW38zUCXoZxdZFiqvHe1Azrkh+NO+tGzZ0rVt\n23avdRYtWrTHpMJly0KpfWzXzM7OolSp0nmuV7ZsadLSYrsCQ+fOnXf/36VLF5599lkWLlzIG2+8\nwb333kvoMR1zzDEAfP3116xfv57HHnsMgGrVqtG9e3feeecd2rdvT/ny5Slbtuzux7xy5Uq++uor\npk2bRlpaGjVq1OCGG25gwoQJnH/++ZQvX56TTz6ZSy+9FICuXbvywgsvkJaWluexIlWoUIFjjz02\npscdMmvWLKI9b7JvVJ/xpzqNL9Vn/KlOCyYzEy64AObNg8ceg5NOWluk6rMoJHbbgJ1AX+fcLuBj\nM5sJnAksAjKA/cLWrwpkOOecmUWWhcrTiZM33tj3bdPTt8b9yhMjR47k6aefZtmyZYD/5bVu3TpW\nrlzJYYcdttf6y5cvZ82aNey///67l2VlZXHKKadE3f+aNWs48MAD94i7bt26e5xuPeigg3b/X6lS\nJTIzM9m1a1e+jyUiIlKcbN8OHTrAjBnw8MNw//0UmZa6kKKQ2H0fZVl4N8QF+IET3wT3mwXLQmUN\nzCwt7HRsM6AA6VjRtXz5cm644QY+/PBDTjzxREqXLs0xxxyDc446derw888/77VNnTp1qF+/Pj/9\n9FPUfUZOFFyrVi02bNiwx+XQVqxYQe3aURtB83UsERGR4mrnTrj8cnj/fbj3XujdO9kRRVeY052U\nMbMKQGmgtJlVMLMywCfACqBnsE5r4DRgerDpSOAuM6ttZrWBHsAIAOfcEmAu0CvY3yXAUfjTuiln\ny5YtmBk1avhBv8OHD+eHH34AoFu3bgwcOJA5c+bgnGPp0qUsX76c448/nrS0NAYMGMC2bdvIysri\nhx9+4NtvvwV8X8Nly5aRnZ0N+OTspJNOomfPnmRmZvL999/z6quvxjQlSl7HEhERKY527YKrroLJ\nk+Gf/4THH4eiegGlwpzu5EH8adf7gKuD/x90zu0ELgLOxfeXGwp0ds4tDrZ7GXgXmB/cpgTLQjoB\nLYGNwOPApc65tQl/NEnQtGlTevTowYknnkjNmjWZP38+rVu3BqBjx4488MADXHnllaSlpXHxxRez\nYcMGSpcuzZQpU5g7dy7169enevXqdOvWjU2bNu3eDnx/uObNmwMwZswYli1bRq1atWjfvj2PPPJI\n1MEPkfI6loiISHGTlQVdusD48XDLLfD000U3qQMwV5QmXylELVu2dOH9xkIWLVpEkyZN4nKM8NOZ\nsrd9qWt1+o0v1Wf8qU7jS/UZf6rT2GVnQ7duMHy4//vyy3sPqCys+jSzOc65lnmtp0uKiYiIiERw\nDm691Sd111wTPakriopBiCIiIiKFxzm480546SU/YGLYsOKR1IESOxEREZHdnIP77oPnnoP27eH1\n16FMUZhDJEZK7ERERETwSV3PnvDEE3DeeTB2rL9QQXGixE5ERERKvFBSN2CAvwbsxIlQrlyyo8o/\nJXYiIiJSokUmdZMmQfnyyY5q3yixExERkRIrlZI6UGInIiIiJVSqJXWgxC6ltW3blldeeaXEHFdE\nRCRWqZjUgRK7lNG7d++YrucqIiJS0qVqUgdK7ERERKQESeWkDpTYFTtr1qyhQ4cO1KhRg/r16zNo\n0CCmTZtGv379GDduHFWqVKFZs2a711++fDmtW7cmLS2NM888k3Xr1u0u69ixIwcddBBVq1alTZs2\nLFiwYHdZly5duPXWWznvvPNIS0vjhBNO4Oeff95d/sUXX3DcccdRtWpVjjvuOL744ovCqQAREZF9\nlOpJHUAxmks5Sa66ChYv3qdNK2Vnx3YNksaN4Y038lwtOzubCy64gIsuuogxY8awatUq2rVrx4sv\nvsj999/P0qVLGTVq1B7bjB49mqlTp1KnTh3OOeccBg4cSP/+/QE455xzGDZsGOXKlePee+/lqquu\nYu7cubu3HTt2LFOnTqV58+Zce+21PPDAA4wdO5YNGzZw3nnnMWjQIK644grGjx/Peeedx9KlS6lW\nrVr+KklERKQQlISkDtRiV6x8++23rF27locffphy5crRoEEDbrjhBsaOHZvjNtdddx2NGjWiYsWK\nXHbZZXskbtdffz1paWmUL1+e3r17M2/ePDZt2rS7vH379hx//PGUKVNmj6Tvvffeo2HDhlxzzTWU\nKVOGK664gsaNG/Puu+8m7sGLiIjsI+fg7rtTP6kDtdjlLYaWtJxsTU8nLS0tbqEsX76cNWvWsP/+\n++9elpWVxSmnnELdunWjbnPQQQft/r9SpUpkZGTs3u6BBx5g/PjxrF27llJBy+K6deuoWrVqrtuu\nWbNmr+PVrVuX1atXx+FRioiIxI9zcMcd8PzzcMEFMH586iZ1oBa7YqVOnTrUr1+fP//8c/ctPT2d\n999/HzPL175Gjx7N5MmTmTFjBps2bWLZsmUAOOfy3LZWrVosX758j2UrVqygdu3a+YpBREQkkbKz\n4aabfFLXoQNMmJDaSR0osStWjj/+eNLS0hgwYADbtm0jKyuLH374gW+//ZaaNWuybNkysrOzY9pX\neno65cuXp1q1amzdupX7778/5jjOPfdclixZwujRo9m1axfjxo1j4cKFnH/++fv60EREROIqKwu6\ndoUhQ6BTJxg7tnhe+zW/lNgVI6VLl2bKlCnMnTuX+vXrU716dbp168amTZvo2LEjANWqVaN58+Z5\n7qtz587UrVuX2rVr07RpU1q1ahVzHNWqVWPKlCk89dRTVKtWjSeeeIIpU6ZQvXr1fX5sIiIi8bJr\nF3TuDCNG+L+jRkGZEtL5rIQ8zNRRq1YtxowZE7Xss88+2+P+rFmz9rjfpUsXunTpAkCVKlWYPHny\nHuWdO3fe/f+IESP2KGvbti2rVq3aff/kk09mzpw5UeOIPK6IiEhh2bkTrrzSn3bt1g1efjm2CSpS\nRQl6qCIiIpLKtm+HSy/1Sd0tt5S8pA6U2ImIiEgK2LYN2reHd96BO++EwYNLXlIHSuxERESkmNu6\nFS68EKZOhfvug6eegnxOFpEylNiJiIhIsZWR4ScdnjEDeveGfv1KblIHGjwRlXMu3/PCSf7EMl+e\niIhIbjZu9EndV1/5hK5nz2RHlHxqsYtQoUIF1q9fr8QjgZxzrF+/ngoVKiQ7FBERKab++ANOO80n\ndU8/raQuRC12EQ455BBWrVrF2rVrC7yvzMxMJS85qFChAoccckiywxARkWJo5Uo44wxYsgSGDvXT\nmoinxC5C2bJlqV+/flz2NWvWLI499ti47EtERETg55/h9NNh9WoYMwYuvzzZERUtSuxERESkWFiw\nwLfUbdgAb70FupLl3gqtj52Z3WZms81su5mNyGGdh83MmVm7sGVmZgPMbH1wG2BhIxvMrJ6ZzTSz\nrWa2OHxbERERSQ2zZ0ObNrB5s5/WRElddIU5eGIN0BcYFq3QzA4DOgL/iyjqDlwMNAOOBi4Abgwr\nHwN8B1QDHgAmmFmNuEYuIiIiSfPpp/CPf0B2tp/W5LTTkh1R0VVoiZ1zbpJz7m1gfQ6r/Bu4F9gR\nsfxa4Cnn3Crn3GpgINAFwMwaAc2BXs65bc65icD3QIcEPAQREREpZNOmwVlnQcWKMGsWtGqV7IiK\ntiLRx87MOgLbnXPvR5k/7ghgXtj9ecGyUNkvzrn0HMojj9Md3wJIzZo1E36x+oyMjIQfo6RRncaX\n6jP+VKfxpfqMv+JUpx9/XJ2+fZtywAE7ePLJeWzcuI2iFnpRq8+kJ3Zmlgb0A87IYZUqwKaw+5uB\nKkE/u8iyUHntaDtyzg0BhgC0bNnStW3bdt8Dj8GsWbNI9DFKGtVpfKk+4091Gl+qz/grLnU6ciQ8\n+ijUrw8ffliBunVPSHZIURW1+iwKExT3Bl53zi3LoTwD2C/sflUgw/kZhCPLQuXpiIiISLH07LNw\n7bXQpInvX1e3brIjKj6KQmJ3OvB/Zvabmf0G1AHeNLN7g/IF+IETIc2CZaGyBkGrX7RyERERKSac\ngwcegDvv9H3pPvkEDj442VEVL4V2KtbMygTHKw2UNrMKwC58Ylc2bNVvgbuAqcH9kcBdZvZ+cL8H\nMAjAObcrIbljAAAgAElEQVTEzOYCvczsQeBc4ChgYoIfjoiIiMRRVhbccgsMGQJnnw0TJkDlysmO\nqvgpzD52DwK9wu5fDTzinOsdvpKZZQEbnXMZwaKXgQbA/OD+K8GykE7ACGAjsAK41DlX8OuBiYiI\nSKHYvh2uvtonc506wWuvQblyyY6qeCq0xC5I4HrHsF69iPsOuCe4RVt/GdC2gOGJiIhIEqSnQ/v2\n8OGHcOutMGgQlCoKHcWKKVWdiIiIJMW6dX7i4Q8/hN694fnnldQVVNKnOxEREZGSZ8UKOPNMWLIE\nBg/2rXVScErsREREpFAtWuSTut9+gzfegCuuSHZEqUOJnYiIiBSab76Bc8+Fbdvg3Xf9CFiJH53J\nFhERkUIxfbrvU5edDTNmKKlLBCV2IiIiknAjR8L558P++/uJh088MdkRpSYldiIiIpIwzkH//v4S\nYY0awZdfwpFHJjuq1KU+diIiIpIQWVlwxx3w73/DKafA5MlwwAHJjiq1KbETERGRuMvM9FeTmDgR\nOnSAUaOgQoVkR5X6dCpWRERE4mrjRj+dycSJfn66ceOU1BUWJXYiIiISNytXwsknw6efwuOP+6tJ\nlC6d7KhKDp2KFRERkbj44Qc/hcnvv8Nrr0HnzsmOqORRYiciIiIF9vHHcNFFsGsXTJkCZ52V7IhK\nJp2KFRERkQIZP973qStf3id4SuqSR4mdiIiI7BPn4Mkn4bLL4NBD4YsvoEWLZEdVsulUrIiIiOTb\nrl1w++3w0kvQujW8/TZUr57sqESJnYiIiORLRgZcfjm8/z507OgvF6bpTIoGnYoVERGRmK1ZA23a\n+KTunntg7FgldUWJWuxEREQkJvPnw7nn+uTuxRfhppuSHZFEUmInIiIiefrgA39psOxsePddn+BJ\n0aNTsSIiIpKrYcN8IpeW5q8ooaSu6FJiJyIiIlE5Bw8+CF27QpMm8NVXcOyxyY5KcqNTsSIiIrKX\nzEzo1g3eeAPOOAMmTID99kt2VJIXtdiJiIjIHv74A04/3Sd1XbvCe+8pqSsulNiJiIjIbgsWwAkn\nwJdfwhNPwNChULZssqOSWOlUrIiIiAAwbZq/PFhWFrz1Flx0UbIjkvxSi52IiIgweDCcd54/5frZ\nZ0rqiisldiIiIiXYrl1w663+uq/Nm8M332jka3FWaImdmd1mZrPNbLuZjQhb3srMPjCzDWa21szG\nm9nBYeVmZgPMbH1wG2BmFlZez8xmmtlWM1tsZu0K6zGJiIgUZ3/+6VvpXngBLr0UPv4YatVKdlRS\nEIXZYrcG6AsMi1h+ADAEqAfUBdKB4WHl3YGLgWbA0cAFwI1h5WOA74BqwAPABDOrEf/wRUREUscv\nv8BJJ8F//gMPPADjxkGlSsmOSgqq0BI759wk59zbwPqI5VOdc+Odc5udc1uBwUDrsFWuBZ5yzq1y\nzq0GBgJdAMysEdAc6OWc2+acmwh8D3RI/CMSEREpnj77DI4/Hn7+GUaOhL59oZQ6Z6WEojgqtg2w\nIOz+EcC8sPvzgmWhsl+cc+k5lO/BzLrjWwCpWbMms2bNilPI0WVkZCT8GCWN6jS+VJ/xpzqNL9Vn\n/L399v688EI2lSrtYuDAH6hTZzOq4n1X1F6jRSqxM7OjgYeB8LE4VYBNYfc3A1WCfnaRZaHy2tH2\n75wbgj/tS8uWLV3btm3jE3gOZs2aRaKPUdKoTuNL9Rl/qtP4Un3Gz65dcPfd8NxzcMQR8M475WjQ\noHmywyr2itprtMgkdmZ2ODAVuMM592lYUQYQPt91VSDDOefMLLIsVJ6OiIiIALBxI1x+OXzwAbRu\nvY6pU6uTlpbsqCQRisQZdTOrC8wA+jjnXo8oXoAfOBHSjL9O1S4AGphZWg7lIiIiJdqiRb4/3Qcf\nwIMPwqOP/qCkLoUV5nQnZcysAlAaKG1mFYJltYGPgMHOuZeibDoSuMvMagfr9gBGADjnlgBzgV7B\n/i4BjgImFsJDEhERKdLee89fHmz1aj/qtU8fDZJIdYV5KvZBoFfY/auBRwAHNAB6m1nvUKFzrkrw\n78tB+fzg/ivBspBO+ERvI7ACuNQ5tzb+4YuIiBQPzvnrvPbsCYccApMna9LhkqLQEjvnXG+gdw7F\nj+SynQPuCW7RypcBbQsUnIiISIrYtg26dYPRo6F1a5g4EWrWTHZUUljUICsiIpIiVq2CNm18Ute1\nK3z0kZK6kkaJnYiISAr46is47jj47jsYNAiGDoVy5ZIdlRQ2JXYiIiLF3JAhvqVu+3aYPh1uvx3+\nuqq6lCRK7ERERIqp7dvhhhvgxhuhaVOYPRtOPz3ZUUkyFZkJikVERCR2q1bBpZfC11/DlVf6U6+V\nKiU7Kkk2tdiJiIgUM598Ai1a+Ba6p5+GUaOU1ImnxE5ERKSYcA6ef96fbs3O9leTuPNO9aeTv+hU\nrIiISDGwbZvvS/f66761btIkOPTQZEclRY1a7ERERIq4Zcv8ZMOvvw5dusCnnyqpk+iU2ImIiBRh\nH34ILVvC/PnwwgswbBhUrJjsqKSoUmInIiJSBGVnw+OPw5lnQtmyMHMm3Hyz+tNJ7tTHTkREpIjZ\nuBE6d4YpU/wp2DffhFq1kh2VFAdqsRMRESlCvvvOD46YMgXuusu31Cmpk1gpsRMRESkiXn0VTjwR\n1q2D8ePhqaf8aViRWOV4KtbMPolxH5nOuTPjFI+IiEiJs20b3HorDB8ORx4JEydCo0bJjkqKo9z6\n2B0H3JTH9gY8F79wRERESpalS/2lwebNg6uvhpdegsqVkx2VFFe5JXZfOOdey2sHZnZlHOMREREp\nMSZPhmuv9S12L77oJyDWqFcpiBz72DnnTo9lBzoNKyIikj+7dsF998HFF0PVqvDZZ3DTTUrqpOA0\n3YmIiEghWr0arrjCXz3i7LNh1CioVi3ZUUmqiCmxM7OVgItStB1YBUwCXnTO7YpjbCIiIill2jS4\n5hrYsAH69oWePaGU5qeQOIq1xW4QcHXwdyVwKHArMB7YAPQA6gD3JCBGERGRYm3XLnjoIejf389J\nN3MmtGmT7KgkFcWa2HUBznDOrQktMLOpwH+cc0eY2UxgBkrsRERE9rBypT/1+vnncNZZ8PrrUKNG\nsqOSVBVrA/DBQEbEsi1AaC7sJcD+8QpKREQkFbz/Phx7LHz5JfTr5+8rqZNEijWxexeYbGbtzKyx\nmbUDJgbLAU4EliUgPhERkWJn506491447zyoUAFmzVJ/OikcsZ6KvRHoDbyMb6X7H/Am8GhQ/gtw\nXryDExERKW5WrPCnXr/4wo96HTlSrXRSeGJK7JxzmcB9wS1a+W/xDEpERKQ4mjLFTzi8aZMfKHH3\n3Wqlk8KVa2JnZg3y2oFz7pf4hSMiIlL8ZGbCPffA88/DIYfAO+9A69bJjkpKorxa7Jbi568LzYUd\n/n/ofukExCUiIlIsLFoEnTrB99/7K0m88oomHJbkybWB2DlXyjlXOvhbCtgc+j9UFuuBzOw2M5tt\nZtvNbERE2elmttjMtprZTDOrG1ZmZjbAzNYHtwFmf110xczqBdtsDfbRLvaHLyIism+cg6FDoUUL\nWLIEXngBJk1SUifJld8z/9GuPhGrNUBfYFj4QjOrjr9yxUPAgcBsYFzYKt2Bi4FmwNHABfjBHCFj\ngO+AasADwAQzUzdVERFJmI0b4bLLoHt3aNAAvv0Wbr5Z13qV5Cu0Lp3OuUnOubeB9RFFlwALnHPj\ng0EavYFmZtY4KL8WeMo5t8o5txoYiJ8wGTNrBDQHejnntjnnJgLfAx0S/oBERKRE+vxzOOYYmDDB\nJ3PffgtHHpnsqES8WKc7SaQjgHmhO865LWa2NFi+OLI8+P+IsG1/cc6l51C+BzPrjm8BpGbNmsya\nNStODyG6jIyMhB+jpFGdxpfqM/5Up/FVlOozKwveeKMur71Wj8qVd9Gnz4+cfPI6vv462ZHlT1Gq\n01RQ1Oozr1Gxr7Pn6dfKZjYyfB3nXOcCxlAFWBuxbDOQFla+KaKsStDPLrIsVF472oGcc0OAIQAt\nW7Z0bdu2LVDgeZk1axaJPkZJozqNL9Vn/KlO46uo1OfKlXD11fDJJ/4ar6NGlaVOneLZTFdU6jRV\nFLX6jGVUbLh+CYghA9gvYllVID2H8qpAhnPOmVle24qIiBTI+PFw442weTP06eOvIFFa80FIEZVr\nYuece6QQYliA70cHgJlVBg4LlofKmwHfBPebRZQ1MLO0sNOxzYA3Eh20iIikts2b4fbb/ZUj6tWD\nd9/V3HRS9MXcxy6YgqQZ/vTnbs650TFuXyY4XmmgtJlVAHYBbwFPmlkH4D2gFzDPObc42HQkcJeZ\nvR/c7wEMCo69xMzmAr3M7EHgXOAo/HVsRURE9slnn8E118CyZf5KEoMGwX6R54dEiqCYEjsz64mf\njmQhsC2syAExJXbAg/ikLeRq4BHnXO8gqRsMjAK+BjqFrfcy0ACYH9x/JVgW0gkYAWwEVgCXOuci\n++yJiIjkaccOeOQRfzmwqlXhzTehY8dkRyUSu1hb7HoALZ1zC/f1QM653vipTKKVzQAa51DmgHuC\nW7TyZUDbfY1LREQEYPFiP0Bizhxo1w5GjIDaUYfiiRRdsc5jtx5YlsA4REREksI5ePFFaN4cfvgB\nnnkGpk9XUifFU6wtdv8EhpjZs8Af4QXOuRVxj0pERKQQ/P47dO0K770HRx0Fb7zh/4oUV7EmduWA\nM4ErI5Y7/GAIERGRYmXyZLjhBli7Fnr0gL59oUKFZEclUjCxnop9AbgfP2dc2bBbuQTFJSIikhB/\n/ulHul58MZQvDzNmwMCBSuokNcTaYlcGGO6cy0pkMCIiIon0wQdw/fWwahV07gzPPQf775/sqETi\nJ9YWu4HAfcFlvERERIqVLVvg1lvhzDP9lCZvvQWvvaakTlJPrC12/wccBNxvZuvDC5xzh8Y9KhER\nkTj5/HN/6vXnn+GSS+Cll6BGjWRHJZIYsSZ2Vyc0ChERkTjLzISHH/b956pWhVGj4MorQeeeJJXF\nlNg55z5OdCAiIiLx8t//+j50CxbAWWfBq69qXjopGWK9pFgZ4ArgWPa+Vmz3BMQlIiKSbzt2wOOP\n+6lLypeHl1/2U5qolU5KilhPxY4CjgKmAr8nLhwREZF989130KULfP89nHKKvyRYgwbJjkqkcMWa\n2J0N1HHOpScyGBERkfzavh369IH+/X0r3aBBfgRsqVjnfRBJIbEmdguAAwEldiIiUmR8+y1cd53v\nS3faafDKK2qlk5It1sTuGuAVM/sPEadinXMj4x6ViIhILjIzoXdvePJJqFQJXnwRundXK51IrIld\nF+AU4ABgW9hyByixExGRQvPll/7qEYsXwxlnwNChULdusqMSKRpiTezuAI51zi1KZDAiIiI52boV\nHnoInnkG0tJ8Qte1q0a8ioSLNbH7HViRyEBERERy8vHHftqSn36Cc87x05jUqZPsqESKnlh7IzwD\njDKzVmbWIPyWyOBERKRky8goQ/fu0LYtrFsHw4fDe+8pqRPJSawtdv8O/l4UsdwBpeMXjoiICDgH\nkyZB9+7HsWEDXHYZPPccHHRQsiMTKdpivaSYxhmJiEihWL0abrsN3n4batRwvPMOXHBBsqMSKR6U\nsImISJGQnQ0vvQRNm8LkyX6S4eHDv1VSJ5IPOSZ2ZvZaLDsws+HxC0dEREqixYvh1FPh5puhdm34\n7DMYPBgqV85KdmgixUpup2IvNbMRQF4DydsD18UtIhERKTF27IAnnvCXBHPOTzp8333+0mAikn+5\nJXZrgWEx7OO3OMUiIiIlyKefwk03wcKFcOKJ/nJgTZsmOyqR4i3HxM45V68Q4xARkRJi3Tq45x4/\ndUnVqvDvf/sET5cDEym4WKc7ERERKRDn4LXX4F//gvXroVMnfxUJTWEiEj/6fSQiIgm3aJGfZPi6\n62D//WH6dBgzRkmdSLwpsRMRkYTZtg0efBCaNYMvv/T/z58PZ56Z7MhEUlORSezMrJ6ZvW9mG83s\nNzMbbGZlgrLTzWyxmW01s5lmVjdsOzOzAWa2PrgNMNMloUVEkm36dDjySHjsMTjpJPj+ez/6tWLF\nZEcmEidTp9Lszjth69ZkR7JbzImdmTU2s4fM7N9h94+OYywv4EfiHgwcA5wK3GJm1YFJwEPAgcBs\nYFzYdt2Bi4FmwNHABcCNcYxLRETyYdUquPxyOPts2LzZ96ubORMaN052ZCJxsnWrn0H73HOptGIF\nLF2a7Ih2iymxM7OOwCdAbeCaYHEV4Ok4xlIfGOecy3TO/QZMA44ALgEWOOfGO+cygd5AMzMLfURc\nCzzlnFvlnFsNDAS6xDEuERGJQWhOusaN4c03oWtXP/Fw586g8yiSMr79Fo49Fl54AS65hNnDhsHR\n8WznKhhzzuW9ktkioJNzbp6ZbXTOHWBmZYE1zrkacQnE7EbgJOBm4ABgOr6Vri1Qzjl3c9i684He\nzrmJZrYJONM593VQ1gKY5ZxLi3KM7vgWPmrWrNli7Nix8Qg9RxkZGVSpUiWhxyhpVKfxpfqMv5Ja\np3Pm7M+gQQ1ZsaIyDRum889//kTTppsLvN+SWp+JpDrdN5aVxaFvvEHdkSPJLl+en26/nd/POouM\nLVsKpT5PO+20Oc65lnmu6JzL8was568kcEPwtwzwRyzbx3iMJsAcYBfggBH4q168CvSPWPdzoEvw\nfxbQOKysYbC95Xa8Fi1auESbOXNmwo9R0qhO40v1GX8lrU5XrnTussucA+cOOMC5F15wbteu+O2/\npNVnYVCd7oOffnKuVSv/Qj/5ZOd+/XV3UWHVJzDbxZBPxdrHbg5/nYIN6QR8E+P2uTKzUvhTr5OA\nykB1fKvdACAD2C9ik6pAevB/ZHlVICOoBBERSYBop11//NFf67V06WRHJxInzsHQoXDMMTBnDvTv\nD7NmQb16yY4sR7FOUPx/wH/MrCtQ2cymA42AeA1YPxA4FBjsnNsObDez4UBfYBC+Hx0AZlYZOAxY\nECxagB84EUoym4WViYhInH34Idx2m+8/17y5v3JEq1bJjkokzv74A264Ad55x1/rbtQo37euiIup\nxc45txhoDPwbeBAYDhzlnPspHkE459YBvwI3mVkZM9sfn8x9D7wFHGlmHcysAtALmBfEBDASuMvM\naptZbaAH/jSuiIjE0YoVfrRru3bw++++7/g33yipkxT07rtw1FE+qfvnP2H27GKR1EEMLXZmVhpY\nAjR1zr2ZwFguAZ4F7sP3m/sIuNM5t9bMOgCDgVHA1/jTwCEvAw2A+cH9V4JlIiISB1u3wpNPwoAB\nfsLh66/3Z6RqxGXonEgRsnGjT+RGjoTateGDD/wvmWIkz8TOOZdlZllARWB7ogJxzs3Fj4CNVjYD\n32IYrcwB9wQ3ERGJE+dg/Hi4+27fWteqFQwaBMcdl+zIRBLg/ff9qdc1a/wcPc8+CwcckOyo8i3W\nwRPPAuPM7FQzO8zMGoRuiQxORESSY948OO00f+p11y54/XX44gsldZKCNm3yzdDnnQfZ2TB5sp9V\nuxgmdRD74InBwd8zIpY7QOOfRERSxLp1/nquQ4dC2bJw//3Qsydo2jNJSdOnQ7du/nIpV17pm6Sr\nVUt2VAUSU2LnnCsy15QVEZH427nTD4bo3Rv+/BPat4eBA6GBzstIKtq8Gf71L/8L5m9/g0mT/Is+\nBcTaYiciIilq+nS46y5YuBCOPBImTIDTT092VCIJMmOGn3hxxQro1Amefx6qV092VHETU2JnZp/i\nT7vuxTnXJq4RiYhIoViwwDdaTJvmuxMNHgw33ghl9JNfUtHmzXDvvfDSSz6RGz8eLr002VHFXaxv\n31ci7h8EdMVPPyIiIsXI779Dr17+LFSpUnDnnb5f3YEHJjsykQR57z246Sbfl65DB9/v4G9/S3ZU\nCRFrH7vXIpeZ2UT8RMWPxjsoERGJv23b/AwOjz8O6em+S9GAAdCwYbIjE0mQtWvhjjtgzBg46CDf\nz6BDh2RHlVAFGRSxGjg6XoGIiEhiOOe/1xo39qNcGzb0l7ucNElJnaQo5/wlwJo08S/+66/3nUhT\nPKmD2PvYXR+xqBL+ShFfxT0iERGJm88/9wMjvvnGT6Q/ciRcdZU/BSuSkpYv96ddp03zw7rHjStR\no4Fi7WN3TcT9LcAXwDPxDUdEROLhp59869yECVC5Mjz6KPToAZUqJTsykQTJyvJ953r29P0OevTw\nL/wS9qKPtY/daYkORERECu733/132ZAhfhL966+HPn2gVq1kRyaSQAsX+omGv/wSjj4aXnmlxF4m\nJabGeDNramY1g/+rmNkjZtbLzEpWGiwiUkRlZPiE7vDDfaPF2Wf7y4K9+qqSOklhmZl+Vu1jj4U5\nc6BvX5g9u8QmdRD7qdgxwGXA78BA4O9AJvAye5+mFRGRQrJzp2+ceOQR31p3/PHwxBNw6qnJjkwk\nwWbMgFtu8f0OTj7Zz9/TuHGyo0q6WBO7es65H83M8IMmmgLbgF8TFpmIiOTIOXjrLd+daMkS31L3\n/PN+vlWzZEcnkkC//eb7z40e7a/rOmwYXHutRgQFYq2FTDNLA44HVjjn1gHbgQoJi0xERKL67DNo\n3drP3PDnn/Dvf/suRh07KqmTFJadDS++6FvlRo+G666DxYv9XyV1u8XaYjca+AhIAwYHy5qjFjsR\nkUIzdy488AC8/74f6dqrl2+4SEtLdmQiCTZ3rp/C5Ouv/dx077wDbXRF02hiHRV7p5mdCex0zs0M\nFmcDdyYsMhERAfyp1l69YOxYKFvWdyt66CE/kb5ISktP9y/+556DcuWgXz//a6ZcuWRHVmTFfKln\n59x/Qv+bWQNgnXNudkKiEhERVq70I12HD/dnoTp39gMA69dPdmQiCRbqRHrHHf76ruecA4MH+wmH\nJVexTncyxsxOCv6/DlgALDCzrokMTkSkJFq71l8tomFDP+L1wgth/nx47TUldVIC/Pijn6+nQwf/\ni2b8eHjvPSV1MYq1t+HpQKh17i6gHX4gxX2JCEpEpCTavNmfdWrQAJ55xnch+uYbf03XI45IdnQi\nCZaeDvfeC0cdBR99BP/6FyxapKHe+RTrqdhyzrkdZlYbONA59zlAaNJiERHZd1u2+JGtAwbAhg3Q\nqpXvSnSarvkjJYFzMGYM3H03rFkD7drBoEF+kITkW6yJ3Vwz6wnUBd4DCJK8zYkKTEQk1W3dCi+9\n5BO6P/7wDRXDh8MFF6iBQkqI77+H22+HTz6BOnX8xY0vuURvgAKI9VRsV+AooCLwYLDsROCNRAQl\nIpLKMjN9g8Rhh/kBftWqwZtv+hkdLrxQ32lSAvz5px8Y0bw5fPUVPPign5OuQwe9AQoo1ulOfgau\njFg2AZiQiKBERFLR9u1+MES/fv6MU6NG8PTTcNllULp0sqMTKQTZ2TBiBNx3nx8ldP75vkPp4Ycn\nO7KUEeuoWDOzG8zsQzP7PljWxswuS2x4IiLF344d8PLLfpTrbbdBxYowciQsWABXXKGkTkqITz/1\nFzPu2hX22w/efdfflNTFVaynYh/Fn44dChwaLFsF3JuIoEREUsGOHb6FrlEjP2l+6dL+spaLF8M1\n10CZmGcSFSnGfv3VN0u3aeOnMnn8cfjhB99aJ3EX68dKF+BY59w6M3sxWPYroEllREQibN/uE7j+\n/WHFCjj0UBgyxF+nXBPmS4mxebNP4p55xv/K6dYN+vTRJVMSLNbErjSQEfzvgr9VwpaJiJR427bB\nxIm1ufpqWL0a6tXzp2CvvRbKl092dCKFJCvLD+9+8EH4/Xdo29Ynd8cck+zISoRYT8VOBZ42s/Lg\n+9wBfYB34xmMmXUys0VmtsXMfjazU4Llp5vZYjPbamYzzaxu2DZmZgPMbH1wGxDEJyJSKLZsgaee\n8leFGDy4IRUr+u+1JUuge3cldVKCzJoFLVrADTdAlSr+smAffaSkrhDFmtjdCRwMbAKq4lvq6hLH\nPnZmdgYwALgOSAPaAL+YWXVgEvAQcCD+ChjjwjbtDlwMNAOOBi4AboxXXCIiOUlP96db69Xzk+Qf\ncADcf/9CFi2CLl2gbNlkRyhSSJYu9fPPnXaa71P35JN+dNDFF2v6kkKWZ2IXtH5VBzriB060Ag5z\nzrV3zqXHMZZHgEedc18557Kdc6udc6uBS4AFzrnxzrlMoDfQzMwaB9tdCzzlnFsVrD8Q3ydQRCQh\n/vzTdxWqVw969vRdhsaN8/3BzzjjDw2KkJLjjz/8BMNNmsDkyXDjjfDTT/6Xjpqqk8Kcc3mvZLYF\nSHPOZSckCLPSwDbgYaAbUAF4G7gb6I+/pNnNYevPB3o75yaa2SbgTOfc10FZC2CWcy4tynG641v4\nqFmzZouxY8cm4uHslpGRQZUqVRJ6jJJGdRpfqs/82bChHOPHH8I779Ri69YyHH54Op07L6d163WU\nCn4mq07jS/UZf/Go01LbtnHIxIkcOmYMZbZuZf0JJ/BL9+5saVDyxlQW1mv0tNNOm+Oca5nnis65\nPG/AZ0DjWNbdlxtQCz8oYzb+lG914HPgMeBVoH/E+p8DXYL/s8JjAxoG+7LcjtmiRQuXaDNnzkz4\nMUoa1Wl8qT5j8/PPzt10k3PlyzsHzrVu7dyUKc5lZ++9ruo0vlSf8VegOt2507mhQ507+GD/ZmjZ\n0rmPPopbbMVRYb1Ggdkuhpwq1hMGs4BpZjYCWMlfI2Nxzg2LcR+52Rb8fd459z8AM3saf/myT4D9\nItavCoROA2dElFcFMoJKEBHZZ/Pn+z5048b5gX7nnONPvZ5ySrIjEylkzsGUKf6KEQsX+pFCY8b4\n+elKxdpdXwpDrIlda/y8dadGLHdAgRM759xGM1tFWMIY9v8CfD86AMysMnBYsDxU3gz4JrjfLKxM\nRCTfvvzST7/17ru+33fHjv777Nhjkx2ZSBJ88w3ccw98/LG/sPGzz/oZt9WHrkiK9VqxpyU6EGA4\ncO8J1isAACAASURBVLuZTQN24kfiTgHeAp40sw7Ae0AvYJ5zbnGw3UjgLjN7P7jfAxhUCPGKSApx\nDqZP9y10H3/sR7R26+a/zxo2THZ0Iknw44/w0EMwfjxUqOCbq++9F6pWTXZkkotcEzszq4Q/HXok\n8F/gcefc9gTF0gfft24JkAm8CTzmnMsMkrrBwCjga6BT2HYv46+AMT+4/0qwTEQkTzt2wNixMHCg\nP/VauTLcdZe/1a6d7OhEkmD5cnj0URgxwt+/7jp//5BDkhqWxCavFrt/Ay3xExRfClQDbk9EIM65\nncAtwS2ybAbQeK+NfJkD7gluIiIx2bzZX+br2Wf9VSL+9jc/hcnNN/uzTSIlzu+/Q79+8NJL/hfP\npZf6hK5Jk2RHJvmQV2J3NtDcOfc/M3seP5AhIYmdiEhhWL0annvOX+pr82Zo1MgneNdc4882iZQ4\nGzf6Juv/b+/Ow6Os7v6Pv78Ji0AIi0AAhQAqokLZXUDZBFFakbrUDXfFpWpb7dOnfVp/tbXtr621\ni1urj0tFUdzQSkWsWxTUIohQRHBhR7aAEExCWJLz/PGdcSaRNZlklnxe13Wuydz3neTMIWE+OfdZ\n/vxnKC2FU0+FX/3Kd5CQtLOvYNcsOks1hLDKzHRjXUTS0ocf+nvX44/Dzp0waBD813/B2LGa1Cf1\nVEkJ3Hkn/P73vur24MHeYzdkSLJrJjWwr2DXwMyGA7aH54QQXq+tyomI1EQI8Oqrvv/4Sy/5DNdx\n43xR/EGDkl07keTI2rED7rrLe+U2bPB9XH/zG++p0/ZfaW9fwW4DlZcz2VTlecAnLoiIpIyyMpg0\nye8sffih32K9+mqfENG9e7JrJ5Ik27fDAw9w3K23wsaNcOSRcPfdcNZZ6rbOIHsNdiGELnVUDxGR\nGlu/Hu69F/76Vygs9D1cb7vNl9xq0ybZtRNJku3b4cEHfXHG1aup6NgRHn4Yxo9HGxtnHv2Likja\nmz/fb7c+8YRP5uvbF+64wxfF1xqqUm9t3w4PPeS3WVevhm7d4OGHea9TJ4aefHKyaye1RMFORNJS\nRQW8+KIHujfe8KFBY8fCD37gY781VEjqre3bvUfuN7+BVas80D30kPfQNWxIKChIdg2lFinYiUha\n2bLF36PuuQeWLvUFhW+4AW68EQ4/PNm1E0miHTs80P361x7ounatFOikflCwE5G0sHChT+R79FFf\naqtrV1++5IoroGXLZNdOJIlKS+GBB+D22/2Wa9euPqbuoosU6OohBTsRSVm7dsHUqR7o3njDj40a\n5T10Y8ZAdnZy6yeSVFu3+myhP/7RZwsdfrgCnSjYiUjq2bTJOyDuvRdWroScHLj+evjud6HHbjcX\nFKlHNm3y7VPuusvHJvTs6c/POUezXEXBTkRSx+zZvlTJE0/4WnTdu/vC+JdcArm5ya6dSJKtXevT\nvf/2N981YuBA+OlP4fTTtQ6dfEXBTkSSqqQEJk/2QPf++z6b9bTT/HbrKafo/UqE5ct926+HHvIZ\nr0OHeqAbOVLTv+VrFOxEJCkWLfKOh0cegaIiaNsWfvxjmDDBx36L1HsffOATIp56CsrL/S+en/7U\n93QV2QMFOxGpMzt2wPPPe+9cdCmtk06Ca6+FM8/UYsIiX21wfPvt8Mor3mV9zjnwox9Bv37Jrp2k\nAQU7Eal1y5b5ZL0HHvBtv3JzfTLE1Vf7uG+Rem/XLnj6ab/lOm8eNGnis4VuuskXGBbZTwp2IlIr\nduyAF16A++/3DogQoE8f37v1/PN9pqtIvVdS4n/1/PGPsGIFHHww3HqrhzptcCzVoGAnIgn1ySfe\nM/f3v/vSWs2a+SLCEybAgAEa6y0CeNf1Pfd4+eIL75W75x649FJo2jTZtZM0pmAnIjVWVgZTpnjv\n3Jtv+rGBA31no/POg+bNk1s/kZQxf75vcPzEE96t3b+/Dzo980ytQScJoZ8iEam2BQv8LtKjj3qn\nQ24uXHcdXHWV33YVEaCiAv75T/jzn30LFTMYOxa+/31fukTd2JJACnYickC2bPHOhocegjlz/Nig\nQX6r9ZxzdBdJ5CvFxT4m4S9/gc8+84GlN97oizQefniyaycZSsFORPapogJef93D3JQpvkZqXh78\n8Idw2WVw9NHJrqFIClm50rf7+t//9UUa8/PhD3/wwaYtWya7dpLhFOxEZI+WLfMOh7//3d+rGjSA\nb30LLr8cTj1V+4yLfCUEeOstuPtueO45X1B48GC/3TpunMbPSZ3RT5qIVFJc7O9Lf/+799IBHHOM\nb1E5fjy0a5fU6omkluJimDTJA92HH3qA+853PNAde2yyayf1kIKdiFBe7jtBTJwIzz7rS2u1aAHX\nXOO9c1qmRKSKTz6Be++Fhx+GrVuhY0f4xS985lCHDsmundRjCnYi9diKFU35yU/gscdg9WrIzobR\no+Hii33SXpMmya6hSAopL4dp07x37l//8mNDhvg2KuPGaWyCpAQFO5F6prAQJk/23rk5c/xWUZ8+\nvnPR+edD+/ZJrqBIqiks9J65v/4Vli/3qd9XX+27Q/TqlezaiVSiYCdSD5SW+jJakyZ5h8OuXR7g\nzj13Jf/zP535xjeSXUORFBOCr7Z9330+FXzHDl+i5E9/8t0hNLtVUlRWsisQz8yOMLMyM3ss7tjJ\nZrbYzErN7A0zy487Z2b2OzPbFCm/M9NIIBHw8Pbyy3DJJb40ybnnwiuv+Lju6dNh1Sq45pqlCnUi\n8TZt8n1bjzoKhg+HZ56B00/3W68ff+yTIhTqJIWlWo/dPcDs6BMzawNMAa4EpgK3AU8Cx0cumQCM\nA3oDAXgFWAb8re6qLJI6QoBZs+Dxx+HJJ2HDBh83N2oUXHghnHGGtvcS+ZoQYOZM75175hlfqLFL\nF98T77LLNBlC0krKBDszOw/YArwDRJfkPhNYGEJ4OnLNrcBGM+sRQlgMXALcEUJYHTn/BzzsKdhJ\nvbJ4sYe5xx+HJUv82AknwC23eA+dligR2Y0vvvCZQ/fdBx995H8FfetbPh181Ch/LpJmLISQ7Dpg\nZrnAHGAE3jt3eAhhvJn9BWgUQrg27toFwK0hhGfNrAg4JYQwK3KuP1AQQthtn4SZTcCDH3l5ef0n\nT55cq6+ruLiYnJycWv0e9Y3aNGbNmoN44412vPFGO5Ys8TbJzy9h5Mj1jBixgY4dy/b5NdSeiac2\nTayEt2dFBa3mzqX9Sy/RdsYMsnbupKxdO9aOGcPaMWPY0bZt4r5XitLPaGLVVXsOHz78/RDCgH1d\nlyo9drcBD4YQVlcZIpcDFFa5divQPO58UZVzOWZmYTeJNYRwP3A/wIABA8KwYcMSU/s9KCgooLa/\nR31T39t01Sp46im/zTo7Mmihc2ff2uvCC6F372aYdQO67dfXq+/tWRvUpomVsPZcscJntj78cGwb\nlbFj4fLLOejUU+manU3Xmn+XtKCf0cRKtfZMerAzsz7ASKDvbk4XA7lVjrUAvtzD+RZA8e5CnUi6\nWrsWnn7aw9w77/ixDh3ge9/zCRHHHQdZKTUNSiRFlJXB88/Dgw/Ca6/5WLqjj9Y2KpLRkh7sgGFA\nF2BlpLcuB8g2s6PxsXKXRC80s2bAYcDCyKGF+MSJ9yLPe8edE0lb69b5tl5PPunbT4YAbdvCtdd6\nmDvxRA3/EdmtEOCDD7xnbtIk2LzZZwxdeSVccYVv86XFEySDpUKwux+IH+z2QzzoRcfV3W5mZwEv\nAj8H5kcmTgBMBG4ys2mR5zcDd9Z6jUVqwerVvlzWs8/CjBn+/tS6tb8XnXsuDBumfcRF9mjNGp8I\nMXEiLIz8fT9kiO+Jd/bZ0KxZcusnUkeS/jYRQigFSqPPzawYKAshFEaenwXcDTwGzALOi/v0+/DB\nRAsizx+IHBNJCytWeJB75hl4910/1qaNh7mzz/ZltBo1Sm4dRVJWSYnfap04EV59FSoqID8ffvYz\n3xfviCOSXUOROpf0YFdVCOHWKs9fBXrs4doA/ChSRNLCkiUe5J55BubM8WN5eX6b9eyzvZNBPXMi\ne1BR4TtCTJzov0TFxX6r9dJLPcyddJIGnUq9prcPkVoWAsyb52Pmnn8eFkT6lw85BG680cPcoEEa\nMyeyVwsX+kKNjz3ms1qzsuCUUzzMnXGG798qIgp2IrVh1y54++1YmFuxwo8fdhjcfDOcdZZms4rs\n04oVdHriCZ8C/p//+LFeveAPf4ALLtCOECK7oWAnkiDbtvkwn+eeg6lTYeNGP96vn4+Z+/a34Zhj\nNCFPZK8KC/0W6+OPw8yZHAa+WON//zecfz707p3sGoqkNAU7kRooLIQXX4QXXvA9wktKvBduyBAY\nN85Lfn6yaymS4r78Ev7xDw9z//oXlJf7LKLrrmNujx70++531b0tsp8U7EQOQAi+peTUqV7efdeP\nNWniW0uOGwenn+7vSSKyFyUlMG2ar779z396l3ezZt4rd8EFMHIkNGzI1oIChTqRA6BgJ7IPO3f6\nunIvvOBhbulSP96hA1x1lQe5k0/2cCcie1FaWjnMlZb6ej6nnuph7vTTNQlCpIYU7ER2o7AQpk/3\n26zTp0NRZEfiPn3gllt8i8l+/dSRILJP27bBSy/5JsdTp8bC3OjR8J3veJhr0SLZtRTJGAp2IvjS\nWHPnemfCiy/C7Nl+i7VRI18keOxY+Na3fAy3iOxDSQm8/LL3zE2d6s8bNoyFubFjFeZEaomCndRb\nRUXwyise5F56Cdav9+OHHgoTJsCYMX6LVTsRieyHLVv89uqUKd7NvW2bh7lTTomFuZYtk11LkYyn\nYCf1Rggwf753JEyfDjNn+npz2dm+QPAPfuBhrmdPLUkisl/WrfPZrFOmwOuv+y/UQQf5mLkzz/Ru\n7latkl1LkXpFwU4yWmGh98q9/LKvorBunR9v29bHao8Z4x0Keu8R2U/Ll/tijVOm+CrcIUBuLpx7\nroe50aPVzS2SRAp2klF27YJZs7xHbvp0eP99f99p0MB75W64wTsT+vTRxAeR/RIdgPqPf/jU8OgO\nEG3bwpVXepgbMcIHpIpI0inYSVoLAT77zHd8eOUVvxsUncGan+9j5U491d93cnOTW1eRtFFW5r9M\n0TV+1qzx4926wfe/79uoDB6sDY5FUpCCnaSdjRvhtddiYS66D2vTpjB0qAe50aOhe3eNlRPZb/Hb\nqLz8si9LYgbHH+9d3WPHwlFH6ZdKJMUp2EnK27bNh/K88go891x/Pv3Uj2dlwYABMH687/pw/PHQ\nuHFy6yqSNkKAefM8zE2bBv/+d+VtVM44A775TcjLS3ZNReQAKNhJytm509eRe/11L++8A9u3+7lD\nDsnm2mt9t6HhwzXpQeSAbN3qXd3TpnlZu9aPd+wIV1zhvXInn6zdH0TSmIKdJF15uS9DEg1yb73l\n65kCtGvnw3lGjPBOhOXL32PYsGFJra9I2ggBPv441is3Y4b/5ZSVBSecANdf71PDe/fWLVaRDKFg\nJ3WuogIWLoSCAnjjDX/cvNnPtWzpAW7ECC9HH135/Wb58iRUWCSdFBX5X0gvv+wl+ktz8MG+UPA3\nv+lr/Bx8cFKrKSK1Q8FOal1Fha+Q8OabHuJmzIBNm/xcs2YwZIiHuOHDfRkSTbQTOQDl5b6uT3Sx\nxnff9WNm0L+/D0L95jdh4ED9conUAwp2knDl5T4mu6DAw9yMGb7bEHiQGzwYhg3zGawDB/quQyJy\nAD7/3EPcyy/7mLnoX0odOniQGz3au77btEluPUWkzinYSY2Vlflkh5kzPcS9/baP0QZo3hxOOslD\n3NCh0K+fgpzIASsq8r+Uouv8LFrkxxs18l+w0aO99OqlsXIi9ZyCnRywLVs8vEWD3OzZsGOHn2vV\nym+tRnvk+vTxXR9E5ABs3+7TwaNBbvZsH9MAvpnx977n4+SGDtX2XSJSid5yZa9CgFWr/D1mxgwP\ncwsW+HGATp3g7LO90+DEE32yg7bqEjlA5eXwwQc+6eG11/yXbds2P9epE1x6qS9DMmIEtG+f1KqK\nSGpTsJNKduzw8XHvvBMrn38eO3/MMb5N10kneencOXl1FUlb8QNRCwp8jZ/o+IVWrXwJkpEjPcwd\nfrhur4rIflOwq+c2bvRJdNEQ9957PmYOfAH6Y4+Fiy+GQYN82SutkCBSDeXlPjU8ur7PW2/FNjXO\nzfW/koYP9zEMmhouIjWgYFeP7NjhCwHPmuW7B82aBZ99FjvfqZPvIjRokJfevTXRQaRadu6EuXP9\nluqMGR7kolPDozOK4oOcBqKKSILof5MMFQKsXBkLcf/+t7/PRLfmatQI+vb15a2ivXGdOiW3ziJp\nq7TUf8neeoveL7zguz2Ulvq55s19AGo0yPXtqyAnIrVG/7tkiC++gDlz/Fbq7Nn+uG5d7HzXrnDm\nmXDccXD88d5J0Lhx8uorktY2bozNKJoxwxcI3rULgGatWsFpp8UGovburVurIlJnUiLYmVlj4F5g\nJNAaWAL8JITwUuT8ycA9QGdgFnBpCGFF5JwBvwWujHy5B4AfhxCdt5l5Skp8Al00wM2eDUuWxM43\nb+4L/15+uYe4447zPVdFpBoqKmDxYg9yb7/tj598EjvftStccIGHuCFDeOfzzxk2fHjy6isi9VpK\nBDu8HquAocBKYAzwlJn1AoqBKXhwmwrcBjwJHB/53AnAOKA3EIBXgGXA3+qw/rWmrMzHXL//vpfZ\ns+HDD2NLWjVu7L1vp53mEx0GDoTu3bXkiEi1FRf7X0zRGUXvvhsbH5ed7b9wN9zgYxhOPBEOPbTy\n569ZU/d1FhGJSIlgF0IoAW6NO/RPM1sG9AcOBhaGEJ4GMLNbgY1m1iOEsBi4BLgjhLA6cv4PeNhL\nu2C3bVvlEPf++7Bw4Vd3eMjK8uVGLr3UA9yxx/papY0aJbXaIumrosLHw8XPKFqwwGexgi89Mnhw\nbEbRwIFaEFhEUpql4h1LM8sDVgB9gGuBRiGEa+POLwBuDSE8a2ZFwCkhhFmRc/2BghBC89183Ql4\n6CMvL6//5MmTa/V1FBcXk5OTs9tzpaXZLFmSw6efevnkk+YsX96MigpfryorK5CfX0L37sV07/4l\n3bt/yWGHFdOkSUWt1jnV7a1N5cDVt/ZsuGULuYsW0XzRInIjpUFJCQDBjNL8fLYedRRFPXuytWdP\nSg899IC7v+tbm9Y2tWfiqU0Tq67ac/jw4e+HEAbs67qU6LGLZ2YNgUnAIyGExWaWAxRWuWwrEA1u\nOUBRlXM5ZmZVx9mFEO4H7gcYMGBAGDZsWC28gpiCggKGDRvGunU+Jm7ePH/84IPKy4xkZ/uODRdf\nDP37e+nd22jaNCfy8rTSfFS0TSUxMro9S0p8KvicOT6GYdYsWLo0dr5dO9/JITIQ1QYOpFluLs2A\nDjX4thndpkmg9kw8tWlipVp7plSwM7Ms4FFgB3B95HAxkFvl0hbAl3s43wIoTvbkiTvugKee6sXK\nlZVnpzZp4pPkRo3yoTp9+/rt1CZNkldXkbQXHYw6e3YsyC1aVHkwar9+MHZsbGp4fr52dBCRjJMy\nwS4yu/VBIA8YE0LYGTm1EB9HF72uGXBY5Hj0fG/gvcjz3nHnkubNN2Hx4lyOPRYuusgDXJ8+PrFB\nKx+I1EBZmc8gmjvXB6LOmeOhLjoYNTsbevWCK67wMXEDBvhfT1ptW0TqgZQJdsBfgaOAkSGEbXHH\nnwNuN7OzgBeBnwPzIxMnACYCN5nZtMjzm4E766jOe/TEE/Dee28zfPiwZFdFJH0VF/t2KXPnxspH\nH8VCnBn06OHLjURDXO/e6gIXkXorJYKdmeUDVwPbgXUWuz1ydQhhUiTU3Q08hq9jd17cp98HdAMW\nRJ4/EDmWVM2a6S6PyAHZsMFD3Lx5XubO9Rmr0VEV2dk+LXz8eO8C79fPQ1zzr82TEhGpt1Ii2EUW\nG95jDAohvAr02MO5APwoUkQk1e3aBZ9+6uFt/vxYmIsfjNqoEXzjG3DVVR7g+vXz26sHHZS8eouI\npIGUCHYikqEKC31duGj5z3/8sawsdk3btt7zNn68P/bu7bdXNSZOROSAKdiJSM2VlvrYt/gQt2AB\nrF8fuyYry2cPjRsXC3B9+kD79hq3ICKSIAp2IrL/ysp83NvChR7kFi708tlnsbFwAB07emjr2dNv\nofbqBUcdpUkNIiK1TMFORL5u27avB7iPPoIlS2JrwwHk5vqEhuHDYwGuZ084+ODk1V1EpB5TsBOp\nxxoWFcHMmb6Y7+LFscflyyv3wEUD3NCh/nj00f7YsaNuo4qIpBAFO5FMt3MnLFsGn3zivXCLF38V\n4gZv2lT52rZtfeLCqFH+eMwxCnAiImlEwU4kE4QAa9Z4eIuWjz/2x6VLobw8dm1WFnTtCscfz8qc\nHDpHQ1yPHrqFKiKS5hTsRNJFRQWsXesTFT791B/jS0lJ5evbtPFZqIMHw5FH+sdHHOElsh7c0oIC\nOqfQ5tUiIlIzCnYiqWTnTli50nvZli6tHNyWLPFJDfFycjyonXZaLLxFS+vWyXkNIiKSNAp2InVt\nyxYPbUuWxAJc9OOVKyvfNgVo0cLD2xlnwOGHVy7t2mnsm4iIfEXBTiTRtm71WaXLlvljtESfFxVV\nvt4MDj0UunXzZUO6dYuVww7zcW8KbyIish8U7EQORAiwcSOsWOFl5crYxytWeHDbvPnrn9eqlU9Y\nOPlk6NKlcnDLz4fGjev6lYiISAZSsBOJV1oKq1bB6tX+uGqVh7dogFu58uvj3ABatoTOnWHIEA9w\nXbp46drVg1uLFnX9SkREpB5SsJP6o6QEPv/cQ9vnn8eCW3yI++KLr3+ema/jlp8Pffv6Y+fO/hj9\nODe37l+PiIhIFQp2kv7Ky6Gw0NdxW7MmFtyiJfq86ti2qLZtfYzbiSdCp05fLx07QqNGdfuaRERE\nqkHBTlLXrl0e2Nat88C2di35b78Nkyf7em7RILd+/ddnkoIvxNuhAxxyCIwY4eHtkENipVMnPxZZ\n001ERCTdKdhJ3QrBl/tYv97LunWxsnZt5Y8LCyvvVwp0jX7Qpo33pPXqBaNH+8cdOvjjIYd4YMvL\ng+zsOn+JIiIiyaJgJzW3c6eHsMJC2LDBy/r1scf4smGDX787DRp4GGvfHgYO9KDWvr0f69gROnbk\n3RUrOOHb39atURERkd1QsJOvKyvzkLZxY6xEg1s0vMU/7m55j6jsbF9ENy8PevaMfRwt7dt76dDB\nd0rIytpr1baXlSnUiYiI7IGCXSYLwZfv2LQpVjZurPw8eiy+VN1ztKrWrX3CQXxYa9u28mNenj/u\nR1gTERGRxFCwSwcVFb6bwZYtvhzH5s3+uK+yaRNs3773r924sYexNm1826oTTog9r1qixxs2rJvX\nLSIiIgdEwa4uhOCL2m7ZEitFRZWfb97sZXcfFxV5uNubrCzvHWvVyh87dfKtqPZVmjbVdlUiIiIZ\nQsGutlx5JcdOn+49ZkVFe54wUFWLFh7OWrXyxW/79Ik9b9nSw1g0vMWX5s11y1NERKSeU7CrLU2b\nsr1tW5p26eJhrWVLL7v7OBrmWrTQ8hwiIiJSbQp2teXOO5lfUMCwYcOSXRMRERGpJ3TvTkRERCRD\nKNiJiIiIZAgFOxEREZEMkRHBzsxam9lzZlZiZivM7IJk10lERESkrmXK5Il7gB1AHtAHeNHM5ocQ\nFia3WiIiIiJ1J+177MysGXAWcEsIoTiEMBP4B3BRcmsmIiIiUrcshJDsOtSImfUF3g4hNI07djMw\nLIRwepVrJwATAPLy8vpPnjy5VutWXFxMTk5OrX6P+kZtmlhqz8RTmyaW2jPx1KaJVVftOXz48PdD\nCAP2dV0m3IrNAbZWObYVaF71whDC/cD9AAMGDAi1vcZcgdaxSzi1aWKpPRNPbZpYas/EU5smVqq1\nZ9rfigWKgdwqx1oAXyahLiIiIiJJkwnB7hOggZkdEXesN6CJEyIiIlKvpH2wCyGUAFOAX5pZMzM7\nERgLPJrcmomIiIjUrbQPdhHXAU2ADcDjwLVa6kRERETqm0yYPEEI4QtgXLLrISIiIpJMab/cSXWZ\nWSGwopa/TRtgYy1/j/pGbZpYas/EU5smltoz8dSmiVVX7ZkfQmi7r4vqbbCrC2Y2Z3/WnJH9pzZN\nLLVn4qlNE0vtmXhq08RKtfbMlDF2IiIiIvWegp2IiIhIhlCwq133J7sCGUhtmlhqz8RTmyaW2jPx\n1KaJlVLtqTF2IiIiIhlCPXYiIiIiGULBTkRERCRDKNjVkJk1NrMHzWyFmX1pZvPM7LS48yeb2WIz\nKzWzN8wsP5n1TRdm9piZrTOzrWb2iZldGXdObVpNZnaEmZWZ2WNxx9Se1WBmBZG2LI6Uj+POqU2r\nwczOM7NFZlZiZkvM7KTIcbXnAYr7uYyWcjO7K+682rQazKyLmU0zs82R96i7zaxB5FxKtKmCXc01\nAFYBQ4EWwM+ApyL/+G3wfWxvAVoDc4Ank1XRNPNboFsIIRff+/dXZtZfbVpj9wCzo0/UnjV2fQgh\nJ1KOBLVpdZnZKOB3wGVAc2AIsFTtWT1xP5c5QHtgG/A06Ge0hu4FCoEOQB/8vf+6VGpTTZ6oBWb2\nH+AXwMHApSGEQZHjzfDVqfuGEBYnsYppxcyOBAqA7wEtUZtWi5mdB5wJfAQcHkIYb2YTUHtWi5kV\nAI+FEB6oclxtWg1m9g7wYAjhwSrH1Z41ZGaXAD8HDgshBLVp9ZnZIuDmEMK0yPPbgVzgfVKkTdVj\nl2Bmlgd0BxYCxwDzo+dCCCXAZ5Hjsg9mdq+ZlQKLgbXANNSm1WJmucAvgZuqnFJ71sz/N7ONZva2\nmQ2LHFObHiAzywYGAG3N7DMzWx25xdUEtWciXAJMDLGeHLVp9f0ZONfMmprZIcBpwHRSqE0V7BLI\nzBoCk4BHIgk9ByiqctlW/DaD7EMI4Tq8rU7Cu7i3ozatrtvw3pDVVY6rPavvv4FuwCH4OlZTf6Zd\nKQAABotJREFUzeww1KbVkQc0BM7Gf9/7AH3xoS1qzxqIjPMaCjwSd1htWn1vAT3x9lqN33J9nhRq\nUwW7BDGzLOBRYAdwfeRwMd5FG68F8GUdVi2thRDKQwgzgUOBa1GbHjAz6wOMBP60m9Nqz2oKIcwK\nIXwZQtgeQngEeBsYg9q0OrZFHu8KIawNIWwE/ojaMxEuAmaGEJbFHVObVkPkfX463tHQDGgDtMLH\nhqZMmyrYJYCZGfAg/lfnWSGEnZFTC4Hecdc1Aw6LHJcD04BY26lND8wwoAuw0szWAT8EzjKzuag9\nEykAhtr0gIUQNuO9H/GDvqMfqz1r5mIq99aB2rS6WgOdgbsjf9BtAh7G/wBJmTZVsEuMvwJHAaeH\nELbFHX8O6GlmZ5nZQfjg1fkanLp3ZtYusuxBjpllm9lo4HzgNdSm1XE//h9Mn0j5G/AiMBq1Z7WY\nWUszG21mB5lZAzO7EJ/FOR21aXU9DNwQ+f1vBfwA+Cdqz2ozs0H4UIGnq5xSm1ZDpCd5GXBN5Pe+\nJT5+8T+kUpuGEFRqUIB8/C/LMrwrNloujJwfiQ/+34bP7OyS7DqnegHaAm8CW/AxCguAq+LOq01r\n1r634rM51Z7Vb8O2+LIxX0Z+Tv8NjFKb1qhNG+JLSWwB1gF3AgepPWvUpvcBj+7hnNq0em3aJ9Je\nm/FZr08BeanUplruRERERCRD6FasiIiISIZQsBMRERHJEAp2IiIiIhlCwU5EREQkQyjYiYiIiGQI\nBTsRERGRDKFgJyKSQczsGjPbZWbFZtZtPz+nReT6nWb2s9quo4jUHgU7EUm6SKiIlgoz2xb3/MJk\n168mzGydmZ1Yx9+2IISQE0JYGlePo81sipltMrMtZjbPzG40MwshFIUQcoBn67ieIpJgCnYiknSR\nEJITCRcr8e35oscmJbt+e2JmDdLhe5hZD+Bd4GPgmBBCS+AC4ESgcU2/voikDgU7EUl5kT2DbzGz\npWa20cwmRfZpxMx6RG49XmFmn0d6pC43s0Fm9mGkd+qPcV/rGjN73czuM7OtZvaRmQ2JO9/azCZG\netpWmdnPzSyryufeY2abgR9Hvn+BmX1hZoVm9oiZNY9c/zTQDvhXpPfxRjM71cw+q/L6vurVM7Pf\nmtnjZvakmX0JnLe317+ffgW8GkL4SQhhHUAI4aMQwndCCGXV+TcRkdSkYCci6eCHwCl4D9OhwE7g\nT3Hns4FvAN2Ay4C7gJuBYZHjl5nZcXHXDwHmAwcDvwWeN7PcyLlJQFHkax0LjAMuqvK584A2wB2R\nY78E2gO9gCOBnwKEEM4BNgCnRHof79zP13sW8AjQAr89uq/Xvy8jgWcO4HoRSVMKdiKSDq4BfhxC\nWBPpYfoFcK6ZWdw1vwwhbA8hvBB5PjGEsDGEsBJ4B+gbd+2qEMK9IYSdIYSJwGpgtJnl48HtphBC\naQhhLb4Z/Xlxn7s0hPC/IYTyEMK2EMLiEMLrIYQdkd6wPwNDa/h63wwhTAshVIQQtu3n698tM8vG\nA+LaGtZJRNJArY8PERGpiUh46QRMM7MQdyoL73EDKA8hbIo7tw1YX+V5Ttzz1VW+zQqgI5APHAQU\nxmWmLCD+1umqKvXrCPwFGAQ0j1xf0xD11ffYj9e/cW9fKIRQbmZFQIca1klE0oB67EQkpYUQAvA5\nMCKE0DKuHBRC2Guo2YtDqzzvDKzBA1Ux0Cru++SGEPrFV6nK594OlAA9Qwi5wJWA7eX6EqBp9ImZ\nNQRaV7nmq89J0Ot/Fb+9KyIZTsFORNLB34DfmlknADNrZ2an1+DrdYpMhGhgZuPxHrF/hRCWAf8G\nfm9mzc0sy8yO2MdyJc3xMLjVzDoDN1U5vx4frxe1CGhtZidHQt0v2Pf/xTV9/bcAo8zs12aWF/ka\nR0YmaBx0AF9HRFKcgp2IpIPf471Or0dmir4D9Nv7p+zVW/iYuy/wiQ5nhhCKIufOB1oCiyPnnwTy\n9vK1/h8+qaEIeI6vrwX3a+DXkdm510d62b6HT9JYDaxjH7dTqeHrDyEswm8VHw0sMrMtwGRgBrB9\nf7+OiKQ+815+EZH6wcyuAc4OIYxMdl1qg5ldiY/52wn0jfRC7utzWuDjDBsCt4UQflu7tRSR2qLJ\nEyIiGSSE8ADwwAF+ThHeSykiaU63YkVEREQyhG7FioiIiGQI9diJiIiIZAgFOxEREZEMoWAnIiIi\nkiEU7EREREQyhIKdiIiISIZQsBMRERHJEP8H+OiYC4l97WgAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "\n", "# Plot pure component vapor pressures\n", "T = np.linspace(20,80)\n", "\n", "plt.plot(T,PsatA(T),'b')\n", "plt.plot(T,PsatB(T),'r')\n", "\n", "plt.ylabel('Pressure [mmHg]')\n", "plt.xlabel('Temperature [C]')\n", "plt.legend([A,B],loc='best')\n", "plt.title(\"Vapor Pressure of the Pure Components from Antoine's Equations\")\n", "plt.grid();" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## A Short Database of Antoine Equations for Selected Compounts" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": false }, "outputs": [], "source": [ "# Psat by Antoine's equation for selected compounds\n", "Psat = dict()\n", "Psat_Tmin = dict()\n", "Psat_Tmax = dict()\n", "\n", "Psat['benzene'] = lambda T: 10**(6.90565 - 1211.033/(T + 220.790))\n", "Psat_Tmin['benzene'] = 8\n", "Psat_Tmax['benzene'] = 103\n", " \n", "Psat['ethanol'] = lambda T: 10**(8.04494 - 1554.3/(T + 222.65))\n", "Psat_Tmin['ethanol'] = -2\n", "Psat_Tmax['ethanol'] = 100\n", " \n", "Psat['methanol'] = lambda T: \\\n", " (T <= 65) and 10**(7.89750 - 1474.08/(T + 229.13)) or \\\n", " (T > 65) and 10**(7.97328 - 1515.14/(T + 232.85))\n", "Psat_Tmin['methanol'] = -14\n", "Psat_Tmax['methanol'] = 110\n", " \n", "Psat['toluene'] = lambda T: 10**(6.95464 - 1344.8/(T + 219.48))\n", "Psat_Tmin['toluene'] = 6\n", "Psat_Tmax['toluene'] = 137\n", "\n", "Psat['water'] = lambda T: \\\n", " (T <= 60) and 10**(8.10765 - 1750.286/(T + 235.0)) or \\\n", " (T > 60) and 10**(7.96681 - 1668.21/(T + 228.0))\n", "Psat_Tmin['water'] = 0\n", "Psat_Tmax['water'] = 150\n", "\n", "\n", "# Tsat functions\n", "Tsat = dict()\n", "Tsat_Pmin = dict()\n", "Tsat_Pmax = dict()\n", "\n", "from scipy.optimize import brentq as fzero\n", "\n", "for s in Psat.keys():\n", " Tsat_Pmin[s] = Psat[s](Psat_Tmin[s])\n", " Tsat_Pmax[s] = Psat[s](Psat_Tmax[s])\n", " Tsat[s] = lambda P,s = s: fzero(lambda T: Psat[s](T)-P,Psat_Tmin[s],Psat_Tmax[s])" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": false }, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python [conda root]", "language": "python", "name": "conda-root-py" }, "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.5.2" } }, "nbformat": 4, "nbformat_minor": 0 }