clear all;
G = tf(1,[1 0.2 7]);
time = [0:0.001:100];
step_response = step(G,time);

figure;
plot(time,step_response);
xlabel('t (sec)');
ylabel('y');