%H = tf([10],[10 2*0.7*sqrt(10) 1 0],'inputdelay',0.35); H = tf([10],[10 2*0.7*sqrt(10) 1 0]) Hd = c2d(H,0.01,'foh') Hd2 = c2d(H,0.1,'foh') %step(H,'-',Hd,'--') bode(H, Hd, Hd2)