[zeros (14,2); Br];
Cz = [C zeros (2,18)];
x = zeros (14,1); xr = zeros (18,1); u = zeros (2,1);
yy = []; uu = []; f = [.0010; .0010]; z = [0, 0];
for i = 1:2000,
y = C * x; e =-z + y;
u = Cr * xr; xr = Ar * xr + Br * e;
y = C * x; x = Ad * x + Bd * (u + f);
yy = [yy; y ']; uu = [uu; u'];
end
x1 = x; xr1 = xr; u1 = u;
figure (2)
subplot (2,2,1); plot (yy (:, 1)); grid; ylabel ('y1, MPa'); title ('Perehod proces Braga');
subplot (2,2,3); plot (yy (:, 2)); grid; ylabel ('y2,');
subplot (2,2,2); plot (uu (:, 1)); grid; ylabel ('u1, M3/c'); title ('Perehod proces Par');
subplot (2,2,4); plot (uu (:, 2)); grid; ylabel ('u2, kg/c');