m=1;
k=3;
b=4;
s=tf('s');
gs=1/(m*s^2+b*s+k)