%

clear all;
close all;
clc;

d1 = 0.35;

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

xbel0 = 0.15;
L  = 30;
ib = 0;
C  = 45;
B  = 0.1;
hb = d1;   % waterdiepte op benedenstroomse rand

%Q = B * hb * C * sqrt( hb * ib );
Q = 0.03;     % input instead of steady flow equivalent
N = 100;      % for Belanger-computation

%function [x,h] = belanger(Q,C,L,B,ib,hb,N);
[xbel,hbel] = belanger(Q,C,L,B,ib,hb,N);
xbel=xbel';
hbel=hbel';
xbel=xbel+xbel0;

plot(xbel,hbel)
