diff options
Diffstat (limited to '2087/CH4/EX4.35/example4_35.sce')
| -rwxr-xr-x | 2087/CH4/EX4.35/example4_35.sce | 26 |
1 files changed, 26 insertions, 0 deletions
diff --git a/2087/CH4/EX4.35/example4_35.sce b/2087/CH4/EX4.35/example4_35.sce new file mode 100755 index 000000000..af3ef0531 --- /dev/null +++ b/2087/CH4/EX4.35/example4_35.sce @@ -0,0 +1,26 @@ +
+
+//example 4.35
+//calculate 200 year flood for stream using gumbel's method
+clc;funcprot(0);
+//given
+sigma=1.1413; //standard deviation
+yn=0.5436;
+T=50;
+y50=-2.303*log10(2.303*log10(T/(T-1)));
+K50=(y50-yn)/sigma;
+T=100;
+y100=-2.303*log10(2.303*log10(T/(T-1)));
+K100=(y100-yn)/sigma;
+x50=878; x100=970; //given peak flood
+A=[K50 1;K100 1];
+B=[x50;x100];
+C=A\B;
+xavg=C(2);
+sigmad=C(1);
+T=200;
+y200=-2.303*log10(2.303*log10(T/(T-1)));
+K200=(y200-yn)/sigma;
+x200=xavg+K200*sigmad;
+x200=round(x200);
+mprintf("200 year flood for stream=%f cumecs.",x200);
|