Merge pull request #1 from prashantsinalkar/masterHEADmaster

Initial commit

14 files changed, 4110 insertions, 0 deletions

diff --git a/Principles_Of_Foundation_Engineering_by_B_M_Das/1-Geotechnical_Properties_of_Soil.ipynb b/Principles_Of_Foundation_Engineering_by_B_M_Das/1-Geotechnical_Properties_of_Soil.ipynb
new file mode 100644
index 0000000..d7dc262
--- /dev/null
+++ b/Principles_Of_Foundation_Engineering_by_B_M_Das/1-Geotechnical_Properties_of_Soil.ipynb
@@ -0,0 +1,291 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 1: Geotechnical Properties of Soil"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 1.1: 1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 1.1\n",
+"clc(); funcprot(0);\n",
+"V=0.25; // ft^3\n",
+"W=30.8; //lb\n",
+"Wd=28.2; // weight dried lb\n",
+"Gs=2.67;\n",
+"Gammaw=62.4; //lb/ft^3\n",
+"Gamma=W/V;\n",
+"disp(Gamma, 'moist unit weight in lb/ft^3')\n",
+"w=(W-Wd)/W;\n",
+"disp(w*100, 'moisture content in %');\n",
+"Gammad=Wd/V;\n",
+"disp(Gammad, 'dry unit weight in lb/ft^3')\n",
+"Vs=Wd/Gs/Gammaw;\n",
+"Vv=V-Vs;\n",
+"e=Vv/Vs;\n",
+"disp(e,'void ratio');\n",
+"n=e/(1+e);\n",
+"disp(n,'porosity');\n",
+"Vw=(W-Wd)/Gammaw;\n",
+"S=Vw/Vv;\n",
+"disp(S*100,'saturation in %');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 1.2: 2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 1.2\n",
+"clc; funcprot(0);\n",
+"e=0.72;\n",
+"w=12/100; //moisture content\n",
+"Gs=2.72;\n",
+"Gammaw=9.81;//kN/m^3\n",
+"Gammad=Gs*Gammaw/(1+e);\n",
+"disp(Gammad,'dry unit weight in kN/m^3');\n",
+"Gamma=Gs*Gammaw*(1+w)/(1+e);\n",
+"disp(Gamma,'moist unit weight in kN/m^3');\n",
+"Gammasat=(Gs+e)*Gammaw/(1+e);\n",
+"wa=Gammasat-Gamma;//water added\n",
+"disp(wa,'water added in kN/m^3');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 1.3: 3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 1.3\n",
+"clc; funcprot(0);\n",
+"gmax=17.1; // Gammadmax\n",
+"Dr=0.7;\n",
+"w=8/100;//moisture content\n",
+"gmin=14.2; //Gammadmin\n",
+"deff('y=f(x)','y=(x-14.2)/(17.1-14.2)*17.1/x-0.7')\n",
+"[x, v, info]=fsolve(16,f);//solving for gammad\n",
+"Gamma=x*(1+w);\n",
+"disp(Gamma,'moist unit weight in kN/m^3');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 1.4: 4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 1.4\n",
+"clc; funcprot(0);\n",
+"F200=30;\n",
+"PI=8;\n",
+"GI=0.01*(F200-15)*(PI-10);\n",
+"if GI<0 then\n",
+"    GI=0;\n",
+"    disp('the soil is A-2-4(0)')\n",
+"end\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 1.5: 5.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 1.5\n",
+"clc; funcprot(0);\n",
+"F200=41;\n",
+"LL=31;\n",
+"PI=12;\n",
+"if 100-F200>50 then\n",
+"    disp('soil is coarse grained')\n",
+"end\n",
+"//from table 1.5 and 1.6\n",
+"disp('soil is SC');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 1.6: 6.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 1.6\n",
+"clc; funcprot(0);\n",
+"Gammad=14.5;\n",
+"Gammasat=17.2;\n",
+"Gammaw=9.81;\n",
+"printf('Point \t Sigma(kN/m^2)\t u(kN/m^2)\t sigmadash=sigma-u (kN/m^2)\n');\n",
+"\n",
+"printf('  A\t   %d\t \t    %d\t\t %d \n',0,0,0);\n",
+"printf('  B\t%.2f\t \t%.2f\t\t %.2f \n',4*Gammad,0,4*Gammad-0);\n",
+"printf('  C\t%.2f\t \t%.2f\t\t %.2f \n',4*Gammad+5*Gammasat,5*Gammaw,4*Gammad+5*Gammasat-5*Gammaw);\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 1.7: 7.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 1.7\n",
+"clc; funcprot(0);\n",
+"//parta\n",
+"e1=0.92;\n",
+"e2=0.86;\n",
+"Hc=2.8;\n",
+"s2=212;//sigma2dash\n",
+"s1=140;//sigma1dash\n",
+"Cc=(e1-e2)/log10(s2/s1);\n",
+"Sc=Cc*Hc/(1+e1)*log10(s2/s1);\n",
+"disp(Sc*1000,'consolidated depth in mm');\n",
+"//part b\n",
+"Sct=40;\n",
+"T50=0.197;\n",
+"t=4.5;\n",
+"Cr=T50*12.7^2/t;\n",
+"U=Sct/Sc*100/1000;\n",
+"H=Hc/2;\n",
+"Tv=%pi/4*U^2/100^2;\n",
+"t=Tv*H^2/Cr*1000^2/60/24;\n",
+"disp(t,'time required in days');\n",
+"\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 1.8: 8.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 1.8\n",
+"clc; funcprot(0);\n",
+"Cv=7.061;\n",
+"tc=15*24*60;\n",
+"H=2.8/2*1000;\n",
+"Scmax=87.5; // consolidation\n",
+"Tc=Cv*tc/H^2;\n",
+"tv=31.6*24*60;\n",
+"Tv=Cv*tv/H^2;\n",
+"//from figure 1.28\n",
+"Sct=Scmax*0.36;\n",
+"disp(Sct,'consolidation in 31.6 days in mm');\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_Foundation_Engineering_by_B_M_Das/10-Braced_Cuts.ipynb b/Principles_Of_Foundation_Engineering_by_B_M_Das/10-Braced_Cuts.ipynb
new file mode 100644
index 0000000..fdbf711
--- /dev/null
+++ b/Principles_Of_Foundation_Engineering_by_B_M_Das/10-Braced_Cuts.ipynb
@@ -0,0 +1,182 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 10: Braced Cuts"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 10.1: 1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 10.1\n",
+"clc; funcprot(0);\n",
+"Gamma=18;\n",
+"H=7;\n",
+"sigmaa=0.3*Gamma*H;\n",
+"disp(sigmaa,'maximum pressure intensity in kN/m^2');\n",
+"//partb\n",
+"A=54.02;\n",
+"B1=1/2*1.75*37.8+37.8*1.75-A;\n",
+"B2=45.2;\n",
+"C=54.02;\n",
+"s=3; //spacing\n",
+"Pa=C*s;\n",
+"disp(Pa,'strut loads in kN');\n",
+"Pb=(B1+B2)*s;\n",
+"disp(Pb,'strut loads in kN')\n",
+"Pc=C*s;\n",
+"disp(Pc,'strut loads in kN');\n",
+"//partc\n",
+"Me=45.2*1.196-37.8*1.196*1.196/2;//Me=Mmax\n",
+"Sall=170e3;//sigmaall\n",
+"S=Me/Sall;\n",
+"disp(S,'section modulus in m^3/m');\n",
+"//partd\n",
+"Mmax=(B1+B2)*s^2/8;\n",
+"S=Mmax/Sall;\n",
+"disp(S,'section modulus in m^3/m');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 10.2: 2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 10.2\n",
+"clc; funcprot(0);\n",
+"phi=32;\n",
+"Gamma=112;\n",
+"s=12;//spacing\n",
+"H=27;\n",
+"Ka=(tan(45*%pi/180-32*%pi/180/2))^2;\n",
+"sigmaa=0.65*Gamma*Ka*H;\n",
+"A=sigmaa*15*15/2/9;\n",
+"B1=603.44*15-A;\n",
+"C=sigmaa*12*12/2/9; \n",
+"B2=sigmaa*s-C;\n",
+"Pa=A*s/1000;\n",
+"Pb=(B1+B2)*s/1000;\n",
+"Pc=C*s/1000;\n",
+"disp(Pa,'strut loads at A in kN');\n",
+"disp(Pb,'strut loads at B in kN');\n",
+"disp(Pc,'strut loads at C in kN');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 10.3: 3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 10.3\n",
+"clc; funcprot(0);\n",
+"//parta\n",
+"s=12;\n",
+"A=1/2*6*3.62;\n",
+"C=1/2*1.81*3;\n",
+"B1=1/2*2.5*1.509;//B'\n",
+"B11=1/2*2.414*4;//B''\n",
+"sigmaall=24;\n",
+"Mmax=max(A,C,B1,B11);\n",
+"S=Mmax*s/sigmaall;\n",
+"disp(S,'section modulus in in^3/ft');\n",
+"//partb\n",
+"\n",
+"A=7.583;\n",
+"Mmax=7.583*s^2/8;\n",
+"S=Mmax*s/sigmaall;\n",
+"disp(S,'section modulus in in^3/ft');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 10.4: 4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 10.4\n",
+"clc; funcprot(0);\n",
+"q=0;\n",
+"Gamma=18;\n",
+"L=25;\n",
+"c=35;\n",
+"T=4;\n",
+"H=7;\n",
+"B1=4;//B'\n",
+"B11=T*sqrt(2);//B''\n",
+"FS=(5.14*c*(1+B11*0.2/L)+c*H/B1)/(Gamma*H+q)\n",
+"disp(FS,'factor of safety is')\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_Foundation_Engineering_by_B_M_Das/11-Pile_Foundations.ipynb b/Principles_Of_Foundation_Engineering_by_B_M_Das/11-Pile_Foundations.ipynb
new file mode 100644
index 0000000..91498ec
--- /dev/null
+++ b/Principles_Of_Foundation_Engineering_by_B_M_Das/11-Pile_Foundations.ipynb
@@ -0,0 +1,527 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 11: Pile Foundations"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.10: 10.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 11.10\n",
+"clc; funcprot(0);\n",
+"Hp=350;\n",
+"vp=0.0016;\n",
+"Sl=0.762e-3;\n",
+"f=115;\n",
+"Qu=(0.746*Hp+98*vp)/(vp+Sl*f);\n",
+"disp(Qu,'ple load capacity in kN');"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.11: 11.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 11.11\n",
+"clc; funcprot(0);\n",
+"Lg=9.92;\n",
+"Bg=7;\n",
+"n1=3;\n",
+"Nc=8.75;\n",
+"n2=4/1000;\n",
+"Ap=14^2/12^2;\n",
+"cup=1775;\n",
+"a1=0.4;//alpha1\n",
+"p=4*14/12;\n",
+"cu1=1050;\n",
+"L1=15;\n",
+"a2=0.54;//alpha2\n",
+"cu2=1775;\n",
+"L2=45;\n",
+"FS=4;\n",
+"Qu=n1*n2*(9*Ap*cup+a1*p*cu1*L1+a2*p*cu2*L2);\n",
+"Qu2=Lg*Bg*cup*Nc+2*(Lg+Bg)*(cu1*L1+cu2*L2);\n",
+"disp(Qu2/1000,'load in kip')\n",
+"Qall=Qu/FS;\n",
+"disp(Qall,'allowed load in kip');\n",
+"\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.1: 1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 11.1\n",
+"clc; funcprot(0);\n",
+"//parta\n",
+"phi=30;\n",
+"pa=2000;\n",
+"q=100*50/1000;\n",
+"Nq=55;\n",
+"Ap=16*16/16/12;\n",
+"Qp=Ap*q*Nq;\n",
+"qp=0.4*pa*Nq*tan(phi*%pi/180)*Ap;\n",
+"disp(qp,'ultimate load in lb');\n",
+"disp(qp/1000,'ultimate load in kip');\n",
+"disp('there is change in answer because of calculation mistake in the book');\n",
+"//partb\n",
+"Nsigma=36;\n",
+"Ap=16*16/12/12;\n",
+"q=110*50/1000;\n",
+"Qp=Ap*q*Nsigma*((1+2*(1-sin(phi*%pi/180)))/3);\n",
+"disp(Qp,'ultimate load in kip');\n",
+"//partc\n",
+"Nq=18.4;\n",
+"Qp=Ap*q*Nq;\n",
+"disp(Qp,'ultimate load in kip');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.12: 12.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 11.12\n",
+"clc; funcprot(0);\n",
+"z1=21/2;\n",
+"Lg=9;\n",
+"Bg=6;\n",
+"Qg=500*1000;\n",
+"Cc1=0.3;\n",
+"Cc2=0.2;\n",
+"Cc3=0.25;\n",
+"H2=12;\n",
+"H3=6;\n",
+"H1=21;\n",
+"e1=0.82;\n",
+"e2=0.7;\n",
+"e3=0.75;\n",
+"s1=Qg/(Lg+z1)/(Bg+z1); //sigma1\n",
+"s2=500*1000/(9+27)/(6+27);//sigma2\n",
+"s3=500*1000/(9+36)/(6+36);//sigma3\n",
+"ss1=6*105+(27+21/2)*(115-62.4);//sigmadash1\n",
+"ss2=6*105+(27+21)*(115-62.4)+(120-62.4)*6;//sigmadash2\n",
+"ss3=6*105+48*(115-62.4)+12*(120-62.4)+3*(122-62.4);//sigmadash3\n",
+"sc1=Cc1*H1/(1+e1)*log10((ss1+s1)/ss1);\n",
+"sc2=Cc2*H2/(1+e2)*log10((ss2+s2)/ss2);\n",
+"sc3=Cc3*H3/(1+e3)*log10((ss3+s3)/ss3);\n",
+"sc=sc1+sc2+sc3;\n",
+"disp(sc*12,'total settlement in inch');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.2: 2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 11.2\n",
+"clc; funcprot(0);\n",
+"//parta\n",
+"K=1.3;\n",
+"f0=0;\n",
+"Delta=0.8*30;\n",
+"D=16/12;\n",
+"L1=50;\n",
+"p=4*16/12;\n",
+"Gamma=110/1000;\n",
+"L=15*D;\n",
+"sigma=Gamma*L;\n",
+"f20=K*sigma*tan(Delta*%pi/180);\n",
+"Qs=(f0+f20)/2*(p*L)+f20*p*(L1-L);\n",
+"disp(Qs,'ultimate load in kip');\n",
+"//partb\n",
+"FS=4;\n",
+"Qp=56.45/3+234.7/3+179.9/3;\n",
+"Qu=Qs+Qp;\n",
+"Qall=Qu/FS;\n",
+"disp(Qall,'allowed load in kip');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.3: 3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 11.3\n",
+"clc; funcprot(0);\n",
+"K=0.25;\n",
+"Ap=16*16/12/12;\n",
+"phi=30*%pi/180;\n",
+"Nq=25;\n",
+"q=110*50/1000;\n",
+"sigmao=q/2;\n",
+"p=4*16/12;\n",
+"L=50;\n",
+"FS=4;\n",
+"Qu=q*Nq*Ap+K*sigmao*tan(0.8*phi)*p*L;\n",
+"Qall=Qu/FS;\n",
+"disp(Qall,'allowed load in kip');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.4: 4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 11.4\n",
+"clc; funcprot(0);\n",
+"FS=4;\n",
+"Ap=0.1295;\n",
+"Nc=9;\n",
+"cu2=100;\n",
+"Qp=Ap*Nc*cu2;\n",
+"D=[5, 10, 30];\n",
+"avgD=[2.5, 7.5,20];\n",
+"sigma=[45, 110.5, 228.5];\n",
+"cu=[30, 30, 100];\n",
+"alpha=[0.6 0.9 0.725];\n",
+"L=[5, 5, 20];\n",
+"p=%pi*0.406;\n",
+"Qs=0;\n",
+"disp(Qp,'bering capacity in kN');\n",
+"printf('depth (m)\t avg Depth(m)\t avgVerticalStress(kN/m^2)\t cu(kN/m^2)\t cu/sigma\t alpha\n' )\n",
+"for i=1:3\n",
+"    cusig(i)=cu(i)/sigma(i);\n",
+"    Qs=Qs+alpha(i)*cu(i)*L(i)*p;\n",
+"    printf('%.2f\t          %.2f\t           %.2f\t\t\t   %.2f\t   %.2f\t\t %.2f\n',D(i),avgD(i),sigma(i),cu(i),cusig(i),alpha(i));\n",
+"end\n",
+"disp(Qs,'bearing capacity in kN');\n",
+"//part2\n",
+"lambda=0.136;\n",
+"L=30;\n",
+"fav=lambda*(178.48+2*76.7);\n",
+"Qs2=p*L*fav;\n",
+"//part3\n",
+"fav1=13;\n",
+"fav2=31.9;\n",
+"fav3=93.43;\n",
+"Qs3=p*(fav1*5+fav2*5+fav3*20);\n",
+"disp(Qs3,'bearing capacity in kN');\n",
+"Qsavg=Qs/3+Qs2/3+Qs3/3;\n",
+"Qu=Qp+Qsavg\n",
+"Qall=Qu/FS;\n",
+"disp(Qall,'allowed bearing capacity in kN');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.5: 5.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 11.5\n",
+"clc; funcprot(0);\n",
+"D=[6, 12, 20];\n",
+"fc=[34.34, 54.94, 70.63];\n",
+"alpha=[0.84, 0.71, 0.63];\n",
+"dL=[6, 6, 8];\n",
+"p=4*0.305;\n",
+"Qs=0;\n",
+"printf('  depth(m)\t fc(kN/m^2)\t alpha \t deltaL(m)\t Q(kN)\n');\n",
+"for i=1:3\n",
+"    Q(i)=alpha(i)*fc(i)*p*dL(i);\n",
+"    Qs=Q(i)+Qs;\n",
+"    printf('%.2f\t\t    %.2f\t  %.2f\t   %.2f\t\t %.2f\n',D(i),fc(i),alpha(i),dL(i),Q(i));\n",
+"end\n",
+"disp(Qs,'bearing force in kN');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.6: 6.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 11.6\n",
+"clc; funcprot(0);\n",
+"L=21;\n",
+"Qwp=502-350;\n",
+"Qws=350;\n",
+"Ap=0.1045;\n",
+"Ep=21e6;\n",
+"epsilon=0.62;\n",
+"Se1=(Qwp+epsilon*Qws)*L/Ap/Ep;\n",
+"//part2\n",
+"Iwp=0.85;\n",
+"qwp=152/Ap;\n",
+"Es=25e3;\n",
+"D=0.356;\n",
+"mus=0.35;\n",
+"Se2=qwp*D/Es*Iwp*(1-mus^2);\n",
+"//part3\n",
+"p=1.168;\n",
+"Iws=2+0.35*sqrt(L/D);\n",
+"Se3=Qws/p/L*D/Es*Iws*(1-mus^2);\n",
+"Se=Se1+Se2+Se3;\n",
+"disp(Se*1000,'settlement in mm');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.7: 7.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 11.7\n",
+"clc; funcprot(0);\n",
+"Ep=207e6;\n",
+"Ip=123e-6;\n",
+"nh=12000;\n",
+"//from table\n",
+"xz=0.008;\n",
+"Ax=2.435;\n",
+"T=(Ep*Ip/nh)^0.2;\n",
+"Qg1=xz*Ep*Ip/Ax/T^3;\n",
+"//part2\n",
+"Fy=248000;\n",
+"d1=0.254;\n",
+"Am=0.772;\n",
+"Mzmax=Fy*Ip*2/d1;\n",
+"Qg2=Mzmax/Am/T;\n",
+"if Qg2>Qg1 then\n",
+"    Qg=Qg1;\n",
+"    disp(Qg,'lateral load in kN');\n",
+"end\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.8: 8.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 11.8\n",
+"clc; funcprot(0);\n",
+"//part1\n",
+"Ep=207e6;\n",
+"Ip=123e-6;\n",
+"nh=12000;\n",
+"//from table\n",
+"xo=0.008;\n",
+"L=25;\n",
+"Fy=248000;\n",
+"D=0.254;\n",
+"Am=0.772;\n",
+"Gamma=18;\n",
+"phi=35;\n",
+"Kp=(tan(%pi/4+phi*%pi/360))^2;\n",
+"My=Fy*Ip*2/D;\n",
+"Qug=140*Kp*D^3*Gamma;\n",
+"\n",
+"//part2\n",
+"Qg1=xo*(Ep*Ip)^0.6*nh^0.4/0.15/L;\n",
+"\n",
+"if Qug>Qg1 then\n",
+"    Qg=Qg1;\n",
+"    disp(Qg,'lateral load in kN');\n",
+"end\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 11.9: 9.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 11.9\n",
+"clc; funcprot(0);\n",
+"Wrh=30*12;\n",
+"E=0.8;\n",
+"Wr=7.5;\n",
+"S=1/8;\n",
+"C=0.1;\n",
+"FS=6;\n",
+"n=0.4;\n",
+"Wp=12/12*12/12*80*150+550;\n",
+"Wp=Wp/1000;\n",
+"Qu=E*Wrh/(S+C)*(Wr+n^2*Wp)/(Wr+Wp);\n",
+"Qall=Qu/FS;\n",
+"disp(Qall,'allowed bearing capacity in kip')\n",
+"//part2\n",
+"He=30*12;\n",
+"L=80*12;\n",
+"Ap=12*12;\n",
+"Ep=3e6/1000;\n",
+"FS=4;\n",
+"Qu=E*He/(S+sqrt(E*He*L/2/Ap/Ep));\n",
+"Qall2=Qu/FS;\n",
+"disp(Qall2,'allowed bearing capacity in kip')\n",
+"\n",
+"//partc\n",
+"a=27;\n",
+"b=1;\n",
+"He=30;\n",
+"FS=3;\n",
+"Qu=a*sqrt(E*He)*(b-log10(S));\n",
+"Qall3=Qu/FS;\n",
+"disp(Qall3,'allowed bearing capacity in kip')\n",
+"\n",
+"\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_Foundation_Engineering_by_B_M_Das/12-Drift_Shaft_Foundations.ipynb b/Principles_Of_Foundation_Engineering_by_B_M_Das/12-Drift_Shaft_Foundations.ipynb
new file mode 100644
index 0000000..7dd4632
--- /dev/null
+++ b/Principles_Of_Foundation_Engineering_by_B_M_Das/12-Drift_Shaft_Foundations.ipynb
@@ -0,0 +1,356 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 12: Drift Shaft Foundations"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 12.1: 1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 12.1\n",
+"clc; funcprot(0);\n",
+"Ap=%pi/4*1.75^2;\n",
+"FS=4;\n",
+"Nq=37.75;\n",
+"L=8;\n",
+"Es=50000;\n",
+"mus=0.265;\n",
+"pu=100;\n",
+"Db=1.75;\n",
+"q=6*16.2+2*19.2;\n",
+"phi=36*%pi/180;\n",
+"Fqs=1+tan(phi);\n",
+"Fqd=1+2*tan(phi)*(1-sin(phi))^2*atan(L/Db);\n",
+"Ir=Es/(2*(1+mus)*q*tan(phi));\n",
+"delta=0.005*(1-phi/20*180/%pi+25/20)*q/pu;\n",
+"Irr=Ir/(1+Ir*delta);\n",
+"Fqc=exp(-3.8*tan(phi)+(3.07*sin(phi)*log10(2*Irr))/(1+sin(phi)));\n",
+"Qp=Ap*(q*(Nq-1)*Fqs*Fqd*Fqc);\n",
+"Qpall=Qp/FS;\n",
+"disp(Qpall,'allowed load in kN');\n",
+"disp('due to rounding off error there is slight change in answer')\n",
+"\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 12.2: 2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 12.2\n",
+"clc; funcprot(0);\n",
+"Ap=%pi/4*1.75^2;\n",
+"q=135.6;\n",
+"w=0.83;\n",
+"FS=4;\n",
+"phi=36;\n",
+"Nq=0.21*exp(0.17*phi);\n",
+"Qp=Ap*q*(w*Nq-1);\n",
+"Qpall=Qp/FS;\n",
+"disp(Qpall,'allowed load in kN');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 12.3: 3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 12.3\n",
+"clc; funcprot(0);\n",
+"Ap=%pi/4*1.5^2;\n",
+"Db=1.5;\n",
+"z=3;\n",
+"p=%pi*1;\n",
+"Li=6;\n",
+"N60=30;\n",
+"sigmazi=16*z;\n",
+"Beta=2.0-0.15*z^0.75;\n",
+"fi=Beta*sigmazi;\n",
+"qp=57.5*N60;\n",
+"qpr=1.27/Db*qp;\n",
+"Qunet=qpr*Ap+fi*p*Li;\n",
+"disp(Qunet,'allowed load in kN');\n",
+"//part b\n",
+"k1=0.315; //from table\n",
+"k2=12/1.5/1000*100;\n",
+"Qunet2=qpr*Ap*k1+fi*p*Li*k2;\n",
+"disp(Qunet2,'allowed load in kN');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 12.4: 4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 12.4\n",
+"clc; funcprot(0);\n",
+"Nc=9;\n",
+"Ap=%pi/4*1.5^2;\n",
+"cu=105;\n",
+"Qpnet=Ap*cu*Nc;\n",
+"disp(Qpnet,'net ultimate bearing point capacity in kN');\n",
+"//part2\n",
+"alpha=0.4;\n",
+"Ds=1.5;\n",
+"p=%pi*Ds;\n",
+"Qs=alpha*p*(50*8+105*3);\n",
+"disp(Qs,'skin resistance in kN');\n",
+"//part3\n",
+"FS=3;\n",
+"Qu=Qpnet/FS+Qs/FS;\n",
+"disp(Qu,'working load in kN');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 12.5: 5.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 12.5\n",
+"clc; funcprot(0);\n",
+"cub=3000;\n",
+"L=20+5;\n",
+"Db=4;\n",
+"Ap=%pi/4*Db^2;\n",
+"alpha=0.55;\n",
+"cu1=800;\n",
+"L1=7;\n",
+"L2=5.5;\n",
+"cu2=1200;\n",
+"p=%pi*2.5;\n",
+"k=alpha*p*(cu1*L1+cu2*L2);//f*p*deltaLi\n",
+"j1=6*cub*(1+0.2*L/Db);\n",
+"j2=9*cub;\n",
+"[qp]=min(j1,j2);\n",
+"Qu=k/1000+qp*Ap/1000;\n",
+"disp(Qu,'allowed load in kN');\n",
+"//part b\n",
+"k1=0.57; //from table\n",
+"k2=0.89;\n",
+"Qunet2=qp*Ap*k1+k*k2;\n",
+"disp(Qunet2/1000,'allowed load in kN');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 12.6: 6.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 12.6\n",
+"clc; funcprot(0);\n",
+"Qws=1005-250;\n",
+"Qwp=250;\n",
+"epsilon=0.65;\n",
+"L=11;\n",
+"Ds=1.5;\n",
+"Es=14000;\n",
+"Ap=%pi/4*1.5^2;\n",
+"Ep=21e6;\n",
+"Cp=0.04;\n",
+"Db=1.5;\n",
+"mus=0.3;\n",
+"p=%pi*1.5;\n",
+"Nc=9;\n",
+"qp=105*Nc;\n",
+"se1=(Qwp+epsilon*Qws)*L/(Ap*Ep);\n",
+"se2=Qwp*Cp/(Db*qp);\n",
+"Iws=2+0.35*sqrt(L/Ds);\n",
+"se3=Qws/p/L*Ds/Es*(1-mus^2)*Iws;\n",
+"se=se1+se2+se3;\n",
+"disp(se*1000,'net settlement in mm');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 12.7: 7.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 12.7\n",
+"clc; funcprot(0);\n",
+"Ds=1;\n",
+"Ep=22e6;\n",
+"Ri=1;\n",
+"cu=100;\n",
+"Ip=%pi*Ds^4/64;\n",
+"Qc=7.34*Ds^2*Ep*Ri*(cu/Ep/Ri)^0.6;\n",
+"disp(Qc,'bearing force in kN');\n",
+"Mc=3.86*Ds^3*Ep*Ri*(cu/Ep/Ri)^0.6;\n",
+"disp(Mc,'bearing moment on kNm')\n",
+"//from figure\n",
+"xoQM=0.0046*1;\n",
+"xoMQ=0.0041*1;\n",
+"xo=0.5*(xoQM+xoMQ);\n",
+"disp(xo*1000,'net ground line deflection in mm');\n",
+"//partb\n",
+"Ip=0.049;\n",
+"Qg=150;\n",
+"Mg=200;\n",
+"deff('y=f(T)','y=338e-6*T^3+300.6e-6*T^2-0.00435');\n",
+"[x]=fsolve(2,f);\n",
+"T=x;\n",
+"k=[0, 0.4, 0.6, 0.8, 1, 1.1, 1.25];//z/T\n",
+"Am=[0, 0.36, 0.52, 0.63, 0.75, 0.765, 0.75];\n",
+"Bm=[1, 0.98, 0.95, 0.9, 0.845, 0.8, 0.73];\n",
+"printf('z/T\t Am\t Bm\t Mz(kN-m)\n');\n",
+"for i=1:7\n",
+"    Mz(i)=Am(i)*Qg*T+Bm(i)*Mg;\n",
+"    printf('%.2f\t %.2f\t %.2f\t %.2f \n',k(i),Am(i),Bm(i),Mz(i));\n",
+"end\n",
+"disp(1*T,'depth in m');\n",
+"//partc\n",
+"Mmax=400;\n",
+"sigma=Mmax*Ds/2/Ip;\n",
+"disp(sigma,'tensile stress in kN/m^2');\n",
+"//partd\n",
+"//from figure\n",
+"k=0.85;\n",
+"L=k*1;\n",
+"disp(L,'length in m');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 12.8: 8.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 12.8\n",
+"clc; funcprot(0);\n",
+"qu=3000;\n",
+"Ds=3*12;\n",
+"L=15*12;\n",
+"FS=3;\n",
+"Ecore=0.36e6;\n",
+"f=min(2.5*qu^0.5,0.15*qu);\n",
+"Qu=%pi*Ds*L*f/1000;\n",
+"Emass=Ecore*(0.266*80-1.66);\n",
+"Ec=17.9*Emass;\n",
+"Ac=%pi/4*Ds^2;\n",
+"If=0.35;\n",
+"se=Qu*L/Ac/Ec+Qu*If/Ds/Emass;\n",
+"Qall=Qu/FS;\n",
+"disp(Qall,'allowed load in kN');\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_Foundation_Engineering_by_B_M_Das/13-Foundations_on_Difficult_Soils.ipynb b/Principles_Of_Foundation_Engineering_by_B_M_Das/13-Foundations_on_Difficult_Soils.ipynb
new file mode 100644
index 0000000..604946f
--- /dev/null
+++ b/Principles_Of_Foundation_Engineering_by_B_M_Das/13-Foundations_on_Difficult_Soils.ipynb
@@ -0,0 +1,134 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 13: Foundations on Difficult Soils"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 13.1: 1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 13.1\n",
+"clc; funcprot(0);\n",
+"Sw=1;\n",
+"Z=2;\n",
+"deltaSf=0.0033*Z*Sw*1000;\n",
+"disp(deltaSf,'free surface swell in mm');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 13.2: 2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 13.2\n",
+"clc; funcprot(0);\n",
+"//from figure\n",
+"deltaS=1/100*1/2*(0.55+0+0.55+1.2+1.2+2+2+3);\n",
+"disp(deltaS*1000,'total swell in mm');\n",
+"//partb\n",
+"D=[5.2, 4.2, 3.2, 2.2, 1.2];\n",
+"deltaS=[0, 0.00275, 0.0115, 0.0275, 0.0525];\n",
+"printf('depth(m)\t  total swell (m) \n');\n",
+"for i=1:5\n",
+"    printf('%.2f\t            %f \n', D(i),deltaS(i));\n",
+"end\n",
+"plot(deltaS*1000,D);\n",
+"xtitle('depth vs total swell','total swell (m)','depth (m)');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 13.3: 3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 13.3\n",
+"clc; funcprot(0);\n",
+"phi=12*%pi/180;\n",
+"Ds=0.8;\n",
+"Z=5;\n",
+"sigmaT=450;\n",
+"U=%pi*Ds*Z*sigmaT*tan(phi);\n",
+"deff('y=f(D)','y=1202-450*6.14/1.25*3.14/4*(D^2-0.8^2)');\n",
+"[x]=fsolve(1,f);\n",
+"Db=x;\n",
+"disp(Db,'diameter of bell in m');\n",
+"//partb\n",
+"D=600;\n",
+"cu=450;\n",
+"Nc=6.14;\n",
+"FS=cu*Nc*%pi/4*(Db^2-Ds^2)/(U-D);\n",
+"if FS>2 then\n",
+"    disp('the structure is compatible with safety measures');\n",
+"end\n",
+"//check bearing capacity\n",
+"L=D+300;//dead+live load\n",
+"Dp=L/%pi*4/Db^2;//downward pressure\n",
+"FS=2763/Dp;\n",
+"if FS>3 then\n",
+"    disp('the structure is safe in bearing ');\n",
+"end\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_Foundation_Engineering_by_B_M_Das/14-Soil_Improvement_and_Ground_Modification.ipynb b/Principles_Of_Foundation_Engineering_by_B_M_Das/14-Soil_Improvement_and_Ground_Modification.ipynb
new file mode 100644
index 0000000..3b93d53
--- /dev/null
+++ b/Principles_Of_Foundation_Engineering_by_B_M_Das/14-Soil_Improvement_and_Ground_Modification.ipynb
@@ -0,0 +1,132 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 14: Soil Improvement and Ground Modification"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 14.1: 1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 14.1\n",
+"clc; funcprot(0);\n",
+"Cc=0.28;\n",
+"Hc=6;\n",
+"eo=0.9;\n",
+"Cv=0.36;\n",
+"H=3;\n",
+"t=9;\n",
+"sigmao=210;\n",
+"sigmap=115; //deltasigmap\n",
+"Sc=Cc*Hc/(1+eo)*log10((sigmao+sigmap)/sigmao);\n",
+"disp(Sc*1000,'primary consolidation in mm');\n",
+"Tv=Cv*t/H^2;\n",
+"//from table\n",
+"k=1.8; //constant\n",
+"sf=k*sigmap;\n",
+"disp(sf,'deltasigmaf in kN/m^2');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 14.2: 2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 14.2\n",
+"clc; funcprot(0);\n",
+"Tv=0.36;\n",
+"sigmap=115;\n",
+"Uv=sqrt(4*Tv/%pi)*100;\n",
+"disp(Uv,'Uv in %');\n",
+"//from table \n",
+"k=0.12; //constant\n",
+"sf=k*sigmap;\n",
+"disp(sf,'deltasigmaf in kN/m^2');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 14.3: 3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 14.3\n",
+"clc(); funcprot(0);\n",
+"Cc=0.31;\n",
+"Hc=15;//ft\n",
+"eo=1.1;\n",
+"n=10;\n",
+"Uv=0.09;\n",
+"sigmao=1000;\n",
+"deltasigma=2000; // deltasigmap+deltasigmaf\n",
+"Sc=Cc*Hc/(1+eo)*log10((sigmao+deltasigma)/sigmao);\n",
+"disp(Sc,'primary consolidation in ft');\n",
+"m=n^2/(n^2-1)*log(n)-(3*n^2-1)/4/n^2;\n",
+"A=2/m;\n",
+"Ur=(0.096-1/A*(1-exp(-A*0.096)))/0.192;\n",
+"Uvf=1-(1-Ur)*(1-Uv);\n",
+"Sc30=Sc*Uvf*12; //settlement after 30 days\n",
+"disp(Sc30,'settlement after 30 days in inch');\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_Foundation_Engineering_by_B_M_Das/2-Natural_Soil_Deposits_and_Subsoil_Exploration.ipynb b/Principles_Of_Foundation_Engineering_by_B_M_Das/2-Natural_Soil_Deposits_and_Subsoil_Exploration.ipynb
new file mode 100644
index 0000000..1227a46
--- /dev/null
+++ b/Principles_Of_Foundation_Engineering_by_B_M_Das/2-Natural_Soil_Deposits_and_Subsoil_Exploration.ipynb
@@ -0,0 +1,78 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 2: Natural Soil Deposits and Subsoil Exploration"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 2.1: 1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 2.1\n",
+"clc; funcprot(0);\n",
+"Distance=[2.5,5,7.5,10,15,20,25,30,35,40,50];\n",
+"Time=10^-3*[11.2,23.3,33.5,42.4,50.9,57.2,64.4,68.6,71.1,72.1,75.5]\n",
+"//part1\n",
+"distance=5.25;\n",
+"time=23e-3;\n",
+"v1=distance/time;\n",
+"disp(v1,'speed in m/s');\n",
+"//part2\n",
+"distance=11;\n",
+"time=13.5e-3;\n",
+"v2=distance/time;\n",
+"disp(v2,'speed in m/s');\n",
+"//part3\n",
+"distance=14.75;\n",
+"time=3.5e-3;\n",
+"v3=distance/time;\n",
+"disp(v3,'speed in m/s');\n",
+"plot(Distance,Time);\n",
+"xtitle('distance vs time','Distance in m','time in s');\n",
+"//part4\n",
+"xc=10.4;\n",
+"Ta=65e-3;\n",
+"Z1=1/2*sqrt((v2-v1)/(v2+v1))*xc;\n",
+"disp(Z1,'thickness of layer 1 in m');\n",
+"Z2=1/2*(Ta-2*Z1*sqrt(v3^2-v1^2)/v3/v1)*v3*v2/sqrt(v3^2-v2^2);\n",
+"disp(Z2,'thickness of layer 2 in m');\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_Foundation_Engineering_by_B_M_Das/3-Shallow_Foundations.ipynb b/Principles_Of_Foundation_Engineering_by_B_M_Das/3-Shallow_Foundations.ipynb
new file mode 100644
index 0000000..0032d1a
--- /dev/null
+++ b/Principles_Of_Foundation_Engineering_by_B_M_Das/3-Shallow_Foundations.ipynb
@@ -0,0 +1,330 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 3: Shallow Foundations"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.1: 1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 3.1\n",
+"clc; funcprot(0);\n",
+"Nc=17.69;\n",
+"Nq=7.44;\n",
+"Ny=3.64;\n",
+"q=3*115;\n",
+"Gamma=115; //lb/ft^3\n",
+"c=320;\n",
+"B=5;//ft\n",
+"FS=4;//factor of safety\n",
+"qu=1.3*c*Nc+q*Nq+0.4*Gamma*B*Ny\n",
+"qall=qu/FS; //q allowed\n",
+"Q=qall*B^2;\n",
+"disp(Q,'allowable gross load in lb');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.2: 2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 3.2\n",
+"clc; funcprot(0);\n",
+"Gamma=105;//lb/ft^3\n",
+"Gammasat=118;//lb/ft^3\n",
+"FS=3;\n",
+"pa=2014.125;//lb/ft^2\n",
+"Depth=[5,10,15,20,25];\n",
+"N60=[4,6,6,10,5];\n",
+"sigmao=[0,0,0,0,0];\n",
+"phi=[0,0,0,0,0]\n",
+"Gammaw=62.4;\n",
+"s=0;\n",
+"printf('depth (ft)\tN60\t  \tstress(lb/ft^2)\t phi(degrees)\n')\n",
+"for i=1:5\n",
+"    sigmao(i)=2*Gamma+(Depth(i)-2)*(Gammasat-Gammaw);\n",
+"    phi(i)=sqrt(20*N60(i)*sqrt(pa/sigmao(i)))+20;\n",
+"    printf(' %.2f\t       %.2f\t\t   %.2f  \t%.2f \n',Depth(i),N60(i),sigmao(i),phi(i));\n",
+"avgphi=phi(i)/5+s;\n",
+"s=avgphi;\n",
+"end\n",
+"disp(round(avgphi),'average friction angle in degrees');\n",
+"//using graph get the values of other terms in terms of B and solve for B\n",
+"deff('y=f(x)','y=-150000/x^2+5263.9+5527.1/x+228.3*x');\n",
+"[x, v, info ]=fsolve(4,f);\n",
+"disp(x,'the width in ft');\n",
+"\n",
+"\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.3: 3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 3.3\n",
+"clc; funcprot(0);\n",
+"phi=25; //degrees\n",
+"Es=620; //kN/m^2\n",
+"Gamma=18;//kN/m^2\n",
+"Df=0.6;\n",
+"B=0.6;\n",
+"L=1.2;\n",
+"Fqc=0.347;\n",
+"Nq=10.66;\n",
+"Nc=20.72;\n",
+"Ngamma=10.88;\n",
+"mu=0.3;\n",
+"Fyd=1;\n",
+"c=48;//kN/m^2\n",
+"q=Gamma*(Df+B/2);\n",
+"Ir=Es/(2*(1+mu)*(c+q*tan(phi*%pi/180)));\n",
+"disp(Ir,'value of Ir');\n",
+"Fcc=Fqc-(1-Fqc)/(Nq*tan(phi*%pi/180));\n",
+"Fcs=1+Nq/Nc*B/L;\n",
+"Fqs=1+B/L*tan(phi*%pi/180);\n",
+"Fys=1-0.4*B/L;\n",
+"Fcd=1+0.4*Df/B;\n",
+"Fqd=1+2*tan(phi*%pi/180)*(1-sin(phi*%pi/180))^2*Df/B;\n",
+"q1=0.6*18;\n",
+"Fyc=Fqc;\n",
+"qu=c*Nc*Fcs*Fcd*Fcc+q1*Nq*Fqs*Fqd*Fqc+1/2*Gamma*Ngamma*Fys*Fyd*Fyc;\n",
+"disp(qu,'ultimate bearing capacity in kN/m^2');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.4: 4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 3.4\n",
+"clc; funcprot(0);\n",
+"q=110*4; //lb/ft^2\n",
+"Nq=33.3;\n",
+"phi=35;\n",
+"Df=4;\n",
+"B=6;\n",
+"Gamma=110;\n",
+"Ngamma=48.03;\n",
+"B1=6-2*0.5;\n",
+"Fqi=1;\n",
+"Fyi=1;\n",
+"Fyd=1;\n",
+"Fqs=1;\n",
+"Fys=1;\n",
+"Fqd=1+2*tan(phi*%pi/180)*(1-sin(phi*%pi/180))^2*Df/B;\n",
+"qu=q*Nq*Fqs*Fqd*Fqi+1/2*B1*Gamma*Ngamma*Fys*Fyd*Fyi;\n",
+"Qult=B1*1*qu;\n",
+"disp(Qult,'ultimate bearing capacity in lb/ft');\n",
+"disp(Qult/2000,'ultimate bearing capacity in ton/ft');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.5: 5.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 3.5\n",
+"clc; funcprot(0);\n",
+"e=0.5;\n",
+"B=6;\n",
+"k=e/B;\n",
+"Gamma=110;\n",
+"q=440;\n",
+"disp('get the values of Nqe and Nye from the figure from the value of e/B');\n",
+"Nye=26.8;\n",
+"Nqe=33.4;\n",
+"Qult=B*1*(q*Nqe+1/2*Gamma*B*Nye);\n",
+"disp(Qult,'ultimate bearing capacity in lb/ft');\n",
+"disp(Qult/2000,'ultimate bearing capacity in ton/ft');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.6: 6.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 3.6\n",
+"clc; funcprot(0);\n",
+"Df=0.7;\n",
+"//from table\n",
+"Nq=18.4;\n",
+"Ny=22.4;\n",
+"q=12.6;\n",
+"phi=30; //angle\n",
+"L=1.5;\n",
+"Fyd=1;\n",
+"Gamma=18;\n",
+"L1=0.85*1.5;\n",
+"L2=0.21*1.5;\n",
+"B=1.5;\n",
+"A=1/2*(L1+L2)*B;\n",
+"B1=A/L1; //B'\n",
+"Fqs=1+B1/L1*tan(phi*%pi/180);\n",
+"Fys=1-0.4*B1/L1;\n",
+"Fqd=1+2*tan(phi*%pi/180)*(1-sin(phi*%pi/180))^2*Df/B;\n",
+"Qult=A*(q*Nq*Fqs*Fqd+1/2*Gamma*B1*Ny*Fys*Fyd);\n",
+"disp(Qult,'ultimate load in kN');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.7: 7.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 3.7\n",
+"clc; funcprot(0);\n",
+"e=0.15;\n",
+"B=1.5;\n",
+"Fqs=1;\n",
+"L=1.5;\n",
+"Gamma=18;\n",
+"q=0.7*18;\n",
+"//from table\n",
+"Nqe=18.4;\n",
+"Nye=11.58;\n",
+"Fys=1+(2*e/B-0.68)*(B/L)+(0.43-3/2*e/B)*(B/L)^2;\n",
+"Qult=B*L*(q*Nqe*Fqs+1/2*Gamma*Nye*Fys);\n",
+"disp(Qult,'ultimate load in kN');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 3.8: 8.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 3.8\n",
+"clc; funcprot(0);\n",
+"q=16;\n",
+"Nqei=14.2;\n",
+"Gamma=16\n",
+"B=1.5;\n",
+"Nyet=20;\n",
+"Qult=B*(Nqei*q+1/2*Gamma*B*Nyet);\n",
+"disp(Qult,'ultimate load in kN/m');\n",
+"\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_Foundation_Engineering_by_B_M_Das/4-Ultimate_Bearing_Capacity_of_Shallow_Foundations.ipynb b/Principles_Of_Foundation_Engineering_by_B_M_Das/4-Ultimate_Bearing_Capacity_of_Shallow_Foundations.ipynb
new file mode 100644
index 0000000..64cce6e
--- /dev/null
+++ b/Principles_Of_Foundation_Engineering_by_B_M_Das/4-Ultimate_Bearing_Capacity_of_Shallow_Foundations.ipynb
@@ -0,0 +1,295 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 4: Ultimate Bearing Capacity of Shallow Foundations"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 4.1: 1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 4.1\n",
+"clc; funcprot(0);\n",
+"FS=4;\n",
+"q=110*2;\n",
+"Nq=90;\n",
+"Ny=50;\n",
+"Gamma=110;\n",
+"m1=0.34;\n",
+"B=2.5;\n",
+"L=2.5;\n",
+"H=1.5;\n",
+"phi=35;\n",
+"m2=0.45;\n",
+"Fqs=1-0.34*B/L;\n",
+"Fys=1-0.45*B/L;\n",
+"qu=q*Nq*Fqs+1/2*Gamma*Ny*Fys*B;\n",
+"Qall=qu*B^2/FS;\n",
+"disp(Qall,'bearing load in lb');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 4.2: 2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 4.2\n",
+"clc; funcprot(0);\n",
+"FS=3;\n",
+"cu=72;\n",
+"q=18;\n",
+"B=1;\n",
+"H=0.25;\n",
+"qu=5.14*(1+(0.5*B/H-0.707)/5.14)*cu+q;\n",
+"qall=qu/FS;\n",
+"disp(qall,'bearing capacity of soil in kN/m^2');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 4.3: 3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 4.3\n",
+"clc; funcprot(0);\n",
+"k=0;//B/L;\n",
+"c2=30;\n",
+"Gamma=17.5;\n",
+"H=1.5;\n",
+"Df=1.2;\n",
+"H=1.5;\n",
+"B=2.0;\n",
+"Ks=2.5;\n",
+"phi=40;\n",
+"pi=%pi;\n",
+"qu=(1+0.2*k)*5.14*c2+(1+k)*Gamma*H^2*(1+2*Df/H)*Ks*tan(phi*pi/180)/B+Gamma*H;\n",
+"Qu=qu*B;\n",
+"disp(Qu,'bearing capacity in kN/m');\n",
+"disp('there is slight variation due to rounding off error')\n",
+"//soil 2\n",
+"Ny=109.4;\n",
+"Nq=64.2;\n",
+"Fqs=1;\n",
+"Fys=1;\n",
+"qt=Gamma*Df*Nq*Fqs+1/2*Gamma*Ny*Fys*B;\n",
+"disp(qt,'bearing capacity in kN/m^2');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 4.4: 4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 4.4\n",
+"clc; funcprot(0);\n",
+"B=1;\n",
+"L=1.5;\n",
+"c2=48;\n",
+"ca=108;\n",
+"D=1;\n",
+"H=1;\n",
+"Gamma=16.8;\n",
+"FS=4;\n",
+"qu=(1+0.2*B/L)*5.14*c2+(1+B/L)*2*ca*H/B+Gamma*D;\n",
+"c1=120;\n",
+"gamma1=16.8;\n",
+"Df=1;\n",
+"qt=(1+0.2*B/L)*5.14*c1+gamma1*Df;\n",
+"disp(qt,'qt in kN/m^2');\n",
+"disp('no need to calculate qt since it is not useful for calculation')\n",
+"disp(qu/FS,'allowable shear stress in kN/m^2');\n",
+"disp(qu/FS*1*1.5,'allowable load in kN');\n",
+"clear()"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 4.5: 5.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 4.5\n",
+"clc; funcprot(0);\n",
+"c=50;\n",
+"//from table\n",
+"Ncq=6.3;\n",
+"FS=4;\n",
+"qu=c*Ncq;\n",
+"qall=qu/4;\n",
+"disp(qall,'allowed shear stress in kN/m^2');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 4.6: 6.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 4.6\n",
+"clc; funcprot(0);\n",
+"Gamma=16.8;\n",
+"B=1.5;\n",
+"//from table\n",
+"Nyq=120;\n",
+"qu=1/2*Gamma*B*Nyq;\n",
+"disp(qu,' shear stress in kN/m^2');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 4.7: 7.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 4.7\n",
+"clc; funcprot(0);\n",
+"phi=35;\n",
+"Df=1.5;\n",
+"B=1.5;\n",
+"Gamma=17.4;\n",
+"A=%pi/4*Df^2;\n",
+"m=0.25;\n",
+"Ku=0.936;\n",
+"Fq=1+2*(1+m*Df/B)*Df/B*Ku*tan(phi*%pi/180);\n",
+"Qu=Fq*Gamma*A*Df;\n",
+"disp(Qu,'bearing capacity in kN');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 4.8: 8.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 4.8\n",
+"clc; funcprot(0);\n",
+"cu=52;\n",
+"B=1.5;\n",
+"L=3;\n",
+"k=0.107*cu+2.5;\n",
+"disp(k,'Df/B of square');\n",
+"A=L*B;\n",
+"Beta=0.2;\n",
+"Gamma=18.9;\n",
+"Df=1.8;\n",
+"Qu=A*(Beta*(7.56+1.44*B/L)*cu+Gamma*Df);\n",
+"disp(Qu,'ultimate shear force in kN');\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_Foundation_Engineering_by_B_M_Das/5-Shallow_Foundations_Allowable_Bearing_Capacity_and_Settlement.ipynb b/Principles_Of_Foundation_Engineering_by_B_M_Das/5-Shallow_Foundations_Allowable_Bearing_Capacity_and_Settlement.ipynb
new file mode 100644
index 0000000..7f78a82
--- /dev/null
+++ b/Principles_Of_Foundation_Engineering_by_B_M_Das/5-Shallow_Foundations_Allowable_Bearing_Capacity_and_Settlement.ipynb
@@ -0,0 +1,433 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 5: Shallow Foundations Allowable Bearing Capacity and Settlement"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.10: 10.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 5.10\n",
+"clc; funcprot(0);\n",
+"Cc=0.32;\n",
+"Hc=2.5;\n",
+"eo=0.8;\n",
+"sigmao=2.5*16.5+0.5*(17.5-9.81)+1.25*(16-9.81);\n",
+"m1=[2, 2, 2];\n",
+"z=[2, 3.25, 4.5];\n",
+"n1=[4, 6.5, 9];\n",
+"Ic=[0.19, 0.085, 0.045];\n",
+"Dsigma=[28.5, 12.75, 6.75];//deltasigma\n",
+"printf('m1\t z(m)\t n1\t Ic\t Dsigma \n');\n",
+"for i=1:3\n",
+"    printf('%.2f\t %.2f\t %.2f\t %.2f\t %.2f\n',m1(i),z(i),n1(i),Ic(i),Dsigma(i));\n",
+"end\n",
+"Dsigmaav=1/6*(Dsigma(1)+4*Dsigma(2)+Dsigma(3));\n",
+"Sc=Cc*Hc/(1+eo)*log10((sigmao+Dsigmaav)/sigmao);\n",
+"disp(Sc*1000,'settlement in mm')\n",
+"//partb\n",
+"B=1;\n",
+"L=2;\n",
+"z=0.5+1.5;\n",
+"B=B+z;\n",
+"L=L+z;\n",
+"A=0.6;\n",
+"//from table\n",
+"kcr=0.36;\n",
+"Sep=kcr*Sc;\n",
+"disp(Sep*1000,'settlement in mm');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.11: 11.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 5.11\n",
+"clc; funcprot(0);\n",
+"N60=(3+7+12+12+16)/5;\n",
+"B=[2, 2.25, 2.3];\n",
+"Fd=[1.248, 1.22, 1.215];\n",
+"Qoac=102000*9.81/1000;//actual Qo\n",
+"Se=25;\n",
+"printf('B(m)\t Fd\t qnet(kN/m^2)\t Qo \n');\n",
+"for i=1:3\n",
+"    qnet(i)=10/0.08*(B(i)+0.3)^2/(B(i))^2*Fd(i)*Se/25;\n",
+"    Qo(i)=qnet(i)*B(i)^2;\n",
+"    printf('%.2f\t %.2f\t %.2f\t\t %.2f\n',B(i), Fd(i),qnet(i),Qo(i))\n",
+"end\n",
+"disp(Qoac,'value of Qo in kN');\n",
+"printf('since Qo is 1000 kN thus B is equal to 2.3m from the table');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.1: 1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 5.1\n",
+"clc; funcprot(0);\n",
+"//first solution\n",
+"B1=2.5;\n",
+"B2=B1;\n",
+"z=12.5;\n",
+"L1=5;\n",
+"L2=L1;\n",
+"m=B1/z;\n",
+"n=B2/z;\n",
+"//from table of tthe values using m,n\n",
+"q=2000;\n",
+"I=0.0328;\n",
+"deltasigma=q*4*I;\n",
+"disp(deltasigma,'change in pressure in lb/ft^2');\n",
+"//second solution\n",
+"Ic=0.131;//from table\n",
+"deltasigma=q*Ic;\n",
+"disp(deltasigma,'change in pressure in lb/ft^2');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.2: 2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 5.2\n",
+"clc; funcprot(0);\n",
+"qo=100;\n",
+"H1=3;\n",
+"H2=5;\n",
+"//from table\n",
+"IaH2=0.126;\n",
+"IaH1=0.175;\n",
+"deltasigma=qo*((H2*IaH2-H1*IaH1)/(H2-H1));\n",
+"disp(deltasigma,'change in pressure in kN/m^2');\n",
+"TS=4*deltasigma;\n",
+"disp(TS,'total change in pressure in kN/m^2');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.3: 3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 5.3\n",
+"clc; funcprot(0);\n",
+"H=7;\n",
+"Gamma=17.5;\n",
+"q0=Gamma*H\n",
+"disp(q0,'pressure change in kN/m^2');\n",
+"//part2\n",
+"//from figure\n",
+"Ileftside=0.445;\n",
+"Irightside=0.445;\n",
+"deltasigma=q0*(Ileftside+Irightside);\n",
+"disp(deltasigma,'change in stress in kN/m^2');\n",
+"//partc\n",
+"//from figure 5.11\n",
+"I=0.24;//I'\n",
+"Dsigma1=43.75*I;//deltasigma1\n",
+"I2=0.495;//I'\n",
+"Dsigma2=I2*q0;//deltasigma2\n",
+"I3=0.335;//I'\n",
+"Dsigma3=I3*78.75;//deltasigma3\n",
+"Dsigma=Dsigma1+Dsigma2-Dsigma3;\n",
+"disp(Dsigma,'total stress increase in A in kN/m^2');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.4: 4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 5.4\n",
+"clc; funcprot(0);\n",
+"zbar=5;\n",
+"mus=0.3;\n",
+"F1=0.641;\n",
+"F2=0.031;\n",
+"z1=2;\n",
+"z2=1;\n",
+"z3=2;\n",
+"Es1=10000;\n",
+"Es2=8000;\n",
+"Es3=12000;\n",
+"qo=150;\n",
+"//from table\n",
+"If=0.709;\n",
+"Es=(Es1*z1+Es2*z2+Es3*z3)/zbar;\n",
+"disp(Es,'modulus of elasticity in kN/m^2');\n",
+"Is=F1+(2-mus)/(1-mus)*F2;\n",
+"Sc=qo*(1/Es-mus^2/Es)*Is*If*2;\n",
+"Scrigid=0.93*Sc;\n",
+"disp(Scrigid*1000,'settelement in mm');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.5: 5.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 5.5\n",
+"clc; funcprot(0);\n",
+"B=5;\n",
+"L=10;\n",
+"Ef=2.3e6;\n",
+"Eo=1400;\n",
+"k=25;\n",
+"t=1.0;\n",
+"mus=0.3;\n",
+"Df=5;\n",
+"qo=5000;\n",
+"Ig=0.69;\n",
+"Be=sqrt(4*B*L/%pi);\n",
+"If=%pi/4+1/(4.6+10*(Ef/(Eo+2*Be/2*k))*(2*t/Be)^3);\n",
+"Ie=1-1/(3.5*exp(1.22*mus-0.4)*(Be/Df+1.6));\n",
+"Se=qo*Be*Ig*If*Ie/Eo*(1-mus^2)/144;\n",
+"disp(Se*12,'settlement in inches');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.6: 6.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 5.6\n",
+"clc; funcprot(0);\n",
+"q=3.06;\n",
+"qbar=25;\n",
+"C1=1-0.5*(q/(qbar-q));\n",
+"Sum=0;\n",
+"C2=1+0.2*log10(10/0.1);\n",
+"L=[1, 2, 3, 4, 5];\n",
+"Dz=[48, 48, 96, 48, 144];\n",
+"Es=[750, 1250, 1250, 1000, 2000];\n",
+"z=[24, 72, 144, 216, 312];\n",
+"Iz=[0.275, 0.425, 0.417, 0.292, 0.125];\n",
+"printf('Layer No.\t deltaz (in)\t Es(lb/in^2)\t z to the middle of the layer (in)  Iz at the middle of the layer  Iz/delta(z) \n');\n",
+"for i=1:5\n",
+"    k(i)=Iz(i)/Es(i)*Dz(i);\n",
+"    printf('%.2f\t \t %.2f\t\t   %.2f\t\t %.2f \t\t\t\t\t   %.2f\t\t    %.4f \n ',L(i),Dz(i),Es(i),z(i),Iz(i),k(i));\n",
+"    Sum=Sum+k(i);\n",
+"end\n",
+"Se=C1*C2*(qbar-q)*Sum;\n",
+"disp(Se,'settlement in inches');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.7: 7.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 5.8\n",
+"clc; funcprot(0);\n",
+"Df=1;\n",
+"B=1.75;\n",
+"qnet=120;\n",
+"N60=10;\n",
+"Fd=1+0.33*Df/B;\n",
+"Se=2*qnet/N60/Fd*(B/(B+0.3))^2;\n",
+"disp(Se,'settlement in mm');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.8: 8.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 5.8\n",
+"clc; funcprot(0);\n",
+"Df=1;\n",
+"B=1.75;\n",
+"qnet=120;\n",
+"N60=10;\n",
+"Fd=1+0.33*Df/B;\n",
+"Se=2*qnet/N60/Fd*(B/(B+0.3))^2;\n",
+"disp(Se,'settlement in mm');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 5.9: 9.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 5.9\n",
+"clc; funcprot(0);\n",
+"Ny=23.76;\n",
+"Nq=16.51;\n",
+"q=3*110;\n",
+"Gamma=110;\n",
+"B=4;\n",
+"Nqe=0.63*Nq;\n",
+"Nye=0.4*Ny;\n",
+"que=q*Nqe+1/2*Gamma*B*Nye;\n",
+"disp(que,'bearing capacity in lb/ft^2');\n",
+"//part 2\n",
+"V=0.4;\n",
+"A=0.32;\n",
+"g=9.81;\n",
+"kh=0.26;\n",
+"k=0.92;//tan(alphae)\n",
+"Seq=0.174*k*V^2/A/g*kh^-4/A^-4;\n",
+"disp(Seq,'settelement in m');\n",
+"disp(Seq*39.57,'settlement in inches')\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_Foundation_Engineering_by_B_M_Das/6-Mat_Foundations.ipynb b/Principles_Of_Foundation_Engineering_by_B_M_Das/6-Mat_Foundations.ipynb
new file mode 100644
index 0000000..e6145ce
--- /dev/null
+++ b/Principles_Of_Foundation_Engineering_by_B_M_Das/6-Mat_Foundations.ipynb
@@ -0,0 +1,259 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 6: Mat Foundations"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.1: 1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 6.1\n",
+"clc; funcprot(0);\n",
+"B=30;\n",
+"L=45;\n",
+"Df=6.5;\n",
+"cu=1950;\n",
+"qunet=5.14*cu*(1+0.195*B/L)*(1+0.4*Df/B);\n",
+"disp(qunet,'allowed force in lb/ft^2');\n",
+"\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.2: 2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 6.2\n",
+"clc; funcprot(0);\n",
+"N60=10;\n",
+"Df=2;\n",
+"B=10;\n",
+"Se=25;\n",
+"qnetall=N60/0.08*(1+0.33*Df/B)*Se/25;\n",
+"disp(qnetall,'allowed pressure in kN/m^2');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.3: 3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 6.3\n",
+"clc; funcprot(0);\n",
+"cu=2800;\n",
+"B=60;\n",
+"L=100;\n",
+"Df=5;\n",
+"B=60;\n",
+"Gamma=120;\n",
+"A=60*100;\n",
+"Q=25e6;\n",
+"FS=5.14*cu*(1+0.195*B/L)*(1+0.4*Df/B)/(Q/A-Gamma*Df);\n",
+"disp(FS,'factor of safety');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.4: 4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 6.4\n",
+"clc; funcprot(0);\n",
+"Cc=0.28;\n",
+"Hc=18*12;\n",
+"e0=0.9;\n",
+"sigmao=11*100+40*(121.5-64)+18/2*(118-62.4);\n",
+"H2=5+40+18;\n",
+"H1=5+40;\n",
+"qo=3567;\n",
+"//from table\n",
+"IaH2=0.21;\n",
+"IaH1=0.225;\n",
+"Dsigma=qo*((H2*IaH2-H1*IaH1)/(H2-H1))*4;\n",
+"Scp=Cc*Hc/(1+e0)*log10(sigmao/sigmao+Dsigma/sigmao);\n",
+"disp(Scp,'settlement in inches');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.5: 5.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 6.5\n",
+"clc; funcprot(0);\n",
+"P=['A','B','C','D','E','F','G','H','I','J','K','L','M','N'];//point\n",
+"k=1.2*ones(1,14);//Q/A\n",
+"x=[-38,-24, -12, 0, 12, 24, 38, 38, 24, 12, 0, -12, -24, -38];\n",
+"x1=0.0017*x;\n",
+"y=[48,48,48,48,48,48,48, -48, -48, -48, -48, -48, -48, -48];\n",
+"y1=-0.0011*y;\n",
+"printf('point\t Q\A (kip/ft^2)\t x(ft)\t 0.0017x(ft)\t\t y(ft)\t 0.0011y(ft)\t q(kip/ft^2)\n')\n",
+"for i=1:14\n",
+"    q(i)=1.2+x1(i)+y1(i);\n",
+"    printf('%s\t %.2f\t\t %.2f\t\t%.2f\t\t %.2f\t %.2f\t\t %.2f\n ',P(i),k(i),x(i),x1(i),y(i),y1(i),q(i))\n",
+"end\n",
+"printf('the soil pressure at all point is less than the given qallnet=1.5 kip/ft^2');\n",
+"\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 6.6: 6.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 6.6\n",
+"clc; funcprot(0);\n",
+"//solving for d\n",
+"deff('y=f(d)','y=(96+2*d)*d-2615.1');\n",
+"[x]=fsolve(19,f);\n",
+"d1=x;\n",
+"deff('y=f(d)','y=(96+4*d)*d-6046.4');\n",
+"[x]=fsolve(28,f);\n",
+"d2=x;\n",
+"d=max(d2,d1);\n",
+"d=round(d)\n",
+"//now coming to design part\n",
+"h=d+3+1;\n",
+"disp(h,'total slab thickness in inches');\n",
+"qa=1.082;\n",
+"qb=1.106;\n",
+"qm=1.212;\n",
+"qn=1.188;\n",
+"q1A=qa/2+qb/2;\n",
+"disp(q1A,'force in strip ABMN in kip/ft^2');\n",
+"q2A=qm/2+qn/2;\n",
+"disp(q2A,'force in strip ABMN in kip/ft^2');\n",
+"q1=1.106/3+1.127/3+1.147/3;\n",
+"disp(q1,'force in strip BCDKLM in kip/ft^2');\n",
+"q2=1.253/3+1.233/3+1.212/3;\n",
+"disp(q2,'force in strip BCDKLM in kip/ft^2');\n",
+"q1=1.147/3+1.167/3+1.188/3;\n",
+"disp(q1,'force in strip DEFIJK in kip/ft^2');\n",
+"q2=1.294/3+1.273/3+1.253/3;\n",
+"disp(q2,'force in strip DEFIJK in kip/ft^2');\n",
+"q1=1.188/2+1.212/2;\n",
+"disp(q1,'force in strip FGHI in kip/ft^2');\n",
+"q2=1.318/2+1.294/2;\n",
+"disp(q2,'force in strip FGHI in kip/ft^2');\n",
+"//checking for force\n",
+"//net soil reaction <load \n",
+"netforce=1/2*(1.094+1.2)*14*96+1/2*(1.127+1.233)*24*96+1/2*(1.167+1.273)*24*96+1/2*(1.2+1.306)*14*96;\n",
+"if netforce<8761 then\n",
+"    disp('forces generated is OK')\n",
+"end\n",
+"//checking for reinforcement requirement\n",
+"Q1=1.4*180+1.7*120;\n",
+"Q2=1.4*360+1.7*200;\n",
+"Q3=1.4*400+1.7*240;\n",
+"Q4=1.4*180+1.7*120;\n",
+"//solving for a\n",
+"deff('y=f(a)','y=0.9*0.51*a*60*(29-a/2)-95.05*12');\n",
+"[x]=fsolve(1.4,f);\n",
+"a=x;\n",
+"As=0.51*a\n",
+"fy=60000;\n",
+"disp(As,'required area in in^2');\n",
+"Asmin=200/fy*12*29;\n",
+"disp(Asmin,'minimum reinforcement required in^2/ft');\n",
+"printf('use No 9 bars at 10 inch centre to centre')\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_Foundation_Engineering_by_B_M_Das/7-Lateral_Earth_Pressure.ipynb b/Principles_Of_Foundation_Engineering_by_B_M_Das/7-Lateral_Earth_Pressure.ipynb
new file mode 100644
index 0000000..e4545a8
--- /dev/null
+++ b/Principles_Of_Foundation_Engineering_by_B_M_Das/7-Lateral_Earth_Pressure.ipynb
@@ -0,0 +1,277 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 7: Lateral Earth Pressure"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.1: 1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 7.1\n",
+"clc; funcprot(0);\n",
+"sigmao=48;\n",
+"phi1=30*%pi/180;\n",
+"phi2=36*%pi/180;\n",
+"Ka1=(tan(%pi/4-phi1/2))^2;\n",
+"Ka2=(tan(%pi/4-phi2/2))^2;\n",
+"sigmaa1=Ka1*sigmao;\n",
+"disp(sigmaa1,'top soil pressure in kN/m^2');\n",
+"sigmaa2=Ka2*sigmao;\n",
+"disp(sigmaa2,'bottom soil pressure in kN/m^2');\n",
+"Po=1/2*3*16+3*12.48+1/3*3*(19.65-12.48)+1/2*3*29.43;\n",
+"zbar=(24*(3+3/3)+37.44*(3/2)+10.76*3/3+44.1*3/3)/Po;\n",
+"disp(zbar,'resultant force acting from the bottom in m');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.2: 2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 7.2\n",
+"clc; funcprot(0);\n",
+"c=14.36;\n",
+"Gamma=17.4;\n",
+"H=6;\n",
+"phi=26*%pi/180;\n",
+"Ka=(tan(%pi/4-phi/2))^2;\n",
+"sigma0=Gamma*H*Ka-2*c*sqrt(Ka);\n",
+"Pa=1/2*Gamma*H^2*Ka-2*c*H*sqrt(Ka);\n",
+"disp(Pa,'active force before which tensile crack appeared in kN/m');\n",
+"zbar=(244.32-323.1)/14.46;\n",
+"disp(zbar,'the line of action on which net force is acting in m');\n",
+"zc=2*c/Gamma/sqrt(Ka);\n",
+"disp(zc,'distance where tensile crack appeared in m');\n",
+"Pa=1/2*(H-zc)*(Gamma*H*Ka-2*c*sqrt(Ka));\n",
+"disp(Pa,'Active force in tensile crack in kN/m');\n",
+"zbar=(H-zc)/3;\n",
+"disp(zbar,'the line of action on which net force is acting in m');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.3: 3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 7.3\n",
+"clc; funcprot(0);\n",
+"pi=%pi\n",
+"H=10;\n",
+"Gamma=110;\n",
+"phi=35*%pi/180;\n",
+"alpha=15*%pi/180;\n",
+"theta=10*%pi/180;\n",
+"zi=asin(sin(alpha)/sin(phi))-alpha+2*theta;\n",
+"disp(zi*180/%pi,'zi in degrees');\n",
+"Ka=cos(alpha-theta)*sqrt(1+(sin(phi))^2-2*sin(phi)*sin(zi))/((cos(theta))^2*(cos(alpha)+sqrt((sin(phi))^2+((sin(alpha))^2))));\n",
+"Pa=1/2*Gamma*H^2*Ka;\n",
+"disp(Pa,'rankine earth pressure in lb/ft');\n",
+"disp('there is slight error in answer due to rounding off error')\n",
+"Beta=atan(sin(phi)*sin(zi)/(1-sin(phi)*cos(zi)));\n",
+"disp(Beta*180/pi,'angle in degrees');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.4: 4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 7.4\n",
+"clc; funcprot(0);\n",
+"H=4.6;\n",
+"Gamma=16.5;\n",
+"Ka=0.297;\n",
+"Po=1/2*Gamma*H^2*Ka;\n",
+"disp(Po,'coulomb active force per unit length in kN/m');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.5: 5.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 7.5\n",
+"clc; funcprot(0);\n",
+"//parta\n",
+"Gamma=105;\n",
+"H=10;\n",
+"Kae=0.474;\n",
+"k1=0;\n",
+"Pae=1/2*Gamma*H^2*Kae*(1-k1)\n",
+"disp(Pae,'active force in lb/ft');\n",
+"//oartb\n",
+"Ka=0.246;\n",
+"Pa=1/2*Gamma*H^2*Ka;\n",
+"disp(Pa,'active force in lb/ft');\n",
+"DPae=Pae-Pa;//deltaPae\n",
+"zbar=(0.6*H*DPae+H/3*(Pa))/Pae;\n",
+"disp(zbar,'the distance of resultant force from bottom in m');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.6: 6.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 7.6\n",
+"clc; funcprot(0);\n",
+"z=[0, 4, 8, 12, 16];\n",
+"Gamma=110;\n",
+"phi=36*%pi/180;\n",
+"H=16;\n",
+"Sa1(1)=0;//sigma(1)\n",
+"Sa2(1)=0;//sigma(2)\n",
+"Sztr(1)=0;//sigma(z)translation\n",
+"printf('z(ft)\t sigma(1)(lb/ft^2)  sigma(2)(lb/ft^2) sigma(z)translation (lb/ft^2)\n');\n",
+"for i=1:5\n",
+"    Sa1(i)=Gamma*(tan(%pi/4-phi*z(i)/2/H))^2*(z(i)-phi*z(i)^2/H/cos(phi*z(i)/H));\n",
+"    Sa2(i)=Gamma*z(i)*(cos(phi)/(1+sin(phi)))^2;\n",
+"    Sztr(i)=Sa1(i)/2+Sa2(i)/2;\n",
+"    printf('%.2f\t %.2f\t\t\t %.2f\t\t\t %.2f\n', z(i),Sa1(i),Sa2(i),Sztr(i));\n",
+"end\n",
+"plot(z,Sztr);\n",
+"xtitle('sigma(z)translation vs z','z(m)','sigma(z)translation (lb/ft^2)')\n",
+"\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 7.7: 7.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 7.7\n",
+"Gammasat=18.86;\n",
+"Gammaw=9.81;\n",
+"clc; funcprot(0);\n",
+"phi1=%pi/180*30;\n",
+"phi2=%pi/180*26;\n",
+"Kp1=(tan(%pi/4+phi1/2))^2;\n",
+"Kp2=(tan(%pi/4+phi2/2))^2;\n",
+"//for top soil\n",
+"c=0;\n",
+"sigma0=31.44;\n",
+"sigmap=sigma0*Kp1+2*c*sqrt(Kp1);\n",
+"disp(sigmap,'passive pressure for top layer in kN/m^2');\n",
+"//for z=2\n",
+"c=10;\n",
+"sigma0=31.44;\n",
+"sigmap=sigma0*Kp2+2*c*sqrt(Kp2);\n",
+"disp(sigmap,'passive pressure for z=2m in kN/m^2');\n",
+"//for z=3\n",
+"c=10;\n",
+"sigma0=15.72*2+(Gammasat-Gammaw)*1;\n",
+"sigmap=sigma0*Kp2+2*c*sqrt(Kp2);\n",
+"disp(sigmap,'passive pressure for z=3m in kN/m^2');\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_Foundation_Engineering_by_B_M_Das/8-Retaining_Walls.ipynb b/Principles_Of_Foundation_Engineering_by_B_M_Das/8-Retaining_Walls.ipynb
new file mode 100644
index 0000000..17b55e3
--- /dev/null
+++ b/Principles_Of_Foundation_Engineering_by_B_M_Das/8-Retaining_Walls.ipynb
@@ -0,0 +1,413 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 8: Retaining Walls"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.1: 1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 8.1\n",
+"clc; funcprot(0);\n",
+"H1=6*tan(10*%pi/180);\n",
+"H2=18;\n",
+"H3=2.75;\n",
+"Gamma1=117;\n",
+"Ka=0.294;//from table\n",
+"H=H1+H2+H3;\n",
+"Pa=1/2*Gamma1*H^2*Ka/1000;\n",
+"Pr=Pa*sin(10*%pi/180);\n",
+"Ph=Pa*cos(10*%pi/180);\n",
+"Mo=Ph*H/3;\n",
+"Sum=0;//sigma Mr\n",
+"S=[1, 2, 3, 4, 5];//section\n",
+"W=[4.05, 1.35, 5.156, 13.01, 1.42];//weight\n",
+"MA=[5.75, 4.67, 6.25, 9.5, 12.5, 12.5];//Moment Arm\n",
+"printf('Section Weight(kip/ft) Moment Arm(ft) Moment (kip-ft/ft)\n');\n",
+"for i=1:5\n",
+"    M(i)=W(i)*MA(i);\n",
+"    Sum=Sum+M(i);\n",
+"    printf('%.2f\t\t %.2f\t\t %.2f\t %.2f\n',S(i),W(i),MA(i),M(i));\n",
+"end\n",
+"FSO=Sum/Mo;\n",
+"if FSO>2 then\n",
+"    printf('safe in overturning with FS=%.2f\n',FSO)\n",
+"end\n",
+"//for sliding\n",
+"phi2=18*%pi/180;\n",
+"V=24.986;\n",
+"B=12.5;\n",
+"c2=0.9;\n",
+"FSS=(V*tan(2/3*phi2)+B*2/3*c2)/Ph;\n",
+"if FSS>2 then\n",
+"    printf('safe in sliding with FS=%.2f\n',FSS)\n",
+"end\n",
+"//for bearing\n",
+"e=B/2-(Sum-Mo)/V;\n",
+"qtoe=V/B*(1+6*e/B);\n",
+"Nc=13.1;\n",
+"Nq=5.26;\n",
+"Ny=4.07;\n",
+"D=0.107;\n",
+"Gamma2=4;\n",
+"B1=B-2*e;//Bdash\n",
+"q=Gamma2*D\n",
+"Fcd=1+0.4*D/B1;\n",
+"Fqd=1+2*tan(phi2)*(1-sin(phi2))^2*(D/B1);\n",
+"Fyd=1;\n",
+"zi=atan(Ph/V);\n",
+"Fci=(1-zi/%pi*2)^2;\n",
+"Fqi=Fci;\n",
+"Fyi=round((1-zi/phi2)^2);\n",
+"qu=c2*Nc*Fcd*Fci+q*Nq*Fqd*Fqi+1/2*Gamma2*B1*Fyd*Fyi;\n",
+"FSB=qu/qtoe;\n",
+"if FSB>3 then\n",
+"    printf('safe in bearing with FS=%.2f\n\n',FSB)\n",
+"end\n",
+"disp('slight changes due to round off error')\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.2: 2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 8.2\n",
+"clc; funcprot(0);\n",
+"c=0.9;\n",
+"B=12.5;\n",
+"Gamma2=4;\n",
+"Fcd=1.138;\n",
+"Fqd=1.107;\n",
+"Nc=7.5;\n",
+"Nq=4;\n",
+"Ny=0;\n",
+"q=0.428;\n",
+"qtoe=2.44;\n",
+"qu=c*Nc*Fcd+q*Nq*Fqd+1/2*Gamma2*B*Ny;\n",
+"FSB=qu/qtoe;\n",
+"if FSB>3 then\n",
+"    printf('safe in bearing with FS=%.2f\n\n',FSB)\n",
+"end\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.3: 3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 8.3\n",
+"clc; funcprot(0);\n",
+"Msum=0;//sum of moment\n",
+"Vsum=0;//sum of force\n",
+"H=15+2.5;//Hdash\n",
+"phi=30*%pi/180;\n",
+"Gamma=121;\n",
+"Ka=(tan(%pi/4-phi/2))^2;\n",
+"Pa=1/2*Gamma*H^2*Ka/1000;\n",
+"Ph=Pa;\n",
+"Mo=Ph*H/3;\n",
+"AN=[1,2,3,4,5,6];//area number\n",
+"W=[0.9,3.375,5.906,3.863,4.764,2.723];//weight\n",
+"MA=[1.783,2.8,5.3,5.15,7.05,9.55];//moment arm\n",
+"\n",
+"printf('AreaNo Weight(kip/ft) Moment Arm(ft) Moment (kip-ft/ft)\n');\n",
+"for i=1:6\n",
+"    M(i)=W(i)*MA(i);\n",
+"    Vsum=Vsum+W(i);\n",
+"    Msum=Msum+M(i);\n",
+"    printf('%.2f\t\t %.2f\t\t %.2f\t %.2f\n',AN(i),W(i),MA(i),M(i));\n",
+"end\n",
+"FSO=Msum/Mo;\n",
+"if FSO>2 then\n",
+"    printf('safe in overturning with FS=%.2f\n',FSO)\n",
+"end\n",
+"//for sliding\n",
+"phi2=20*%pi/180;\n",
+"V=Vsum\n",
+"B=10.3;\n",
+"c2=1.0;\n",
+"FSS=(V*tan(2/3*phi2)+B*2/3*c2)/Ph;\n",
+"printf('safe in sliding with FS=%.2f\n',FSS)\n",
+"e=B/2-(Msum-Mo)/V;\n",
+"qtoe=V/B*(1+6*e/B);\n",
+"disp(qtoe,'soil pressure at toe in kip/ft^2')\n",
+"qheel=V/B*(1-6*e/B);\n",
+"disp(qheel,'soil pressure at heel in kip/ft^2')\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.4: 4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 8.4\n",
+"clc; funcprot(0);\n",
+"Msum=0;//sum of moment\n",
+"Vsum=0;//sum of force\n",
+"H=5+1.5;//Hdash\n",
+"phi=32*%pi/180;\n",
+"Gamma=18.5;\n",
+"Ka=0.424;\n",
+"Pa=1/2*Gamma*H^2*Ka;\n",
+"Ph=Pa*cos(15*%pi/180+2/3*phi);\n",
+"Mo=Ph*H/3;\n",
+"AN=[1,2,3,4,5];//area number\n",
+"A=[4.36, 3.42, 0.77, 2.8, 2.8];//area\n",
+"W=[102.81, 80.64, 18.16, 66.02, 93.14];//weight\n",
+"MA=[2.18, 1.37, 0.98, 1.75, 2.83];//moment arm\n",
+"printf('AreaNo Area(m^2) Weight(kN/m) Momwnt Arm(m) Moment (kN-m/m)\n');\n",
+"for i=1:5\n",
+"    M(i)=W(i)*MA(i);\n",
+"    Vsum=Vsum+W(i);\n",
+"    Msum=Msum+M(i);\n",
+"    printf('%.2f\t\t %.2f\t %.2f\t\t %.2f\t %.2f\n',AN(i),A(i),W(i),MA(i),M(i));\n",
+"end\n",
+"FSO=Msum/Mo;\n",
+"if FSO>2 then\n",
+"    printf('safe in overturning with FS=%.2f\n',FSO)\n",
+"end\n",
+"//for sliding\n",
+"phi2=24*%pi/180;\n",
+"V=Vsum\n",
+"B=3.5;\n",
+"c2=30;\n",
+"Pp=1/2*2.37*18*1.5^2+2*30*1.54*1.5;\n",
+"FSS=(V*tan(2/3*phi2)+B*2/3*c2+Pp)/Ph;\n",
+"printf('safe in sliding with FS=%.2f\n',FSS)\n",
+"disp('if Pp is ignored then FS=1.37');\n",
+"e=B/2-(Msum-Mo)/V;\n",
+"qtoe=V/B*(1+6*e/B);\n",
+"disp(qtoe,'soil pressure at toe in kN/m^2')\n",
+"qheel=V/B*(1-6*e/B);\n",
+"disp(qheel,'soil pressure at heel in kN/m^2')\n",
+"disp('there is difference in answer due to rounding off error')\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.5: 5.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 8.5\n",
+"clc; funcprot(0);\n",
+"//from graph\n",
+"D15=0.04;\n",
+"D50=0.13;\n",
+"D85=0.25;\n",
+"printf('grain size  from grain size distribution D15B in mm =%.2f, D85B in mm = %.2f, D50B in mm = %.2f\n',D15,D85,D50);"
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.6: 6.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 8.6\n",
+"clc; funcprot(0);\n",
+"Sv=2;\n",
+"Sh=3;\n",
+"w=3/12;\n",
+"fy=35000*144;\n",
+"FSb=3;\n",
+"pi=%pi;\n",
+"phi=36*pi/180;\n",
+"Gamma1=105;\n",
+"H=30;\n",
+"t=Gamma1*H*Sv*Sh*FSb/w/fy*(tan(pi/4-phi/2))^2;\n",
+"t=t*12;\n",
+"disp(t,'thiskness in inches');\n",
+"t=t+0.001*50;\n",
+"disp('so take thickness=0.2 inches');\n",
+"//for tie length\n",
+"z=[5,10,15,20,25,30];\n",
+"TL=[38.45, 35.89, 33.34, 30.79, 28.25, 25.7];//tie length\n",
+"printf('z(ft)\t Tie Length (ft)\n');\n",
+"for i=1:6\n",
+"    printf('%.2f\t %.2f\n', z(i),TL(i))\n",
+"end\n",
+"disp('use tie length=40ft')\n",
+"//check for over turning\n",
+"\n",
+"z=30/3;\n",
+"x1=20;\n",
+"L=40;\n",
+"Ka=0.26;\n",
+"Pa=1/2*Gamma1*Ka*H^2;\n",
+"W1=Gamma1*H*L;\n",
+"FSO=W1*x1/(Pa*z);\n",
+"disp(FSO,'factor of safety is');\n",
+"disp('since FS>3 structure is safe')\n",
+"//check for sliding\n",
+"k=2/3;\n",
+"Pa=12285;\n",
+"FSS=W1*tan(k*phi)/Pa;\n",
+"if FSS>3 then\n",
+"    printf('safe in sliding with FS=%.2f\n',FSS)\n",
+"end\n",
+"//check for bearing\n",
+"Mr=126000*20;\n",
+"Mo=12285*10;\n",
+"V=126000;\n",
+"e=L/2-Mr/V+Mo/V;\n",
+"L1=L-2*e;//Ldash\n",
+"c2=1000;\n",
+"Nc=25.8;\n",
+"Gamma2=110;\n",
+"Ny=16.72;\n",
+"qult=c2*Nc+1/2*Gamma2*L1*Ny\n",
+"sigma0=Gamma1*H;\n",
+"FSB=qult/sigma0;\n",
+"if FSB>5 then\n",
+"    printf('safe in bearing with FS=%.2f\n\n',FSB)\n",
+"end\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 8.7: 7.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 8.7\n",
+"clc; funcprot(0);\n",
+"pi=%pi;\n",
+"phi=36*pi/180;\n",
+"Ka=(tan(pi/4-phi/2))^2;\n",
+"z=[8, 12, 16];\n",
+"sigmaG=80*12;\n",
+"Gamma1=110;\n",
+"FS=1.5;\n",
+"for i=1:3\n",
+"    Sv(i)=sigmaG/Gamma1/z(i)/Ka/FS*12;\n",
+"    printf('for z=%.2f ft Sv = %.2f inches\n',z(i),Sv(i));\n",
+"end\n",
+"z=[16,56,76,96,112,144,176];\n",
+"zf=z/12;//z in ft\n",
+"Sv=[1.67,1.67,1.67,1.67,1.33,1.33,1.33];\n",
+"k=[7.48,5.78,4.93,4.08,3.4,2.04,0.68];//0.51(H-z)\n",
+"printf('z(in)    z(ft)  Sv(ft)  0.51(H-z)(ft)  0.438Sv(ft) L(ft) \n')\n",
+"for i=1:7\n",
+"    k2(i)=0.438*Sv(i);//0.438Sv\n",
+"    L(i)=k(i)+k2(i);\n",
+"    printf('%.2f\t %.2f\t %.2f\t %.2f\t\t %.2f\t   %.2f\n',z(i),zf(i),Sv(i),k(i),k2(i),L(i));\n",
+"end\n",
+"Sv=20/12;\n",
+"Ka=0.26;\n",
+"FS=1.5;\n",
+"l1=Sv*Ka*FS/4/tan(2/3*phi);\n",
+"if l1<3 then\n",
+"    l1=3;\n",
+"    disp(l1,'length in ft')\n",
+"end\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_Foundation_Engineering_by_B_M_Das/9-Sheet_Pile_Piles.ipynb b/Principles_Of_Foundation_Engineering_by_B_M_Das/9-Sheet_Pile_Piles.ipynb
new file mode 100644
index 0000000..b6c7841
--- /dev/null
+++ b/Principles_Of_Foundation_Engineering_by_B_M_Das/9-Sheet_Pile_Piles.ipynb
@@ -0,0 +1,403 @@
+{
+"cells": [
+ {
+		   "cell_type": "markdown",
+	   "metadata": {},
+	   "source": [
+       "# Chapter 9: Sheet Pile Piles"
+	   ]
+	},
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.1: 1.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 9.1\n",
+"clc; funcprot(0);\n",
+"sall=30;//sigma allowed\n",
+"pi=%pi;\n",
+"zbar=12.1;\n",
+"L1=10;\n",
+"L2=20;\n",
+"Gamma=0.12;\n",
+"Gamma1=0.1294-0.0624;\n",
+"phi=40*pi/180;\n",
+"Ka=(tan(pi/4-phi/2))^2;\n",
+"Kp=(tan(pi/4+phi/2))^2;\n",
+"s1=Gamma*L1*Ka;//sigma1\n",
+"s2=Gamma*L1*Ka+Gamma1*L2*Ka;//sigma2\n",
+"L3=s2/(Gamma1*(Kp-Ka));\n",
+"disp(L3,'length in ft');\n",
+"P=1/2*s1*L1+s1*L2+1/2*(s2-s1)*L2+1/2*s2*L3;\n",
+"disp(P,'force in kip/ft');\n",
+"s5=Gamma*L1*Kp+Gamma1*L2*Kp+Gamma*L3*(Kp-Ka);//sigma5\n",
+"disp(s5,'pressure in kip/ft');\n",
+"A1=s5/(Gamma1*(Kp-Ka));\n",
+"A2=8*P/(Gamma1*(Kp-Ka))\n",
+"A3=6*P*(2*zbar*(Gamma1*(Kp-Ka))+s5)/(Gamma1*(Kp-Ka))^2\n",
+"A4=P*(6*zbar*s5+4*P)/(Gamma1*(Kp-Ka))^2\n",
+"printf('A1,A2,A3,A4 respectively in %.2f,%.2f,%.2f,%.2f \n',A1,A2,A3,A4);\n",
+"disp('slight error due to rounding off error')\n",
+"//partb\n",
+"deff('y=f(x)','y=x^4+41.7*x^3-270.5*x^2-13363*x-106863');\n",
+"[x]=fsolve(20,f);\n",
+"D=1.88+x;\n",
+"disp(D,'value of D, ft')\n",
+"TL=L1+L2+1.3*D;\n",
+"disp(TL,'total length in ft');\n",
+"//partc\n",
+"z=sqrt(2*P/(Gamma1*(Kp-Ka)));//zdash\n",
+"Mmax=P*(z+zbar)-1/2*(Gamma1*(Kp-Ka))*z^2*z/3;\n",
+"S=Mmax*12/sall;\n",
+"disp(S,'section modulus in in^3/ft')\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.2: 2.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 9.2\n",
+"clc; funcprot(0);\n",
+"sall=172.5*1000;//sigma allowed\n",
+"pi=%pi;\n",
+"c=47;\n",
+"zbar=1.78;\n",
+"L1=2;\n",
+"L2=3;\n",
+"Gamma=15.9;\n",
+"Gamma1=19.33-9.81;\n",
+"phi=32*pi/180;\n",
+"Ka=(tan(pi/4-phi/2))^2;\n",
+"Kp=(tan(pi/4+phi/2))^2;\n",
+"s1=Gamma*L1*Ka;//sigma1\n",
+"s2=Gamma*L1*Ka+Gamma1*L2*Ka;//sigma2\n",
+"P=1/2*s1*L1+s1*L2+1/2*(s2-s1)*L2;\n",
+"disp(P,'force in kN/ft');\n",
+"deff('y=f(x)','y=127.4*x^2-104.4*x-357.15');\n",
+"[x]=fsolve(2,f);\n",
+"D=x;\n",
+"disp(D,'value of D, m')\n",
+"disp(D*1.5,'actual D in m')\n",
+"L4=D*(4*c-(Gamma*L1+Gamma1*L2)-P/D)/4/c;\n",
+"disp(L4,'length in m');\n",
+"s6=4*c-(Gamma*L1+Gamma1*L2);//sigma6\n",
+"s7=4*c+(Gamma*L1+Gamma1*L2);//sigma7\n",
+"z=P/s6;//zdash\n",
+"Mmax=P*(z+zbar)-1/2*s6*z^2;\n",
+"S=Mmax*12/sall;\n",
+"disp(S,'section modulus in m^3/m')\n",
+"disp('slight error due to rounding off error')\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.3: 3.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 9.3\n",
+"clc; funcprot(0);\n",
+"pi=%pi;\n",
+"zbar=2.23;\n",
+"L1=2;\n",
+"L2=3;\n",
+"Gamma=15.9;\n",
+"Gamma1=19.33-9.81;\n",
+"phi=32*pi/180;\n",
+"Ka=(tan(pi/4-phi/2))^2;\n",
+"Kp=(tan(pi/4+phi/2))^2;\n",
+"s1=Gamma*L1*Ka;//sigma1\n",
+"s2=Gamma*L1*Ka+Gamma1*L2*Ka;//sigma2\n",
+"L3=s2/(Gamma1*(Kp-Ka));\n",
+"disp(L3,'length in m');\n",
+"P=1/2*s1*L1+s1*L2+1/2*(s2-s1)*L2+1/2*s2*L3;\n",
+"disp(P,'force in kN/m');\n",
+"deff('y=f(x)','y=x^3+6.99*x^2-14.55');\n",
+"[x]=fsolve(1.4,f);\n",
+"D=L3+x;\n",
+"disp(D,'value of D in m')\n",
+"AL=1.4*D;\n",
+"disp(AL,'actual length in m');\n",
+"//partb\n",
+"L4=1.4;\n",
+"F=P-1/2*(Gamma1*(Kp-Ka)*L4^2);\n",
+"disp(F,'anchor force in kN/m');\n",
+"//partc\n",
+"deff('y=f(x)','y=x^2+6.682*x-14.44');\n",
+"[x]=fsolve(1.7,f);\n",
+"z=x+2;\n",
+"Mmax=-1/2*s1*L1*(x+2/3)+F*(x+1)-s1*x*x/2-1/2*Ka*Gamma1*x^3/3;\n",
+"disp(Mmax,'maximum moment in kN-m/m')\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.4: 4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 9.4\n",
+"clc; funcprot(0);\n",
+"Gamma=15.9;\n",
+"Gamma1=19.33-9.81;\n",
+"GD=0.23;\n",
+"CDL1=1.172;\n",
+"L1=2;\n",
+"L2=3;\n",
+"Dth=(L1+L2)*GD*CDL1;\n",
+"disp(Dth,'theoritical depth in m');\n",
+"Dac=1.4*Dth;\n",
+"disp(Dac,'actual depth in m');\n",
+"printf('approximate it as 2 m');\n",
+"//partb\n",
+"CFL1=1.073;\n",
+"GF=0.07;\n",
+"Gammaa=(Gamma*L1^2+Gamma1*L2^2+2*Gamma*L1*L2)/(L1+L2)^2;\n",
+"F=Gammaa*(L1+L2)^2*GF*CFL1;\n",
+"disp(F,'force in kN/m');\n",
+"//partc\n",
+"GM=0.021;\n",
+"CML1=1.036;\n",
+"Mmax=Gammaa*(L1+L2)^3*GM*CML1;\n",
+"disp(Mmax,'maximum moment in kN-m/m');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.5: 5.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 9.5\n",
+"clc; funcprot(0);\n",
+"Mmax=43.72;\n",
+"sp=['PSA-31','PSA-23'];//sheet pile\n",
+"H=[7.9,7.9]\n",
+"I=[4.41e-6,5.63e-6];\n",
+"p=[0.00466,0.00365];\n",
+"S=[10.8e-5,12.8e-5];\n",
+"Md=[18.63,22.08];\n",
+"printf('SheetPile    I(m^4/m)  H(m)  p\t     Logp  S(m^3/m)  Md(kN-m/m) Md/Mmax \n')\n",
+"for i=1:2\n",
+"    Logp(i)=log10(p(i));\n",
+"    k(i)=Md(i)/Mmax;\n",
+"    printf('%s\t   %.2e    %.2f %f %.2f %.2e   %.2f\t %.3f\n',sp(i),I(i),H(i),p(i),Logp(i),S(i),Md(i),k(i));\n",
+"    \n",
+"end\n",
+"plot(Logp,k);\n",
+"xtitle('','LogP','Md/Mmax');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.6: 6.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 9.6\n",
+"clc; funcprot(0);\n",
+"pi=%pi;\n",
+"R=0.6;\n",
+"L1=10;\n",
+"L2=20;\n",
+"Gammasat=122.4;\n",
+"l1=5;\n",
+"Gamma=110;\n",
+"C=0.68;\n",
+"L=L1+L2;;\n",
+"Gammaw=62.4;\n",
+"Gamma1=Gammasat-Gammaw;//gammadash\n",
+"Gammaav=(Gamma*L1+Gamma1*L2)/(L1+L2);\n",
+"phi=35*pi/180;\n",
+"Ka=(tan(pi/4-phi/2))^2;\n",
+"sa=C*Ka*Gammaav*L;//sigmaa\n",
+"sp=R*sa;//sigmap\n",
+"deff('y=f(x)', 'y=x^2+50*x-1000');\n",
+"[x]=fsolve(15,f);\n",
+"D=x;\n",
+"disp(D,'depth in ft');\n",
+"R=L/D*(L-2*l1)/(2*L+D-2*l1);\n",
+"printf('value of R=%.2f is OK\n',R);\n",
+"//partb\n",
+"F=sa*(L-R*D);\n",
+"disp(F,'Force in lb/ft')\n",
+"//partc\n",
+"Mmax=0.5*sa*L^2*((1-R*D/L)^2-2*l1/L*(1-R*D/L));\n",
+"disp(Mmax,'maximum moment lb-ft/ft');\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.7: 7.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 9.7\n",
+"clc; funcprot(0);\n",
+"pi=%pi;\n",
+"zbar=3.2;\n",
+"c=41;\n",
+"L1=3;\n",
+"L2=6;\n",
+"Gamma=17;\n",
+"Gamma1=20-9.81;\n",
+"phi=35*pi/180;\n",
+"Ka=(tan(pi/4-phi/2))^2;\n",
+"Kp=(tan(pi/4+phi/2))^2;\n",
+"s1=Gamma*L1*Ka;//sigma1\n",
+"s2=Gamma*L1*Ka+Gamma1*L2*Ka;//sigma2\n",
+"P=1/2*s1*L1+s1*L2+1/2*(s2-s1)*L2;\n",
+"disp(P,'Force in kN/m');\n",
+"s6=4*c-(Gamma*L1+Gamma1*L2);//sigma6\n",
+"disp(s6,'pressure in kN/m^2');\n",
+"deff('y=f(x)','y=x^2+15*x-25.43');\n",
+"[x]=fsolve(1.6,f);\n",
+"D=x;\n",
+"disp(D,'depth in m');\n",
+"F=P-s6*D;\n",
+"disp(F,'force in kN/m');\n",
+"printf('slight error due to rounding off')\n",
+""
+   ]
+   }
+,
+{
+		   "cell_type": "markdown",
+		   "metadata": {},
+		   "source": [
+			"## Example 9.8: 4.sce"
+		   ]
+		  },
+  {
+"cell_type": "code",
+	   "execution_count": null,
+	   "metadata": {
+	    "collapsed": true
+	   },
+	   "outputs": [],
+"source": [
+"//example 9.8\n",
+"clc; funcprot(0);\n",
+"pi=%pi;\n",
+"Gamma=105;\n",
+"Cov=14;\n",
+"B=15/12;\n",
+"Ka=0.26;\n",
+"phi=35*pi/180;\n",
+"H=37.5/12;\n",
+"h=15/12;\n",
+"t=6/12;\n",
+"Gc=150;//gamma concrete\n",
+"W=H*t*Gc;\n",
+"k=4.5;//kp*cos(delta)\n",
+"Pu=1/2*Gamma*H^2*(k-Ka*cos(phi));\n",
+"disp(Pu,'force in lb/ft');\n",
+"Pus=[(Cov+1)/(Cov+H/h)]*Pu;\n",
+"disp(Pus,'force in lb/ft');\n",
+"Be=0.227*(H+h)+B;\n",
+"Pu=Pus*Be;\n",
+"disp(Pu,'resistance of anchor plate i lb/ft')\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
+}