Merge pull request #1 from prashantsinalkar/masterHEADmaster

Initial commit

15 files changed, 3086 insertions, 0 deletions

diff --git a/Principles_Of_Geotechnical_Engineering_by_B_M_Das/10-stress_in_a_soil_mass.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/10-stress_in_a_soil_mass.ipynb
new file mode 100644
index 0000000..23b49a9
--- /dev/null
+++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/10-stress_in_a_soil_mass.ipynb
@@ -0,0 +1,96 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 10: stress in a soil mass"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 10.1: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"sx= 2000 //lb/ft^3\n",
+"sy= 2500 //lb/ft^3\n",
+"T= 800 //lb/ft^3\n",
+"t= 0.348//radians\n",
+"//calculations\n",
+"s1= (sx+sy)/2+sqrt(((sy-sx)/2)^2+T^2)\n",
+"s2= (sx+sy)/2-sqrt(((sy-sx)/2)^2+T^2)\n",
+"sn= (sx+sy)/2+(sy-sx)*cos(2*t)/2-T*sin(2*t)\n",
+"Tn= (sy-sx)*sin(2*t)/2+T*cos(2*t)\n",
+"//results\n",
+"printf ('principle stress s1 = % 2f lb/ft^3 ',s1)\n",
+"printf ('principle stress s2 = % 2f lb/ft^3 ',s2)\n",
+"printf ('normal stress = % 2f lb/ft^3 ',sn)\n",
+"printf ('shear stress = % 2f lb/ft^3 ',Tn)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 10.3: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"x= 3 //m\n",
+"y= 4 //m\n",
+"P= 5 //kN\n",
+"z= 2 //m\n",
+"//calculations\n",
+"r= sqrt(x^2+y^2)\n",
+"k= r/z\n",
+"I= 3/(2*%pi*((r/z)^2+1)^2.5)\n",
+"s= P*I/z^2\n",
+"//results\n",
+"printf ('verticle stress increase at 2m = % 4f kN/m^3 ',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/Principles_Of_Geotechnical_Engineering_by_B_M_Das/11-Compressibility_of_Soil.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/11-Compressibility_of_Soil.ipynb
new file mode 100644
index 0000000..3bab158
--- /dev/null
+++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/11-Compressibility_of_Soil.ipynb
@@ -0,0 +1,439 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 11: Compressibility of Soil"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.10: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"To=3000 // lb/ft^2\n",
+"eo=1.1\n",
+"e1=0.9\n",
+"e=eo-e1\n",
+"ea=(eo+e1)/2\n",
+"T1=3000 // lb/ft^2\n",
+"T=1 // in\n",
+"t = 2 // min\n",
+"m=(e/T1)/(1+ea)\n",
+"U=50\n",
+"Tv=0.197\n",
+"Gw=62.4 //lb/ft^3\n",
+"Cv=Tv*(T/(2*12)^2)/t\n",
+"k=Cv*m*Gw *10^7\n",
+"printf('a)k = %f x10^-7 ft/min\n',k)\n",
+"\n",
+"\n",
+"U=60\n",
+"Tv=0.286\n",
+"H=6\n",
+"t60=Tv*H^2/(Cv*60*24)\n",
+"printf(' b)t60 =%f days',t60)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.11: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"t50=19\n",
+"Hd=2.24/2\n",
+"Cv=0.197*Hd^2/t50\n",
+"printf('Cv = %f cm^2/min',Cv)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.12: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"LL=40\n",
+"Cc=0.009*(LL-10)\n",
+"H=10*12\n",
+"eo=1.0\n",
+"Gss=120\n",
+"Gsc=110\n",
+"Gd=100\n",
+"To=10*Gd +10*(Gss-62.4)+10*(Gsc-62.4)/2\n",
+"\n",
+"Tt=0.408\n",
+"Tm=0.232\n",
+"Tb=0.019\n",
+"Tav= (Tt+4*Tm+Tb)/6\n",
+"Sc=Cc*H*log10((To+Tav*1000)/To)/(1+eo)\n",
+"printf('Sc = %f in',Sc)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.13: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"H = 6\n",
+"Cc = 0.28\n",
+"eo = 0.9\n",
+"Cv = 0.36\n",
+"To=210\n",
+"Tp=115\n",
+"Sc= Cc*H*log10((To+Tp)/To)/(1+eo)\n",
+"t2=9\n",
+"Hd=3\n",
+"Tv=Cv*t2/Hd^2\n",
+"U=0.67\n",
+"Tf=0.677*Tp\n",
+"printf('Tf = %f kN/m^2',Tf)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.1: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"Tz=150\n",
+"b=1\n",
+"l=2\n",
+"z=5*b\n",
+"Es= (10000*2 + 8000*1 +12000*2)/5\n",
+"a=4\n",
+"H=z\n",
+"m=l/b\n",
+"n=2*H/b\n",
+"F1=0.641 //from tables 11.1 and 11.2\n",
+"F2=0.031\n",
+"u=0.3\n",
+"Is= F1 + ((2-u)/(1-u))*F2\n",
+"If=0.71 //from  table 11.3\n",
+"Sef= Tz *a*b/l *(1-u^2)*Is*If/Es\n",
+"Ser=0.93*Sef\n",
+"printf('The elastic settlement at the centre of foundation = %f m',Ser)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.2: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"// one value of e is done \n",
+"Gs=2.75\n",
+"A=30.68\n",
+"Ms=128\n",
+"p=1\n",
+"Hs=Ms/(A*Gs*p)\n",
+"H=2.540\n",
+"Hv=H-Hs\n",
+"e=Hv/Hs\n",
+"printf('the value of e for give values = %f',e)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.3: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"e11=0.9\n",
+"e21=0.8\n",
+"T2=4\n",
+"T1=2\n",
+"Cc= (e11-e21)/log10(T2/T1) // from loading branch\n",
+"e1=0.67\n",
+"e2=0.655\n",
+"Cs=(e1-e2)/log10(T2/T1)\n",
+"k=Cs/Cc\n",
+"T3=12\n",
+"e3=e11-Cc*log10(T3/T1)\n",
+"printf('Compression index Cc= %f\n',Cc)\n",
+"printf(' Cs/Cc = %f\n',k)\n",
+"printf(' e3 = %f',e3)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.4: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"Gd=14\n",
+"Gss=18\n",
+"Gsc=19\n",
+"Gw=9.81\n",
+"To= 2*Gd+4*(Gss-Gw)+2*(Gsc-Gw)\n",
+"LL=40\n",
+"Cc=0.009*(LL-10)\n",
+"H=4\n",
+"T=100\n",
+"e=0.8\n",
+"Sc= Cc*H*log10((To+T)/To)/(1+e)\n",
+"printf('a)Primary Consolidation Sc = %f m\n',Sc)\n",
+"\n",
+"Tc=190\n",
+"Cs=Cc/6\n",
+"Sc= Cs*H*log10((To+T)/To)/(1+e)\n",
+"printf(' b)Primary Consolidation Sc = %f m\n',Sc)\n",
+"\n",
+"Tc=170\n",
+"Sc= Cc*H*log10((To+T)/Tc)/(1+e)+ Cs*H*log10(Tc/To)/(1+e)\n",
+"printf(' c)Primary Consolidation Sc = %f m\n',Sc)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.5: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"Gs=18\n",
+"Gw=9.81\n",
+"H=10\n",
+"eo=1.1\n",
+"To=5*(Gs-Gw)\n",
+"T1=48\n",
+"T=To+T1\n",
+"e1=1.045 // void ratio corresponding to T \n",
+"e=eo-e1\n",
+"Sc=H*e/(1+eo)\n",
+"printf('The settlement in the field Sc = %f m',Sc)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.6: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"T=8.5\n",
+"eo=0.8\n",
+"Cc=0.28\n",
+"To=2650\n",
+"T1=970\n",
+"C1=0.02\n",
+"t2=5\n",
+"t1=1.5\n",
+"H=8.5*12\n",
+"epr=Cc*log10((To+T1)/To)\n",
+"ep=eo-epr\n",
+"C2=C1/(1+ep)\n",
+"Sc=epr*H/(1+eo)\n",
+"Ss=C2*H*log10(t2/t1)\n",
+"TS=Sc+Ss\n",
+"printf('Total consolidation settlement of the clay = %f in',TS)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.7: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//T50 = Cvtlab /H^2 lab = Cvtfield?H^2 fiels\n",
+"tl=140\n",
+"Hf=3\n",
+"Hd=0.025/2\n",
+"tf=tl*Hf^2/Hd^2\n",
+"k=tf/(3600*24)\n",
+"printf('t field = %f days',k)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.8: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//Tv is directly proportional to U^2\n",
+"t1=93.333\n",
+"U2=30\n",
+"U1=50\n",
+"t2=t1*U2^2/U1^2\n",
+"printf('t2 = %f days',t2)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.9: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"t90=75*24*60*60 // time in sec\n",
+"T90=0.848\n",
+"Hd=1.5*100 //in cm\n",
+"Cv=T90*Hd^2/t90\n",
+"printf('Cv = %f cm^2/sec',Cv)"
+   ]
+   }
+],
+"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/Principles_Of_Geotechnical_Engineering_by_B_M_Das/12-Shear_Strength_of_Soil.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/12-Shear_Strength_of_Soil.ipynb
new file mode 100644
index 0000000..3921e14
--- /dev/null
+++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/12-Shear_Strength_of_Soil.ipynb
@@ -0,0 +1,252 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 12: Shear Strength of Soil"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 12.1: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"D=50 // in mm\n",
+"A= %pi/4 *(D/1000)^2\n",
+"// solving for test 1 \n",
+"N=150\n",
+"Sp=157.5\n",
+"Sr=44.2\n",
+"Tf=Sp/A\n",
+"Tr=Sr/A\n",
+"// from graph\n",
+"k=tand(27)\n",
+"k1=tand(14.6)\n",
+"\n",
+"printf('Peak strength Tf = 40+ t*%f\n',k )\n",
+"printf(' Residual strength Tr = t*%f',k1)\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 12.2: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"T3=16 // lb/in^2\n",
+"Tf=25 // lb/in^2\n",
+"T1=T3+Tf\n",
+"a= asind((T1-T3)/(T1+T3)) // Mohr's circle\n",
+"printf('a)Angle of friction,a = %f\n',a)\n",
+"b= 45+ a/2\n",
+"printf(' b)Angle b that the failure plane makes with the major principal plane = %f',b)\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 12.3: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"T1=41\n",
+"T3=16\n",
+"a=58\n",
+"T=(T1+T3)/2 + (T1-T3)*cosd(2*a)/2\n",
+"tf=(T1-T3)*sind(2*a)/2\n",
+"printf('a)the normal stress T = %f lb/in^2',T)\n",
+"printf(' and the shear stress tf = %f lb/in^2',tf)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 12.4: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//For normally consolidated clay, c' = 0.\n",
+"a=30\n",
+"T3=10\n",
+"T1=T3*(tand(45+a/2))^2\n",
+"Tf=T1-T3\n",
+"printf('The deviator stress at failure = %f lb/in^2',Tf)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 12.5: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"T13=70\n",
+"T1f=130\n",
+"T11=T13+T1f\n",
+"\n",
+"T23=160\n",
+"T2f=223.5\n",
+"T21=T23+T2f\n",
+"\n",
+"a= 2*(atand(((T11-T21)/(T13-T23))^0.5)-45)\n",
+"c= (T11-T13*(tand(45+a/2))^2)/(2*tand(45+a/2))\n",
+"printf('the shear strength parameter c = %f kN/m^2',c)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 12.6: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"T3=12\n",
+"Tf=9.1\n",
+"T1=T3+Tf\n",
+"u=6.8\n",
+"a= asind((T1-T3)/(T1+T3))\n",
+"\n",
+"a1= asind((T1-T3)/(T1+T3-2*u))\n",
+"\n",
+"printf('a)Consolidated-undrained angle of shearing resistance = %f degrees\n',a)\n",
+"printf(' b)Drained friction angle = %f degrees',a1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 12.7: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"T3=12\n",
+"a=27.8\n",
+"T1=T3*(tand(45+a/2))^2\n",
+"Tf=T1-T3\n",
+"printf('the deviator stress at failure = %f lb/in^2',Tf)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 12.8: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"PI=28\n",
+"OCR=3.2\n",
+"To=160\n",
+"Kn=0.11+0.0037*PI\n",
+"Ko=OCR^0.8 * Kn\n",
+"Cu=Ko*To\n",
+"printf('the average undrained shear strength of the clay = %f kN/m^2',Cu)"
+   ]
+   }
+],
+"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/Principles_Of_Geotechnical_Engineering_by_B_M_Das/13-Lateral_Earth_Pressure.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/13-Lateral_Earth_Pressure.ipynb
new file mode 100644
index 0000000..9da380b
--- /dev/null
+++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/13-Lateral_Earth_Pressure.ipynb
@@ -0,0 +1,350 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 13: Lateral Earth Pressure"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 13.10: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"H=28\n",
+"C=210\n",
+"b=10\n",
+"G=118\n",
+"c=20\n",
+"kh=0.1\n",
+"Ka=tand(45-c/2)\n",
+"zo=2*C/(G*(Ka))\n",
+"n=zo/(H-zo)\n",
+"Nac=1.60\n",
+"Nav=0.375\n",
+"L=1.17\n",
+"Pae= G*(H-zo)^2*(L*Nav)-C*(H-zo)*Nac\n",
+"printf('The magnitude of the active force, Pae = %f lb/ft',Pae)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 13.1: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"OCR=2\n",
+"a=30\n",
+"Ko=(1-sind(a))*(OCR)^sind(a)\n",
+"//at z=0\n",
+"To1=0\n",
+"Th1=0\n",
+"u1=0\n",
+"//at z=10\n",
+"To2=10*100\n",
+"Th2=Ko*To2\n",
+"u2=0\n",
+"//at z=15\n",
+"To3= 10*100+5*(122.4-62.4)\n",
+"Th3=Ko*To3\n",
+"u3=5*62.4\n",
+"//Lateral force Po =Area 1 + Area 2+ Area3+ Area 4\n",
+"Po =(1/2)*10*707+5*707+(1/2)*5*212.1+(1/2)*5*312\n",
+"z=((3535)*(5+10/3)+3535*(5/2)+530.3*(5/3)+780*(5/3))/Po\n",
+"printf('z = %f ft',z)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 13.2: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//c=0\n",
+"a=36\n",
+"G=16\n",
+"Ka=(1-sind(a))/(1+sind(a))\n",
+"//at z=0 Tp=0\n",
+"z=6\n",
+"To=G*z\n",
+"Ta=Ka*To\n",
+"Pa=z*Ta/2\n",
+"\n",
+"printf('a)Rankine active force per unit length of the wall = %f kN/m',Pa)\n",
+"printf(' and the location of the resultant is z = 2m\n')\n",
+"\n",
+"\n",
+"p=36\n",
+"G=16\n",
+"Kp=(1+sind(a))/(1-sind(a))\n",
+"//at z=0 Tp=0\n",
+"z=6\n",
+"To=G*z\n",
+"Tp=Kp*To\n",
+"Pp=z*Tp/2\n",
+"\n",
+"printf(' b)Rankine passive force per unit length of the wall = %f kN/m',Pp)\n",
+"printf(' and the location of the resultant is z = 2m')\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 13.3: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"H=12\n",
+"a=20\n",
+"b=20\n",
+"G=115\n",
+"c=30\n",
+"Oa= asind(sind(a)/sind(c))-a+2*b\n",
+"Ka= (cosd(a-b)*sqrt(1+(sind(c))^2-2*sind(c)*cosd(Oa)))/((cosd(b))^2*(cosd(a)+sqrt((sind(c))^2-(sind(a))^2)))\n",
+"Pa=G*H^2*Ka/2\n",
+"B= atand((sind(c)*sind(Oa))/(1-(sind(c)*cosd(Oa))))\n",
+"printf('The active force Pa per unit length of the wall = %f lb/ft\n',Pa)\n",
+"printf( ' The resultant will act a distance of 12/3 = 4 ft above the bottom of the wall with B = %f degree',B)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 13.4: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"a=30\n",
+"Ka1=(1-sind(a))/(1+sind(a))\n",
+"a=35\n",
+"Ka2=(1-sind(a))/(1+sind(a))\n",
+"//at z=0 so T0=0\n",
+"//atz=3\n",
+"To=3*16\n",
+"Ta1=Ka1*To\n",
+"Ta2=Ka2*To\n",
+"\n",
+"// At z=6\n",
+"To=3*16+3*(18-9.81)\n",
+"Ta2=Ka2*To\n",
+"\n",
+"Pa =(1/2)*3*16+3*13.0+ (1/2)*3*36.1\n",
+"z= (24 *(3+3/3)+39.0*(3/2)+54.15*(3/3))/Pa\n",
+"printf('The force per unit length of the wall = %f kN/m\n',Pa)\n",
+"printf (' The location of the resultant = %f m ',z)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 13.5: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"Ka= (tand(1))^2\n",
+"G=16.5\n",
+"cu=10\n",
+"H=6\n",
+"//at z=0\n",
+"z=0\n",
+"Ta=G*z-2*cu\n",
+"//zt z=6\n",
+"z=6\n",
+"Ta=G*z-2*cu\n",
+"\n",
+"zo=2*cu/G\n",
+"// Before the tensile crack occurs\n",
+"Pa= G*H^2/2 - 2*cu*H\n",
+"printf('Pa before the tensile crack occurs = %f kN/m\n',Pa)\n",
+"//After the tensile crack occurs\n",
+"Pa=(H-zo)*Ta/2\n",
+"printf(' Pa after the tensile crack occurs = %f kN/m',Pa)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 13.6: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"H=15\n",
+"a=10\n",
+"G=118\n",
+"b=20\n",
+"C=250\n",
+"Zo=2*C*sqrt((1+sind(b))/(1-sind(b)))/G\n",
+"//at z=0 Ta=0\n",
+"//at z=15 \n",
+"z=15\n",
+"K=0.3\n",
+"Ta=G*z*K*cosd(a)\n",
+"Pa=(H -Zo)*Ta/2\n",
+"printf('The Rankine active force Pa on the retaining wall after the tensile crack occurs = %f lb/ft',Pa)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 13.7: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"c=30\n",
+"b=15\n",
+"a=10\n",
+"Ka=0.3872 // from table 13.8\n",
+"H=4\n",
+"G=15\n",
+"Pa=G*H^2*Ka/2\n",
+"printf('The active force per unit length Pa = %f kN/m\n',Pa)\n",
+"printf(' The resultant will act at a vertical distance equal to H/3 = 4/3 = 1.33 m above \n the bottom of the wall and will be inclined at an angle of 15to the back face of the wall.')"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 13.9: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"kh=0.2\n",
+"kv=0\n",
+"H=4\n",
+"a=0\n",
+"b=0\n",
+"c=15\n",
+"d=30\n",
+"G=15.5\n",
+"B= atand(kh/(1-kv))\n",
+"b1=b+B\n",
+"a1=a+B\n",
+"Ka=0.452\n",
+"Pa=G*H^2*Ka/2\n",
+"Pae=Pa*(1-kv)*((cosd(b1))^2/((cosd(b))^2*(cosd(B))^2))\n",
+"Ka=0.3014\n",
+"Pa=G*H^2*Ka/2\n",
+"P1=Pae-Pa\n",
+"z= ((Pa*H/3)+P1*0.6*H)/Pae\n",
+"printf('Pae = %f kN/m\n',Pae)\n",
+"printf(' Z = %f m',z)"
+   ]
+   }
+],
+"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/Principles_Of_Geotechnical_Engineering_by_B_M_Das/14-Lateral_Earth_Pressure_Curved_Failure_Surface.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/14-Lateral_Earth_Pressure_Curved_Failure_Surface.ipynb
new file mode 100644
index 0000000..73ef301
--- /dev/null
+++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/14-Lateral_Earth_Pressure_Curved_Failure_Surface.ipynb
@@ -0,0 +1,104 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 14: Lateral Earth Pressure Curved Failure Surface"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 14.1: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"G=15.7\n",
+"a=0\n",
+"b=15\n",
+"c=30\n",
+"H=3\n",
+"Kp=4.977 // from table 13.9\n",
+"Pp=Kp*G*H^2/2\n",
+"printf('a)the passive force = %f kN/m\n',Pp)\n",
+"// for part b\n",
+"Kp=4.53\n",
+"Pp=Kp*G*H^2/2\n",
+"printf(' b)the passive force = %f kN/m\n',Pp)\n",
+"// for part c\n",
+"Kp=4.13\n",
+"Pp=Kp*G*H^2/2\n",
+"printf(' c)the passive force = %f kN/m\n',Pp)\n",
+"//for part d\n",
+"Kp=4.56\n",
+"Pp=Kp*G*H^2/2\n",
+"printf(' d)the passive force = %f kN/m\n',Pp)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 14.2: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"G=16\n",
+"H=7\n",
+"c=30\n",
+"Ta=0.65*G*H*(tand(45-c/2))^2\n",
+"A=Ta*3*3/4\n",
+"B1=Ta*3-54.61\n",
+"C=Ta*4*4/4\n",
+"B2=Ta*4-97.08\n",
+"s=2\n",
+"As=A*s\n",
+"Bs=(B1+B2)*s\n",
+"Cs=C*s\n",
+"printf( 'The strut loads at level A = %f kN\n',As)\n",
+"printf( ' The strut loads at level B = %f kN\n',Bs)\n",
+"printf( ' The strut loads at level C = %f kN',Cs)"
+   ]
+   }
+],
+"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/Principles_Of_Geotechnical_Engineering_by_B_M_Das/15-Slope_Stability.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/15-Slope_Stability.ipynb
new file mode 100644
index 0000000..6968ca3
--- /dev/null
+++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/15-Slope_Stability.ipynb
@@ -0,0 +1,333 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 15: Slope Stability"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 15.11: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"C=20\n",
+"G=18.5\n",
+"H=21.62\n",
+"c=25\n",
+"r=0.25\n",
+"Fs=3.1*tand(c)\n",
+"printf('Fs = %f',Fs)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 15.1: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"Gs=17.8\n",
+"Gw=9.81\n",
+"C=10\n",
+"c=20\n",
+"b=15\n",
+"H=6\n",
+"G=Gs-Gw\n",
+"Fs= C/(Gs*H*cosd(b)*cosd(b)*tand(b))+G*tand(c)/(Gs*tand(b))\n",
+"printf('a)The factor of safety = %f \n',Fs)\n",
+"Fs=2\n",
+"H=2.247/(Fs-0.61)\n",
+"printf(' b)H= %f m',H)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 15.2: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"G=105\n",
+"c=15\n",
+"C=600\n",
+"b=45\n",
+"Fs=3\n",
+"Cd=C/Fs\n",
+"c1= atand(tand(c)/Fs)\n",
+"\n",
+"H= 4*Cd*(sind(b)*cosd(c1)/(1-cosd(b-c1)))/G\n",
+"printf('The depth of the cut slope = %f ft',H)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 15.3: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"Cu=40\n",
+"G=17.5\n",
+"b=60\n",
+"a=35\n",
+"c=72.5\n",
+"m=0.195\n",
+"Hc=Cu/(G*m)\n",
+"r=Hc/(2*sind(a)*sind(c/2))\n",
+"BC=Hc*((1/tand(a))-(1/tand(b)))\n",
+"printf('a)The maximum depth Hc = %f m\n',Hc)\n",
+"printf(' b)The radius, r = %f m\n',r)\n",
+"printf(' c)The distance BC.= %f m',BC)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 15.4: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"Gs=17.29\n",
+"d=9.15\n",
+"d1=6.1\n",
+"D=d/d1\n",
+"a=40\n",
+"m=0.175\n",
+"b=40\n",
+"H=6.1\n",
+"Cu=H*Gs*m\n",
+"printf('a)The undrained cohesion of the clay Cu = %f kN/m^2\n',Cu)\n",
+"printf(' b)The nature of the critical circle is midpointcircle\n')\n",
+"d=1.5\n",
+"b=40\n",
+"n=0.9\n",
+"D1=n*H\n",
+"printf(' c)Distance = %f m',D1)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 15.5: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"Fs=1\n",
+"b=56\n",
+"Kh=0.25\n",
+"M=3.66\n",
+"Cu=500\n",
+"G=100\n",
+"Hc=Cu*M/G\n",
+"printf('a)The maximum depth = %f ft\n',Hc)\n",
+"Fs=2\n",
+"H=Cu*M/(G*Fs)\n",
+"printf(' b)H= %f ft',H)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 15.6: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"b=45\n",
+"c=20\n",
+"C=24\n",
+"G=18.9\n",
+"m=0.06\n",
+"Hc=C/(G*m)\n",
+"Cd=G*Hc*m\n",
+"Fc=C/Cd\n",
+"printf('a)Critical height of slope = %f \n',Hc)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 15.7: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"FSs=1\n",
+"c=20\n",
+"G=18.9\n",
+"C=24\n",
+"Hcr=C/(G*tand(c)*0.17)\n",
+"printf('a)Critical height Hc = %f m\n',Hcr)\n",
+"H=10\n",
+"k=C/(G*H*tand(c))\n",
+"Fs=4*tand(c)\n",
+"printf(' b)Fs = %f',Fs)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 15.8: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"W=22.4\n",
+"C=20\n",
+"a=70\n",
+"s=sind(a)\n",
+"c=cosd(a)\n",
+"l=2.924\n",
+"Wn=W*s\n",
+"Wn1=W*c\n",
+"//doing this to all values\n",
+"F1=30.501\n",
+"F2=776.75\n",
+"F3=1638\n",
+"Fs=(F1*C+F3*tand(C))/F2\n",
+"printf('Fs = %f',Fs)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 15.9: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"C=20\n",
+"G=18.5\n",
+"r=0.25\n",
+"H=21.62\n",
+"C=25\n",
+"b= atand(0.5)\n",
+"//from table 15.3 \n",
+"m=1.624\n",
+"n=1.338\n",
+"Fs=m-n*r\n",
+"printf(' The value of Fs for D= 1 is %f',Fs)"
+   ]
+   }
+],
+"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/Principles_Of_Geotechnical_Engineering_by_B_M_Das/16-Soil_Bearing_Capacity_for_Shallow_Foundations.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/16-Soil_Bearing_Capacity_for_Shallow_Foundations.ipynb
new file mode 100644
index 0000000..a56a48e
--- /dev/null
+++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/16-Soil_Bearing_Capacity_for_Shallow_Foundations.ipynb
@@ -0,0 +1,263 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 16: Soil Bearing Capacity for Shallow Foundations"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 16.1: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"c=20\n",
+"// from table 16.1\n",
+"Nc=17.69\n",
+"Nq=7.44\n",
+"Ng=3.64\n",
+"\n",
+"Df=3\n",
+"G=110\n",
+"q=G*Df\n",
+"\n",
+"C=200\n",
+"B=4\n",
+"\n",
+"Qu= C*Nc+q*Nq+G*B*Ng/2\n",
+"\n",
+"Fs=3\n",
+"Qall=Qu/Fs\n",
+"printf('Qa = %f lb/ft^2',Qall)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 16.2: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"G=18.15\n",
+"qa=30000*9.81/1000\n",
+"\n",
+"Nc=57.75\n",
+"Nq=41.44\n",
+"Ng=45.41\n",
+"C=0\n",
+"q=G*1\n",
+"B=1\n",
+"(1.3*C*Nc+q*Nq+0.4*G*B*Ng)*B^2/3 == qa\n",
+"B= sqrt(294.3/(250.7+109.9))\n",
+"printf(' B = %f m',B)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 16.3: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"B=1.2\n",
+"L=1.2\n",
+"c=32\n",
+"C=0\n",
+"Df=1\n",
+"G=16\n",
+"Nq=23.18\n",
+"Ng=22.02\n",
+"Nc=1\n",
+"Lqs=1+0.1*B*(tand(45+c/2))^2/L\n",
+"Lgs=Lqs\n",
+"Lqd=1+0.1*Df*tand(45+c/2)/B\n",
+"Lgd=Lqd\n",
+"Lcs=1\n",
+"Lcd=1\n",
+"Gs=19.5\n",
+"q=0.5*G+0.5*(Gs-9.81)\n",
+"Qu= C*Lcs*Lcd*Nc+q*Lqs*Lqd*Nq+(Gs-9.81)*Lgs*Lgd*B*Ng/2\n",
+"Qa=Qu/3\n",
+"Q=Qa*B^2\n",
+"printf('the gross load = %f kN',Q)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 16.4: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"e=0.1\n",
+"B=1\n",
+"X=B-2*e\n",
+"Y=1.5\n",
+"B1=0.8\n",
+"L1=1.5\n",
+"c=30\n",
+"Df=1\n",
+"Nq=18.4\n",
+"Ng=15.668\n",
+"q=1*18\n",
+"G=18\n",
+"Lqs=1+e*(B1/L1)*(tand(45+c/2))^2\n",
+"Lgs=Lqs\n",
+"Lqd=1+e*(Df/B1)*tand(45+c/2)\n",
+"Lgd=Lqd\n",
+"qu=q*Lqs*Lqd*Nq+Lgs*Lgd*G*B1*Ng/2\n",
+"Qu=qu*B1*L1\n",
+"printf('The magnitude of the gross ultimate load = %f kN',Qu)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 16.5: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"B=1.5\n",
+"Df=0.75\n",
+"e=0.1*B\n",
+"G=17.5\n",
+"c=30\n",
+"C=0\n",
+"q=G*Df\n",
+"Nq=18.4\n",
+"Ng=15.668\n",
+"Lqd=1+0.1*(Df/B)*tand(45+c/2)\n",
+"Lgd=Lqd\n",
+"Quc=q*Nq*Lqd+Lgd*B*Ng/2\n",
+"k=0.8\n",
+"a=1.754\n",
+"Qua=Quc*(1-a*(e/B)^k)\n",
+"printf('The gross ultimate load per unit length = %f kN',Qua)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 16.6: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"Qup=280\n",
+"Bp=0.7 // in m\n",
+"Bf=1.5\n",
+"Quf=Qup*Bf/Bp\n",
+"printf('The ultimate bearing capacity = %f kN/m^2',Quf)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 16.7: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"a=2500\n",
+"//doing for the first values only\n",
+"Bf=4\n",
+"Bp=0.305\n",
+"q=a/Bf^2\n",
+"Sep=4\n",
+"Sef=Sep*(2*Bf/(Bf+Bp))^2\n",
+"printf('Sef = %f mm',Sef)"
+   ]
+   }
+],
+"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/Principles_Of_Geotechnical_Engineering_by_B_M_Das/18-Subsoil_Exploration.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/18-Subsoil_Exploration.ipynb
new file mode 100644
index 0000000..1c6c804
--- /dev/null
+++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/18-Subsoil_Exploration.ipynb
@@ -0,0 +1,105 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 18: Subsoil Exploration"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 18.1: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//solving for z=5 only\n",
+"To=0.275\n",
+"Cn=To^(-0.5)\n",
+"N60=8\n",
+"N160=Cn*N60\n",
+"printf('(N1)60 = %f',N160)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 18.2: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"z=5\n",
+"To=0.275\n",
+"Cn=2/(1+To)\n",
+"N60=8\n",
+"N160=Cn*N60\n",
+"printf('(N1)60 = %f',N160)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 18.3: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"pa=1 // 14.7 lb/in^2 = 1ton/ft^2\n",
+"To=0.275 // ton/ ft^2\n",
+"N60=8\n",
+"c= atand((N60/(12.2+20.3*(To/pa)))^0.34)\n",
+"printf('The average soil friction angle = %f',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/Principles_Of_Geotechnical_Engineering_by_B_M_Das/2-weight_volume_relationships.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/2-weight_volume_relationships.ipynb
new file mode 100644
index 0000000..9da5bd1
--- /dev/null
+++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/2-weight_volume_relationships.ipynb
@@ -0,0 +1,61 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 2: weight volume relationships"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.2: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"//from graph\n",
+"d= 0.15 //mm\n",
+"w= 0.17 //mm\n",
+"a= 0.27 //mm\n",
+"//calculations\n",
+"C= a/d\n",
+"c= w^2/(a*d)\n",
+"//results\n",
+"printf ('uniformity coefficient = % f ',C)\n",
+"printf ('coefficient of gradation = % 2f ',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/Principles_Of_Geotechnical_Engineering_by_B_M_Das/3-weight_volume_relationships.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/3-weight_volume_relationships.ipynb
new file mode 100644
index 0000000..7064ae0
--- /dev/null
+++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/3-weight_volume_relationships.ipynb
@@ -0,0 +1,164 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 3: weight volume relationships"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.2: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"V= 1.2 //m^3\n",
+"M= 2350 //Kg\n",
+"w= 0.086\n",
+"G= 2.71\n",
+"W= 1000 //kg/m^3\n",
+"//calculations\n",
+"R= M/V\n",
+"D= M/((1+w)*V)\n",
+"e= (G*W/D)-1\n",
+"n= e/(e+1)\n",
+"S= (w*G/e)*100\n",
+"v= (M-(M/(1+w)))/W\n",
+"//results\n",
+"printf ('moist density = % f kg/m^3 ',R)\n",
+"printf ('dry density = % f kg/m^3 ',D)\n",
+"printf ('void ratio = % 3f ',e)\n",
+"printf ('porosity = % 3f ',n)\n",
+"printf ('Degree of saturation = % 3f ',S)\n",
+"printf ('volume of water in soil sample = % 3f m^3 ',v)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.3: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"n= 0.4\n",
+"G= 2.68\n",
+"w= 0.12\n",
+"R= 1000 //kg/m^3\n",
+"V= 10 //m^3\n",
+"//calculations\n",
+"d= G*R*(1-n)*(1+w)\n",
+"s= ((1-n)*G+n)*R\n",
+"M= s-d\n",
+"m= M*V\n",
+"//results\n",
+"printf ('mass of water to be added for full saturation = % f kg ',m)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.4: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"d= 16.19 //kN/m^3\n",
+"w= 0.23\n",
+"W= 9.81 //kN/m^3\n",
+"//calculations\n",
+"R= d*(1+w)\n",
+"G= d/(W-d*w)\n",
+"e= w*G\n",
+"//results\n",
+"printf ('satuarated unit weight = % 2f kN/m^3 ',R)\n",
+"printf ('specific gravity = % 2f ',G)\n",
+"printf ('void ratio = % 2f ',e)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.5: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"G= 2.68\n",
+"w= 0.12\n",
+"d= 1794.4 //kg/m^3\n",
+"W= 1000 //kg/m^3\n",
+"emax= 0.75\n",
+"emin= 0.4\n",
+"//calculation\n",
+"e= (G*W*(1+w)/d)-1\n",
+"D= ((emax-e)/(emax-emin))*100\n",
+"//results\n",
+"printf ('relative density of compaction in percentage = % f ',D)"
+   ]
+   }
+],
+"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/Principles_Of_Geotechnical_Engineering_by_B_M_Das/4-plasticity_and_structure_of_soil.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/4-plasticity_and_structure_of_soil.ipynb
new file mode 100644
index 0000000..6fa3c99
--- /dev/null
+++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/4-plasticity_and_structure_of_soil.ipynb
@@ -0,0 +1,61 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 4: plasticity and structure of soil"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 4.1: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"V1= 24.6 //cm^3\n",
+"V2= 15.9 //cm^3\n",
+"M1= 44 //g\n",
+"M2= 30.1 //g\n",
+"W= 1 //g/cm^3\n",
+"//calculations\n",
+"SL= (((M1-M2)/M2)*100)-(((V1-V2)/M2)*W*100)\n",
+"//results\n",
+"printf ('shrinkage limit of the soil = % 2f ',SL)\n",
+""
+   ]
+   }
+],
+"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/Principles_Of_Geotechnical_Engineering_by_B_M_Das/6-soil_compaction.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/6-soil_compaction.ipynb
new file mode 100644
index 0000000..ecda060
--- /dev/null
+++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/6-soil_compaction.ipynb
@@ -0,0 +1,127 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 6: soil compaction"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.2: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"G= 2.6\n",
+"LL= 20\n",
+"P= 20\n",
+"//calclations\n",
+"R= (4804574*G-195.55*(LL)^2+156971*(P)^0.5-9527830)^0.5\n",
+"n= (1.195e-4)*((LL)^2)-1.964*G-(6.617e-5)*(P)+7.651\n",
+"w= %e^n\n",
+"//results\n",
+"printf ('maximum dry density = % f kg/m^3 ',R)\n",
+"printf ('optimum moisture content = % 2f ',w)\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.3: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"do= 1570 //kg/m^3\n",
+"mo= 0.545 //kg\n",
+"M1= 7.59 //kg\n",
+"M2= 4.78 //kg\n",
+"M3= 3.007 //kg\n",
+"w= 0.102 //\n",
+"dmax= 19 //KN/m^3\n",
+"//calculations\n",
+"Ms= M1-M2\n",
+"Mc= Ms-mo\n",
+"Vh= Mc/do\n",
+"Dc= M3/Vh\n",
+"Du= Dc*9.81/1000\n",
+"f= Du/(1+w)\n",
+"Rc= f*100/dmax\n",
+"//results\n",
+"printf ('dry unit weight of compaction in the field = % 2f kN/m^3 ',f)\n",
+"printf ('relative compaction in the field = % f ',Rc)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.4: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"D1= 0.36 //mm\n",
+"D2= 0.52 //mm\n",
+"D5= 1.42 //mm\n",
+"//calculations\n",
+"Sn= 1.7*(sqrt((3/(D5)^2)+(1/(D2)^2)+(1/(D1)^2)))\n",
+"//results\n",
+"printf ('sustainabilty number = % f ',Sn)"
+   ]
+   }
+],
+"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/Principles_Of_Geotechnical_Engineering_by_B_M_Das/7-permiability.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/7-permiability.ipynb
new file mode 100644
index 0000000..b89d2d3
--- /dev/null
+++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/7-permiability.ipynb
@@ -0,0 +1,355 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 7: permiability"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.10: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"k1= 0.302e-7 //cm/sec\n",
+"k2= 0.12e-7 //cm/sec\n",
+"e1= 1.1\n",
+"e2= 0.9\n",
+"e= 0.75\n",
+"//calcualtions\n",
+"n= (log10((k1/k2)*((1+e1)/(1+e2))))/log10(e1/e2)\n",
+"C= k1/(e1^n/(1+e1))\n",
+"k= C*(e^n/(1+e))\n",
+"//results\n",
+"printf ('hydraulic conductivity = %e cm/sec ',k')"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.11: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"H1= 2 //m\n",
+"H2= 3 //m\n",
+"H3= 4 //m\n",
+"k1= 1e-4 //cm/sec\n",
+"k2= 3.2e-2 //cm/sec\n",
+"k3= 4.1e-5 //cm/sec\n",
+"//calculations\n",
+"H= H1+H2+H3\n",
+"Kh= (1/H)*((k1*H1)+(k2*H2)+(k3*H3))\n",
+"Kv= H/((H1/k1)+(H2/k2)+(H3/k3))\n",
+"P= Kh/Kv\n",
+"//results\n",
+"printf ('ration of equivalent hydraulic conductivity = % 2f ',P)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.12: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"H= 450 //mm\n",
+"h= 150 //mm\n",
+"k1= 1e-2 //cm/sec\n",
+"k2= 3e-3 //cm/sec\n",
+"k3= 4.9e-4 //cm/sec\n",
+"h1= 300 //mm\n",
+"//calculations\n",
+"Kv= H/(h*(1/k1+1/k2+1/k3))\n",
+"i= h1/H\n",
+"q= Kv*i*100*3600\n",
+"//results\n",
+"printf ('rate of water supply = % 2f cm/hr ',q)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.1: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"L= 30 //cm\n",
+"A= 177 //cm^2\n",
+"h= 50 //cm\n",
+"Q= 350 //cm^3\n",
+"t= 300 //sec\n",
+"//claculations\n",
+"k=Q*L/(A*h*t)\n",
+"//results\n",
+"printf ('hydraulic conductivity = % 3f cm/sec ',k)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.2: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"L= 203 //mm\n",
+"A= 10.3 //cm^2\n",
+"a= 0.39 //cm^2\n",
+"h0= 508 //mm\n",
+"h180= 305 //mm\n",
+"t= 180 //sec\n",
+"//calculations\n",
+"k= 2.303*a*L*log10(h0/h180)/(A*t)\n",
+"//results\n",
+"printf ('hydraulic conductivity of sand = % 2f in/sec ',k)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.3: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of varilables\n",
+"k= 3e-7 //cm/sec\n",
+"n= 0.0911e-4 //g*sec/cm^2\n",
+"dw= 1 //g/cc\n",
+"//calculations\n",
+"K= k*n/dw\n",
+"//results\n",
+"printf ('absolute premeability = % 4f cm^2 ',K)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.4: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"k= 5.3e-5 //m/sec\n",
+"H= 3 //m\n",
+"a= 0.139 //radians\n",
+"//calculations\n",
+"A= H*cos(a)\n",
+"i= sin(a)\n",
+"q= k*i*A*3600\n",
+"//results\n",
+"printf ('rate of seepage = % 4f m^3/hr/m ',q)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.5: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"L= 50 //m\n",
+"k= 0.08e-2//m/sec\n",
+"h= 4 //m\n",
+"H1= 3 //m\n",
+"H= 8 //m\n",
+"a= 0.139 //radians\n",
+"//calculations\n",
+"i= h*cos(a)/L\n",
+"A= H1*cos(a)\n",
+"q= k*i*A\n",
+"//results\n",
+"printf ('flow rate = % 2f m^3/sec/m ',q)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.6: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"k1= 0.02 //cm/sec\n",
+"e1= 0.5 \n",
+"e2= 0.65\n",
+"//calculations\n",
+"k2= k1*(e2^3/(1+e2))/(e1^3/(1+e1))\n",
+"//results\n",
+"printf ('hydraulic conductivity at void ratio of 0.65 = % 2f cm/sec ',k2)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.8: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"e= 0.6\n",
+"D10= 0.09 //mm\n",
+"//calculations\n",
+"k= 2.4622*(D10^2*(e^3/(1+e)))^0.7825\n",
+"//results\n",
+"printf ('hydraulic conductivity = % 4f cm/sec ',k)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.9: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"e= 0.6\n",
+"D10= 0.09 //mm\n",
+"D60= 0.16 //mm\n",
+"//calculations\n",
+"Cu=D60/D10\n",
+"k= 35*(e^3/(1+e))*(Cu^0.6)*(D10^2.32)\n",
+"//results\n",
+"printf ('hydraulic conductivity = % 3f cm/sec ',k)"
+   ]
+   }
+],
+"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/Principles_Of_Geotechnical_Engineering_by_B_M_Das/8-seepage.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/8-seepage.ipynb
new file mode 100644
index 0000000..99bf4d9
--- /dev/null
+++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/8-seepage.ipynb
@@ -0,0 +1,200 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 8: seepage"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.1: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"H1= 12 //in\n",
+"H2= 20 //in\n",
+"z= 8 //in\n",
+"h1= 24 //in\n",
+"h= 20 //in\n",
+"k1= 0.026 //in/sec\n",
+"D= 3 //in\n",
+"//calculations\n",
+"k2= H2*k1/((z/(1-h/h1))-H1)\n",
+"i= h1/(H1+H2)\n",
+"A= %pi/4*D^2\n",
+"keq= (H1+H2)/((H1/k1)+(H2/k2))\n",
+"q= keq*A*i*3600\n",
+"//results\n",
+"printf ('rate of water flow = % 2f in^3/hr ',q)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.2: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"Nd= 6\n",
+"H1= 5.6 //m\n",
+"H2= 2.2 //m\n",
+"k= 5e-5 //cm/sec\n",
+"dL= 4.1 //m\n",
+"//calculations\n",
+"H= (H1-H2)/Nd\n",
+"h1= 5.61-H\n",
+"h2= 5.61-5*H\n",
+"q= 2.38*(H1-H2)*k/Nd\n",
+"i= H/dL\n",
+"//results\n",
+"printf ('at point a,water will rise to height of = % 3f m ',h1)\n",
+"printf ('at point b,water will rise to height of = % 3f m ',h2)\n",
+"printf ('total rate of seepage per unit lenghth = %e m^3/sec/m ',q)\n",
+"printf ('average hydraulic gradient at c = % 3f ',i)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.3: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"k1= 5.67 //ft/day\n",
+"k2= 11.34 //ft/day\n",
+"//from graph\n",
+"Nd= 8\n",
+"Nf= 2.5\n",
+"H= 20\n",
+"//calculations\n",
+"q= sqrt(k1*k2)*H*Nf/Nd\n",
+"//results\n",
+"printf ('average rate of flow = % 2f ft^3/day/ft ',q)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.4: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"B= 6 //m\n",
+"L= 120 //m\n",
+"s= 3 //m\n",
+"T= 6 //m\n",
+"x= 2.4 //m\n",
+"H= 5 //m\n",
+"k= 0.008 //cm/sec\n",
+"//calculations\n",
+"b=B/2\n",
+"a1= b/T\n",
+"a2= s/T\n",
+"a3= x/b\n",
+"Q= 0.378*k*H*L*36*24\n",
+"//results\n",
+"printf ('seepage under the dam = % 2f m^3/day ',Q)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.5: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"b= %pi/4 //degrees\n",
+"a= %pi/6.//degrees\n",
+"B= 10 //ft\n",
+"H= 20 //ft\n",
+"h= 25 //ft\n",
+"k= 2e-4 //ft/min\n",
+"//calculations\n",
+"r= H/tan(b)\n",
+"d= 0.3*r+(h-H)/tan(b)+B+h/tan(a)\n",
+"L= d/cos(a)-sqrt((d/cos(a))^2-(H/sin(a))^2)\n",
+"q= k*L*tan(a)*sin(a)*24*60\n",
+"//results\n",
+"printf ('seepage rate = % 4f ft^3/day/ft ',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/Principles_Of_Geotechnical_Engineering_by_B_M_Das/9-in_situ_stresses.ipynb b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/9-in_situ_stresses.ipynb
new file mode 100644
index 0000000..5423a71
--- /dev/null
+++ b/Principles_Of_Geotechnical_Engineering_by_B_M_Das/9-in_situ_stresses.ipynb
@@ -0,0 +1,176 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 9: in situ stresses"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.1: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"Ds= 16.5 //kN/m^3\n",
+"S= 19.25 //kN/m^3\n",
+"g= 9.8 //kN/m^3\n",
+"h1= 6 //m\n",
+"h2= 13 //m\n",
+"//at point A\n",
+"Sa= 0\n",
+"Ua= 0\n",
+"Sa1= 0\n",
+"//at point B\n",
+"Sb= h1*Ds\n",
+"Ub= 0 \n",
+"Sb1= Sb-Ub\n",
+"//at point C\n",
+"Sc= h1*Ds+h2*S\n",
+"Uc= h2*g\n",
+"Sc1= Sc-Uc\n",
+"//results\n",
+"printf ('total pressure at C= % 2f kN/m^3 ',Sc)\n",
+"printf ('pore water pressure at C = % 2f kN/m^3 ',Uc)\n",
+"printf ('effective stress at point C= % 2f kN/m^3 ',Sc1)\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.2: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"h= 20 //ft\n",
+"g= 120 //kg/ft^3\n",
+"h1= 12 //ft\n",
+"w= 62.4 //kg/ft^3\n",
+"//calculations\n",
+"H= h-(h1*w/g)\n",
+"//results\n",
+"printf ('maximu depth that can be made in clay = % 2f ft ',H)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.3: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"G= 2.68\n",
+"e= 0.52\n",
+"g= 9.81 //kN/m^3\n",
+"h1= 0.7 //m\n",
+"h2= 1 //m\n",
+"h3= 1.5 //m\n",
+"h4= 2 //m\n",
+"//calculations\n",
+"//for soil A\n",
+"sa= (G+e)*g/(1+e)\n",
+"//point a\n",
+"Sa= h1*g+h2*sa\n",
+"u= (h2+h1+h3/2)*g\n",
+"Es= Sa-u\n",
+"//point b\n",
+"sb= h1*g+h4*sa\n",
+"ub= (h4+h1+h3)*g\n",
+"Eb= sb-ub\n",
+"i= h3/2\n",
+"s= i*g\n",
+"//results\n",
+"printf ('effective stress at point a= % 2f kN/m^2 ',Es)\n",
+"printf ('effective stress at point b= % 2f kN/m^2 ',Eb)\n",
+"printf ('seepage force per unit voume = % 2f kN/m^3 ',s)"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.4: solved.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"clc\n",
+"//initialisation of variables\n",
+"C0= 0.357\n",
+"H1= 30.5 //ft\n",
+"H2= 5 //ft\n",
+"w= 62.4 // lb/ft^3\n",
+"D= 20\n",
+"g= 112 // lb/ft^3\n",
+"//calculations\n",
+"G= g-w\n",
+"FS= D*G/(C0*w*(H1-H2))\n",
+"//results\n",
+"printf ('safety factor = % 2f ',FS)"
+   ]
+   }
+],
+"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
+}