m = 1e6;

wp = 2*pi*m * [ 0.7 1.3 ];
ws = 2*pi*m * [ 0.9 1.1 ];
Rp = 5;
Rs = 40;

[N, Wn] = cheb2ord(wp, ws, Rp, Rs, 's');
N

return;

[b, a] = cheby2(N, Rs, Wn, 'stop', 's');
h = tf(b, a);

bodehz(h);

return;

[z, p, k] = zpkdata(h);
z = z{1};  p = p{1};
z 
p 
k

figure; bode3(h, 2e6, 4, 50); fixfig;

figure; pzmap(h); fixfig;
