Matlab reference code

num=1;
den=[1 10 20];
Kp=350;
Ki=300;
Kd=50;
contr=tf([Kd Kp Ki],[1 0]);
plant=tf(num,den);
sys_cl=feedback(contr*plant,1);
t=0:0.01:2;
step(sys_cl,t);