% obiekt
A1=4;
A2=5;
g=9.81;
Aw1=0.5;
Aw2=0.5;
hdiffmax=2;
hmax2=6;
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];
Cz3 = [0 1];
Dz3 = 0;
Z3=ss(Az3,Bz3,Cz3,Dz3);
figure(1)
step(Z3)
% model kupfmullera
tau3=0.7;
T3=15.9-tau3;
k3=1.1;
G3=tf([k3],[T3 1]);
figure(2)
step(G3, 'r')
% zestawienie
figure(3)
step(Z3)
hold on
step(G3, 'r')