// 2ŽŸŒn (ƒÖn^2)/(s^2+2ƒÄƒÖns+ƒÖn^2) ‚̃{�[ƒh�ü�}

s=poly(0,'s');

zeta = 2^(-0.5);

fn = 1;               // [Hz]
wn = 2*%pi*fn;         // [rad/sec]

h=syslin('c',wn^2/(s^2+2*zeta*wn*s+wn^2));

bode(h);