%regulator
Kr=0.2;
Ti=25*Kr/(10*Kr+1)+0.5

%proces
num=[10];
den=[5 2 1];
[A,B,C,D]=tf2ss(num,den)
C1=[C;1 0;0 1]
D1=[D;0;0]

%ci
x0=linsolve(C,10);
x30=1;

%simulare+grafice
tf=50;
sim('labS5_simulare')
plot(w.time,w.data,y.time,y.data,c.time,c.data)
hold on
plot(x1.time,x1.data,':k',x2.time,x2.data,':b',x3.time,x3.data,':g')