%% Zbiornik 4
A1=4;
A2=5;
g=9.81;
Aw1=2;
Aw2=2;
hdiffmax=9;
hmax2=10;
a1z4=Aw1*sqrt(2.*g).*hdiffmax.^-0.5;
Az4 = [-a1z4/A1, a1z4/A1; a1z4/A2, -a1z4/A2];
Bz4 = [1/A1; 1/A2];
Cz4 = [0 1];
Dz4 = 0;
Z4=ss(Az4,Bz4,Cz4,Dz4);
figure(1)
step(Z4)
figure(2)
impulse(Z4)
figure(3)
nyquist(Z4)
figure(4)
bode(Z4)
figure(5)
pzmap(Z4)