s=%s;
num=10*s+4;
den=4+4*s+s^2;
g=syslin('c',num/den);
t=0:0.01:10;
gs=csim('step',t,g);
gr=csim(t,t,g);
subplot(2,1,1);
xgrid;
xtitle('Response to unit-ramp','t(seg)','Amplitude')
plot2d(t,gr,3);
plot2d(t,t,4);
legends(['Input','Output'],[4,3],opt=4);
subplot(2,1,2);
plot2d(t,gs);
xgrid;
xtitle('Response to unit-step','t(seg)','Amplitude')