Merge pull request #1 from prashantsinalkar/masterHEADmaster

Initial commit

15 files changed, 3010 insertions, 0 deletions

diff --git a/Water_and_Wastewater_Engineering_by_G_M_Fair/10-Groundwater_collection.ipynb b/Water_and_Wastewater_Engineering_by_G_M_Fair/10-Groundwater_collection.ipynb
new file mode 100644
index 0000000..bb90556
--- /dev/null
+++ b/Water_and_Wastewater_Engineering_by_G_M_Fair/10-Groundwater_collection.ipynb
@@ -0,0 +1,98 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 10: Groundwater collection"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 10.1: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"w1=1000//ft\n",
+"w2=2000//ft\n",
+"r=700//gpm\n",
+"d=10//days\n",
+"q=2//ft\n",
+"u=1.87*[(3.4*10^-5)/(3.2*10^4)]*(d^6/d)//ft\n",
+"W=7.94//ft\n",
+"p=114.6*(7*10^2)*W/(3.2*10^4)//ft\n",
+"U=1.87*[(3.4*10^-5)/(3.2*10^4)]*(4*d^6/d)//ft\n",
+"Wu=6.55//ft\n",
+"P=114.6*(7*10^2)*Wu/(3.2*10^4)//ft\n",
+"R=54//ft\n",
+"//CALCULATIONS\n",
+"W1=R+p+P//ft\n",
+"D=R+q*p//ft\n",
+"//RESULTS\n",
+"printf('the expected drawndown the first well is pumped at a rate=% f ft',W1)\n",
+"printf('the drawdown in each well all the three are pupped at a rate=% f ft',D)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 10.3: Example_2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"g=20//ft\n",
+"k=10^-1//cm/sec\n",
+"g1=3.28*10^-3//fps\n",
+"w=2//ft\n",
+"w1=30//ft\n",
+"//CALCULATIONS\n",
+"Q=(1/2)*(g1)*[(g^2)-(2^2)]/(w1)//cfs\n",
+"//RESULTS\n",
+"printf('the flow into a foot of gallery=% f cfs',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/Water_and_Wastewater_Engineering_by_G_M_Fair/11-Surface_water_collections.ipynb b/Water_and_Wastewater_Engineering_by_G_M_Fair/11-Surface_water_collections.ipynb
new file mode 100644
index 0000000..e395a9d
--- /dev/null
+++ b/Water_and_Wastewater_Engineering_by_G_M_Fair/11-Surface_water_collections.ipynb
@@ -0,0 +1,182 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 11: Surface water collections"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.1: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"s=20//mph\n",
+"t=90//min\n",
+"w=1.31//ft\n",
+"h=7.5//miles\n",
+"h1=0.22//ft\n",
+"t1=1100//min\n",
+"t2=6.0//min\n",
+"p=32.2//ft\n",
+"l=5.12//length\n",
+"l1=2.8//length\n",
+"p1=1400//ft\n",
+"d=73//depth\n",
+"h3=2.06//ft\n",
+"e=173.0//ft\n",
+"hi=0.2//ft\n",
+"//CALCULATIONS\n",
+"W=s*w//mph\n",
+"hs=h1*[(W)^2/p]^0.53*h^0.47//ft\n",
+"Ts=t2*(W/p)^0.44*(h/p)^0.28//sec\n",
+"Td=t1*h/(p*Ts)//min\n",
+"Ls=l1/(l*(Ts)^2)//ft\n",
+"D=d/(l*(Ts)^2)//ft\n",
+"H=(W)^2*[h*(1/(p1*d))]//ft\n",
+"hr=h3*l1//ft\n",
+"M=e+hi+hr//ft\n",
+"//RESULTS\n",
+"printf('the overwater wind speed=% f mph',W)\n",
+"printf('the significant wave height=% f ft',hs)\n",
+"printf('the significant wave period=% f sec',Ts)\n",
+"printf('the minimum wind duration required to reach the significant wave height=% f min',Td)\n",
+"printf('the significant wave lenght adn steepness=% f ft',Ls)\n",
+"printf('the reservoir depth ratio=% f ft',D)\n",
+"printf('the wind tide or set up=% f ft',H)\n",
+"printf('the run up =% f ft',hr)\n",
+"printf('the maximum elevation reached by the waves=% f ft',M)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.2: Example_2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"g=264//quartz\n",
+"p=0.39//percent\n",
+"//CALCULATIONS\n",
+"S=(1-p)*(g-1)//in\n",
+"//RESULTS\n",
+"printf('the hydralic gradient and seepage velocity=% f in',S)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.3: Example_3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"w=40//ft\n",
+"k=2*10^-3//cm/sec\n",
+"p=3.28*10^-3//cfs\n",
+"h=6.47*10^5//gpd\n",
+"p1=0.433//ft\n",
+"m=9//ft\n",
+"delh=w/(18*9)//in\n",
+"k1=4.94*10^-4//cm/sec\n",
+"//CALCULATIONS\n",
+"Q=k*p*w*(9/18)//cfs\n",
+"Q1=Q*h//gpd/ft width\n",
+"P=(1-8/18)*w*p1//Psig\n",
+"H=k1/k//in\n",
+"//RESULTS\n",
+"printf('the seepage through each foot width of the foundation=% f gpd/ft/ width',Q1)\n",
+"printf('the excess hydrostatic pressure on the upstream  side of the bottom of the sheet pilling=% f Psig',P)\n",
+"printf('the maximum hydraulic gradient and its relations to the coeeficent=% f in',H)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.4: Example_4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"d=120//ft\n",
+"w=16//ft\n",
+"d1=120/0.8//ft\n",
+"p=60*0.8//ft\n",
+"h=2//ft\n",
+"v=18.74*0.8//ft\n",
+"s=95.23//ft\n",
+"s1=0.8//ft\n",
+"//CALCULATIONS\n",
+"W=d-h*p//ft\n",
+"S=s*s1//ft\n",
+"//RESULTS\n",
+"printf('in succession from the intersection of the upstream slop=% f ft',S)"
+   ]
+   }
+],
+"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/Water_and_Wastewater_Engineering_by_G_M_Fair/12-Water_transmission.ipynb b/Water_and_Wastewater_Engineering_by_G_M_Fair/12-Water_transmission.ipynb
new file mode 100644
index 0000000..d108103
--- /dev/null
+++ b/Water_and_Wastewater_Engineering_by_G_M_Fair/12-Water_transmission.ipynb
@@ -0,0 +1,188 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 12: Water transmission"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 12.1: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"c=100//in\n",
+"a=10//in\n",
+"Q=0.976//ft\n",
+"//CALCULATIONS\n",
+"G=a*Q//ft\n",
+"//RESULTS\n",
+"printf('the graphical basic =% f ft',G)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 12.2: Example_2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"a=27.6//sq ft\n",
+"h=1.37//ft\n",
+"d=1.53*(27.9)^0.38*(1.36)^0.24//ft\n",
+"//CALCULATIONS\n",
+"R=d/4//ft\n",
+"A=(%pi*d^2)/4//sq ft\n",
+"//RESULTS\n",
+"printf('The diameter hydraulics radius and area of the hydraulically equivalent circular conduit=% f sq ft',A)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 12.3: Example_3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"h1=13.5//ft\n",
+"h2=19.0//ft\n",
+"h3=27.5//ft\n",
+"c1=2.0*10^4//ft\n",
+"c2=2.1*10^4//ft\n",
+"c3=2.2*10^4//ft\n",
+"//CALCULATIONS\n",
+"H=h1+h2+h3//ft\n",
+"C=c1+c2+c3//ft\n",
+"//RESULTS\n",
+"printf('the most economical distributions of the available head=% f ft',C)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 12.4: Example_4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"p=60//in\n",
+"h=20//percent\n",
+"a=1000//ft\n",
+"h1=40//percent\n",
+"c=0.5//ft\n",
+"p1=14.3//ft\n",
+"p2=6.1//ft\n",
+"d=11.7*10^-2//ft\n",
+"//CALCULATIONS\n",
+"P=p2/p1//ft\n",
+"D=d*p//ft\n",
+"//RESULTS\n",
+"printf('the air valve with a discharge the change in slop=% f ft',D)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 12.5: Example_5.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"p=90//deg\n",
+"h=48//in\n",
+"p1=100//psig\n",
+"P=(1/2*%pi)*h^2*p1*0.7071//lb\n",
+"r=3000/54-31//ft\n",
+"s=9000//psi\n",
+"l=170//in\n",
+"b=6.5*10^-6//ft\n",
+"w=46//ft\n",
+"w1=1000//ft\n",
+"//CALCULATIONS\n",
+"D=(1/4*%pi)*h^2*p1//lb\n",
+"P=[r]*h^2//lb\n",
+"T=%pi*h*(1/4)*s//lb\n",
+"T1=(1/2)*l//tons\n",
+"Del=b*w*w1//ft per\n",
+"//RESULTS\n",
+"printf('the accorance with unless otherwise stated=% f ft per',Del)"
+   ]
+   }
+],
+"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/Water_and_Wastewater_Engineering_by_G_M_Fair/13-Water_distribution.ipynb b/Water_and_Wastewater_Engineering_by_G_M_Fair/13-Water_distribution.ipynb
new file mode 100644
index 0000000..5b3e2f0
--- /dev/null
+++ b/Water_and_Wastewater_Engineering_by_G_M_Fair/13-Water_distribution.ipynb
@@ -0,0 +1,171 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 13: Water distribution"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 13.2: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"p1=7.8//ft\n",
+"p2=6.0//ft\n",
+"p3=7.4//ft\n",
+"p4=6.5//ft\n",
+"p=7.6//ft\n",
+"h=1.0//ft\n",
+"h1=6.7//ft\n",
+"p5=3.3//ft\n",
+"//CALCULATIONS\n",
+"D=p1-p2//mgd\n",
+"D1=p1-p3//mgd\n",
+"D2=p-p4//mgd\n",
+"D3=p4+h//mgd\n",
+"D4=h1-p5//mgd\n",
+"//RESULTS\n",
+"printf('the demand is taken =% f mgd',D3)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 13.3: Example_2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"w=500//ft\n",
+"p=20//psig\n",
+"h=40//psig\n",
+"h1=1000//in\n",
+"q=1250//ft\n",
+"g=2.308/0.75//ft\n",
+"g1=2.308/1.00//ft\n",
+"s=5200//gpm\n",
+"a=250//gpm\n",
+"//CALCULATIONS\n",
+"H=[h1-(1/2)*(w)]//ft\n",
+"H1=(h-p)*g//percent\n",
+"Q=[q-(1/2)*(w)]//ft\n",
+"Q1=(h-p)*g1//percent\n",
+"S=s/a//gpm\n",
+"//RESULTS\n",
+"printf('the number of standard fire streams=% f gpm',S)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 13.6: Example_3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"h1=2.1*3//ft\n",
+"h2=2.1//ft\n",
+"h=8.4//ft\n",
+"p=1000//ft\n",
+"h3=5.7//ft\n",
+"h4=4.2*3//ft\n",
+"q=4.2//ft\n",
+"s=1.68//ft\n",
+"q1=1.33//ft\n",
+"//CALCULATIONS\n",
+"A=p*h/h2//ft\n",
+"B=p*(h3+h4)/q//ft\n",
+"C=p*(h1+h2)/s//ft\n",
+"//RESULTS\n",
+"printf('the equilent pipe for the Hazen willians coefficent=% f ft',A)\n",
+"printf('the equilent pipe for the Hazen willians coefficent=% f ft',B)\n",
+"printf('the equilent pipe for the Hazen willians coefficent=% f ft',C)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 13.8: Example_4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"d=10//hr\n",
+"p=50000//in\n",
+"a=7.5//mgd\n",
+"w=0.75//mg\n",
+"s=5.03//mg\n",
+"//CALCULATIONS\n",
+"S=s/w//mg\n",
+"P=S-s//mg\n",
+"//RESULTS\n",
+"printf('a steady gravity supply equal to maximum daily=% f mg',P)"
+   ]
+   }
+],
+"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/Water_and_Wastewater_Engineering_by_G_M_Fair/14-Wastewater_flows.ipynb b/Water_and_Wastewater_Engineering_by_G_M_Fair/14-Wastewater_flows.ipynb
new file mode 100644
index 0000000..deaf80a
--- /dev/null
+++ b/Water_and_Wastewater_Engineering_by_G_M_Fair/14-Wastewater_flows.ipynb
@@ -0,0 +1,427 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 14: Wastewater flows"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 14.10: Example_10.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"Q1=30//cfs\n",
+"Q2=16//cfs\n",
+"a=32//sq ft\n",
+"r=1.6//ft\n",
+"i=10^-4//ft\n",
+"n=1.25*10^-2//ft\n",
+"h2=0.50//ft\n",
+"c=3.33//ft\n",
+"h1=5.20//ft\n",
+"l=72//ft\n",
+"s=12320//ft\n",
+"//CALCULATIONS\n",
+"L=s-l//ft\n",
+"//RESULTS\n",
+"printf('the forchheimer s methos =% f ft',L)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 14.11: Example_11.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"q=1.0//cfs\n",
+"g=2.0//percent\n",
+"g1=5.6//percent\n",
+"r=0.015//cfs\n",
+"w=90//percent\n",
+"Q=10*0.9*q//ft\n",
+"p=0.10//ft\n",
+"h=(3.48*g1^1/3)//ft\n",
+"d=p^2/3*100//ft\n",
+"//CALCULATIONS\n",
+"D=h*d//in\n",
+"//RESULTS\n",
+"printf('The maximum depth of flow in the gutter=% f in',D)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 14.1: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"n=0.013//ft\n",
+"s=4.90//ft\n",
+"v=0.590//ft\n",
+"d=0.463//ft\n",
+"w=3.9*10^-2//ft\n",
+"p=1.696//ft\n",
+"//CALCULATIONS\n",
+"V=s*v//fps\n",
+"Q=s*d//cfs\n",
+"N=(w*p)^2*1000//percent\n",
+"//RESULTS\n",
+"printf('the velocity of flow and rate of discharge=% f percent',N)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 14.2: Example_2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"v=1.34//fps\n",
+"s=3.7*10^-3//fps\n",
+"k=0.8//ft\n",
+"r=20//ft\n",
+"k1=0.04//ft\n",
+"v=3.0//fps\n",
+"v1=5.0//fps\n",
+"d=10^-1//ft\n",
+"d1=1.34//ft\n",
+"//CALCULATIONS\n",
+"K=r*k1//ft\n",
+"V=sqrt(r)//times\n",
+"D=d*(v/d1)^2//cm\n",
+"D1=d*(v1/d1)^2//cm\n",
+"//RESULTS\n",
+"printf('the minimum velocity and the gradient at the which coarse quartz=% f cm',D1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 14.3: Example_3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"\n",
+"\n",
+"clc\n",
+"//initialisation of variables\n",
+"v=2.5//fps\n",
+"q=0.873//cfs\n",
+"s=5.20//percent\n",
+"a=0.252//ft\n",
+"r=0.684//ft\n",
+"r1=1.46//ft\n",
+"v1=0.776//ft\n",
+"q1=0.196//ft\n",
+"n=0.78//ft\n",
+"R=0.939//ft\n",
+"//CALCULATIONS\n",
+"V=v1*v//fps\n",
+"Q=q1*q//cfs\n",
+"R1=r1*s//percent\n",
+"Vs=R*v//ft\n",
+"N=n*Vs//fps\n",
+"Qs=a*R*q//cfs\n",
+"N1=n*Qs//cfs\n",
+"//RESULTS\n",
+"printf('the required grades and associated velocity and rates=% f cfs',V)\n",
+"printf('the depth and a grade=% f cfs',Q)\n",
+"printf('the self cleaning flow=% f cfs',N1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 14.4: Example_4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"Q=0.873//cfs\n",
+"s=5.20//percent\n",
+"d=0.161//cfs\n",
+"q1=0.185//ft\n",
+"d2=2.5//ft\n",
+"v=0.91//ft\n",
+"s1=1.70//ft\n",
+"s3=1.46//ft\n",
+"w=0.185//ft\n",
+"d1=0.30//ft\n",
+"v1=0.732//ft\n",
+"//CALCULATIONS\n",
+"q=d/Q//cfs\n",
+"Vs=v*d2//fps\n",
+"Ss=s1*s//percent\n",
+"Va=v1*d2//fps\n",
+"Ss1=s3*s//percent\n",
+"//RESULTS\n",
+"printf('the depth and velocity of flow and the required slop=% f percent',Ss1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 14.5: Example_5.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"d1=0.67//ft\n",
+"h1=2.00//ft\n",
+"h2=4.04//ft\n",
+"hv1=0.062//ft\n",
+"hv2=0.254//ft\n",
+"d=0.19//ft\n",
+"h=0.2//ft\n",
+"h1=0.04//ft\n",
+"q=0.644//ft\n",
+"q1=0.65//ft\n",
+"v=0.92//ft\n",
+"d2=6.5//ft\n",
+"v1=3.69//ft\n",
+"d3=0.542//ft\n",
+"hv3=0.21//ft\n",
+"delv=0.15//ft\n",
+"d4=0.02//ft\n",
+"//CALCULATIONS\n",
+"H=d1+hv1//ft\n",
+"H1=d1+hv2//ft\n",
+"he=h*d//ft\n",
+"hi=d+h1//ft\n",
+"H2=d3+hv3//ft\n",
+"he1=h*delv//ft\n",
+"S=d4+h1//ft\n",
+"//RESULTS\n",
+"printf('the required slope=% f ft',hi)\n",
+"printf('the lower sewer and the invert drop in the transition=% f ft',S)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 14.6: Example_6.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"q=60//cfs\n",
+"D=4//ft\n",
+"w=0.177//ft\n",
+"s=0.59//ft\n",
+"h=4.0//ft\n",
+"d1=1.0//ft\n",
+"v=0.90//ft\n",
+"d1=0.42//ft\n",
+"h1=6.0//ft\n",
+"h2=1.5//ft\n",
+"dl=1.3//ft\n",
+"p=0.41//ft\n",
+"u=0.8//ft\n",
+"u1=3.2//ft\n",
+"y=0.45//ft\n",
+"//CALCULATIONS\n",
+"H=s*D//ft\n",
+"d2=d1*D//ft\n",
+"V=v*D//ft\n",
+"P=p*D//ft\n",
+"D1=y*D//ft\n",
+"//RESULTS\n",
+"printf('the critical depth=% f ft',H)\n",
+"printf('the alternate stages for an energy =% f ft',V)\n",
+"printf('the alternate stages for an energy head=% f ft',P)\n",
+"printf('the lower alternate stage with upper alternate stage=% f ft',D1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 14.7: Example_7.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"d=106//cfs\n",
+"q=400//cfs\n",
+"d1=0.40//cfs\n",
+"w=10//ft\n",
+"//CALCULATIONS\n",
+"D=d/q//cfs\n",
+"D1=d1*w//cfs\n",
+"//RESULTS\n",
+"printf('the water level in this well rises=% f cfs',D1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 14.8: Example_8.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"Q=8.07*10^-2//ft\n",
+"N=0.012//ft\n",
+"d=0.47//ft\n",
+"q=10//ft\n",
+"//CALCULATIONS\n",
+"D=d*q//ft\n",
+"//RESULTS\n",
+"printf('teh water surface in the sewer when it is flowing at maximum capacity=% f ft',D)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 14.9: Example_9.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variable\n",
+"g=sqrt(3)//ft\n",
+"d=5.67//ft\n",
+"//CALCULATIONS\n",
+"C=g*d//ft\n",
+"//RESULTS\n",
+"printf('The rate of propagation of a discontinuous surge=% f ft',C)"
+   ]
+   }
+],
+"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/Water_and_Wastewater_Engineering_by_G_M_Fair/15-Wastewater_collection.ipynb b/Water_and_Wastewater_Engineering_by_G_M_Fair/15-Wastewater_collection.ipynb
new file mode 100644
index 0000000..04024a1
--- /dev/null
+++ b/Water_and_Wastewater_Engineering_by_G_M_Fair/15-Wastewater_collection.ipynb
@@ -0,0 +1,185 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 15: Wastewater collection"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 15.1: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"q=0.25//in\n",
+"Q=0.34//in\n",
+"r=0.76//in\n",
+"v=0.83//in\n",
+"n=0.78//in\n",
+"r1=0.84//in\n",
+"v1=0.70//in\n",
+"w=2//in\n",
+"q1=0.056//in\n",
+"d=0.16//in\n",
+"v2=0.53//in\n",
+"n1=0.80//in\n",
+"d1=0.18//in\n",
+"n2=0.46//in\n",
+"//CALCULATIONS\n",
+"V=v*w//fps\n",
+"N=v1*w//fps\n",
+"V1=v2*w//fps\n",
+"V2=n2*w//fps\n",
+"//RESULTS\n",
+"printf('The one fourth their full flow=% f fps',N)\n",
+"printf('The one enghteenth their full flow=% f fps',V2)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 15.2: Example_2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"v=2.5//fps\n",
+"N=0.015//fps\n",
+"a=(40+27)//in\n",
+"b=(40*27+27*19)/a\n",
+"c=0.440//cfs\n",
+"w=49*0.09/100//cfs\n",
+"g=0.008//percent\n",
+"Q=0.82//cfs\n",
+"r=0.795//cfs\n",
+"t=2.35*1.16//fps\n",
+"d1=113.20-113.03//ft\n",
+"d2=12//ft\n",
+"//CALCULATIONS\n",
+"R=r/Q//cfs\n",
+"D=g*r//in\n",
+"D2=d1*d2//in\n",
+"//RESULTS\n",
+"printf('The required capacity and find the slope size and hydraulic characteristics of the system=% f in',D2)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 15.3: Example_3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"p=20//min\n",
+"N=0.012//in\n",
+"k=2.19//min\n",
+"l=k+1.97//min\n",
+"q=340/(60*3.94)//min\n",
+"r=2.56*0.508//min\n",
+"del=0.42//min\n",
+"j=84.28//min\n",
+"w1=0.92//min\n",
+"//CALCULATIONS\n",
+"r1=r*k//cfs\n",
+"w=p+q//min\n",
+"G=j-del//min\n",
+"S=(G-w1)//min\n",
+"//RESULTS\n",
+"printf('The required capacity and find the slop size and hydraulic=% f min',S)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 15.4: Example_4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"a=42//in\n",
+"d=45//mgd\n",
+"d1=0.75//in\n",
+"s=60//ft\n",
+"p1=9//in\n",
+"p2=8.4//in\n",
+"p3=9//in\n",
+"c1=13*63.6//sq in\n",
+"c2=9*55.4//sq in\n",
+"c3=9.21//sq ft\n",
+"M=d*1.547//cfs\n",
+"v=M/c3//fps\n",
+"g=0.025*32.2//ft/sec^2\n",
+"//CALCULATIONS\n",
+"F=v/sqrt(g*(p1/12))//ft\n",
+"S=s/d1//in\n",
+"//RESULTS\n",
+"printf('the port near the end of the diffuser pipe=% f in',F)"
+   ]
+   }
+],
+"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/Water_and_Wastewater_Engineering_by_G_M_Fair/16-Machinery_and_equipment.ipynb b/Water_and_Wastewater_Engineering_by_G_M_Fair/16-Machinery_and_equipment.ipynb
new file mode 100644
index 0000000..d30d090
--- /dev/null
+++ b/Water_and_Wastewater_Engineering_by_G_M_Fair/16-Machinery_and_equipment.ipynb
@@ -0,0 +1,62 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 16: Machinery and equipment"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 16.2: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"p=500//ft\n",
+"p1=6//in\n",
+"t=500//cfm\n",
+"p2=7//psig\n",
+"P=p2+14.7//psia\n",
+"T=520*(P/14.7)^0.283//F\n",
+"f=0.048*p1^0.027/(t)^0.148//in\n",
+"//CALCULATIONS\n",
+"delP=20*10^-3*p*T*(t)^2/(38*10^3*P*p1^5)//psia\n",
+"//RESULTS\n",
+"printf('the pressure drop=% f psia',delP)"
+   ]
+   }
+],
+"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/Water_and_Wastewater_Engineering_by_G_M_Fair/2-Water_system.ipynb b/Water_and_Wastewater_Engineering_by_G_M_Fair/2-Water_system.ipynb
new file mode 100644
index 0000000..5393ee1
--- /dev/null
+++ b/Water_and_Wastewater_Engineering_by_G_M_Fair/2-Water_system.ipynb
@@ -0,0 +1,301 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 2: Water system"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.1: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"w=3000//sq ft\n",
+"w1=2000//sq ft\n",
+"w2=1000//sq ft\n",
+"r=15//in\n",
+"a=12//in\n",
+"h=7.5//in\n",
+"//CALCULATIONS\n",
+"G=w*(r/a)*h//gal\n",
+"g=w1*(r/a)*h//gal\n",
+"g1=w2*(r/a)*h//gal\n",
+"//RESULTS\n",
+"printf('the normally be stored to tide the supply over dry spells=% f gal',G)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.2: Example_2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"m=17.378//mg\n",
+"h=20//in/sq mile \n",
+"d=365//in\n",
+"s=0.75//percent\n",
+"a=100//sq miles\n",
+"p=750000//gpd per sq mile\n",
+"t=180//in\n",
+"c=150//gpcd\n",
+"n=64699 //gpd per sq mile\n",
+"//CALCULATIONS\n",
+"R=h*m//mg\n",
+"A=R/d//mgd\n",
+"S=s*a*t//billion gal\n",
+"Q=a*p/c//gpd\n",
+"P=a*n/c//people against\n",
+"//RESULTS\n",
+"printf('the surface water sheds and storage requirements=% f mg',R)\n",
+"printf('the well watered sections of north america=% f billion gal',S)\n",
+"printf('the average consumpition populations=% f gpd',Q)\n",
+"printf('the presence of proper storage=% f people against',P)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.3: Example_3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"w=20//ft\n",
+"r=3//ft a day\n",
+"g=500//ft\n",
+"g1=1000//ft\n",
+"h=7.5/1440//ft\n",
+"p=7.5/1000000//ft\n",
+"r1=2//ft a day\n",
+"//CALCULATIONS\n",
+"W1=w*g1*r*h//gpm\n",
+"W2=w*g1*r1*r*p//mgd\n",
+"//RESULTS\n",
+"printf('the ground water laterally =% f gpm',W1)\n",
+"printf('the water from both sides=% f mgd',W2)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.4: Example_4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"p1=10//mgd\n",
+"p2=6940//gpm\n",
+"w=67000//people\n",
+"d=2//min\n",
+"v=d*p2/d//gal\n",
+"v1=98.2//cu ft each\n",
+"q=30//min\n",
+"q1=q*p2/d//gal\n",
+"q2=13900//cu ft\n",
+"a=1390//sq ft\n",
+"s=2//hr\n",
+"s1=120*p2/d//gal\n",
+"s2=55700//cu ft \n",
+"s3=s2/8//sq ft\n",
+"r=3//gpm/sq ft\n",
+"r1=6//rapid\n",
+"//CALCULATIONS\n",
+"D=sqrt(v1*4/%pi)//ft\n",
+"S=p2/s3//gpm/sq ft\n",
+"A=p2/(r1*r)//sq ft\n",
+"//RESULTS\n",
+"printf('the capacity of the components of a rapid sand filtration plant=% f sq ft',A)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.5: Example_5.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"r=10000//ft\n",
+"l=400000//people\n",
+"q=1000000//mgd\n",
+"w=100//gpcd\n",
+"w1=150//gpcd\n",
+"m=50//percent\n",
+"g=1.5//ft\n",
+"h1=2.32//cfs\n",
+"h2=139//cfs\n",
+"d=12//ft\n",
+"c=100//ft\n",
+"l=10.8//ft\n",
+"l2=0.85//ft\n",
+"l1=1000//ft\n",
+"//CALCULATIONS\n",
+"a=r*w/q//mgd\n",
+"b=l*w1/q//mgd\n",
+"a1=a*g//mgd\n",
+"b1=b*g//mgd\n",
+"D=d*sqrt(h1/%pi)//in\n",
+"D1=d*sqrt(h2/%pi)//in\n",
+"L=l/l1//ft\n",
+"L1=l2/l1//ft\n",
+"//RESULTS\n",
+"printf('the higher loss of head in small conduits at equal velocity=% f ft',L1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.6: Example_6.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"a=12//in\n",
+"b=24//in\n",
+"r=500//gpm\n",
+"d=200//gpcd\n",
+"d1=150//gpcd\n",
+"p1=113//sq in\n",
+"p2=425//sq in\n",
+"v1=3//fps\n",
+"v2=2.35//cfs\n",
+"v3=9.42//cfs\n",
+"h=646000//gpd\n",
+"w=720000//gpd\n",
+"//CALCULATIONS\n",
+"D1=v2*h//gpd\n",
+"D2=v3*h//gpd\n",
+"W1=D1-w//gpd\n",
+"W2=D2-w//gpd\n",
+"R1=D1/d//people\n",
+"R2=D2/d//people\n",
+"S=W1/d1//people\n",
+"S1=W2/d1//people\n",
+"//RESULTS\n",
+"printf('the absence of fire service for a maximum draft=% f gpd',D2)\n",
+"printf('The residential fire flow requirements=% f gpd',W2)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.7: Example_7.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"w=100000//ft\n",
+"c=250//per capita\n",
+"p1=0.3//percent\n",
+"p2=0.1//percent\n",
+"p3=0.60//percent\n",
+"w1=15//mgd\n",
+"//CALCULATIONS\n",
+"T=c*w//$\n",
+"W=p1*T//$\n",
+"W1=p2*T//$\n",
+"W2=p3*T//$\n",
+"W3=((w1)^2/3)*T//$\n",
+"//RESULTS\n",
+"printf('the replacement cost of the water of a city=% f $',W3)"
+   ]
+   }
+],
+"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/Water_and_Wastewater_Engineering_by_G_M_Fair/3-Wastewater_systems.ipynb b/Water_and_Wastewater_Engineering_by_G_M_Fair/3-Wastewater_systems.ipynb
new file mode 100644
index 0000000..82f2568
--- /dev/null
+++ b/Water_and_Wastewater_Engineering_by_G_M_Fair/3-Wastewater_systems.ipynb
@@ -0,0 +1,273 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 3: Wastewater systems"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.1: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"v=2.5//fps\n",
+"c=0.013//gpd\n",
+"p=300//gpd\n",
+"d=50//per care\n",
+"m=5.20//ft\n",
+"a=1000//ft\n",
+"//CALCULATIONS\n",
+"C=[(%pi*64)/(4*144)]*v*646000//gpd\n",
+"M=m/a//ft\n",
+"P=C/p\n",
+"A=P/d//acres\n",
+"//RESULTS\n",
+"printf('the acres will it drain if the population density=% f acres',A)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.2: Example_2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"a=37.4//acres\n",
+"r=2//in\n",
+"p=30//min\n",
+"v=3//fps\n",
+"r1=0.6//in\n",
+"h=1.0//cfs\n",
+"p1=50//percent\n",
+"q=646000//gpd\n",
+"//CALCULATIONS\n",
+"R=r*r1*a//cfs\n",
+"A=R/v//sq ft\n",
+"D=12*sqrt(4*A/%pi)//in\n",
+"P=r*r1*q/p1//gpcd\n",
+"//RESULTS\n",
+"printf('the per capita storm runoff for a population density=% f gpcd',P)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.3: Example_3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"w=1.0//cfs\n",
+"w1=3.0//cfs\n",
+"w2=45.0//cfs\n",
+"v=3.0//fps\n",
+"h=144//ft\n",
+"D=12*sqrt(4*w/(%pi*w1))//in\n",
+"d1=1.95//cfs\n",
+"D1=12*sqrt(4*d1)/(%pi*v)//in\n",
+"d2=41.6//cfs\n",
+"D2=12*sqrt(4*d2)/(%pi*w1)//ins\n",
+"//CALCULATIONS\n",
+"C=%pi*(D)^2*3/(4*h)//cfs\n",
+"C1=%pi*(1/4)*3//cfs\n",
+"V=(d2*4)/(%pi*4^2)//fps\n",
+"//RESULTS\n",
+"printf('The minimum dry-weather flow =% f cfs',C)\n",
+"printf('The maximum dry-weather flow in excess actual capacity=% f cfs',C1)\n",
+"printf('the storm flow in axcess of maximum dry-weather flow=% f fps',V)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.4: Example_4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"t=0.8//mgd\n",
+"d=8000//people\n",
+"a=2//hr\n",
+"v=800000//ft\n",
+"h=10//ft\n",
+"a1=4//in\n",
+"d1=1//sq ft per capita\n",
+"a3=3//mgad\n",
+"//CALCULATIONS\n",
+"V=v*(a/24)/a//gal\n",
+"S=V/h//sq ft\n",
+"S1=(v/h)/S//gpd per sq ft\n",
+"V1=a*d/h//cu ft\n",
+"D=d/S//ft\n",
+"E=d1*d/a1//sq ft\n",
+"A=t/a3//acre\n",
+"//RESULTS\n",
+"printf('the capacity of the components of a small trickling-filter =% f acre',A)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.5: Example_5.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"w=2000//sq miles\n",
+"r=0.1//cfs\n",
+"d=80000//ft\n",
+"p=100//gpd\n",
+"p1=80//ft\n",
+"p2=340//percent\n",
+"h=646000//ft\n",
+"//CALCULATIONS\n",
+"L=r*w//cfs\n",
+"R=6*p1/1.4//cfs\n",
+"P=p1*(p2-L)/p2//percent\n",
+"D=(d*p)//gpd\n",
+"D1=(L*h)//gpd\n",
+"//RESULTS\n",
+"printf('the percent of removal of pollutional load needed=% f percent',P)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.6: Example_6.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"p=10000//people\n",
+"p1=4//ft\n",
+"w=10//in\n",
+"s=80//gpcd\n",
+"h=43560//ft\n",
+"p2=20//ft\n",
+"//CALCULATIONS\n",
+"R=[(w/12)*7.5*h]/365//gpad\n",
+"A=p*s/R//acres\n",
+"A1=1.7//sq miles\n",
+"P=p/500//acres\n",
+"D=p2*h*4*7.48/(p*s)//days\n",
+"//RESULTS\n",
+"printf('the detention period in ponds =% f days',D)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.7: Example_7.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"p=100000//people\n",
+"a=75//$\n",
+"a2=47//in\n",
+"b=10//in\n",
+"//CALCULATIONS\n",
+"P=a*p//people\n",
+"S=((a)*(b^5))/(b)^1/4//$\n",
+"//RESULTS\n",
+"printf('the money is inversed in the sanitary sewerage system=% f $',S)"
+   ]
+   }
+],
+"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/Water_and_Wastewater_Engineering_by_G_M_Fair/4-Information_analysis.ipynb b/Water_and_Wastewater_Engineering_by_G_M_Fair/4-Information_analysis.ipynb
new file mode 100644
index 0000000..e139167
--- /dev/null
+++ b/Water_and_Wastewater_Engineering_by_G_M_Fair/4-Information_analysis.ipynb
@@ -0,0 +1,202 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 4: Information analysis"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 4.1: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"y=19.5//in\n",
+"x=396.8//in\n",
+"n=6//in\n",
+"y1=2.20//in\n",
+"x1=51.14//in\n",
+"p=5.64//in\n",
+"//CALCULATIONS\n",
+"Beta=(x-n*(y)*(y1))/(x1-n*(y1)^2)\n",
+"X=p+Beta//minimum\n",
+"//RESULTS\n",
+"printf('the method of leate squares =% f minimum',X)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 4.3: Example_2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"a=12//in\n",
+"h=121//in\n",
+"p=11//in\n",
+"s=220//in\n",
+"//CALCULATIONS\n",
+"B={a/[p*(h-1)]}*s//per unit\n",
+"//RESULTS\n",
+"printf('the interval of time a noted before=% f per unit',B)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 4.5: Example_3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"a=4404//ft\n",
+"q=9//ft\n",
+"mu=12//ft\n",
+"//CALCULATIONS\n",
+"F=sqrt(a/q)//ft\n",
+"CF=F/mu*100//percent\n",
+"//RESULTS\n",
+"printf('the coefficient of fluctuation is =% f percent',CF)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 4.7: Example_4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"h2=5//in\n",
+"x=3.72//in\n",
+"x1=1.28//in\n",
+"//CALCULATIONS\n",
+"H=h2*x1/x//in\n",
+"//RESULTS\n",
+"printf('the either side of the center of the scale=% f in',H)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 4.8: Example_5.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"p=80//in\n",
+"q=20//in\n",
+"//CALCULATIONS\n",
+"K=p+q//ft\n",
+"//RESULTS\n",
+"printf('the moments of the arithmetically normal frequency curve=% f ft',K)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 4.9: Example_6.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"g=3.2541//in\n",
+"g1=3.46//in\n",
+"m=0.5390//ft\n",
+"h=2/99//ft\n",
+"p=1.52//ft\n",
+"//CALCULATIONS\n",
+"L=sqrt(g*h)//in\n",
+"mu=g1*p//in\n",
+"M=g1/p//percent\n",
+"//RESULTS\n",
+"printf('the points necessary to plot the straigt line of fit on log probability=% f percent',M)"
+   ]
+   }
+],
+"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/Water_and_Wastewater_Engineering_by_G_M_Fair/5-Water_and_wastewater_volume.ipynb b/Water_and_Wastewater_Engineering_by_G_M_Fair/5-Water_and_wastewater_volume.ipynb
new file mode 100644
index 0000000..e11b7c0
--- /dev/null
+++ b/Water_and_Wastewater_Engineering_by_G_M_Fair/5-Water_and_wastewater_volume.ipynb
@@ -0,0 +1,249 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 5: Water and wastewater volume"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.1: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"t1=5.25//yr\n",
+"t2=10.00//yr\n",
+"yi=171000//in\n",
+"ye=111000//in\n",
+"yt=5.23300//in\n",
+"yl=5.04532//in\n",
+"yn=31500//in\n",
+"ym=0.09853//in\n",
+"tm=9.25//yr\n",
+"tn=10.00//yr\n",
+"//CALCULATIONS\n",
+"T=t1/t2//yr\n",
+"T1=tm/tn//yr\n",
+"Y=yi-ye//in\n",
+"Yt=yt-yl//in\n",
+"//RESULTS\n",
+"printf('the fifth intercensal year =% f yr',T)\n",
+"printf('the ninth postcensal year =% f yr',T1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.2: Example_2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"y0=30000//in\n",
+"y1=172000//in\n",
+"y2=292000//in\n",
+"a=172//ft\n",
+"p=30//ft\n",
+"y=292//ft\n",
+"q=322000//ft\n",
+"g=313//ft\n",
+"n=0.05//ft\n",
+"d=-2.442//ft\n",
+"//CALCULATIONS\n",
+"L=[2*p*a*y2-(a)^2*q]/[p*y-(a)^2]//moreover\n",
+"m=(g-p)/p//ft\n",
+"N=n*d//in\n",
+"Y=g/[1+m*(N)]//in\n",
+"//RESULTS\n",
+"printf('the saturation populations=% f moreover',L)\n",
+"printf('the coefficients=% f in',N)\n",
+"printf('the equation of a logistic curve=% f in',Y)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.4: Example_3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"p=100000//in\n",
+"d=150//in\n",
+"h=1000000//in\n",
+"a1=2.0//draft\n",
+"a2=3.0//draft\n",
+"a3=1.6//draft\n",
+"m=1.5//in\n",
+"q=2.5//in\n",
+"v=1020//in\n",
+"w=100//in\n",
+"t=0.01//in\n",
+"v1=13.2//mgd\n",
+"//CALCULATIONS\n",
+"A=d*p/h//mgd\n",
+"M=m*A//mgd\n",
+"M1=q*A//mgd\n",
+"V=v*sqrt(w)*(1-t*sqrt(w))//gpm\n",
+"D=M+v1//mgd\n",
+"L=a1*A//mgd\n",
+"L1=(4/3)*M//max\n",
+"H=a2*A//mgd\n",
+"H1=(4/3)*M1//max\n",
+"F=a3*A//mgd\n",
+"F1=(4/3)*M//max\n",
+"//RESULTS\n",
+"printf('the resulting capacities of the four system =% f max',F1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.6: Example_4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"r=48//in\n",
+"A=450//gpd/acre\n",
+"B=8000//gpd/mile\n",
+"S=5280/350//manholes/mile\n",
+"//CALCULATIONS\n",
+"C=(B-S*100)/12//gpd/mile\n",
+"//RESULTS\n",
+"printf('the ground a quarter of it eventually =% f gpd/mile',C)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.7: Example_5.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"p1=20//ft\n",
+"p2=30//ft\n",
+"w=5//person\n",
+"s=17800//in\n",
+"h=1200//in\n",
+"q=100//in\n",
+"i=1//in\n",
+"//CALCULATIONS\n",
+"S=p1*p2*i*s/(h*w)//gpcd\n",
+"P=(q*p1*p2/S)//percent\n",
+"//RESULTS\n",
+"printf('the degree of separation of stormwater=% f percent',P)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.8: Example_6.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"s=105//gpcd\n",
+"m=315//gpcd\n",
+"m1=35//gpcd\n",
+"Q1=360//gpcd\n",
+"Q2=30//gpcd\n",
+"p1=20//pecent\n",
+"p2=15//persons/acer\n",
+"D=21//persons/acre\n",
+"I=2000//gpd/acre\n",
+"//CALCULATIONS\n",
+"A=D*(s+Q2)+I//gpd/acre\n",
+"R=D*(m+Q2)+I//gpd/acre\n",
+"L=D*(m1+Q2)+I//gpd/acre\n",
+"//RESULTS\n",
+"printf('the average peak and low rates of flow =% f gpd/acre',L)"
+   ]
+   }
+],
+"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/Water_and_Wastewater_Engineering_by_G_M_Fair/6-Elements_of_hydrology.ipynb b/Water_and_Wastewater_Engineering_by_G_M_Fair/6-Elements_of_hydrology.ipynb
new file mode 100644
index 0000000..609375c
--- /dev/null
+++ b/Water_and_Wastewater_Engineering_by_G_M_Fair/6-Elements_of_hydrology.ipynb
@@ -0,0 +1,137 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 6: Elements of hydrology"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.1: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"H=1360//ft\n",
+"t=60//f\n",
+"a=(10^3)*5.5*(10^-3)//f\n",
+"q=(1.36*10^3)*5.5*(10^-3)//f\n",
+"s=(4-1.36)*(10^3)*(3.2*10^-3)//f\n",
+"//CALCULATIONS\n",
+"T=t-q-s//F\n",
+"T1=T+3*a//F\n",
+"//RESULTS\n",
+"printf('the temperature at the mountain top=% f F',T)\n",
+"printf('the temperature on the plain beyond the mountain=% f F',T1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.2: Example_2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"t=60//f\n",
+"v=0.52//in\n",
+"t1=80//F\n",
+"p=40//percent\n",
+"v1=1.03*0.40//in\n",
+"w=8//mph\n",
+"pa=29.0//in\n",
+"p1=0.497//ft\n",
+"q=1.32*10^-2//ft\n",
+"r=0.268//ft\n",
+"//CALCULATIONS\n",
+"E=p1*(1-q*pa)*(1+r*w)*(v-v1)//in\n",
+"//RESULTS\n",
+"printf('the evaporation for the a day during=% f in',E)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.3: Example_3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"t=47//f\n",
+"q=8000//ft\n",
+"a=100//ft\n",
+"d=0.10//in\n",
+"d1=7//degree days\n",
+"s1=14000//ft\n",
+"s2=7000//ft\n",
+"s=1000//ft\n",
+"g=32//ft\n",
+"h=17.37//ft\n",
+"h1=1.547//ft\n",
+"//CALCULATIONS\n",
+"T=q+s*(t-g)/3//ft\n",
+"T1=t-3*1//F\n",
+"T2=(T1+g)/2//F\n",
+"T3=d1*d*a//sq mile in\n",
+"M=h*T3//mgd\n",
+"M1=M*h1//cfs\n",
+"//RESULTS\n",
+"printf('the upper boundary of the melting zone and temperature at the snow line=% f F',T1)\n",
+"printf('The average temperature of =% f cfs',M1)"
+   ]
+   }
+],
+"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/Water_and_Wastewater_Engineering_by_G_M_Fair/7-Rainfall_and_runoff.ipynb b/Water_and_Wastewater_Engineering_by_G_M_Fair/7-Rainfall_and_runoff.ipynb
new file mode 100644
index 0000000..9a10bfe
--- /dev/null
+++ b/Water_and_Wastewater_Engineering_by_G_M_Fair/7-Rainfall_and_runoff.ipynb
@@ -0,0 +1,129 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 7: Rainfall and runoff"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.5: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"n=20//ft\n",
+"s=sqrt(12676/19)//ft\n",
+"c=45.5//ft\n",
+"q=551400//ft\n",
+"q1=12700//ft\n",
+"h=8.5//ft\n",
+"w=s/c//ft\n",
+"//CALCULATIONS\n",
+"D=q/(2*s*q1)//cfs\n",
+"D1=D*(1+h/n)//cfs\n",
+"//RESULTS\n",
+"printf('the record runoff of a stream draining=% f cfs',D1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.6: Example_2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"i=16/(62)^0.66//in hr\n",
+"q=(16*10^0.31)/(62)^0.66//in hr\n",
+"c=1.0//max\n",
+"C1=c*(0.01)^0.31//in\n",
+"C2=c*(0.1)^0.31//in\n",
+"x1=640//cfs\n",
+"//CALCULATIONS\n",
+"Y1=C1*i*c*x1//cfs\n",
+"Y2=C2*q*c*x1//cfs\n",
+"//RESULTS\n",
+"printf('the time of concentration=% f cfs',Y2)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.8: Example_3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"d=163*48.5//cfs\n",
+"a=48.5//ft\n",
+"q=100//cfs\n",
+"Q=45.5*a//cfs\n",
+"c=0.57//cfs\n",
+"v=1.8//cfs\n",
+"p=0.45//ft\n",
+"//CALCULATIONS\n",
+"P=d/(q*sqrt(a))//percent\n",
+"C=Q/(a^0.8*(1+2*a^-0.3))//cfs\n",
+"d1=2.6//cfs\n",
+"T=(1-p*c+v*c*2)//cfs\n",
+"//RESULTS\n",
+"printf('the meyers rating =% f percent',P)\n",
+"printf('the magnitude of the maximum peak flood =% f cfs',T)"
+   ]
+   }
+],
+"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/Water_and_Wastewater_Engineering_by_G_M_Fair/8-Storage_and_runoff_control.ipynb b/Water_and_Wastewater_Engineering_by_G_M_Fair/8-Storage_and_runoff_control.ipynb
new file mode 100644
index 0000000..5996112
--- /dev/null
+++ b/Water_and_Wastewater_Engineering_by_G_M_Fair/8-Storage_and_runoff_control.ipynb
@@ -0,0 +1,167 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 8: Storage and runoff control"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.2: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"a=0.75//ft\n",
+"p=123//mg\n",
+"v=100//ft\n",
+"s=33//mg\n",
+"s1=67//mg\n",
+"d=26.6//mgd\n",
+"d1=0.0477//mgd\n",
+"q=0.750//gpd/sq mile\n",
+"d2=365//days\n",
+"//CALCULATIONS\n",
+"S=p/a//mg per sq mile\n",
+"Cv=v*s/s1//percent\n",
+"M=d*d1//mgd per sq mile\n",
+"D=v*q/M//MAF\n",
+"D1=(v*p)/(M*d2)//MAF\n",
+"R=p/q//days\n",
+"//RESULTS\n",
+"printf('the use monthly averages rather then daily stream flow=% f days',R)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.3: Example_2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"d=750000//gpd per sq mile\n",
+"v=0.22//ft\n",
+"a=1.27//ft\n",
+"q=0.30//ft\n",
+"d1=365//days\n",
+"p=0.25//ft\n",
+"//CALCULATIONS\n",
+"Q=q*a*d1//mg/sq mile\n",
+"H=p*a*d1//mg/sq mile\n",
+"//RESULTS\n",
+"printf('the results obtained by normal analytical procedures and by Hazen s=% f mg/sq mile',H)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.4: Example_3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"d=30.0//mgd\n",
+"a=40.0//sq miles\n",
+"a1=1500//acres\n",
+"r1=47.0//in\n",
+"r2=27.0//in\n",
+"q=0.9//in\n",
+"k=640//in\n",
+"h=0.052//gpd/sq mile\n",
+"//CALCULATIONS\n",
+"Q=r2-(r2+a-r1)*q*a1/(k*a)//in\n",
+"D=d+a*h//mgd\n",
+"A=a-(q*a1/k)*[1-(r1-a)/(r2)]//sq miles\n",
+"R=r2+a-r1//in\n",
+"S=R*q//in\n",
+"//RESULTS\n",
+"printf('the revised mean annual runoff=% f in',Q)\n",
+"printf('the equivalent mean draft=% f mgd',D)\n",
+"printf('the equivalent land area=% f sq miles',A)\n",
+"printf('the adjusted flowline=% f in',S)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.6: Example_4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"p=100//ft\n",
+"q=27000//acre-ft\n",
+"p1=10//ft\n",
+"s=8250//acre-ft\n",
+"//CALCULATIONS\n",
+"R=p*s/q//percent\n",
+"//RESULTS\n",
+"printf('the ratio of peak inflow from fuller values=% f percent',R)"
+   ]
+   }
+],
+"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/Water_and_Wastewater_Engineering_by_G_M_Fair/9-Groundwater_flow.ipynb b/Water_and_Wastewater_Engineering_by_G_M_Fair/9-Groundwater_flow.ipynb
new file mode 100644
index 0000000..a4316bc
--- /dev/null
+++ b/Water_and_Wastewater_Engineering_by_G_M_Fair/9-Groundwater_flow.ipynb
@@ -0,0 +1,239 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 9: Groundwater flow"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.1: Example_1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"t=10//C\n",
+"s=74.2//days\n",
+"c=0.01//mm\n",
+"d=245//mm\n",
+"//CALCULATIONS\n",
+"h=s/(d*c)//cm\n",
+"//RESULTS\n",
+"printf('the high will water at a temperature =% f cm',h)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.2: Example_2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"p1=1000//ft\n",
+"p2=50//ft\n",
+"g=20//ft/mile\n",
+"v=5280//ft\n",
+"q=7.5*10^-6//ft\n",
+"t=60//F\n",
+"k=2835//ft/days\n",
+"p=7.5//ft\n",
+"//CALCULATIONS\n",
+"S=g/v//ft\n",
+"W=k*(g/v)//ft/day\n",
+"Q=W*p1*p2*q//mgd\n",
+"P=k*p//ft\n",
+"P1=P*p2//mgd\n",
+"//RESULTS\n",
+"printf('the velocity of flow =% f mgd',Q)\n",
+"printf('the standard coefficient pf permeability=% f mgd',P1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.3: Example_3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"p=40//ft\n",
+"d=56//ft\n",
+"d1=140//ft\n",
+"p1=30//ft\n",
+"w=3.28*10^-4//fps\n",
+"//CALCULATIONS\n",
+"Q=w*(p/d1)*2*d*p//cfs\n",
+"q=Q/p//cfs\n",
+"K=w*(p/d1)//fps\n",
+"x0=q/(2*%pi*K)//ft\n",
+"Z=2*%pi*x0//ft\n",
+"//RESULTS\n",
+"printf('the yield of the well if the coefficient of permeability=% f ft',x0)\n",
+"printf('the distance of the point of stagnation =% f ft',Z)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.4: Example_4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"p=5*10^6//ft\n",
+"Q=350//gpm\n",
+"x=225//ft\n",
+"u=10^-2//ft\n",
+"g=1.87//ft\n",
+"p2=7*10^2//ft\n",
+"p3=10.9//ft\n",
+"w=4.0//ft\n",
+"t=114.6//ft\n",
+"d=10//ft\n",
+"p1=5//ft\n",
+"w1=3.2*10^4//ft\n",
+"W=21.75//ft\n",
+"//CALCULATIONS\n",
+"T=t*Q*4/p1//gpd/ft\n",
+"S=u*(w1)/[g*(p)]//ft\n",
+"U=g*[(S)/(T)]*(x^2/d)//ft\n",
+"P=t*(p2)*p3/(T)//ft\n",
+"U1=g*[(S)/(T)]*(1/d)//ft\n",
+"P1=t*(p2)*W/(T)//ft\n",
+"//RESULTS\n",
+"printf('the type curve as if a transparency of the observed data had moved into place over the type=% f ft',P1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.5: Example_5.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"Q=350//gpm\n",
+"x=225//ft\n",
+"t=1//min\n",
+"p=1.6//ft\n",
+"t2=10//min\n",
+"p2=4.5//ft\n",
+"p3=700//gpm\n",
+"T=3.2*10^4//gpd/ft\n",
+"t0=0.3//min\n",
+"u=1.15*10^-5\n",
+"//CALCULATIONS\n",
+"S=t0*(T)*t0/[(x)^2*1440]//ft\n",
+"P=[(114.6*p3)/(T)]*(-0.5772*2.3*log(u))//ft\n",
+"//RESULTS\n",
+"printf('A straight line being drawn through the ppints for the higher=% f ft',P)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.6: Example_6.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"h=4.8//ft\n",
+"m=13.4//ft\n",
+"k=10^-1//cm/sec\n",
+"k1=3.28*10^-3//fps\n",
+"n=7//ft\n",
+"n1=11//ft\n",
+"q=1.0*10^-2\n",
+"//CALCULATIONS\n",
+"Q=k1*h*n/n1//cfs/ft\n",
+"Q1=2*q*10^2//cfs\n",
+"//RESULTS\n",
+"printf('A satisfactory orthogonal system the flow of into the collector =% f cfs',Q1)"
+   ]
+   }
+],
+"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
+}