Added(A)/Deleted(D) following books

 A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/Chapter3.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/Chapter4.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/Chapter5.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/Chapter7.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter10.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter12.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter13.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter14.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter15.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter16.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter17.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter19.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter20.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter21.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter22.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter23.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter24.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter25.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter26.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter27.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter28.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter29.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter30.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter31.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter32.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter33.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter34.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter35.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter36.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter6.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/chapter9.ipynb
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/screenshots/chapter32.png
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/screenshots/chapter331.png
A  A_Textbook_of_Applied_Electronics_by_R_S_Sedha/screenshots/chapter332.png
A  Chemical_Engineering_Thermodynamics___by_S._Sundaram/ch10_1.ipynb
A  Chemical_Engineering_Thermodynamics___by_S._Sundaram/ch11_1.ipynb
A  Chemical_Engineering_Thermodynamics___by_S._Sundaram/ch12_1.ipynb
A  Chemical_Engineering_Thermodynamics___by_S._Sundaram/ch13_1.ipynb
A  Chemical_Engineering_Thermodynamics___by_S._Sundaram/ch14_1.ipynb
A  Chemical_Engineering_Thermodynamics___by_S._Sundaram/ch15_1.ipynb
A  Chemical_Engineering_Thermodynamics___by_S._Sundaram/ch1_1.ipynb
A  Chemical_Engineering_Thermodynamics___by_S._Sundaram/ch2_1.ipynb
A  Chemical_Engineering_Thermodynamics___by_S._Sundaram/ch3_1.ipynb
A  Chemical_Engineering_Thermodynamics___by_S._Sundaram/ch4_1.ipynb
A  Chemical_Engineering_Thermodynamics___by_S._Sundaram/ch5_1.ipynb
A  Chemical_Engineering_Thermodynamics___by_S._Sundaram/ch7_1.ipynb
A  Chemical_Engineering_Thermodynamics___by_S._Sundaram/ch8_1.ipynb
A  Chemical_Engineering_Thermodynamics___by_S._Sundaram/ch9_1.ipynb
A  Chemical_Engineering_Thermodynamics___by_S._Sundaram/screenshots/TaVSN2_1.png
A  Chemical_Engineering_Thermodynamics___by_S._Sundaram/screenshots/TempVSMoleFraction13_1.png
A  Chemical_Engineering_Thermodynamics___by_S._Sundaram/screenshots/Tempvsequlibrum14_1.png

1 files changed, 631 insertions, 0 deletions

