%% Zbiornik 3
A1=4;
A2=5;
g=9.81;
Aw1=2;
Aw2=2;
hdiffmax=9;
hmax2=10;
a1z3=Aw1*sqrt(2.*g).*hdiffmax.^-0.5;
a2z3=Aw2*sqrt(2.*g).*hmax2.^-0.5;
Az3 = [-a1z3/A1, a1z3/A1; a1z3/A2, (-a1z3/A2-a2z3/A2)];
Bz3 = [1/A1; 0];
Cz5 = [0 1];
Dz5 = 0;
Z3=ss(Az3,Bz3,Cz5,Dz5);
figure(1)
step(Z3)
figure(2)
impulse(Z3)
figure(3)
nyquist(Z3)
figure(4)
bode(Z3)
figure(5)
pzmap(Z3)