wfalp=3;  % współczynnik falistości w paśmie przepustowym
wfalz=20; % współczynnik falistości w paśmie zaporowym
for n=2:9
[z,p,k]=ellipap(n,wfalp,wfalz);
[L,M]=zp2tf(z,p,k)
H(n)=tf(L,M)
end