1 files changed, 28 insertions, 0 deletions
diff --git a/2858/CH3/EX3.2/Ex3_2.sce b/2858/CH3/EX3.2/Ex3_2.sce
new file mode 100755
index 000000000..b40e27e83
--- /dev/null
+++ b/2858/CH3/EX3.2/Ex3_2.sce
@@ -0,0 +1,28 @@
+//example 3.2
+clc; funcprot(0);
+Gamma=105;//lb/ft^3
+Gammasat=118;//lb/ft^3
+FS=3;
+pa=2014.125;//lb/ft^2
+Depth=[5,10,15,20,25];
+N60=[4,6,6,10,5];
+sigmao=[0,0,0,0,0];
+phi=[0,0,0,0,0]
+Gammaw=62.4;
+s=0;
+printf("depth (ft)\tN60\t  \tstress(lb/ft^2)\t phi(degrees)\n")
+for i=1:5
+    sigmao(i)=2*Gamma+(Depth(i)-2)*(Gammasat-Gammaw);
+    phi(i)=sqrt(20*N60(i)*sqrt(pa/sigmao(i)))+20;
+    printf(" %.2f\t       %.2f\t\t   %.2f  \t%.2f \n",Depth(i),N60(i),sigmao(i),phi(i));
+avgphi=phi(i)/5+s;
+s=avgphi;
+end
+disp(round(avgphi),"average friction angle in degrees");
+//using graph get the values of other terms in terms of B and solve for B
+deff('y=f(x)','y=-150000/x^2+5263.9+5527.1/x+228.3*x');
+[x, v, info ]=fsolve(4,f);
+disp(x,"the width in ft");
+
+
+