zd=[0.05 0.1 0.2 0.707];
wn=1;
w=logspace(-0.1,0.1,1e4)
for k=zd
sysH=tf(wn^2,[1 2*k*wn wn^2])
bode(sysH,w)
grid on
hold on
end