1 files changed, 41 insertions, 0 deletions
diff --git a/2087/CH14/EX14.26/example14_26.sce b/2087/CH14/EX14.26/example14_26.sce
new file mode 100755
index 000000000..c014a8673
--- /dev/null
+++ b/2087/CH14/EX14.26/example14_26.sce
@@ -0,0 +1,41 @@
+
+
+//example 14.26
+//calculate bed width of channel;also check depth using Kutter equation
+clc;funcprot(0);
+//given
+Q=5;              //discharge
+S=0.2/1000;       //bed slope
+m=0.8;            //critical velocity ratio
+s=1/2;            //side slope of chanel
+C=30;
+//assuming
+D=1;
+Vo=0.55*m*D^0.64;
+A=Q/Vo;
+B=A/D-(s*D);
+P=B+2.43*D;
+R=A/P;
+V=C*(R*S)^0.5;
+//Vo>V
+//hence take second trial
+D=0.8;     //assume
+Vo=0.55*m*D^0.64;
+A=Q/Vo;
+B=A/D-(s*D);
+P=B+2.43*D;
+R=A/P;
+V=C*(R*S)^0.5;
+//again Vo>V
+//hence we take third trial
+D=0.7;
+Vo=0.55*m*D^0.64;
+A=Q/Vo;
+B=A/D+(s*D);
+P=B+2.43*D;
+R=A/P;
+V=C*(R*S)^0.5;
+B=round(B*100)/100;
+//Vo is almost equal to V;
+mprintf("Width of channel section=%f m.",B);
+mprintf("\nDepth of channel section=%f m.",D);