From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 2087/CH4/EX4.28/example4_28.sce | 32 ++++++++++++++++++++++++++++++++
 1 file changed, 32 insertions(+)
 create mode 100755 2087/CH4/EX4.28/example4_28.sce

(limited to '2087/CH4/EX4.28/example4_28.sce')

diff --git a/2087/CH4/EX4.28/example4_28.sce b/2087/CH4/EX4.28/example4_28.sce
new file mode 100755
index 000000000..4f8e9641f
--- /dev/null
+++ b/2087/CH4/EX4.28/example4_28.sce
@@ -0,0 +1,32 @@
+
+
+//example 4.28
+//calculate ordinates of storm hydrograph
+clc;funcprot(0);
+//given
+fi=2.5;          //infiltration index
+B=10;            //Base flow
+O=[0 110 365 500 390 310 250 235 175 130 95 65 40 22 10 0 0 0];   //ordinates of unit hydrograph
+R1=2;R2=6.75;R3=3.75;
+r1=(R1*10-(fi*3)-5)/10;             //rainfall excess in first three hour
+r2=(R2*10-(fi*3))/10;               //rainfall excess in second three hour
+r3=(R3*10-(fi*3))/10;               //rainfall excess in third three hour
+
+for i=1:18
+    s1(i)=r1*O(i);                
+end
+for i=2:18
+    s2(i)=r2*O(i-1);
+end
+for i=3:18
+    s3(i)=r3*O(i-2);
+end                                 //surface run-off from rainfall excess during succesive unit periods
+mprintf("ordinates of storm hydrograph");
+for i=1:18
+    T(i)=s1(i)+s2(i)+s3(i);
+    t(i)=T(i)+B;
+    t(i)=round(t(i)*10)/10;
+    mprintf("\n%f",t(i));
+end
+
+
-- 
cgit