{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "\n", "*This notebook contains course material from [CBE20255](https://jckantor.github.io/CBE20255)\n", "by Jeffrey Kantor (jeff at nd.edu); the content is available [on Github](https://github.com/jckantor/CBE20255.git).\n", "The text is released under the [CC-BY-NC-ND-4.0 license](https://creativecommons.org/licenses/by-nc-nd/4.0/legalcode),\n", "and code is released under the [MIT license](https://opensource.org/licenses/MIT).*" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "\n", "< [Products: Product Design and Analysis](http://nbviewer.jupyter.org/github/jckantor/CBE20255/blob/master/notebooks/A.00-Projects-Product-Design-and-Analysis.ipynb) | [Contents](toc.ipynb) | [Pyrotechnic Design for Airbags](http://nbviewer.jupyter.org/github/jckantor/CBE20255/blob/master/notebooks/A.02-Pyrotechnic-Design-for-Airbags.ipynb) >
"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Diesel Engine Emissions Control"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Summary\n",
"\n",
"This [Jupyter notebook](http://jupyter.org/notebook.html) demonstrates an analysis of NOx emissions control in diesel engines."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Background"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Diesel engines offer greater energy efficiency than equivelently sized gasoline engines by operating at higher compression and combustion temperatures. The problem, however, is that diesel engines produce more particulates because the fuel is injected late in the compression cycle, and produces greater NOx emissions because of the higher combustion temperatures. \n",
"\n",
"Here we look into a particular technology for treating NOx that has found wide spread use in heavy trucks and has been recently introduced in the automobile and light truck marketplace. "
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"image/jpeg": "/9j/4AAQSkZJRgABAQAAAQABAAD/2wCEABALDA4MChAODQ4SERATGCgaGBYWGDEjJR0oOjM9PDkz\nODdASFxOQERXRTc4UG1RV19iZ2hnPk1xeXBkeFxlZ2MBERISGBUYLxoaL2NCOEJjY2NjY2NjY2Nj\nY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY//AABEIAWgB4AMBIgACEQED\nEQH/xAAbAAEAAgMBAQAAAAAAAAAAAAAAAwQBAgUGB//EAEYQAAICAgAEBAQDBgIHBQkBAAABAgME\nEQUSITETQVFhBiIycRQjgRUzQlKRsaHBNVNicnPR8BY0Q4LhJCUmNmN1krLxB//EABkBAQADAQEA\nAAAAAAAAAAAAAAABAgMEBf/EACQRAQEAAgICAgICAwAAAAAAAAABAhEDEiExBBNBUSIyI2Fx/9oA\nDAMBAAIRAxEAPwD5+AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA\nAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA\nAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAXVgSwxrrIc8K5Sj6pGY4t09ctcnt6XT\nzPRcJdy4LjOiKlvJ5bF6xZ0MdRruqjFLSzZRX6oDx7wclWKt0zUpdlruZjgZMrZVxpnzRW2tdkel\nrnOvi9TVV0lyTT8To+z7FyqvWZZN878XFbUO0lpgeLsxbq+fng48mubfls1opnfaq64uUn2SPSyx\nLMjH4jXGmxTlySipvrrsUOA1/h+MOmzUbHCcYt+T5XoChPhuVXa651NSUebXsP2dkecGukW/17Ho\n4Rtqlw1ZM140Zzdj7/J0/wDU2yJw/EZlsbIuqyEJ1+yUkB5vK4ZlYmvHrcdvS+5JXwfMsvVKr1Nx\nUkn07nouIVwc71bfBwtyoyrW96WzXOyqLsrFyse7U6b/AA59e8drX+YHnJ8KyoWuvk3JSUenq1sl\nlwTLjk10JKUrN617f/w6+FmVUytdlm9Zalv200Zouq4ZbQp5MbZeLKSkn2Ti0B5zJxLMVQ8Ra509\nfoQHb+IsrHyq8OdH1csnNej2itgcHvzo1OrWrG1/QDmg6ePwa7IuhVHpKVcpv20V7MKVeCsl9vE8\nPQFQFnMxnjSrT/jgpf1Wy5VwuuXC3lzu5ZNtKOgOUDorhrfEoYu+ktNv0Wi5LgdU14lF3PW0+V/7\nS7oDhA7PFuDR4cpfPzNSSXuWHwXGX4Vbnu5Jtvt1jsDzwOv+z8ariONRa24Wxi+j9S/+ycPeVy1T\naps5PqA8yD0N/CsVcI/EKDjNVqfPzdG/Q88+4AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA\nAAAAAABaxuIZOLXKumxxjJ7aRqs2+L6Wy+rn/UrgC1PiGTZNSnbJtdi3wz8Rn5vzXuChFuUt9kjl\nF7hGTXjZMvGT8OyDg2vLYHQ/Czs4tVXTl+JC2O/ET7Ikx+FVPiGRVfdJWVJzTX8UUu5rRm4HDspP\nHj4nJU48z/ibMy43TPNrynUlLwJVSSffo0v7gRy4c7+GXZ8LZOSk1FesV3IbsGrHzMaq21+HdBS5\nn5bLWPx6rFxacaNKnV4coWbffZBmZ+HkYdcXRLxq4ckZ79OwFuPCaaeIU4uRNu1wc5LfT26macLB\nrz76LqrFqLsj19EUp8YjZxSGTZW3FVKtrfto1v4zzZ3jwq+VVOpLfk1oDmXuDun4aahvomRmZPbb\n9TAA9DwLiNeJwvK556sh81S9d9zzwA9nfxfh9WTdkQfPuEIxSTX8XMynfZgZmHk1xyVBSv8AFimn\n5p9P6nmNmdv1A7mbPFsqxrlapWVQjFw130aX8ZjZw50VpV81jlKKXlpf8jjbCA70s7BhfRlQlOVq\njGMotdO3Ukq4zh4lmLCiMpVVzlOe132cjOx40V40o/8Ai18zKYHW4lxZZuHXTrUo2Sk393s1hxVx\nyMWTXy0xS19lo5YA7GRxTHtx6fyfzqtJT9dM3/bdTll89O4X2+Io+hxAB1J8YlLHljqGq3Xya/sc\ntgAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA0Z5X6AYR1OO1QoyaY1xUY\numEmvfRzowm+0W/0OvxTGtswuGuuEpfkNPp/tMDigtR4flS7Uz/obS4ZlwqlZKpqMVt79AKYAAAu\n8KxK8zJlXbNxioOXT2LLXB4v6rZfowOSDq+NwuL2qZy+5t+0sGH0YUX9wOTo2jXKTSUZP9Dovi0d\n7hi1R/Qft3I2nGNcdekQN+I49k8DhzVcm/DafT3KUOHZc+2PPXrotW8ezbFrxOn2IP2plf6xgZ/Z\nGZ51cv3ZIuDXvvKEf1Kss3Il9Vsv6kbusfebf6gdBcFnv5sipfqV8zAniVwm5RlGXTaKjk33bL9r\n5uB1dd6vl/ZAc8AAAAAAAAAAAAABkwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAdDgVUL+L\n41Vq3Cc9NFy/i1dF064YNK5Za6r/ANDn8JsdfFMWa6NWL+5vxuHh8ZzYeStl/cCd8dt/hoqj9kaz\n4/mySXNFKPRLXY5YAuz4tmz73P8AQu8LyrsiGdC2cp7xpNbfocU6nw/82fOt9p0zj/gBzH3MGZdJ\nMwB0/h7/AErCP80ZR/wZzrFy2SXo2i9wKXLxfGe+9iX9WVcxay7lrXzv+4EIAAAAAATUY12Q9U1y\nn9kBFrb0jbwrOVy5JaXsen4J8PWV3RyMuCeuqr9T0tmBjOhq+FVUJ9XHzCZLXzHll6M6ccWx8Dsl\nyy+W5dNex69/sHE/gjN/2M5HFeGW4NlUK4wTkunqRbpbpk+eNNGC3k0p5Nig04cz5X7GqxdruNo6\nVWBc/Cx85G6oq9B2T0qho2UJPsmX1XXHtFEi0l0iVuTTHhl/LnKib/hN1i2PyL+36GNsr3bT4+P7\nU/wc/U3WGl3kWNsxtkdqv9PHPwjWNWl22ZVFa/hNxtjtV/rwn4a+FD+U1lRB+RIB2qLhjfwqyxfR\nkE4OL0zo6NZ1qa6ovMmGXB+nOBPdS4dfIgLSubLG4+wAEqgAAAAAAAAAAAAAAAAAAAAAAAAAA3pk\n4WwknpqSaOj8SL/33ky/mlzHMOr8ROFmfGyuSkpVQb166A5IJIUznJKMX1PQ43w/iRphZlZSTfVx\nQTpQ4LwO3iynKElGMPXzO/w34Yli5NdsbeaSbTWvZl/Eu4ZiVcsbo1QXlHuza34nxcXpjQ5td2/M\nja848r4cOHwflWWydkkuptm/B8sfElZGzmkuyJbvifInNyh0TK8uO5Nso88pOKe+XfcbW+mz3XAw\nvyeJUOfy8lsd78tM24tyvimS4PcXY9NE86vHypTUfmnLaLcOC5VnXwtfci5aU6OGot+Rnw5/ys9F\nDgUv47YRJ1wnFr/eZG/sZ/dit9bzCpm/I3WLL1PTLE4bDp1mWcerDl9FC6eo+2X0dZHlasTr2nJ+\nyOzgX5WJWq8WjXXe2jrSshDpXXCLXsQPIypPSj/REXkv6TJEULOMSe4y5ffsRz4fm3vmychb+5a8\nLOm+smZ/ZuTJpuRllyVrLfwprhdK62ZO2c3ilVVNkI0Tctrqejjwdv8AeS0cbjeNDGyYRg9/L3L8\nee77Z57chQNlAk0NG7GLEcVNbbN/w0TR5MEl8xFLNin6kabeFuNMIo21BehQeamu3+Jr47fmOtO+\nLotw12RVtjFPaZX55epq235jqtObXpKnszoihLRKnspZp0cecyAZBVrpgGdDRIAzoaCdNXFSWmUL\n6/DmdAqZr+ZF8XN8jGa2qgA0cIAAAAAAAAAAAAAAEtWPZb9MQmTaIHRr4fHpzS2/Qswwa0vp2Rtp\nOK1xTblb8mdr8BF9YrRRtg42uPkhsvH191Wrx52PSWvdll4cI4bsdi8XeuVehlQb8zdU+w2p/FCu\nH26Ta0mbR4fKT0tts9Fh8NUqKJZNjUbflrR04cOriq66ZR8XFlzWt+hS5L+P085RwXaj4lc9vy0a\n5HCbY2qMape2z28I2xyb7lYpV2x/Kfv7EfhO7Fh+Ilq2D2/XRS56m0x42HCsvncVW00WFwPMm+u/\n6np53VVXPKcm4SWkRT41VF/Ktsp9mV9La/dcSv4bvl9XQtf9muWqTnLbSb0XZ8am/wB3BET4llT+\nmC/oU+zKJ1L+XAp4XbObrfytf2OhVwylWRm+tWtP7nQVqVbyJ1/mv5WZdVeljpvcvmK3lyqfrjnR\n4fGuE5afiKW4G0YZlnnI6DcHPx1J8tXRr1JPxsZ188I9CcbcvZZI5y4dkz7y/wASWHCJfxyEc++U\npPl16GK8jKvtSe0iLMtollWIcJqj1lMmhj49L+VrqQWynHUW31Iv2dbY+ZT0idWedm5+ks4U1zc5\nS6ElefQ5KMI/4EMeFJ/vLGWasGijrzdSMspaaqDJzXXPl5e5C8rKnPlgnpl2yWKusmmbVZWPLagl\ntE3rSbjnWLNsaW2cbitVlWRFWvbaO/Zmy8dpQ6LzOFxm535fNJa6GvHjq+lc7uOeDVySRo59To0x\n2hs6zf3NUbvq9mNDSGNGdsGSUtq5PzJdb6ohRJFgGiSqXkaNBFMo148tVODEXtGxlp6Mu4AAJAAA\nKecvmg/YuFLM+pGmHtzfJ/qrAA0cAAAAAAAAAAAABbwcV3z2/pQTjjcrqM4WI7JKc1qB1YVL6YLU\nfYRq6qK7LyLUeWC1vqZXJ3cfFJGsaIxRvypIRk5PsTKPQp2b48aDl32OS489km/NndaSRyUkXw9u\nX5WOtRHGGiRRUpJdjJtXHnthFfxPRpXG9dWoO3wJ0fl4sVKD9WSQnj+G7+RqWW+R+wVF9dNNamnN\nPc/sbRhNTyFyxlXBc1X3Ma0jMqaPy6I2NPFfO/sYnNK2d8Hz1XPlT9DVtqNdsqvmvXJM2hGNXi48\ndcsO2/Uzz9L4+0V1FV0lidlDqmI8LxYvcpJ/c0yFdHFjJL8599dyisbNuXXm6meO9L+HU8DBpX8B\nh5eJWtR5N+xQjwjIn9UtE0eC8vWU+pW/9TB7svc9rwtf4ke7FS7XDdiel9hKvlUsaLe112Iz5pxs\nc/kgtNGdXHCLkqeTStW5NepPXiVUV8qa/UrS8WNdi3ub6w+xVWPmW93L9S2Ns9VGUjouzFr6toz+\nNxUtx0ilXwm2T+ewsrhFXTcxufmo1+o1tyqrdPXXyJfGlyLlNq8PFqfdN/c3lKqL0jbDKVnlNK0p\n2SfZkU6Lbeza/UZPE66ZquuPNN+Qxsyy6TUq3HRpZKrutY8Km188uhZqwKaVzOxI0yp5DgvCIZ4W\nTbrb8vNmWc1+WmK3KeJX1lKLPL8furnmc1S0uU9BDhG+tk9M898QUQpz1CD2uUtwX+XtHJ6cptsc\npsom2jscyIG0lpmAMaGjICWUZMIyBJDr00Ya0yfCr8SUtrsjS2Om16Fatr8tYPqSES6MlXYzyd3D\ndxkAFW4YAAyUsz6kXCnmfUi+PtzfI/qrAA1cAAAAAAAAAAANoJykkj0OJTGiiMddfM5HDq07Odrp\nE7Ve5dWUyrt+Px+N1tvkW33ZpF80tmt0+Z6FXRmO3bjisQWmWF2K67k67ENpGLHqDfscnZ1L3qqR\nyjXjeX8z+8Z2S4vXJrSW3zEJZ4bGf4+lxW3zdjS+nI9f4kJP8ZLmXirw+Ux4Uo1xxYWPnqlzv7Fm\nW3kyq8JOmMeeOu3MaRuhyrJsq+e5+G0vJGFXjKcpZE7VNOqcdQb9SGNbuoUrel8Htv1RNKqluOIp\nSXgfmcxHKb3PJqXPXY+Vr0KZ+l8fbSVsa7VlSl+W+mjbH4g8vJVGPVtvv7IxdVC66OK4OMNb2djh\nHD6cVOyK3JrWyvHjMrpOWWlyvFhGK5ltkWVw+F8HyycJeqLoOr68da0x7V47PxrcJuc480t8ql7F\nXw4cyo6pSXNtHseIY8cjFnFx20to8h4dsKnr97zdPscPLhMa6cMu0a81XN4+3qvp9zX9qqdijCGl\n6klle7I1cnyTXzP0N/wuHT31sznX8rXf4QZOZf1VSe/YrR/G29HzdTpPLxK1002K86FibgtaLz+V\n8RW/7VqMC7nUpy6fct/hl5sryzXKxwins3i7ZrszWbntS+VbIwfDy431vmfmmS0ydfNO7XX0LCon\nLuZeCpLU5dCMs4mY1WlxKuL+SG2aPiORZ0rrLixMOlfNrfuI5OLF8sOX9EYyzK+mnmOdP8db6o4P\nGq7K8pKyXNLlR6nKzfDklGPMmeb49J2ZybWvkR2cWMlY55bctLRkGToYNJI0JX2IyQAMkAjbRhIy\nBf4Zr8wjy48trXqYwLPDuafaRNnL5k/VFXRfPGo+ZJHsR92bx7FavwbbAwNlHZtkwNmNhG2Spl/U\ni1sqZXdFsfbn+R/VXABq4QAAAAAAAALuDeqPNZGPqwmea7PDsdLGTfmy5LUIPRtXFQphFeSIb5fN\no57XrYY6kiJLb2ySPcxGPQz2K10YRPHyLEexXj2RaphKUfli2QugynqpnMR1c+qcafmjor4WFZkS\n6aUV3l6F8ctPO+Rx9uT/AEr1UztlywW2zv8ABeG148PxeRZ2+VL0ZbxMGWDiuVcYyuck17otupSu\nePKv8tx5+n8xa5uXLHGeIkljzjTXj13bsjLnf2N4wnLItfSVXL8vpzEEbVyrLnFqb/K16L1EYare\nHGbUoPxNv0KKpOeax677KvzbHySfsSR5KJTx4rkilvb82zMbXK93qxOmUeWO/wCZCNMsiqCu6XKW\n9+xjyZzHG2pxVsiF8MRcr3cvPz0drhknTwmNlm3JRcpbIPBhC3xG3za0JXvkcP4WtNHJxfPwxyTn\nNuFwzinGviB5F+HbVRTW9RjJbbLXw78SzvtysbifLC3HfWS7PqV8PguTh3XfszN8GNr24yW1stYv\nwtPDxMjwbk8vJa57Zrevsethy48k3jWNjv4uZRm47tompR21s8lOyMrnkc2lFuJ3a6f2NwV1SlGV\ni31S+ps4UnCco1OPyWfNJ+5zfIs26OHwishONUoeJuc+sUQQ4bk2dZS19yy5wcne006uiT8ySrNn\nbFyUdIwx3fTS6ntDDhC7zmWo41GPHuirdLJs/dt/oaR4flWfW+/uT5xvmni+onduLVY5+ZmfFofw\nQ2aw4TrrbYSrCxK+supGWUqdaVJcRvm9Qhr9CzCV0quabaNpZmHR0jHf2RpLPhctKPQth5vpXJEs\nR5T2pEsOG01tNz6ledlqXLWuVE1Nc/C3KW2y1xyt9qyprJ41S+fTPL8ethdnOda1HlO21CxuPLvR\nWuwqJ/VDbNcf4e1cv5PN1xc5qKW2zoLhU+Te1s6ePw+iu9SitaLGTTJxXhPrs3x5JWNxseVuqlVN\nxl3RA1pnqsjhUJvmf1a6nHyeFXQm+SPNEv2hquYZJ5Yd8O9bI3TZF9YNfoO0/aNU17DRajBcq2Zd\nLfVRf9CU9arR2pJryJp2uxLfkaTXLJo1KWtsePLTXzNl2AI23wx6smBsxsrpptkGNmNkI7MlbK7o\nsFbJe5Ivix57vFAADRxgAAAAAAAMxW2WsOCeTD7laH1HRwcexWqco6RF9NeLXabdje0itL5rHs3i\n5a7DwpN7MOtehlzYxhNGspLfQkVPqzeNKSHSk+TJ6axs0l0JFm2wr5YS5UbKK1rRHZXuJPRW/ItQ\n+NOy2PNKUlvqtnqMfHqhVGEU08mP9DzeJVy3xc+2+57CErIyvtdcZqv9z9iNac+ee2i0ra7Yz+Wl\nckn6Mx+fXjWRT3Y5cy9dEjVUlXTKtx/E/PJryZmM8d2yyOaUVUuTXqQybeJvJjW4J1cnNpLzNar6\npVvKth87l4bNvAnGhVRtXi83P/5TdVSsvn0TpUeb/wAxCGs6KfCWJXPU4vxCbGXiZVl8ZbraSRXU\npKqOTKv8yT5JfYmxeWmdmPHsupyfNn+Gr4p7Hsi0c/I43j08SWFYnGX83kRrjuN+GvufMo0y5X7s\n8jH4/JJ4i1rqx3Fprui8+IVUY/PdNJ+mzz2LxfxcyGPdRKp2R5ob80XsnFpufNd5dls7PjZZ8HJ1\ny/KvTsqZ+es2zxYzapj017lFV2xol25m9x+xvZVGNsqFH5H1TNVOO43y6KHytHZle13W8moy4vxa\n65R+WS+YtqFNcHGOkUmpRhOMZblPrEgjVbrbnt+aJx8K5eViu5Y7bmaz4pY+lUSNw2/m2wtLql0G\nWrdk3G1NmRdZ88nokuaivnfRkUJuMuhm2cLtKS0/I0w1FMttHRWnzaLFVUFXvRA5baSLEpKurZbC\neUZI3OMpuK7iF3KtEXNCTc4dzTb8xl4pF1KKhtLv3KzsjPfL10FdpcrXQgdlNKluSiTbuElbqXmj\neM9TTfVHPs4lRBNR22VbOKza1CKXuVkq/XbtZdjnD8t9StLIqhBc80mvc4lmXdZ3sZA233ezTX7W\n6OxdxKlL5FzMo2Z9k+mopfYqAmSRaYxPT889s99j149WLVKVcF8ifVHg8SPNZH3aPV8UulCNdce0\nak2a4qc/8daeV4pZG3iN84RUYuXRIqbM2S5rJN+bNdkVpMvDOxsxsxshG2dg12Y2Sjs22Y2Y2a7G\nkXJtshyEt7N2yG17ZbFjyXcRgAu5wAAZjFyektkscW6XaDLnBIRnlT5lv5P80d1VLyiv6Aeer4Zf\nPutFmvhH88jtKv2NlV7Acyrh1MH22y3GC7IteEbKr2ArKs2VfsWlWbKv2IWioqn6G3hFxVG6q9iq\n8UVSbKku+F7G0a/YirKSo3pep3KqlDwqoXf936z36MrKrovl6llxacNw3K1csn7FajaRK+NdsoyU\npSe6fsbcspWxodf5c47s16kW4Rip7klivlfubLngrIRs3Zb88PZFEtoW1tLLkpRdb8LXqZVLhCVE\nLd2OSs/QkUuayqvkjKlx3P8A3jWF0dTyZ16kn4aXboQhNz2SyeZSUqeTpv1MYvNZXOd0fzN637Gj\nqrcY4kLOWSfibZJB82TK2NqdXLyvr2Zh8nC58dxntaVzM34fqzeITyLW3GUOVL0fqVsf4VUcC3Hs\nubcrOeMvQ7rtsrolOVfzb6RT309ST8QlZCt7Tkt9jzrn8vCa0nWLnYPCJ15UMnKvd1lceWHT6URc\nQl4l+ozel6FvMybJ1T8GXVPWjnyxrHOuO/qW3I24OPkuXfkT6aObjHUX+pjnTSg4rkfV/cw4T+eW\nnqHTZpLmUYya6S7HZ0ie1iWUocym49YdEjTlhyygpadnX7EW3/F0foY3/XyI6J7ppVt+GoyXyfWY\njFuVqaTgusfci3re+m/cjsyoVRSlalr3H11PeJpLVEbZV6betGHXGNqr31ktr2KlnGoQflP2II8e\nS3KVKlZ/C/QfXktJK6Ea+7hLrF/NsxmZFdVcVa9Rl2aOLPiltkp8qUFPukdbglUcum53Q5+RdObr\no1wxu0cmExx3VWWdjY8OWDlNPt0KlnF5NahBL7nW4zhwr4bZd4Si46S0jyvY0uKuGqs2Z19nebX2\nK7k29ttmNmNjTRtsGuxslO2djZqNkjbY2a7Mx6kWJxdPhNLvy6oL+ZHf4tbCEb5vyjyoqfC+K5Xu\n1x6RRj4rmqa41LpKT2zTGOfmu83l9jZrsxzEaW235ka7NdmOYaR2bbMNmnOauZbSlySbRq5aI3Nm\nOYnTO5pHPoRN7GzBKluwABAAAOv8OR5syzf8j/uj0iqOB8Kx5s+xf/Tf90erUAK0ajfwi1yGVW32\nRGxU8PRlVlxY85PSiyeGBY4/TobHOVfU2VTOxVwx66lqHD649wlw66JPyZZrwLJ9ondppqgvsSSy\nKa/Qhbbjw4RNrrotVcGgus2T2cSgl8iKlvEbX9PQJ8r0MLHr10X6kPEeWFP5XLzL6Sgsi2e22HNS\nW7X0j17lKmSo3zapg69xt/emqsh4k75wa8D5Yr2NJ5OPCrIlHISna9xXoU7uL0qcYxXPDl1Na1tm\nd8NOlvp0lWlR4Ndvz2vnj7Eii7Lkk4yrUdS+55az4hnLJUvDUeVNL7Ea41PlcFNpN7Ei/wBNerdk\n66fHsr/MT5P0IZOiG8bmUOf5+bZ5x8TtuXW9s5fEMicrU1Jsnr2V5MOk293DIrnbC+Ny5K1yyTfd\nmysurqulLlla3uv7HziORdHopySfuWauMZ1U4y8aTceiLdNMpXv5TWqanVpXfvGvU0bpsslBc0I4\nvV/7R4uv4izoRa5k9vfXuWqvia5NKdMNPpPS7or1WkyvqPVJVShK6U0qrlqtP1NLaZQx63OMZWRf\nSC8keYzOP+LSqceuMIRe477nOyOJ5WRbzzulvWuhHVrOOvWZl2JjXqdrTnKO2k+xyMjimLVTZGnm\nlY30k/I4Fljn3bbNNiYtJxz8rVmbfY/mmyCU5Pu2/uabBOk6kZ2xswAtElfWR674bcK+G32S7bPI\n1dz1OBBw4DJdnZIvh7U+Tf8AHIg+IuIq7CdUOkWzy+zrcc+SFUTjtk32z4/GLI2a79xsrpfbbZjZ\nrsbJR2bbMbNWzGwjs32WMapzkku7fQrRTbPUfDXDXbP8TatVw6r3JkW7am3oeGVV8L4eufrNx2zw\nvHs+efxKyyXZdEem+JOJxxcaUYy/MmtJeiPBue2231LOTt52k5jVyI3I1bJ0m5pHM1c9mgGlLlWd\nsbMAlXYAAAAAAAAAAPRfBdMruJ2Rj/qn/dHua+GykuvQ8j//AJ1r9rXt/wCpf90fQnbGJFFOHDYr\nuWI4lUf4TMsmKI5Zi8lsJWYwhFdIozzRRQnmT8uhDLIsl5gdKV0I/wASIZ5kEu5zm2+7MNb7kC1Z\nkuzSi9EUm9Nzl0XmRw7kPFpeHw6yYaYeS3Oxal1sTfsUreNVJfJFs89K1y6pGE5P2KbdExjqWcYu\nl2fKvYqW59k1p2S6+5V5fNmeaEQtJG7slLsjMW0upF4yXYxK35G0VaS6VW9ybMABJvXYa82A+w9I\nsl9ktaNNGWYG0dYdPQwADchsAxsG2TBjY2EbZMbMBslG2dmUzTZmPcJlWaFuaR7OrGl+zsevetJy\nZ5zgGG8vPhHXyrqz22VGGNjWXTa0o6ii+M8MfkZbsj5/xy5WZnIu0Fo5rZvkXeLfOx/xPZE2QiXU\nZ2NmmxsnStybbMbNdmvNsaR2b7G9voYrhOyWoptnd4TwKd8lO75Ie40mX9ouDcLtzr0uVqC7s9pk\n5GHwnhz5pJQgukf5mcjiHFMXg+N4GJpz82jx2dxHIz57um2l2RaKZ578M8T4hZxDLndPom/lW+yK\newCWQAAAAAAAAAAAAAAAAAAPSfBDceJ3NP8A8L/NHtnJt9WeI+Cf9JW/8L/NHtyAAMEgYZnZhgAN\nmHIhLOzn8dyFDh8ot65nouua8zznxBa8i5VRfSJDTCOb4sEuhHK9t6iYVUY/UzKlFfSijbYnbIyq\nfVhSk36GeWXn1CZa3ShFd9kN819MRbuK9CBvfcheMgxsyQvtkPsZXYxIJaMwG+pq2EWs7GzXZjZK\nm2dgxsxsIrOw2a7GxpDI2a7MbJ0jbYkqXzES7na4Fw55uVGLXyLrJk6XmUk3Xc4BjvDxPF1+Za9R\nXsafF2dKjDjS5fmWeXojvaqxqndPUa6Y9NnzXjXEZ8S4hbe38u9RXoi7jyy3dqmzVy9zTmZjY0js\n3czNcZ2y5YLbZEZjJxe4vRKNu1i/D19yTsthWn7onyOCV4dkJeJGyHm+ZHC/E3f62f8AUxK+2S07\nJP8AUG3qvxfC8Fc0dTl7I5+d8S3XxcKEq4nB2ANrLJ2ScpybbNQAgAAAAAAAAAAAAAAAAAAAAAej\n+Cv9I3f8L/NHtjxHwY0uIW7/ANW/7o9m59AN2zDZXtyqq0+ecV+pzcnj2JQvr5vsQOzzGkrUu7PK\nZPxUntUQ/VnLyOO5l+1z8q9gPcW5lcPqmkVbOJ1pdHs8G8m5vbm2WcXMnzqM3tMiteKYXL+T02Rx\nWfK+RHLnZZbJzk+rEJ8y0ZlHT9jLddnSWtFD1NJLXYlNZRfcbWnHI1gWq0n3K8ejJ4PqFusa5UOn\nY577nWsXNA5tkdSaClxRmUDKBIyuxiRk1t+WPUJqFvqYb6mNmGyZGW2Qa7MbLaV7NmY2a7MNoaUu\nTbYcjTmNebZOkdm7ZjezEYSk9RTZ1uG8ItvlzTXLH3JP+oeH4U8u6MIR3s91gY1HC8TkckptbnIo\nUxo4XRuOubXdnmOL8ZsyZuFc/k82Tpnlnvwv/FXxB+L/APY8Wf5Mfqa8zywBLMAAAAAAAAAAAAAA\nAAAAAAAAAAAAAAAAAAAAF/hHEf2dfO1LfNHlLWT8R5du+X5DjDYE9uZfc9zsk/1IW9vqYAAAADKe\nnswAO3w26Ntem/miXPqetdDzuPc6bVNfqeix768irmh38zLOPQ+Pl3mr7RNNSD7FnlU1rzIZ1uPQ\nzdXVB2JYsjktG8OwUsWYvaKeRD5tot1dTS+O0DTn8pnRL4bJK6GxsmCCEdsq5E+ael5F/MlGivlj\n9TOS35vuXkY8mUnhs2a8xq59DXmLyOW5t3JaNXJEbbMFtM7k3cxHc5KKXVkZlNp7XcaV27GLwKy5\nKVlsYJ+50J/D1Mcbdd0ZWLr1ejzaybl2sl/Uy8q9rTtlr7jRuvTY8OH4VClfOHiehBkfEcK4uGNW\neccpN7bbNRo2tZfEMjLk3ZN69NlUAlAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA\nAAAWMXKnjy3F9PNFcBMyuN3Hp6LIXR5oMmj8/Ro8xRkTpluMno6+HxGE+jepGGWGnq8PyMc/GXtf\ndHXsYVOn7EkLVJb2iWMYy7SRm6fr35iGMdGWtosrHkZ/DT80l9xs+vSoqFrZjInDFq5pPr5I3zcz\nHw62lNTs9EzzOXm25M25vp5I0xxc3NzY4TUbZOR4trkVnPZrswbaeblncmW9mACWYAAAAAAAAAAA\nAAAAAAAAAAAAAAAAA0AAM6AwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAZT12MACerLtr1yzZbhxa\nyPVpHNBW4ytsefPH1XX/AG5al0WiC/i2Rctc7SOeB1ib8jkvutpTlJ7k22agFmNtvsOpwDhseI5z\njbLloqg7LX6RSOWeh+EmrJcQxE9WX4sow936BCeq7gnEHZiQxvwr5ZeHdKW1tLzPPSxbI434jl3V\nz8nN7lvh3B8jOzfw7i60tucpLSikX51S/wCx3LD5+TOfVf7rA5K4blSvppVTdl0OeC9V/wBIsv4f\n4hHkboepvSZ6StSq+J+BKS1KOGtr35ZFTGzciXw3xyUrpbV8Nde229gcC3hWXXmxxJVPxpLaXqiZ\n/D/EFGt+A/zHpex6XFk5cY+HbJ9ZPF6v17nPx8/Kl8OcYsldNzd0Ou/VgcrH4ddicXpx8vGlNyf7\ntd5L2IqOFZWZObx6W4qfL9meopk7OJ/DU5NuTo6t933Of49lHwxnuqTi5ZvK2vTT/wCQHDu4bk05\nkcWdbV0vpj6m2ZwnLwa1ZfU4wb1s9JK+2HFPh/IhX49qxduPnL6v8iLitNWRwTKvwb7VTCxStptX\naT9GBw+BY9eVxnFouXNXOepI79L4Jl8UfDVhWVylY61Pm7P17nF+GP8A5iwv+IdjN+IMXDz7/wAN\nw+uORCySVm+z21sDiR4Pk35uTRjw5/Am4yfl3aIreFZdObDEnU/Gs+lep2eDxslwfLy8rJnXjWXr\nmVa3Kcu/+Z2LVGXGPhyUOdRcHpz+rXlsDgcK+HMi7idNGXBwhOEp/ok/8yrLF8Lh2bz40nZXbGPi\nb6Q79DrcEzcjK+MJu22Unu7Sb/2ZdCHE3/2a4r4if/equbf3YHLr4JnXYyvhQ3Bra90VJYtscaGQ\n4Pw5ycU/dHquKW8Sj8T0wwG+VVwVK/h1y9f8dlHJrsl8KYsYpyccyaaXXyQFfg3DFPingZlb06JW\nJb/2dplTC4Tl5/NKipyinpvyPVcrj8TVRktNcO0//wADnYNaq4DTZmZFscey2Xh11Lq2uj2BQ4Tw\nK3J43HByYOHKm5r20yrDhWTdm3Y9Vbcq29r0PXXScfjnB5eZc2Kuj7/RLuU+DumPAeJyyHapPJSs\ncPqUfL9O4HlszDvwrfCyIOEtb+51vhTExsrJypZdbshTRKxRT1tol+I7qreG8PVdduoqSjZZ3lHZ\nv8E2eFkcQsSUuTEm9PswK3EM3hduJOGPw+VVj1yzcuxYwuH0ZPGKKbcaVcHic7i33et7KHEeNPPx\n/BeNTV829wjpnpa3/wDEuF/9tX/6sDzy4NKXw5PiOtyVvKv93XcrUcFzr8ZX10twa2vdHRrlN/Bl\niUpaWZ116cp3sq7Fo4pg8qyHOFdfhRq+mS0B4zB4Vl57ksepy5HqXsytfVOi2dVi1OD00er/ABU6\neFcetx+apyyIJLzXVnkpzdk3KTbk+rbA1AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACXHvsx\n7o21ScZxe00RADs5PxNnZGPOl8kFNalKEdNlfh3Gcnh1c66uWVc3txmtrfqc4AdKfHMyziFObOW7\nqo8sX26df+ZWjnXRxr8dS/LukpTXq0VgBehxbKhbjWRn82NHkr9kRQzroYl+Mn+XfJSl90VgB0Ic\nYyYW4lifzYkeWv2RC8+54dmLv8uyzxGvcqgC8uK5MbcS2M2p4sVGt+3/AEyXP47lZ+P4FnJCtvcl\nBa5n7nMAE+HlWYeVXkVdJ1vaI7bHbdOyXecnJ/qaADo8O4zk8Opspq5ZVzabhNbW/U3t4/m25GNd\nOac8dt1v02csAWsPPuw8+OZS9Wxk5J/fuS38VyLq8mt6UcianNL1RQAHXq+Is6rEjjxlH5YuMbGv\nmS9mR8O45lcOplTVySg5c2pren6nMAHSfG8uWd+LlJO3w3Xv21o2wuO5eHifh4ckq1Lmipx3yv2O\nWAOjPjeZPiNGc5/n0xUYy9v+mb08dy6M27Jhy7v/AHkGtxf6HLAFziPEsjiV0bL2tRWoxj0UV7Ix\ng8QuwHc6Xrxa3XL7MqADOzoLjOUsuGSpfPCrwV/u60c4AXsbimRjYN+HBp03dZJrfUs4/wAR52Pj\nRpi4ScFqE5R3KK9mcgAWo59yxb8fmbjfJSnt92iqAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA\nAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA\nAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAH//Z\n",
"text/html": [
"\n",
" \n",
" "
],
"text/plain": [
" "
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.7.3"
}
},
"nbformat": 4,
"nbformat_minor": 2
}