Initial commit

12 files changed, 3570 insertions, 0 deletions

diff --git a/Elementary_Heat_Power_by_H_L_Solberg/1-Matter_and_energy.ipynb b/Elementary_Heat_Power_by_H_L_Solberg/1-Matter_and_energy.ipynb
new file mode 100644
index 0000000..7a4ff72
--- /dev/null
+++ b/Elementary_Heat_Power_by_H_L_Solberg/1-Matter_and_energy.ipynb
@@ -0,0 +1,388 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 1: Matter and energy"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 1.10: Example_10.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"g=32.2 //ft/s^2\n",
+"p1=100 //psia\n",
+"p2=25 //psia\n",
+"v2=2 //cu ft\n",
+"//calculations\n",
+"W=p1*144*v2 - p2*144*v2\n",
+"//results\n",
+"printf('Work done = %d ft-lb',W)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 1.11: Example_11.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"g=32.2 //ft/s^2\n",
+"n=100 //rpm\n",
+"p1=100 //psia\n",
+"p2=25 //psia\n",
+"v2=2 //cu ft\n",
+"//calculations\n",
+"W=p1*144*v2 - p2*144*v2\n",
+"Hp=W*n/33000\n",
+"//results\n",
+"printf('Horsepower developed = %.1f hp',Hp)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 1.12: Example_12.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"P=50 //hp\n",
+"m=30 //lb\n",
+"E=19000 //Btu/lb\n",
+"//calculations\n",
+"eta= P*2545/(m*E) *100\n",
+"//results\n",
+"printf('Efficiency = %.1f percent',eta)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 1.1: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"g=32.2 //ft/s^2\n",
+"m=500 //lb\n",
+"rate=10 //ft/s^2\n",
+"//calculations\n",
+"F1=m/g *rate\n",
+"ms=m/g\n",
+"F2=ms*rate\n",
+"//results\n",
+"printf('Force in case 1 = %.1f lbf',F1)\n",
+"printf('\n Force in case 2 = %.1f lbf',F2)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 1.2: Example_2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"g=32.2 //ft/s^2\n",
+"g2=32.0 //ft/s^2\n",
+"rate=10 //ft/s^2\n",
+"w1=500 //lbf\n",
+"//calculations\n",
+"fd1=w1*g2/g\n",
+"F=fd1/g2 *rate\n",
+"ms=w1/g\n",
+"F2=ms*rate\n",
+"//results\n",
+"printf('Net weight of body in case 1 = %.1f lbf',F)\n",
+"printf('\n Force in case 2 = %.1f lbf',F2)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 1.3: Example_3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"g=32.174 //ft/s^2\n",
+"m=500 //lbm\n",
+"rate=10 //ft/s^2\n",
+"//calculations\n",
+"F=1/g *m*rate\n",
+"//results\n",
+"printf('Force required = %.1f lbf',F)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 1.4: Example_4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"g1=32.174 //ft/s^2\n",
+"gc=g1\n",
+"g2=30 //ft/s^2\n",
+"m=100 //lbm\n",
+"//calculations\n",
+"w1=g1/gc *m\n",
+"w2=g2/gc *m\n",
+"//results\n",
+"printf('Weight in case 1 = %d lbf',w1)\n",
+"printf('\n Weight in case 2 = %.1f lbf',w2)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 1.5: Example_5.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"ge=32.174 //ft/s^2\n",
+"gm=5.47 //ft/s^2\n",
+"we=50 //lbm\n",
+"//calculations\n",
+"wm=we*gm/ge\n",
+"//results\n",
+"printf('In case a, it will weigh the same, weight  = %d lbm',we)\n",
+"printf('\n In case b, weight = %.1f lbf',wm)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 1.6: Example_6.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"g=32.2 //ft/s^2\n",
+"p1=100 //psig\n",
+"p2=29.0 //in of Hg\n",
+"//calculations\n",
+"BP=p2*0.491\n",
+"AP=BP+p1\n",
+"//results\n",
+"printf('Absolute pressure = %.2f psia',AP)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 1.7: Example_7.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"g=32.2 //ft/s^2\n",
+"Pb=29.5 //in of Hg\n",
+"Pv=10 //in of Hg\n",
+"//calculations\n",
+"AP=(Pb-Pv)*0.491\n",
+"//results\n",
+"printf('Absoulte pressure = %.2f psia',AP)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 1.8: Example_8.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"g=32.2 //ft/s^2\n",
+"v1=1 //cu ft\n",
+"p1=100 //psia\n",
+"//calculations\n",
+"v2=2*v1\n",
+"W=144*p1*(v2-v1)\n",
+"//results\n",
+"printf('Work done = %d ft-lb',W)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 1.9: Example_9.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"g=32.2 //ft/s^2\n",
+"v1=1 //cu ft\n",
+"p1= 100 //psia\n",
+"p2=50 //psia\n",
+"v2=3 //cu ft\n",
+"//calculations\n",
+"pa=(p1+p2)/2\n",
+"W=pa*(v2-v1)*144\n",
+"//results\n",
+"printf('Work done = %d ft-lb',W)"
+   ]
+   }
+],
+"metadata": {
+		  "kernelspec": {
+		   "display_name": "Scilab",
+		   "language": "scilab",
+		   "name": "scilab"
+		  },
+		  "language_info": {
+		   "file_extension": ".sce",
+		   "help_links": [
+			{
+			 "text": "MetaKernel Magics",
+			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+			}
+		   ],
+		   "mimetype": "text/x-octave",
+		   "name": "scilab",
+		   "version": "0.7.1"
+		  }
+		 },
+		 "nbformat": 4,
+		 "nbformat_minor": 0
+}
diff --git a/Elementary_Heat_Power_by_H_L_Solberg/10-Drafts_fans_blowers_and_compressors.ipynb b/Elementary_Heat_Power_by_H_L_Solberg/10-Drafts_fans_blowers_and_compressors.ipynb
new file mode 100644
index 0000000..14df2c4
--- /dev/null
+++ b/Elementary_Heat_Power_by_H_L_Solberg/10-Drafts_fans_blowers_and_compressors.ipynb
@@ -0,0 +1,303 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 10: Drafts fans blowers and compressors"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 10.1: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"hb=29 //in of Hg\n",
+"sg=0.491\n",
+"Ra=53.3\n",
+"Ta=460+40 //R\n",
+"Tg=540+460 //R\n",
+"H=300 //ft\n",
+"gam=62.4 //lb/cu ft\n",
+"//calculations\n",
+"pb=hb*sg*144\n",
+"rhoa=pb/(Ra*Ta)\n",
+"rhog=pb/(Ra*Tg)\n",
+"dp=H*(rhoa-rhog)\n",
+"D=dp/(gam)\n",
+"//results\n",
+"printf('Theoretical draft = %.1f psf',dp)\n",
+"printf('\n Draft = %.2f ft H2O',D)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 10.2: Example_2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"md=15 //lb per lb of coal\n",
+"x=0.1\n",
+"mss=1 //basis\n",
+"rea=29 //in of Hg\n",
+"sg=0.491\n",
+"R=53.3\n",
+"T=540+460 //R\n",
+"V=25 //fps\n",
+"gam=0.038 //lb/ft^3\n",
+"//calculations\n",
+"m=mss-mss*x+md\n",
+"ms=m\n",
+"rhog=rea*0.491*144/(R*T)\n",
+"A=ms/(gam*V)\n",
+"//results\n",
+"printf('stack area = %.1f sq ft',A)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 10.3: Example_3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"p=144*29*0.491 //psf\n",
+"R=53.3\n",
+"T=70+460 //R\n",
+"gamw=62.4 //lb/ft^3\n",
+"gama=0.073 //lb/ft^3\n",
+"hw=3/12 //ft\n",
+"hw2=3.5/12 //ft\n",
+"hv=32.2 //ft/s^2\n",
+"ms=9 //lb\n",
+"g=32.2 //ft/s^2\n",
+"//calculations\n",
+"rhoa=p/(R*T)\n",
+"hs=hw*gamw/gama\n",
+"ht=hw2*gamw/gama\n",
+"hv=ht-hs\n",
+"V=sqrt(2*g*hv)\n",
+"msv=ms*V*60\n",
+"mm=msv*gama\n",
+"airhp= ht*mm/33000\n",
+"//results\n",
+"printf('Velocity head = %d ft of air',hv)\n",
+"printf('\n velocity of air in the duct = %.1f fps',V)\n",
+"printf('\n volume = %d cu ft per min',msv)\n",
+"printf('\n Mass flow rate = %d lb/min',mm)\n",
+"printf('\n Air hp = %.1f hp',airhp)\n",
+"disp('The answers in the textbook are a bit different due to rounding off error in the textbook. Please use a calculator')"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 10.4: Example_4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"A2=9 //sq ft\n",
+"p2=3/12 *62.4 //psf\n",
+"p1=-1/12 *62.4 //psf\n",
+"ms=20000 //cfm\n",
+"A1=16 //sq ft\n",
+"gam=0.075 //lb/ft^3\n",
+"g=32.2 //ft/s^2\n",
+"inp=17 //hp\n",
+"//calculations\n",
+"V2=ms/60 *1/A2\n",
+"V1=ms/60 *1/A1\n",
+"ht=(p2-p1)/gam +(V2^2 -V1^2)/(2*g)\n",
+"airhp=ht*ms*gam/33000\n",
+"eta=airhp/inp *100\n",
+"//results\n",
+"printf('Total head = %.1f ft of air',ht)\n",
+"printf('\n Air hp = %.1f hp',airhp)\n",
+"printf('\n Effifciency = %.1f percent',eta)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 10.5: Example_5.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"n1=400 //rpm\n",
+"mv1=10000 //lb\n",
+"mv2=15000 //lb\n",
+"h1=2 //in of water\n",
+"hp1=4 //hp\n",
+"//calculations\n",
+"n2=mv2/mv1 *n1\n",
+"h2=h1*(n2/n1)^2\n",
+"hp2=hp1 *(n2/n1)^3\n",
+"//results\n",
+"printf('The speed = %d rpm',n2)\n",
+"printf('\n The pressure = %.1f in of water',h2)\n",
+"printf('\n Power = %.1f hp',hp2)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 10.6: Example_6.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"m=100000 //lb/hr\n",
+"p1=1 //psia\n",
+"x=0.8\n",
+"p2=14.7 //psia\n",
+"t2=300 //F\n",
+"//calculations\n",
+"disp('from table A3 and A2')\n",
+"h2=1192.8 //Btu/lb\n",
+"hf=69.7 //Btu/lb\n",
+"hfg=1036.3 //Btu/lb\n",
+"h1=hf+x*hfg\n",
+"W=h2-h1\n",
+"power=m*W\n",
+"hp=power/2545\n",
+"//results\n",
+"printf('Power required = %d hp',hp)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 10.7: Example_7.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"p1=14.7 //psia\n",
+"t1=60 //F\n",
+"p2=60 //psia\n",
+"t2=440 //F\n",
+"m=10 //lb/sec\n",
+"//calculations\n",
+"disp('From mollier charts,')\n",
+"h2=216.3 //Btu/lb\n",
+"h1=124.3 //Btu/lb\n",
+"W21=h2-h1\n",
+"power=W21*m\n",
+"hp=power*3600/2545\n",
+"cp=0.237\n",
+"W212=cp*(t2-t1)\n",
+"power2=W212*m\n",
+"hp2=power2*3600/2545\n",
+"//results\n",
+"printf('Power required = %d hp',hp)\n",
+"printf('\n Power required = %d hp',hp2)\n",
+"printf('\n Work done = %.1f Btu/lb',W212)"
+   ]
+   }
+],
+"metadata": {
+		  "kernelspec": {
+		   "display_name": "Scilab",
+		   "language": "scilab",
+		   "name": "scilab"
+		  },
+		  "language_info": {
+		   "file_extension": ".sce",
+		   "help_links": [
+			{
+			 "text": "MetaKernel Magics",
+			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+			}
+		   ],
+		   "mimetype": "text/x-octave",
+		   "name": "scilab",
+		   "version": "0.7.1"
+		  }
+		 },
+		 "nbformat": 4,
+		 "nbformat_minor": 0
+}
diff --git a/Elementary_Heat_Power_by_H_L_Solberg/11-Feed_water_heaters_and_condensers.ipynb b/Elementary_Heat_Power_by_H_L_Solberg/11-Feed_water_heaters_and_condensers.ipynb
new file mode 100644
index 0000000..d19e3c7
--- /dev/null
+++ b/Elementary_Heat_Power_by_H_L_Solberg/11-Feed_water_heaters_and_condensers.ipynb
@@ -0,0 +1,162 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 11: Feed water heaters and condensers"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.1: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"m1=1000 //lb/hr\n",
+"m2=5000 //lb/hr\n",
+"m3=3000 //lb/hr\n",
+"//calculations\n",
+"disp('From mollier charts,')\n",
+"h5=196.16 //Btu/lb\n",
+"h1=38.04 //Btu/lb\n",
+"h2=67.97 //Btu/lb\n",
+"h3=117.89 //Btu/lb\n",
+"h4=1156.3 //Btu/lb\n",
+"m4=(m1*h1+m2*h2+m3*h3-(m1+m2+m3)*h5)/(h5-h4)\n",
+"//results\n",
+"printf('Pounds of steam entering the heater = %d lb/hr',m4)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.2: Example_2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"P1=100 //psia\n",
+"T1=400 //F\n",
+"T2=70 //F\n",
+"//calculations\n",
+"disp('From mollier charts,')\n",
+"h1=1227.6 //Btu/lb\n",
+"h2=298.4 //Btu/lb\n",
+"h3=279.9 //Btu/lb\n",
+"h4=38.04 //Btu/lb\n",
+"m1=(h3-h4)/(h1-h2)\n",
+"//results\n",
+"printf('Mass of steam required = %.2f lb steam per lb water',m1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.3: Example_3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"h0=1260 //Btu/lb\n",
+"msr=15 //lb\n",
+"m4=15 //lb per hr per kw\n",
+"t2=80 //F\n",
+"t3=60 //F\n",
+"//calculations\n",
+"h1=h0-3413/msr\n",
+"disp('from mollier charts,')\n",
+"h4=58 //Btu/lb\n",
+"dt=t2-t3\n",
+"m3=m4*(h1-h4)/dt\n",
+"//results\n",
+"printf('enthalpy of steam entering the condenser = %d Btu/lb',h1)\n",
+"printf('\n mass of cooling water = %d lb per hr per kw',m3)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.4: Example_4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"m4=8*1000000 //lb per hr\n",
+"dt=12 //F\n",
+"//calculations\n",
+"disp('from mollier charts,')\n",
+"dh4=950 //Btu/lb\n",
+"m3=m4*(dh4)/dt\n",
+"//results\n",
+"printf('\n mass of cooling water = %.3e lb per hr',m3)"
+   ]
+   }
+],
+"metadata": {
+		  "kernelspec": {
+		   "display_name": "Scilab",
+		   "language": "scilab",
+		   "name": "scilab"
+		  },
+		  "language_info": {
+		   "file_extension": ".sce",
+		   "help_links": [
+			{
+			 "text": "MetaKernel Magics",
+			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+			}
+		   ],
+		   "mimetype": "text/x-octave",
+		   "name": "scilab",
+		   "version": "0.7.1"
+		  }
+		 },
+		 "nbformat": 4,
+		 "nbformat_minor": 0
+}
diff --git a/Elementary_Heat_Power_by_H_L_Solberg/12-The_Gas_turbine_power_plant.ipynb b/Elementary_Heat_Power_by_H_L_Solberg/12-The_Gas_turbine_power_plant.ipynb
new file mode 100644
index 0000000..54491f3
--- /dev/null
+++ b/Elementary_Heat_Power_by_H_L_Solberg/12-The_Gas_turbine_power_plant.ipynb
@@ -0,0 +1,113 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 12: The Gas turbine power plant"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 12.1: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"T1=80 //F\n",
+"T2=460 //F\n",
+"T3=1300 //F\n",
+"T4=780 //F\n",
+"//calculations\n",
+"disp('from mollier charts,')\n",
+"h1=129.1 //Btu/lb\n",
+"h2 = 221.2 //Btu/lb\n",
+"h3= 438.8 //Btu/lb\n",
+"h4 = 301.5 //Btu/lb\n",
+"wcom=h2-h1\n",
+"wcob=h3-h2\n",
+"wtur=h3-h4\n",
+"eta=(wtur-wcom)/wcob *100\n",
+"//results\n",
+"printf('\n work done by compressor = %.1f btu input as work per lb of air compressed',wcom)\n",
+"printf('\n Heat supplied in the combustor = %.1f Btu supplied per lb of air ',wcob)\n",
+"printf('\n work done in the turbine = %.1f Btu output as work per lb of air',wtur)\n",
+"printf('\n Cycle efficiency = %.1f percent',eta)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 12.2: Example_2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"T1=80 //F\n",
+"T2=460 //F\n",
+"T=700 //F\n",
+"T3=1300 //F\n",
+"T4=780 //F\n",
+"//calculations\n",
+"disp('from mollier charts,')\n",
+"h1=129.1 //Btu/lb\n",
+"h2 = 221.2 //Btu/lb\n",
+"h3= 438.8 //Btu/lb\n",
+"h4 = 301.5 //Btu/lb\n",
+"wcom=h2-h1\n",
+"wcob=h3-h2\n",
+"wtur=h3-h4\n",
+"output=-wcom+wtur\n",
+"h=281.1 //Btu/lb\n",
+"Q=h3-h\n",
+"eff=output/Q *100\n",
+"//results\n",
+"printf('\n Heat supplied in the combustor = %.1f Btu supplied per lb of air ',Q)\n",
+"printf('\n Cycle efficiency = %.1f percent',eff)"
+   ]
+   }
+],
+"metadata": {
+		  "kernelspec": {
+		   "display_name": "Scilab",
+		   "language": "scilab",
+		   "name": "scilab"
+		  },
+		  "language_info": {
+		   "file_extension": ".sce",
+		   "help_links": [
+			{
+			 "text": "MetaKernel Magics",
+			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+			}
+		   ],
+		   "mimetype": "text/x-octave",
+		   "name": "scilab",
+		   "version": "0.7.1"
+		  }
+		 },
+		 "nbformat": 4,
+		 "nbformat_minor": 0
+}
diff --git a/Elementary_Heat_Power_by_H_L_Solberg/13-Mechanical_Refrigeration.ipynb b/Elementary_Heat_Power_by_H_L_Solberg/13-Mechanical_Refrigeration.ipynb
new file mode 100644
index 0000000..17f0ba3
--- /dev/null
+++ b/Elementary_Heat_Power_by_H_L_Solberg/13-Mechanical_Refrigeration.ipynb
@@ -0,0 +1,134 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 13: Mechanical Refrigeration"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 13.1: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"\n",
+"disp('From mollier diagram from ammonia, values are found')\n",
+"disp('part a')\n",
+"h1=65 //Btu/lb\n",
+"printf('enthalpy in case a = %d Btu/lb',h1)\n",
+"h2=99 //Btu/lb\n",
+"v2=0.93 //ft^3/lb\n",
+"printf('\n In case 2, enthalpy and specific volume are %d Btu/lb and %.2f ft^3/lb respectively',h2,v2)\n",
+"h3=583 //Btu/lb\n",
+"v3=8.8 //ft^3/lb\n",
+"s3=1.275\n",
+"printf('\n In case 3, enthalpy, specific volume and entropy are %d Btu/lb, %.2f ft^3/lb and %.3f respectively',h3,v3,s3)\n",
+"h4=720 //Btu/lb\n",
+"v4=10.4 //ft^3/lb\n",
+"s4=1.50\n",
+"printf('\n In case 4, enthalpy, specific volume and entropy are %d Btu/lb, %.2f ft^3/lb and %.3f respectively',h4,v4,s4)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 13.2: Example_2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"mr=3 //lb\n",
+"mj=5 //lb\n",
+"t2=67 //F\n",
+"t1=60 //lb\n",
+"ihp=7.25\n",
+"//calculations\n",
+"disp('From mollier charts,')\n",
+"h4=709 //Btu/b\n",
+"h3=618 //Btu/lb\n",
+"energyin=ihp*2545/60\n",
+"energyout=mr*(h4-h3) + mj*(t2-t1)\n",
+"//results\n",
+"printf('Energy in = %.1f Btu/min',energyin)\n",
+"printf('\n Energy out = %.1f Btu/min',energyout)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 13.3: Example_3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"mr=3 //lb\n",
+"hp=10 //hp\n",
+"//calculations\n",
+"h3=618 //Btu/lb\n",
+"h1=131 //Btu/lb\n",
+"Qe=mr*(h3-h1)\n",
+"work=hp*2545/60\n",
+"cop=Qe/work\n",
+"//results\n",
+"printf('Coefficient of performance = %.2f',cop)"
+   ]
+   }
+],
+"metadata": {
+		  "kernelspec": {
+		   "display_name": "Scilab",
+		   "language": "scilab",
+		   "name": "scilab"
+		  },
+		  "language_info": {
+		   "file_extension": ".sce",
+		   "help_links": [
+			{
+			 "text": "MetaKernel Magics",
+			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+			}
+		   ],
+		   "mimetype": "text/x-octave",
+		   "name": "scilab",
+		   "version": "0.7.1"
+		  }
+		 },
+		 "nbformat": 4,
+		 "nbformat_minor": 0
+}
diff --git a/Elementary_Heat_Power_by_H_L_Solberg/2-Fuels_and_Combustion.ipynb b/Elementary_Heat_Power_by_H_L_Solberg/2-Fuels_and_Combustion.ipynb
new file mode 100644
index 0000000..59b6918
--- /dev/null
+++ b/Elementary_Heat_Power_by_H_L_Solberg/2-Fuels_and_Combustion.ipynb
@@ -0,0 +1,547 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 2: Fuels and Combustion"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.10: Example_10.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"Gf=11.57 //lb per lb of fuel\n",
+"H=4.4/100\n",
+"M=13.5/100\n",
+"mr=700\n",
+"mf=10000\n",
+"mc=1 //lb\n",
+"//calculations\n",
+"pro=M+9*H\n",
+"mrf=mr/mf\n",
+"Aa=Gf+pro+mrf-mc\n",
+"At=8.83\n",
+"ea=(Aa-At)/At *100\n",
+"//results\n",
+"printf('Excess air = %.1f percent',ea)\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.11_a: Example_11.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clear\n",
+"//Initialization of variables\n",
+"Gf=11.57 //lb per lb of fuel\n",
+"tg=500 //F\n",
+"ta=70 //F\n",
+"//calculations\n",
+"Q1=0.24*Gf*(tg-ta)\n",
+"//results\n",
+"printf('Heat loss = %d Btu per lb of fuel',Q1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.11_b: Example_12.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"Co=0.1\n",
+"Co2=14.1\n",
+"Cb=0.646\n",
+"//calculations\n",
+"Q2=Co/(Co+Co2) *Cb*10160\n",
+"//results\n",
+"printf('Heat loss = %d Btu per lb of fuel',Q2)\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.11c: Example_13.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"mf=10000 //lb \n",
+"mr=700 //lb\n",
+"Cr=0.2\n",
+"//calculations\n",
+"Q3=mr*Cr/mf *14600\n",
+"//results\n",
+"printf('Heat loss = %d Btu per lb of fuel',Q3)\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.11d: Example_14.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"M=0.135\n",
+"tg=500 //F\n",
+"ta=70 //F\n",
+"//calculations\n",
+"Q4=M*(1089+0.46*tg-ta)\n",
+"//results\n",
+"printf('Heat loss = %.1f Btu per lb of fuel',Q4)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.11e: Example_15.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"Per=0.044 //percentage\n",
+"tg=500 //F\n",
+"ta=70 //F\n",
+"//calculations\n",
+"Q5=9*Per*(1089+0.46*tg-ta)\n",
+"//results\n",
+"printf('Heat loss = %.1f Btu per lb of fuel',Q5)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.1: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"x1=0.135\n",
+"x2=0.056\n",
+"veca=[32.5 48.4 5.6 13.5]\n",
+"B1=11788\n",
+"//calculations\n",
+"vecb=veca/(1-x1)\n",
+"vecc=veca/(1-x1-x2)\n",
+"B2=B1/(1-x1)\n",
+"B3=B1/(1-x1-x2)\n",
+"vecb(4)=0\n",
+"vecc(3)=0\n",
+"vecc(4)=0\n",
+"//results\n",
+"printf('In Moisture free case, ')\n",
+"format('v',6);vecb\n",
+"disp(vecb)\n",
+"printf('In Moisture and Ash free case, ')\n",
+"format('v',6);vecc\n",
+"disp(vecc)\n",
+"printf('Energy in Moisture free case  = %d Btu per lb',B2)\n",
+"printf('\n Energy in Moisture and ash free case  = %d Btu per lb',B3)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.2: Example_2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"y1=13.5\n",
+"x1=0.135\n",
+"x2=0.056\n",
+"veca=[66 1.5 1.1 5.6 5.9 19.9]\n",
+"//calculations\n",
+"vecb=[veca y1]\n",
+"vecb(5) = vecb(5) - 1/9*y1\n",
+"vecb(6) = vecb(6) - 8/9*y1\n",
+"vecc=vecb/(1-x1)\n",
+"vecd=vecb/(1-x1-x2)\n",
+"vecd(4)=0\n",
+"vecd(7)=0\n",
+"vecc(7)=0\n",
+"s1=sum(vecc)\n",
+"s2=sum(vecd)\n",
+"//results\n",
+"printf('With moisture as a separate item, ')\n",
+"format ('v',6);vecb\n",
+"disp(vecb)\n",
+"printf('In Moisture free case, ')\n",
+"format('v',4);vecc\n",
+"disp(vecc)\n",
+"printf('In Moisture and Ash free case, ')\n",
+"format('v',5);vecd\n",
+"disp(vecd)\n",
+"printf('Total Mositure free content = %.1f percent',s1)\n",
+"printf('\n Total Mositure and ash free content = %.1f percent',s2)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.3: Example_3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"H=5.9\n",
+"O=19.9\n",
+"H2=4.4\n",
+"O2=7.9\n",
+"//calculations\n",
+"Ha1=H-O/8\n",
+"Ha2=H2-O2/8\n",
+"//results\n",
+"printf('Available hydrogen in case 1 = %.1f percent by weight',Ha1)\n",
+"printf('\n Available hydrogen in case 1 = %.1f percent by weight',Ha2)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.4: Example_4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"H1=0.059\n",
+"O1=0.199\n",
+"H2=0.044\n",
+"O2=0.079\n",
+"C=0.66\n",
+"S=0.011\n",
+"//calculations\n",
+"Qh1= 14600*C+62000*(H1-O1/8)+4050*S\n",
+"Qh2=14600*C+62000*(H2-O2/8)+4050*S\n",
+"//results\n",
+"printf('Heating value in case 1 = %d Btu/lb ',Qh1)\n",
+"printf('\n Heating value in case 2 = %d Btu/lb ',Qh2)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.5: Example_5.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"H1=0.059\n",
+"O1=0.199\n",
+"C=0.66\n",
+"S=0.011\n",
+"//calculations\n",
+"Qh1= 11.52*C+34.56*(H1-O1/8)+4.32*S\n",
+"//results\n",
+"printf('Theoretical air required = %.2f lb of air per lb of coal ',Qh1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.6: Example_6.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"mf=10000 //lb\n",
+"mr=700 //lb\n",
+"Cr=0.20\n",
+"Cco2=14.1\n",
+"Co2=5.1\n",
+"Cco=0.1\n",
+"Cf=0.66\n",
+"//calculations\n",
+"Cn2=100-(Cco2+Co2+Cco)\n",
+"Ci=mf*Cf\n",
+"Ca=mr*Cr\n",
+"Cb=(Ci-Ca)/mf\n",
+"Cb2=((mf*Cf)-mr*Cr)/(mf)\n",
+"veca=[Cco2 Co2 Cco Cn2]\n",
+"vecb=veca\n",
+"vecb(1)=vecb(1) *44\n",
+"vecb(2)=vecb(2) *32\n",
+"vecb(3)=vecb(3) *28\n",
+"vecb(4)=vecb(4) *28\n",
+"sumvec=sum(vecb)\n",
+"Lbc=Cco2*12 + Cco*12\n",
+"Gc=sumvec/Lbc\n",
+"Gf=Gc*Cb\n",
+"//results\n",
+"printf('Carbon in the dry products combustion = %.3f lb per lb of fuel',Cb)\n",
+"printf('\n In case 2, Carbon in the dry products combustion = %.3f lb per lb of fuel',Cb2)\n",
+"printf('\n Dry gaseous products of combstion per lb of coal = %.2f lb ',Gf)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.7: Example_7.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"mf=10000 //lb\n",
+"mr=700 //lb\n",
+"Cr=0.20\n",
+"Cco2=14.1\n",
+"Co2=5.1\n",
+"Cco=0.1\n",
+"Cf=0.66\n",
+"//calculations\n",
+"Cn2=100-(Cco2+Co2+Cco)\n",
+"Ci=mf*Cf\n",
+"Ca=mr*Cr\n",
+"Cb=(Ci-Ca)/mf\n",
+"Cb2=((mf*Cf)-mr*Cr)/(mf)\n",
+"veca=[Cco2 Co2 Cco Cn2]\n",
+"vecb=veca\n",
+"vecb(1)=vecb(1) *44\n",
+"vecb(2)=vecb(2) *32\n",
+"vecb(3)=vecb(3) *28\n",
+"vecb(4)=vecb(4) *28\n",
+"Cbb1=Cb*Cco*12/(Cco2*12 + Cco*12)\n",
+"Cbb2= Cb*(veca(3) /(veca(3) + veca(1)))\n",
+"//results\n",
+"printf('In case 1, Carbon burned per lb of fuel = %.5f lb per lb of fuel',Cbb1)\n",
+"printf('\n In case 2, Carbon burned per lb of fuel = %.5f lb per lb of fuel',Cbb2)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.8: Example_8.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"H=4.4/100\n",
+"M=13.5/100\n",
+"H2=0.059\n",
+"//calculations\n",
+"pro=M+9*H\n",
+"pro2=9*H2\n",
+"//results\n",
+"printf('In case 1, watervapor present in products = %.3f lb',pro)\n",
+"printf('\n In case 2, watervapor present in products = %.3f lb',pro2)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.9: Example_9.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"Gf=11.57 //lb per lb of fuel\n",
+"H=4.4/100\n",
+"M=13.5/100\n",
+"mr=700\n",
+"mf=10000\n",
+"mc=1 //lb\n",
+"//calculations\n",
+"pro=M+9*H\n",
+"mrf=mr/mf\n",
+"Aa=Gf+pro+mrf-mc\n",
+"//results\n",
+"printf('Actual air supplied = %.2f lb of air supplied per lb of fuel',Aa)"
+   ]
+   }
+],
+"metadata": {
+		  "kernelspec": {
+		   "display_name": "Scilab",
+		   "language": "scilab",
+		   "name": "scilab"
+		  },
+		  "language_info": {
+		   "file_extension": ".sce",
+		   "help_links": [
+			{
+			 "text": "MetaKernel Magics",
+			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+			}
+		   ],
+		   "mimetype": "text/x-octave",
+		   "name": "scilab",
+		   "version": "0.7.1"
+		  }
+		 },
+		 "nbformat": 4,
+		 "nbformat_minor": 0
+}
diff --git a/Elementary_Heat_Power_by_H_L_Solberg/3-Internal_Combustion_Engines.ipynb b/Elementary_Heat_Power_by_H_L_Solberg/3-Internal_Combustion_Engines.ipynb
new file mode 100644
index 0000000..f982fbe
--- /dev/null
+++ b/Elementary_Heat_Power_by_H_L_Solberg/3-Internal_Combustion_Engines.ipynb
@@ -0,0 +1,420 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 3: Internal Combustion Engines"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.10: Example_10.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"bore=3 //in\n",
+"str=4 //in\n",
+"rpm=3000 //rpm\n",
+"air=110 //cu ft per min\n",
+"//calculations\n",
+"pdv=bore*bore*%pi*str*2*bore/4\n",
+"pde=pdv*rpm /2\n",
+"req=air*1728\n",
+"eff=req/pde *100\n",
+"//results\n",
+"printf('Volumetric efficiency = %.1f percent',eff)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.11: Example_11.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"x1=11.5\n",
+"x2=4.1\n",
+"x3=0.4\n",
+"x4=2.3\n",
+"x5=0.2\n",
+"x6=81.5\n",
+"yc=0.842\n",
+"yh=0.158\n",
+"basis=1\n",
+"bhp=100\n",
+"burn=8.9 //gal/hr\n",
+"sg=0.731\n",
+"Qh=20750 //Btu/lbm\n",
+"rate=66 //gpm\n",
+"ex=1100 //F\n",
+"air=70 //F\n",
+"cp=0.254\n",
+"h2=4330\n",
+"h4=62000\n",
+"h5=23700\n",
+"//calculations\n",
+"c1=x1*44\n",
+"c2=x2*28\n",
+"c3=x3*32\n",
+"c4=x4*2\n",
+"c5=x5*16\n",
+"c6=x6*28\n",
+"summ=c1+c2+c3+c4+c5+c6\n",
+"carbon=x1*12 + x2*12+x5*12\n",
+"hydrogen=x4*2+x5*4\n",
+"lbdrygas=summ/carbon *yc\n",
+"lbfuel=carbon/yc\n",
+"lbH=lbfuel*yh\n",
+"lbH2=lbH-hydrogen\n",
+"lb3=lbH2*9\n",
+"lbwater=lb3/lbfuel\n",
+"lbair=lbdrygas+lbwater-basis\n",
+"bsfc=burn*sg*8.33/bhp\n",
+"fuelmin=bsfc*bhp/60\n",
+"energy=2545/bsfc\n",
+"per1=energy/Qh\n",
+"Ec=rate*8.33*10\n",
+"Eclb=Ec/fuelmin\n",
+"per2=Eclb/Qh\n",
+"dryloss=(ex-air)*cp*lbdrygas\n",
+"per3=dryloss/Qh\n",
+"hv2=h2*c2/lbfuel\n",
+"hv4=h4*c4/lbfuel\n",
+"hv5=h5*c5/lbfuel\n",
+"hv=hv2+hv4+hv5\n",
+"per4=hv/Qh\n",
+"eh2=lbwater*(1066+0.5*ex-air)\n",
+"per5=eh2/Qh\n",
+"rad=1017 \n",
+"per6=rad/Qh\n",
+"//results\n",
+"printf('Air supplied per lb of fuel = %.1f lb air per lb fuel',lbair)\n",
+"printf('\n Percentage of energy supplied utilized in Btu = %.2f percent',per1*100)\n",
+"printf('\n Percentage of energy absorbed by coolant = %.2f percent',per2*100)\n",
+"printf('\n Energy lost in sensible heat = %.2f percent',per3*100)\n",
+"printf('\n Energy supplied in combustiles in exhaust = %.2f percent',per4*100)\n",
+"printf('\n Energy supplied in water formed by combustion = %.2f percent',per5*100)\n",
+"printf('\n Energy supplied unaccounted for = %.2f percent',per6*100)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.1: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"re=6\n",
+"k=1.4\n",
+"//calculations\n",
+"nt=1-1/re^(k-1)\n",
+"ntt=nt*100\n",
+"//results\n",
+"printf('Thermal efficiency = %.1f percent',ntt)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.2: Example_2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"d=3.25 //in\n",
+"str=4 //in\n",
+"v=6 //cu in\n",
+"//calculations\n",
+"Dp=d^2 *%pi*str/4\n",
+"r=(Dp+v)/v\n",
+"//results\n",
+"printf('Compression ratio = %.2f',r)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.3: Example_3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"per=20\n",
+"Dp=100\n",
+"//calculations\n",
+"r=Dp/per +1\n",
+"//results\n",
+"printf('Compression ratio = %d ',r)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.4: Example_4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"r=16\n",
+"rc=4\n",
+"k=1.4\n",
+"//calculations\n",
+"etat=1-1/r^(k-1) *((rc^k -1)/(k*(rc-1)))\n",
+"eta=etat*100\n",
+"//results\n",
+"printf('Thermal efficiency = %.1f percent',eta)\n",
+"disp('The answer is a bit different due to rounding off error in the textbook')"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.5: Example_5.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"F=200 //lb\n",
+"area=1.65 //sq. in\n",
+"len=3.5 //in\n",
+"//calculations\n",
+"ord=area/len\n",
+"mep=ord*F\n",
+"//results\n",
+"printf('MEP of an engine = %.1f psi',mep)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.6: Example_6.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"Pi=90 //psi\n",
+"L=5/12 //ft\n",
+"r=5 //in\n",
+"x=1.5 //ft\n",
+"rpm=1500 //rpm\n",
+"//calculations\n",
+"A=%pi*x*x\n",
+"N=rpm*4/2\n",
+"Ihp=Pi*L*A*N/33000\n",
+"//results\n",
+"printf('IHP of cylinder = %d',Ihp)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.7: Example_7.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"r=4 //ft\n",
+"n=300 //rpm\n",
+"F=60 //lb\n",
+"//calculations\n",
+"Bhp=2*%pi*r*F*n/33000\n",
+"//results\n",
+"printf('Bhp of the engine = %.1f',Bhp)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.8: Example_8.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"C=1/4000\n",
+"F=125 //lb\n",
+"n=3500 //rpm\n",
+"//calculations\n",
+"Bhp=F*n*C\n",
+"//results\n",
+"printf('Bhp developed by the engine = %.1f',Bhp)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.9: Example_9.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"r=1.75 //ft\n",
+"F1=72 //lb\n",
+"F2=24 //lb\n",
+"n=500 //rpm\n",
+"m=1.8 //lb\n",
+"mi=15 //min\n",
+"Qh=20000 //Btu/lb\n",
+"//calculations\n",
+"Bhp=2*%pi*r*F1*n/33000\n",
+"Fhp=2*%pi*r*F2*n/33000\n",
+"Ihp=Bhp+Fhp\n",
+"Fc=m*60/mi\n",
+"Bsfc=Fc/Bhp\n",
+"Isfc=Fc/Ihp\n",
+"etam=Bhp/Ihp *100\n",
+"etabt=Bhp*2545/(Fc*Qh)\n",
+"etait=Ihp*2545/(Fc*Qh)\n",
+"//results\n",
+"printf('Thermal efficiency = %d percent',etam)\n",
+"printf('\n Brake thermal effficiency = %.1f percent',etabt*100)\n",
+"printf('\n Indicated thermal effficiency = %.1f percent',etait*100)"
+   ]
+   }
+],
+"metadata": {
+		  "kernelspec": {
+		   "display_name": "Scilab",
+		   "language": "scilab",
+		   "name": "scilab"
+		  },
+		  "language_info": {
+		   "file_extension": ".sce",
+		   "help_links": [
+			{
+			 "text": "MetaKernel Magics",
+			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+			}
+		   ],
+		   "mimetype": "text/x-octave",
+		   "name": "scilab",
+		   "version": "0.7.1"
+		  }
+		 },
+		 "nbformat": 4,
+		 "nbformat_minor": 0
+}
diff --git a/Elementary_Heat_Power_by_H_L_Solberg/5-Steam_Generation.ipynb b/Elementary_Heat_Power_by_H_L_Solberg/5-Steam_Generation.ipynb
new file mode 100644
index 0000000..0c02a5b
--- /dev/null
+++ b/Elementary_Heat_Power_by_H_L_Solberg/5-Steam_Generation.ipynb
@@ -0,0 +1,496 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 5: Steam Generation"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.10: Example_10.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"m=150000 //lb\n",
+"P1=1000 //psia\n",
+"Ts=900 //F\n",
+"Tf=200 //F\n",
+"//calculations\n",
+"disp('From mollier charts,')\n",
+"h2=1448.2 //Btu/lb\n",
+"hf=167.99 //Btu/lb\n",
+"correc=2.2 //Btu/lb\n",
+"hc=hf+correc\n",
+"Q=m*(h2-hc)\n",
+"//results\n",
+"printf('Heat absorption = %d Btu/hr',Q)\n",
+"disp('The answer is a bit different due to rounding off error in textbook')"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.11: Example_11.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"m=150000 //lb\n",
+"P1=1000 //psia\n",
+"Ts=900 //F\n",
+"Tf=200 //F\n",
+"//calculations\n",
+"disp('From mollier charts,')\n",
+"h2=1448.2 //Btu/lb\n",
+"hf=167.99 //Btu/lb\n",
+"correc=2.2 //Btu/lb\n",
+"hc=hf+correc\n",
+"Q=m*(h2-hc)\n",
+"output=Q/1000\n",
+"//results\n",
+"printf('Output of the steam generating unit = %d kB/hr',output)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.12: Example_12.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"m=150000 //lb\n",
+"P1=1000 //psia\n",
+"Ts=900 //F\n",
+"Tf=200 //F\n",
+"m2=21000 //lb\n",
+"HV=12000 //Btu/lb\n",
+"//calculations\n",
+"disp('From mollier charts,')\n",
+"h2=1448.2 //Btu/lb\n",
+"hf=167.99 //Btu/lb\n",
+"correc=2.2 //Btu/lb\n",
+"hc=hf+correc\n",
+"Q=m*(h2-hc)\n",
+"output=Q\n",
+"inpu=m2*HV\n",
+"eta=output/inpu\n",
+"//results\n",
+"printf('Efficiency of the steam generating unit = %.1f percent',eta*100)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.13: Example_13.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"hv=11780 //Btu/lb\n",
+"steam=55000 //lb/hr\n",
+"coal=6480 //lb\n",
+"x1=0.66\n",
+"x2=0.044\n",
+"x3=0.079\n",
+"x4=0.015\n",
+"x5=0.11\n",
+"z1=14.5\n",
+"z2=0.2\n",
+"z3=4.4\n",
+"z4=80.9\n",
+"xash=0.076\n",
+"xmois=0.115\n",
+"yc=0.21\n",
+"refuse=622 //lb/hr\n",
+"cp=0.24\n",
+"tg=400 //F\n",
+"ta=70 //F\n",
+"Qco=10160 //Btu/lb\n",
+"Qc=14600 //Btu/lb\n",
+"//calculations\n",
+"disp('From steam tables,')\n",
+"hf=269.6 //Btu/lbm\n",
+"hfg=1.5 //Btu/lbm\n",
+"h1=hf+hfg\n",
+"h2=1196.5\n",
+"Qb=h2-h1\n",
+"h3=1407.7 //Btu/lbm\n",
+"Qs=h3-h2\n",
+"h4=h3-h1\n",
+"out=steam*h4/1000\n",
+"eff=steam*h4/(coal*hv)\n",
+"//Energy balance\n",
+"Ci=coal*x1\n",
+"Cr=refuse*yc\n",
+"Cb=(Ci-Cr)/coal\n",
+"lbt= z1*44+z2*28+z3*32+z4*28\n",
+"lbC=z1*12+z2*12\n",
+"dry=lbt/lbC *Cb\n",
+"loss1=dry*cp*(tg-ta)\n",
+"loss2=z2*12/(lbC) *Cb*Qco\n",
+"loss3=Cr*Qc/coal\n",
+"loss4=xmois*(1089+0.46*tg-ta)\n",
+"loss5=x2*9*(1089+0.46*tg-ta)\n",
+"loss6=steam*h4/coal\n",
+"//results\n",
+"printf('Heat absorbed in the boiler = %.2f Btu per lb of steam generated',Qb)\n",
+"printf('\n Heat absorbed in the superheater = %.2f Btu/lb of steam',Qs)\n",
+"printf('\n Heat absorbed in steam generating = %.2f Btu/lb of steam generated',h4)\n",
+"printf('\n Output of steam generating unit = %d kB',out)\n",
+"printf('\n Efficiency of steam generating unit = %.1f percent',eff*100)\n",
+"printf('\n Carbon burned to CO and CO2 = %.2f lb of C per lb of fuel',Cb)\n",
+"printf('\n Dry products of combustion = %.2f lb per lb of fuel',dry)\n",
+"printf('\n Loss due to sensible heat in dry gaseous products of combustion = %d Btu/lb of fuel',loss1)\n",
+"printf('\n Loss due to CO in dry products of combustion = %.1f Btu/lb of fuel',loss2)\n",
+"printf('\n Loss due to C in refuse = %.1f Btu/lb of fuel',loss3)\n",
+"printf('\n Loss due to evaporating moisture in fuel = %.1f Btu/lb of fuel',loss4)\n",
+"printf('\n Loss due to water vapor formed from H = %.1f Btu/lb of fuel',loss5)\n",
+"printf('\n Energy absorbed in generating steam = %d Btu/lb of fuel',loss6)\n",
+"disp('The answers are a bit different due to rounding off error in the textbook')"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.1: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"x=0.98\n",
+"vg=26.80\n",
+"vf=0.01672\n",
+"//calculations\n",
+"vx=x*vg+(1-x)*vf\n",
+"//results\n",
+"printf('Specific volume of wet steam = %.6f cu ft per lb',vx)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.2: Example_2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"hf=167.99 //Btu/lb\n",
+"hg=4.5 //Btu/lb\n",
+"//calculations\n",
+"hc=hf+hg\n",
+"//results\n",
+"printf('Enthalpy of water = %.1f Btu/lb',hc)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.3: Example_3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"x=0.97\n",
+"hg=1187.2 //Btu/lb\n",
+"hf=298.40 //Btu/lb\n",
+"hfg=888.8 //Btu/lb\n",
+"//calculations\n",
+"hx1=x*hg+(1-x)*hf\n",
+"hx2=hf+x*hfg\n",
+"hx3=hg-(1-x)*hfg\n",
+"//results\n",
+"printf('\n In case 1, enthalpy = %.1f Btu/lb',hx1)\n",
+"printf('\n In case 2, enthalpy = %.1f Btu/lb',hx2)\n",
+"printf('\n In case 3, enthalpy = %.1f Btu/lb',hx3)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.4: Example_4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"h1=1172 //Btu/lb\n",
+"hf1=355.36 //Btu/lb\n",
+"hfg1=843 //Btu/lb\n",
+"//calculations\n",
+"h2=h1\n",
+"x1= (h2-hf1)/hfg1\n",
+"//results\n",
+"printf('Quality of steam = %.1f percent',x1*100)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.5: Example_5.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"P=200 //psia\n",
+"x=0.95\n",
+"m=1//lb\n",
+"//calculations\n",
+"disp('From mollier chart,')\n",
+"hx=1156 //Btu/lb\n",
+"sx=1.495 //Btu/lb F\n",
+"//results\n",
+"printf('Enthalpy = %d Btu/lb',hx)\n",
+"printf('\n entropy = %.3f Btu/lb F',sx)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.6: Example_6.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"P=200 //psia\n",
+"T=600 //F\n",
+"m=1 //lb\n",
+"//calculations\n",
+"disp('From mollier chart,')\n",
+"hx=1322 //Btu/lb\n",
+"sx=1.676 //Btu/lb F\n",
+"//results\n",
+"printf('Enthalpy = %d Btu/lb',hx)\n",
+"printf('\n entropy = %.3f Btu/lb F',sx)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.7: Example_7.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"P=200 //psia\n",
+"T=260 //F\n",
+"//calculations\n",
+"disp('From mollier chart,')\n",
+"hx=1174 //Btu/lb\n",
+"x1=2.8\n",
+"y1=100-x1\n",
+"//results\n",
+"printf('Quality = %.1f percent',y1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.8: Example_8.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"P=200 //psia\n",
+"T=500 //F\n",
+"//calculations\n",
+"disp('From mollier chart,')\n",
+"hi=1269 //Btu/lb\n",
+"hf=1063 //Btu/lb\n",
+"dh=hi-hf\n",
+"y1=91\n",
+"//results\n",
+"printf('Quality = %.1f percent',y1)\n",
+"printf('\n Change in enthalpy = %d Btu/lb',dh)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.9: Example_9.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"P=200 //psia\n",
+"Ts=260 //F\n",
+"Tf=220 //F\n",
+"m=10000 //lb\n",
+"Pc=20 //psia\n",
+"//calculations\n",
+"disp('From mollier charts,')\n",
+"hf=188 //Btu/lb\n",
+"h2=1172 //Btu/lb\n",
+"Q=m*(h2-hf)\n",
+"//results\n",
+"printf('Heat absorption = %d Btu/hr',Q)"
+   ]
+   }
+],
+"metadata": {
+		  "kernelspec": {
+		   "display_name": "Scilab",
+		   "language": "scilab",
+		   "name": "scilab"
+		  },
+		  "language_info": {
+		   "file_extension": ".sce",
+		   "help_links": [
+			{
+			 "text": "MetaKernel Magics",
+			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+			}
+		   ],
+		   "mimetype": "text/x-octave",
+		   "name": "scilab",
+		   "version": "0.7.1"
+		  }
+		 },
+		 "nbformat": 4,
+		 "nbformat_minor": 0
+}
diff --git a/Elementary_Heat_Power_by_H_L_Solberg/6-Steam_power_plant_cycles.ipynb b/Elementary_Heat_Power_by_H_L_Solberg/6-Steam_power_plant_cycles.ipynb
new file mode 100644
index 0000000..8c343df
--- /dev/null
+++ b/Elementary_Heat_Power_by_H_L_Solberg/6-Steam_power_plant_cycles.ipynb
@@ -0,0 +1,189 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 6: Steam power plant cycles"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.1: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"P1=200 //psia\n",
+"T1=600 //F\n",
+"P2=2 //psia\n",
+"J=778\n",
+"//calculations\n",
+"disp('from mollier charts,')\n",
+"h1=1322 //Btu/lb\n",
+"h2=974 //Btu/lb\n",
+"vf2=0.01623 //cu ft per lb\n",
+"hf2=94 //Btu/lb\n",
+"t2=126 //F\n",
+"Wtj=h1-h2\n",
+"Qout=h2-hf2\n",
+"Wp=(P1-P2)*vf2\n",
+"Wpj=Wp/J\n",
+"h3=hf2+Wpj\n",
+"Qin=h1-h3\n",
+"etat=((h1-h2)-Wpj)/(h1-(hf2+Wpj))\n",
+"eta=((h1-h2))/(h1-(hf2))\n",
+"//results\n",
+"printf('Efficiency of rankine cycle = %.1f percent',etat*100)\n",
+"printf('\n Efficiency of rankine cycle neglecting boiler feed pump = %.1f percent',eta*100)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.2: Example_2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"B=70 //F\n",
+"P1=140 //psia\n",
+"x=0.986\n",
+"P2=14.7 //psia\n",
+"ms=2000 //lb/hr\n",
+"Ihp=80\n",
+"//calculations\n",
+"disp('From mollier charts,')\n",
+"hc=38 //Btu/lb\n",
+"hf=324.82 //Btu/lb\n",
+"hfg=868.2 //Btu/lb\n",
+"h1=hf+x*hfg\n",
+"Qin=ms*(h1-hc)\n",
+"eta=Ihp*2545*100/(Qin)\n",
+"Qw=Ihp*2545\n",
+"Qr=Qin-Qw\n",
+"per=Qr/Qin *100\n",
+"//results\n",
+"printf('Heat input to the boiler = %d Btu/hr',Qin)\n",
+"printf('\n Cycle efficiency = %.1f percent',eta)\n",
+"printf('\n Heat rejected to waste = %d Btu/hr or %.1f percent of Qin',Qr,per)\n",
+"disp('The answer is a bit different due to rounding off error in textbook')"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.3: Example_3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"B=70 //F\n",
+"P1=140 //psia\n",
+"x=0.986\n",
+"P2=14.7 //psia\n",
+"ms=2000 //lb/hr\n",
+"Ihp=80\n",
+"//calculations\n",
+"disp('From mollier charts,')\n",
+"hc=180 //Btu/lb\n",
+"hf=324.82 //Btu/lb\n",
+"hfg=868.2 //Btu/lb\n",
+"h1=hf+x*hfg\n",
+"Qin=ms*(h1-hc)\n",
+"eta=Ihp*2545*100/(Qin)\n",
+"Qw=Ihp*2545\n",
+"Qr=Qin-Qw\n",
+"per=Qr/Qin *100\n",
+"//results\n",
+"printf('Heat input to the boiler = %d Btu/hr',Qin)\n",
+"printf('\n Cycle efficiency = %.2f percent',eta)\n",
+"printf('\n Heat rejected to waste = %d Btu/hr or %.2f percent of Qin',Qr,per)\n",
+"disp('The answer is a bit different due to rounding off error in textbook')"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.4: Example_4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"m=1.24 //lb\n",
+"HV=11300 //Btu/lb\n",
+"//calculations\n",
+"HR=m*HV\n",
+"eff=3413/HR\n",
+"//results\n",
+"printf('Plant heat rate = %d Btu/kw hr',HR)\n",
+"printf('\n Overall efficiency = %.1f percent',eff*100)"
+   ]
+   }
+],
+"metadata": {
+		  "kernelspec": {
+		   "display_name": "Scilab",
+		   "language": "scilab",
+		   "name": "scilab"
+		  },
+		  "language_info": {
+		   "file_extension": ".sce",
+		   "help_links": [
+			{
+			 "text": "MetaKernel Magics",
+			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+			}
+		   ],
+		   "mimetype": "text/x-octave",
+		   "name": "scilab",
+		   "version": "0.7.1"
+		  }
+		 },
+		 "nbformat": 4,
+		 "nbformat_minor": 0
+}
diff --git a/Elementary_Heat_Power_by_H_L_Solberg/7-Steam_turbines.ipynb b/Elementary_Heat_Power_by_H_L_Solberg/7-Steam_turbines.ipynb
new file mode 100644
index 0000000..aad6b7a
--- /dev/null
+++ b/Elementary_Heat_Power_by_H_L_Solberg/7-Steam_turbines.ipynb
@@ -0,0 +1,343 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 7: Steam turbines"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.10: Example_10.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"P=200 //psia\n",
+"T=540 //F\n",
+"pow=1000 //kw\n",
+"ms=16000 //lb/hr\n",
+"//calculations\n",
+"disp('From mollier charts,')\n",
+"h1=1290 //Btu/hr\n",
+"h2=940 //Btu/hr\n",
+"dh=h1-h2\n",
+"hf2=83 //Btu/lb\n",
+"etat=(h1-h2)/(h1-hf2)\n",
+"act=pow*3413/(ms*(h1-hf2))\n",
+"etae=act/etat\n",
+"//results\n",
+"printf('Ideal thermal efficiency = %.1f percent',etat*100)\n",
+"printf('\n Actual thermal efficiency = %.1f percent',act*100)\n",
+"printf('\n Engine efficiency = %.1f percent',etae*100)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.1: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"P1=200 //psia\n",
+"T1=500 //psia\n",
+"m=1 //lb /s\n",
+"P4=140 //psia\n",
+"P11=1 //psia\n",
+"x=0.808\n",
+"//calculations\n",
+"disp('From mollier charts,')\n",
+"h1=1268.9 //Btu/lb\n",
+"h4=1234.7 //Btu/lb\n",
+"V4=223.8*sqrt(h1-h4)\n",
+"v4=3.584 //cu ft/lb\n",
+"A4=m*v4/V4\n",
+"h11=907.4 //Btu/lb\n",
+"V11=223.8*sqrt(h1-h11)\n",
+"vf=0.01614 //cu ft/lb\n",
+"vg=333.6 //cu ft/lb\n",
+"vfg=vg-vf\n",
+"v11=x*vg\n",
+"A11=m*v11/V11\n",
+"//results\n",
+"printf('Area of nozzle = %.5f sq ft',A4)\n",
+"printf('\n Area of nozzle = %.4f sq ft',A11)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.4: Example_4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"P1=200 //psia\n",
+"T1=500 //F\n",
+"P2=1 //psia\n",
+"alpha=20 //degrees\n",
+"n=3600\n",
+"g=32.2 //ft/s^2\n",
+"//calculations\n",
+"disp('From mollier charts,')\n",
+"V1=4240 //fps\n",
+"Vb=V1*cosd(alpha) /2\n",
+"R=Vb*60/(n*2*%pi)\n",
+"work=1/32.2 *(V1*cosd(alpha))*Vb\n",
+"eff=work/(V1^2 /(2*g)) *100\n",
+"//results\n",
+"printf('Blade velocity = %d fps',Vb)\n",
+"printf('\n Blade radius = %.1f ft',R)\n",
+"printf('\n Work done = %d ft-lb per lb of steam',work)\n",
+"printf('\n Blade efficiency = %.1f percent',eff)\n",
+"disp('The answers are a bit different due to rounding off error in textbook.')"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.5: Example_5.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"P1=200 //psia\n",
+"T1=500 //F\n",
+"P2=1 //psia\n",
+"alpha=20 //degrees\n",
+"n=3600\n",
+"g=32.2 //ft/s^2\n",
+"Vb=1200 //fps\n",
+"//calculations\n",
+"disp('From mollier charts,')\n",
+"V1=4240 //fps\n",
+"V1x=3980 //fps\n",
+"V2x=-1580 //fps\n",
+"work=1/32.2 *(V1x - V2x)*Vb\n",
+"eff=work/(V1^2 /(2*g)) *100\n",
+"//results\n",
+"printf('\n Work done = %d ft-lb per lb of steam',work)\n",
+"printf('\n Blade efficiency = %.1f percent',eff)\n",
+"disp('The answers are a bit different due to rounding off error in textbook.')"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.6: Example_6.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"P1=200 //psia\n",
+"T1=500 //F\n",
+"P2=1 //psia\n",
+"alpha=20 //degrees\n",
+"n=3600\n",
+"g=32.2 //ft/s^2\n",
+"//calculations\n",
+"disp('From mollier charts,')\n",
+"V1=2450 //fps\n",
+"Vb=V1*cosd(alpha) /2\n",
+"R=Vb*60/(n*2*%pi)\n",
+"work=1/32.2 *(V1*cosd(alpha))*Vb\n",
+"w3=3*work\n",
+"//results\n",
+"printf('Blade velocity = %d fps',Vb)\n",
+"printf('\n Blade radius = %.2f ft',R)\n",
+"printf('\n Work done = %d ft-lb per lb of steam',w3)\n",
+"disp('The answers are a bit different due to rounding off error in textbook.')"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.7: Example_7.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"P1=200 //psia\n",
+"T1=500 //F\n",
+"P2=1 //psia\n",
+"alpha=20 //degrees\n",
+"n=3600\n",
+"g=32.2 //ft/s^2\n",
+"stage=2\n",
+"//calculations\n",
+"disp('From mollier charts,')\n",
+"V1=4240 //fps\n",
+"Vb=V1*cosd(alpha) /(2*stage)\n",
+"R=Vb*60/(n*2*%pi)\n",
+"V1x=3980 //fps\n",
+"V2x=-1990 //fps\n",
+"work1=1/g *(V1x-V2x)*Vb\n",
+"work2=1/g *(-V2x)*Vb\n",
+"wt=work1+work2\n",
+"//results\n",
+"printf('Blade velocity = %d fps',Vb)\n",
+"printf('\n Blade radius = %.2f ft',R)\n",
+"printf('\n Total Work done = %d ft-lb per lb of steam',wt)\n",
+"disp('The answers are a bit different due to rounding off error in textbook.')"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.8: Example_8.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"alpha=20 //degrees\n",
+"n=3600\n",
+"g=32.2 //ft/s^2\n",
+"V1=500 //fps\n",
+"//calculations\n",
+"Vb=V1*cosd(alpha)\n",
+"V1x=Vb\n",
+"work=1/32.2 *(V1x)*Vb\n",
+"//results\n",
+"printf('Blade velocity = %d fps',Vb)\n",
+"printf('\n Work done = %d ft-lb per lb of steam',work)\n",
+"disp('The answers are a bit different due to rounding off error in textbook.')\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.9: Example_9.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"pow=1000 //kw\n",
+"ms=16000 //lb/hr\n",
+"P=200 //psia\n",
+"T=540 //F\n",
+"//calculations\n",
+"disp('From mollier charts,')\n",
+"h1=1290 //Btu/hr\n",
+"h2=940 //Btu/hr\n",
+"dh=h1-h2\n",
+"rate=3413/dh\n",
+"act=ms/pow\n",
+"//results\n",
+"printf('Ideal steam rate = %.2f lb per kw hr',rate)\n",
+"printf('\n Actual steam rate = %d lb per kw hr',act)"
+   ]
+   }
+],
+"metadata": {
+		  "kernelspec": {
+		   "display_name": "Scilab",
+		   "language": "scilab",
+		   "name": "scilab"
+		  },
+		  "language_info": {
+		   "file_extension": ".sce",
+		   "help_links": [
+			{
+			 "text": "MetaKernel Magics",
+			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+			}
+		   ],
+		   "mimetype": "text/x-octave",
+		   "name": "scilab",
+		   "version": "0.7.1"
+		  }
+		 },
+		 "nbformat": 4,
+		 "nbformat_minor": 0
+}
diff --git a/Elementary_Heat_Power_by_H_L_Solberg/8-Steam_engines.ipynb b/Elementary_Heat_Power_by_H_L_Solberg/8-Steam_engines.ipynb
new file mode 100644
index 0000000..b87ae6b
--- /dev/null
+++ b/Elementary_Heat_Power_by_H_L_Solberg/8-Steam_engines.ipynb
@@ -0,0 +1,122 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 8: Steam engines"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.1: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"area1=2.7\n",
+"len=3.4\n",
+"scale=60\n",
+"area2=2.75\n",
+"dia=12 //ft\n",
+"d2=2.5 //ft\n",
+"L=15/12 //ft\n",
+"n=250 //rpm\n",
+"F=600 //lb\n",
+"r=3 //ft\n",
+"//calculations\n",
+"Ah=dia^2 *%pi/4\n",
+"Ac=(dia^2 -d2^2)*%pi/4\n",
+"Pih=area1/len *scale\n",
+"Pic=area2/len *scale\n",
+"Hihp=Pih*L*Ah*n/33000\n",
+"Cihp=Pic*L*Ac*n/33000\n",
+"Tihp=Hihp+Cihp\n",
+"Bhp=2*%pi*r*F*n/33000\n",
+"Fhp=Tihp-Bhp\n",
+"eff=Bhp/Tihp *100\n",
+"//results\n",
+"printf('Ihp = %.1f ihp',Tihp)\n",
+"printf('\n Bhp = %.1f bhp',Bhp)\n",
+"printf('\n Fhp = %.1f fhp',Fhp)\n",
+"printf('\n Efficiency = %.1f percent',eff)\n",
+"disp('The answer is a bit different due to rounding off error in the textbook.')"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.2: Example_2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"Ihp=101.1\n",
+"Bhp=85.7\n",
+"md=3000 //lb/hr\n",
+"h1=1172 //Btu/hr\n",
+"h2=180 //Btu/hr\n",
+"h3=1025 //Btu/hr\n",
+"//calculations\n",
+"eta1=Ihp*2545/(md*(h1-h2)) *100\n",
+"eta2=Bhp*2545/(md*(h1-h2)) *100\n",
+"etat=(h1-h3)/(h1-h2) *100\n",
+"engeff=eta1/etat *100\n",
+"rate1= md/Ihp\n",
+"rate2=md/Bhp\n",
+"h22=h1-2545/rate1\n",
+"//results\n",
+"printf('Actual thermal efficiency based upon Ihp = %.2f lb per ihp hr',eta1)\n",
+"printf('\n Actual thermal efficiency based upon Bhp = %.2f lb per ihp hr',eta2)\n",
+"printf('\n Ideal thermal efficiency = %.2f percent ',etat)\n",
+"printf('\n Engine efficiency = %.1f percent',engeff)\n",
+"printf('\n Steam rate = %.2f lb per ihp hr',rate1)\n",
+"printf('\n Steam rate = %.2f lb per bhp hr',rate2)\n",
+"printf('\n Enthalpy of exhaust steam = %d Btu/lb of steam',h22)"
+   ]
+   }
+],
+"metadata": {
+		  "kernelspec": {
+		   "display_name": "Scilab",
+		   "language": "scilab",
+		   "name": "scilab"
+		  },
+		  "language_info": {
+		   "file_extension": ".sce",
+		   "help_links": [
+			{
+			 "text": "MetaKernel Magics",
+			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+			}
+		   ],
+		   "mimetype": "text/x-octave",
+		   "name": "scilab",
+		   "version": "0.7.1"
+		  }
+		 },
+		 "nbformat": 4,
+		 "nbformat_minor": 0
+}
diff --git a/Elementary_Heat_Power_by_H_L_Solberg/9-Pumps.ipynb b/Elementary_Heat_Power_by_H_L_Solberg/9-Pumps.ipynb
new file mode 100644
index 0000000..b61182a
--- /dev/null
+++ b/Elementary_Heat_Power_by_H_L_Solberg/9-Pumps.ipynb
@@ -0,0 +1,353 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 9: Pumps"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.10: Example_10.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"sr=2\n",
+"//calculations\n",
+"hr=sr^2\n",
+"capr=sr\n",
+"hpr=sr^3\n",
+"//results\n",
+"printf('head is %d times the original',hr)\n",
+"printf('\n capacity is  %d times the original',capr)\n",
+"printf('\n power is  %d times the original',hpr)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.1: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"h=200 //ft\n",
+"gam=64 //lb per cu ft\n",
+"//calculations\n",
+"P=h*gam/144\n",
+"//results\n",
+"printf('Pressure = %.1f psi',P)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.2: Example_2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"P=20 //psi\n",
+"gam=62.4 //lb per cu ft\n",
+"//calculations\n",
+"h=P*144/gam\n",
+"//results\n",
+"printf('height = %.1f ft',h)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.3: Example_3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"h=3/12 //ft\n",
+"gam=63.4 //lb per cu ft\n",
+"gam2=0.075 //lb per cu ft\n",
+"//calculations\n",
+"P=h*gam\n",
+"h2=P/gam2\n",
+"//results\n",
+"printf('Air height required = %d ft of air',h2)\n",
+"disp('The answer is a bit different due to roundoff error in textbook.')"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.4: Example_4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"dif=4 //in\n",
+"gam=62.4 //lb per cu ft\n",
+"density=13.6 //g/cc\n",
+"//calculations\n",
+"pv=dif*1/12 *density*gam/144 - dif/12 *gam/144\n",
+"hv=pv*144/gam\n",
+"//results\n",
+"printf('velocity pressure = %.2f psi',pv)\n",
+"printf('\n velocity head = %.1f ft of water ',hv)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.5: Example_5.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"dif=4 //in\n",
+"gam=62.4 //lb per cu ft\n",
+"gam2=0.08 //lb per cu ft\n",
+"//calculations\n",
+"pv=dif*1/12 *gam/144\n",
+"hv=pv*144/gam2\n",
+"//results\n",
+"printf('velocity pressure = %.3f psi',pv)\n",
+"printf('\n velocity head = %.1f ft of air ',hv)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.6: Example_6.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"hw=3/12 //ft\n",
+"gam1=62.4 //lb/ft^3\n",
+"gam2=0.07 //lb/ft^3\n",
+"g=32.2 //ft/s^2\n",
+"//calculations\n",
+"p=hw*gam1\n",
+"hg=p/gam2\n",
+"V=sqrt(2*g*hg)\n",
+"//results\n",
+"printf('velocity of gas = %.1f fps',V)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.7: Example_7.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"h=4 //in\n",
+"den=13.6 //g/cc\n",
+"Ar=1/9\n",
+"A1=12 //sq in\n",
+"gam=62.4 //lb/ft^3\n",
+"g=32.2 //ft/s^2\n",
+"//calculations\n",
+"dh=(h*den-h)/12\n",
+"Vr=1/Ar\n",
+"V22=2*g*dh/(1-Ar^2)\n",
+"V2=sqrt(V22)\n",
+"A2=A1*Ar\n",
+"v2=1/gam\n",
+"ms=A2*V2/(v2*144)\n",
+"//results\n",
+"printf('Flow rate of water = %.1f lb/sec',ms)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.8: Example_8.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"mdot=8000 //lb/min\n",
+"A1=1 //sq ft\n",
+"A2=3/4 //sq ft\n",
+"P2=50 //psi\n",
+"P1=10 //psi\n",
+"gam=62.4 //lb/ft^3\n",
+"y2=-2 //ft\n",
+"y1=-4 //ft\n",
+"g=32.2 //ft/s^2\n",
+"eff=0.7\n",
+"//calculations\n",
+"v=1/gam\n",
+"cap=mdot/8.33\n",
+"V1=mdot*v/A1 /60\n",
+"V2=mdot*v/A2 /60\n",
+"ht= (y2-y1) + (V2^2 -V1^2)/(2*g) + (P2-P1)*144/gam\n",
+"Hhp=mdot*ht/33000\n",
+"Php=Hhp/eff\n",
+"//results\n",
+"printf('Capacity = %d gpm',cap)\n",
+"printf('\n Total dynamic head = %.1f ft',ht)\n",
+"printf('\n Hydraulic hp = %.1f hp',Hhp)\n",
+"printf('\n pump hp = %.1f hp',Php)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.9: Example_9.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"clear\n",
+"//Initialization of variables\n",
+"z=12 //ft\n",
+"gam1=62.4 //lb/ft^3\n",
+"sg=0.8\n",
+"P2=100 //psia\n",
+"P1=-10 //psia\n",
+"mm=10000 //lb/min\n",
+"//calculations\n",
+"gam2=sg*gam1\n",
+"p2g=P2*144/(gam2) +z\n",
+"p1g=P1*144*0.491/(gam2)\n",
+"ht=p2g-p1g\n",
+"Hhp=mm*ht/33000\n",
+"//results\n",
+"printf('Total dynamic head = %.1f ft of oil',ht)\n",
+"printf('\n Hydraulic hp = %.1f hp',Hhp)"
+   ]
+   }
+],
+"metadata": {
+		  "kernelspec": {
+		   "display_name": "Scilab",
+		   "language": "scilab",
+		   "name": "scilab"
+		  },
+		  "language_info": {
+		   "file_extension": ".sce",
+		   "help_links": [
+			{
+			 "text": "MetaKernel Magics",
+			 "url": "https://github.com/calysto/metakernel/blob/master/metakernel/magics/README.md"
+			}
+		   ],
+		   "mimetype": "text/x-octave",
+		   "name": "scilab",
+		   "version": "0.7.1"
+		  }
+		 },
+		 "nbformat": 4,
+		 "nbformat_minor": 0
+}