% model obiektu
A1=4;
A2=5;
g=9.81;
Aw1=0.5;
Aw=0.5;
h1max=6;
h2max=6;
a1=Aw1*sqrt(2.*g).*h1max.^-0.5;
a2=Aw*sqrt(2.*g).*h2max.^-0.5;
A = [-a1/A1, 0; a1/A2, -a2/A2];
B = [1/A1; 0];
C = [0 1];
D = 0;
sys=ss(A,B,C,D);
figure(1)
step(sys)
% model kupfmullera
tau=0.3;
T=15.2-tau;
k=1.11;
G1=tf([k],[T 1]);
figure(2)
step(G1, 'r')
% zestawienie
figure(3)
step(sys)
hold on
step(G1, 'r')