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)