s=%s;
num=2;
den=s*(s+1);
gr=syslin('c',num/den);
cr=gr/. 1;
t=0:0.01:10;
r1=(1/2)*t^2;
c1=csim(r1,t,cr);
clf;
plot2d(t,c1,2);
plot2d(t,r1,3);
xtitle('Close-loop system response to input r(t)=(1/2*t^2) ','t(seg)','Amplitude');
legends(['Input','Output'],[3,2],opt=4)
xgrid;