L2_AN_5_equations_5_inconnus.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "b824cdef",
   "metadata": {},
   "source": [
    " # Détermination de 5 inconnus avec 5 équations par bilan matière  \n",
    "\n",
    "Nous considérerons un séparateur avec 1 entrée (numérotée 1) et deux sorties (numérotée 2 et 3) fonctionnant avec un mélange binaire de constituants A et B dont les compositions sont données par des fractions, $x_A$ et $x_B$. Les débits globaux sont notés $F_T$ et les partiels, $F_A$ et $F_B$. \n",
    "\n",
    ">Le bloc de code ci-dessous permet de dessiner le flowsheet et de faire apparaître les débits globaux et partiels ainsi que les compositions des trois courants."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "id": "ec990393",
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQIAAAEeCAYAAABhWxuwAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAWD0lEQVR4nO3de7BdZXnH8e8vAQIBElIsVW5JAJUGwdQBQbzEqqDIrajDxanilTIiOOpQsQUJIBSjSIt0CJa2I7QKiI4zEBATEYKVOwTlklhQRBBMCeFmLpDk6R/v2rDY7HPO3mevtdfee/0+M2vOzrvf/a43JznPedbtfRQRmFm9Tah6AmZWPQcCM3MgMDMHAjPDgcDMcCAwMxwIzAwHAjPDgcDMcCAwMxwIzAwHAjPDgcDMcCAwMxwIzAwHAjPDgcDMcCAwMxwIbIBIepOkhyVtV/Vcho0DgQ2SecB2wFeqnsiwkRcvtUEg6U3Az4HNgDXALhHxaLWzGh7OCGxQzAM2zV5PwFlBoZwRWN9rygYanBUUyBmBDYJ8NtDgrKBAzgisr42QDTQ4KyiIMwLrd62ygQZnBQVxRmB9a4xsoMFZQQGcEVg/Gy0baHBWUABnBNaX2swGGpwVdMkZgfWrdrKBBmcFXXJGYH2nw2ygwVlBF5wRWD/qJBtocFbQBWcE1lfGmQ00OCsYp42qnoBZk+OBicDTLd4TMGWE9wAmAx8HvlrO1IaXMwLrK5J2AN44wtuvAb4NHDzKEDdHxBOFT2zIORDYwMiCxMMRoarnMmx8stDMHAjMzIHAbESSnpK0OtseGea1Eh0IzEa2e0RsFhGbAX8Evlf1hMriQGC1JelOSdFiWw4QEb/P+k0k3dfQ1Zl1SbtIulDS3ZLWS7p+lL6zJP1U0ipJf5B0ejaPjvq0y/cRWJ1NB54BvtDUvqzxIgsKr8r6/XWX+9sNeD9wM7DJSJ0kTQMWAfcBhwI7A+eQfnGf3G6fjkSEN28DsQE7pP+yhY23Hri9jX4bA7cAC7rc34Tc6yuA60fo92VgJTAl1/b3wKpGWzt9Otl8aGC1JGkW6bfnL8fqGxEvAGcC7+5mnxGxoc2uBwDXRsQzubZLSYcnczro0zYHAqurA7Ovt0qalNsEIGmGpN1y/T9HOmFI9r4kbTTWNs657QoszTdExMOk3/a7dtCnbQ4EVldvy75eQHpYqbHNz9qnk4LEaklrgNeRju8bjgZeaGMbj2nAUy3aV2bvtdunbT5ZaHU1i/SD+smm9kUAEXEDsPkon78S2KucqQGtr1Coqb2dPm1xILC62hZ4PCIuGefnn2TkpyC7tRLYqkX7VF7KAtrp0zYfGljtZNfaJ5O7TDgOZR4aLKXpOD974GpzXjov0E6ftjkjsDraL/t6axdjlHlocA1woqQtI+LZrO0IYDVwQwd92uZAYHX0nuzrT8Y7QESsAFZ08hlJk3nphON2wBRJH8r+fHVErMpezwdOAH4o6WvATsBc4Ju5y4Xt9Gl/btmNCGZ9r6j1CLJbe+eQbvDp2Q+ApBnAb0d4e2ZEPJTrOws4H3gL6Zj/ImBuRKzvpE/bc3MgsEHhhUnK45OFZuZAYGYOBGaGA4GZ4UBgZjgQmBkOBGaGA4GZ4UBgZvhZA7NRSbqB9LivgF8DnxjPvfz9zhmB2egOiYg3RsQewMPAiVVPqAwOBGajiIinASRNID3rP5QP5zgQmI1B0tWkhUtfD8yreDqlcCAwG0NEvB94NWkhk89UPJ1SOBBYLUmak5U3OyDXNlPScknnNffPnvH/DvDRXs6zVxwIrJayVYp/BpwCIGkqcBXpt/7ns7Zpkv4i97EPAvf0eKo94cuHVmenAosl7Q98kbTY6JG5FX6mAZdL2oR0+fB+4PhKZloyr1BkA6OMFYokLQT2JS31tXdEPFLU2IPEhwZWdw+QljY/ta5BABwIrMYkHQN8Argb+FTF06mUA4HVkqT9SCsAfxo4Ftg7fwWhqe9x2RWGnXo5x15yILDayaocfx+YFxEXR8TNpJqHp7XoO4V0IvEOYPeeTrSHfLLQBkYRJwslbQPcAtwOHN6oayDpHaQKQQdFxIJc/7OA50nnEZ6NiDO6+Cv0LV8+tFqJiOXAzBbti0mXCF8kaXvgSGAP4DDg4F7MsQo+NDAb2ZnA2RHxHOlGoj0qnk9pfGhgA6OXlY4kzQZuApaTnjgUqZT6lhGxpuz995ozArPWvgEcFRHTI2JGREwn1S3creJ5lcKBwKxJdhlxSkT8qOmtZQzp4YEPDWxguAhqeZwRmJkDgZk5EJgZDgRmhu8sNBuV6xqYGbiugZm5roGZAa5rYGa4roHZ0HJdg5dzILBacl2Dl/PlQ6sz1zXI+KEjGxiua1AeHxpY3bmuAQ4EVmOua/ASBwKrpXbqGkg6TNIKSUsk3Sfpbkmzqphv2RwIrHY6qGswG/hWRMyOiFnAbcBRPZ1sjzgQWK1kdQ2uAhaSXTrMnAHsJenAXNtsUmETJG0FvAG4sScT7TFfNbCB0eulyiT9DngOWA+8DjglIr7ei333mu8jMGtB0jRgi2z1YiS9Frhd0oV+DNmsPmYD9zX+EBH/S3rycGpVEyqTA4FZa7OBext/kHQ48EhE/L6yGZXIhwZmrb0ReK+kfUi3Hj8KHFLtlMrjk4U2MFzXoDw+NDAzBwIzcyAwMxwIzAyfLBxYM05acAzw4arn0UsbXlg76fnHfr3PpjvufkPVc6nAdx86+8BvlzW4M4LB9WHSte7a0EabrN34VTveXfU8KjCbkoO+7yMYbEseOvvAd1Y9CSvXjJMWXF/2PpwRmJkDgZk5EJgZDgRmhgOBmeFAYGY4EJgZDgRmhgOBmeFAYGY4EJgZDgRmhgOBmeFAYGY4EJgZDgRmhgOBmeFAYGY4EJgZDgRmhgNBpSTtVPUczMCBoDKS3gk8KOk2SXMqno7VnANBdTYBngH2BBY4IFiVHAiq1SgztTkOCFYhB4L+4oBglXAg6E8OCNZTDgT9zQHBesKBYDA4IFipCimCKmkWcG8RY9moGgHhmrWPP7B20qt3qWNlYCtBUdWQlwPXAn9e0Hh1MAuYCGzc4efWAkzYeNNnC5+R1VYhgSAingDeV8RYdSFpf+ByYGqbH1kLbAD+DThr4623v6ysuVn9+BxB/1sLrAYuBGZGxOci4o8Vz8mGTFGHBla8l2UA/uG3MlWSEUi6R9J6STF279pxBmA9V9WhwXxgj4r23a8cAKwypQQCSXdKihbbcoCIOD8iCrvcKOlDkn4haYWkNZKWSTpZ0iYt+s6S9FNJqyT9QdLpkiZ22qcgG+EAYH2grHME00lP1n2hqX1ZSfvbGvgZ8HXgKeDNwFzg1cBnG50kTQMWAfcBhwI7A+eQAuLJ7fYpyBOkh44uxOcArGJlBYKtgLsi4t9LGv9lIuLCpqafSZoCHCfp+IhonIs4FtgM+EBEPAMszPrNlTQva2unTxFzvlPSVhGxvojxzLpR+KFBdpfhBOCXRY/doRWkZ/7zDgCubfphvpT0gz+ngz6FcBCwflHGOYIDs6+3SpqU29TOh5VsNNY2wmcnSpos6W3ACcAFuWwAYFdgaf4zEfEwsCp7r90+ZkOljEDwtuzrBcCa3Da/0UHSUknrstfrJOV/8I4GXmhja+VP2XYjcANwYtP700jnEJqtzN5rt4/ZUCnjHMEs0g/qJ5vaFzVeRMRov1mvBPYa5773BSaTThZ+BTgf+ExTn1b3LqipvZ0+ZkOjjECwLfB4RFwyzs8/CTw9ng9GxJ3Zy59LegL4jqRzIuLBrH0l6URms6m8lAW008dsqBR6aJBda59Md5cJuzk0yGsEhZm5tqU0HedL2oH0eO/SDvqYDZWiM4L9sq+3djFGN4cGeW/Nvv4213YNcKKkLSOi8RjvEaQbem7ooI/ZUCk6ELwn+/qT8Q4QEStIl/7aJunHpHMQ9wLrSUHgi8BlucMCSCcsTwB+KOlrwE6kG4++mbtc2E4fs6FSdCDYM/u6uOBxx3Ib8DFgBrAO+A3wZXJXKgAiYqWkd5NOIl5JOuY/l/SD3nYfs2FTaCCIiHcWOV4H+z0FOKXNvvcB7+q2j9kw8cIkZuZAYGYOBGaGA4GZ4UBgZjgQmBkOBGaGA4GZ4UBgZlQYCCRdIOlR1zYwq16VGcH3gDdVuH8zy1RW8iwiFgO0uZShvdIcgBknLbi+4nn0TGzYoHVPP77jxtO2/V3Vc+mx2cCSMnfgcwQ2MGL9C5usf3bFjKrnUYElwHfL3IFevshv70mKiHBaYGPKVop62P9fildGXYM5WXmzA3JtMyUtl3Re0fszs+6VkhFIug7YNCL2lTQV+AVpybBDm4t6OCOwdjkjKE9Z5whOBd4iaX/gctJio0fmg4CkiyQ9kr1+RNJFJc3FrGOSdsgK4d4v6V5J89ot0jOISjtHIGkhqc7AU8DeEfFIKTuy2uhlRiDpNcB2EXF7VlV7IXBeRPyg7H1XocyrBg+QljY/1UHA+pGkc7PzWc3bgoh4LCJuB4iI50m1PHeodsblKescwTHAt4D7gTURsU/hO7HaKTojkLQImMQrS+M9HhEP5fptTbqEt39E3F/EvvtN4TcUSdqPtALwp4BfAzdJOiAirmnR97is784R8Zui52I2ht2B/4qIm0fqIGkScAXwz8MaBKD4Ske7Ad8H5kXExdk3eBFwWou+U0i1B+4g/YOY9YykbYBtSLUwRuozEfhv4K6IOKdXc6tCYYEg+8ZeRTqpkl9a/AxgL0kHNn3kJOBi4Dpgj6LmYdamxv+5ZZI2ym35w44LgWdJv7CGWmGHBhGxnJfXGWy0LyZVEn6RpO2BI0n/GIcBBxc1D7M2NQLBz5vajwYulvRWUkXve4C7svjwHxExlDfFVfXQ0ZnA2RHxnKR7gH+saB5WX7uTrmwd1dS+DCAi/oemX2DDrOfPGkiaDdwELAeC9M3eFtgyItb0dDI2UIq8aiDpduDBiDii+5kNviqePvwGcFRETI+IGRExnXT78W4VzMVqSNIEYBbwq6rn0i96GgiyB5GmRMSPmt5ahk8YWu+8FtgMB4IXVf4Yslm7ijo0kPQh0mVu37+S8cIkVke7A38iHZIazghsgPgx5PI4IzAzBwIzcyAwM6otcHKPpPUucGJWvSozgvn43gGzvlBZIIiI8yNixEdAOyVpVrbG3CpJf5B0evYY6Vif+9gIq9QcW8T4ZoOgskpHRZI0jbTuwX3AocDOwDmkQHdym8O8C1id+/OLN5oUNL5Z3xqKQAAcS7pl9AMR8QywMFv4ZK6keVnbWG6LiOdKHN+sb5VR4OSELLU+Jdf2dkkbJC0pen+ZA4Brm34gLyX98M4ZgPHNKlV4IMgWblhJtiBkdjfYQuD/gL3aGUPJRmNtuY/sCixtmsfDwKrsvXY8KGmdpGWS/q7pvSLGtwFTp9oGZZ0sPAXYUtJJwF2kdQf2iIgXGh0kLZW0Lnu9TlL+B+1oUlGUsbaGaaT6Cc1WZu+N5rFsvh8hrZR0CzBf0ucLGt8G1zrgSxHxl8BfAXsDH6h2SuUo5RxBRPyrpNOBfwI2kAqc/LGpz2i/Sa+kzewhP2SLNo3Qnp/HtcC1uaZrspVrT5b0LxGxoZvxbXBFxGOkXxRExPOShra2QZknCx8D/gy4pFEoogNPAk930H8lsFWL9qm0/k0+liuAw4EZpKsHRY9vAyarbfA3wP4VT6UUpRwaSLqEtOLQGuCD4xii00ODpTQdq2fnJjan6di+Q43f9mWNbwOgDrUNyrhq8CXgb4FvA58FtshfQcj6nJ1dRVgtaW32Nb+ScePQYKyt4RrgvZK2zLUdQbov4IZx/DU+CDwB/K6k8a1ikuZkV7cOyLXNlLRc0nm5tnrUNoiIwjbgENJv0RtzbSuA55r6XQ9cn/vzUmBxF/udRjoUWQi8BzgGeA74alO/j5JOAE3Ptf0A+BLpEuFBwCXZ3+H4Tsf3Vu5GOj6PAse7DvhF9noqqdjJVcDEXJ+LgP8kW7tjWLci/5FmkdL13+e/aaSsIIC5ubbHgNOz19NJRSROKmD/15F+Sz9GKqwysanPx7K5zMi1nUVaM3FV9tk7gI+MZ3xv5W4lBIK3Z/8f9iedMF4CbJF7/63Z+7/K3lsCnFD196GMrZIVirLLhuuyb/KmwIKIOKjnE7GBUsYKRZIWAvuSTvruHTWt3N3zh44kzQAmRMSmEbEZsB9woKTtej0XM1KRk8nAqXUNAlDN04eHkQ4FAIiIRdlLBwLrKUnHAJ8A7iZV766tKgLBHODFyCvpXGBtRNxawVyspiTtB5wPfJr0UNne+SsITX2Py64w7NTLOfZSFU8fvgGYKWk16a7DZ4D3VjAPqylJu5HqGsyLiIuztkXAaaRLxfm+U0jVkO8gLYM+lHUQvJy5DYwiThZK2ob0PMntwOGR/QBIegfpnpCDImJBrv9ZwPOk8wjPRsQZXfwV+tawrEdg1paIWA7MbNG+mKbqx5K2B44kLal3GOmhtKHkVYzNRnYmcHakBWvuYYjX2PShgQ2MXlY6kjQbuAlYTrrfRcC2wJYRsabs/feaMwKz1r4BHBUR0yNiRkRMJ9VK3K3ieZXCgcCsSXYZcUpE/KjprWUM6eGBDw1sYLgIanmcEZiZA4GZORCYGQ4EZoYDgdmIJL1Z0pPZcnprJN3iugZm9bMW+FxETAK2Bl4PzKt2SuVwIDAbQUTcHRGXZK//RFrMdpduxuy0qrak60eo1h2S3tLN2Hl+6MisDZJ2Id1V+A9djDGeqtqfAaY0tZ1Oqrx0W5djv8iBwGwM2TL2twFX5x9RHoeOq2pHxH1Nc9kE2BO4LCLWdTN2ng8NrJbardotaWPgfuA3EXFIl7stoqr2+0jL63+vyLEdCKyWov2q3b8iLWG/Z/MYFVXtPhJ4FLixqb2rsR0IrM5Grdot6VjSlYIdgFVZRa4rcp/vZdVuJE0mLY5yWbzyIaGuxvY5AqutGKNqd0TMB+aPMkTPqnZnDga24JWHBV2P7UBgdTdIVbuPBB4YYZ5dje1DA6utQaraLWkq6YTgSNlAVxW7HQislgawavdhwCRGDgTdVeyuuviiN2/tbhRUBJXBrNr9Y2BJt2OPtDkjsFqRNAv4Aana1jtyb50KbC5pbq7t9cDi7HPTSWX5rh7vviNiJfBuYCIpmzgNODfbd96ErI+yfb8q+9ylBYzdkpcqs4HR66XK6lS12xmBWQt1q9rtQGDWWq2qdjsQmLVWq6rdvqHIrLVaVe12IDBrISK6WoBk0PjQwMwcCMzMgcDMcCAwM3yy0PqMpH0Y+ez8Nlmf0W6bvTwi7i98YkPOgcD6zceBTzF6tjp3hPb1pKW5HAg65GcNrK9I2hFYRrq3v1PPAdtGxLNj9rSX8TkC6yuRFty8lJcv6NGO1cDXHATGxxmB9Z1xZgXOBrrgjMD6zjiyAmcDXXJGYH2pw6zA2UCXnBFYX+ogK3A2UABnBNa32swKnA0UwBmB9a02sgJnAwVxRmB9bYyswNlAQZwRWF8bJStwNlAgZwTW90bICpwNFMgZgfW9FlmBs4GCOSOwgdCUFTgbKJgzAhsIuaxgA84GCueMwAZGlhWcB3zEgaBYDgRm5kMDM3MgMDMcCMwMBwIzw4HAzHAgMDMcCMwMBwIzw4HAzHAgMDMcCMwMBwIzw4HAzHAgMDMcCMwMBwIzw4HAzHAgMDPg/wHyJOqIRiSW6AAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 360x360 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "xs, ys = 40, 40       #position du separateur\n",
    "lx, ly = 40, 20       #taille du separateur\n",
    "la, hla, hwa=25, 5,5  #taille des fleches\n",
    "\n",
    "#tracé du flowsheet\n",
    "def plot_flowsheet_sd():\n",
    "    fig, ax = plt.subplots(figsize=(5, 5)) \n",
    "    ax.set_aspect( 1 )\n",
    "    plt.axis('off')\n",
    "    plt.plot([xs,xs+lx,xs+lx,xs,xs],[ys,ys,ys+ly,ys+ly,ys])\n",
    "    #plt.plot([xs,xs+lx],[ys,ys+ly],'b--')\n",
    "    plt.arrow(xs-la-hla,ys+ly/2,la,0, head_width=hwa, head_length=hla, fc='k', ec='k')\n",
    "    plt.arrow(xs+lx/2,ys,0,-la, head_width=hwa, head_length=hla, fc='k', ec='k')\n",
    "    plt.arrow(xs+lx/2,ys+ly,0,la, head_width=hwa, head_length=hla, fc='k', ec='k')\n",
    "    #plt.text(45,48,'Séparateur',c='b', size=15)\n",
    "    plt.text(xs-la-hla,ys+ly/2-6, r\"$F^1$\", size=15)\n",
    "    plt.text(xs-la-hla,ys+ly/2-14, r\"$x^1_A$\", size=15)\n",
    "    plt.text(xs-la-hla,ys+ly/2-22, r\"$x^1_B$\", size=15)\n",
    "    plt.text(xs+lx/2+5,ys-6, r\"$F^2$\", size=15)\n",
    "    plt.text(xs+lx/2+5,ys-14, r\"$x^2_A$\", size=15)\n",
    "    plt.text(xs+lx/2+5,ys-22, r\"$x^2_B$\", size=15)\n",
    "    plt.text(xs+lx/2+5,ys+ly+25-6, r\"$F^3$\", size=15)\n",
    "    plt.text(xs+lx/2+5,ys+ly+25-14, r\"$x^3_A$\", size=15)\n",
    "    plt.text(xs+lx/2+5,ys+ly+25-22, r\"$x^3_B$\", size=15)\n",
    "\n",
    "#Les 4 paramètres fixés\n",
    "FT1=300\n",
    "XB1=0.5\n",
    "XB2=0.7\n",
    "FT3=100\n",
    "\n",
    "#Analyse de la cohérence du jeu de données\n",
    "if ((FT1*XB1)-(FT1-FT3)*XB2)/FT3 >1:\n",
    "    print ('Attention le jeu de données ne correspond pas à une solution physique')\n",
    "    print ('XB3>1, il y a soit :')\n",
    "    print ('      i)   trop de B en entrée (essayez de réduire XB1 ou FT1)' )\n",
    "    print ('      ii)  pas assez de B en sortie 2 (essayez d augmenter XB2 ou FT2=FT1-FT3)' )\n",
    "if ((FT1*XB1)-(FT1-FT3)*XB2)/FT3 <0:\n",
    "    print ('Attention le jeu de données ne correspond pas à une solution physique')\n",
    "    print ('XB3<0, il y a soit :')\n",
    "    print ('      i)   pas assez de B en entrée (essayez d augmenter XB1 ou FT1)' )\n",
    "    print ('      ii)  trop de B en sortie 2 (essayez de réduire XB2 ou FT2=FT1-FT3)' )    \n",
    "#Tracé des 4 paramètres fixés sur le flowsheet\n",
    "def plot_flowsheet_data():\n",
    "    plt.text(xs-la-hla,ys+ly/2-6, r\"$F^1$\"+'={}'.format(FT1), size=15)\n",
    "    plt.text(xs-la-hla,ys+ly/2-22, r\"$x^1_B$\"+'={:.2f}'.format(XB1), size=15)\n",
    "    plt.text(xs+lx/2+5,ys-22, r\"$x^2_B$\"+'={:.2f}'.format(XB2), size=15)\n",
    "    plt.text(xs+lx/2+5,ys+ly+25-6, r\"$F^3$\"+'={}'.format(FT3), size=15)\n",
    "def plot_flowsheet():\n",
    "    plot_flowsheet_sd()\n",
    "    plot_flowsheet_data()  \n",
    "plot_flowsheet()\n",
    "plt.show()"
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "85b16940",
   "metadata": {},
   "source": [
    "1- Qu’est ce qui vous semble manquer sur ce schéma ? Faites des hypothèses qui vous semblent cohérentes pour le compléter. \n",
    "\n",
    "> Il manque les unités. Si les fractions sont molaires alors les débits devront être exprimés en mole par unité de temps. Si les fractions sont massiques alors les débits devront être exprimés en masse par unité de temps. \n",
    ">\n",
    "> Considérons par exemple que les fractions sont molaires (nous utilisons dans ce cours x pour les fractions molaires) et que les débits sont en moles par heure.\n",
    "\n",
    "2- Ecrire 5 équations reliant les inconnues.\n",
    "\n",
    ">Les équations, permettant de décrire les bilans de matière, sont les suivantes :\n",
    ">\n",
    ">**Bilan sur le débit global**\n",
    ">1. $F^1_T=F^2_T+F^3_T$\n",
    ">\n",
    ">**Bilan sur les débits partiels de A et de B**\n",
    ">\n",
    ">2. $F^1_A=F^2_A+F^3_A$ soit $F^1_T x^1_A=F^2_T x^2_A+F^3_T x^3_A$\n",
    ">3. $F^1_B=F^2_B+F^3_B$ soit $F^1_T x^1_B=F^2_T x^2_B+F^3_T x^3_B$\n",
    ">\n",
    ">**Lien entre les différentes fractions molaires**\n",
    ">\n",
    ">4. $x^1_A+x^1_B=1$\n",
    ">5. $x^2_A+x^2_B=1$\n",
    ">6. $x^3_A+x^3_B=1$\n",
    ">\n",
    ">Il y a donc 6 équations pour 9 inconnues \n",
    ">\n",
    ">$F^1_T,F^2_T,F^3_T,x^1_A,x^2_A,x^3_A,x^1_B,x^2_B,x^3_B$. \n",
    ">\n",
    ">Il faut cependant noter que ces équations ne sont pas indépendantes puisque la première est la somme >des équations 2 et 3 en considérant les équations 4, 5 et 6.\n",
    ">\n",
    ">Il y a donc seulement **5 équations indépendantes** qui forment un système comportant **9 inconnues**. Il faut donc **fixer 4 variables** pour obtenir une solution unique et déterminer les 5 inconnues à l'aide de 5 équations.\n",
    ">\n",
    ">On a fixé dans cet exercice les 4 paramètres suivants : $F^1_T,F^3_T,x^1_B,x^2_B$. Les 5 équations nous permettront alors de déterminer les 5 inconnues : $F^2_T,x^1_A,x^2_A,x^3_A,x^3_B$\n",
    ">\n",
    ">Les lignes de codes ci-dessous permettent de visualiser les 4 données fixées et de vérifier si une solution \"physique\" existe. Si ce n'est pas le cas, un message s'affiche pour vous dire dans quel sens modifier les 4 données pour retomber sur un problème ayant une solution physique. Vous pouvez donc vous exercer en changeant le jeu des 4 paramètres fixés et refaire les calculs. \n",
    ">\n",
    ">A noter que si vous choisissez $x^1_B=x^2_B$, il s'agit d'un simple diviseur de courant qui ne modifie pas la composition des courants de sortie."
   ]
  },
  {
   "attachments": {},
   "cell_type": "markdown",
   "id": "e3c7a587",
   "metadata": {},
   "source": [
    "3- Déterminer les inconnues.\n",
    "\n",
    "> Les blocs suivant permettent d'appliquer les différents bilans (équations par équations) pour déterminer l'ensemble des inconnues. Vous pouvez donc exécuter les blocs successivement pour comprendre comment on arrive à déterminer l'ensemble des inconnues. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "id": "63c8e5a8",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAARYAAAExCAYAAAC55I3BAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAb8UlEQVR4nO3dedRcVZnv8e8viQECJARYtM1gAtgOiaTTCwRlMOkWRGRq0A5gK9IOXJYC3tsuGnAlAmEURboxaNrFvQ6xu1Ghly46DAaQSWQSwpwgIEIggoYwhCSQ4bl/7FOkUrzJW/XWrvH8Pmud9dZ7atc+O4R68px9ztmPIgIzs5yGdXoAZtZ/HFjMLDsHFjPLzoHFzLJzYDGz7BxYzCw7BxYzyy5vYJG+i/QsUudvjummsZiVjLLeICd9CFgI/JEI5eu4x8diVjJ5A8ubvSq65svcyrFI7wROAT4AvA+4lYipLTmWWQ8Z0ekB9LiJwMeAO4CRHR6LWdfojslbaWqPzoVcRcRORPwD8HCnB2PWLeoLLNIUpEA6qGrfzkgvIF3SqsF1/Vgi1rb1eGY9ov45FulGYFMi9kYaA9wO/B44nIg1NW0bm9eQpgK/qvszOcciCRg+6DEjVg8ypiuAbT3HYtbYqdAZwAeRPgL8FFgFHL3eF1m6DGlR8XoR0mUD9iQJacSbW+WLXb0v7W/9WOAzxecH28ysTo1dFZLmAXsDLwF7EbFoaEfVccD3B2238Uwj11i2AXauYyz3DNKPMxazQqNXhR4H9ge+POQvcnIV8P6q33cHZtfsa9dYXgRebuLzZlajkTmW44FvA48CK4n4QL5RaCqNzbHkG0uO7Cn144zFrFBfxiIdAMwCPg88BvwG6SAirhmg7ZeKtrsS8WS+oTYwFukI4DLgGdL9JauAY4h4ZIAea7MnM2vS4BmLNBH4NTCLiOnFvnnAGCL2rGk7GphPOr04m4hf1DcKTaWejKXesUhnAUHEmcXvlwGLiZhR13jqJY0i3SAH8BVgNGliGeBqIpZnPZ5Zj9h4YJG2A+4E7gGmUWmcnsO5GTiEiLlV7c8D3gBGAa8ScXZ9o9BUBgssjYxF+gVwGRFXIW0FXAt8jYhf1jWeeknjSZe5B7IzEU9lPZ5Zj8j3rJC0I3ALMAk4AjiUiGl5Om94LH8AlgFrgHcBM4j4RkfGYlZCOW/pPxe4gIhlwEOkANN+0lhgCyImEjEJ2A2YXpymmVkb5HkIUZoMTAOmIn0VELA90qZErMxyjPpNBtZN0kb8rngOaQzwSpvHYlZKuTKWb5KuuowjYjwR40hzDxMz9d+IyVQ/EChNAxYR8UwHxmJWSs1nLOlhwNFE/LzmnYWk06HfNn2Mxvw1cCDSB0iXmZ8FDmvzGMxKrTULPZlZqXXHeixm1lccWMwsOwcWM8vOgcXMsnNgMbPsHFjMLDsHFjPLzoHFzLJzYDGz7BxYzCw7BxYrLUm7Slok6b2dHku/cWCxMjsH+Evg/E4PpN/4IUQrJUm7khYk2xRYAeweEY92dlT9wxmLldU5rFs2ZCTOWrJyxmKlU5OtVDhrycgZi5VRdbZS4awlI2csViobyFYqnLVk4ozFymagbKXCWUsmzlisNAbJViqctWTgjMXKZGPZSoWzlgycsVgp1JmtVDhraZIzFiuLerKVCmctTXLGYn2vwWylwllLE5yxWBk0kq1UOGtpgjMW62tDzFYqnLUMUZ6i8Gbd63hgOPDyBt4fs5H3NgO+CJzUgnH1NWcs1tck/QXw/o00uQr4B2DlBt6/NyKeyz6wPufAYqUmKYBREbGi02PpJ568NbPsHFjMLDsHFrM2kfSSpBXFtkjSDp0eU6s4sJi1z24RsVlEbAY8D/xXpwfUKg4sZplIuldSDLC9ABARzxTthpMuZTd15UTSOyX9u6T7Ja2RdNNG2k6QdIOk5ZKekzSzGEdDberl+1jM8hkHvAL8c83+hZUXRZDZtmj3t00ebyLwMeAO0p3CA5I0FrgeeAQ4HNgVuIiUWEyvt00jfLnZSi3n5WZJa4D7ImKPQdq9DbgN+HNEHNzE8YZFxNri9RXAthExdYB2pwP/AoyLiFeKff8CnAm8PSJeqadNI2PzqZBZBpImkL5PDwzWNiJWAecCH27mmJWgUoeDgOtqgsPlpNOxKQ20qZsDi1kelczjLkmbVG0CkDRe0sSq9l8mTeBSvC9JIwbbhji29wALqndExNPA8uK9etvUzYHFLI99i5/fJT0eUNlmF/vHkYLOCkkrgXeR5kcqPgOsqmMbirHASwPsX1q8V2+bunny1iyPCaQv/udq9l8PEBE3A5tv5PNXsfFnmpo10GSqavbX06YuDixmeWwP/DEi5gzx8y+y4aesm7UU2GqA/WNYl6XU06ZuPhUya1Jxr8coqi4rD0ErT4UWUDNPImknUga1oIE2dXPGYta8A4qfdzXRRytPha4BTpG0ZUS8Wuw7irSQ1c0NtKmbA4tZ8/Yvfv5yqB1ExBJgSSOfkTSKdRPAOwCjJX2i+P3qiFhevJ4NnAz8t6SvA7uQ7k/5VtXl5Xra1D823yBnZZbjBrniVvopwLBo4xdK0njg9xt4e+eIeKqq7QRgFvBB0pzJZcCZEbGmkTZ1j82BxcrMCz21hidvzSw7BxYzy86Bxcyyc2Axs+wcWMwsOwcWM8vOgcXMsnNgMbPsHFjMLDsHFjPLzg8hmrWRpJtJ654IeAz47FAe8ut2flbISq3dzwpJGhMRLxevvwW8FhEz2nHsdvKpkFkbVQWVYaRFlPryX3YHFrM2k3Q1aYX+dwMXdng4LeHAYtZmEfEx4O2kFee+2OHhtIQDi1kGkqYUdZoPqtq3s6QXJF1S275YPOmHwLHtHGe7ePLWSi1zidUbgU0jYm9JY4DbSSu8HR4Ra4r6yCMj4vmi/deACRFxdLPH7ja+3GyWzxnALZI+AnyFtKr+0VVLO44FfippJOly86PASR0ZaYs5Y7FSy325WdI8YG/SmrF7RcSiHP32Gs+xmOX1OKnG0BllDSrgwGKWjaTjgc8C9wOf7/BwOsqBxSwDSQeQSmd8ATgB2Kv6ClFN2y8VV5B2aecY28lzLFZqmeoKTQR+DcyKiOnFvnnAmIjYs6btaGA+qVbz2RHxi6Eet5s5sFipNRtYJG0H3AncA0yrFCyT9CFSadJDImJuVfvzgDdI8zCvRsTZTf4RupIDi5VaOx9ClLQjcAswCTgCODQiprX6uJ3gORaz9jkXuCAilgEPkQJMX3LGYqXWroxF0mTgN8ALpCeaBWwPbBkRK1t57E5wxmLWHt8EjomIcRExPiLGkW73n9jhcbWEA4tZixWXnUdHxM9r3lpIn54O+VTISq3dK8iVhTMWM8vOgcXMsnNgMbPsHFjMLDsv9GTWRq4rZFYCrivUGj4VMmsj1xUys5ZwXSEzy851hcysLq4rtD5P3lqpua5Qa/hys1k+ritUcMZipea6Qq3hORazvFxXCAcWs2xcV2gdBxazDOqpKyTpCElLJM2X9Iik+yVN6MR4W82BxaxJRV2hnwEXRsSPIuIO4HrgrJqmk4FvR8TkiJgA3A0c09bBtokDi1kTirpC/wPMA6qf+TkbeL+kg6v2TQZ+W3xuK+B9wK1tGWib+aqQlVqb6wr9AVgGrAHeBcyIiG+0+rid4MBipdbG8h9jgccjYpvi978iVU/cqR+XTfCpkFl7TAYeqfwSEb8jPdk8plMDaiUHFrP2mAw8XPlF0jRgUUQ807ERtZBv6Tdrj78GDpT0AdKt/s8Ch3V2SK3jORYrNdcVag2fCplZdg4sZpadA4uZZec5FgNg/Glzjwc+2elxtNvKpx+cssmOE27VsOFrOz2WNvvPpy44+Hut6twZi1V8knRJtFTetu07HkDDyhZUJtPif0R8udmqzX/qgoOndnoQ1lrjT5t7U6uP4YzFzLJzYDGz7BxYzCw7BxYzy86Bxcyyc2Axs+wcWMwsOwcWM8vOgcXMsnNgMbPsHFjMLDsHFjPLzoHFzLJzYDGz7BxYzCw7BxYzy86Bxcyyc2Axs+wcWMwsOwcWM8vOgaWPSNpFkjo9DjMHlj4haRLwBPCgpI86wFgnObD0j7cBrwATgStwgLEOcmDpL5WylpvjAGMd5MDS3xxgrCMcWMrBAcbayoGlXBxgrC0cWMrpLQEmIgb5iFn9urIovKTxwGOkKx3WOpUAc+Ubix+LTbZ/9z2dHpD1h64MLMASYC7wjk4PpIeMA7YZwudWAcM0ctOlmcdjJdaVgSUiXgWO6PQ4eomk3YEbgDF1fmQVsAb4D+CskduOm9OqsVn5eI6lfFYBK4EfAe+KiM9HxDMdHpP1ma7MWKwl1stQHEyslfoiY5H0kKQ1knxp462coVjb9UVgAWYDkzo9iC7jgGId0xOBRdK9kmKA7QWAiJgVEQ9nPN4nJN0uaYmklZIWSpouaeQAbSdIukHScknPSZopaXijbTIZgQOKdYFemWMZR3py959r9i9s0fG2AX4FfAN4CdgTOBN4O3BipZGkscD1wCPA4cCuwEWkgD293jaZvASsJQUUz6FYR/VKYNkKuC8i/m87DhYR/16z61eSRgNfknRSrLtN9QRgM+DIiHgFmFe0O1PShcW+etrkGPMTksZGxJoc/Zk1o+tPhSRNII3zgQ4PZQlQeyp0EHBdTXC4nBRIpjTQJgsHFesWXR9YgIOLn3dJ2qRqq+vhOSUjBts28NnhkkZJ2hc4GfhurP9QzXuABdWfiYingeXFe/W2MesrvRBY9i1+fpc0KVnZZlcaSFogaXXxerWk6i/yZ0gTmoNtA3mt2G4FbgZOqXl/LGluo9bS4r1625j1lV6YY5lA+uJ/rmb/9ZUXEbGxf/mvAt4/xGPvDYwiTd5+DZgFfLGmzUD3zqhmfz1tzPpGLwSW7YE/RsRQn2V5EXh5KB+MiHuLl7dJ+jPwQ0kXRcQTxf6lpInlWmNYl6XU08asr3T1qVBxr8comrus3MypULVKkNm5at8CauZJJO1EWo5gQQNtzPpKt2csBxQ/72qij2ZOhartU/z8fdW+a4BTJG1ZPJENcBSwgjQnU28bs77S7YFl/+LnL4faQUQsIV0qrpuka0lzOA+THtzbB/gK8JOq0yBIE8gnA/8t6evALqQb6b5VdXm5njZmfaXbA8sexc9b2nzcu4HjgPHAauBJ4HSqrkQBRMRSSR8mTepeRZozuZgUOOpuY9ZvujqwRMTUDh13BjCjzraPAH/XbBuzftLVk7dm1pscWMwsOwcWM8vOgcXMsnNgMbPsHFjMLDsHFjPLzoHFzLJzYDGz7BxYzCy7vgksc6WHFktrcNEys47r6meFGvEEzL4ffvVVeKjTY+lRUwDGnzb3pg6Po61Wvfjs+BFjt3+qziWU+8VkYH4rD9A3GcvJEbO+mrFoGdIEpBuQliM9hzSTeoqMScchxQDbCVn6t6zWLHtxHLG2b74HdZoP/GcrD9A3GUtWeYqM/R1pMaeKJzP3n9VTFxxcqn+yK6RDAjgwIlYM2tjq5sAysDeLjFEUGaMoMoZ0IfUt0HQ3Ecta2L9Z1+r6FPAS6WSkuFR6c32Ur0v7/Ulae5U0v0WHPQi4jtYVGWt1/2Yd1fWB5eSIS+6GpfsUNX3OkXY6EuY9Bn+6t961bCUhjRh0W+ctRcZovMjYE0irkRYi/a+a93L0bz1E0k6SbpD0qKSHJV1Yb9G9XtT1gQXgdpgxGbacJZ22P9y3GuJGmHRGxJur618rLXiuKFr2nLT62uaKljVTZGwxafW5TwOHAncCs5H+T6b+rTetBk6NiPcCfwPsBRzZ2SG1Tk/MsXw54tI7pZn/BOe/Amt/AHvNiHi+us1Hq4qWbV9sVYayUv/QioxFXAdcV7XnGqRNgOlI/0bE2qb6t64l6WLgfw/w1tURcTDpHx0i4g1JDwA7tXF4bdUTgQXgz7B4c9j6pzDn9Ih7Gvx4o0XLchcZuwKYRlqc+8kW9G/dYTfgNt5aiveP1b9I2gb4e+Aj7RlW+/VEYPmBNOcYmPgorNwdPk5aQb8RnwG+X0e7yjnvW4qMkafIWCUbaVX/1lm7AT+OiDs21EApe70C+NeIeLRtI2uzrp9jmSWd+kn41I/he9fBiZNgi+orRNV+JP0EKb4u1V5ZqZwKDbZVXAMciLRl1b5miox9HPgz8IcW9W8dJmk7YDtSLaoNtRkO/AdwX0Rc1K6xdUREdO32DTjsJYifwa2VfXfAkgdgWW3bs2GHP8CqB+G178C5TR0bxgYsDpgXsH/A8QHLAs6paXdswOqAcVX7rgw4NeCggEMC5gREwEkN9++t5Rspi9wsQz/7F33tQzoTqGyqanMZKXNWs8fr9q1rM5bzpAnT4Mr7YNFD8KHK/jvhjN1g8+9IZ1a3fy/89E749e/goa3TjPvQRSwFPgwMJ2U7Z5GKjJ1R03JY0ab6suFC4LPAlcDPgAnAsUR8ewj9W++YVPy8jfWvNH4aQNI+wOdIRfjukzRf0smdGGg7qIikPe18aY9/hN/8EHb8SzhrVzj4byP6dsbd8lF6Gn5UNHlLv6TvA/sCx9S8tTDW1ewujZ6YvB3M7nD5jXD5jIjnL5Juegf8U6fHZKWzG3BvNH7Fsi917alQvS6SjtoPdt0fjnlGWn0U/PgdMHKmNKbTY7NykDSMdMr7YKfH0i16PrDsB9/5IVywY8SInSJG7BgxYhGs2gIO6fTYrDT+ivSclwNLoacDy6XSjM1gkxMiTq/e/xy8uJWLsFv77Fb8dGAp9MXkrdlQ5Zi8lXQW8BVgy/AXCnBgsZLLdVXI1tfTp0Jm1p0cWMwsOwcWM8uubwKLpIckrZHrCpl1XN8EFmA2657XMLMO6pvAEhGzImNdIUkTijVKl0t6TtJM1VH3R9JxkmKA7YSadkPq36wX9MWzQrmpxXWFMvVv1rUcWAb2Zt2fKOr+qKj7I+nCqLOuUNRRV6iJ/s26VtefCkk6uTiVmFG1bz9Ja9XiukLR4rpCLezfrKO6PrBExCWkxadPgVSfBZgH/Ik6V95XMmKwreojb6n7E0OoKyRptaSFqqOu0BD6N+taXR9YCjOALSWdBtxHWgJwUlTVFZK0QEVdoeIL3VV1heS6QqVXpqJlPTHHEhGXSpoJnA+sBfaKmrpCUVVXaABtqysUA9QVKlZmny7p38J1hcqsUrTsHkkjSZn3kaRlTPtKTwSWwmJga2DOEFbpcl0h67iIWExJipb1xKmQpDnARGAlqZRGoxo9FXpL3R+1uK5Qpv6tR1QVLbtukKY9qesDi6RTgU8B3wNOBLZQTV0hSRcUV4lWSHq9+HloVZMh1RVSi+sKZezfekgpipZ1uv7IxjbgMNK/8rdW7VtCTV0h4CbgpqrfFwC3NHHcsaSUdR6pXszxwDJq6v4Ax5LOm8dV7bsSOJV0SfkQYE7xZzip0f69teX/sVx1haYUfR1UtW9n4AXgkqp9w0lB5Vud/rO39L9rpwewkb+oCaTTk2dYv+jTicVf4JlV+xYDM4vX44BXgdMyHP9GUhaxGDgbGF7T5rhiLOOr9p1Hqi20vPjsb4FPD6V/b235/yxLYCn6uhG4vXg9hlQV8X+q/14pSdGyvlhBrrjMvJr0P8mmwNyI8GLaNqicK8hJ2g+4BTiQtFTlXwD7RnEHdlG07DbgIWBN8bH/F+lerb7S84FF0njgyYgYVvy+P+kUY8eIeLaTY7Pul3tpSknzgL1JV/f2iohFOfrtNV0/eVuHI0inPgBExPXFyx06MxwruceBUcAZZQ0q0B+BZQrw5l+gpIuB1yPirs4NycpI0vGkut33A5/v8HA6qh8Cy/uA9xSXmF8Djiad45q1jaQDgFnAF0hPr+8l6aANtP1S8WDtLu0cYzv1/ByLWTMy1RWaCPwamBUR04t984AxEbFnTdvRwHzS3eBnR8QvhnrcbubAYqXWbGCRtB3pQdN7gGlRfKEkfYh0s+MhETG3qv15wBukeZhXI+LsJv8IXcmBxUqtnQXLJO1Iuhw9iXTR4dCImNbq43ZCP8yxmPWKc4ELivtaHqKPF393xmKl1q6MRdJk4DekW/yDtETG9qR6zytbeexOcMZi1h7fBI6JiHERMT4ixgG/Jz2133ccWMxarLjsPDoifl7z1kL69HTIp0JWau2cvC0TZyxmlp0Di5ll58BiZtk5sJhZdg4sZm0iaU9JLxbrMq+UdGe/1hVyYDFrn9eBL0fEJsA2wLuBCzs7pNZwYDFrk4i4PyLmFK9fI1VteGczfUqaUFRXXC7pOUkzJQ3fSPubiiUbBto+2Ezf1XqpYJlZ35D0TtJdt19too+xwPXAI8DhwK7ARaSEYfoGPvZFYHTNvpnA3wB3N9n3mxxYzNqsqCd1N3B19ZIKQ3ACsBlwZES8Aswr1ns5U9KFxb71RMQjNWMZCewB/CQiVjfTdzWfCpllIOnk4nRiRtW+/YpCevOr9r0NeJS0APxhTR72IOC6mi/55aSAMKXOPj5KqnP1Xzn7dmAxy6Ao4bEUOAXeLJk7D/gT61fZfJBUS2qP2j6UjBhsq/rIe6gpyRsRT5NqWq1XwncjjgaeBW6t2d9U3w4sZvnMALaUdBpwH2l5hEkRsQpA0gmkK0E7AcuLdZqvqPp8ozXGx5LKjNRaWry3UZJGAYeSToNqHxpsqm/PsZhlEhGXSpoJnA+sJdUVer7q/dnA7I10Uakx3tBhB9inDeyvdSiwBW89DWq6bwcWs7wWA1sDcyLingY/+yLwcgPtlwJbDbB/DANnG7WOBh7fwDib6tunQmaZSJpDuoS8Evj4ELpo9FRoATXzHcXczubUzI8MMNYxpAnaDWUrQ+4bHFjMspB0KvAp4HvAicAW1VeIijYXFFeJVhS39a+QdGhVk8qp0GBbxTXAgcXl64qjSJPDNw8y5COATdhwYGmmbzpeld6bt05upPmCzZrs47Cin1ur9i0BltW0uwm4qer3BcAtTRx3LOnUax6wP3A8sAw4p6bdscBqYFzVvmuB+c32vaHNGYtZEyRNAK4klfn9UNVbZwCbSzqzat+7SeU/kDSOVF/86qEeOyKWAh8GhpOynbOAi4tjVxtWtFFx7G2Lz12eoe8BeWlKK7U21xVaTcocAtgUmBsRh7T6uJ3gjMWsDSSNB4ZFxKYRsRlwAHCwpB06O7LWcGAxa48jgFcrv0TE9cVLBxYzG7IppHkYACRdDLweEXd1bkit4xvkzNrjfcDOklaQ7sp9BTiws0NqHU/eWqm5rlBr+FTIzLJzYDGz7BxYzCw7T95aX5M0iXSpd2NmSHp9A+9dFRH3Zh5W33NgsX73CdKC1RtbXf70DexfS1p42oGlQb4qZH1N0jbA08CoIXx8BbBzVC3WZPXxHIv1tYhYAswirZHSiNeByxxUhsYZi/W9IWYtzlaa4IzF+t4QshZnK01yxmKl0GDW4mylSc5YrBQayFqcrWTgjMVKo86sxdlKBs5YrDTqyFqcrWTijMVKZZCsxdlKJs5YrFQ2krU4W8nIGYuVzgayFmcrGTljsdIZIGtxtpKZMxYrpZqsxdlKZs5YrJSqspbA2Up2zlistCRtDfwA+IIDS14OLGaWnU+FzCw7BxYzy86Bxcyyc2Axs+wcWMwsOwcWM8vOgcXMsnNgMbPsHFjMLDsHFjPLzoHFzLJzYDGz7BxYzCw7BxYzy86Bxcyyc2Axs+wcWMwsOwcWM8vu/wO1Fxwsw4CthQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 360x360 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot_flowsheet()\n",
    "XA1=1-XB1\n",
    "plt.text(xs-la-hla,ys+ly/2-14, r\"$x^1_A$\"+'={:.2f}'.format(XA1), color='r',size=15)\n",
    "plt.text(0,95, r\"$x^1_A+x^1_B=1$\", color='r',size=15)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "id": "bb48b385",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAARcAAAEvCAYAAABv+oUUAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAf5ElEQVR4nO3debhcVZnv8e8vCQQCJAR4HICYILRoAjF9QaICggooMtig4tiIQ3NpROyrjZArgYR5aBwgIHq5VyAOyGBrRyYDiIBKgGYMEEYZApEhBCEkwQzv/WPtIjuVM1Sd2ruqTtXv8zz7OXV2rVp7Bc55z7vX3nu9igjMzIo2pNUDMLPO5OBiZqVwcDGzUji4mFkpHFzMrBQOLmZWisaCizQK6TakxUjbFTSm9jummdWt0cxlCbAPcHkBY2nnY5pZnRoLLhHLiXihoLG07zHNrG6ecwGQjkC6A+l1pAtbPRyzTjCs6UeUpgEQMa3px+7ds8BJwEeA9Vs8FrOOUFvmIp2B9J+5789Euh5pnVJG1ezjRfyKiF8DC0vp36wL1Zq5nA48hjQJeC/wUWAXIpYjXQVMArZF+hERFxYwrt6PB/R7TOm3wC699H0LEfsWMEYz60NtwSViIdL3gYuBUaRf9L9l732s8FH1dbxajungYdZy9Uzo3gVsD0wh4um6jiL9FullpJeBY4Bj3vg+ZRnFHs/MWk41recibQ9cC/wO2JKIParenwDcC4wlYn4/fU0D+p7Q7et40mXA5tl3rwMfI2JZ1eevBnbtpfebidi7l+OelB3vkD7/DWbWr/5Pi6QtgFnAYcB1wF+Qdifixlyro4GZwLuAvoNL48d7B7ADESuQfg1MBG5bo4/egkfvxxxG+m8xFBiKtB6wgogVA/53mHW5vk+LpJHAVcB3ifgvIpYAZwIn59pMBBaQMo13NTSa/o4nrQsMyQLLhsBbgYcbOmZyLLCUdMr2hez1sQX0a9a1ajst6rMH/QT4FjAGOIyIfy1gXL0d692kbOZ+YFvgaCIuLu14ZjZgjQWXdKn4UuBxKqcUER8qZGQ9H+/zwLZEHIc0CriViMayJTMrRaN36B4DTCZikcSYydw6d454EFgFXAkcHUGRK4BvB9yevR4N/K2nRhJjgAtJE79ljcXM+jDwZ4ukHYClRCzK9qw4lSmPBdoF+EdgMnBg40NcwwTgO0jXA/+HNOnbkxWkYPKuEsdiZn1ofM6lt47FOcBjEXy/lAO0wVgktgGOIt1FvB1wcwS7F3kMs8GqlAcXJTYF/gnYq4z+61HyWCYAHwNuBdYtoX+zQavwJRckhpMWcvp+BA/W+JndpeLnQwYyljrNimBMBJ8iXcEys0xNwUViN4mQ2Du3byuJ5yXOzu0bCvwMuCuCs4ofbnuNJYJVZfRr1glqnnORuAFYL4L3S4wC/gT8Bfh4BCuzNheQLkl/uZ4rMxK7A7+PQM0ei4Sydn2KoM+7dSUuBzbznItZUs9p0fHA+yT2It3bshz4TO6XeWfgK8COwF0Sd0sc2VNHEpIYVtnIfrnz+7L9pY8F+GL2+f42M6tDXVeLJGYD7wdeBiZHDOw5IolDgJ/0166vTKbAsWwKbFXDWO7opx9nLmY59V4tehTYA/jGQH+ZM7OA9+S+3wE4v2pfs8byEr3cjGdmA1fPnMuhwDnAg8CyCN5b2CDqn3MpbCxFZFFZP85czHJqylwk9gRmAF8lPYX8Z4m9I7i6h7Zfy9puHcHjRQ621rFIHABcADxNuv9kOfDZCB7oocvqLMrMCtBv5iIxAfgjMCMiLUOQzXeMimCnqrYjgbtJpxonRvCbmgZRY+ZS61gkpgMRwbTs+wuABRFMrWU8tZIYQbqJDtKT4SNJk80AV0WwpMjjmQ0mfQYXiTcBc4A7gIMql3QlPgD8Adg3gitz7U8B/g6MAF6N4MSaBlFDcKlnLBK/AS6IYJbExsA1wHER/K6W8dRKYhzpEnhPtorgiSKPZzaYFPZskcSWwE2kleEOAPaL4KBCOq9/LE8Ci4GVpJXrpkZwZivGYtatirz9/2TgtAgWA3NJQabpJEYDG0YwIYKJpEW+j81O2czeIGlrSfMleU2gEhQSXCQmAQcB35F4Avg1sLXEekX0X6dJsHriNoJHgCCVKDHLO4m0VOqprR5IJyoqc/kP0tWYsRGMi2AsaS5iQkH912MSuYcIJQ4C5kfg8iT2Bklbk56WHwLs5eyleA0Hl+wBwpER/LrqrYdozanRu4GPZ7f83w58Dti/BeOw9nYSq2/FWBdnL4UrbbEos3aVZS1zYY3T9qXADhFRxtIcXanw9VzMBoF81lLh7KVgzlysq/SStVQ4eymQMxfrNj1lLRXOXgrkzMW6Rj9ZS4Wzl4I4c7Fu0lfWUuHspSDOXKwr1Ji1VDh7KYAzF+sWtWQtFc5eCuDMxTpenVlLhbOXBjlzsW5QT9ZS4eylQc5crKMNMGupcPbSgFLKuZq1kUNJpWt6W4R9VB/vrQ8cDny9hHF1PGcu1tEkvZm+10ieBXwKWNbL+3dGxLOFD6wLOLhYV5MUwIiIWNrqsXQaT+iaWSkcXMysFA4uZk0i6WVJS7NtvqQtWj2mMjm4mDXP9hGxfkSsDzwH/KLVAyqTg4tZQSTdKSl62J4HiIins3ZDSZe5G7qaImkbST+SdI+klZJu7KPteEnXS1oi6VlJJ2TjqKtNPXyfi1lxxgKvAN+s2v9Q5UUWaDbL2n2wweNNIFX8vJV0R3GPJI0GriNVxfg4sDVwFim5OLbWNvXypWjrakVeipa0ErgrInbsp906wC3AixGxTwPHGxIRq7LXlwObRcTuPbSbAnwbGBsRr2T7vg1MA94SEa/U0qbe8fm0yKwAksaTfp/u7a9tRCwnFRH8cCPHrASWGuwNXFsVIC4hnZrtVkebuji4mBWjkoHcJml4bhOApHGS8nW8vkGa1CV7X5KG9bcNcGzvBObld0TEU8CS7L1a29TFwcWsGLtkX39IepSgsp2f7R9LCjxLJS0j1TD/WO7zXwSW17ANxGjg5R72L8req7VNXTyha1aM8aRf/q9U7b8OICL+AGzQx+dn0fczUI3qaXJVVftraVMzBxezYmwO/DUiZg7w8y/R+9PZjVoEbNzD/lGszlZqaVMXnxaZNSi7F2QEuUvOA1DmadE8quZNJI0hZVLz6mhTF2cuZo3bM/t6WwN9lHladDVwlKSNIuLVbN+nSYth/aGONnVxcDFr3B7Z198NtIOIWAgsrOczkkawelJ4C2CkpE9m318VEUuy1+cDRwK/knQ68HbS/SvfzV16rqVNXXwTnXW1Im6iy2673w0YEk38hZI0DvhLL29vFRFP5NqOB2YA7yPNoVwATIuIlfW0qWt8Di7WzbxYVHk8oWtmpXBwMbNSOLiYWSkcXMysFA4uZlYKBxczK4WDi5mVwsHFzErh4GJmpfCzRWZNJOkPpKUNBDwMfHmgz+60O9/+b12t2bf/SxoVEX/LXn8XeC0ipjbj2M3m0yKzJsoFliGktVI69q+7g4tZk0m6irQ497bAGS0eTmkcXMyaLCI+BryFtLjU4S0eTmkcXMwKIGm3rHTr3rl9W0l6XtLZ1e2zNVIuAg5u5jibyRO61tUKrrh4A7BeRLxf0ijgT6TFnD4eESuzkqnrRsRzWfvjgPER8ZlGj92OfCnarDjHAzdJ2gv4FmlB7c/kVnIbDVwqaV3SpegHga+3ZKRN4MzFulrRl6IlzQbeT1omcnJEzC+i38HIcy5mxXqUVGbk+G4OLODgYlYYSYcCXwbuAb7a4uG0nIOLWQEk7UlaOf9fgMOAyfkrR1Vtv5ZdWXp7M8fYbJ5zsa5WUGmRCcAfgRkRcWy2bzYwKiJ2qmo7EribVL71xIj4zUCP2+4cXKyrNRpcJL0JmAPcARxUqVsk6QOkSoX7RsSVufanAH8nzcu8GhEnNvhPaFsOLtbVmvngoqQtgZuAicABwH4RcVDZx20Vz7mYNc/JwGkRsRiYSwoyHcuZi3W1ZmUukiYBfwaeJz0JLWBzYKOIWFbmsVvFmYtZc/wH8NmIGBsR4yJiLOnRgAktHldpHFzMSpZdkh4ZEb+ueushOvjUyKdF1tVciL48zlzMrBQOLmZWCgcXMyuFg4uZlcKLRZk1kesWmXUJ1y0qj0+LzJrIdYvMrDSuW2RmpXDdIjOrmesWrc0TutbVXLeoPL4UbVYc1y3KceZiXc11i8rjORezYrluUcbBxawgrlu0JgcXswLUUrdI0gGSFkq6W9IDku6RNL4V420GBxezBmV1iy4DzoiIiyPiVuA6YHpV00nAORExKSLGA7cDn23qYJvIwcWsAVndot8Cs4H8M0InAu+RtE9u3yTgv7PPbQxsB9zclIG2gK8WWVdrct2iJ4HFwErgHcDUiDiz7OO2ioOLdbUmlhYZDTwaEZtm3/8DqUrjmE5dcsGnRWbNMQl4oPJNRDxCeiJ6VKsGVDYHF7PmmATcX/lG0kHA/Ih4umUjKplv/zdrjncDH5H0XtJjAc8A+7d2SOXynIt1NdctKo9Pi8ysFA4uZlYKBxczK4WDi5mVwhO6BsC4Y648FPhcq8fRbMueum+34VuOv1lDhq5q9Via7OdPnLbPj8s8gDMXq/gc6V6MrrLOZm+7Fw3ptsAyiSb8IfF9LpZ39xOn7bN7qwdh5Rp3zJU3NuM4zlzMrBQOLmZWCgcXMyuFg4uZlcLBxcxK4eBiZqVwcDGzUji4mFkpHFzMrBQOLmZWCgcXMyuFg4uZlcLBxcxK4eBiZqVwcDGzUji4mFkpHFzMrBQOLmZWCgcXMyuFg4uZlcLBpYNIersktXocZuDg0jEkTQQeA+6T9FEHGWs1B5fOsQ7wCjABuBwHGWsxB5fOUimfuQEOMtZiDi6dzUHGWsbBpTs4yFjTObh0FwcZaxoHl+60VpCJiH4+YlaftixEL2kc8DDpCoiVpxJkrvj7godj+Obb3tHqAVnnaMvgAiwErgTe1uqBDCJjgU0H8LnlwBCtu96igsdjXa4tg0tEvAoc0OpxDCaSdgCuB0bV+JHlwErgZ8D0dTcbO7OssVl38pxL91kOLAMuBt4REV+NiKdbPCbrQG2ZuVgp1shUHFCsbB2RuUiaK2mlJF/yWJszFWuJjgguwPnAxFYPos04qFhLDYrgIulOSdHD9jxARMyIiPsLPN4nJf1J0kJJyyQ9JOlYSev20Ha8pOslLZH0rKQTJA2tt01BhuGgYm1isMy5jCU98fvNqv0PlXS8TYHfA2cCLwM7AdOAtwBHVBpJGg1cBzwAfBzYGjiLFLSPrbVNQV4GVpGCiudUrOUGS3DZGLgrIv5vMw4WET+q2vV7SSOBr0n6eqy+nfUwYH3gwIh4BZidtZsm6YxsXy1tihjzY5JGR8TKIvoza1TbnxZJGk8a570tHspCoPq0aG/g2qoAcQkpmOxWR5tCOLBYO2n74ALsk329TdLw3FbTA3dKhvW39fLZoZJGSNoFOBL4Yaz5EM47gXn5z0TEU8CS7L1a25h1nMEQXHbJvv6QNFFZ2c6vNJA0T9KK7PUKSflf5i+SJjn723ryWrbdDPwBOKrq/dGkuY5qi7L3am1j1nEGw5zLeNIv/1eq9l9XeRERfWUAs4D3DPDY7wdGkCZ0jwNmAIdXtenp3hpV7a+ljVlHGQzBZXPgrxEx0GdfXgL+NpAPRsSd2ctbJL0IXCTprIh4LNu/iDTZXG0Uq7OVWtqYdZy2Pi3K7gUZQWOXnBs5LcqrBJqtcvvmUTVvImkMaSmDeXW0Mes47Z657Jl9va2BPho5LcrbOfv6l9y+q4GjJG2UPckN8GlgKWmOptY2Zh2n3YPLHtnX3w20g4hYSLqMXDNJ15DmdO4nPey3M/At4Je5UyJIk8pHAr+SdDrwdtLNdt/NXXqupY1Zx2n34LJj9vWmJh/3duAQYBywAngcmELuChVARCyS9GHSRO8s0hzK90jBo+Y2Zp2orYNLROzeouNOBabW2PYB4EONtjHrNG09oWtmg5eDi5mVwsHFzErh4GJmpXBwMbNSOLiYWSkcXMysFA4uZlYKBxczK0XHBJcrpbkLpJW4dpFZW+iY4PIYnP8T1y4yaxtt/WxRPY6MmAFAbUvr9i8tDH4O8D7Sw4YXANPpbxFs6RDgJz28869EnJ9rN7D+y7MbwLhjrryxRcdvieUvPTNu2OjNn6hxSeZOMQm4u+yDdExwKVQxtYY+RFqzpeLxgvu3Aqxc/NLYYRu/5Sk0dFWrx9JEdwM/L/sgWnMx+w4gBRGN/RmSpgDfBsZSWXNF+jaVwmh9rcOyOnPZiIjFhfdvhcrqi4+IiKX9Nra6tP2cy9nSkUhxrvTGEginS7u+IK2aJd1d0mH3Bq6lvFpDZfdv1nJtH1yOjDj7dli0c1bW4yRpzIEw+2F44c5al6+UhDSs3221tWoNUX+toceQViA9hPQ/q94ron+zttb2wQXgTzB1Emw0QzpmD7hrBcQNMPH4iDcW1r5GmvdsVrvoWWnFNY3VLmqk1tAC0kJT/wzsB8wBzkf6XwX1b4OQpDGSrpf0oKT7JZ1Ra2G/wWpQTOh+I+LcOdIJX4JTX4FVF8LkqRHP5dt8NFe7aPNsyxnIIt0DqzUUcS1wbW7P1UjDgWORfkBEZeLQtYy6ywrg6Ii4Q9K6wGzgQOCK1g6rPIMiuAC8CAs2gE0uhZlTIu6o8+P11i4qutbQ5cBBpDV5Hy+hf2sDkr4H/FsPb10VEfuQsloi4u+S7gXGNHF4TTcogsuF0szPwoQHYdkO8AnS4tn1+CI933tSrZKmrlVriGJqDVWykrL6t9baHriFtcv+/jX/jaRNgX8C9mrOsFqj7edcZkhHfw6+8FP48bVwxETYMH/lKO9i6ZdIcbpUfcWlclrU31ZxNfARpI1y+xqpNfQJ4EXgyZL6t/awPXBbRNxatT1RaaB0inw58P2IeLBVA22KiGjb7UzY/2WIy+Dmyr5bYeG9sLi67YmwxZOw/D547Tw4uaFjw+iABQGzA/YIODRgccBJVe0ODlgRMDa374qAowP2Dtg3YGZABHy97v69lb6Rssn1C+jnTVlfX+6jzVBSYPluq//dzdjaNnM5RRp/EFxxF8yfCx+o7J8Dx28PG5wnTcu3fxdcOgf++AjM3QQmN3TwiEXAh0k/DLOA6aRaQ8dXtRyStcnP+j8EfJk0UXcZMB44mIhzBtC/DR6V59oekjQst+V/Nn4EvEoqsNfxOuIO3VOlHT8Pf74ItnwrTN8a9vlgREdPllkxirpDV9I3SY9wVPtiRFwsaWfSfMxcUhVPgP8XEWc3ctx2NigmdPuzA1xyA1wyNeK5s6Qb3wZfavWYrOtsDzwKfLZq/0MAEfFH1sxwO96gz1zOkj59OFyyEFYG6f/em2HoqbDxcRH1XH62LlRg5nIH8FhEfLqYkQ1+bTvnUqtd4byL4LQtI4aNiRi2ZcSw+bB8Q9i31WOz7iBpCGlu7b5Wj6WdDOrgcq40dX0YfljElPz+Z+GljV2b2ZrnH0gPnTq45Az60yKzRhRxWiTpk6Qrg1tHxOP9te8WgzpzMWsT2wOvAX9p9UDaiTMX62peLKo8zlzMrBQOLmZWCgcXMytFxwQXSXMlrZSLopm1hY4JLsD5uCiaWdvomOASETMi4v6i+pM0PlvzdImkZyWdIGloDZ87RFL0sB1WRP9mg0VHPLhYNJVcFK2g/s3amoNLzw4j3c59YKTaQrMljQSmSTojaitadnv0VhStmP7N2lrbnxZJOjI7rZia27erpFUquShalFwUrcT+zVqu7YNLtpjOIrJFj5UWsp4NvECN5UKUDOtvy31kraJlMYCiaJJWSHpINRRFG0D/NshIp+0k3f6S9Ojr0iPLpP+cI03r2DVe2j64ZKYCG0k6BriLtFbpxMgVRZM0T1lRtOyXuq2KoslF0Yx4HW7+RsQ2w+EXm8LbtoW3ntHqUZVlUMy5RMS5kk4ATgVWAZOjqiha5Iqi9aBpRdGih6Jo2Yrvx0r6QbgoWteKmHIPcE96fdxr0qwnYYNtGulTYjxwDvA+0h+sC4DpEW8spdnb5z4JfBPYllTS5klgJnBGBH9vtH8YJMElswDYBJgZLopmg5x08jbwpQnws/898D5o5KrjpsDvgTNJP3M7AdOAtwBHFNB/e5cWqWykiBqkS7uvDuDzh2Sf73PLtb8J+EVVH2OydvsN4Pifyj67VRn9e2voZ6uQ0iL1HXP6RnD7IvjpfzXWT0yBWAQxMrfv2xBL8vvq6O9kiJchVET/bT/nIulo4AvAj0kRdUNVFUWTdFp29WippNezr/vlmgyoKJpKLopWYP/WYtLZR0qEdG7uqubpu0ovrJJm3b163/R1YLcH4ZnHIz6/f4OHza46UtRVx4XAuoX13+q/HP38Vdmf9Jfl5ty+hVQVRQNuBG7MfT8PuKmB444mnYbNBvYADgUWU1W0DDiYVGB8bG7fFcDR2f+YfVmddX293v69NeVnrLDMBW57Ce56Jb0+cQw8vAxueQ6mrbO6zdXz4JpH4Hit/fkQxLD+tlz75yGm9dDPaxBH1TbmGAoxAmIXiCcgziqq/5b/z+3jf/p40hWcp8kWtcr2H5H9QEzL7VsAnJC9HksqPHVMAce/gZRNLABOBIZWtTkkG8u43L5TSOUklmSf/W/gnwfSv7em/JwVGFy+/7VUXPOcY+DPL8IDS+GEN+fePyy9P29Zeu+BpXDx5avfj0PS+31vufbLIf5t7XHEfIhTahtzLMv1fRHEkKL674iV6LJL0CtIPyjrAVdGhFf/t34VvRKdNGchbLcJvLIKLpwcMaXmiw8SmwJb9dcugjuy9suBf4/gB1X9PANcGMF3ajjm/wBGkCZ0jwN+HsHhRfQ/mK4W9UjSOGBIRKyXfb8H6Xb6LSLimZYOzrrQiwtgg03g0pn1BJZM069qRnBn9vIWiReBiyTOiuCxRvtv+wndGhxAOg0CICKuy15u0ZrhWLeSLpwJe0yAB5fBDp8YQBf13uw5j6o7uiXGkO5bWeMO8BpVAk0le2qo/04ILrsB8yvfSPoe8HpE3Na6IVm3kWYcDZ/7Avz0x3DtETBxw/yVozXbXvzLdGXp9OorLgO8qklRVx13zr5Wqhg01n+rJ9QKmJB7lFTYeympvMMCYLdWj8vb4NgoYEIXztwfXg64LHdV89aFcO/itdueuAU8uRzuew3OO7mx48ZoiAUQsyH2gDgUYjFE1VXNOBhiBcTY3L5rIP4dYm+IvSCmZ5+9pN7+e9sGfeYSEdtExNCIWD8iNoiIt0aE7xWxppBOGQ8HXQF3zYe5H1j9zpzjYfsNpPOmrfmJd10Kc/4Ij8yFTSY3cuwIFgEfBoaSsp7pwPeA46uaDsna5B+SvJ10tfMy4FLSc3BTSM/E1dt/jzriapHZQDWzbpF06o7w+T/DRVvCW6fD1vtEfHBM2cdtlUGfuZgNHjtcAjdcEjH1OfjbjfC2N7V6RGVy5mJdrVmZi3TWp+HwS2DhyjTNI+DNQ+HUjSOOq+fy86DhzMWsKXY9Dy46LWLLYRFjhkVsOQzmL4cNO/ZmTwcXs5KlS9LrD484bMqa7zz7Emz8odaMqnw+LbKu5kL05XHmYmalcHAxs1I4uJhZKRxczKwUDi5mTSJpJ0kvZUuxLpM0R5LrFplZw14HvhERw0mr728LdGzdIgcXsyaJiHsiYmb2+jXSgu0N1i3SeEnXS1oi6VlJJ0ga2kf7G7PyyD1t72uk72qDfiU6s8FI0jbABKCBukUaSF2hw4GRVftOAP6R9KR0I32vwcHFrMmykjK3A1dFxJUNdHUYqczHgRHxCml515HANElnZPvWEBEPVI1lXWBH4JcRsaKRvqv5tMisAJKOzE4tcnWLtGtWT+vu3L51gAeBxyOioLpF0Ujdoo+SSt38oui+HVzMChARZ5MWtD4KQNIYUl2qF1hzacr7SKsm7ljdh5Jh/W25j7yTqrVsI+IpUlmbvmqn530GeAa4uWp/w307uJgVZyqwkaRjgLtIaytMjIjlAJIOI10hGgMsySqDXp77fL0LdI+m51X4F2Xv9UnSCNIKdL+MtR8ybKhv8JyLWWEi4lxJJwCnAquAyRHxXO7984Hz++iiskB3XYftYZ962V9tP2BD1j4lKqJvBxezgi0ANgFmRkS71y36DPBoL+NsuCaST4vMCiJpJuny8jKgRXWLVFNdIUmjSJO2vWUtA+67wsHFrACSjga+APyYVM98w/yVo6zNadnVo6XZIwBLJe2XazLAukUaSF2hA4Dh9B5cGuk7aXXdGG/eWrlRSN0i9s/6ydUtYiGwuKrdjcCNue/nATc1cNzRpNOw2cAewKHAYqCqbhEHk2qpj83tuwa4u9G++9qcuZg1QNJ44ApS1c9c3SKOBzaQNC23b1vgpuxzY0klh68a6LEjYkB1iyRtln3ukgL67pWXubSu1ty6RVpByiACWA+4MiK8QLeZDZykccCQiFgvItYH9gT2kbRFa0dWHgcXs+Y4AHi18k1EXJe9dHAxs4bsRpqXAUDS94DXI+K21g2pXL6Jzqw5tgO2krSUdPfuK8BHWjukcnlC17qa6xaVx6dFZlYKBxczK4WDi5mVwsHFzErhq0XW0SRNJN1j0pepkl7v5b1ZEXFnwcPqCg4u1uk+SVphv6+SGFN62b+KtFK+g8sA+FK0dTRJmwJPASMG8PGlwFaRW03Oauc5F+toEbEQmEFawKkerwMXOLAMnDMX63gDzF6ctTTImYt1vAFkL85aCuDMxbpCndmLs5YCOHOxrlBH9uKspSDOXKxr1Ji9OGspiDMX6xo1ZC/OWgrkzMW6Sj/Zi7OWAjlzsa7SR/birKVgzlys6/SSvThrKZgzF+s6PWQvzlpK4MzFulJV9uKspQTOXKwr5bKXwFlLKZy5WNeStAlwIfAvDi7Fc3Axs1L4tMjMSuHgYmalcHAxs1I4uJhZKRxczKwUDi5mVgoHFzMrhYOLmZXCwcXMSuHgYmalcHAxs1I4uJhZKRxczKwUDi5mVgoHFzMrhYOLmZXCwcXMSuHgYmal+P80iT6m/CCKDQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 360x360 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot_flowsheet()\n",
    "XA1=1-XB1\n",
    "plt.text(xs-la-hla,ys+ly/2-14, r\"$x^1_A$\"+'={:.2f}'.format(XA1), color='r',size=15)\n",
    "plt.text(0,95, r\"$x^1_A+x^1_B=1$\", color='r',size=12)\n",
    "XA2=1-XB2\n",
    "plt.text(xs+lx/2+5,ys-14, r\"$x^2_A$\"+'={:.2f}'.format(XA2), color='b',size=15)\n",
    "plt.text(0,90, r\"$x^2_A+x^2_B=1$\", color='b',size=15)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "id": "25c0e51a",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAARcAAAEvCAYAAABv+oUUAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAjMklEQVR4nO3deZgcVb3/8fcnCQECSQx4lT1BEDQsRlFQWQVEkEVBL6goIipyEfG6IMRLmDSLbILKJvB4fwaigizKlU1kMQIqm2wGCCjIHllCWEIWsnx/f1Q1qXRmerqnq7p7uj+v56lnpqtPn3MGMt85dU7V+SoiMDPL25BWd8DMOpODi5kVwsHFzArh4GJmhXBwMbNCOLiYWSEaCy7SaKQ7kOYgbZpTn9qvTTOrW6Mjl7nA7sBlOfSlnds0szo1FlwiFhLxQk59ad82zaxunnMBkA5DugtpAdKUVnfHrBMMa3qL0mQAIiY3ve2+PQscD3wMWLnFfTHrCLWNXKRTkH6beX0q0o1IKxTSq2a3F/EbIq4AZhVSv1kXqnXkcjLwKNIE4IPArsA2RCxEugaYAGyMdB4RU3LoV9/tAf22KV0FbNNH3bcSsUcOfTSzKmoLLhGzkH4MXAiMJvlFfyV97+O596pae7W06eBh1nL1TOjeA2wGTCTiqbpaka5CehnpZeAo4Kg3XyejjHzbM7OWU037uUibAdcBfwDWIWLnivc3Ae4HxhLxdD91TQaqT+hWa0+6FFgrfbUA+DgR8ys+fy2wbR+130LEbn20e3za3oFVfwYz61f/l0XS2sCVwCHADcC/kHYgYlqm1JHAVODdQPXg0nh7GwFbELEI6Qpgc+COZeroK3j03eYwkv8WQ4GhSCsBi4hYNOCfw6zLVb8skkYB1wCnE/E7IuYCpwInZMpsDswkGWm8u6He9NeeNBwYkgaWVYE1gUcaajNxNDCP5JLt8+n3R+dQr1nXqu2yqGoN+jnwHWBd4BAi/iuHfvXV1ntIRjMPABsDRxJxYWHtmdmANRZckqXiS4DHKF9SROyYS896b29/YGMijkEaDdxGRGOjJTMrRKN36B4FbEXEbACkO6oXb9imwJ3p92OAV6qUNbMWGnhwkbYA5r0ZWBLzkVYnoqg7XTcBdkH6OrCEZNLXzNpQQ5dFEh8CTgfeIHk+54AIFubUt7Zp08zq1+hT0U8AO0awPcm8yyca71JbtmlmdWpoziWCZzMvF5FcqhSqFW2aWf1y2c9FYn1gN6CvW/lzl2ebEodJ3CWxQGJKw50zs8b3c5EYBVwAfCGCN2ooPxkgIvnajDZr4P1czHJW08hF4hSJ32Zenypxo8QKwEXA5AgezqtT1dqTGJZ3mxH8JoIr8H4uZrmp9bLoZOAjEhMkDiHZX2Uf4DPAVsAxEtMk9supX722l64Kfba/NiWukni5j6Npl25m3aymy6IIZkn8mMz+KhG8QvKw4tS8O1WlPSL6bzMC7+di1mID2s8lgrr2V8mOJEj3c6lhJDHg9sys9Wq6iU5imf1VIti54v0393OJqL7lQi0TutXak1huP5cI5ld8vup+LhH0uiWDxPFpewdW+xnMrH/9XhZJLLe/isQOEUzLFMttP5ca2tsI2CKCRRJX0Mt+Ln0FjyptLrOfi8RKwKIIvJ+L2QBVvSxKl3yvAU6P4HcRLLefi0Ru+7n0157EcGBIGli8n4tZG2t4PxeJZfZziaCw/VwkltvPJQLv52LWhhp9cHECFfu5RFDYfi4S+wMbR3CMxGjgtogGd7+zriVpA+BPwEcj4qFW96fTNHr7/1HAVhHsGsFHgVVz6FM1mwL3pt97Pxdr1PEkl9YntrojnWjAIxeJLVjl39/iiDU3AsYDH2RynAPsHVHMna7a6Opb2Ol/3ksMWcjrb1vAFT/fI15b864i2nqzzZLeDvwWWAgsBvaPnphZZJtWvHTUMh1YiWSObQuPXvI14JFLBH/jiDW/BOwOXJae266owALA/nvsxxr3/Uece/cYNrzuf/nuWuMKa2upF4Ftoie2J7mp78tNaNOKdzxLV0uH49FL7hrbcqEnFgIvqKScutNve83f4qEnFmdejiSZTLZBLB21fJKl//6HArtIerdHL/mpKbiopDHAS8DrmdOvAOtET6PpA+pvTyWVt1s4obfPF9DeBOA84C3ALnm0aS2VHbWUlUcvn2x6bzpUrSOXCcBL0ROrF9iXmtpTSUu3W+iJqtstqJRkd4yeKtkd+2kv/fy9wFYqaV9gIt67d9DqZdRS5tFLzmqdc5kAPFhgP2pqTyUt3W6hJ3Lb4qGv9tI2V8y8fAWYm2O71ny9jVrKPPeSo1pHLu+l71++a0h+OTdWSedFT0ypeP8qYJs+6r01eqK3J5j7am/pdgslHQP8NHri1zX9BNX1+fMB71NJJ5OsFM0HDsqhPWuBKqOWMo9eclTbg4sl3Q9sCG/u+vZU9MRmhXUqp/ZqvSxq9s9nrSHpIuDTVP+juhi4KiI+2ZROdbB+g0t6WTAH+FD0RKH3lOTRXsVIaaX0a/mp6eVGSs3++aw1Ku5r6Y/ve8lBLZdFm5Is+d6fPamSvgF8ChhLMg/xAnBV9MQPK8pV3/6gJyqfYO61vVrbzAaPGkcuDbVng0a1uZZKXjnKQS3/sd8LTK9cmYmeOBM4UyWdB1wYPfHn3j7cS/AYUHv1tNnm7VmT1TDXUslzLzmoZbVoAlDtcmET8l1J6q+9vNtsdnvWfPWMWsq8ctSgxrdcKGl69MSmOfWnLdtsxc9o+ahzrqWS514a0NDt/yppbVgmA2Lhmt1mK35Gy9XBJJc5fT1BP7rKeysDhwLfKKBfHa/RpGib0vxnbZrdZit+RsvP6cAtVd6/EvgKLLsPc8bdufeoSzQaXFoxF9HsNj3fMohFxHNUSfkrCeDqiJjXtE51iYbnXMwGM0kBjHBwyV8uiejNzCo5uJg1iaSXJc1Lj6clrd3qPhXJwcWseTaLiJUjYmXgOZIn/DuWg4tZTiTdLSl6OZ4HiIin0nJDSZa5G5rwlLShpPMk3SdpsaRpVcqOl3SjpLmSnpV0bNqPusrUo9HVIjNbaizwKvDtivNv7j2UBpq3puU+0mB7mwAfB24juaO4V5LGkOT7ehD4BLABcBrJ4OLoWsvUy6tF1tXyXC2StBi4JyLe30+5FYBbgRcjYvcG2hsSEUvS7y8D3hoRO/RSbiLwPWBsRLyanvseMBlYIyJeraVMvf3zZZFZDiSNJ/l9Wu7p+koRsZBk/+edGmmzHFhqsBtwXUWAuJjk0mz7OsrUxcHFLB/lEcgdklbMHAKQNE7SJpny3ySZ1CV9X5KG9XcMsG/vAmZkT0TEkyTbiLyrjjJ1cXAxy0d5g7KfkjxKUD7OTc+PJQk88yTNBzYimS8p+yJJ4r3+joEYA7zcy/nZ6Xu1lqmLJ3TN8jGe5Je/MmneDQAR8SdglSqfvxL4QDFdA3pfmVLF+VrK1MzBxSwfawH/joipA/z8SxSX+3w2Sc6tSqNZOlqppUxdfFlk1qD0XpARZJacB6DIy6IZVMybSFqXZCQ1o44ydfHIxaxxH02/3tFAHUVeFl0LHCFpZES8lp7bj2QzrD/VUaYuDi5mjds5/fqHgVYQEbOAWfV8RtIIlk4Krw2MkvTp9PU1EVFO4HcucDjwG0knA+8guX/l9MzScy1l6uKb6Kyr5XETXXrb/fbAkGjiL5SkccC/+nh7/Yh4PFN2PHAW8CGSOZSfAZMjYnE9Zerqn4OLdTPv51IcT+iaWSEcXMysEA4uZlYIBxczK4SDi5kVwsHFzArh4GJmhXBwMbNCOLiYWSH8bJFZE0n6E8nWBgIeAQ4a6LM77c63/1tXa/bt/5JGR8Qr6fenA69HxKRmtN1sviwya6JMYBlCsldKx/51d3AxazJJ15Bszr0xcEqLu1MYBxezJouIjwNrkGwudWiLu1MYBxezHEjaPk3dulvm3PqSnpd0RmX5dI+UC4ADmtnPZvKErnW1nDMu3gSsFBEfljQa+AvJZk6fiIjFacrU4RHxXFr+GGB8RHym0bbbkZeizfLTA9wsaRfgOyQban8ms5PbGOASScNJlqIfAr7Rkp42gUcu1tXyXoqWdD3wYZJtIreKiKfzqHcw8pyLWb7+SZJmpKebAws4uJjlRtLBwEHAfcBXWtydlnNwMcuBpI+S7Jz/VeAQYKvsylFF2a+nK0vvaGYfm81zLtbVckotsgnwZ+CsiDg6PXc9MDoitqwoOwq4lyR963ER8X8DbbfdObhYV2s0uEh6G3A7cBewbzlvkaTtSDIV7hERV2fK/wB4g2Re5rWIOK7BH6FtObhYV2vmg4uS1gFuBjYH9gb2jIh9i263VTznYtY8JwAnRcQcYDpJkOlYHrlYV2vWyEXSBOCvwPMkT0ILWAsYGRHzi2y7VTxyMWuOHwKfjYixETEuIsaSPBqwSYv7VRgHF7OCpUvSoyLiioq3HqaDL418WWRdzYnoi+ORi5kVwsHFzArh4GJmhXBwMbNCeLMosyZy3iKzLuG8RcXxZZFZEzlvkZkVxnmLzKwQzltkZjVz3qLleULXuprzFhXHS9Fm+XHeogyPXKyrOW9RcTznYpYv5y1KObiY5cR5i5bl4GKWg1ryFknaW9IsSfdKelDSfZLGt6K/zeDgYtagNG/RpcApEXFhRNwG3ACUKopOAM6MiAkRMR64E/hsUzvbRA4uZg1I8xZdBVwPZJ8ROg74gKTdM+cmAH9LP/cWYFPglqZ0tAW8WmRdrcl5i54A5gCLgY2ASRFxatHttoqDi3W1JqYWGQP8MyJWT1+/kyRL47qduuWCL4vMmmMC8GD5RUT8g+SJ6NGt6lDRHFzMmmMC8ED5haR9gacj4qmW9ahgvv3frDneA3xM0gdJHgt4BtirtV0qludcrKs5b1FxfFlkZoVwcDGzQji4mFkhHFzMrBCe0DUAxh119cHA51rdj2ab/+Tft19xnfG3aMjQJa3uS5P96vGTdj+/yAY8crGyz5Hci9FVVnjrevejId0WWCbQhD8kvs/Fsu59/KTdd2h1J6xY4466eloz2vHIxcwK4eBiZoVwcDGzQji4mFkhHFzMrBAOLmZWCAcXMyuEg4uZFcLBxcwK4eBiZoVwcDGzQji4mFkhHFzMrBAOLmZWCAcXMyuEg4uZFcLBxcwK4eBiZoVwcDGzQji4mFkhHFw6iKR3SFKr+2EGDi4dQ9LmwKPA3yXt6iBjrebg0jlWAF4FNgEuw0HGWszBpbOU02eugoOMtZiDS2dzkLGWcXDpDg4y1nQOLt3FQcaaxsGlOy0XZCKin4+Y1actE9FLGgc8QrICYsUpB5nL35j5SKy41sZ3tbpD1jnaMrgAs4CrgfVa3ZFBZCyw+gA+txAYouErzc65P9bl2jK4RMRrwN6t7sdgImkL4EZgdI0fWQgsBn4JlIa/dezUovpm3clzLt1nITAfuBDYKCK+EhFPtbhP1oHacuRihVhmpOKAYkXriJGLpOmSFkvyksfyPFKxluiI4AKcC2ze6k60GQcVa6lBEVwk3S0pejmeB4iIsyLigRzb+7Skv0iaJWm+pIclHS1peC9lx0u6UdJcSc9KOlbS0HrL5GQYDirWJgbLnMtYkid+v11x/uGC2lsd+CNwKvAysCUwGVgDOKxcSNIY4AbgQeATwAbAaSRB++hay+TkZWAJSVDxnIq13GAJLm8B7omI/21GYxFxXsWpP0oaBXxd0jdi6e2shwArA/tExKvA9Wm5yZJOSc/VUiaPPj8qaUxELM6jPrNGtf1lkaTxJP28v8VdmQVUXhbtBlxXESAuJgkm29dRJhcOLNZO2j64ALunX++QtGLmqOmBOyWG9Xf08dmhkkZI2gY4HPhpLPsQzruAGdnPRMSTwNz0vVrLmHWcwRBctkm//pRkorJ8nFsuIGmGpEXp94skZX+Zv0gyydnf0ZvX0+MW4E/AERXvjyGZ66g0O32v1jJmHWcwzLmMJ/nl/3LF+RvK30REtRHAlcAHBtj2h4ERJBO6xwBnAYdWlOnt3hpVnK+ljFlHGQzBZS3g3xEx0GdfXgJeGcgHI+Lu9NtbJb0IXCDptIh4ND0/m2SyudJolo5Wailj1nHa+rIovRdkBI0tOTdyWZRVDjTrZ87NoGLeRNK6JFsZzKijjFnHafeRy0fTr3c0UEcjl0VZW6df/5U5dy1whKSR6ZPcAPsB80jmaGotY9Zx2j247Jx+/cNAK4iIWSTLyDWT9HuSOZ0HSB722xr4DvDrzCURJJPKhwO/kXQy8A6Sm+1Ozyw911LGrOO0e3B5f/r15ia3eydwIDAOWAQ8Bkwks0IFEBGzJe1EMtF7Jckcyo9IgkfNZcw6UVsHl4jYoUXtTgIm1Vj2QWDHRsuYdZq2ntA1s8HLwcXMCuHgYmaFcHAxs0I4uJhZIRxczKwQDi5mVggHFzMrhIOLmRWiY4LL1dL0mdJinLvIrC10THB5FM79uXMXmbWNtn62qB6HR5wFQG1b6/Yv2Rj8TOBDJA8b/gwo0d8m2NKBwM97eee/iDg3U25g9Rdne4BxR109rUXtt8TCl54ZN2zMWo/XuCVzp5gA3Ft0Ix0TXHKVT66hHUn2bCl7LOf6LQeL57w0dthb1ngSDV3S6r400b3Ar4puRMtuZt8BpCCisT9D0kTge8BYynuuSN+jnBit2j4sS0cuI4mYk3v9lqs0v/iIiJjXb2GrS9vPuZwhHY4UZ0tvboFwsrTtC9KSK6V7C2p2N+A6iss1VHT9Zi3X9sHl8Igz7oTZW6dpPY6X1t0Hrn8EXri71u0rJSEN6/dYarlcQ9Sfa+hRpEVIDyN9reK9POq3QUQlrauSblRJD6mkB1TSKSp19kRP2wcXgL/ApAkw8izpqJ3hnkUQN8HmPRFvbqz9e2nGs2nuomelRb9vLHdRI7mGZpJsNPUFYE/gduBcpG/lVL8NRjfxNqYwksnAcYgX+QxL2KfV3SrSoJjQ/WbE2bdLx34JTnwVlkyBrSZFPJcts2smd9Fa6ZExkE26B5ZrKOI64LrMmWuRVgSORvoJEeWJQ+cy6jAq6UfAf/fy1jXczFeAQyPiLknDeZZ/8RK7A5c30N5/kvwR24IkVc3DwA+jJy7qpWyvq5PRs3R1spYy9RgUwQXgRZi5Cqx2CUydGHFXnR+vN3dR3rmGLgP2JdmT97EC6rd28Dz7MA+4vuL8AlaOiJkko1qYzEjmMYrf8lxlFXX6Nkk2im8BLwIfB36lkt4aPXFmuZBK/a9O1lKmXoMiuEyRpn4WNnkI5m8BnyLZPLseX6T3e08qla+Bl8s1RD65hsqjkqLqt1YayRo8zTM8zacr3vl3+RuVtCKLuIK/sJhHuLDBFveMnngx8/omlbQWSdA5M3P+EJLFgn2iJ14FrldJo4DJKumU9FwtZerS9nMuZ0lHfg4+/ws4/zo4bHNYNbtylHWh9GukOFmqXHEpXxb1d5RdC3wMaWTmXCO5hj5F8pfliYLqtxZTSW9jZYbzOg9ExG0Vx+NpmaEs4SKmM45bOC4iHmqkzYrAUnYP8LaKc7sB11UEiMrVyVrK1KWtRy4/lPb6Kpz0O7j1yxFfA7hdOmU7OBI4Llv2eGntA2Cf6TB3FOxC9pe0/txFb+YaoiLX0DLLx9IBwP8DNiDiifTc5SRJ3O4HhpIEjf2AwzPzLbXVb4PHTLZjTeB1/qbSMiuPi6MnvZksOJ9/8j6u4NKIOC378XTlaGh/zURPLOqnyIdJLm2y3gXcVFHPkyqpvDp5ZY1l6tK2I5cfSOP3hcvvgaenw3bl87dDz2awyjnS5Gz5d8Mlt8Of/wHTV4OtGmo8YjawE8n/7CuBEkmuoZ6KkkPSMtklxYeBg0gm6i4FxgMHEEuvgeuo3waLBewEwK5MJLsCGXwBQCVtjTiItzOKHnZUSfeqpMMzNTScdlgl7UQyX3J2xVu1rE7mvoLZtiOX7ye5flZYD9ghcz59huisQzPnTpTevz9seQGssyaUNoDdG+5ALbmGIqYAUyrOfR/4fi712+AxjC2YxfJrP3NZiclkU+A9s8z7S/+cNJR2WCWNI7ml//+iJ6b0UqSW1clcVzDbNrjUYwu4+Ca4eFLEc6dJ09aDL7W6T9ZlRvAOnmcOz/KRinceBoiIP7PsCLdSvSuab1JJq5HM4z0JfL6XIrWsTua+gjnog8tp0n6HwgazYNxT0mf3A94OQ4+VRh8TMaD/WWb1UElDGMkYHuXBqP82ibJ6VzTLbY8ArgKGA7tHT7zey2eWW51UabnVyVrK1GXQB5dt4ZwL4KRDIiaWz/1LemNV2AP4ZQu7Zt3iDTZiOEN4lft6e1slHcjSm+veTXJvynxgMfDB6ImFDOCyKJ04vhR4J7B19MTzfRS9FjhCJY2MnngtPVe5OllLmboM6qeiz5YmbQdHbhaxavb8n6V/PwxXHxTx5Vb1zQaHPJ6K1kE6nLH8hEv4WjwQ5/dZrqQhwGvAetET9axe9lXf+cBXgW+SrFBm3RM9sSAtN4ZkBWk6UF6dPB34cfTE0bWWqdegHrl8PeI4KpakAbaOWGPrFvTHulSwDW8AD3NDPyXfCczOI7Ckdkm//qSX99YHHgeInpidriSdRTJCeplkdXJyuXAtZeo1qEcuZo3KZeQilYDvACOjyi+UStoX+GL0ROOrmYNA297nYjZYRERPRKxaLbCk3kNyc2VXcHAxa57NofdJ307k4GLWPO/BwcXM8pSuxvwH8Eir+9IsHTOhK2k6yT0EQ6LRDbqta3iD7uJ00sjlXJwUzaxtdExwiYizIuKBvOqTNF7SjZLmSnpW0rGS+n0kXtKBkqKX45A86jcbLAb1TXRFUcFJ0XKq36ytObj0bumWf5Fu+ad0yz/plKhtQ6c7o6+kaPnUb9bW2v6ySNLh6WXFpMy5bSUtUcFJ0aLgpGgF1m/Wcm0fXCLiDJK9Jo4AULKR9fXAC9T4FKkSw/o7Mh9ZLmlZDCApmqRFkh5WDUnRBlC/DTLSSVtKd74k/XOB9I/50m9vlyZ37Mpm2weX1CRgpKSjSDYgDmDzyCRFkzRDaVK09Je6rZKiyUnRjFgAt3wzYsMV4aLVYb2NYc1TWt2rogyKOZeIOFvSscCJwBJgq6hIihaZpGi9aFpStOglKZrSpGiSfhJOita1IibeR3qHbsQxr0tXPgGrbNhInRK9JzILqiYyk/g0SQqSjUk2hHoCmAqcEsEbjdYPgyS4pGYCqwFTB7Dbl5OiWVuRTtgQvrQJ/LL//Zb7rINGVh1XB/4InEryb25Lku0V1gAOy6F+iIi2P0giapAs7b42gM8fmH6+6pEpfzNwUUUd66bl9hxA+/+Zfnb9Iur30dC/rQBWbm6bpZFw52z4xe8aqycmQsyGGJU59z2IudlzddR3AsTLEMqj/rafc5F0JMmmw+eTRNRVVZEUTdJJ6erRPEkL0q97ZooMKCmaCk6KlmP91mLSGYdLhHR2ZlXz5G2lF5ZIV9679FxpBdj+IXjmsYj992qw2XTVkbxWHWeR7MWbT/2t/svRz1+VvUj+stySOTcLmFNRbhowLfN6BnBzA+2OIbkMux7YGTgYmAMcX1HuAGARMDZz7nKSpG27kezjWx51faPe+n005d9YbiMXuOMluOfV5Pvj1oVH5sOtz8HkFZaWuXYG/P4f0KPlPx+CGNbfkSn/PMTkXup5HeKI2vocQyFGQGwD8TjEaXnV3/L/uVX+p48nWcF5ivQBy/T8Yek/iMmZczOBY9Pvx5LsU3pUDu3fRDKamEmynebQijIHpn0Zlzn3A5J0EnPTz/4N+MJA6vfRlH9nOQaXH38dIuDMo+CvL8KD8+DYt2fePyR5f8b85L0H58GFly19Pw5M3q9+ZMovhPjv5fsRT0P8oLY+x/xM3RdADMmr/o54Kjpdgl5E8g9lJeDqiNijtb2ywSDvp6Kl22fBpqvBq0tgylYRE2tefJBYnWTv26oiuCstvxD4bsSye+hKPANMieB/amjzfcAIkgndY4BfRXBoHvUPptWiXkkaR7LNwkrp651JbqdfOyKeqfphs9y9OBNWWQ0umVpPYEk1fVUzgrvTb2+VeBG4QOK0CB5ttP62n9Ctwd4kl0EARER5B/a1W9Md61bSlKmw8ybw0HzY4lMDqKLemz2XT2QmGklkVg405dFTQ/V3QnDZHni6/ELSj4AFEVGZx8WsMNJZR8LnPg+/OB+uOww2XzW7crRs2Qt/nawsnVy54jLAVU3yWnUsZ+T5Vy71t3pCLYcJuX+SZK6bB7xOMjm6fav75WNwHOQwoQun7gUvB1yaWdW8bRbcP2f5ssetDU8shL+/Duec0Fi7MQZiJsT1EDtDHAwxB6JiVTMOgFgEMTZz7vcQ34XYDWIXiFL62Yvrrb+vY9CPXCJiw4gYGhErR8QqEbFmRPheEWsK6QfjYd/L4Z6nYfp2S9+5vQc2W0U6Z/Kyn3j3JXD7n+Ef02G1rRppO4LZwE7AUJJRT4kkkVlPRdEhaZnsQ5J3kqx2XgpcQvIc3ESSZ+Lqrb9XHbFaZDZQzdxDVzrx/bD/X+GCdWDNEmywe8RH1i263VYZ9CMXs8Fji4vhposjJj0Hr0yD9d7W6h4VySMX62rNGrlIp+0Hh14MsxYn0zwC3j4UTnxLxDH1LD8PGh65mDXFtufABSdFrDMsYt1hEesMg6cXwqode7Ong4tZwZIl6ZVXjDhk4rLvPPsSvGXH1vSqeL4ssq7mpGjF8cjFzArh4GJmhXBwMbNCOLiYWSEcXMyaRNKWkl5Kt2KdL+l2Sc5bZGYNWwB8MyJWJNl9f2OgY/MWObiYNUlE3BcRU9PvXyfZsL3BvEUaL+lGSXMlPSvpWElDq5SflqZH7u34UCN1Vxr0O9GZDUaSNgQ2ARrIW6SB5BU6FBhVce5Y4L0kT0o3UvcyHFzMmixNKXMncE1EXN1AVYeQpPnYJyJeJdnedRQwWdIp6bllRMSDFX0ZDrwf+HVELGqk7kq+LDLLgaTD00uLTN4ibZvm07o3c24F4CHgsYjIKW9RNJK3aFeSVDcX5V23g4tZDiLiDJINrY8AkLQuSV6qF1h2a8q/k+ya+P7KOpQY1t+R+ci7qNjLNiKeJElrUy13etZngGeAWyrON1y3g4tZfiYBIyUdBdxDsrfC5hGxEEDSISQrROsCc9PMoJdlPl/vBt1j6H0X/tnpe1VJGkGyA92vY/mHDBuqGzznYpabiDhb0rHAicASYKuIeC7z/rnAuVWqKG/QXVezvZxTH+cr7QmsyvKXRHnU7eBilrOZwGrA1Iho97xFnwH+2Uc/G86J5Msis5xImkqyvDwfaFHeItWUV0jSaJJJ275GLQOuu8zBxSwHko4EPg+cT5LPfNXsylFa5qR09Whe+gjAPEl7ZooMMG+RBpJXaG9gRfoOLo3UnWh13hgfPlp5kEveIvZK68nkLWIWMKei3DRgWub1DODmBtodQ3IZdj2wM3AwMAeoyFvEASS51Mdmzv0euLfRuqsdHrmYNUDSeOBykqyfmbxF9ACrSJqcObcxcHP6ubEkKYevGWjbETGgvEWS3pp+7uIc6u6Tt7m0rtbcvEVaRDKCCGAl4OqI8AbdZjZwksYBQyJipYhYGfgosLuktVvbs+I4uJg1x97Aa+UXEXFD+q2Di5k1ZHuSeRkAJP0IWBARd7SuS8XyTXRmzbEpsL6keSR3774KfKy1XSqWJ3StqzlvUXF8WWRmhXBwMbNCOLiYWSEcXMysEF4tso4maXOSe0yqmSRpQR/vXRkRd+fcra7g4GKd7tMkO+xXS4kxsY/zS0h2yndwGQAvRVtHk7Q68CQwYgAfnwesH5nd5Kx2nnOxjhYRs4CzSDZwqscC4GcOLAPnkYt1vAGOXjxqaZBHLtbxBjB68aglBx65WFeoc/TiUUsOPHKxrlDH6MWjlpx45GJdo8bRi0ctOfHIxbpGDaMXj1py5JGLdZV+Ri8eteTIIxfrKlVGLx615MwjF+s6fYxePGrJmUcu1nV6Gb141FIAj1ysK1WMXjxqKYBHLtaVMqOXwKOWQnjkYl1L0mrAFOCrDi75c3Axs0L4ssjMCuHgYmaFcHAxs0I4uJhZIRxczKwQDi5mVggHFzMrhIOLmRXCwcXMCuHgYmaFcHAxs0I4uJhZIRxczKwQDi5mVggHFzMrhIOLmRXCwcXMCuHgYmaF+P97eWGu4sIc6gAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 360x360 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot_flowsheet()\n",
    "XA1=1-XB1\n",
    "plt.text(xs-la-hla,ys+ly/2-14, r\"$x^1_A$\"+'={:.2f}'.format(XA1), color='r',size=15)\n",
    "plt.text(0,95, r\"$x^1_A+x^1_B=1$\", color='r',size=12)\n",
    "XA2=1-XB2\n",
    "plt.text(xs+lx/2+5,ys-14, r\"$x^2_A$\"+'={:.2f}'.format(XA2), color='b',size=15)\n",
    "plt.text(0,90, r\"$x^2_A+x^2_B=1$\", color='b',size=12)\n",
    "FT2=FT1-FT3\n",
    "plt.text(xs+lx/2+5,ys-6, r\"$F^2_T$\"+'={}'.format(FT2), color='g',size=15)\n",
    "plt.text(0,85, r\"$F^1_T=F^2_T+F^3_T$\", color='g',size=12)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "id": "77b1bb53",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAARcAAAEvCAYAAABv+oUUAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAApXUlEQVR4nO3deZhcVZnH8e8vCSQhJCHgsIeExSABQiRIQLZiEcEICCqgaERQiMjmOEgyQ+gUoECQoOzkEQhkhh1lDIsxLCWgwx6EsINsgQhkgZiVLO/8cW7R1ZXu6qq691Z1db+f56mnu26de87J0m+fe86975GZ4ZxzSetW7w445zonDy7OuVR4cHHOpcKDi3MuFR5cnHOp8ODinEtFvOAi9Ud6AmkR0g4J9anjtemcq1jckcsSYBRwRwJ96chtOucqFC+4mK3A7KOE+tJx23TOVcznXACkk5GeQlqONKXe3XGuM+hR8xalCQCYTah52217HzgP+CrQu859ca5TKG/kIk1E+kPB+4uQHkBaK5Ve1bo9s99jdhcwL5X6neuCyh25XAi8gTQc2A04CNgTsxVI9wLDgW2RrsFsSgL9ars9oN02pbuBPduo+1HMvp5AH51zJZQXXMzmIf0GuBHoT/hB/yT67GuJ96pUe+W06cHDubqrZEJ3JrAjMA6zdytqRbob6WOkj4GxwNjP3odRRrLtOefqTmXlc5F2BKYDfwY2x+yAos+3B54DBmE2u526JgClJ3RLtSfdDmwavVsOfA2zZUXn3wfs1Ubtj2B2cBvtnhe1d2zJP4Nzrl3tXxZJmwHTgDHA/cCbSBnMcgWlzgSmAtsBpYNL/PaGACMwW4l0FzAMeKJFHW0Fj7bb7EH4u+gOdEfqBazEbGXVfw7nurjSl0VSP+BeYBJmf8RsCXAR8MuCMsOAOYSRxnaxetNee9LaQLcosKwLbAK8GqvN4CxgKeGS7XvR92clUK9zXVZ5l0Ula9D1wM+BgcAYzH6SQL/aamsnwmjmBWBb4EzMbkytPedc1eIFl7BUfBvwD/KXFGb7JdKz1ts7BtgWs7OR+gOPYRZvtOScS0XcO3THAiMxWwCA9ETp4rHtADwZfT8A+KREWedcHVUfXKQRwNLPAkuwDGkDzNK603V74ECknwKrCZO+zrkOKNZlkcTuwCTgU8LzOaPNWJFQ3zpMm865ysV9KvptYD8z9iHMuxwWv0sdsk3nXIVizbmY8X7B25WES5VU1aNN51zlEsnnIrElcDDQ1q38iUuyTYmTJZ6SWC4xJXbnnHPx87lI9ANuAL5vxqdllJ8AYBa+1qLNMng+F+cSVtbIRWKixB8K3l8k8YDEWsDNwAQzXkmqU6Xak+iRdJtm/N6Mu/B8Ls4lptzLoguBfSWGS4wh5Fc5AjgaGAmcLZGTOCqhfrXaXrQq9J322pS4W+LjNl41u3Rzrisr67LIjHkSv6Egv4oZnxAeVpyadKdKtIdZ+22a4flcnKuzqvK5mFFRfpXCkQRRPpcyRhJVt+ecq7+ybqKTaJFfxYwDij7/LJ+LWemUC+VM6JZqT2KNfC5mLCs6v2Q+FzNaTckgcV7U3rGl/gzOufa1e1kksUZ+FYmMGbmCYonlcymjvSHACDNWStxFK/lc2goeJdpskc9Fohew0gzP5+JclUpeFkVLvvcCk8z4oxlr5HORSCyfS3vtSawNdIsCi+dzca4Di53PRaJFPhczUsvnIrFGPhczPJ+Lcx1Q3AcXh1OUz8WM1PK5SBwDbGvG2RL9gcfMYma/c12WpK2BvwBfMbOX6t2fzibu7f9jgZFmHGTGV4B1E+hTKTsAz0bfez4XF9d5hEvr8+vdkc6o6pGLxAj6/PNnnLHJEGAosBsT7ErgcLN07nTVkHseYf//+iLWbQWLN1zOXdd/3f61yVNptPVZm1ltBPwBWAGsAo6xJpuTZpsufdGoZRbQizDHNsJHL8mqeuRixtOcsckPgVHAHdGxvdMKLAAc8/Wj2Pjv/2ZXPzOAbaZfy39sOji1tprNBfa0JtuHcFPf8TVo06XvPJpXS9fGRy+Ji5dyoclWAB8pq4S60257tU/x0GSrCt72JUwmuwYWjVq+QfP//+7AgZK289FLcsoKLspqADAfWFxw+BNgc2uKu31A5e0pq3y6hV+2dn4K7Q0HrgHWAw5Mok1XV4Wjlrz86OUbNe9NJ1XuyGU4MN+abIMU+1JWe8qqOd1Ck5VMt6Bs2N3Rmkrs7thOe9H5zwIjldWRwDg8d2/DamXUkuejl4SVO+cyHHgxxX6U1Z6yak630GSJpXhoq72ozZ4Fbz8BliTYrqu91kYteT73kqByRy5fpO0fvnsJP5zbKqtrrMmmFH1+N7BnG/U+ak3W2hPMbbXXnG4hq7OBq6zJbi3rT1Bam38+YGdldSFhpWgZcFwC7bk6KDFqyfPRS4LKe3Axq+eAbeCzrG/vWpPtmFqnEmqv3MuiWv/5XH1Iuhn4FqV/qa4C7jazb9SkU51Yu8EluixYBOxuTZbqPSVJtFc0UuoVfc0/Nb3GSKnWfz5XH0X3tbTH73tJQDmXRTsQlnyfKzyorE4BvgkMIsxDfATcbU3266JypdMfNFnxE8yttldum4XBo8yRS6z2XMMoNddSzFeOElDOX/YXgVnFKzPWZJcBlymra4Abrcn+2trJrQSPqtqrpM0O3p6rsTLmWor53EsCylktGg6UulzYnmRXktprL+k2a92eq71KRi15vnIUU/yUC1nNsibbIaH+dMg26/FndMmocK6lmM+9xBDr9n9ltRm02AExdbVusx5/RpeoEwiXOW09Qd+/xGe9gZOAU1LoV6cXd1O0Haj9sza1brMef0aXnEnAIyU+nwb8CFrmYS7wTOI96iLiBpd6zEXUuk2fb2lgZvYBJbb8lQRwj5ktrVmnuoh4mehyuf7ADKJ8LpbJzErzvGopl9ud8BvsU8IlzmjLZFYkfU7MPq6ZNyaTaTdvTLXnuUCSAet4cEle3Ex0SyjI51KD86r1NrCfZTL7EFJyHpbSOXGEvDGhvUryxlR7nnOpipfPJfwm/0i5XE3Oq5ZlMhXnganmnDgsk6kqb0y157nak/QxkH8Qdh4w0szeq1+P0lVePpdcru18J5lM4vlcqlFOH5XLVZwHpppzqu0nsBNV5I1RLje8mvNcze1oZu8CSHqa8IT/3vXtUnoqy+eSydQqn0s1hlOij8rlmvPAZDIl88BUc45yuQkAlslMiNHPZ4GRyuUqyhtjmUxV57lkSXqGcMd3sY/MbMOCwNKdsMy9KGZ72wBnALsRVjUfMbNMG2WHApcBuwMfA78DsmbNmRbLKVOJSoJLRSsmyuUmAp+3TObw6P1FwM7AQSlNjA6nrbQQuVxzHphM5pWC4232EbDWzkmrn8rlelomszx62yJvTDv97NbWea7mBgELgX8vOt78f076EPhcVG7fmO1tD3wNeIxwR3GrJA0g7Pf1ImHucGvgYsKc61nllqlU/HwuuVxzPpdc7hrLZKZEH10IvBEN2Xcj/CDsmQ8sJc7L11s6D0wmU5wHplROluY8MLlcyAOTydxaqo/K5b7fxjlxtdXPnZXLtZU3plQ/dy9xnqut9YCZZnZtWwXMbENJawGPAtcRFjaqNc3M/hdA0h2EoNWaMYSR0hFmthCYIakfMEHSxOhYOWUqUl4+l1xuzXwnmUy7+U6iS4UjCHdB7mmZzLuVdrBc9e5juZdF9e6naymppejokuIF4HozazfASzoUuM3MqnksobX67gA+19plkaSHgffN7OiCY1sQVkQPNbNp5ZSptE/tb0Sfy/Uk7AG9u2UyleY7mQk0Ee69SDOw1KWPRaOrXtGx06P3a4yuGuHv0lUtPwJ5QmqRGvVTMzNJg4E+ZpZfzTsN+CBfSOFuvu7tNWJmK6vo2xeAB4vqeUfSkuizaWWWqUj1+VxyudZznWQyv44+3xG4ijAhehxwU9H520d1DrJMZnZxo8rlSueByWQKUzm0nZOlRD/b62N7/SwMHmWOXOryd+lqIv9L5qrolTcZOJHwb3uvpG6ACP/GBxWU+wFwfRntVLOPzwDCBG2xBdFn5ZapSPn5XIpWSyyTCblOcrmQ6yST+SzXiXK5zQiRbgxhkuhN5XIZy2RyBVWcCUwl/CZf4weiKHhU1cdS/Syzj+32s0J1+bt0NTGUcJd08U2M9wOY2V+APiXOnwZ8KZ2uAWGBopiKjpdTpmzlBJfhVJDPJVq+vReYZJnMH6NjFxHuE9kjej8MmENYet2O8ChAHO31sUU/y+ljnfrZEf4uXXU2Bf5pZlOrPH8+6e19voAw2VysP82jlXLKVMbMYr146KFZVZxzPQ89tD4PPbQTDz10Vdw+dJZ+NkIfO9uL8Fu5d8w6ukf1zIhRx7FRHSVfJc6/A8i18dnDwM1FxwZGdR5SbplKX/HyuYQhe0W5TqLl1D0I8wbdKWMSK65G6Gcj9NG16SvR1ydi1JHmZdF9wBmS+prZv6JjRxGSYf2lgjIVqUc+l7HASMtkFgAol4vzD1KuRuhnI/TRte6A6Oufq63AzOYRnjcqm6R1CDfRAWwG9JP0rej9vWaWv6HyauBU4PeSLgS2AiYAk6z5/pVyylQk7lPRFeU6US43Alia/2GILFMul/ZjBY3Qz0boo2vdLtHXh2vc7obA7dFrN8Kkcv79hvlCZrYA2J8wsp0GZIFLCLc2lF2mUrFz6DrXyDyfS3rijlycc65VHlycc6nw4OKcS4UHF+dcKjy4OOdS4cHFOZcKDy7OuVR4cHHOpcKDi3MuFXGfLXLOVUCTJ39Mr14hU938+fN4/vmRdu21nXLvIh+5OFdLjz++o40e3dtGj+7N4sUfsMUWN9e7S2nxkYtzNWTXXhv2Ljr22O5st11vVq6Mt3dRLtf6XkMtd+Js69xvAOcA2xLSfVxmmcykpOr3kYtzNaYrrviQww9fwXrrbcrbbx9ZdT1h9877CQmdDiMEip8Tnmhu79w9gN8TctAcQtjm5MKCBPOx6gd/Ktp1cfV6Klo//OFabLXVo3z66Vw799yq9i5SLjcO+AUhMfvC6NgvCHlYNs4fa+Pc6UBvy2T2Ljg2iZARb2PLZD6NUz/4yMW5ROiUU05VLmc6/fTxnx078cS9dNddq/WrXz1bXN6uv34F77//S4YM2T9GswcD04t+yG8hbG62TzvnDidKHl7gz4RM/7snUL/PuTiXBLvsskt19dUT2GGHM4BzdfzxA9l33xnMnv0Rr732JQAdf/xgpD72u9+FjIMbbngac+c2712Uy5W3d1Emk9+7aM29hjKZd5TLlbPXUC+aN+bLy28LvB0htWWc+n3k4lxiXnhhPNts01ennTaWESNmsmqV8cwzw+z668Pe6D16DGKPPZ7QjTcu1dSpy1hvvSG88MLXCmr4AWF7kvZeeXH2GnqdNXP27hp9XT+B+n3k4lxS7LLLrtBVV53DwQefz+LFq5k+faRdd91nIxO75po09i6qdq+hq4GrlMv9mLBzwK6EyVoI+47Hrd+Di3OJ+uSTOfTqtT4PPTTVJk+udMveSvcuirPX0HXAToTdIScTdvo8k7DsnA+IsfYy8ssi5xKisWOnsvPO2/P228sYMuSbVVRR6WXRy4S5j+Y+5HIDCaOjl0s1ZJnMKstkTgb+DRgGbAQ8Fn2c/1p1/eDBxblE6LTTzmT//b/HjBmTefLJk9l663ULV45alB037lblcqYTTyxecclfFrX3yrsP+Kpyub4Fxyraa8gymQWWyTxvmcwi4CTgb5bJ5ANHrPr9PhfXpSVxn4t+8pNDGTXqf3n66UetqWkvAF111Tx69uxpxx23bouyxx+/GQce+BZLlnzKc8/9xi655L+qbjfc5PYiMAvI7zU0CfiNZTJnFZQbTbgM2toymbejY7sBexK2Ae4HfAf4KrCnZTLPVVJ/W3zk4lwM+tGPhrLvvnfy+uuzefPNz25I48UXm9hqqz46/fQJLU4YNOg2Xnrpr8yePYt+/UbGaTvas6qcvYa6RWVUcGwFYRRyFzAFWAfYIx9YKqy/VT5ycV1aLe/Q1Qkn7ML++/8f06dvzgYbZNl001H2s58NTLvdevGRi3O1MmTILcyceYtdd90HLF6cY6ONNmz/pMblIxfXpdVq5KKTTjqKww67hYULV2EGEgwY0J2bblrPrruukuXnhuEjF+dqYccdr2T69Avs29/uYUce2cO+/e0efPTRCnr3/nq9u5YWDy7OpUynnz6enj172qRJ41p8MHfufPr02a9O3UqdXxa5Ls03ok+Pj1ycc6nw4OKcS4UHF+dcKjy4OOdS4SkXnKshSX8hpDEQ8CpwnJmVzEXbqHy1yHVptV4tktTfzD6Jvp8ELDazVp+ebnR+WeRcDRUElm6EvCid9re7BxfnakzSvYRsb9sCE+vcndR4cHGuxszsa8DGhA3JTqpzd1LjwcW5BEjaR5JJOrjg2JaSPpR0aXF5M1sF3ACMrmU/a8kndF2XluSErqQHgV5m9mVJ/YG/AW8Ch5nZKkkDgLXN7IOo/NnAUDM7Om7bHZEvRTuXnCbgYUkHErbpWAEcHY1SIOz1c5uktQlL0S8Bp9SlpzXgIxfXpSW9FC1pBvBlwtYbI81sdhL1NiKfc3EuWa8T8tE2deXAAh5cnEuMpBOA44C/Az+qc3fqzoOLcwmQ9BXgcuDHwBhgZOHKUVTmcEnzJD0r6UVJf5c0tB79rQUPLs7FJGl74HZgopndaGaPAfcTtuIoNBy4zMyGm9lQ4EnCfkGdkgcX52KQtCFwNzADKHxG6FzgS5JGFRwbDjwdnbcesAPwSE06Wge+WuS6tJruWyS9DSwCVgFDgPFmdlHa7daLBxfXpdVsa5FwA93rZrZB9P7zwFPAwM6acsEvi5yrjeGEfZcBMLPXCE9E969Xh9LmwcW52hgOvJB/I+lIYLaZvVu3HqXMb/93rjZ2Ar4qaTfCYwHvAYfWt0vp8jkX16X5vkXp8csi51wqPLg451LhwcU5lwoPLs65VPiErgNg8Nh7TgC+W+9+1Nqyd57fp+fmQx9Rt+6r692XGrvprQtGTU6zAR+5uLzvEu7F6FLW+twWz6FuXS2wDKcGv0j8PhdX6Nm3LhiVqXcnXLoGj70nV4t2fOTinEuFBxfnXCo8uDjnUuHBxTmXCg8uzrlUeHBxzqXCg4tzLhUeXJxzqfDg4pxLhQcX51wqPLg451LhwcU5lwoPLs65VHhwcc6lwoOLcy4VHlycc6nw4OKcS4UHF+dcKjy4OOdS4cHFOZcKDy6diKStJKne/XAOPLh0GpKGAW8Az0s6yIOMqzcPLp3HWsBCYHvgDjzIuDrz4NK55LfP7IMHGVdnHlw6Nw8yrm48uHQNHmRczXlw6Vo8yLia8eDSNa0RZMysnVOcq0yH3Ihe0mDgVcIKiEtPPsjc+emcV63npts+Ve8Ouc6jQwYXYB5wD7BFvTvSQAYBG1Rx3gqgm9butSDh/rgurkMGFzP7F3B4vfvRSCSNAB4A+pd5ygpgFfA/QHbtzw2amlbfXNfkcy5dzwpgGXAjMMTMfmRm79a5T64T6pAjF5eKFiMVDygubZ1i5CJplqRVknzJY00+UnF10SmCC3A1MKzenehgPKi4umqI4CLpGUnWyutDADO73MxeSLC9b0n6m6R5kpZJekXSWZLWbqXsUEkPSFoi6X1J50jqXmmZhPTAg4rrIBplzmUQ4Ynffy86/kpK7W0APARcBHwM7ApMADYGTs4XkjQAuB94ETgM2Bq4mBC0zyq3TEI+BlYTgorPqbi6a5Tgsh4w08yurUVjZnZN0aGHJPUDfirpFGu+nXUM0Bs4wswWAjOichMkTYyOlVMmiT6/IWmAma1Koj7n4urwl0WShhL6+VyduzIPKL4sOhiYXhQgbiEEk30qKJMIDyyuI+nwwQUYFX19QlLPgldZD9wp6NHeq41zu0taR9KewKnAVdbyIZwvAC8XnmNm7wBLos/KLeNcp9MIwWXP6OtVhInK/OvqfAFJL0taGX2/UlLhD/MPCJOc7b1aszh6PQL8BTij6PMBhLmOYguiz8ot41yn0whzLkMJP/zHFx2/P/+NmZUaAUwDvlRl218G1iFM6J4NXA6cVFSmtXtrVHS8nDLOdSqNEFw2Bf5pZtU++zIf+KSaE83smejbRyXNBW6QdLGZvREdX0CYbC7Wn+bRSjllnOt0OvRlUXQvyDrEW3KOc1lUKB9otiw49jJF8yaSBhJSGbxcQRnnOp2OPnL5SvT1iRh1xLksKrRH9PXNgmP3AWdI6hs9yQ1wFLCUMEdTbhnnOp2OHlwOiL7+udoKzGweYRm5bJL+RJjTeYHwsN8ewM+BWwsuiSBMKp8K/F7ShcBWhJvtJhUsPZdTxrlOp6MHl12irw/XuN0ngWOBwcBK4B/AOApWqADMbIGk/QkTvdMIcyiXEIJH2WWc64w6dHAxs0yd2h0PjC+z7IvAfnHLONfZdOgJXedc4/Lg4pxLhQcX51wqPLg451LhwcU5lwoPLs65VHhwcc6lwoOLcy4VHlycc6noNMHlHmnWHGkVvneRcx1Cpwkub8DV1/veRc51GB362aJKnGp2OQDlpdZtX0gMfhmwO+Fhw98BWdpLgi0dC1zfyic/wezqgnLV1Z+efQAGj70nV6f262LF/PcG9xiw6VtlpmTuLIYDz6bdSKcJLolKZq+h/Qg5W/L+kXD9LgGrFs0f1GO9jd9B3VfXuy819CxwU9qNqGUy+05AMszi/RqSxgG/AAaRz7ki/YL8xmil8rA0j1z6YrYo8fpdoqL9xdcxs6XtFnYV6fBzLpdKpyLZFdJnKRAulPb6SFo9TXo2pWYPBqaT3l5DadfvXN11+OByqtmlT8KCPaJtPc6TBh4BM16Fj54pN32lJKQe7b6arbHXEJXvNfQG0kqkV5BOLPosifpdA1FWA5XVA8rqJWX1grKaqGznnujp8MEF4G8wfjj0vVwaewDMXAn2IAxrMvsssfafpJffj/Yuel9a+ad4exfF2WtoDiHR1PeBQ4DHgauRfpZQ/a4RPciGTKEvE4BzEXM5mtUcUe9upakhJnRPM7vicemcH8L5C2H1FBg53uyDwjIHFexdtGn0KlBNku7q9hoymw5MLzhyH1JP4Cyk32KWnzj0vYw6GWV1CXB6Kx/dy8P8CDjJzJ6StDbv8ybzGQXcGaO9bxN+iY0gbFXzCvBra7KbWynb6uqkNTWvTpZTphINEVwA5sKcPrD+bTB1nNlTFZ5e6d5FSe81dAdwJCEn7z9SqN91BB9yBEuBGUXHl9PbzOYQRrUwgb4spR9/4IPiKir074TdKH4GzAW+BtykrD5nTXZZvpCy7a9OllOmUg0RXKZIU78D278Ey0bANwnJsyvxA1q/96RY/hp4jb2GSGavofyoJK36XT31ZWNm8x6z+VbRJ//Mf6OserKSu/gbq3iVG2O2eIg12dyC9w8qq00JQeeyguNjCIsFR1iTLQRmKKt+wARlNTE6Vk6ZinT4OZfLpTO/C9/7b5g8HU4eBusWrhwVulG6FckulIpXXPKXRe298u4DvorUt+BYnL2Gvkn4zfJ2SvW7OlNWG9KbtVnMC2b2WNHrrahMd1ZzM7MYzCOca2YvxWmzKLDkzQQ2LDp2MDC9KEAUr06WU6YiHXrk8mvp0B/DBX+ER483OxHgcWni3nAmcG5h2fOkzUbDEbNgST84kMIf0sr3LvpsryGK9hpqsXwsjQauA7bG7O3o2J2ETdyeA7oTgsZRwKkF8y3l1e8axxz2ZhNgMU8r22LlcZU1RTeTGZN5nZ25i9vN7OLC06OVo+7tNWNNtrKdIl8mXNoU+gLwYFE97yir/OrktDLLVKTDjlx+JQ09Eu6cCbNnwd75449D047Q50ppQmH57eC2x+Gvr8Gs9WFkrMbNFgD7E/6xpwFZwl5DTUUlu0VlCpcUXwGOI0zU3Q4MBUZjzdfAFdTvGsVy9gfgIMZRuAJpfB9AWe2BOI6N6EcT+ymrZ5XVqQU1xN52WFntT5gvuaLoo3JWJxNfweywI5f/DHv9rLUFkCk4Hj1DdPlJBcfOl3Y5Bna9ATbfBLJbw6jYHShnryGzKcCUomP/CfxnIvW7xtGDEcxjzbWfJfRiAoVb4L3X4vPmXyexth1WVoMJt/T/rzXZlFaKlLM6megKZocNLpUYAbc8CLeMN/vgYim3Bfyw3n1yXcw6bMWHLOJ99i365BUAM/srLUe4xSpd0fyMslqfMI/3DvC9VoqUszqZ+ApmwweXi6WjToKt58Hgd6XvHAVsBN3PkfqfbVbVP5ZzlVBW3ejLAN7gRav8Nom8Slc0822vA9wNrA2MsiZb3Mo5a6xOKrvG6mQ5ZSrS8MFlL7jyBrhgjNm4/LE3pU/Xha8D/1PHrrmu4lOGsDbdWMjfW/tYWR1L88112xHuTVkGrAJ2syZbQRWXRdHE8e3A54E9rMk+bKPofcAZyqqvNdm/omPFq5PllKlIQz8VfYU0fm84c0ezdQuP/1X65ytwz3Fmx9erb64xJPFUtI7TqQzit9zGifaCTW6zXFbdgH8BW1iTVbJ62VZ9k4EfA6cRVigLzbQmWx6VG0BYQZoF5FcnJwG/sSY7q9wylWrokctPzc6laEkaYA+zjfeoQ39cF2XsyafAK9zfTsnPAwuSCCyRA6Ovv23lsy2BtwCsyRZEK0mXE0ZIHxNWJyfkC5dTplINPXJxLq5ERi5SFvg50NdK/EApqyOBH1iTxV/NbAAd9j4X5xqFmTWZ2bqlAktkJ8LNlV2CBxfnamcYtD7p2xl5cHGudnbCg4tzLknRasy/Aa/Wuy+10mkmdCXNItxD0M3iJuh2XYYn6E5PZxq5XI1viuZch9FpgouZXW5mLyRVn6Shkh6QtETS+5LOkdTuI/GSjpVkrbzGJFG/c42ioW+iS4tS3hQtofqd69A8uLSuOeWfRSn/FKX8kyZaeQmdnrS2NkVLpn7nOrQOf1kk6dTosmJ8wbG9JK1WypuiWcqboqVYv3N11+GDi5ldSsg1cQaAQiLrGcBHlPkUqYIe7b0KTllj0zKrYlM0SSslvaIyNkWron7XYKQLdpWenC+9vlx6bZn0h8elCZ12ZbPDB5fIeKCvpLGEBMQGDLOCTdEkvaxoU7Toh7pDbYom3xTNYcvhkdPMtukJN28AW2wLm0ysd6/S0hBzLmZ2haRzgPOB1cBIK9oUzQo2RWtFzTZFs1Y2RVO0KZqk35pvitZlmY37O9EdumZnL5amvQ19tolTp0TrG5kZJTcyk/gWYQuSbQkJod4GpgITzfg0bv3QIMElMgdYH5haRbYv3xTNdSjSL7eBH24P/9N+vuU26yDOquMGwEPARYT/c7sS0itsDJycQP1gZh3+RYioRlja/VcV5x8bnV/yVVD+YeDmojoGRuUOqaL9b0fnbplG/f6K9X/LgN61bTPbF55cAP/9x3j12DiwBWD9Co79AmxJ4bEK6vsl2MdgSqL+Dj/nIulMQtLhyYSIuq6KNkWTdEG0erRU0vLo6yEFRaraFE0pb4qWYP2uzqRLT5Uw6YqCVc0L95I+Wi1Ne7b5WHYt2OcleO8fZsccGrPZaNWRpFYd5xFy8SZTf71/c7TzW+VQwm+WRwqOzQMWFZXLAbmC9y8DD8dodwDhMmwGcABwArAIOK+o3GhgJTCo4NidhE3bDibk8c2Puk6ptH5/1eT/WGIjF3hiPsxcGL4/dyC8ugwe/QAmrNVc5r6X4U+vQZPWPN8E1qO9V0H5D8EmtFLPYrAzyuuzdQdbB2xPsLfALk6q/rr/45b4Rx9KWMF5l+gBy+j4ydF/iAkFx+YA50TfDyLkKR2bQPsPEkYTcwjpNLsXlTk26svggmO/ImwnsSQ692ng+9XU76+a/D9LMLj85qdgBpeNhf+bCy8uhXM2Kvh8TPj85WXhsxeXwo13NH9ux4bPS78Kyq8AO33NfthssF+V12dbVlD3DWDdkqq/UzwVHS1BryT8R+kF3GNmX69vr1wjSPqpaOnxebDD+rBwNUwZaTau7MUHiQ0IuW9LMuOpqPwK4D/MWubQlXgPmGLGf5XR5s7AOoQJ3bOBm8w4KYn6G2m1qFWSBhPSLPSK3h9AuJ1+MzN7r+TJziVu7hzosz7cNrWSwBKp+aqmGc9E3z4qMRe4QeJiM96IW3+Hn9Atw+GEyyAAzCyfgX2z+nTHdVXSlKlwwPbw0jIY8c0qqqj0Zs81NzITcTYyywea/OgpVv2dIbjsA8zOv5F0CbDczIr3cXEuNdLlZ8J3vwf/PRmmnwzD1i1cOWpZ9sZbw8rShcUrLlWuapLUqmN+R543E6m/3hNqCUzIvU7YuW4psJgwObpPvfvlr8Z4kcCELlx0KHxscHvBquZj8+C5RWuWPXczeHsFPL8YrvxlvHZtANgcsBlgB4CdALYIrGhV00aDrQQbVHDsT2D/AXYw2IFg2ejcWyqtv61Xw49czGwbM+tuZr3NrI+ZbWJmfq+IqwnpV0PhyDth5myYtXfzJ483wY59pCsntDxju9vg8b/Ca7Ng/ZFx2jZjAbA/0J0w6skSNjJrKiraLSpT+JDkk4TVztuB2wjPwY0jPBNXaf2t6hSrRc5Vq5Y5dKXzd4Fj/g9u2Bw2ycLWo8z2HZh2u/XS8CMX5xrHiFvgwVvMxn8An+Rgiw3r3aM0+cjFdWm1GrlIFx8FJ90C81aFaR4BG3WH89czO7uS5eeG4SMX52piryvhhgvMNu9hNrCH2eY9YPYKWLfT3uzpwcW5lIUl6d49zcaMa/nJ+/Nhvf3q06v0+WWR69J8U7T0+MjFOZcKDy7OuVR4cHHOpcKDi3MuFR5cnKsRSbtKmh+lYl0m6XFJvm+Rcy625cBpZtaTkH1/W6DT7lvkwcW5GjGzv5vZ1Oj7xYSE7TH3LdJQSQ9IWiLpfUnnSOpeonwu2h65tdfuceou1vCZ6JxrRJK2AbYHYuxbpGr2FToJ6Fd07Bzgi4QnpePU3YIHF+dqLNpS5kngXjO7J0ZVYwjbfBxhZgsJ6V37ARMkTYyOtWBmLxb1ZW1gF+BWM1sZp+5iflnkXAIknRpdWhTsW6S9ov20ni04thbwEvAPM0to3yKLs2/RQYStbm5Oum4PLs4lwMwuJSS0PgNA0kDCvlQf0TI15fOErIm7FNehoEd7r4JTvkBRLlsze4ewrU2pvdMLHQ28BzxSdDx23R5cnEvOeKCvpLHATEJuhWFmtgJA0hjCCtFAYEm0M+gdBedXmqB7AK1n4V8QfVaSpHUIGehutTUfMoxVN/ici3OJMbMrJJ0DnA+sBkaa2QcFn18NXF2iinyC7oqabeWY2jhe7BBgXda8JEqibg8uziVsDrA+MNXMOvq+RUcDr7fRz9h7IvllkXMJkTSVsLy8DKjTvkUqa18hSf0Jk7ZtjVqqrjvPg4tzCZB0JvA9YDJhP/N1C1eOojIXRKtHS6NHAJZKOqSgSJX7FqmafYUOB3rSdnCJU3dQ731j/OWver5IZN8iDo3qKdi3iHnAoqJyOSBX8P5l4OEY7Q4gXIbNAA4ATgAWAUX7FjGasJf6oIJjfwKejVt3qZePXJyLQdJQ4E7Crp8F+xbRBPSRNKHg2LbAw9F5gwhbDt9bbdtmVtW+RZI+F513SwJ1t8nTXLourbb7FmklYQRhQC/gHjPzBN3OuepJGgx0M7NeZtYb+AowStJm9e1Zejy4OFcbhwP/yr8xs/ujbz24OOdi2YcwLwOApEuA5Wb2RP26lC6/ic652tgB2FLSUsLduwuBr9a3S+nyCV3Xpfm+RenxyyLnXCo8uDjnUuHBxTmXCg8uzrlU+GqR69QkDSPcY1LKeEnL2/hsmpk9k3C3ugQPLq6z+xYhw36pLTHGtXF8NSFTvgeXKvhStOvUJG0AvAOsU8XpS4EtrSCbnCufz7m4Ts3M5gGXExI4VWI58DsPLNXzkYvr9KocvfioJSYfubhOr4rRi49aEuAjF9clVDh68VFLAnzk4rqECkYvPmpJiI9cXJdR5ujFRy0J8ZGL6zLKGL34qCVBPnJxXUo7oxcftSTIRy6uSykxevFRS8J85OK6nDZGLz5qSZiPXFyX08roxUctKfCRi+uSikYvPmpJgY9cXJdUMHoxfNSSCh+5uC5L0vrAFODHHlyS58HFOZcKvyxyzqXCg4tzLhUeXJxzqfDg4pxLhQcX51wqPLg451LhwcU5lwoPLs65VHhwcc6lwoOLcy4VHlycc6nw4OKcS4UHF+dcKjy4OOdS4cHFOZcKDy7OuVR4cHHOpcKDi3MuFf8P6zaEpgjlkTYAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 360x360 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot_flowsheet()\n",
    "XA1=1-XB1\n",
    "plt.text(xs-la-hla,ys+ly/2-14, r\"$x^1_A$\"+'={:.2f}'.format(XA1), color='r',size=15)\n",
    "plt.text(0,95, r\"$x^1_A+x^1_B=1$\", color='r',size=12)\n",
    "XA2=1-XB2\n",
    "plt.text(xs+lx/2+5,ys-14, r\"$x^2_A$\"+'={:.2f}'.format(XA2), color='b',size=15)\n",
    "plt.text(0,90, r\"$x^2_A+x^2_B=1$\", color='b',size=12)\n",
    "FT2=FT1-FT3\n",
    "plt.text(xs+lx/2+5,ys-6, r\"$F^2_T$\"+'={}'.format(FT2), color='g',size=15)\n",
    "plt.text(0,85, r\"$F^1_T=F^2_T+F^3_T$\", color='g',size=12)\n",
    "XA3=(FT1*XA1-FT2*XA2)/FT3\n",
    "plt.text(xs+lx/2+5,ys+ly+25-14, r\"$x^3_A$\"+'={:.2f}'.format(XA3), color='c',size=15)\n",
    "plt.text(0,80, r\"$F^1_T x^1_A=F^2_T x^2_A+F^3_T x^3_A$\", color='c',size=12)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "id": "0dd407ed",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAARcAAAEvCAYAAABv+oUUAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAtuUlEQVR4nO2deZhU5ZWH3wMoIIuADi6IEGUwohIijujgclFjVKKJxmgmKiFiFJeIScYRZsTuMhoVoyZBBXmiIkzinhhxCUHtGzWZuEIUcI8bihugyCoNZ/74btHVRVfVrbr31tbnfZ77VNd3v+X0dup8y/0dUVUMwzDipkOlDTAMoz4x52IYRiKYczEMIxHMuRiGkQjmXAzDSARzLoZhJEI05yKyLSJPI7IKkb1jsqn6xjQMo2iiRi5rgNHAPTHYUs1jGoZRJNGci+oGVD+OyZbqHdMwjKKxNRcAkfMQeRaR9YjMrLQ5hlEPdCr7iCKNAKg2ln3s3LwPXAZ8HehaYVsMoy4IF7mITEHkDxnvr0bkUUS2SsSqco+n+ntU7wOWJdK/YbRDwkYuVwFvIDIMOAA4CjgI1Q2IPAQMA/ZA5CZUZ8ZgV+7xgIJjijwAHJSj7ydR/UYMNhqGkYdwzkV1GSK/BGYB2+L+0T8L7h0Tu1X5xgszpjkPw6g4xSzozgf2ASah+m5Ro4g8gMiniHwKTAQmbn7voox4xzMMo+JIKD0XkX2AucCfgV1QPSLr/l7AC8AAVJcU6KsRyL+gm288kbuBnYN364FjUF2X1f5h4OAcvT+B6tE5xr0sGG9s3u/BMIyCFJ4WifQD5gDjgUeANxHxUPUzal0EzAb2BPI7l+jjDQaGo9qMyH3AUODpVn3kch65x+yE+1l0BDoi0gVoRrW55O/DMNo5+adFIj2Bh4BrUb0f1TXA1cDlGXWGAktxkcaekawpNJ7I1kCHwLF0B3YCXo00puNiYC1uynZq8PXFMfRrGO2WcNOivD3IrcBPgf7AeFTPjsGuXGN9BRfNLAL2AC5CdVZi4xmGUTLRnIvbKr4L+CfpKYXqYbFY1vZ4pwB7oHoJItsCf0c1WrRkGEYiRD2hOxEYgeoKAESezl89MnsDzwRf9wY+y1PXMIwKUrpzERkOrN3sWBzrENkO1aROuu4FHInIucAm3KKvYRhVSKRpkQgHAtcCX+CezxmjyoaYbKuaMQ3DKJ6oT0W/DRymyqG4dZdvRjepKsc0DKNIIq25qPJ+xttm3FQlUSoxpmEYxROLnosIXwKOBnId5Y+dOMcU4TwRnhVhvQgzIxtnGEZ0PRcRegK3Aaep8kWI+o0Aqu61HGOGwPRcDCNmQkUuIkwR4Q8Z768W4VERtgJuBxpVeSUuo/KNJ0KnuMdU5feq3IfpuRhGbISdFl0FjBJhmAjjcfoqJwDfBUYAl4jgi3ByTHa1OV6wK/QfhcYU4QERPs1xlW3qZhjtmVDTIlWWifBLMvRVVPkM97Di7LiNyjMeqoXHVMX0XAyjwpSk56JKUfoqmZEEgZ5LiEii5PEMw6g8oQ7RidBKX0WVI7Lub9ZzUc0vuRBmQTffeCJsoeeiyrqs9nn1XFRpU5JBhMuC8cbm+x4MwyhMwWmRCFvoq4jgqeJnVItNzyXEeIOB4ao0i3Afbei55HIeecZspeciQhegWRXTczGMEsk7LQq2fB8CrlXlflW20HMRITY9l0LjibA10CFwLKbnYhhVTGQ9FxFa6bmokpieiwhb6LmoYnouhlGFRH1wcRhZei6qJKbnIsIpwB6qXCLCtsDfVSOq3xntFhHZHfgL8DVVfanS9tQbUY//TwRGqHKUKl8DusdgUz72BhYEX5ueixGVy3BT6ysqbUg9UnLkIsJwun3wYy7caTAwBDiARr0ROF41mZOuMvjBJzj8f76KdtjA6r7rue/Wb+jnOz2bxFibx0zJDsAfgA3ARuAUbdClSY5pJE8QtSwEuuDW2IZb9BIvJUcuqjzHhTv9ABgN3BOUHZKUYwHglG+czI7/+Bed/nxvBs29mf/ceWBiY7XwCXCQNuihuEN948owppE8l9GyW7o1Fr3ETjTJhQbdAHwsKYnJnILjlV/ioUE3ZrztgVtMNmqYIGr5Fi1//x2BI0VkT4te4iOUc5GU9AaWA6szij8DdtGGqOkDih9PUpKWW7i8rfYJjDcMuAnoBRwZx5hGRcmMWtKko5dvld2aOiVs5DIMWK4Nul2CtoQaT1LSIrfQoHnlFiTlsjtqQ57sjgXGC9ovAEZISk4CJmHavTVLG1FLGoteYibsmsswYHGCdoQaT1LSIrfQoLFJPOQaLxizc8bbz4A1MY5rlJ+2opY0tvYSI2Ejl6+S+5/vIdw/5x6Skpu0QWdm3X8AOChHv09qg7b1BHOu8VrkFlJyCTBNG/TOUN9BfnJ+f8C+kpKrcDtF64DTYxjPqAB5opY0Fr3ESLgHF1PyAjAINqu+vasNuk9iRsU0XthpUbm/P6MyiMjtwInk/1DdCDygqt8qi1F1TEHnEkwLVgEHaoMmeqYkjvGyIqUuwWv6qektIqVyf39GZcg611IIO/cSA2GmRXvjtnxfyCyUlPwI+DYwALcO8THwgDboL7Lq5Zc/aNDsJ5jbHC/smJnOI2TkEmk8o2bIt9aSje0cxUCYH/ZXgYXZOzPaoFOBqZKSm4BZ2qB/batxG86jpPGKGbPKxzPKTIi1lmxs7SUGwuwWDQPyTRf2It6dpELjxT1mucczyk8xUUsa2zmKSHTJhZQs1AbdOyZ7qnLMSnyPRjwUudaSja29RCDS8X9JST9olQExcco9ZiW+RyNWzsRNc3I9Qb9tnntdgXOAHyVgV90TNSna3pT/WZtyj1mJ79GIj2uBJ/LcnwOcAa11mDN4PnaL2glRnUsl1iLKPaatt9QwqvoheVL+igjAg6q6tmxGtROiKdH5/rbAPAI9F/W8hUm2KxXx/QNxn2Bf4KY4Y9TzNsTdJqKNW+rGeF5B3ZhS2xkOEVFgG3Mu8RNViW4NGXouZWhXKm8Dh6nnHYqT5PxmQm2i4HRj3HjF6MaU2s4wEiWanov7JP9YfL8s7UpFPa9oHZhS2kRBPa8k3ZhS2xnlR0Q+BdIPwi4DRqjqe5WzKFnC6bn4fm69E8+LXc+lFMLYKL5ftA5MKW1KtRP4CiXoxojvDyulnVF29lHVdwFE5DncE/6HVNak5ChOz8XzyqXnUgrDyGOj+H6LDozn5dWBKaWN+H4jgHpeYwQ7FwAjxPeL0o1RzyupnREvIvI87sR3Nh+rat8Mx9IRt829KuJ4g4ALgQNwu5pPqKqXo+4QYCpwIPAp8BsgpdqitBimTjEU41yK2jER358C/Kt63vHB+6uBfYGjEloYHUYuWQjfb9GB8bxXMspz2ghoW22SslN8v7N63vrgbSvdmAJ2dsjVzig7A4CVwE+yylv+5kQ+ArYP6o2KON5ewDHA33EnittERHrj8n0txq0d7g5cg1tzvThsnWKJrufi+y16Lr5/k3rezODWVcAbQch+AO4f4aC0Y8nTLt1vfh0Yz8vWgcmnydKiA+P7TgfG8+7MZ6P4/mk52kQll537iu/n0o3JZ+eBedoZ5aUXMF9Vb85VQVX7ishWwJPALbiNjVKZo6p/BBCRe3BOqy3G4yKlE1R1JTBPRHoCjSIyJSgLU6cowum5+P6WeieeV1DvJJgqnIA7BXmQet67xRoYlkrbGHZaVGk7jdbEtRUdTCkWAbeqakEHLyLHAXepaimPJbTV3z3A9m1Ni0TkceB9Vf1uRtmuuB3R41R1Tpg6xdpUOBG973fG5YA+UD2vWL2T+UAD7uxFko6lIjZmRVddgrILgvdbRFe18LM0SiYdgTwt0koa9QtVVREZCHRT1fRu3gTgw3Qlcaf5OhYaRFWbS7Dty8BjWf28IyJrgntzQtYpitL1XHy/ba0Tz/tFcH8fYBpuQfR04HdZ7fcK+hygnrcke1Dx/fw6MJ6XKeWQW5Mlj52FbCxkZ6bzCBm5VORnaZSF9IfMtOBKMwM4C/e7fUhEOgCC+x0flVHv+8CtIcYpJY9Pb9wCbTYrgnth6xRFeD2XrN0S9TyndeL7TuvE8zZrnYjv98N5uvG4RaI3xfc99Tw/o4uLgNm4T/It/iGynEdJNuazM6SNBe0skor8LI2yMAR3Sjr7EOMjAKr6F6BbnvZzgH9LxjTAbVBkI1nlYeqEJoxzGUYRei7B9u1DwLXqefcHZVfjzomMDN4PBZbitl73xD0KEIVCNrayM4yNFbKzGn6WRmnsDHygqrNLbL+c5HKfr8AtNmezLS3RSpg6xaGqkS6amhaW0OZWmpr60NT0FZqapkW1oV7srAUb6+3CfSp3jdhHx6CfeRH6GBv0kffK0/4ewM9x73Hg9qyy/kGfx4atU+wVTc/FhexFaZ0E26kjcesGHQmxiBWVWrCzFmw0cvK14PXpCH0kOS16GLhQRHqo6udB2ck4May/FFGnKCqh5zIRGKGetwJAfD/KLyQstWBnLdhotM0RweufS+1AVZfhnjcKjYhsgztEB9AP6CkiJwbvH1LV9IHK6cD5wO9F5CpgN6ARuFZbzq+EqVMUkZ6K3uEDht/zbb7li7/KF7+gDKT4/nBgbdMoNvniP+2Lv2rPxXQQ30/0sYJjHuSYu0/kOF/8v/ji3+6Lv1UhOw96gu5No3go3abTBtYnaedubzDi7hM5NhjvMV/8nQrZiPtZbu2L/zdf/L9MG8/gnxznD0nKRiMn+wWvj5d53L7A3cF1AG5ROf2+b7qSqq4ADsdFtnOAFHAd7mhD6DrFEknPJfgn7QVcDfzC03C6LKW2KxVf/J2BFZ56a33xLwfme+rllXsopU1EGzsC6qm3yRd/LLCLp95lSbUzHKbnkhyRpkWeOukEX/yytCsVT4uXTyilTRQ8LU06odR2hpE0oZyLL7llAjytDsmFMDb6Urx8QiltSrWTEiUXfDHJBaP6KEpywdPql1zIZaMvLfIJnoaTXCimjS/uhK6n4SQXcti5ABjhS3HSCZ56JbUzjCRJTHLBFycT4KmTCfClRSYgmBbFzTBy2OhLi+SCpy3yCflsJENyIbNNUnb64nf2tG3phAJ2dsjVzjAqSWTJBV9apBN88W/ytLXkQhCyb5YJSDuWPO3S/eaVXPC0NMkFX5x8gqctkgtt2ehLi+RCVpuo5JRc8CW/5EIOOw/M084wKkao3SJftpQJ8LSwTEAwVdgsE+Bpck/zVtrGsNOiSttptMZ2i5KjYOTiS4tMgKelywQk7FgqYmNWdNUlKLsgeL9FdFULP0vDiIuSJRd8aVsmwFMnE+BLbpkAX/y7cQ96AawHjvHUa5Xxzpf8kguehpNcyGdnPhvD2JnpPEJGLhX5WRpGJQgtuZC9W+KpkwnwxckEeNoiE+DLljIBvviep5tlAgYDwz31mn3x7wOGkvVcRpbzKMnGfHaGsDGUnUVSkZ+lYVSC2CUXgu3bh4BrPXUyAcHuxuXASF/8rXE7HM2++N2BnYBXSzM/tI2t7CxkY/C+EnZWw8/SSBCZMeNTunRxSnXLly/jxRdH6M0312XuokjH/wF88Rd66hV8riij/ldwn8CLgD2Aizz1ZkUyIty4VW9nLdhYb5R7QVfGjeuvN9/sUoxcfvlzNDev1oaGusxdFOn4fxCyFyUTgFt3mOapd4kv/ra4tAhJ/9NWvZ21YKMRnc2OZezYjuy5Z1eam6PlLvL9tnMNtc7Emavtt4BLcR9M7wNT1fOujav/qLmiS5EJ2Bt3EhWcNmdS6lvZY1a7nbVgoxEDcsMNH3H88Rvo1Wtn3n77pJL7cdk7H8Ed+PwmzlH8FPdEc6G2I4Hf49bnjsWlObkqQ2A+Uv8QXc+l1RpBEW2O9MU/F7dzUo6j6rVgZy3YaMSAnntuX/nBD7Zit92epF+/KLmLWnINeZ7LNeSkURvF96cEZbm4BJeh4ozg/Z8DZ3KJ+P6Ngc5zlP6jy1zaZVctX8Qgc6mqcN5559PUpEyYMHlz2ZlnHswf/rCJyy9f0Gab8eOPY9asdSWP2dT0OE1Nd2SV7UpTk9LUlFeakqamD2lquiSr7Kig7aFR+1eNKHNpGIZDp079tUyf3sjee18I/EzGjevPqFHzWLLkY1577d8AZNy4gYh009/8xk1/+/adwCeftOQu8v1wuYs8L527aMtcQ573jvh+mFxDXWg5JZ4m/Yzanjhpyyj9R15zMQwjzaJFkxk0qIdMmDCR4cPns3Gj8vzzQ/XWW92Dup06DWDkyKdl1qy1Mnv2Onr1GsyiRcdk9PB9XHqSQleaKLmGXmdLzd79g9c+MfRvkYthxIVOnXqDTJt2KUcffQWrV29i7twResstmyMTvemmJHIXlZpraDowTXz/h7jMAfvjFmvBPQQbtX9zLoYRK599tpQuXfrQ1DRbZ8wo9vmxYnMXRck1dAtOnGwaLivkGlxyvam0pJmNlMvIpkWGERMyceJs9t13L95+ex2DB3+7hC6KnRa9jFv7aLHB9/vjoqOX8w2knrdRPe884F9wj4zsgDsnRcZryf2DORfDiAWZMOEiDj/8VObNm8Ezz5zH7rt3lwsumNxm3UmT7hTfVznrrEOzbqWnRYWuNA8DXxff75FRVlSuIfW8Fep5L6rnrQLOAf6mnpd2HJH6j3z83zBqmTiO/8vZZx/H6NF/5LnnntSGhoMBZNq0ZXTu3FlPP717q7rjxvXjyCPfYs2aL3jhhV/qddf9T8njunMpi4GFOEGx3YBrgV+q512cUW8Mbhq0u3re20HZATi5kAVAT5yg2teBg9TzXiim/1xY5GIYEZAzzhjCqFH38vrrS3jzzZZnhBYvbmC33brJBRc0tmowYMBdvPTSX1myZCE9e46IMnaQDC9MrqEOQR3JKNuAi0LuA2YC2wAj046lyP7bxCIXo11TzgcX5cwz9+Pww/+PuXN3YbvtUuy882j98Y/7Jz1upbDIxTDKxeDBdzB//h16yy0fsnq1zw479C3cqHaxyMVo15QrcpFzzjmZb37zDlau3IgqiEDv3h353e966S231OUDpxa5GEY52GefG5k790r9znc66UknddLvfKcTH3+8ga5ds7NY1A3mXAwjYeSCCybTuXNnvfbaSa1ufPLJcrp1O6xCZiWOTYuMdo2lFkkOi1wMw0gEcy6GYSSCORfDMBLBnIthGIlgkguGUUZmyIxPu+DyFi1n+bIXeXHEzVqfeYsscjGMMvIUT+0zRsd0HaNjuq5m9Ye7suvtlbYpKSxyMYwycrO6vEVjZWzHPdmzazPR8hb50nZeIU/z5xXyxR8EXAgcgEtR84SnnhdX/2DOxTDKzg1yw0fHc/z2H/DByqd4alSp/fiyOa/QYlxeod2Ba3AzkkKSCHsBx+CEobZOoH+bFhlGuTlXz+17H/d1/oiPXulHv1sidLU5r5Cn3jxPvek4WYSfBHnG8zHHU6+/p953yJ2ML0r/5lwMIw5+JD863xdfL5AW9bmz5KyD75P7Nv1cfr4gu/6teuuG93n/8sEMPjzCsEcDcz1tlZzsDpxDyFa5a4Wn3qYk+webFhlGLEzVqb+eLtMb98blLRon4/qPYtS8JSz5+DVc3qJxMm6gIN1+oy5vUV/6TviElrxFvoTLW+Rp7rxFnnrv+BIur1AIIvVvkYthxMQiFk0exKAeE2TCxOEMn7+Rjfo8zw+9VV3eok50GjCSkU/PkllrZ8vsdb3oNXgRFctbFAbLW2QY1cBUnXrDNJl26dEcfcVqVm+ay9wRt2hL3qKbtKryFiXevzkXw4iRz/hsaRe69GmiafYMreq8RYn3b9Miw4iJiTJx9r7su9fbvL1uMJXJW+RL+LxCIYjUvzkXw4iBCTLhosM5/NR5zJvxDM+ctzu7d8/cOQL4ifzkyvvl/k23yW1rfyu/XT9LZq09W84+NqNKSXmLfCk9b1EBIvVvYlFGuyYOsaiz5ezjRjP6j8/x3JMN6vIWTZNpyzrTufPp2pK3qFEafYBGbfQApsiUl9ey9qMGbTikrX4LERxyazOvkKcteYV8aclb5KnLW+SLvw3uEB24HNE9aUkZ8pCn3pqw/efCIhfDiMAZcsaQUYy693VeX/ImLXmLFrO4YTd263aBtOQt6kOfPT7js8cBxsm4AduxXb/lLH+o1LE9jZS3qC9wd3AdAAzJeN+3yP7bxCIXo11TTpnLu+Xu5vWsb97EJt2RHbs8zuMPXqaXmUC3YRilM07GDexM5w6n6qldxuiYrvdy79cO5MDR42Rcv0rblhTmXAyjDHSj2/Ef8uHn6fc36U2PKEpHOppzMQyjdHrQ49BlLFuSfv9T+el1K1ixfobOeLqSdiWJHaIzjDLQi157D2LQl26T29ZuYtOmAQxY+SRPfv00Tqu0aYlhC7pGu8byFiWHTYsMw0gEcy6GYSSCORfDMBLBnIthGIlgC7oGAAMnPngm8L1K21Fu1r3z4qGddxnyhHToGEb2sZ743VtXjp6R5AAWuRhpvgcMq7QR5War7Xd9AenQ3hzLMMrwQWLnXIxMFrx15Wiv0kYYyTJw4oN+OcaxyMUwjEQw52IYRiKYczEMIxHMuRiGkQjmXAzDSARzLoZhJII5F8MwEsGci2EYiWDOxTCMRDDnYhhGIphzMQwjEcy5GIaRCOZcDMNIBHMuhmEkgjkXwzASwZyLYRiJYM7FMIxEMOdiGEYimHMxDCMRzLkYhpEI5lzqCBHZTUSk0nYYBphzqRtEZCjwBvCiiBxlTsaoNOZc6oetgJXAXsA9mJMxKow5l/oinT6zG+ZkjApjzqW+MSdjVAxzLu0DczJG2THn0r4wJ2OUDXMu7ZMtnIyqFmhiGMVRlYnoRWQg8CpuB8RIjrSTufeLpa9q5533eLbSBhn1Q1U6F2AZ8CCwa6UNqSEGANuV0G4D0EG27rIiZnuMdk5VOhdV/Rw4vtJ21BIiMhx4FNg2ZJMNwEbgt0Bq6+0HzE7KNqN9Ymsu7Y8NwDpgFjBYVc9Q1XcrbJNRh1Rl5GIkQqtIxRyKkTR1EbmIyEIR2SgituWxJRapGBWhLpwLMB0YWmkjqgxzKkZFqQnnIiLPi4i2cX0EoKrXq+qiGMc7UUT+JiLLRGSdiLwiIheLyNZt1B0iIo+KyBoReV9ELhWRjsXWiYlOmFMxqoRaWXMZgHvi9ydZ5a8kNN52QBNwNfApsD/QCOwInJeuJCK9gUeAxcA3gd2Ba3BO++KwdWLiU2ATzqnYmopRcWrFufQC5qvqzeUYTFVvyipqEpGewLki8iNtOc46HugKnKCqK4F5Qb1GEZkSlIWpE4fNb4hIb1XdGEd/hhGVqp8WicgQnJ0vVNiUZUD2tOhoYG6Wg7gD50wOLaJOLJhjMaqJqncuwOjg9WkR6ZxxhXrgThydCl052nYUkW1E5CDgfGCatn4I58vAy5ltVPUdYE1wL2wdw6g7asG5HBS8TsMtVKav6ekKIvKyiDQHXzeLSOY/8/dxi5yFrrZYHVxPAH8BLsy63xu31pHNiuBe2DqGUXfUwprLENw//7is8kfSX6hqvghgDvBvJY7978A2uAXdS4DrgXOy6rR1tkayysPUMYy6ohacy87AB6pa6rMvy4HPSmmoqs8HXz4pIp8At4nINar6RlC+ArfYnM22tEQrYeoYRt1R1dOi4CzINkTbco4yLcok7Wi+lFH2MlnrJiLSHydl8HIRdQyj7qj2yOVrwevTEfqIMi3KZGTw+mZG2cPAhSLSI3iSG+BkYC1ujSZsHcOoO6rduRwRvP651A5UdRluGzk0IvIn3JrOItzDfiOBnwJ3ZkyJwC0qnw/8XkSuAnbDHba7NmPrOUwdw6g7qt257Be8Pl7mcZ8BxgIDgWbgn8AkMnaoAFR1hYgcjlvonYNbQ7kO5zxC1zGMeqSqnYuqehUadzIwOWTdxcBhUesYRr1R1Qu6hmHULuZcDMNIBHMuhmEkgjkXwzASwZyLYRiJYM7FMIxEMOdiGEYimHMxDCMRzLkYhpEIdeNcHhRZuFRkI5a7yDCqgrpxLm/A9Fstd5FhVA1V/WxRMZyvej0A4aR1C+OEwacCB+IeNvwNkKKQCLbIWODWNu6cjer0jHql9Z8chwIMnPigX6HxK8KG5e8N7NR757dCSjLXC8OABUkPUjfOJVbiyTV0GE6zJc0/Y+7fiIGNq5YP6NRrx3eQjpsqbUsZWQD8LulBpLWYfR0goqhG+xgSmQT8FzCAtOaKyH+RToyWT4elJXLpgeqq2Ps3YiXIL76Nqq4tWNkoiqpfc/m1yPmI6A0imyUQrhI5+GORTXNEFiQ07NHAXJLLNZR0/4ZRcareuZyv+utnYMXIIK3HZSL9T4B5r8LHz4eVrxQRRDoVvFrYItcQxecaegORZkReQeSsrHtx9G/UEJKS/pKSRyUlL0lKFklKpkiqvhd6qt65APwNJg+DHteLTDwC5jeDPgZDG1Q3C2v/SeTl94PcRe+LNP8pWu6iKLmGluKEpk4DjgWeAqYj8uOY+jdqkcfoy0x60Aj8DOETvssmTqi0WUlSEwu6E1RveErk0h/AFSth00wYMVn1w8w6R2XkLto5uDIoRaS7tFxDqnOBuRklDyPSGbgYkV+hml44tFxGdYak5DrggjZuPcTjnAGco6rPisjWvM+bLGc0cG+E8b6D+xAbjktV8wrwC23Q29uo2+bupDa07E6GqVMMNeFcAD6Bpd2gz10we5Lqs0U2LzZ3Udy5hu4BTsJp8v4zgf6NauAjTmAtMC+rfD1dVXUpLqqFRnqwlp78gQ+zuyiSn+CyUfwY+AQ4BvidpGR7bdCp6UqSKrw7GaZOsdSEc5kpMvs/YK+XYN1w+DZOPLsYvk/bZ0+ySc+Bt8g1RDy5htJRSVL9G5WkBzuyhPdYwolZdz5IfyEp6Uwz9/E3NvIqsyKOeKw26CcZ7x+TlOyMczpTM8rH4zYLTtAGXQnMk5T0BBolJVOCsjB1iqLq11yuF7noe3Dq/8KMuXDeUOieuXOUySyROxHRq0Syd1zS06JCV5qHga8j0iOjLEquoW/jPlneTqh/o8JISvrSla1ZzSJV/XvW9VZQpyObuJ2FDOQJfqaqL0UZM8uxpJkP9M0qOxqYm+Ugsncnw9QpiqqOXH4hctwP4cr74clxqmcBPCUy5RC4CPhZZt3LRPqNgRMWwpqecCSZ/6TF5y7anGuIrFxDrbaPRcYAtwC7o/p2UHYvLonbC0BHnNM4GTg/Y70lXP9G7bCUQ9gJWM1zkmq187hRG4LDZMoMXmdf7uNuVb0ms3mwc9Sx0DDaoM0Fqvw7bmqTyZeBx7L6eUdSkt6dnBOyTlFUbeTyc5EhJ8G982HJQjgkXf4UNOwD3W4Uacysvyfc9RT89TVY2AdGRBpcdQVwOO6XPQdI4XINNWTV7BDUydxSfAU4HbdQdzcwBBiDtsyBi+jfqBXWczgARzGJzB1I5TQASclIhNPZgZ40cJikZIGk5PyMHiKnHZaUHI5bL7kh61aY3cnYdzCrNnL5b5frZ6tdAS+jPHiG6PpzMsquENnvFNj/NthlJ0jtDqMjGxAm15DqTGBmVtl/A/8dS/9G7dCJ4Sxjy72fNXShkcwUeO+1ut/ycRIp7bCkZCDuSP8ftUFntlElzO5krDuYVetcimE43PEY3DFZ9cNrRPxd4QeVtsloZ2zDbnzEKt5nVNadVwBU9a+0jnCzKXZHczOSkj64dbx3gFPbqBJmdzL2Hcyady7XiJx8Duy+DAa+K/IfJwM7QMdLRba9RLWkX5ZhFIOkpAM96M0bLNbij0mkKXZHMz32NsADwNbAaG3Q1W202WJ3UlJb7E6GqVMUNe9cDoYbb4Mrx6tOSpe9KfJFd/gG8NsKmma0F75gMFvTgZX8o63bkpKxtByu2xN3NmUdsBE4QBt0AyVMi4KF47uBfwVGaoN+lKPqw8CFkpIe2qCfB2XZu5Nh6hRFTT8VfYPI5EPgon1Uu2eW/1Xkg1fgwdNVx1XKNqM2iOOpaDldzmcAv+IuztJFOiNnvZR0AD4HdtUGLWb3Mld/M4AfAhNwO5SZzNcGXR/U643bQVoIpHcnrwV+qQ16cdg6xVLTkcu5qj8ja0saYKTqjiMrYI/RTlEO4gvgFR4pUPNfgRVxOJaAI4PXX7Vx70vAWwDaoCuCnaTrcRHSp7jdycZ05TB1iqWmIxfDiEoskYtICvgp0EPz/ENJSk4Cvq8NGn03swao2nMuhlErqGqDqnbP51gCvoI7XNkuMOdiGOVjKLS96FuPmHMxjPLxFcy5GIYRJ8FuzL8Ar1balnJRNwu6IrIQd4agg0YV6DbaDSbQnRz1FLlMx5KiGUbVUDfORVWvV9VFcfUnIkNE5FERWSMi74vIpSJS8JF4ERkrItrGNT6O/g2jVqjpQ3RJIQknRYupf8Ooasy5tE2L5J8Gkn8SSP6JTNFwgk7PaK6kaPH0bxhVTdVPi0Tk/GBaMTmj7GAR2SQJJ0XThJOiJdi/YVScqncuqvprnNbEhQDihKznAR8T8ilScXQqdGU02SJpmZaQFE1EmkXkFQmRFK2E/o0aQ+TK/UWeWS7y+nqR19aJ/OEpkca63dmseucSMBnoISITcQLECgzVjKRoIvKyBEnRgn/qqkqKJpYUzUDXwxMTVAd1htu3g133gJ2mVNqqpKiJNRdVvUFELgWuADYBIzQrKZpmJEVrg7IlRdM2kqJJkBRNRH6llhSt3aI66R8EJ3RVL1ktMudt6DYoSp8itJ3ITMmbyEyEE3EpSPbACUK9DcwGpqjyRdT+oUacS8BSoA8wuwS1L0uKZlQVIpcPgh/sBb8trLecsw+i7DpuBzQBV+P+5vbHySvsCJwXQ/+gqlV/4Tyq4rZ2Py+h/digfd4ro/7jwO1ZffQP6h1bwvjfCdp+KYn+7Yr0t6VA1/KOmeoBz6yA/70/Wj86CXQFaM+Msv8CXZNZVkR/l4N+Cipx9F/1ay4ichFOdHgGzqN2l6ykaCJyZbB7tFZE1gevx2ZUKSkpmiScFC3G/o0KI/Lr80VQkRsydjWvOljk400icxa0lKW2gkNfgvf+qXrKcRGHDXYdiWvXcRlOizee/iv9yVHgU+U43CfLExlly4BVWfV8wM94/zLweIRxe+OmYfOAI4AzgVXAZVn1xgDNwICMsntxSduOxun4pqOuHxXbv11l+RuLLXKBp5fD/JXu65/1h1fXwZMfQuNWLXUefhn+9Bo0yJbtVUA7Fboy6n8E2thGP6tBLwxns3YE3Qb0INC3QK+Jq/+K/3Lz/NKH4HZw3iV4wDIoPy/4g2jMKFsKXBp8PQCnUzoxhvEfw0UTS3Fymh2z6owNbBmYUfZzXDqJNUHb54DTSunfrrL8ncXoXH55LqjC1Inwf5/A4rVw6Q4Z98e7+y+vc/cWr4VZ97Tc17Hufv4ro/4G0Au2tEOXgP48nM26LqPv20A7xNV/XTwVHWxBN+P+ULoAD6rqNyprlVELxP1UtMhTy2DvPrByE8wcoTop9OaDCNvhtG/zosqzQf0NwH+qttbQFeE9YKYq/xNizH2BbXALupcAv1PlnDj6r6XdojYRkYE4mYUuwfsjcMfp+6nqe3kbG0bsfLIUuvWBu2YX41gCyr6rqcrzwZdPivAJcJsI16jyRtT+q35BNwTH46ZBAKhqWoG9X2XMMdorIjNnwxF7wUvrYPi3S+ii2MOeWyYyE6IkMks7mnT0FKn/enAuhwJL0m9E5Dpgvapm53ExjMQQuf4i+N6p8L8zYO55MLR75s5R67qz7nQ7S1dl77iUuKtJXLuO6Yw8b8bSf6UX1GJYkHsdl7luLbAatzh6aKXtsqs2LmJY0IWrj4NPFe7O2NX8+zJ4YdWWdX/WD97eAC+uhhsvjzau9gZdCjoP9AjQM0FXgWbtauoY0GbQARllfwL9T9CjQY8ETQVt7yi2/1xXzUcuqjpIVTuqaldV7aaqO6mqnRUxyoLIz4fASffC/CWw8JCWO081wD7dRG5sbN1iz7vgqb/Cawuhz4goY6uyAjgc6IiLelK4RGYNWVU7BHUyH5J8BrfbeTdwF+45uEm4Z+KK7b9N6mK3yDBKpZwauiJX7Aen/B/ctgvslILdR6uO6p/0uJWi5iMXw6gdht8Bj92hOvlD+MyHXftW2qIkscjFaNeUK3IRueZkOOcOWLbRLfMIsENHuKKX6iXFbD/XDBa5GEZZOPhGuO1K1V06qfbvpLpLJ1iyAbrX7WFPcy6GkTBuS7prZ9Xxk1rfeX859DqsMlYlj02LjHaNJUVLDotcDMNIBHMuhmEkgjkXwzASwZyLYRiJYM7FMMqEiOwvIssDKdZ1IvKUiFjeIsMwIrMemKCqnXHq+3sAdZu3yJyLYZQJVf2Hqs4Ovl6NE2yPmLdIhojIoyKyRkTeF5FLRaRjnvp+kB65revAKH1nU/NKdIZRi4jIIGAvIELeIiklr9A5QM+sskuBr+KelI7SdyvMuRhGmQlSyjwDPKSqD0boajwuzccJqroSJ+/aE2gUkSlBWStUdXGWLVsD+wF3qmpzlL6zsWmRYcSAiJwfTC0y8hbJwUE+rQUZZVsBLwH/VNWY8hZplLxFR+FS3dwed9/mXAwjBlT11zhB6wsBRKQ/Li/Vx7SWpnwRp5q4X3Yf4uhU6Mpo8mWytGxV9R1cWpt8udMz+S7wHvBEVnnkvs25GEZ8TAZ6iMhEYD5OW2Goqm4AEJHxuB2i/sCaIDPoPRntixXo7k3bKvwrgnt5EZFtcAp0d+qWDxlG6htszcUwYkNVbxCRS4ErgE3ACFX9MOP+dGB6ni7SAt1FDdtGmeQoz+ZYoDtbToni6Nuci2HEzFKgDzBbVas9b9F3gddz2Bk5J5JNiwwjJkRkNm57eR1QobxFEiqvkIhsi1u0zRW1lNx3GnMuhhEDInIRcCowA5fPvHvmzlFQ58pg92ht8AjAWhE5NqNKiXmLpJS8QscDncntXKL07ah03hi77KrkRSx5izgu6CcjbxHLgFVZ9XzAz3j/MvB4hHF746Zh84AjgDOBVUBW3iLG4HKpD8go+xOwIGrf+S6LXAwjAiIyBLgXl/UzI28RDUA3EWnMKNsDeDxoNwCXcvihUsdW1ZLyFonI9kG7O2LoOycmc2m0a8qbt0iacRGEAl2AB1XVBLoNwygdERkIdFDVLqraFfgaMFpE+lXWsuQw52IY5eF44PP0G1V9JPjSnIthGJE4FLcuA4CIXAesV9WnK2dSstghOsMoD3sDXxKRtbjTuyuBr1fWpGSxBV2jXWN5i5LDpkWGYSSCORfDMBLBnIthGIlgzsUwjESw3SKjrhGRobgzJvmYLCLrc9ybo6rPx2xWu8Cci1HvnIhT2M+XEmNSjvJNOKV8cy4lYFvRRl0jItsB7wDblNB8LfAlzVCTM8Jjay5GXaOqy4DrcQJOxbAe+I05ltKxyMWoe0qMXixqiYhFLkbdU0L0YlFLDFjkYrQLioxeLGqJAYtcjHZBEdGLRS0xYZGL0W4IGb1Y1BITFrkY7YYQ0YtFLTFikYvRrigQvVjUEiMWuRjtijzRi0UtMWORi9HuyBG9WNQSMxa5GO2ONqIXi1oSwCIXo12SFb1Y1JIAFrkY7ZKM6EWxqCURLHIx2i0i0geYCfzQnEv8mHMxDCMRbFpkGEYimHMxDCMRzLkYhpEI5lwMw0gEcy6GYSSCORfDMBLBnIthGIlgzsUwjEQw52IYRiKYczEMIxHMuRiGkQjmXAzDSARzLoZhJII5F8MwEsGci2EYiWDOxTCMRDDnYhhGIphzMQwjEf4fiWvz3Uc3mKAAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 360x360 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot_flowsheet()\n",
    "XA1=1-XB1\n",
    "plt.text(xs-la-hla,ys+ly/2-14, r\"$x^1_A$\"+'={:.2f}'.format(XA1), color='r',size=15)\n",
    "plt.text(0,95, r\"$x^1_A+x^1_B=1$\", color='r',size=12)\n",
    "XA2=1-XB2\n",
    "plt.text(xs+lx/2+5,ys-14, r\"$x^2_A$\"+'={:.2f}'.format(XA2), color='b',size=15)\n",
    "plt.text(0,90, r\"$x^2_A+x^2_B=1$\", color='b',size=12)\n",
    "FT2=FT1-FT3\n",
    "plt.text(xs+lx/2+5,ys-6, r\"$F^2_T$\"+'={}'.format(FT2), color='g',size=15)\n",
    "plt.text(0,85, r\"$F^1_T=F^2_T+F^3_T$\", color='g',size=12)\n",
    "XA3=(FT1*XA1-FT2*XA2)/FT3\n",
    "plt.text(xs+lx/2+5,ys+ly+25-14, r\"$x^3_A$\"+'={:.2f}'.format(XA3), color='c',size=15)\n",
    "plt.text(0,80, r\"$F^1_T x^1_A=F^2_T x^2_A+F^3_T x^3_A$\", color='c',size=12)\n",
    "XB3=(FT1*XB1-FT2*XB2)/FT3\n",
    "plt.text(xs+lx/2+5,ys+ly+25-22, r\"$x^3_B$\"+'={:.2f}'.format(XB3), color='m',size=15)\n",
    "plt.text(0,75, r\"$F^1_T x^1_B=F^2_T x^2_B+F^3_T x^3_B$\", color='m',size=12)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "5b0bf706",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.9.12"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}