clear all;
close all;

%%% Parametry

g = 9.81;

A1 = 25;
A2 = 50;
Aw1 = 2;
Aw2 = 1;
hmax1 = 20;
hmax2 = 15;

t0 = 5;
czas = 50;

f_wy_max1 = Aw1*sqrt(2*g*(hmax1-hmax2));
f_wy_max2 = Aw2*sqrt(2*g*hmax2);

fwe1=0;
df1=0.1*f_wy_max1;

df2=0;

fwe2=0;

h20=((fwe1+fwe2)^2)/((Aw2^2)*2*g);
h10=(fwe1^2)/((Aw1^2)*2*g)+h20;

model='kaskada1a';
[z]=sim(model,czas);

figure(1);
hold on;
grid on;
plot(t,h1,'r')

figure(2)
hold on;
grid on;
plot(t,h2,'b')

figure(3)
hold on;
grid on;

fwy1=Aw1*sqrt(2*g*(h1-h2));
plot(h1-h2,fwy1,'g')

fwe1=0.5*f_wy_max1;

h20=((fwe1+fwe2)^2)/((Aw2^2)*2*g);
h10=(fwe1^2)/((Aw1^2)*2*g)+h20;

[z]=sim(model,czas);


figure(1)
plot(t,h1,'r')

figure(2)
plot(t,h2,'b')

figure(3)
hold on;
grid on;
fwy1=Aw1*sqrt(2*g*(h1-h2));
plot(h1-h2,fwy1,'g')

fwe1=0.9*f_wy_max1;
h20=((fwe1+fwe2)^2)/((Aw2^2)*2*g);
h10=(fwe1^2)/((Aw1^2)*2*g)+h20;
[z]=sim(model,czas);
%hold on;
%grid on;
%plot(t,h2)
figure(1)
plot(t,h1,'r')

figure(2)
plot(t,h2)

figure(3)
hold on;
grid on;
fwy1=Aw1*sqrt(2*g*(h1-h2));
plot(h1-h2,fwy1);