%

%function [x,h] = belanger(Q,C,L,B,ib,hb,N);

dh = 0.012;
x0 = 0.3;
L  = 29.7;
ib = dh/L;
C  = 45;
B  = 0.1;
hb = 0.402;   % waterdiepte op benedenstroomse rand

Q = B * hb * C * sqrt( hb * ib )

N = 300;      % for Belanger-computation

[x,h] = belanger(Q,C,L,B,ib,hb,N);
x=x+x0;

close all;
figure;

plot(x,h);

mean(h)
