clear all; close all; t = 0:0.0000001:0.00005; s = tf('s'); sys = (1.000000000*10^10)/(s^2 + 200000.*s + 1.000000000*10^10); sys2 = 1/(6.3196375*10^(-13)*s^2 + 1.124*10^(-6)*s + 1); y = step(sys,t); y2 = step(sys2,t); %damp(sys); %damp(sys2); plot(t,y,'-r') hold on plot(t,y2,'-b') xlabel('t'); ylabel('y(t)'); grid on;