K = 4;

%

G_cl = tf([9 * K], [1, 7, 21, 36 + 9 * K])

%

pzmap(G_cl)

%

pole(G_cl)

%

[x, t] = step(G_cl);
plot(t, x)

%

stepinfo(x, t)