//example 14.28
//design a regime channel using Laecy's theory
clc;funcprot(0);
//given
Q=35;           //discharge
f=0.9;          //silt factor
s=1/2;          //side slope

V=(Q*f/140)^(1/6);
A=Q/V;
P=4.75*Q^0.5;
D=(P-(P^2-6.944*A)^0.5)/3.472;
B=P-2.236*D;

R=5*V^2/(2*f);
S=f^(5/3)/(3340*Q^(1/6));
D=round(D*100)/100;
mprintf("Bed slope=%f.",S);
mprintf("\nWidth of channel section=%i m.",B);
mprintf("\nDepth of channel section=%f m.",D);