w=1;
Zr = 1.38; Ti = 2.93; Td = 0.80;
P = Zr; I = Zr/Ti; D = Zr*Td;
sim('sys4pid');

clear m6zn;
t=res(:,1);
w=res(:,2);
u=res(:,3);
y=res(:,4);
plot(t,w,'k:',t,y,'k-',t,u,'k--');
v=axis;
%axis([v(1) v(2) v(3) 1.5]);
legend(' w', ' y', ' u');
ylabel('w,y,u');

xlabel('t');
set(1, 'PaperUnits', 'centimeters')
set(1, 'PaperPosition', [3 10 11 8.25])
print -deps 'obr6str.eps'