diff --git a/Chemical_Engineering_Thermodynamics___by_S._Sundaram/ch2_1.ipynb b/Chemical_Engineering_Thermodynamics___by_S._Sundaram/ch2_1.ipynb
new file mode 100644
index 00000000..f0ca2482
--- /dev/null
+++ b/Chemical_Engineering_Thermodynamics___by_S._Sundaram/ch2_1.ipynb
@@ -0,0 +1,631 @@
+{
+ "metadata": {
+  "name": "",
+  "signature": "sha256:3938cdb77cb1cf91b5da41037c46803e1323bc2d84e6ed7eb8e8d2a220b8bc11"
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+  {
+   "cells": [
+    {
+     "cell_type": "heading",
+     "level": 1,
+     "metadata": {},
+     "source": [
+      "Chapter 2 : P V T Relations"
+     ]
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 2.1  Page No : 29"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "\n",
+      "#Given\n",
+      "m = 140 #m is the mass of N2 in Kg\n",
+      "P = 4.052*(10**5) #P is the pressure of the system in Pa\n",
+      "V = 30 #V is the volume of the system in m**3\n",
+      "R = 8314.4 # R is the gas consmath.tant\n",
+      "\n",
+      "#To determine temperature required\n",
+      "T = P*V/float(((m/28)*R)) #T is the temperature of the system in K\n",
+      "print \"Temperature of the system is \",\n",
+      "print \"%.6f\"%T,\n",
+      "print \"K\"\n",
+      "#end\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Temperature of the system is  292.408352 K\n"
+       ]
+      }
+     ],
+     "prompt_number": 1
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 2.5  Page No : 33"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math\n",
+      "import numpy\n",
+      "from numpy.polynomial import Polynomial as poly\n",
+      "\n",
+      "\n",
+      "#Given\n",
+      "n = 1 #n is Kg moles of methane\n",
+      "T = 423 #T is the temperatue of the system in kelvin\n",
+      "P = float(100) #P is the pressure of the system in atm\n",
+      "Tc = 191 #Tc is the critical temperature of the system in K\n",
+      "Pc = 45.8 #Pc is the critical pressure of the system in atm\n",
+      "R = 0.08206 #R is the gas consmath.tant in (m**3 atm/Kg mole K)\n",
+      "\n",
+      "#To calculate the volume of methane\n",
+      "#(i)Umath.sing ideal gas equation\n",
+      "V1 = (n*R*T)/P #V1 is the volume of the gas in m**3\n",
+      "print \"i)Volume of the gas umath.sing ideal gas equation is \",\n",
+      "print \"%.6f\"%V1,\n",
+      "print \"cubic meter\"\n",
+      "\n",
+      "#(ii)Umath.sing Vander Waals' equation\n",
+      "a = (27*(R**2)*(Tc**2))/(64*Pc) #Vander Waais consmath.tant\n",
+      "b = (R*Tc)/(8*Pc) #Vander Waais consmath.tant\n",
+      "q = numpy.poly1d([1,-(((R*T)+(b*P))/P),-((a*b)/P)+(a/P),0])\n",
+      "r = numpy.roots(q)\n",
+      "print \" ii)Volume of the gas umath.sing Vander Waals equation is\",\n",
+      "print \"%.6f\"%r[0],\n",
+      "print \"cubic meter\"\n",
+      "\n",
+      "#(iii)Umath.sing generalized Z chart\n",
+      "Tr = T/Tc #Tr is the reduced temperatue\n",
+      "Pr = P/Pc #Pr is the reduced pressure\n",
+      "#From the figure A.2.2,\n",
+      "Z = 0.97 #Z is the compressibility factor\n",
+      "V = (Z*R*T)/P \n",
+      "print \" iii)Volume of the gas umath.sing Z chart is \",\n",
+      "print \"%.6f\"%V,\n",
+      "print \"cubic meter\"\n",
+      "\n",
+      "#(iv)Umath.sing molar polarisation method\n",
+      "#From Table 2.2\n",
+      "Pmc = 6.82 #Pmc is the molar polarisation for methane\n",
+      "#From figure A.2.4\n",
+      "Z0 = .965 \n",
+      "Z1 = 14.8*(10**-4)\n",
+      "Z = Z0+(Z1*Pmc)\n",
+      "V = (Z*R*T)/P\n",
+      "print \" iv)Volume of the gas umath.sing molar polarisation method is \",\n",
+      "print \"%.6f\"%V,\n",
+      "print \"cubic meter\"\n",
+      "\n",
+      "\n",
+      "#(v)From experiment\n",
+      "#Given\n",
+      "Z = 0.9848\n",
+      "V = (0.9848*n*R*T)/P\n",
+      "print \" v)Volume of the gas calculated by experimental Z value is \",\n",
+      "print \"%.6f\"%V,\n",
+      "print \"cubic meter\"\n",
+      "#end\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "i)Volume of the gas umath.sing ideal gas equation is  0.347114 cubic meter\n",
+        " ii)Volume of the gas umath.sing Vander Waals equation is"
+       ]
+      },
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        " 0.322787 cubic meter\n",
+        " iii)Volume of the gas umath.sing Z chart is  0.336700 cubic meter\n",
+        " iv)Volume of the gas umath.sing molar polarisation method is  0.338468 cubic meter\n",
+        " v)Volume of the gas calculated by experimental Z value is  0.341838 cubic meter\n"
+       ]
+      }
+     ],
+     "prompt_number": 2
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 2.6  Page No : 35"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "%matplotlib inline\n",
+      "import math\n",
+      "import matplotlib.pyplot as plt\n",
+      "import numpy\n",
+      "\n",
+      "#Given\n",
+      "P1 = 266\n",
+      "T1 = 473.16#Initial temperature in Kelvin\n",
+      "T2 = 273.16#Final temperature in Kelvin\n",
+      "V1 = 80\n",
+      "V2 = 80 #Initial & final volume in litres\n",
+      "N1 = (14.28/28)\n",
+      "N2 = (14.28/28) #Initial and final Kg moles are equal\n",
+      "Tc = 126 #Critical temperature of N2 in K\n",
+      "Pc = 33.5#Critical pressure of N2 in atm\n",
+      "\n",
+      "#To calculate the final pressure achieved\n",
+      "#(i)Umath.sing ideal gas law\n",
+      "p2 = (P1*V1*N2*T2)/(V2*N1*T1);\n",
+      "print \"i)Final pressure of N2  using ideal gas law is \",\n",
+      "print \"%.6f\"%p2,\n",
+      "print \"atm\"\n",
+      "\n",
+      "#(ii)Umath.sing generalized Z chart\n",
+      "Tr1 = T1/Tc#reduced initial temp in k\n",
+      "Pr1 = P1/Pc #reduced initial press in K\n",
+      "#From the Z-chart compressibility factor coressponding to the above Tr1 &Pr1 is\n",
+      "Z1 = 1.07\n",
+      "P2 = [125,135,150]\n",
+      "Z2 = [0.95, 0.96, 0.98]\n",
+      "F = [];\n",
+      "for i in range(0,3):\n",
+      "    F.append((P2[i]/(Z2[i]*T2))-(P1/(Z1*T1)));\n",
+      "plt.plot(P2,F)\n",
+      "plt.ylabel(\"F\")\n",
+      "plt.xlabel(\"P2\")\n",
+      "plt.title(\"P2 vs F\")\n",
+      "plt.show()\n",
+      "P3 = numpy.interp(0,F,P2);\n",
+      "print \" ii)Final pressure of N2 from Z chart is \",\n",
+      "print \"%.6f\"%P3,\n",
+      "print \" atm\"\n",
+      "#(iii)Umath.sing Pseudo reduced density chart\n",
+      "R = 0.082 #gas consmath.tant\n",
+      "v = V1/N1 #Volume per moles of nitrogen in m**3/Kg mole\n",
+      "Dr = (R*Tc)/(Pc*v)\n",
+      "Tr2 = T2/Tc #final reduced temp in K\n",
+      "#From figure A.2.1, reduced pressure coressponding to this Dr and Tr2 is\n",
+      "Pr2 = 4.1#final reduced pressure in atm\n",
+      "p2_ = Pr2*Pc\n",
+      "print \" iii)Final pressure achieved umath.sing Dr chart is \",\n",
+      "print \"%.6f\"%p2_,\n",
+      "print \"atm\"\n",
+      "#end\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "i)Final pressure of N2  using ideal gas law is  153.564460 atm\n"
+       ]
+      },
+      {
+       "metadata": {},
+       "output_type": "display_data",
+       "png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAEZCAYAAABfKbiYAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd4VFX+x/E3BhDXhiyPoIDGgjQVEBewz2JZBCXYZVcF\n7AU7SrOMbRV7wS4qrgUbP2VVEFAGWUEEpBOqRpo0AeklZH5/fG/MMEwgyZQzc+fzep55csuZyeFy\nk2/OOfd7DoiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiISCoVABuBdcBS4C1gb+/ck8AcYC2QD1ye\n4vqs87537RR8XxERSYBfgDbe9sHANOBRbz8IHOVttwRWASeksD4iaWEP1xUQyVBLgKHA0d5+EGup\nAPwIjKb0oJIPtI/YrwysAJoB1YB3gZXAau+zDkxgvUWSSkFFpHwqeV/rAWcDP8UosxfwN2B6KZ/x\nPtApYv8fwHJgMtAZ2A+oC9QArgM2laE+IiKSYQqwsYvV3nY/YM8Y5QYAX+3ic47Axj+qefvvAfd4\n212B74Fjylmf1cCgMrxHRETSRFnGMJ4AxgP77KbcOOBi4C9YgDncO14ZuA+YASwG+nrHKlofERFJ\nU7v7Jf4AMBU4oAyfdRvwf1g32NhSyhyKBZcrK1gfkZQr7S8gESmfXliAOAXritqdgcC/sXGT9yOO\nB4DfgZlY19Y2YHsiKyoiIulhVy2DImxAfV3Eq+duPm8EsJUdn+66FJgFrMdyYZ6l9Adq1FIRidIW\n+wGaC/Qopczz3vkpQPOocznAJOC/yaqgiIhkhhxgHpALVMEep2wUVaYdJU/RtAJ+iDp/B/bkzOCk\n1VJERMrMZZ5KSyyoFGD9xgOBvKgyHbDHM8GelqkO1PL262JB5w30rL6ISFpwGVTqAAsj9hd5x8pa\n5hngLqwvW0RE0oDLoBIuY7noVkgl4BwsA3lSjPMiIuKIy0eKF2NTXRSrh7VEdlWmrnfsAqxrrB2W\nlbwf8A5wReSbjzjiiPD8+fMTW2sREf+bDxzpuhLlVRmreC5Qld0P1Ldm54F6gNMo/emvsJj777/f\ndRXShq5FCV2LEroWJSh7T1LMX+yuFALdgK+xJ8H6Y7O3XuedfxULKO2wAf0N2LxIsVT4AoiISOK4\nzqgf4r0ivRq13203nzHKe4mIiGOa+j5LBAIB11VIG7oWJXQtSuhaJIbfn5zyugdFRKSsKlWqBBWM\nD2qpiIhIwiioiIhIwiioiIhIwiioiIhIwiioiIhIwiioiIjIn2bNiu/9CioiIsLKlXDzzXDyyfF9\njoKKiEgW27IFnnoKGnkzL8bbUnE9TYuIiDgQDsOgQXD33dC4MYweDQ0bxv+5CioiIllm/Hi44w5Y\nuxZefRXOOCNxn63uLxGRLLFgAVx2GeTlQZcu8NNPiQ0ooKAiIuJ769bBPfdA8+Zw+OEwZw5cdRXk\n5CT+eymoiIj41Pbt8MYb0KCBtVImT4YHH4R99kne99SYioiIDw0fDnfeCdWrw+DBcPzxqfm+rlsq\nbYFZwFygRyllnvfOTwGae8eqAeOwJYhnAo8mt5oiIpkhPx/at4cbboBgEEaNSl1AAbdBJQfohwWW\nxkAnYq9RfyRQH7gWeNk7vhn4O9AMONbbjjNlR0Qkc61YATfdBKeeCqefDjNmwPnnQ6UUr5rlMqi0\nxNaeLwC2AQOBvKgyHYAB3vY4oDpQy9vf6H2tigWoVUmsq4hIWtqyBZ54wpIXK1e25MU77oA993RT\nH5dBpQ6wMGJ/kXdsd2Xqets5WPfXMmAk1g0mIpIVwmH4+GMLJv/7H3z/PTz3HPz1r27r5XKgvqzr\n/EY33orftx3r/tof+BoIAKHoNweDwT+3A4GA1qEWkYw3bpy1RjZutKe72rSJ7/NCoRChUCghdXO5\nRn1rIIiNqQD0AoqAvhFlXsECxUBvfxZwGtY6iXQvsAl4Muq41qgXEd/49Vfo1csG3x9+GK64Ijm5\nJpm6Rv0EbAA+FxsXuQQYHFVmMHCFt90aWIMFlJrY+ArAXsCZwKTkVldExI21a6F3bzjuOKhfH2bP\nhq5dkxNQ4uWy+6sQ6IZ1XeUA/YF84Drv/KvAV9gTYPOADUBX79xB2AD+Ht7rP8A3qaq4iEgqFBbC\nm2/C/ffDP/4BU6ZA3bq7f59LLru/UkHdXyKSkYYNs+TFv/7VpqZv0SJ13zue7i9l1IuIpJGZM6F7\nd5g71x4VzstLfa5JPFxn1IuICLB8uWXBBwJw1lmWvNixY2YFFFBQERFxavNm6NvXFsrac09LXrzt\nNqha1XXNKkbdXyIiDoTD8NFH0LMnNGsGY8bAUUe5rlX8FFRERFLshx/g9tttipW33rIuL79Q95eI\nSIoUFECnTnDhhXD99TBhgr8CCiioiIgk3dq11s3VogU0bGjJi507wx4+/A3sw3+SiEh6KCyEV16x\nsZJly2DqVEtk3Htv1zVLHo2piIgkwdChlrx44IHw1Vc2xUo2UFAREUmg6dMtefHnn+HJJ+HcczMv\n1yQe6v4SEUmAZcts8L1NG2jXzoJLhw7ZFVBAQUVEJC6bN8Njj0GTJvCXv1jy4i23ZG7yYrzU/SUi\nUgHhMAwcaOubHHccjB1r09JnOwUVEZFyGjPGVl4sLIR33oFTT3Vdo/ShoCIiUka//GL5JmPGwL//\nDf/6lz9zTeKhyyEisht//AE9esDxx8PRR1vy4uWXK6DE4vqStMXWnZ8L9CilzPPe+SlAc+9YPWAk\nMAOYDtyS3GqKSDYqLISXXoIGDWDlSpg2De691wbkJTaX3V85QD/gDGAxMB5bkz4/okw74EhsLftW\nwMvYWvXbgNuBycA+wERgeNR7RUQqJByGIUMs3+SggyyRsVkz17XKDC6DSkts7fkCb38gkMeOgaED\nthY9wDigOlALWOq9ANZ77zkYBRURidO0aZYJv2CBrbx4zjnZl2sSD5fdX3WAhRH7i7xjuytTN6pM\nLtYtNi7B9RORLLJ0KVx7LZxxhiUtTpuWfdnwieCypRIuY7no/9LI9+0DfALcirVYdhIMBv/cDgQC\nBPw2z7SIxGXTJnjmGXj6aejSxZIXDzjAda1SKxQKEQqFEvJZLmNwayCIDdYD9AKKgL4RZV4BQljX\nGNig/mnAMqAK8AUwBHi2lO8RDofLGrtEJJsUFcEHH1jyYsuWtqTvEUe4rlV6qGTNswrFB5ctlQnY\nAHwusAS4BOgUVWYw0A0LKq2BNVhAqQT0B2ZSekAREYnp++8tebGoCN57D045xXWN/MNlUCnEAsbX\n2JNg/bGB9uu8868CX2FPgM0DNgBdvXMnAZcBU4FJ3rFewNBUVFxEMtPPP1u+ybhxlrz4z38q1yTR\n/D4Epe4vEWHNGnjkEXjzTWuh3H67ck12JZ7uL8VoEfGtbdvgxRcteXHNGpuOvk8fBZRk0txfIuI7\n4bCttti9O9SpA8OGQdOmrmuVHRRURMRXpkyx5MXFi23lxXbtlGuSSur+EhFf+O03uPpqOOssOP98\nmDoV2rdXQEk1BRURyWgbN8LDD9vswTVq2AzCN94IVaq4rll2UveXiGSkoiJ4/33o3Rtat4bx4+Hw\nw13XShRURCTjjB5tjwZXqmRZ8Sed5LpGUkxBRUQyxrx5lrw4YQI8+ihceqmSF9ON/jtEJO2tXm1P\ndLVqBS1a2KSPyoZPT/ovEZG0tW0bvPACNGwI69fDjBk2hrLXXq5rJqVR95eIpJ1wGL74wpIXDz0U\nRoyAY45xXSspCwUVEUkrkydbV9dvv8Gzz0Lbtso1ySTq/hKRtLBkCVx1lQWRCy+05MWzz1ZAyTQK\nKiLi1IYN8OCD1r1Vs6YlL95wA1RWP0pG0n+biDhRVATvvmuzBp90kj0mfNhhrmsl8VJQEZGUGzXK\nkherVIEPP4QTT3RdI0kU191fbbF15+cCPUop87x3fgrQPOL4m9jSwtOSWUERSZy5c22yx86d4a67\nYOxYBRS/cRlUcoB+WGBpjK1P3yiqTDvgSGwt+2uBlyPOveW9V0TS3KpVttriCSdAy5aQn2/Z8BqE\n9x+XQaUltvZ8AbANGAjkRZXpAAzwtscB1YHa3v5oYHXSaykiFbZ1Kzz3nCUvbtpkyYs9eyp50c9c\njqnUARZG7C8CWpWhTB1gaXKrJiLxCIdh8GDr4jr8cPj2W5uaXvzPZVAJl7FcdAO5rO8DIBgM/rkd\nCAQIBALlebuIlNNPP1ny4ooV8Pzzlnci6S0UChEKhRLyWS57NFsDQUrGRXoBRUDfiDKvACGsawxs\nUP80bIAeIBf4L1DaBA7hcLhcMUhEKmjxYrjnHhgyBB54wBIZlWuSmSrZYFeF4oPLMZUJ2AB8LlAV\nuAQYHFVmMHCFt90aWENJQBGRNLBhAwSDcOyxULs2zJkD112ngJKtXAaVQqAb8DUwE/gQyAeu814A\nXwE/YwP6rwI3Rrz/A2AMcBQ27tI1JbUWEcCSF99+Gxo0sCz4iRNtjZP99nNdM3HJ7w/0qftLJAlG\njrRxk2rV4OmnbTlf8Y94ur/UQBWRMpszB+6+G6ZMgb594aKLlGsiO3KdUS8iGWDVKrjtNst+P/FE\nS168+GIFFNmZgoqIlGrrVlvTpEED254501oq1aq5rpmkK3V/ichOwmH47DMLIPXrQygETZq4rpVk\nAgUVEdnBxIk2g/CqVfDii3DWWa5rJJlE3V8iAsCiRTZ78DnnwGWXwaRJCihSfgoqIllu/Xq4/35o\n2hTq1rWck2uuUfKiVIyCikiW2r4d3nzTBuHnzbM5ux55RMmLEh/9LSKShb791sZN9t4bBg2CVtHz\ng4tUkIKKSBaZPdumo58+HR5/HC64QLkmkljq/hLJAr//DrfcAiefDKeeasmLF16ogCKJp6Ai4mNb\nttjcXA0b2gSQM2dC9+6w556uayZ+pe4vER8Kh22spEcPCyjffQeNGrmulWQDBRURnxk/3gbh//gD\nXn4ZzjzTdY0km6j7S8QnFi6Eyy+HvDxLYpw0SQFFUs91UGmLLRE8F+hRSpnnvfNTgOblfK+I761f\nD/feC82aQW6uPeF19dWQk+O6ZpKNXAaVHKAfFhwaA52A6F7fdsCR2LLD1wIvl+O9Ir62fTv07w9H\nHQUFBTB5Mjz0EOy7r+uaSTZzOabSElsmuMDbHwjkYUsKF+sADPC2xwHVgdrAYWV4r4hvjRhhKy/u\ntx98/jn87W+uayRiXAaVOtja8sUWAdF5vbHK1AEOLsN7RXwnP9+SF/PzLXnx/POVayLpxWX3V1kX\nj9ePjGS9lSuhWzdLXPz73y3fRNnwko5ctlQWA/Ui9uthLY5dlanrlalShvcCEAwG/9wOBAIEAoGK\n1lck5bZsgRdesPXgO3WyFkrNmq5rJX4TCoUIhUIJ+SyXf+dUBmYDpwNLgB+xAffIcZF2QDfva2vg\nWe9rWd4LEA6Hy9ogEkkf4TB8+qklLzZpYl1dDRu6rpVki0rWBK5QfHDZUinEAsbX2NNc/bGgcJ13\n/lXgKyygzAM2AF13816RjPfjj5a8uH49vPYanH666xqJlJ3fe2TVUpGMsWAB9Opl68E/9JAlMCrX\nRFyIp6XiOvlRJOutWwd9+kDz5nDkkZa8eOWVCiiSmRRURBzZvh1ef92SFxctgilT4IEHYJ99XNdM\npOI0oaSIA8OHW/LiAQfAf/8Lxx/vukYiiaGgIpJCM2da8uLs2fDEE9Cxo3JNxF/U/SWSAitWwI03\nwmmnwRlnWHA57zwFFPEfBRWRJNq82XJMGjWCKlVg1iy4/XaoWtV1zUSSQ91fIkkQDsPHH1vyYtOm\nMGaMDciL+N2ugsohwIJUVUTEL374wZIXN22CN9+0ubpEssWuur8+j9j+NNkVEcl0v/4K//ynTfR4\n7bUwYYICimSfso6pHJ7UWohksLVroXdvOO446+KaMwe6dFHyomQnDdSLVFBhIbz6KjRoAEuWwNSp\nEAzC3nu7rpmIO7saUzkWWOdt7xWxDbYWyn7JqpRIuvv6a0terFkTvvzSWikisuugosa7SJQZM6B7\nd5g/35IXO3RQrolIJHV/iZTB8uVwww028N62LUyfDnl5Cigi0RRURHZh82Z47DFo3BiqVbPkxVtv\nVfKiSGmU/CgSQzgMH34IPXvaeMnYsVC/vutaiaQ/Vy2VGsBwYA4wDKheSrm2wCxgLtAj4vhFwAxg\nO6AhUkmosWPhxBNtepW334ZBgxRQRMrKVVDpiQWVo4BvvP1oOUA/LLA0xtagb+SdmwacB3yX9JpK\n1igogEsvhYsuguuvt+TFQMB1rUQyi6ug0gEY4G0PADrGKNMSW5u+ANgGDATyvHOzsFaOSNz++MO6\nuVq0sLGT2bNtKd89NOIoUm6ufmxqAcu87WXefrQ6wMKI/UXeMZGEKCyEl1+25MXly2HaNLjvPiUv\nisQjmQP1w4HaMY73idoPe69osY6JxC0chqFDLd+kVi0YMsTWhxeR+CUzqJy5i3PLsICzFDgIWB6j\nzGKgXsR+Pay1Ui7BYPDP7UAgQECd5Flt5kxbz6SgwJIXzz1XuSYioVCIUCiUkM9y9eP0OPA70Bcb\npK/OzoP1lYHZwOnAEuBHbLA+P6LMSKA7MLGU7xMOh9XgEdiwAR56CPr3h3vusVUYq1RxXSuR9FTJ\n/tKqUHxwNabyGNaSmQO08fYBDga+9LYLgW7A18BM4ENKAsp52HhLa6/8kJTUWjLS4MHQpAksWGCT\nPt56qwKKSLL4veGvlkoW+/VXuOUWy4J/8UVbG15Edi8TWyoiSbNtG/Tta48IH3+8tU4UUERSQ9O0\niK98951N/HjIITBuHBxxhOsaiWQXBRXxhRUr4K67YMQIePZZW9JXT3WJpJ66vySjFRXBa6/ZQHyN\nGpCfDxdeqIAi4opaKpKxJk+2ri6A4cOhaVO39RERtVQkA61bZwmMZ50FV14J33+vgCKSLhRUJGOE\nw/Dxx9CoEaxZY0v7XnONJn4USSfq/pKMMH8+3HQTLFoEH3wAp5ziukYiEov+xpO0tmULPPggtGoF\nbdrApEkKKCLpTC0VSVsjRtgcXY0bw8SJcOihrmskIrujoCJp57ff4I47bFnfF16wmYRFJDOo+0vS\nxvbtFkSOPRZyc20gXgFFJLOopSJpYfx4Wxd+331h1Cjr8hKRzKOWiji1Zo2Nm5x7rk1JP3KkAopI\nJlNQESfCYXj3Xcs5KSqyFRmvuELTq4hkOnV/ScrNmmWtk9Wr4bPP7HFhEfEHly2VGsBwbPXHYdiS\nwrG0BWYBc4EeEcefwFaCnAIMAvZPWk0lITZuhD594OSTIS/PxlEUUET8xWVQ6YkFlaOAb9h5jXqA\nHKAfFlgaY2vUN/LODQOaAE2xwNQryfWVOHz5JRx9NMybV7Kkb2W1k0V8x+WPdQfgNG97ABBi58DS\nEpgHFHj7A4E8rIUyPKLcOOCCJNVT4rBwoQWQadPglVdsEkgR8S+XLZVawDJve5m3H60OsDBif5F3\nLNqVwFcJrZ3EZds2ePJJaN7c8k6mTVNAEckGyW6pDAdqxzjeJ2o/7L2ixToW67O2Au/HOhkMBv/c\nDgQCBAKBMnykxOP7722dk9q1LSu+fn3XNRKRXQmFQoRCoYR8lssHOGcBAWApcBAwEmgYVaY1EMTG\nVMDGTYqAvt5+F+Aa4HRgc4zvEQ6HyxKXJBFWroQePWDoUHjmGbjoIj0iLJKJKtkPboV+el12fw0G\nOnvbnYHPYpSZANQHcoGqwCXe+8ACzV3YGEusgCIpUlQE/fvbkr777mtL+l58sQKKSDZy+WNfA/gI\nOAQbiL8YWAMcDLwOtPfKnQ08iz0J1h941Ds+Fws0q7z9scCNUd9DLZUkmzbNplcpLLSB+ObNXddI\nROIVT0vF739LKqgkyfr1EAzCO+/YeifXXAM5Oa5rJSKJkKndX5KBwmEYNMjm51q+vKSlooAiIqBp\nWqQcfv4Zbr4ZfvnFWih6kE5EoqmlIru1ZQs88gi0bGlTrEyerIAiIrGppSK79O23Nvlj/fowYYIt\nniUiUhoFFYlp2TK4804YPRqee84mgNQjwiKyO+r+kh1s3w4vvWSTP9apY+ucdOyogCIiZaOWivxp\n4kSbXqVaNVuB8eijXddIRDKNWirCH3/YU13t29v4yahRCigiUjEKKlksHIYPPrCcky1bYMYM6NJF\nXV0iUnHq/spSc+ZYq2TFCvjkEzjhBNc1EhE/UEsly2zaBPfdByeeaN1dEycqoIhI4qilkkWGDoWb\nboLjjrMExrp1XddIRPxGQSULLF4Mt90GP/0E/frB2We7rpGI+JW6v3yssNAWy2raFBo1gunTFVBE\nJLnUUvGpsWMt56RmTVvet0ED1zUSkWygoOIzq1ZBz57wxRfw1FNw6aV6RFhEUsdV91cNYDgwBxgG\nVC+lXFtsLfu5QI+I4w8BU4DJwDdAvaTVNEOEw/D225ZzsueetqRvp04KKCKSWq5+5TwOrPS+9gAO\nAHpGlckBZgNnAIuB8UAnIB/YF1jnlbsZaApcHeP7ZMXKjzNmWFfXpk22pG+LFq5rJCKZLBNXfuwA\nDPC2BwAdY5RpCczD1q/fBgwE8rxz6yLK7YMFqKyzYQP06GFrm1x6KfzwgwKKiLjlakylFrDM217m\n7UerAyyM2F8EtIrYfwS4HNgItE5CHdPa55/DLbfAKafYkr61a7uukYhIcoPKcCDWr7o+Ufth7xVt\nd/1WfbxXT+AZoGusQsFg8M/tQCBAIMOXLCwosGAyZw689Ra0aeO6RiKS6UKhEKFQKCGf5WpMZRYQ\nAJYCBwEjgYZRZVoDQWywHqAXUAT0jSp3CPAVEGteXd+MqWzdCk8/DU8+CbffDt2724C8iEiiZeKY\nymCgs7fdGfgsRpkJQH0gF6gKXOK9D+94sTxgUlJqmSZGjYJmzWwVxh9/hD59FFBEJD25aqnUAD7C\nWhkFwMXAGuBg4HWgvVfubOBZ7Emw/sCj3vFPgAbAdmA+cAOwPMb3yeiWyvLlcNddtk78c8/Beefp\nEWERSb54Wip+/xWVkUGlqAhefx3uvReuuALuvx/23dd1rUQkW8QTVJRRn2YmTbKck5wcGDECjj3W\ndY1ERMpOE0qmibVrbSbhtm3hmmts/EQBRUQyjYKKY+EwfPSRTa+ybp1lx191Feyh/xkRyUDq/nJo\n3jxbNGvJEhg4EE4+2XWNRETio7+HHdi8GR54AFq3hjPPtMWzFFBExA/UUkmx4cPhxhvhmGNsUL5e\n1s+vLCJ+oqCSIkuWWCb8jz/akr7t2+/+PSIimUbdX0lWWGiJi8ceC0ceaQPxCigi4ldqqSTRuHGW\nc7L//vaIcKNGrmskIpJcaqkkwerVFkw6doQ77rBpVhRQRCQbKKgkUDgM//mP5ZxUqgQzZ8Jll2m+\nLhHJHur+SpD8fHuqa+1aW0CrZUvXNRIRST21VOK0cSP07g2nngrnn29PdymgiEi2UlCJwxdfQJMm\n8MsvMHUq3HyzTQQpIpKt1P1VAQsWwK232uPBr71mWfEiIqKWSrls2waPPw7HHQfNm1vrRAFFRKSE\nq6BSAxgOzAGGAdVLKdcWW89+LtAjxvk7sXXrayShjjsYPdoCybffwg8/wH33QbVqyf6uIiKZxVVQ\n6YkFlaOAb7z9aDlAPyywNAY6AZHZHvWAM4Ffk1nRFSuga1fo1MlWYBwyxDLjRURkZ66CSgdggLc9\nAOgYo0xLYB62hv02YCCQF3H+aeDuZFWweEnfJk2genXLObnoIuWciIjsiquB+lrAMm97mbcfrQ6w\nMGJ/EdDK287z9qcmo3JTplhGfDgMw4ZBs2bJ+C4iIv6TzKAyHKgd43ifqP2w94oW6xjAXkBvrOur\nWKnth2Aw+Od2IBAgEAiUVpR166yL69134eGH4eqrtQKjiPhfKBQiFAol5LNcdebMAgLAUuAgYCTQ\nMKpMayCIjakA9MIG5b/ExmE2esfrAoux7rLlUZ8RDodLi02RheDTT21q+tNPtye8DjywnP8iERGf\nqGT9/BWKD666vwYDnYG+3tfPYpSZANQHcoElwCXYYH0+O3aX/QK0AFZVpCLz50O3bpZ78t57lhkv\nIiIV46pz5zGs+2oO0MbbBzgYa4kAFALdgK+BmcCHWECJtvumSAxbtsBDD0GrVhAI2CqMCigiIvHx\n+7NMMbu/vvnGJn9s2BCefx4OPdRBzURE0lQmdn85sXSprW8yZowFkw4dXNdIRMRfsuLZpu3bbV34\nY46BQw6xObsUUEREEs/3LZUJE+D662HvvSEUsmRGERFJDt+PqdSqFebxx+Hyy5UNLyJSFvGMqfj9\n12z499/D1Ej6dJMiIv6hoFK6MiU/iohIiXiCSlYM1IuISGooqIiISMIoqIiISMIoqIiISMIoqIiI\nSMIoqIiISMIoqIiISMIoqIiISMIoqIiISMK4Cio1sDXs5wDDgOqllGuLLT08F+gRcTwILAImea+2\nO71TRERSzlVQ6YkFlaOw9eZ7xiiTA/TDAkZjbCnhRt65MPA00Nx7DU1yfTNeKBRyXYW0oWtRQtei\nhK5FYrgKKh2AAd72AKBjjDItgXlAAbANGAjkRZz3+7xlCaUfmBK6FiV0LUroWiSGq6BSC1jmbS/z\n9qPVARZG7C/yjhW7GZgC9Kf07jMREUmhZAaV4cC0GK/oNRfD3ivarqYXfhk4DGgG/AY8FW9lRUQk\nc80CanvbB3n70Vqz41hJL3YcrC+WiwWrWOZRErT00ksvvfQq22seGeZxSgJET+CxGGUqA/OxoFEV\nmEzJQP1BEeVuB95PSi1FRCQj1ABGsPMjxQcDX0aUOxuYjUXNXhHH3wGmYmMqnxF7TEZERERERMSN\nN7GnxyLHVJ4A8rFWzCBgf+94LrCJkoTJl1JWy9SIdS0ewq7DZCwfqF7EuV5YUuks4KwU1TFVynMt\ncvHvfRHrOhS7EyjCeg2KZds9USz6WuTi33sCYl+LIDsmlJ8dcc7P98VOTsGSHyMvzpmUPNX2GCXj\nNbmUPqDvB7Guxb4R2zcDb3jbjbFfrlWw6zIPf03ZU55rkYt/74tY1wEsoA4FfqHkF2k23hMQ+1rk\nxijnJ7FuVoRtAAAC1UlEQVSuxf3AHTHKlvu+yPSbZjSwOurYcOyvDoBxQN2U1sidWNdiXcT2PsBK\nbzsP+ABLKi3AbpSWSa5fKpXnWvhZrOsANhvF3VHHsvGegNjXwu9KuxaxEsrLfV9kelDZnSuBryL2\nD8OadiHgZBcVcuARYAHQBXjUO3Yw1tQtFp1Y6lfF16IzOz5xmE33RR72/z016ng23hOlXQvIrnui\nWKyE8nLfF34OKn2ArZQ8brwEa+o2x5p577Njl4hf9QEOAd4Cnt1FuXBqquNU8bV4G3jGO5ZN98Vf\ngN5YV0exXU135Od7YlfXIpvuiWLlSSjf5X3h16DSBWgH/Cvi2FZKmnw/YTkw9VNbLafeB/7mbS9m\nx0H7ut6xbBF5LbLpvjgC6xefgo0h1AUmYo/kZ9s9Udq1OJDsuieKLack8fENSrq4su2+AHYeVGsL\nzABqRpWric18DHA41ozz25xhuex4LSJ/EG4G/uNtFw++VcX+OpmP/ybozKVs18Lv90UupQ86xxqo\nz6Z7IlLktfD7PQE7X4vSEsqz4b7YwQdYU3UrNvnkldijb7+y8+OAFwDTvWMTgfaprmySxboWn2A3\nzmTgU+yvsGK9sUG3WcA/UlrT5CvPtTgf/94XxddhC3Ydukad/5kdHynOhnuiLNfCz/cExP752FVC\nuZ/vCxERERERERERERERERERERERERERERGRTLMdy3mYBnwE7IVlKY/EEnWnA7c4q52IiGSUyBmS\n38Uylmtj8yyBzZo8m5LlskUyll/n/hJJV/8DjgSWYtn9AOuxheUOdlUpkURRUBFJncrYinrRU63n\nYjPijkt1hUREJPMUUjIX3XNYcCm2DzAB6OigXiIikoHWlXK8CvA1cFsK6yIiIhkuVlCphM0M+0yM\ncyIiIqVaG+PYyUARNlhf3DXWNpWVEhERERERERERERERERERERERERERERERERERERHxvf8HPRlu\n9dgw3UoAAAAASUVORK5CYII=\n",
+       "text": [
+        "<matplotlib.figure.Figure at 0x1108b2610>"
+       ]
+      },
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        " ii)Final pressure of N2 from Z chart is  138.489629  atm\n",
+        " iii)Final pressure achieved umath.sing Dr chart is  137.350000 atm\n"
+       ]
+      }
+     ],
+     "prompt_number": 4
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 2.7  Page No : 36"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "\n",
+      "#Given\n",
+      "n = 1 #n is the Kg mole of methane gas\n",
+      "T = 298 #T is the consmath.tant temperature in K\n",
+      "P1 = 1 #P1 is the initial pressure of the system\n",
+      "P2 = 100 #P2 is the final pressure of the system\n",
+      "R = 8314.4 #R is the gas consmath.tant in Nm/Kgmole deg K\n",
+      "\n",
+      "#To compute the work required\n",
+      "#(i)Umath.sing ideal gas law\n",
+      "W = R*T*math.log(P1/float(P2))\n",
+      "print \"i)Work done by the system if the gas obeys ideal gas law is \",\n",
+      "print \"%.2e\"%W,\n",
+      "print \n",
+      "\n",
+      "#(ii)Umath.sing Vander Waals' equation\n",
+      "#Given\n",
+      "#For methane\n",
+      "a = 2.32*(10**5) #Vander Wals' consmath.tant a in N/m**2\n",
+      "b = 0.0428 #Vanderwaals' consmath.tant b in m**3\n",
+      "#V1 and V2 are evaluated by trial and error umath.sing Vanderwaals' equation as P1 and P2 are known\n",
+      "V1 = 11.1 #initial volume of the gas in m**3\n",
+      "V2 = 0.089 #final volume of the gas in m**3\n",
+      "W = (R*T*math.log((V2-b)/(V1-b)))+(a*((1/V2)-(1/V1)))\n",
+      "print \" ii)Work done by the system if the gas obeys Vander Waals equation is \",\n",
+      "print \"%.2e\"%W,\n",
+      "print \"Nm\"\n",
+      "#end\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "i)Work done by the system if the gas obeys ideal gas law is  -1.14e+07\n",
+        " ii)Work done by the system if the gas obeys Vander Waals equation is  -1.10e+07 Nm\n"
+       ]
+      }
+     ],
+     "prompt_number": 5
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 2.8  Page No : 38"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "%matplotlib inline\n",
+      "import math\n",
+      "import matplotlib.pyplot as plt\n",
+      "import numpy\n",
+      "\n",
+      "#Given\n",
+      "V = 27*(10**-3) #Volume of the container in m**3\n",
+      "n = (15/70.91) #n is the Kg moles of chlorine\n",
+      "T = float(293)#T is the temperature in K\n",
+      "R = 0.08206\n",
+      "P = pow(10,(4.39-(1045/293))) #P is the vapour pressure of chlorine\n",
+      "Pc = 76.1 #Critical pressure of Chlorine\n",
+      "Tc = float(417) #Critical temperature of Chlorine \n",
+      "Pr = P/Pc #Reduced pressure of Chlorine\n",
+      "Tr = T/Tc #Critical temperature of Chlorine\n",
+      "M = 70.91 #Molecular weight of the Chlorine\n",
+      "\n",
+      "#To determine the vapour pressure of chlorine, amount of liquid Cl2 and temperature required\n",
+      "#(i)Specific volume of liquid Chlorine\n",
+      "#From figure A.2.2\n",
+      "Zg = 0.93\n",
+      "#From figure A.2.6\n",
+      "Zl = 0.013\n",
+      "vl = ((Zl*R*T)/P)\n",
+      "print \"i)Specific volume of liquid Chlorine from compressibility chart is \",\n",
+      "print \"%.6f\"%vl,\n",
+      "print \"cubic meter /Kgmole\"\n",
+      "\n",
+      "#From Francis relation, taking the consmath.tants from Table 2.3\n",
+      "D = (1.606-(216*(10**-5)*20)-(28/(200-20)))*10**3 #Density of liq Cl2 in Kg/m**3\n",
+      "Vl = M/D;\n",
+      "print \"   Specific volume of liquid Chlorine from Francis relation is \",\n",
+      "print \"%.6f\"%Vl,\n",
+      "print \"cubic meter /Kgmole\"\n",
+      "\n",
+      "#(ii)Amount of liquid Cl2 present in the cylinder\n",
+      "vg = ((Zg*R*T)/P)\n",
+      "V1 = V-vg #V1 is the volume of liquid Chlorine\n",
+      "Vct = 0.027 #volume of the container\n",
+      "Vg = (0.212-(Vct/vl))/((1/vg)-(1/vl)) #By material balance\n",
+      "W = ((V-Vg)*70.9)/vl \n",
+      "print \" ii)Weight of Chlorine at 20deg cel is \",\n",
+      "print \"%.6f\"%W\n",
+      "print \"Kg\"\n",
+      "\n",
+      "#(iii)Calculation of temperature required to evaporate all the liquid chlorine\n",
+      "#math.log P' = 4.39 - 1045/T (given)\n",
+      "#Assume the various temperature\n",
+      "Ng = 0.212#total Kg moles of gas\n",
+      "Ta = [413,415,417]\n",
+      "N = [0,0,0]\n",
+      "for i in range(0,3):\n",
+      "    Tr = Ta[i]/Tc #reduced temperature in K\n",
+      "    P = pow(10,(4.39-(1045/Ta[i])))\n",
+      "    Pr = P/Pc #reduced pressure in K\n",
+      "#From the compressibility factor chart,Z values coressponding to the above Tr &Pr are given as\n",
+      "    Z = [0.4,0.328,0.208]\n",
+      "    N[i] = (P*Vct)/(Z[i]*R*Ta[i])\n",
+      "#end\n",
+      "\n",
+      "\n",
+      "plt.plot(N,Ta)\n",
+      "plt.ylabel(\"Ta\")\n",
+      "plt.xlabel(\"N\")\n",
+      "plt.title(\"Ta vs N\")\n",
+      "plt.show()\n",
+      "T1 = numpy.interp(0.212,N,Ta) #in K\n",
+      "print \" iii)The temperature required to evaporate all the liquid chlorine is \",\n",
+      "print \"%.6f\"%(T1-273),\n",
+      "print \"deg celsius\"\n",
+      "#end\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "i)Specific volume of liquid Chlorine from compressibility chart is  0.012733 cubic meter /Kgmole\n",
+        "   Specific volume of liquid Chlorine from Francis relation is  0.045374 cubic meter /Kgmole\n",
+        " ii)Weight of Chlorine at 20deg cel is  13.112602\n",
+        "Kg\n"
+       ]
+      },
+      {
+       "metadata": {},
+       "output_type": "display_data",
+       "png": "iVBORw0KGgoAAAANSUhEUgAAAZIAAAEZCAYAAAC99aPhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4VOX5//E3RhCsIFIUaaEiFUQQRKhAqZqp4IaIphSl\n+AO+1bq0fmlEQCiLRBEEVLa6IbQqteK3FuUSEQUlI+LGFiAkRhBcQDZFWRQQkszvj/vEDGEmmSQz\nc+bMfF7XNZfnzDxz5j4mzJ37eZ7zHBARERERERERERERERERERERERERERERERERke+AA86jGDgY\ntP+HOMfiBw4BTYKe6w58Guc4RESkij4FLnPx8/3A18DMoOeUSMR1J7gdgIhHdQLeB74FtgN/B2qG\nabsIuLPMc+uA653tqcAuYB+wHmgT5jgBYAZWCTWvauAi0aZEIlI1hUAm8FPg10A34C9h2j7Psd1g\nrYFfAAuBK4FLgBbAqUAfYE85n/slMAu4rxqxi0SVEolI1awBVmDjJp8DTwHpYdrOB9oDTZ39m4B5\nwFHnURc4D/v3+DGws5zPDQAPAtdiCUnEdUokIlXTEngV2IF1SY3HqpNQDmDVR0lV0hf4t7O9FHgU\neAzr3pqJJZbyfO28534ssYi4SolEpGqeAPKBc7AuqVGU/+9pLpZIfg3UBrKDXvs78CuswmgJDIvg\n8x8Cfgt0rGzgItGmRCJSNadglcZBoBXw5wravwachY1tvBD0/K+AzthA/UHgMFBUznFqOP/dBzwC\nDK9s4CLRpkQiUjVDgX7Afmx85AXK72Y6AryEDco/H/R8Pef93wCfYd1WD5VznODPmI4N+qt7S5Je\nGpADLHD2+wB52F9dwWV5P6ddyaMIaBfieA2AJcBGYDFQPyZRi4hIwrgbG1h8xdlvhfUDZwMdwrzn\nfGBTmNcmA/c428OBidEJU0REElET4E1sUHBBmdfKSyQTgHFhXisAGjnbZzr7IiLikliPkUzFZqAU\nV/J9N2CzXEJphE2TxPlvozDtREQkDmKZSHoCu7HxjhoVtA3WGZu9kh9B2wAaaBQRcdWJMTx2V6AX\n0AObN18PmAMMqOB9fTl2VktZu7AurZ1AYyxZHeeXv/xlYPPmzZUMWUQk5W3Gro+KWCwrkpHYkhBn\nY8lhKccnkbKVygnYrK4XCO8VYKCzPRBbfuI4mzdvJhAIJO1j7Nixrsegc9P56fwqfixfHqB9+wCX\nXhpg3Tr346/oAfyyUt/0xPc6kpIuqAxgK9AFWzZiUVCbS4EvsPn0wWZROlV4InA5Nv33MjRrS0QS\n0I4dMGAA3Hgj3HMP+P3QLtQFDUkgll1bwd52HgAvO49Q/FiXWFm3Bm1/g92DQUQk4Rw9CjNmwIMP\nwq23QkEBnHKK21HFVrwSiUSZz+dzO4SYSeZzA52f15V3fkuWwF//Cs2awXvvQcuWcQvLVZWZTeU1\nAae/T0Qkpj7/HIYMgTVrYNo0uPZaqOHRb9caFnilotdaWyIiVXT4MIwbBx06wAUXQF4e9Orl3SRS\nVeraEhGppEAAFiyAu+6yJLJmDZx1lttRuUeJRESkEjZuhMxM+OwzmDkTLr/c7Yjcp64tEZEIfPcd\njBgBXbtC9+6wbp2SSAklEhGRcgQC8MILcN55sH075ObawHqtWm5HljjUtSUiEsb69Tadd98+Sya/\n+Y3bESUmVSQiImXs3WsJpHt3uzJ91SolkfIokYiIOIqL4Z//hFat4MgRyM+HP/8Z0tLcjiyxqWtL\nRARYsQIGDbKksXAhdOxY8XvEqCIRkZT21Vfwpz/B9dfDnXfC8uVKIpWlRCIiKamwEP7+d2jTBurV\ng48+stV6T9C3YqWpa0tEUs6yZfC//wunnw7Z2ZZMpOqUSEQkZXz5JQwbZt1XU6ZA796pty5WLKiI\nE5Gk98MPMGmSLazYvLl1Y/3+90oi0aKKRESS2uuv29pYLVvCBx/AOZW6G7lEQolERJLSli1w9922\ntPu0aXDNNW5HlLzUtSUiSeXgQRg7Fjp1gs6dYcMGJZFYi0ciSQNygAXOfh8gDygCOpRp2w54H9gA\nrAdOCnG8LGCbc8wc4KqoRywinhMIwEsvQevW8PHHkJMDf/sbnBTqW0SiKh5dW5lAPlDX2c8FMoCZ\nIWL5F/D/nDanAUdDHC8ATHEeIiJ89JGtjbVjBzz9NPz2t25HlFpiXZE0AXoAsym9B3ABsDFE2yuw\nKiTX2f8WKA5zXM21EBH277fpvJdeCj17WhWiJBJ/sU4kU4FhhE8IwVpg1cbrwGrnfeEMAtYB/wDq\nVzNGEfGYQACee87uEbJnj42DZGZCzZpuR5aaYplIegK7sXGMSCqImsDFQD/nvxnAZSHaPQGcDbQH\ndgCPRCNYEfGGtWvhkktsJta8ebZab6NGbkeV2mI5RtIV6IV1bdUG6gFzgAFh2m8FlgHfOPuvYYPx\nS8u02x20PZvSQfzjZGVl/bjt8/nw+XyRxi4iCeabb2D0aBtQHzcObr5Zy7tHg9/vx+/3V+sY8Rpr\nSAeGAtcGPZftPLfa2a8PvIVVI0eBRdiA+qIyx2qMVSIAg4GLsCqmrEAgEIhG7CLioqIimD0b7r0X\n+vSB+++HBg3cjip51bDL/SuVG+J5QWLJt3oGMANoCCzEur6uBvZiiWOl03YhpUlkFtaltQaYhHVr\nBYBPgdvjE76IxNv779viiiefDIsX2xInkniSefaTKhIRj9q5E0aMgCVLYPJk6NdP62LFS1UqEl3Z\nLiIJ4+hRmDoV2raFM86AggK46SYlkUSntbZEJCEsXWq3um3SBN55x+6bLt6gRCIirvriCxg6FFau\ntGrkuutUgXiNurZExBWHD8P48dChg62PlZ9v901XEvEeVSQiEnevvgp33QXt2sGqVdCsmdsRSXUo\nkYhI3HzyiSWQTZvgscfgyivdjkiiQV1bIhJz338Po0ZBly6Qng65uUoiyUSJRERiJhCA//zHFlf8\n/HNYv95W661Vy+3IJJrUtSUiMbFhg90jZM8e+Pe/baFFSU6qSEQkqvbtg8GD4bLLoHdvWL1aSSTZ\nKZGISFQUF9vdCVu1sjGRvDy48044Uf0eSU8/YhGptlWrbHFFgFdegYsucjceiS9VJCJSZV9/Dbfd\nBtdeC7ffDu+9pySSipRIRKTSCgvtOpDWrW2J948+gj/+EU7QN0pKUteWiFTK8uXWjXXaabbQ4vnn\nux2RuE2JREQisn073HMPvP02PPww3HCD1sUSo0JURMp15Ag89JCti/WLX1g31o03KolIKVUkIhLW\n4sV2UWHz5nbb2xYt3I5IEpESiYgc57PP4O67Yd06mDYNevZUBSLhxaNrKw3IARY4+32APKAI6FCm\nbTvgfWADsB44KcTxGgBLgI3AYqB+9EMWSU2HDsF990HHjnafkLw8m9qrJCLliUciyQTygYCznwtk\nAMvKtDsR+BdwG3A+kA4cDXG8EVgiaQm85eyLSDUEAjB/PrRpY2tk5eTA6NFQu7bbkYkXxLprqwnQ\nAxgP3O08VxCm7RVYFZLr7H8bpl0vLMkAPAv4UTIRqbKPP4bMTLvl7axZ0K2b2xGJ18S6IpkKDAOK\nI2jbAqtaXgdWO+8LpRGwy9ne5eyLSCUdOADDh8PFF9u9QdatUxKRqollRdIT2I2Nj/giaF8TuBj4\nFXAI67ZaDSwt5z0BSrvMjpOVlfXjts/nw+eLJAyR5BYIwNy5dk1I9+52k6kzz3Q7KnGL3+/H7/dX\n6xixHEKbAPQHCoHaQD1gHjDAeT0bGAKscfZvBK4G/sfZHw0cBh4uc9wCLDHtBBo7x2kV4vMDgUDY\nHCOSktavt6vSv/sOHn0UunZ1OyJJNDVsZkWlckMsu7ZGAk2Bs4G+WGUxoEyb4GDfANoCdbBKKR2b\n3VXWK8BAZ3sgMD96IYskp2+/hUGD4PLL4aabYOVKJRGJnnhe2V5SHmQAW4EuwEJgkfP8XmAKsBLr\nDlsd9NosoKOzPRG4HJv+e5mzLyIhFBfD7Nl2q9uiIsjPt1V609LcjkySSTLPDlfXlqS0Dz+0bqxa\ntawb68IL3Y5IvCDRurZExAW7d8PNN8PvfmfLmyxfriQisaVEIpIkCgth+nS7qLBBA1tcsX9/XZUu\nsae1tkSSgN9vg+lnngnLltmYiEi8KJGIeNi2bTB0qK3MO2WKdWepApF4U9eWiAf98AM8+CC0bw8t\nW1o3Vu/eSiLiDlUkIh7z2mu2Nlbr1rBihd0rRMRNSiQiHrF5MwweDAUFMGMGXH212xGJGHVtiSS4\ngwdhzBjo3Bl+8xtbG0tJRBKJKhKRBBUIwLx5MGSILWeydi00aeJ2VCLHUyIRSUD5+XYx4a5d8Oyz\noIWrJZGpa0skgezfbxVIejpcd53dqVBJRBKdEolIAiguhjlzoFUr2LfP7pU+aBCcqD4D8QD9moq4\nbM0aW1yxsNDum96pk9sRiVSOKhIRl+zZA3fcAT16wC23wAcfKImINymRiMRZURE8+aSth1Wrll2V\nfsstcIL+NYpHqWtLJI7ee8+6serWhTffhHbt3I5IpPqUSETiYMcOGD4cli6Fhx6Cvn21LpYkDxXT\nIjF09Kitytu2LfzsZ7a8yR/+oCQiyUUViUiMvPmmXVR41lnw7rtw7rluRyQSG/GoSNKAHGCBs98H\nyAOKgA5B7ZoBh5y2OcDjYY6XBWwLandVtAMWqY7PP4ff/x5uuw0mTrTVepVEJJnFI5FkAvlAwNnP\nBTKAZSHafgJc6Dz+EuZ4AWBKULvXoxmsSFUdPgzjxkGHDjaInpcHvXqpG0uSX6wTSROgBzAbKPnn\nVABsrOZx9U9TEkYgAK+8YvdKX7sWVq+Ge++FOnXcjkwkPmKdSKYCw4DiCNufjXVX+YGLy2k3CFgH\n/AOoX434RKpl0ya45hqbkfXkk7Zab7NmbkclEl+xHGzvCezGEoMvgvbbgabAt9jYyXygDXCgTLsn\ngPud7XHAI8AtoQ6YlZX147bP58On1e8kSr77DsaPh1mzYMQIG1SvVcvtqEQqz+/34/f7q3WMWHYR\nTQD6A4VAbaAeMA8Y4LyeDQwB1oR5f0Wvgw3QLwDahngtEAgEQjwtUnWBAPzf/8GwYbYq7+TJ0Lix\n21GJRE8NG9SrVG6IZUUy0nkApANDKU0iJYKDbYhVI0VAc6AFsCXEcRsDO5ztDGzwXiTmcnNtRd69\ne2HuXLi4vM5XkRQSzwsSS8qDDGAr0AVYCCxynk/Hxj1ygBeB24G9zmuzKJ0qPAlY77RNBwbHOnBJ\nbXv3QmYmdOsGN9xgg+lKIiKlknn2k7q2pFqKi+GZZ2DkSLvJ1Pjx0LCh21GJxFaidW2JeNbKlba4\n4gknwMKF0LGj2xGJJC6ttSUS5Kuv4NZb7ULCv/zFljZREhEpnxKJCHZ3wkcftYsK69a1xRUHDtQ9\nQkQioa4tSXnLllk31umnQ3a2JRMRiZwSiaSsL7+060GWL4dHHrGFFrUulkjlqXCXlHPkCEyaBBdc\nAM2b261u+/RREhGpKlUkklLeeMOWM2nZEj74AM45x+2IRLxPiURSwqefwuDBsGEDTJ9uCy2KSHSo\na0uS2sGDMHYsXHQRdOpkiURJRCS6VJFIUgoE4OWX4e67oXNnyMmBpk3djkokOSmRSNIpKLBxkO3b\n4emn4be/dTsikeSmri1JGvv323TeSy6BHj2sClESEYk9JRLxvEAAnnsOzjsPvv7axkHuugtq1nQ7\nMpHUoK4t8bS1a+0eIYcO2W1uu3RxOyKR1KOKRDzpm2/gzjvhyiuhf3/48EMlERG3KJGIpxQVwVNP\nWTdWjRp2Vfptt0FamtuRiaQudW2JZ7z/vnVj1a5tV6i3b+92RCICSiTiAbt2wYgRsHgxTJ4M/fpp\nXSyRRKKuLUlYR4/CtGlw/vm2xHtBAdx0k5KISKKJRyJJA3KABc5+HyAPKAI6BLVrBhxy2uYAj4c5\nXgNgCbARWAzUj3rE4rqlS+HCC+G11+Cdd6wSqVvX7ahEJJRIurbqALcAbYDaznMB4OYIPyMTyAdK\nvgZygQxgZoi2nwAXVnC8EVgimQwMd/ZHRBiLJLitW2HIEFixAqZOheuvVwUikugiqUj+BTQCrgT8\nQBPguwiP3wToAcwGSr4OCrBqoqp6Ac86288C11fjWJIgDh+G8eOtCmndGvLzISNDSUTEC8pLJCXV\nyjnAGCx5PIslhs4RHn8qMAwojrD92Vi3lh+4OEybRsAuZ3uXsy8etnixjYOsWgUrV0JWFpx8sttR\niUikyuvaWoGNYRxx9vcBbYGdwOkRHLsnsBtLDL4I2m8HmgLfOp87H+tOO1DOewLOI6SsrKwft30+\nHz5fJGFIPG3bBn37wvPPw1VXuR2NSOrx+/34/f5qHaO8joMcbLziT8BLWBJ5GjgFuBd4soJjTwD6\nA4XY2Eo9YB4wwHk9GxgCrAnz/nCvF2CJaSfQ2GnXKsT7A4FA2BwjCeKOO+DUU+3WtyLivhrWn1yp\nTuXyGm8DppTT5pFKfE46MBS4Nui5bOe51c5+Q6waKQKaA8uA84G9ZY41GdgDTMIG2esTerBdiSTB\nffKJLWvy8cfw05+6HY2IQNUSSXljJGnYTKtTwjwqq+RbPQPYCnQBFgKLnOfTgXVYJfQicDulSWQW\n0NHZnghcjg3YX+bsiwdlZUFmppKIiNdF0rXlVapIEtiGDdCtm1Uluj5EJHFEuyIRiZkxY2D4cCUR\nkWRQXtb5KTYW4VWqSBLUihXQuzds3Ah16rgdjYgEi/Zgu9cpkSSoK66wRHL77W5HIiJlqWtLEl52\nNmzZAjdHusCOiCQ8JRKJm0AARo2C++7T/dRFkokSicTNwoWwf79dyS4iyUOJROKiuBhGj4YHHtBt\ncUWSjRKJxMWLL0KtWnDddW5HIiLRpllbEnOFhdCmDTz2GHTv7nY0IlIezdqShDRnDvzsZ3Ylu4gk\nH1UkElM//AAtW8LcudC1q9vRiEhFVJFIwpk5E9q1UxIRSWaqSCRmvv8ezjkHFi2C9u3djkZEIqGK\nRBLKjBmQnq4kIpLsVJFITOzdCy1awLvv2hiJiHiDKhJJGA8/DL16KYmIpAJVJBJ1u3ZB69awZg2c\ndZbb0YhIZWgZ+WMpkbjkrrtsgcbp092OREQqS4nkWEokLvjiC7jwQsjPh0aN3I5GRCorUcdI0rD7\nvy9w9vsAeUAR0CFE+18A3wFDwhwvC9jmHDMHuCqKsUo1jRtnN6xSEhFJHSfG4TMygXyg5O7cuUAG\nMDNM+ynAwnKOF3DaTIlWgBIdGzfC/Pn2XxFJHbGuSJoAPYDZlJZKBUC4r5rrgS1Y4ilPMnfJedbY\nsTB4MJx2mtuRiEg8xTqRTAWGAcURtD0FuAfruqrIIGAd8A+gflWDk+hZt85uo/vXv7odiYjEWyy7\ntnoCu7FxDF8E7bOwxHOQ8iuOJ4D7ne1xwCPALSEPmJX147bP58PniyQMqYoxY+Bvf4NTTnE7EhGp\nDL/fj9/vr9YxYtlFNAHoDxQCtYF6wDxggPN6NjagvsbZXwY0dbbrY1XMGODxcj6jGTaI3zbEa5q1\nFScffAA33GBjI7Vrux2NiFRHVWZtxbIiGek8ANKBoZQmkRLBwV4atD0WOEDoJNIY2OFsZ2CD9+Ki\nUaPg3nuVRERSVTyXSCkpDzKArUAXbHbWogjeO4vSqcKTgPXYGEk6MDi6YUplvPUWbN0KAwe6HYmI\nuCWZZz+payvGAgHo0sVmavXt63Y0IhINiXpBoiSpBQvg8GEbHxGR1BWPCxIlCRUX29jIhAlwgv4c\nEUlp+gqQKnnhBZvq27On25GIiNs0RiKVdvSoLRM/cyZcdpnb0YhINGmMROLimWfsPiNKIiICqkik\nkg4ftlvo/ve/0Lmz29GISLSpIpGYe/JJ6NhRSURESqkikYgdOGDVyJIl0DbUojQi4nmqSCSmpk+H\nbt2URETkWKpIJCLffAPnngvvvw/nnON2NCISK6pIJGYeeggyMpREROR4qkikQjt3Qps2sHYtNG1a\ncXsR8a6qVCRKJFKhQYOgZk2YMsXtSEQk1pRIjqVEEgWffw4dOsBHH8EZZ7gdjYjEmhLJsZRIouDm\nm+HnP4dx49yORETiIdHukCgeV1BgS8Vv2uR2JCKSyDRrS8IaOxaGDIH69d2OREQSmbq2JKScHLjm\nGqtGfvITt6MRkXjRdSQSNaNHw8iRSiIiUrF4JJI0IAdY4Oz3AfKAIqBDiPa/AL4DhoQ5XgNgCbAR\nWAyo4yXK3n0X8vLg1lvdjkREvCAeiSQTyAdK+plygQxgWZj2U4CF5RxvBJZIWgJvOfsSJYGAVSJj\nx8JJJ7kdjYh4QawTSROgBzCb0j63AqyaCOV6YAuWeMLpBTzrbD/rvEeiZMkS2LUL+vd3OxIR8YpY\nJ5KpwDCgOIK2pwD3AFkVtGsE7HK2dzn7EgWBAIwaZdeMnKiJ4SISoVh+XfQEdmPjI74I2mdhiecg\nkc8YCFDaZXb8AbOyftz2+Xz4fJGEkbrmz4fCQujd2+1IRCRe/H4/fr+/WseI5fTfCUB/oBCoDdQD\n5gEDnNezsQH1Nc7+MqBkScD6WBUzBni8zHELsMS0E2jsHKdViM/X9N9KKCqCdu1sld8ePdyORkTc\nkmjTf0diieFsoC+wlNIkUiI42EudtmcD04DxHJ9EAF4BBjrbA4H50Qs5dc2dC6edBldf7XYkIuI1\n8byOpKQ8yAC2Al2w2VmLInjvLKCjsz0RuBwbsL/M2ZdqOHLEZmmNHw81kvkSVRGJiWT+2lDXVoSe\nfBJefhneeMPtSETEbVr991hKJBE4dAhatLCB9l/9yu1oRMRtiTZGIh7w+OPQubOSiIhUnSqSFLZ/\nv1UjS5farXRFRFSRSKVMnQpXXqkkIiLVo4okRe3ZA+eeCytWQPPmbkcjIolCFYlEbNIk6NNHSURE\nqk8VSQravh3atoX16+1+7CIiJTT991hKJGHceSecfLIthyIiEkyJ5FhKJCFs2QIXXQQffwwNG7od\njYgkGo2RSIXuuw8GDVISEZHo0V0nUkh+Prz+Omza5HYkIpJMVJGkkHvvhaFDoV49tyMRkWSiMZIU\nsXo19Opl1cjJJ7sdjYgkKo2RSFijR9tDSUREok2JJAUsW2aztG65xe1IRCQZKZEkuUAARo2CrCyo\nVcvtaEQkGSmRJLk33rB1tW66ye1IRCRZKZEkseJiGDkSHngA0tLcjkZEklU8EkkakAMscPb7AHlA\nEdAhqF0np10OsB64MczxsoBtQW2vinrESeKll+CEEyAjw+1IRCSZxeOCxEwgH6jr7OcCGcDMMu1y\ngY5AMXAmsAH4L5ZwggWAKc5DwigqgjFjYNo0qJHMk7xFxHWxrkiaAD2A2ZTOSy4ANoZoewhLIgB1\ngH0cn0RK6KuxAs89B2ecAVdc4XYkIpLsYp1IpgLDKE0QFemEdXvlAXeX024QsA74B1C/OgEmox9+\nsFla48erGhGR2ItlIukJ7MbGMSL9OlsBtMHGTqYDp4Zo8wRwNtAe2AE8Uu1Ik8zs2XDeeXDxxW5H\nIiKpIJZjJF2BXljXVm2gHjAHGBDBewuAzcA5wOoyr+0O2p5N6SD+cbKysn7c9vl8+Hy+CD7a2w4e\ntErk1VfdjkREvMDv9+P3+6t1jHh1fKQDQ4Frg57Ldp4rSRTNsNlYhcBZwDvA+cD+MsdqjFUiAIOB\ni4B+IT4zJdfamjwZVq6EF190OxIR8aKqrLUVz2XkS77VM4AZQENgIdb1dTVwCTAcOOo8bqM0iczC\nurTWAJOwbq0A8Clwe3zCT3z79sHDD8Pbb7sdiYikkmQeik25imTsWPjiC3j6abcjERGv0q12j5VS\nieSrr2yAfdUqaNbM7WhExKuUSI6VUolkyBCb9vvoo25HIiJepkRyrJRJJNu2wQUXwIYN0Lix29GI\niJcpkRwrZRLJHXdA/fowcaLbkYiI1ymRHCslEsnmzdC5M2zcCA0auB2NiHidbrWbgrKyIDNTSURE\n3KOKxMM2bIDu3WHTJqhbt+L2IiIVUUWSYsaMgXvuURIREXepIvGoFSugd28bG6lTx+1oRCRZqCJJ\nIaNHW0WiJCIiblMi8aDsbNiyBf74R7cjERFRIvGcQABGjYL77oOaNd2ORkREicRzXnsNDhyAvn3d\njkRExCiReEhxsVUjDzwAaWluRyMiYpRIPOTFF+Gkk6BXL7cjEREppem/HlFYCG3awGOP2UWIIiKx\noOm/SWzOHPj5z6FbN7cjERE5lioSD/jhB2jZEl54AX79a7ejEZFkpookST31FLRrpyQiIokpHokk\nDcgBFjj7fYA8oAjoENSuk9MuB1gP3BjmeA2AJcBGYDFQP/ohJ47vv4cJE2ymlohIIopHIskE8oGS\nfqZcIANYVqZdLtARuBC4AngMS0JljcASSUvgLWc/ac2YAenpdgfEYH6/35V44iGZzw10fl6X7OdX\nFbFOJE2AHsBsSvvcCrBqoqxDQLGzXQfYh1UtZfUCnnW2nwWuj1awiWbvXpgyBe6///jXkvmXOZnP\nDXR+Xpfs51cVsU4kU4FhlCaIinTCur3ygLvDtGkE7HK2dzn7Senhh+G662ygXUQkUZ0Yw2P3BHZj\nYx6+CN+zAmgDtAJeB/xYZRJOgNIus6Syezc88QTk5LgdiYiIeyYAW4FPgR3A98CcoNezOXawvay3\nsDGTsgqAM53txs5+KJ9Qmmj00EMPPfSI7PEJCSqd0llbJbI5NlE0o7RCOgv4AqgX4liTgeHO9ghg\nYtSiFBGRhJUOvOJsZ2CVyiFgJ7DIeb4/sAHrClsBXBX0/lmUJp0GwJukyPRfERERERFJYFdhYySb\nKO3uCuUioBD4XTyCiqKKzs+HTUYouZBzdNwii45Ifn4+7Nw2YJMvvKSi8xtK6c8uF/sd9VKFXdH5\nNcQmzazFfn7/E7fIqq+iczsNeBlYB3yITRLyin9iM15zy2kzAzv3ddi1fUkrDRsUagbUxH5ZzwvT\nbinwKtA7XsFFQSTn56O0y9BrIjm/+thU8CbOfsN4BRcFkf5+luiJddl6RSTnlwU86Gw3BPYQ25mi\n0RLJuT2n795OAAADSUlEQVQEjHG2z8VbP7tLsOQQLpH0AF5ztjsDH1R0QC+vtdUJ+2F/BhwFXgCu\nC9FuEPBf4Ku4RRYdkZ6fVxfejOT8+gHzgG3O/tfxCi4KIv35legHzI19WFETyfntoHTCTD0skRTG\nKb7qiOTczsMmDAF8jCWd0+MTXrW9A3xbzuvBF31/iP1BV+71el5OJD/HBu1LbHOeK9vmOuAJZz8Q\nh7iiJZLzCwBdsfLzNaB1fEKLikjOrwU2uSIbWIVNyPCKSM6vxMnAlVjS9IpIzm8W1uWzHfsdzYxP\naNUWybmto7SrvBM207QJySHU+Zd7bl4oM8OJJClMw6YIB7C/3L3013sk57cGaAocBK4G5mNrkHlB\nJOdXE7vWqBv2Zfs+VmZvimFc0VKZP1quBZYDe2MUSyxEcn4jsW4hH/BLbI28C4ADsQsrKiI5t4nA\ndErHt3IIvaSTV5X9riz3/4mXE8mX2JdoiaaUdoGU6IiVpWB9tFdjpaoXxhUiOb/gf5CLgMexv+C/\niW1oURHJ+W3FurMOOY9l2BeRFxJJJOdXoi/e6taCyM6vKzDe2d6MXZx8LlZdJrJI/+3dHLT/KbAl\nxnHFS9nzb+I8l5ROxH45mwG1qHgw82m8NWsrkvNrROlfDp2wPl2viOT8WmGDmGlYRZKLd7rvIv39\nPBUbO6gTt8iiI5LzmwKMdbYbYV/GDeIUX3VEcm6nOq8B3Ao8E6fYoqUZkQ22dyGCwXavuxob6PoE\n+Jvz3O3OoyyvJRKo+PzuxKZVrgXew37oXhLJz28oNnMrF/hrXKOrvkjObyDwfJzjipaKzq8htqLF\nOuzn1y/eAVZDRef2a+f1Amwyz6nxDrAa5mLjVkewqv9mjv+9fBQ793WUv5SViIiIiIiIiIiIiIiI\niIiIiIiIiIiIeEua2wGIpIhioC52Mzaw62MuB952LSKRKPHyoo0iXnIEuzvoT519Ly0gKlIuJRKR\n+DgKPAUMdjsQkWhTIhGJn8eBmyi9R4dIUlAiEYmfA8AcvLdmmEi5lEhE4msacAvwE7cDEYkWJRKR\n+PoW+A+WTDTgLklBiUQkPoKTxiPYEusiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIqnh/wM7\nGvJ1b2hx8QAAAABJRU5ErkJggg==\n",
+       "text": [
+        "<matplotlib.figure.Figure at 0x1108b5a50>"
+       ]
+      },
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        " iii)The temperature required to evaporate all the liquid chlorine is  140.000000 deg celsius\n"
+       ]
+      }
+     ],
+     "prompt_number": 6
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 2.9  Page No : 40"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math\n",
+      "import numpy\n",
+      "\n",
+      "\n",
+      "#Given\n",
+      "N1 = 0.7 #Kg mole of CH4\n",
+      "N2 = 0.3 #Kg mole of N2\n",
+      "R = 0.08206 #Gas consmath.tant\n",
+      "T = float(323) #Temperature in Kelvin\n",
+      "V = 0.04 #Volume in m**3\n",
+      "a1 = 2.280\n",
+      "b1 = 0.0428 #Vanderwaals consmath.tants for CH4\n",
+      "a2 = 1.345\n",
+      "b2 = 0.0386 #Vanderwaals consmath.tants for N2\n",
+      "Tc1 = 191\n",
+      "Pc1 = 45.8 #Critical temperature in K and pressure of CH4 in atm\n",
+      "Tc2 = 126\n",
+      "Pc2 = 33.5 #Critical temperature in K and pressure of N2 in atm\n",
+      "\n",
+      "#To find Approx Value\n",
+      "def approx(V,n):\n",
+      "    A = numpy.around([V*pow(10,n)])/pow(10,n)\n",
+      "    return A[0]\n",
+      "\n",
+      "#To calculate the pressure exerted by the gas mixture\n",
+      "#(i)Umath.sing ideal gas law\n",
+      "P = (N1+N2)*((R*T)/V)\n",
+      "print \"i) Pressure exerted by the gas mixture umath.sing ideal gas law is \",\n",
+      "print \"%d\"%P,\n",
+      "print \"atm\"\n",
+      "\n",
+      "#(ii)Umath.sing Vander waal equation\n",
+      "P1 = ((N1*R*T)/(V-(N1*b1)))-((a1*(N1**2))/(V**2)) #Partial pressure of CH4\n",
+      "P2 = ((N2*R*T)/(V-(N2*b2)))-((a2*(N2**2))/(V**2)) #Partial pressure of N2\n",
+      "Pt = P1+P2 \n",
+      "print \"ii) Pressure exerted by the gas mixture umath.sing Vander waal equation is \",\n",
+      "print \"%.6f\"%Pt,\n",
+      "print \"atm\"\n",
+      "\n",
+      "#(iii)Umath.sing Zchart and Dalton's law\n",
+      "Tra = T/Tc1 #reduced temperature of CH4\n",
+      "Trb = T/Tc2 #reduced temperature of N2\n",
+      "#Asssume the pressure\n",
+      "P = [660,732,793,815,831]\n",
+      "Pa =[]\n",
+      "Pb = []\n",
+      "Pra = []\n",
+      "Prb = []\n",
+      "for i in range(0,5):\n",
+      "    Pa.append(N1*P[i]) # partial pressure of CH4 for the ith total pressure\n",
+      "    Pb.append(N2*P[i]) # partial pressure of N2 for the ith total pressure\n",
+      "    Pra.append(Pa[i]/Pc1) #reduced pressure of CH4 for the ith total pressure\n",
+      "    Prb.append(Pb[i]/Pc2) #reduced pressure of N2 for the ith total pressure\n",
+      "#end\n",
+      "\n",
+      "#For the above Pr and Tr values compressibility factors  from the figure A.2.3 are given as\n",
+      "Za = [1.154,1.280,1.331,1.370,1.390] #Z values of CH4  \n",
+      "Zb = [1,1,1,1,1]#Z values of N2\n",
+      "V3 = 0.0421\n",
+      "for i in range(0,5):\n",
+      "    Pa[i] = Za[i]*N1*((R*T)/V);#partial pressure of CH4 coressponding to the ith total presure\n",
+      "    Pb[i] = Zb[i]*N2*((R*T)/V);#partial pressure of N2 coressponding to the ith total pressure\n",
+      "    Pt = Pa[i]+Pb[i] #total pressure of the gas mixture\n",
+      "    if Pt-P[i] < 15:\n",
+      "        print \"iii) pressure exerted by the gas mixture umath.sing Z chart and Dalton Law is \",\n",
+      "        print \"%d\"%Pt,\n",
+      "        print \"atm\"\n",
+      "    #end\n",
+      "#end\n",
+      "\n",
+      "#(iv)Umath.sing Amagat's law and Z chart\n",
+      "P = [1000,1200,1500,1700]\n",
+      "for i in range(0,4):\n",
+      "    Pra[i] = P[i]/Pc1\n",
+      "    Prb[i] = P[i]/Pc2\n",
+      "#end\n",
+      "#For the above Pr and Tr values compressibility factors  from the figure A.2.3 are given as\n",
+      "Za = [1.87,2.14,2.52,2.77]\n",
+      "Zb = [1.80,2.10,2.37,2.54]\n",
+      "Va = []\n",
+      "Vb = []\n",
+      "V1 = []\n",
+      "for i in range(0,4):\n",
+      "    Va.append(approx((N1*Za[i]*((R*T)/P[i])),4))\n",
+      "    Vb.append(approx((N2*Zb[i]*((R*T)/P[i])),4))\n",
+      "    V1.append(approx((Va[i]+Vb[i]),4))\n",
+      "    if V1[i]-V <= 0.003:\n",
+      "        print \"iv) Pressure exerted by the gas mixture umath.sing Amagat law and Zchart is \",\n",
+      "        print \"%d\"%P[i],\n",
+      "        print \"atm\"\n",
+      "#end\n",
+      "#end\n",
+      "#end\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "i) Pressure exerted by the gas mixture umath.sing ideal gas law is  662 atm\n",
+        "ii) Pressure exerted by the gas mixture umath.sing Vander waal equation is  1353.867785 atm\n",
+        "iii) pressure exerted by the gas mixture umath.sing Z chart and Dalton Law is  843 atm\n",
+        "iv) Pressure exerted by the gas mixture umath.sing Amagat law and Zchart is  1700 atm\n"
+       ]
+      }
+     ],
+     "prompt_number": 8
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 2.10  Page No : 41"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "\n",
+      "#Given\n",
+      "yN2 = 1.0/4 #mole faction of N2 in the mixture\n",
+      "yH2 = 3.0/4 #mole fraction of H2 in the mixture\n",
+      "V = 5.7 #V is the rate at which mixture enters in m**3 in 1 hour\n",
+      "P = float(600) #P is in atm\n",
+      "T = float(298) #T is in K\n",
+      "TcN2 = float(126) #critical temp of N2 in K\n",
+      "TcH2 = 33.3 #critical temp of H2 in K\n",
+      "TcNh3 = 406.0 #critical temp of NH3 in K\n",
+      "PcN2 = 33.5 #critical pressure of N2 in atm\n",
+      "PcH2 = 12.8 #critical pressure of H2 in atm\n",
+      "PcNH3 = 111.0 #critical pressure of NH3 in atm\n",
+      "R = 0.082 #gas consmath.tant\n",
+      "\n",
+      "#To calculate the amount of ammonia leaving the reactor and the velocity of gaseous product leaving the reactor\n",
+      "#(i)Calculation of amount of NH3 leaving the reactor\n",
+      "Tcm = (TcN2*yN2)+(TcH2*yH2) #critical temperature of the mixture\n",
+      "Pcm = (PcN2*yN2)+(PcH2*yH2) #critical pressure of the mixture\n",
+      "Trm = T/Tcm \n",
+      "Prm = P/Pcm \n",
+      "#From figure A.2.3\n",
+      "Zm = 1.57 #compressibility factor of the mixture\n",
+      "N = (P*V)/(Zm*R*T) #Kg mole of the mixture \n",
+      "N1 = 0.25*N #Kg mole of N2 in feed\n",
+      "#N2+3H2 - 2NH3\n",
+      "W = 2*0.15*N1*17 \n",
+      "print \"i)Ammonia formed per hour is \",\n",
+      "print \"%.6f\"%W,\n",
+      "print \"Kg\"\n",
+      "\n",
+      "#(ii)Calculation of velocity\n",
+      "N1 = 0.25*N-(0.25*N*0.15) #Kg mole of N2 after reactor\n",
+      "N2 = 0.75*N-(0.75*N*0.15) #Kg mole of H2 after reactor\n",
+      "N3 = 0.25*N*2*0.15 #Kg mole of NH3 after reactor\n",
+      "Nt = N1+N2+N3 #total Kg moles after reactor\n",
+      "y1NH3 = N3/Nt #mole fraction of NH3 after reactor\n",
+      "y1N2 = N1/Nt #mole fraction of N2 after reactor\n",
+      "y1H2 = N2/Nt #mole fraction of H2 after reactor\n",
+      "T1cm = (TcN2*y1N2)+(TcH2*y1H2) \n",
+      "P1cm = (PcN2*y1N2)+(PcH2*y1H2) \n",
+      "T1 = 448 #in K\n",
+      "P1 = 550 #in atm\n",
+      "T1rm = T1/T1cm \n",
+      "P1rm = P1/P1cm\n",
+      "#From Figure A.2.2\n",
+      "Zm1 = 1.38\n",
+      "V1 = (Zm1*Nt*R*T1)/P1\n",
+      "d = 5*(10**-2)#diameter of pipe\n",
+      "v = V1/((math.pi/4)*(d**2)*3600)\n",
+      "print \"ii)Velocity in pipe is \",\n",
+      "print \"%.6f\"%v,\n",
+      "print \"m/s\"\n",
+      "#end\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "i)Ammonia formed per hour is  113.659704 Kg\n",
+        "ii)Velocity in pipe is  1.075261 m/s\n"
+       ]
+      }
+     ],
+     "prompt_number": 9
+    }
+   ],
+   "metadata": {}
+  }
+ ]
+}
\ No newline at end of file