function manning = manning(A)
%M=[.035,.0001,10,2;.020,.0002,8,1;.015,.0010,20,1.5;030,.0007,24,3;.022,.0
%003,15,2.5]
for I = 1:length(A),
    n(I)=A(I,1);
    S(I)=A(I,2);
    B(I)=A(I,3);
    H(I)=A(I,4);
end

for I=1:length(n),
    U(I) = (S(I)^.5)*(B(I)*H(I)/(B(I)+2*H(I)))^(2/3)/n(I);
end

for I = 1:length(A),
    A(I,5)=U(I);
end
fprintf('      n          S        B         H         U\n')
disp(A)
end