{ "cells": [ { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": false, "deletable": true, "editable": true }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Contre-courant\n", "X|z=0 = 1.0 X|z=1 = 0.602\n", "Y|z=0 = 0.806 Y|z=1 = 0.01\n", "A(X0-X1) 7.96e-06 S(Y1-Y0) -7.96e-06\n", "Efficacité dans le raffinat 0.4 d'extraction 0.4\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEKCAYAAAD9xUlFAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl8VPXVx/HPyUZIgIQs7EtCcEMBWYUAIqIVbdVq9RG3\naquiIi59Wqt16WYXW9un7vte17prK26IG2FHBATRhDWIQMJOCGT5PX/cyxgpy5DMzJ0k3/frNS8y\nd+7cOYxXTu79nd/5mXMOERERgISgAxARkfihpCAiIiFKCiIiEqKkICIiIUoKIiISoqQgIiIhSgoi\nIhKipCAiIiFKCiIiEpIUdAAHKicnx+Xl5QUdhohIozJ79uwy51zu/vZrdEkhLy+PWbNmBR2GiEij\nYmbLw9lPt49ERCRESUFEREKUFEREJKTRjSmIiAShqqqK0tJSKisrgw5ln1JTU+nSpQvJycn1er+S\ngohIGEpLS2ndujV5eXmYWdDh7JFzjvLyckpLS8nPz6/XMaJ2+8jMHjWztWa2YC+vm5ndaWbFZjbP\nzPpHKxYRkYaqrKwkOzs7bhMCgJmRnZ3doKuZaI4pPA6M2cfrJwIH+Y9xwH1RjEVEpMHiOSHs0tAY\no5YUnHMfAev3scupwJPOMw3INLOO0Ypnefk2/vb2YopKyqisqonWx4iINGpBVh91BlbWeV7qb/sv\nZjbOzGaZ2ax169bV68M+K93EvR8Uc85D0+n7u3c47+Hp3PtBMZ+t3EhNrdapFpH4tnLlSvLz81m/\n3vtde8OGDeTn57Ns2bKIfk6jGGh2zj0IPAgwcODAev0LfkrfThxzSC4zlqxnSkkZRcXl/PWtxcBi\n2qQmMaRHNoUF2QzrmUPPdq0axWWiiDQfXbt25fLLL+f666/nwQcf5Prrr2fcuHFEuu1PkElhFdC1\nzvMu/raoaZOazHG92nNcr/YArNuyg6lLyikqLmNKSRnvLFwDQLvWLSgsyKawZw7DeubQObNlNMMS\nEQnLz372MwYMGMDtt9/OJ598wt133x3xzwgyKbwOTDCz54CjgE3OudWxDCC3dQtO6duJU/p2AmDl\n+gqmFJcxpaScT4rLeHXu1wDkZacxtCCHYT2zKSzIISs9JZZhikic+d0bn7Pw680RPWavTm34zcmH\n73Of5ORkbrvtNsaMGcM777xT77kI+xK1pGBmzwLHADlmVgr8BkgGcM7dD7wJnAQUAxXAT6IVS7i6\nZqUxdnA3xg7uhnOOL9dsZUpxGUUlZbzx2dc8O2MFAId1bMMw/1bToPwsWrVoFHfhRKQJmDhxIh07\ndmTBggUcf/zxET++Ode4BlkHDhzoguiSWl1Ty7xVmygqLqOopJxZyzews7qWpASjb9dMhvm3m/p1\ny6RFUmLM4xOR6Fq0aBGHHXZYoDHMnTuXc889l4kTJzJ8+HCmT59Ox47/XbS5p1jNbLZzbuD+PkNJ\noZ4qq2qYvXxD6HbT/NKN1DpITU5gUF4Whf7tpsM7ZZCYoEFrkcYu6KTgnKOwsJDf//73HH/88dx1\n111MmzaNp59++r/2bUhS0H2PekpNTmSYPxANsGl7FdOXlFNUUk5RSRl/eesLADJaJjOkRxbDeuZQ\nWJBDQW66KptE5IA99NBDdOvWLXTLaPz48Tz22GN8+OGHjBw5MmKfoyuFKFm7uZKpS8q9K4niclZt\n3A5A+zYtKCzICZW/dlJlk0ijEPSVwoHQlUIcatcmlVOP7MypR3bGOcfK9duZUlLGlOIyPvxyHa98\n6lXf5uekhxLEkB7ZqmwSkUApKcSAmdEtO41u2d04e3A3amsdi9ds8W41FZfx2tyveXr6CszgsA5t\nvNLXnjkMzssiXZVNIhJD+hcnAAkJxmEd23BYxzZcNDyfqppa5pVuCk2ie6JoOQ99vJSkBKNft0xv\njkRBNv26tSUlSesiiUj0KCnEgeTEBAZ0b8uA7m25cvRBbN9Zw6zl6ykq8cYk7nr/K+6c9BUtkxMZ\nlJ8VmiPRq2MbElTZJCIRpKQQh1qmJDLioFxGHJQLwKaKKqYt3dWOo5w/T/QqmzLTkhmSn82wntkM\nVWWTiESAkkIjkJGWzAmHd+CEwzsAsGZzJUUlXlXT1JJy3vr8GwA6tEmt07Mpm44ZqmwSkQOjpNAI\ntW+Tymn9unBavy4451heXhHq/Dp58Vpe9iubeuSkU+j3axraI5u2qmwSabScc4wYMYIbb7yRE088\nEYAXXniBRx55hLfeeitin6Ok0MiZGXk56eTlpHPuUd2prXUs+mYzU/3xiFfmrOKpaV5lU6+ObfxJ\ndNkMzs8iLUX/+UUaCzPj/vvv58wzz2TUqFFUV1dzww03RDQhgCavNXlVNbV8tnIjU4rLmVJSxqcr\nNlBV40hONI7smum348jhyK6ZqmwS2Yd4mbz2y1/+kvT0dLZt20br1q25+eab/2sf9T6SsG3fWcPM\nZd8uNLTg6004B2kpiQzKywq1B1dlk8h3fecf2onXwzfzI/sBHXrDibfud7dt27bRv39/UlJSmDVr\nFi1atNh3rD7NaJY9apmSyNEH53L0wV5l08aKnUzzezZNKS7jT296y51mpiUztIc/aF2QTX6OKptE\n4kF6ejpnnXUWrVq12mNCaCglhWYuMy2FMUd0ZMwRXvvdupVNRcVlTFzgVTZ1zEhlaEE2w/zbTR0y\nUoMMWyRYYfxGH00JCQkkJETndq+SgnzH7pVNy8orKNpV2fTFWl6e41c25ab7CSKbIT2yyUxTZZNI\nU6CkIHtlZuTnpJO/W2VTkT9o/dKcUv45bTlmcHinNgwryKGwZw6D8tqqskmkkdJAs9TbzupaPivd\n6C9ZWv6dyqZ+3dr6SSKbI7tmkpyoyiZp3OKl+igcGmiWQKQkeavMDcrL4prjoGJnNTOXbQg19rt9\n0pf84z2vsmlwflYoSRzWQZVNIvFKSUEiJi0liZEH5zJyt8qmXXMkPli8CIC2ackMLcgOzZHIy05T\nZZNInFBSkKjZvbLpm02V/prW3sD1m/O9yqZOGamhfk2FBTm0b6PKJolPzrm4/wWmoUMCSgoSMx0y\nUvnRgC78aIBX2bS0bBtT/IWG3lu0hhdnlwJQkJseWtN6aI9sMtKSA45cBFJTUykvLyc7OztuE4Nz\njvLyclJT6/+LlQaaJS7U1joWrt4cmiMxY+l6tlfVkGBwROeM0ByJQXlZtExJDDpcaYaqqqooLS2l\nsrIy6FD2KTU1lS5dupCc/N1fptTmQhq1ndW1zF25q7KpjE9XbKS61pGSmEC/bpkM82839emiyiaR\ncCgpSJOybUc1M5d9uxrdwtWbcQ7Sd1U2+bebDu3QWpVNInugklRpUtJbJHHMIe045pB2AGzYtpOp\nS8pDcyQm+5VN2ekpDAm148imW5Yqm0QOhJKCNEpt01M4qXdHTurtVTZ9vXE7RSXloTkS/5m3GoDO\nmS1DVU2FBdm0U2WTyD7p9pE0Oc45StZtY+quJUuXlLNpexUAB7VrFVpo6Kge2WS0VGWTNA8aUxDx\n1dQ6Fn69mSklZUwpLmPmsvVUVtWSYNC7SyaF/u2mgXltSU1WZZM0TUoKInuxo7qGuSs2huZIzF3p\nVzYlJTCgW1vvdlPPHPp0ziBJlU3SRCgpiIRp645qZi5dHxq0Xrh6MwCtWiRxVH5WaLb1Ie1ba9Ba\nGq24qD4yszHAHUAi8LBz7tbdXs8BngI6+rH8zTn3WNQCcg70P7XsplWLJEYd2o5Rh3qVTeu37WRq\nSbnfjqOMSV+sBSCnVQpDC7yV6Ib1zKFrVlqQYYtERdSuFMwsEfgSOB4oBWYCZzvnFtbZ57dAS+fc\ndWaWCywGOjjndu7tuPW+Ulg9D14dD8feCAePUXKQsK3auJ0i/ypiSnEZa7fsAKBrVksKe3idXwsL\ncshtHfmlEUUiJR6uFAYDxc65JX5AzwGnAgvr7PMN0Me8a/JWwHqgOirR7NgMVdvg2bHQeSCM/jX0\nGBmVj5KmpXNmS84c2JUzB3b1K5u2ep1fi8t4c8Fqnp+1EoCD27cKdX49qkcWbVJV2SSNTzSvFM4A\nxjjnLvafnw8c5ZybUGefRGAScAjQGjjLOfefPRxrHDAOoFu3bgOWL19ev6BqqmDuM/DhX2DzKsgf\n6SWHLvtNniJ7VFPrWLBqkzdHouS7lU19umQyrKdX2dS/uyqbJFiBDzSHmRRuAtoBVwMFwLtAX+fc\n5r0dNyIDzVWVMOtR+PjvUFEGh5wEx94E7Q9v2HGl2dtRXcOc5Rv9xn5lfFa6iZpaR4ukBAbmtQ1N\nouutyiaJsXhICkOB3zrnTvCf/wrAOffnOvtMBP7knPvYf/4+cL1zbsbejhvR6qMdW2H6fTDlLu/2\nUu8z4JhfQXZBZI4vzd6WyipmLP22Z9MX32wBoHWLJI7qke3NkeiZw8HtW6mySaIqHpJCEt5A82hg\nFd5A8znOuc/r7PMPYJNz7rdm1h6Yg3elULa340alJLViPRTdBdPvh+od0O88GPlLyOgS2c+RZq9s\n647QanRFJWUsL68AIKdVCwoLvk0SqmySSAs8KfhBnATcjleS+qhz7o9mdhmAc+5+v+LoMaAbkADc\n6px7al/HjOo8hS1rvFtKsx8DDAZdBMP/F1rlRufzpNkr3VBBkb9caVFJOevqVDZ5a1p7t5tyWqmy\nSRomLpJCNMRk8trGFd5g9NxnIKklDB0PQydAy8zofq40a845itdu9ZcsLWfaknK2VHrFeId2aB0a\njziqRxatVdkkB0hJIRLKvoLJf4LPX4bUTBh2NRx1KaSkx+bzpVmrrqllwdfeanRFxeXMXLaeHdW1\nJCYYfbpk+FcS2fTvpsom2T8lhUhaPQ/e/wN89Takt4MRP4eBP4EkXdJL7FRW1TBnxYbQ7aZ5e6hs\nGtYzh96dM0jUQkOyGyWFaFgxHSb9HpZ/AhldYeR10PdsSNSyFBJ7uyqbdg1ahyqbUpMY0iM71I6j\nZztVNomSQvQ4B0smw6Rb4Os5kN0TRt0AvU6DBNWdS3DWbdnB1CXloXUkVqz3KptyW7cItQcv7JlN\nl7aqbGqOlBSizTn44j/ebaV1i6B9bxh9Mxz0PfVVkriwcn2FP4munKKScsq2epVN3bPTQoPWhQXZ\nZKuyqVlQUoiV2hqY/yJ88CfYsAy6DPZaZ+SPCDoykRDnHF+u2RpqDz59STlbdnxb2TTMbw8+OD+b\nVi10O7QpUlKItZoq+PSf8OFtsOVr6DHKu3LoPCDoyET+S3VNLfP9nk1TisuYtXwDO6trSUow+nbN\n9K8icujfPZMWSapsagqUFIJStR1mPgKf/B9UlMOhP4BRN0L7XkFHJrJXlVU1zFm+wV+ytJx5pRup\ndZCanMCgvCy/simbwzupsqmxUlII2o4tMO0+r33Gji3Q+0wY9SvI6hF0ZCL7tbmyiulL1ofmSCxe\n41U2tdlV2eTfbirIVWVTY6GkEC8q1sOU22H6g1BbBf3O9/oqtekUdGQiYVu3ZUcoQUwpKaN0w3YA\n2rdpQWFBDkP98tfOmS0DjlT2Rkkh3mz5Bj76G8x+HCwBBl8Cw38G6TlBRyZywFaUV4T6NRUVl1G+\nzVssMS87zVvT2k8UWekpAUcquygpxKsNy+DDv8Jnz0JyGgy9wnukZgQdmUi91NY6Fq/ZEkoQ05eu\nZ6tf2XRYxzahSXSD87NIV2VTYJQU4t26xTD5j7DwNWjZFoZdA4PHQYomFknjVl1Ty7xVmygq9gat\nZy/fwM4ar7LpyK6Zoc6v/bqpsimWlBQai6/nehPgit+FVu3h6Guh/wWQpMtuaRoqq2qYtWxD6HbT\n/N0qm4b5t5t6dWqjyqYoUlJobJZP9foqrSiCzG4w8nroc5b6KkmTs2l7FdOXlIfmSHy1disAGS2T\nGdLDSxKFBTkU5KarsimClBQaI+egZJKXHFZ/BjkHe3McDjtFfZWkyVq7uZKpS7wEMaW4nFUbv61s\n2rXQ0LCe2XTMUGVTQygpNGbOwaLX4f0/Qtli6NDHa53R8zj1VZImzTnHivUVTPFLX6eVlIcqm/Jz\n0kPLlQ7tkU1bVTYdECWFpqC2Bub9y+urtHEFdBsKx94MecOCjkwkJnZVNtXt2bRtZw1m0KtjGy9B\nFGQzOE+VTfujpNCUVO+ET5/0+ipt/QYKRnt9lTr1CzoykZiqqqllXukm/1ZTGZ+u2MjOmlqSE/3K\nJn+hoSO7ZpKSpFuudSkpNEU7K2DmQ/DJP2D7BjjsZBh1E7Q7NOjIRAKxfWcNM5et9+ZIlJQxf9Um\nnIOWyYkMys8KzZHo1bENCc28sklJoSmr3ART74Wpd0NVhVeldMz10DYv6MhEArWpooqpS8r9dSTK\nKFm3DYDMtGSG9sgOzZHokdP8KpsinhTMrDPQHQjduHPOfVTvCOtJSaGObeVeN9aZD0NttTe/4ehr\noU3HoCMTiQtrNld+u9BQcRlfb6oEoGNGqtevyb/d1CEjNeBIoy+iScHM/gKcBSwEavzNzjl3SoOi\nrAclhT3Y/DV8dBvMeRISkryZ0cN/BmlZQUcmEjeccywrr2BKcRlT/dtNGyqqAOiRm+6VvxZkM7Qg\nm8y0plfZFOmksBjo45zbEYngGkJJYR/WL4UPboV5z0NKKyicAEPGQ2qboCMTiTu1tY5F32wOdX6d\nsXQ9FX5l0+Gd2oTmSAzKa0taSuOvbIp0UpgInOmc2xqJ4BpCSSEMaxd5rTO++De0zPKuGgZfAsma\n/COyNzura5lXujE0R+LTFRuoqnEkJxr9uralsGd2qLIpObHxVTZFOim8BPQFJgGhqwXn3FUNCbI+\nlBQOwKrZMOkWWDIZWnf0xhv6na++SiJhqNhZzcxlGyjy50gs+NqrbEpLSWRwflaoPXhjqWyKdFK4\nYE/bnXNP1CO2BlFSqIelH8P7t8DK6ZDZHUbd4K0El6AOlSLh2lixk2lLykNXEkv8yqa2ackM9de0\nHtYzh7zstLisbIpG9VEKcLD/dLFzrqoB8dWbkkI9OQdfveMlh2/mQ+6hcOxN3hrScXgCi8S71Zu2\nM7XESxJFJWWs9iubOmWkhvo1FRbk0L5NfFQ2RfpK4RjgCWAZYEBX4AKVpDZCtbWw6DWvr1L5V96s\n6GNvhoJjlRxE6sk5x9KybUzxFxqauqScjX5lU892rSj0rySG9sgmIy05kBgjnRRmA+c45xb7zw8G\nnnXODWhwpAdISSFCaqph3nPwwV9g0wroPtxrndFtSNCRiTR6tbWOhas3h+ZIzFi6nu1VNSQYHNE5\nw7/VlM3A7lm0TInNbdxIJ4V5zrk++9u2h/eNAe4AEoGHnXO37mGfY4DbgWSgzDk3cl/HVFKIsOod\nMPsJb57DtrVw0Pe820od+wYdmUiTsbO6lrkrN4bmSHy60qtsSklMoF+3TH8NiWz6RrGyKdJJ4VGg\nFnjK33QukOic++k+3pMIfAkcD5QCM4GznXML6+yTCRQBY5xzK8ysnXNu7b5iUVKIkp3bYPoDMOUO\nqNwIh5/mreWQc1DQkYk0ORU7q5mxdH1ooaGFqzfjHKTvqmzyFxo6tEPriFU2RToptACuAIb7mz4G\n7t3XZDYzGwr81jl3gv/8VwDOuT/X2Wc80Mk5d9N+g/ApKUTZ9o1eT6Wp90L1duh7DhxznbcanIhE\nxYZtfmVTSRlFxeUsKfMqm7LSU/yeTV5Lju4NqGwKvCGemZ2BdwVwsf/8fOAo59yEOvvsum10ONAa\nuMM59+S+jqukECNb1/l9lR4BVwsDfwIjfgGt2wcdmUiTt3rT9lC/piklZazZ7P3+fcmIfG78fq96\nHTPcpLDPudtm9i/n3P+Y2Xzgv7LH/sYUwpAEDABGAy2BqWY2zTn35W5xjAPGAXTrpt9YY6JVLoz5\nMwy9Aj78q5cc5vwThlwGhVepr5JIFHXMaMkZA7pwxoAuOOcoWbeNqSVlHNIh+i1r9tfQ42r/zx/U\n49ir8EpXd+nib6urFCh3zm0DtpnZR3gzp7+TFJxzDwIPgnelUI9YpL4yusApd8Kwq2Hyn+CT270E\nUXiVlyBatA46QpEmzczo2a4VPdu1isnn7XOY2zm32v9xvHNued0HMH4/x54JHGRm+f7Et7HA67vt\n8xow3MySzCwNOApYdOB/DYm67AI44xG47BPIGw6T/wB3HOmNPVRVBh2diERIuLVPx+9h24n7eoNz\nrhqYALyN9w/9v5xzn5vZZWZ2mb/PIuAtYB4wA69sdUG4wUsAOhwBZz8LF7/v/fz2r+Cu/jD7cagJ\nZJK7iETQPgeazexyvCuCHkBJnZdaA1Occ+dFN7z/poHmOLPkQ691RulMaJvvlbEe8SNIaHxdJEWa\nsohUH5lZBtAW+DNwfZ2Xtjjn1jc4ynpQUohDzsGXb3vJYc0CaHe4NwHukBPVOkMkToSbFPY3prDJ\nObfMOXe2P46wHa8KqZWZqQxIPGZwyBi49GP40SPe/IbnzoaHR8OSD4KOTkQOQFjX+GZ2spl9BSwF\nPsRrjDcxinFJY5SQAL3PgCtmwMl3wpZv4MlT4YmTYeXMoKMTkTCEe+P3D8AQ4EvnXD7evIJpUYtK\nGrfEZBhwAVw5B8bcCmsWwiPHwTNj4RvVEYjEs3CTQpVzrhxIMLME59xkYL/3pqSZS06FIZfD1Z95\n7blXFMH9w+HFn0J5yf7fLyIxF25S2GhmrYCPgKfN7A5gW/TCkialRSs4+hdechjxv7B4Itw9CF6/\nEjaVBh2diNQRbkO8dLxB5gS8DqkZwNP+1UNMqfqoCdi6Fj7+O8x61Hs+8CIY8XOvtYaIREXEGuL5\nLbDfc86NilRwDaGk0IRsXAkf3gpzn4Gklt/2VWqZGXRkIk1OREpSAZxzNUCtP2dBJHIyu8Kp93jV\nSgef4F093NHH+3On7k6KBCHc20evAf2Ad6kzluCcuyp6oe2ZrhSasNXzYPIf4cu3IL2dNw4x4EJI\nahF0ZCKNXqQX2blgD5vd/tY+iAYlhWZg5QyY9HtY9jFkdIWR10HfsyFxf019RWRvInb7yJfpnHui\n7gOv/YVI5HUdDBe8Aee/Cum58PoEuPcoWPAS1NYGHZ1IkxZuUtjTlcKFEYxD5LvMoGAUXPI+nPU0\nJCR78xseONrrsxSlFQNFmrv9rbx2NnAOkG9mdddCaA0E0hBPmhkzOOwHXnO9+S/CB3+CZ/4HugyG\n0b+G/BFBRyjSpOzvJm0RsBrIAf5eZ/sWvDUQRGIjIRH6ngVHnA6fPuUtEfrED6DHMXDsr6HLgKAj\nFGkSwhpojicaaBYAqrZ7k98+/jtUlMMh34djb4T2hwcdmUhciuhAs5mdbmZfmdkmM9tsZlvMbHPD\nwxSpp+SWMPQKr3XGqJu8SqX7hsFLl6ivkkgDhDvQ/FfgFOdchnOujXOutXOuTTQDEwlLi9Yw8lov\nOQy/Bha94fVVeuNq2LQq6OhEGp1wk8Iafz1lkfiUlgXH/RaunguDLoJPn4Y7+8HbN8K2sqCjE2k0\nwp28dgfQAXgV2LFru3Pu5eiFtmcaU5CwbFjuDUZ/9gwkp8GQ8VA4AVLVrUWap0hPXmsDVADfA072\nHz+of3giUda2O/zwHhg/HXoeBx/9FW7vA5/cDjsrgo5OJG6p+kiah9WfwaRboPhdaNUejr4W+l8A\nSSlBRyYSE5GuPjrYzCaZ2QL/eR8zu6mhQYrETMe+cN6L8JO3ILsnvPkLuGuAN/ZQUx10dCJxI9zb\nRw8BvwKqAJxz84Cx0QpKJGq6D4UL/wPnvewNTr82Hu4bCp+/qr5KIoSfFNKcczN226Zfr6RxMoOe\no2HcB/A//wRLgBcugAdHwlfvqq+SNGvhJoUyMysAHICZnYHX/kKk8TKDXqfA5UVw2gNQuQmePgMe\nOxGWTQk6OpFAhJsUrgAeAA41s1XANcDlUYtKJJYSEqHvWJgwC77/d1i/FB4/Cf55Onz9adDRicTU\nAVUfmVk6kOCc2xK9kPZN1UcSdVXbYcZD8Mk/YPt6OOwUGHUjtDs06MhE6i3S1Ud/MrNM59w259wW\nM2trZn9oeJgicSi5JQy7ymudMfJ6KJnsDUa/chlsWBZ0dCJRFe7toxOdcxt3PXHObQBOik5IInEi\ntQ2M+pWXHIZeAZ+/AncNhH//L2zWkJo0TeEmhUQzC62ebmYtAa2mLs1DejZ87w9w1Vzo/2OY8wTc\neSS8cxNsKw86OpGICjcpPA1MMrOLzOwi4F3gif29yczGmNliMys2s+v3sd8gM6v2q5pE4lObjvCD\n//MGpHv9EIruhjv6wge3QqU6yUvTEPZAs5mdCIz2n77rnHt7P/snAl8CxwOlwEzgbOfcwj3s9y5Q\nCTzqnHtxX8fVQLPEjbWLYPIfvXbdLbNg+M9g8CXemIRInIl0QzyccxOdc7/wH/tMCL7BQLFzbolz\nbifwHHDqHva7EngJWBtuLCJxod1hcNZTcMlk6NQP3r0Z7jgSZj4M1TuDjk6kXqK58lpnYGWd56X+\ntrrH7QycBtx3IEGLxJXO/eH8l732GW3z4D8/h7sHwmfPQW1N0NGJHJCgV167HbjOObfPpjNmNs7M\nZpnZrHXr1kXgY0WiIG84/PQtOOcFb92GVy6F+wph4etqnSGNRjRXXlsFdK3zvIu/ra6BwHNmtgw4\nA7jXzH64+4Gccw865wY65wbm5uYeYBgiMWQGB38Pxn0IZz7uXSn863x4aBQUT1JykLiXFOZ+s8zs\neQ5s5bWZwEFmlo+XDMYC59TdwTmXv+tnM3sc+Ldz7tUwYxKJXwkJcPhpcOjJMO95r0LpqdOh+zA4\n9mavW6tIHAo3KdRdeW0XB+w1KTjnqs1sAvA2kIhXWfS5mV3mv35//UIWaUQSk6DfudD7DJjzJHx0\nGzw2BnoeD6Nv9tZ5EIkjWnlNJJZ2VsCMB7xlQSs3evMdRt0IuQcHHZk0cZHufdTFzF4xs7X+4yUz\n69LwMEWamZQ0bz7DNfPg6F9C8Xtw71Hw6njYsDzo6ETCHmh+DHgd6OQ/3vC3iUh9pGbAsTd6fZWG\njIf5L3rLg755LWxZE3R00oyFmxRynXOPOeeq/cfjgMqARBoqPQdO+CNc9ak39jDzEa91xru/gYr1\nQUcnzVCwfy5yAAAPSUlEQVS4SaHczM4zs0T/cR6gTmAikZLRGU6+AybMhMNOhil3eLOjP7wNdgS2\nfIk0Q+EmhZ8C/wN8g7cM5xnAhVGKSaT5yi6AHz0El0/xJsNN/oOXHKbeA1WVQUcnzUC4SeH3wAXO\nuVznXDu8JPG76IUl0sy1PxzOfgYungQdjoC3b4C7+sPsx6GmKujopAkLNyn08RfWAcA5tx7oF52Q\nRCSky0D48WtwwRvQpjO8cTXcMxjmvQC1++wOI1Iv4SaFBDNru+uJmWUR/sQ3EWmo/KPhonfg7Och\nOQ1evhjuHw5f/EetMySiwk0KfwemmtktZnYLUITXJE9EYsUMDhkDl34MP3oEqivhuXPg4dGw5IOg\no5MmIqyk4Jx7EjgdWOM/TnfO/TOagYnIXiQkeG0zrpgBp9zlzWt48lR4/AewcmbQ0UkjpzYXIo1d\n9Q5vAPqj22DbOjh4DBx7E3ToHXRkEkcivvKaiMSppBZw1KXe7OjRv4YVU73xhhd/CmXFQUcnjYyS\ngkhTkZIOI34OV8+DEb+AxW95lUqvTYCNK/f/fhGUFESanpaZXlvuqz/zriDmPe/NcZh4HWzVUuiy\nb0oKIk1Vq1wY82e4cg70OQtmPOT1VZr0e9i+Yf/vl2ZJSUGkqcvsCqfe7VUrHXIifPx3Lzl89DfY\nsTXo6CTOKCmINBc5PeGMR+GyT6DbUHj/FrjzSJh2v1fBJIKSgkjz06E3nPM8XPQu5B4Kb10Hd/b3\nlgutqQ46OgmYkoJIc9V1MFz4b6+3Uuv28PqV3ipwC15SX6VmTElBpLnrcYzXjXXss5CY4s1veOBo\nr6S1kU1ulYZTUhARr6/SoSfBZVPg9Iehahs8exY8cjws/Sjo6CSGlBRE5FsJCdDnTK9S6eQ7YNMq\neOJkr7dS6eygo5MYUFIQkf+WmAwDLvTWjj7hT/DNfHj4WHj2HFjzedDRSRQpKYjI3iWnwtArvNnR\no26CZR/DfcPgpUugvCTo6CQKlBREZP9atIaR13rJYfg1sOgNr6/SG1d7t5ikyVBSEJHwpWXBcb/1\nksPAn8KnT8Od/eCtG2BbWdDRSQQoKYjIgWvdHk66Da6cDb3PhOn3ea0z3v8jVG4KOjppACUFEam/\ntt3hh/fA+Olw0PHw0V/h9j7wyT9gZ0XQ0Uk9KCmISMPlHgxnPg6XfuTNlH7vt15fpekPqq9SI6Ok\nICKR07EvnPsC/PRtyO4JE6+FuwZ6Yw/qq9QoKCmISOR1GwIX/gfOexnSs+G18XDvEPj8FfVVinNR\nTQpmNsbMFptZsZldv4fXzzWzeWY238yKzKxvNOMRkRgyg56j4ZLJcNZTkJAIL1wID46Er95VX6U4\nFbWkYGaJwD3AiUAv4Gwz67XbbkuBkc653sAtwIPRikdEAmIGh50MlxfBaQ941UlPnwGPjoFlU4KO\nTnYTzSuFwUCxc26Jc24n8Bxwat0dnHNFzrld6wJOA7pEMR4RCVJCIvQdCxNmwff/DhuWweMnwT9P\nh1Vzgo5OfNFMCp2BlXWel/rb9uYiYGIU4xGReJCUAoMuhqvnwvG3wNefwkOj4PnzYO0XQUfX7MXF\nQLOZjcJLCtft5fVxZjbLzGatW7cutsGJSHQkt4RhV3mzo4/5FZR84A1Gv3wprF8adHTNVjSTwiqg\na53nXfxt32FmfYCHgVOdc+V7OpBz7kHn3EDn3MDc3NyoBCsiAUltA8dcD9fMg8IJsPBVuHsg/Pt/\nYfPqoKNrdqKZFGYCB5lZvpmlAGOB1+vuYGbdgJeB851zX0YxFhGJd2lZ8L0/wFVzof8FMOcJbwLc\nOzfBtj3+vihRELWk4JyrBiYAbwOLgH855z43s8vM7DJ/t18D2cC9ZjbXzGZFKx4RaSTadIQf/J83\nIH34aVB0t9dX6YNboXJz0NE1eeYaWa3wwIED3axZyh0izcbaL2DyH7x23S2zYPjPYPAl3piEhM3M\nZjvnBu5vv7gYaBYR2at2h3qT3y6ZDJ36wbs3e+26Zz4M1TuDjq7JUVIQkcahc384/2W48E1omwf/\n+bk3IP3Zc1BbE3R0TYaSgog0LnnD4CcT4dwXITUDXrkU7iuEha+rdUYEKCmISONj5q3fMO5Dr2W3\nq4V/ne9NgiuepOTQAEoKItJ4JSR4FUqXT4Uf3gcV5fDU6fD492H51KCja5SUFESk8UtMgiPPgQmz\n4aS/QXkxPDYGnj4TVn8WdHSNipKCiDQdSSleuepVc+G430HpTHjgaPjXBbBO82PDoaQgIk1PShoM\nv8brqzTyOih+D+49Cl69AjauCDq6uKakICJNV2oGjLrBSw5DxsP8F+DO/vDmtbBlTdDRxSUlBRFp\n+tJz4IQ/wlWfQr/zYNajXuuMd38DFeuDji6uKCmISPOR0RlOvh0mzIRep8CUO7zk8OFtsGNL0NHF\nBSUFEWl+snrA6Q96S4TmH+31VrqjL0y9B6oqg44uUEoKItJ8te8FY5+Gi9+HDr3h7Ru8vkqzHoOa\nqqCjC4SSgohIlwHw49fggjcgowv8+xq4exDMewFqa4OOLqaUFEREdsk/Gi56B85+HlJawcsXw/3D\nYNG/m03rDCUFEZG6zOCQMXDpR3DGY1CzE54/Fx4eDSWTm3xyUFIQEdmThAQ44nQYPx1Oudub1/DP\nH8ITJ8PKGUFHFzVKCiIi+5KYBP3Ph6vmwJi/wLov4JHj4Zmz4Jv5QUcXcUoKIiLhSGoBQy7zZkeP\n/jWsmAr3D4cXfgJlXwUdXcQoKYiIHIiUdBjxc7h6Hoz4BXz5NtwzGF5rGn2VlBREROqjZSaMvtm7\nchh8Kcz7F9w1AN78ZaPuq6SkICLSEK1y4cRb4co50HcszHwY7jwS3vsdbN8QdHQHTElBRCQSMrvC\nKXfBFTPgkJPgk/+D2/vCR7fBjq1BRxc2JQURkUjK6QlnPAKXfQLdC+F9v6/StPsaRV8lJQURkWjo\n0BvOeQ4ues/rsfTW9XBXf5j9RFz3VVJSEBGJpq6DvJ5KP34NWneAN67yqpXmvxiXfZWUFEREYqHH\nMXDxJBj7LCSlwksXwQMjYPHEuGqdoaQgIhIrZnDoSXDZFDj9YaiqgGfHejOkl3wYdHSAkoKISOwl\nJECfM71KpZPvgM1fw5OnwBOnQOmsYEML9NNFRJqzxGQYcKE3x+GEP8Oaz71urM+eDd8sCCQkJQUR\nkaAlp8LQ8d7s6FE3wbIpXl+lFy+C8pKYhhLVpGBmY8xssZkVm9n1e3jdzOxO//V5ZtY/mvGIiMS1\nFq1g5LVw9VwYfg0sftNbAe71q2BTaUxCiFpSMLNE4B7gRKAXcLaZ9dpttxOBg/zHOOC+aMUjItJo\npGXBcb+Fq+bCoIth7jNwZ38oujvqHx3NK4XBQLFzbolzbifwHHDqbvucCjzpPNOATDPrGMWYREQa\nj9bt4aS/wpWzofeZ0LZ71D8yKYrH7gysrPO8FDgqjH06A6ujGJeISOPStjv88J6YfFSjGGg2s3Fm\nNsvMZq1bty7ocEREmqxoJoVVQNc6z7v42w50H5xzDzrnBjrnBubm5kY8UBER8UQzKcwEDjKzfDNL\nAcYCr++2z+vAj/0qpCHAJuecbh2JiAQkamMKzrlqM5sAvA0kAo865z43s8v81+8H3gROAoqBCuAn\n0YpHRET2L5oDzTjn3sT7h7/utvvr/OyAK6IZg4iIhK9RDDSLiEhsKCmIiEiIkoKIiISYi6PFHcJh\nZuuA5fV8ew5QFsFwIiVe44L4jU1xHRjFdWCaYlzdnXP7relvdEmhIcxslnNuYNBx7C5e44L4jU1x\nHRjFdWCac1y6fSQiIiFKCiIiEtLcksKDQQewF/EaF8RvbIrrwCiuA9Ns42pWYwoiIrJvze1KQURE\n9qHJJIWGLP25v/dGOa5z/Xjmm1mRmfWt89oyf/tcM5sV47iOMbNN/mfPNbNfh/veKMd1bZ2YFphZ\njZll+a9F8/t61MzWmtkeV1MP8PzaX1xBnV/7iyuo82t/ccX8/DKzrmY22cwWmtnnZnb1HvaJ3fnl\nnGv0D7yGeyVADyAF+Azotds+JwETAQOGANPDfW+U4yoE2vo/n7grLv/5MiAnoO/rGODf9XlvNOPa\nbf+Tgfej/X35xz4a6A8s2MvrMT+/wowr5udXmHHF/PwKJ64gzi+gI9Df/7k18GWQ/341lSuFhiz9\nGc57oxaXc67IObfBfzoNb02JaGvI3znQ72s3ZwPPRuiz98k59xGwfh+7BHF+7TeugM6vcL6vvQn0\n+9pNTM4v59xq59wc/+ctwCK8FSjritn51VSSwt6W9Qxnn3DeG8246roI77eBXRzwnpnNNrNxEYrp\nQOIq9C9VJ5rZ4Qf43mjGhZmlAWOAl+psjtb3FY4gzq8DFavzK1yxPr/CFtT5ZWZ5QD9g+m4vxez8\nimrrbAmfmY3C+592eJ3Nw51zq8ysHfCumX3h/6YTC3OAbs65rWZ2EvAqcFCMPjscJwNTnHN1f+sL\n8vuKazq/DljMzy8za4WXhK5xzm2O1HEPVFO5UmjI0p9hLQkaxbgwsz7Aw8CpzrnyXdudc6v8P9cC\nr+BdKsYkLufcZufcVv/nN4FkM8sJ573RjKuOsex2aR/F7yscQZxfYQng/NqvgM6vAxHT88vMkvES\nwtPOuZf3sEvszq9ID5oE8cC74lkC5PPtYMvhu+3zfb47UDMj3PdGOa5ueCvPFe62PR1oXefnImBM\nDOPqwLfzWAYDK/zvLtDvy98vA+++cHosvq86n5HH3gdOY35+hRlXzM+vMOOK+fkVTlxBnF/+3/tJ\n4PZ97BOz86tJ3D5yDVj6c2/vjWFcvwaygXvNDKDaeQ2v2gOv+NuSgGecc2/FMK4zgMvNrBrYDox1\n3lkY9PcFcBrwjnNuW523R+37AjCzZ/EqZnLMrBT4DZBcJ66Yn19hxhXz8yvMuGJ+foUZF8T+/BoG\nnA/MN7O5/rYb8BJ6zM8vzWgWEZGQpjKmICIiEaCkICIiIUoKIiISoqQgIiIhSgoiIhKipCAiIiFK\nCiIiEqKkINJAZnZZnR78S81sctAxidSXJq+JRIjfv+Z94K/OuTeCjkekPnSlIBI5d+AtyqKEII1W\nk+h9JBI0M7sQ6A5MCDgUkQbR7SORBjKzAcATwAj37SpnIo2Sbh+JNNwEIAuY7A82Pxx0QCL1pSsF\nEREJ0ZWCiIiEKCmIiEiIkoKIiIQoKYiISIiSgoiIhCgpiIhIiJKCiIiEKCmIiEjI/wNwzwS1P9aB\nqwAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "Co-courant\n", "X|z=0 = 1.0 X|z=1 = 0.714660643197\n", "Y|z=0 = 0.01 Y|z=1 = 0.580678713606\n", "A(X0-X1) 5.70678713606e-06 S(Y1-Y0) 5.70678713606e-06\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEKCAYAAAD9xUlFAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl8XXWd//HXJ/u+J02bJk2670AJZUcQKi1SUEYWxXXU\njiij46gzjDMjo6MooiOOMw6D/EAYBRRFQKTsW9lKWyjd9zVpmzRLkzRLs31/f3xv05C26U2bm5Pl\n/Xw8zuPee+659356H6d533O+3/P9mnMOERERgKigCxARkcFDoSAiIl0UCiIi0kWhICIiXRQKIiLS\nRaEgIiJdFAoiItJFoSAiIl0UCiIi0iUm6AL6KicnxxUXFwddhojIkLJixYoq51zuibYbcqFQXFzM\n8uXLgy5DRGRIMbOd4Wyn00ciItJFoSAiIl0UCiIi0mXItSmIiAShra2NsrIyWlpagi6lVwkJCYwd\nO5bY2NiTer1CQUQkDGVlZaSmplJcXIyZBV3OMTnnqK6upqysjJKSkpN6j4idPjKze82s0szWHOd5\nM7P/NLMtZrbKzOZEqhYRkVPV0tJCdnb2oA0EADMjOzv7lI5mItmm8Gtgfi/PLwAmhZZFwP9EsBYR\nkVM2mAPhsFOtMWKh4Jx7FajpZZOrgQec9xaQYWajI1XPvroWfvjUenbXNEXqI0REhrwgex8VALu7\nPS4LrTuKmS0ys+Vmtnz//v0n9WFv76jhnte284E7XmLRA8t5Y2sVmp9aRIaK3bt3U1JSQk2N/61d\nW1tLSUkJO3bs6NfPGRJdUp1zdzvnSp1zpbm5J7xK+5iuOm0MS/7hEm66eALLdtTwiV8t5fI7X+X/\n3trJwUPt/VyxiEj/Kiws5KabbuKWW24B4JZbbmHRokX097A/QYZCOVDY7fHY0LqIGZORyLcun8qb\n/3QpP/7YbOJiovjXx9Zwzm0v8J3H17CpoiGSHy8ickq+/vWv89Zbb3HnnXfy2muv8c1vfrPfPyPI\nLqlPADeb2cPA2UCdc27vQHxwQmw015UWcu2ZY3l39wEeeGMHD7+9mwfe3Mnckiw+ec44Lp8xiviY\n6IEoR0SGmO/+eS3r9tT363tOH5PGrQtn9LpNbGwsd9xxB/Pnz+fZZ5896WsRehOxUDCzh4CLgRwz\nKwNuBWIBnHN3AU8BVwBbgCbgc5GqpZcamVOUyZyiTP71ykM8sqKMB5fu4qsPvUtWchzXnjmWG+YW\nUZKTPNCliYgc0+LFixk9ejRr1qxh3rx5/f7+NtQaW0tLS10kR0nt7HS8tqWKB5fu4rn1FXR0Os4Z\nn8XH5xZx+Yx8EmJ19CAyEq1fv55p06YFWsPKlSu58cYbWbx4MRdccAFLly5l9OijO20eq1YzW+Gc\nKz3RZ+iK5h6iooyLJudy0eRcKutbeGRFGQ8v28XXHl5JemIsHz2jgOtKC5k+Ji3oUkVkBHHOcdNN\nN3HnnXdSVFTEt771Lb75zW/y29/+tl8/Z0j0PgpKXloCX7lkIq988xJ+8/mzuWhyLg8u3cUV/7mE\nhb94jQfe3MGBptagyxSREeBXv/oVRUVFXaeMvvzlL7N+/XpeeeWVfv0cnT7qowNNrTz2bjm/X17G\nur31xEVHMW/6KP7qzAIumpRLTLRyVmQ4Ggynj8Kl00cDKCMpjs+eX8Jnzy9h7Z46HllexuMry/nL\n6r3kpsbzkdPH8NEzxur0kogMSQqFUzBjTDozrkrn21dM46WNlfxhRRn3vb6DXy3ZztT8VK6ZU8BV\npxWQn54QdKkiImFRKPSDuJgoLp+Rz+Uz8qlpbOXJVXt49J1ybntqAz9cvIFzx2fzkdMLuHxmPumJ\n/d+vWESkvygU+llWchyfPreYT59bzPaqRh57t5zHV5bzD39cxb88toaLp+Ry9ekFfHBqHolx6t4q\nIoOLQiGCSnKS+fq8yfzdZZN4r6yOJ1bu4c+r9vDsugqS4qKZN30UC2eP4cLJObp6WkQGBYXCADAz\nTi/M4PTCDP75w9NYur2aP7+3l8Vr9vL4yj2kxscwb8Yorpw9mvMnKiBEJDjqPznAoqOM8ybk8MNr\nZrHsny/j1587i/kz83l+XQV//evllH7/ef7+dyt5bl0FLW0dQZcrIoOEc44LLriAxYsXd6175JFH\nmD+/t7nM+k7XKQwSre2dvLZlP0+t3seza/dR39JOclw0l0zNY8HM0Vw8JZfkeB3YiQRlMFynsGbN\nGq699lreffdd2tvbOeOMM3j66aeZMGHC+7bTdQrDQFxMFB+cOooPTh1F60dn8ea2ap5es5dn11bw\n5Kq9xMVEceHEHC6fkc+l0/LITokPumQRGWAzZ85k4cKF3H777TQ2NvLpT3/6qEA4VQqFQSguJooP\nTM7lA5Nz+f5HHMt21PDM2n08u7aCFzZUEmVQOi6LedNHMW/6KIo1iqvIwFp8C+xb3b/vmT8LFvzo\nhJvdeuutzJkzh7i4OCJx1kShMMhFRxnnjM/mnPHZfOfK6azdU8+z6yp4bl0FP3hqPT94aj0T81K4\nbNooLpuWxxlFmURHDf7JxUXk5CQnJ3P99deTkpJCfHz/nzFQKAwhZsbMgnRmFqTz9/Mms7umiefX\nV/D8+gruWbKNu17ZSmZSLBdPyeOSqXl8YFIu6Um6WE6k34Xxiz6SoqKiiIqKTD8hhcIQVpiVxOfO\nL+Fz55dQ39LGkk1VvLChgpc37udP75YTHWWcWZTJJVPzuHhKLlPzUzHTUYSIHJ9CYZhIS4jlw7NH\n8+HZo+nodKzcfYCXN1by4oZKbn96A7c/vYH8tATfVjEll/Mn5mjIDRE5irqkjgCV9S28vGk/L2+s\nZMnmKhpa2okyOL0wg4sm53LhpFxOG5uuYb9FejEYuqSGS11SpVd5aQlcV1rIdaWFtHd08l7ZAV7e\nuJ9XN1fx8xc2c+fzm0lNiOG8CdlcMCmXCyfmMC47SaeaREYghcIIExMdxZnjsjhzXBbf+NAUahtb\neW1LFa9truK1LVU8s7YCgIKMRM6fmM35E3M4d0I2eaka/ltkJFAojHCZyXEsPG0MC08bg3OO7VWN\nvL6lite3VPP0mn38fnkZAJPyUjhvQjbnTsjm7JJsMpPjAq5cZOA55wb9EfSpNgkoFKSLmTE+N4Xx\nuSl86txiOjoda/fU8cbWal7fUsXvlu/m/jd3YgZT89M4Z3wWZ5dkc3ZJlkJChr2EhASqq6vJzs4e\ntMHgnKO6upqEhJM/sldDs4Sttb2TVWUHeHNrNW9uq+adXbW0tHUCMDU/lbklPiTOKsnU6SYZdtra\n2igrK6OlpSXoUnqVkJDA2LFjiY19f+/CcBuaFQpy0g61d7CqrI6l26pZur2GFTtraWr1I7sWZydx\nVnEWZxVnUVqcSUlO8qD9dSUyEigUZMC1dXSydk89y7bX8PaOGpbvqKG2qQ2A7OQ4zhyXyZnjMikt\nzmRmQbrmjRAZQAoFCVxnp2Nb1UGW7ahl2Q5/JLGzugmAuOgoZhakMacokznjMplTlEl+uk45iUSK\nQkEGpf0Nh1ixs5Z3dtXyzs5aVpXX0dru2yVGpydwRlFGaJa6TGYVpGsea5F+oovXZFDKTY1n/sx8\n5s/MB3y7xPq9Dby7q5Z3dh1g5e5anlq9D/AjxE4elcrphenMHpvBaWMzmDwqRVdei0SQjhRk0Kk6\neIiVuw7wXtkBVu4+wHu7D1Df0g5AQmwUM8akM6sgndlj/VKSk6LhwkVOQKePZNhwzrGjuolVZQdY\nVVbHqrIDrCmvpzk0h3VSXDTTR6d1DSs+qyCdCbnJOqIQ6UahIMNaR6dj6/6DrCqrY015HavL61i3\n50hQxMdEMTU/lelj0pkxJo3pY9KYmp9KUpzOmMrINChCwczmAz8HooF7nHM/6vF8DvAbYDS+feMn\nzrn7entPhYIcT0enY3vVQdaU17OmvI61e+pZu6eu69STGZTkJDNtdBrT8lP97eg0Rqcn6BoKGfYC\nDwUziwY2AfOAMmAZ8HHn3Lpu2/wbkOic+0czywU2AvnOudbjva9CQfrCOUdZbTPr99azbm896/bU\ns35fPbtrmru2SUuIYepofyQxJT+VqfmpTB6VSmqC5puQ4WMw9D6aC2xxzm0LFfQwcDWwrts2+4DZ\n5n+mpQA1QHsEa5IRxswozEqiMCuJD83I71rf0NLGhn0NbNhbz/rQ7R9XlNEYuiIb/EixU/JTmTQq\nhSmjfFBMzEshIVbdZGX4imQoFAC7uz0uA87usc2vgBeAPUAqcL1zrrPnG5nZImARQFFRUUSKlZEl\nNSG2axiOwzo7HeUHmtmwr4FNFQ1s3OeXJZv309bhj6jNoCgriUl5KUzMSw3dpjAhL4WUeLVXyNAX\n9F78T8Aq4BJgAvCcmS1xztV338g5dzdwN/jTRwNepYwIUVFHjirmTR/Vtb6to5Od1Y1s3HeQzZUN\nbK44yKaKBl7ZdCQswF98NzEvhQm5PiQm5CYzMTeF3NR4tVnIkBHJUCgHCrs9Hhta1935wG3ON2xs\nMbPtwFTg7QjWJdInsdFRTMxLZWJeKr5PhOfDooktlQ1s3d/IlsqDbKk8yO+X7+4aGBAgNT6Gktxk\nxuckMz43hZKc5K4lWUcXMshEco9cBkwysxJ8GNwAfKLHNhuAS4ElZjYKmAJsi2BNIv3Gh4U/fdSd\nc4599S1srWxk6/6DbNt/kG1VjSzbUctjK/e8b9u81HiKc3xgFOckU5ydTHFOEuOykjXEhwQiYqHg\nnGs3s5uBZ/BdUu91zq01sy+Fnr8LuA24z8xWAVHAPzrnqiJVk8hAMDNGpycyOj2RCyblvO+5lrYO\ntlc1sm1/IzuqG9le1ciOqkaeX19B1cH3d7oblRbPuKxkxmUnMS47iaLsZMZl+fvpibE6JSURoYvX\nRAaJ+pY2dlQ1srO6iZ3VjWyvamJXjX9c2XDofdumJsRQlJVEUagNpDAricLMRAqzkijISFQPKTnK\nYOiSKiJ9kJYQy+yxGcwem3HUc02t7eyuaWZndSO7aprYXdPErpomNlY08MKGyq6RZg/LS42nMCuJ\nsZmJjM1MpCDD3y/ITFRoSK8UCiJDQFJcDFNCF9f11NnpqGw4xK6aJspqmyirbWZ3jb99Z1ctf1m1\nl/bO958RyEmJoyAjkTEZiV23fklgdHoiOSlxOj01QikURIa4qCgjPz2B/PQE5pZkHfV8e0cnlQ2H\nKKttpvxAE2U1zeypa6astpmNFQ28tLGya67tw+JioshPS2B0egJjMhIZnZ4QWhLJD93PSlZwDEcK\nBZFhLiY6qutIAI4ODeccB5raKD/QzJ7QsreuhT11Lew90Mzb22uoqG856mgjLiaKUWnx5KclMCq0\n5KclMCo9gVGp8eSn+3U6VTW0KBRERjgzIzM5jszkOGYWpB9zm45OR/XBQ+ypa2FfXQv76prZW99C\nRV0L++pbWFNex/PrK4464gA/tlReWgJ5qfGMCt0efpwbWvJS40mJj9GRxyCgUBCRE4qOMv+HPC3h\n/ZekduOco76lnYr6ltByiIr6FipD9ysbWli2o4bK+kO0dhwdHgmxUT4kUnxQ5HS79ffjuu4nxUUr\nQCJEoSAi/cLMSE+MJT0xlsmjjm4QP8w5R11zG5UNh9jf4MNif9d9f7s9dLFfTeOxB0xOjI0mOyWO\n7JR4cpLjyE6JIys5npyUI/ezk+PICi06hRU+hYKIDCgzIyMpjoykuF7DA/xQIjWNrT40Dh6i+mAr\n1QcPURW6vz90SmvNnjpqGlvfNxZVd8lx0WSlxJGV5EMiM9nfzwyFRubh9UmxodpiiR2hM/cpFERk\n0IqNjupqxD4R5xz1ze1UNx6iprGV6sZWakJL9cFWahoPUdPUxv6Dh9hUcZCaxtaumfqOJTUhhoyk\nWDJDAZaZFEtG4pHQyEiKJSMxjvRu69MSYob8NLAKBREZFsyM9KRY0pNiGZ8b3muaWzuobWr1S2Nb\n1/2axlYONLVxoKmV2tDtzupGahtbu2byO56U+Jiu02jdl7TEmG73Q0tCLOmJMaQl+MfxMVGBt5Uo\nFERkxEqMiyYx7nB33fB0dDrqm9s40OxDpK65jbpQcNQ1t3Og2a+rb26jrrmNrfsPUt/i7x+rd1Z3\ncdFRpCbEkJYYS2pCjF/ifaCkJsRywaQcLpmSd6r/7F4pFERE+iA66kgX3hKS+/TalrYOGlrau0Ki\nvrmN+pZ26prbaGhpo77ZP9fQ0h56ro2K+kM0hNYlxUUrFEREhouE2GgSYqPJTY0/qdcPxACmQ7tF\nRERkBBmI9gaFgoiIdFEoiIhIF4WCiIh0USiIiEgXhYKIiHRRKIiISBeFgoiIdAn74jUzKwDGdX+N\nc+7VSBQlIiLBCCsUzOx24HpgHXB4WEEHKBRERIaRcI8UPgJMcc4dimQxIiISrHDbFLYBsZEsRERE\nghfukUITsNLMXgC6jhacc1+NSFUiIhKIcEPhidAiIiLDWFih4Jy738zigMmhVRudc22RK0tERIIQ\nbu+ji4H7gR2AAYVm9hl1SRURGV7CPX30U+BDzrmNAGY2GXgIODNShYmIyMALt/dR7OFAAHDObSKM\n3khmNt/MNprZFjO75TjbXGxmK81srZm9EmY9IiISAeEeKSw3s3uA34Qe3wgs7+0FZhYN/DcwDygD\nlpnZE865dd22yQB+Ccx3zu0ys8hOPioiIr0K90jhJvzVzF8NLetC63ozF9jinNvmnGsFHgau7rHN\nJ4BHnXO7AJxzleEWLiIi/S/c3keHgP8ILeEqAHZ3e1wGnN1jm8lArJm9DKQCP3fOPdCHzxARkX7U\nayiY2e+dc9eZ2Wr8WEfv45yb3Q+ffyZwKZAIvGlmb4XaLLrXsQhYBFBUVHSKHykiIsdzoiOFr4Vu\nrzyJ9y4HCrs9Hhta110ZUO2cawQazexV4DTgfaHgnLsbuBugtLT0qHASEZH+0WubgnNub+jul51z\nO7svwJdP8N7LgElmVhK68O0Gjr4q+nHgAjOLMbMk/Oml9X3/Z4iISH8It/fRPOAfe6xbcIx1XZxz\n7WZ2M/AMEA3c65xba2ZfCj1/l3NuvZk9DawCOoF7nHNr+vqPEBEZNlqboG63Xw7shrqy0OMymPFR\nmPvFiH78idoUbsIfEYw3s1XdnkoFXj/RmzvnngKe6rHurh6P7wDuCLdgEZEhrbUJDuwKLTtDyy4f\nAAd2QVPV+7e3aEgrgPSxEJMQ8fJOdKTwILAY+CHQ/eKzBudcTcSqEhEZqjo7oWEP1O44xrITGnv0\nvI+Oh4xCyCiC0bMhPXQ/vdCvTx0NUdEDVn6voeCcqwPqgI8DhC4uSwBSzCzl8PUFIiIjSnur/1Vf\ns80vtduhZru/f2AndLQe2daiIb0AMoth8uWQOQ4yiv0f/sxxkJwHUeFeMhZ54Q6ItxB/jcIYoBI/\nV/N6YEbkShMRCVBHO9TtguqtUL3F39ZsDf3h3wWu88i2cSmQWQJ5U2HKAh8AWSX+Nr0QoofOHGXh\nNjR/HzgHeN45d4aZXQJ8MnJliYgMkMZqqNoE1ZuhanMoBDb7X/6d3WYIiE+DrPEwZg7MutbfP7wk\n54JZcP+GfhRuKLQ556rNLMrMopxzL5nZnRGtTESkv3R2+F/3VZtg/0Z/W7XZ3zZ3ax6NjoOsCZAz\nGaZcATmTIHuiX5ecM2z+8Pcm3FA4YGYpwKvAb82sEmiMXFkiIieho9036O5fD/s3+ADYv8EHQHvL\nke2ScyFnCky/yt8e/uOfUTSgjbqDUbihcDXQDHwdP0JqOvC9SBUlItKrzk7fd79yXWhZD5Ub/C//\njkNHtksvgtwpUPIB/+s/d4q/TcoKrvZB7oShEBoC+0nn3CX4C8zuj3hVIiKHNR+AirV+qTx8ux5a\nDx7ZJm0s5E2DCRdD7jTf4JszBeJTAit7qDphKDjnOsys08zSQ11URUT6X2en79q5bzVUrAndrvVH\nBIclZMComXD6JyBvemiZCgnpwdU9zIR7+uggsNrMnqNbW4Jz7qsRqUpEhrf2Vn+uf+97sG8V7F3l\ng+Dwr3+L9uf5C8+Gsz7vgyBvOqSNGRGNvUEKNxQeDS3dabRSETmxthZ/2mfPSh8Ce1f60z+HL/CK\nSzny6z9/ll9yp0Fs5Id0kKOFGwoZzrmfd19hZl873sYiMkJ1tPmG3/J3YM+7fqlcB53t/vnETBh9\nGpxzE+TPhtGn+37+g+iK3pEu3FD4DPDzHus+e4x1IjJSOOev7i1fcWTZt/pI18+EDBhzBpz3VRhz\nug+AjCKd/hnkTjRK6sfx8yiXmFn3uRBSAQ2IJzKSNNf6P/xlK6BsGZQv9+sAYpP8H/2zvuCDoGCO\nH/ZBATDknOhI4Q1gL5AD/LTb+gb8HAgiMhx1dvo+/7uXwu63fQhUbQw9ab7759QrYWwpFJRC7lSI\nDvfEgwxmJxoldSewEzh3YMoRkUC0NvmjgN1vwa6lUPY2tIR6oCdmwdizYPa1/nbMHEhIC7ZeiZhw\nR0m9BrgdyAMstDjnnPYMkaGoqQZ2vRla3vI9gw4P/pY7FaZ/xHcHLTwbsifoNNAIEu7x3o+Bhc45\nzZ8sMhQ17IOdr8OO12HnG35sIPADwBWcCefdDIXnQOFcDQExwoUbChUKBJEhpGEf7HgNdizxt9Vb\n/Pq4FP/rf9bHYNx5/lSQrgeQbsINheVm9jvgMaBrtCnnXM8L2kQkCI3VsONV2L4Etr/q5wMAPwfA\nuPNgzmeg+HzIP00NwtKrcPeONKAJ+FC3dY6jr3IWkYHQ2gg734RtL8G2V6BitV8flwJF58KcT0Hx\nhf5CsRE+FLT0TVih4Jz7XKQLEZFedHb4xuBtL8LWl31X0c423yZQeDZ88F/88NBjzhhSUz/K4BNu\n76PJwP8Ao5xzM81sNnCVc+77Ea1OZCSrK4etL8CWF2Dby9BywK/PnwXnfhnGX+wbh+OSAixShptw\nTx/9CvgW8L8AzrlVZvYgfu5mEekP7Yd8F9Etz8Pm54/0EEodDVM/DBM+6I8GUnKDrVOGtXBDIck5\n97a9v69yewTqERlZ6spg87M+BLa9DG2N/pTQuPP8qKETL/NXD+s6ARkg4YZClZlNIDRctpl9DD/8\nhYj0RWeHv3J409Ow6Rk/hwD4aSNPuwEmzfMNxJoxTAISbih8BbgbmGpm5cB24JMRq0pkODl0ELa+\nCBsX+6OCpio/iUzRuTDvezDpcj93sI4GZBAIt/fRNuAyM0sGopxzDZEtS2SIa9gHG5+CDU/B9lf8\nhDIJ6TDpQzB5Pky81M8tIDLIhNv76Dbgx865A6HHmcA3nHP/EsniRIaU6q2w/s+w4Uk/qihAZjGc\n9UWYsgCKzlF3URn0wj19tMA59+3DD5xztWZ2BaBQkJHLOT+r2LonfBhUrvXrR58Gl/yL7zGkRmIZ\nYsINhWgzi3fOHQIws0QgPnJliQxSzvl5htc97pearYD59oH5P/JBkFEUdJUiJy3cUPgt8IKZ3Rd6\n/Dng/hO9yMzm46fsjAbucc796DjbnQW8CdzgnPtDmDWJDAznYN8qWPsnWPsY1G73DcUlF/rRRade\nCSl5QVcp0i/CbWi+3cxWAZeGVv27c+6Z3l5jZtHAfwPzgDJgmZk94Zxbd4ztbgee7WvxIhFVuR7W\n/BHWPOqPCCwaxn8ALvi6D4Lk7KArFOl3YQ+X6JxbDCzuw3vPBbaEei5hZg8DVwPremz3t8AfgbP6\n8N4ikVG7A1b/wQdB5VqwKH/dwPlfhWlXaa4BGfYiOfNaAbC72+My4Owe71sAfBS4BIWCBKWxyp8a\nWvV7Pw0l+EHmFtwBMz6iU0MyogQ989qdwD865zqtlx4aZrYIWARQVKRGPOkHbc3+OoL3fucHnets\nh7wZcOmtMPOvIHNc0BWKBCKSM6+VA4XdHo8NreuuFHg4FAg5wBVm1u6ce6z7Rs65u/FXVFNaWur6\nWIeI55wfcG7lg77n0KF6SB0D534FZl0H+TODrlAkcJGceW0ZMMnMSvBhcAPwie4bOOdKDt83s18D\nT/YMBJFTVrsT3nsY3nvQtxnEJsP0q/1YQ8UXQlRU0BWKDBoRm3nNOdduZjcDz+C7pN7rnFtrZl8K\nPX/XyZUsEoa2Zn9B2bv/56enxKDkIrj4n2DaQohLDrpCkUHJnBtaZ2NKS0vd8uXLgy5DBiPnYO9K\neOf/fA+iQ3V+mInTb/RHBbqoTEYwM1vhnCs90Xbh9j4aC/wCOD+0agnwNedc2cmXKNJPWup8z6F3\n7od9qyEmwZ8eOuNTMO58nR4S6YNwTx/dBzwIXBt6/MnQunmRKErkhJzz8xIsv89fYNbe7KepvOIn\nMOtaSMwIukKRISncUMh1zt3X7fGvzezvIlGQSK8ONcDqR2D5vf6oIC4FTrse5nzGT1qvwedETkm4\noVBtZp8EHgo9/jhQHZmSRI6hYh0s/3++F1HrQX9UcOXP/FFBfGrQ1YkMG+GGwl/j2xR+hu919Abw\n2QjVJOJ1tMGGv8Dbd8PO1yE6HmZeA6Wfh7GlOioQiYBwQ+F7wGecc7UAZpYF/AQfFiL9q7EKVtwH\ny+6Fhj2+19Bl3/UNxxqETiSiwg2F2YcDAcA5V2NmZ0SoJhmp9q2Gt+7ybQYdh2D8JXDlf/gpLKOi\ng65OZEQINxSizCyzx5FC2COsihxXZwdsegbe+iXsWAKxSXDGJ2HuIsibGnR1IiNOuH/Yfwq8aWaP\nhB5fC/wgMiXJiNDa6McgeuuXULMN0sbCvO/BnE9rQnuRAIU7yc4DZrYc+GBo1TU9J8sRCcvBSlj6\nv74nUXMtFJwJH7sXpl0N0Tr4FAlaXybZWcfRE+SIhKdqM7zxC9+ltKPVz2V83t/6eQvUi0hk0NBP\nM4ms8hXw2s9g/ZMQHQenf8KHQfaEoCsTkWNQKEj/cw62vQyv/YcfoTQhHS78Bpz9N5rFTGSQUyhI\n/+nshE2LYclP/RFCSj586Ptw5md11bHIEKFQkFPX2QHrHoNXf+onu88shivv9KeKYuKDrk5E+kCh\nICevo92PUPrqHVC9GXImw0fv9nMcqyeRyJCk/7nSdx3tsOYPoTDY4ie8v/bXvlup5i4QGdIUChK+\nzg5/ZPAp6vayAAAMBUlEQVTK7T4MRs2C638DUz6sMBAZJhQKcmKdnb7N4OUfQdVGGDVTYSAyTCkU\n5Picg01Pw4s/gIrVkDsVrr0fpl2lMBAZphQKcmzbX4UXvgdlyyCzBK75lW9A1milIsOaQkHeb8+7\n8Px3YdtLkDoGFv4cTr8RomODrkxEBoBCQbzqrfDi92Hto5CYBR/6AZz1BYhNCLoyERlACoWRrrHK\n9yZafq8fm+iib/mxiRLSg65MRAKgUBip2pr9XAZLfgZtTX4eg4tvgdT8oCsTkQApFEaazk4/3eUL\n34X6ct+t9LJ/g9zJQVcmIoOAQmEk2fkGPPNt35g8+nTfo6j4/KCrEpFBRKEwEtTugOdu9RegpRX4\n8YlmXatrDUTkKAqF4ezQQT+nwRv/5a8vuPjbvhE5LinoykRkkFIoDEfOweo/wHP/Cg17Yfb1cOmt\nkF4QdGUiMsgpFIabvatg8T/Arjd9u8F1D0Dh3KCrEpEhIqInlc1svpltNLMtZnbLMZ6/0cxWmdlq\nM3vDzE6LZD3DWnMt/OWbcPcHoGozXPUL+OJLCgQR6ZOIHSmYWTTw38A8oAxYZmZPOOfWddtsO/AB\n51ytmS0A7gbOjlRNw1JnJ7z3kD9V1Fzrr0K+5NuQmBl0ZSIyBEXy9NFcYItzbhuAmT0MXA10hYJz\n7o1u278FjI1gPcNPxVr4yzf8qaLCs+GKn8Do2UFXJSJDWCRDoQDY3e1xGb0fBXweWBzBeoaP1kZ4\n+Yfw5i/9cBRX/ZcftE5dTEXkFA2KhmYzuwQfChcc5/lFwCKAoqKiAaxsENq4GJ76FtTt9kNTXPZd\nSMoKuioRGSYiGQrlQGG3x2ND697HzGYD9wALnHPVx3oj59zd+PYGSktLXf+XOgTU7/W9itY/AbnT\n4HNPw7hzg65KRIaZSIbCMmCSmZXgw+AG4BPdNzCzIuBR4FPOuU0RrGXo6uyEFffB8/8G7Yfgg/8K\n530VYuKCrkxEhqGIhYJzrt3MbgaeAaKBe51za83sS6Hn7wK+A2QDvzQzgHbnXGmkahpyqjbDE1+F\nXW9AyUVw5Z2QPSHoqkRkGDPnhtbZmNLSUrd8+fKgy4isjjZ44xfw8o/8JDeX3+Ybkn1wioj0mZmt\nCOdH96BoaJZu9q2Bx78Me9+D6VfDgjsgdVTQVYnICKFQGCzaW2HJT2HJT/yFZ9c94ENBRGQAKRQG\ng32r4U83QcVqmHUdLLhd3UxFJBAKhSB1tMFrP/NzJCdmwQ0PwdQrgq5KREYwhUJQ9m+EP/2NnwVt\n1rWw4Mc6OhCRwCkUBlpnJyz9H3j+uxCXDNfeDzM+EnRVIiKAQmFg1ZXBn74EO5bA5AWw8OfqWSQi\ng4pCYaCs/gP85e+hox0W/qcft0jXHYjIIKNQiLSWOj/5zerfw9i5cM3/Qtb4oKsSETkmhUIk7XoL\nHv0i1JXDxd+GC78B0frKRWTw0l+oSOho9xehvXI7pBfCXz+taTFFZEhQKPS3ujL44xf9IHazroMP\n/xQS0oKuSkQkLAqF/rT+SXj8K/6itI/+L5x2Q9AViYj0iUKhP7Qfgue+A0vvgtGnwcfu0xDXIjIk\nKRROVc12eOSzsHclnP0lmPc9iIkPuioRkZOiUDgV657wp4vM4PrfwrQrg65IROSUKBRORnsrPH8r\nvPVLGDMHrv01ZI4LuioRkVOmUOirunJ45DNQtgzm/g186N91ukhEhg2FQl9sexn+8Ne+Yflj98HM\na4KuSESkXykUwuEcvPYf8OL3IXsSXP8byJ0cdFUiIv1OoXAiLfXw2E2w4UmYcQ1c9QuITwm6KhGR\niFAo9Gb/Rnj4RqjZBpffBud8WSObisiwplA4ng1/gUf/xjcif/pxKLkw6IpERCJOodBTZye8ege8\nfBuMOcO3H6SPDboqEZEBoVDo7tBBP2/yhidh9g2w8E6ITQy6KhGRAaNQOKx2Bzz0Cdi/Hi7/IZxz\nk9oPRGTEUSgA7HgNfvcpcB3wyT/ChA8GXZGISCCigi4gcCvuhweuhqRs+MKLCgQRGdFG7pFCZ4cf\n7vrN/4IJl8LH7oXEjKCrEhEJ1MgMhUMN8McvwKan/fhFl9+muZNFRBiJoVBXDg9eD5Xr/FSZZ30h\n6IpERAaNiLYpmNl8M9toZlvM7JZjPG9m9p+h51eZ2ZxI1sOelXDPpb6n0Y2/VyCIiPQQsVAws2jg\nv4EFwHTg42Y2vcdmC4BJoWUR8D+RqoctL8B9V0BUDHz+GZh4WcQ+SkRkqIrkkcJcYItzbptzrhV4\nGLi6xzZXAw847y0gw8xGR6SazGIoOge+8DyMmhGRjxARGeoiGQoFwO5uj8tC6/q6Tf/IngCfehRS\n8yPy9iIiw8GQuE7BzBaZ2XIzW75///6gyxERGbYiGQrlQGG3x2ND6/q6Dc65u51zpc650tzc3H4v\nVEREvEiGwjJgkpmVmFkccAPwRI9tngA+HeqFdA5Q55zbG8GaRESkFxG7TsE5125mNwPPANHAvc65\ntWb2pdDzdwFPAVcAW4Am4HORqkdERE4sohevOeeewv/h777urm73HfCVSNYgIiLhGxINzSIiMjAU\nCiIi0kWhICIiXcyf1h86zGw/sPMkX54DVPVjOf1lsNYFg7c21dU3qqtvhmNd45xzJ+zTP+RC4VSY\n2XLnXGnQdfQ0WOuCwVub6uob1dU3I7kunT4SEZEuCgUREeky0kLh7qALOI7BWhcM3tpUV9+orr4Z\nsXWNqDYFERHp3Ug7UhARkV4Mm1A4lak/T/TaCNd1Y6ie1Wb2hpmd1u25HaH1K81s+QDXdbGZ1YU+\ne6WZfSfc10a4rm91q2mNmXWYWVbouUh+X/eaWaWZrTnO80HtXyeqK6j960R1BbV/naiuAd+/zKzQ\nzF4ys3VmttbMvnaMbQZu/3LODfkFP+DeVmA8EAe8B0zvsc0VwGLAgHOApeG+NsJ1nQdkhu4vOFxX\n6PEOICeg7+ti4MmTeW0k6+qx/ULgxUh/X6H3vgiYA6w5zvMDvn+FWdeA719h1jXg+1c4dQWxfwGj\ngTmh+6nApiD/fg2XI4VTmfoznNdGrC7n3BvOudrQw7fwc0pE2qn8mwP9vnr4OPBQP312r5xzrwI1\nvWwSxP51wroC2r/C+b6OJ9Dvq4cB2b+cc3udc++E7jcA6zl6BsoB27+GSyicytSfkZwStK/v/Xn8\nr4HDHPC8ma0ws0X9VFNf6jovdKi62MwOT2w9KL4vM0sC5gN/7LY6Ut9XOILYv/pqoPavcA30/hW2\noPYvMysGzgCW9nhqwPaviA6dLeEzs0vw/2kv6Lb6AudcuZnlAc+Z2YbQL52B8A5Q5Jw7aGZXAI8B\nkwbos8OxEHjdOdf9V1+Q39egpv2rzwZ8/zKzFHwI/Z1zrr6/3revhsuRwqlM/RnWlKARrAszmw3c\nA1ztnKs+vN45Vx66rQT+hD9UHJC6nHP1zrmDoftPAbFmlhPOayNZVzc30OPQPoLfVziC2L/CEsD+\ndUIB7V99MaD7l5nF4gPht865R4+xycDtX/3daBLEgj/i2QaUcKSxZUaPbT7M+xtq3g73tRGuqwg/\n89x5PdYnA6nd7r8BzB/AuvI5ch3LXGBX6LsL9PsKbZeOPy+cPBDfV7fPKOb4DacDvn+FWdeA719h\n1jXg+1c4dQWxf4X+3Q8Ad/ayzYDtX8Pi9JE7hak/j/faAazrO0A28EszA2h3fsCrUcCfQutigAed\nc08PYF0fA24ys3agGbjB+b0w6O8L4KPAs865xm4vj9j3BWBmD+F7zOSYWRlwKxDbra4B37/CrGvA\n968w6xrw/SvMumDg96/zgU8Bq81sZWjdt/GBPuD7l65oFhGRLsOlTUFERPqBQkFERLooFEREpItC\nQUREuigURESki0JBRES6KBRERKSLQkHkFJnZl7qNwb/dzF4KuiaRk6WL10T6SWj8mheBHzvn/hx0\nPSInQ0cKIv3n5/hJWRQIMmQNi7GPRIJmZp8FxgE3B1yKyCnR6SORU2RmZwL3Axe6I7OciQxJOn0k\ncupuBrKAl0KNzfcEXZDIydKRgoiIdNGRgoiIdFEoiIhIF4WCiIh0USiIiEgXhYKIiHRRKIiISBeF\ngoiIdFEoiIhIl/8PKIo78yhM9/cAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Transfert in a counter-current column \n", "# Shooting method to solve the problem for boundary conditions at both end\n", "from numpy import *\n", "from scipy import integrate\n", "from scipy import optimize\n", "kx,ky=1.e-4, 1.e-4\n", "K=2.\n", "sect, a=0.001, 100.\n", "A, S = 0.00002, 0.00001\n", "X_0, Y_1= 1., 0.01\n", "n=100\n", "L=2.\n", "z = linspace(0, L, n) \n", "#Set of ordinary differential equations\n", "def dX_dz(X, z):\n", " kgl=1./((1./kx)+(1./(ky*K)))\n", " N=-kgl*(X[0]-X[1]/K)*sect*a\n", " return array([N/A,N/S])\n", "Y0_try=1.\n", "# function giving the gap between dc/dz|z=1 and the targeted boundary conditions\n", "def cible(Y_0):\n", " X0 = array([X_0, Y_0]) # boundary conditions at t=0\n", " X, infodict = integrate.odeint(dX_dz, X0, z, full_output=True)\n", " return X[n-1,1]-Y_1\n", "#find the dc/dz|z=0 to have a cible=0\n", "Y0_found=optimize.newton(cible, Y0_try)\n", "X0 = array([X_0, Y0_found]) \n", "X, infodict = integrate.odeint(dX_dz, X0, z, full_output=True)\n", "print \"Contre-courant\"\n", "print \"X|z=0 = \", X[0,0], \" X|z=1 =\", X[n-1,0]\n", "print \"Y|z=0 =\", X[0,1], \" Y|z=1 =\", X[n-1,1]\n", "print \"A(X0-X1)\",A*(X[0,0]-X[n-1,0]),\"S(Y1-Y0)\", -S*(X[0,1]-X[n-1,1])\n", "print \"Efficacité dans le raffinat\", (X[0,0]-X[n-1,0])/(X[0,0]-(X[n-1,1]/K)), \"d'extraction\", (X[0,1]-X[n-1,1])/(K*X[0,0]-X[n-1,1])\n", "import matplotlib.pyplot as plt\n", "[X, Y]=plt.plot (z, X )\n", "plt.legend([X, Y],[\"X\",\"Y\"], loc='best')\n", "plt.xlabel ('z')\n", "plt.ylabel ('concentration')\n", "plt.show()\n", "# Transfert in a co-current column \n", "# Shooting method to solve the problem for boundary conditions at both end\n", "#Set of ordinary differential equations\n", "def dX_dz_co(X, z):\n", " kgl=1./((1./kx)+(1./(ky*K)))\n", " N=-kgl*(X[0]-X[1])*sect*a\n", " return array([N/A,-N/S])\n", "Y_0=Y_1\n", "X0 = array([X_0, Y_0]) \n", "X, infodict = integrate.odeint(dX_dz_co, X0, z, full_output=True)\n", "print \"Co-courant\"\n", "print \"X|z=0 = \", X[0,0], \" X|z=1 =\", X[n-1,0]\n", "print \"Y|z=0 =\", X[0,1], \" Y|z=1 =\", X[n-1,1]\n", "print \"A(X0-X1)\",A*(X[0,0]-X[n-1,0]),\"S(Y1-Y0)\", -S*(X[0,1]-X[n-1,1])\n", "import matplotlib.pyplot as plt\n", "[X, Y]=plt.plot (z, X )\n", "plt.legend([X, Y],[\"X\",\"Y\"], loc='best')\n", "plt.xlabel ('z')\n", "plt.ylabel ('concentration')\n", "plt.show()\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true, "deletable": true, "editable": true }, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] } ], "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.13" } }, "nbformat": 4, "nbformat_minor": 2 }