s=tf('s')
wfal=1;
for n=2:10
e=sqrt(10^(wfal/10)-1);
u=(1/n)*asinh(1/e) ;
omc=1/(cosh(1/n*acosh(sqrt((2*e^2+1)/e^2))));
if mod(n,2)==0
    zwieksz=1;
    for k=1:(n/2)
    alp=((2*k-1)/(2*n))*pi;
    B=1/(cosh(u)^2-cos(alp)^2);
    A=2*B*sinh(u)*cos(alp);
    a=A/omc;
    b=B/(omc.^2);
    L=1;
    M=[b a 1]
    h=tf(L,M)
    zwieksz=zwieksz*h
    end
    H(n)=zwieksz
else
    zwieksz=1;
    for k=2:((n+1)/2)
        bet=((k-1)/n)*pi;
        B=1/(cosh(u)^2-cos(bet)^2)
        A=2*B*sinh(u)*cos(bet)
        a=A/omc
        b=B/(omc.^2)
        L=1;
        M=[b a 1]
        h=tf(L,M)
        zwieksz=zwieksz*h
    end
    H(n)=zwieksz*(1/((1/sinh(u))*(s+1)))
    end
end