1 files changed, 33 insertions, 0 deletions
diff --git a/1379/CH10/EX10.1.4/example10_4.sce b/1379/CH10/EX10.1.4/example10_4.sce
new file mode 100755
index 000000000..5cf847dbb
--- /dev/null
+++ b/1379/CH10/EX10.1.4/example10_4.sce
@@ -0,0 +1,33 @@
+
+
+//example 10.4
+clc; funcprot(0);
+//exapple 10.4
+// Initialization of Variable
+rho=998;
+w0=40;//density of slurry
+mu=1.01/1000;
+g=9.81;
+rho1=2660;//density quartz
+h=0.25;
+t=18.5*60;
+mp=[5 11.8 20.2 24.2 28.5 37.6 61.8];
+d=[30.2 21.4 17.4 16.2 15.2 12.3 8.8]/1000000;
+u=h/t;
+d1=sqrt(18*mu*u/g/(rho1-rho));
+function[a]=inter(d,f,g,b);//interpolation linear
+    for i=1:b
+        if d>f(i+1)& d<=f(i)  then
+            break
+        else 
+            continue
+        end
+        break
+    end
+    
+    a=-(d-f(i+1))/(f(i)-f(i+1))*(g(i+1)-g(i))+g(i+1);
+endfunction
+[a]=inter(d1,d,mp,6);
+phi=1-a/100;
+rhot=phi*(rho1-rho)/rho1*w0+rho;
+disp(rhot,"the density of suspension at depth 25cm in kg/m^3 is")