clear all;
n = 2;     %II rząd
for k = 1:1:(n/2)
    a(k) = 2*cos((2*k-1)/(2*n)*pi);
    b(k) = 1;
    num = 1;
    den = [b(k),a(k),1]
    B(k) = tf(num,den);
end
b2 = B(1);
n = 3;     %III rząd
for k = 1:1:((n+1)/2)
    a(k) = 2*cos(((k-1)/n)*pi);
    b(k) = 1;
    b(1) = 0;
    a(1) = 1;
    num = 1;
    den = [b(k),a(k),1];
    B(k) = tf(num,den);
end
b3 = B(1)*B(2);
n = 4;     %IV rząd
for k = 1:1:(n/2)
    a(k) = 2*cos((2*k-1)/(2*n)*pi);
    b(k) = 1;
    num = 1;
    den = [b(k),a(k),1];
    B(k) = tf(num,den);
end
b4 = B(1)*B(2);
n = 5;     %V rząd
for k = 1:1:((n+1)/2)
    a(k) = 2*cos(((k-1)/n)*pi);
    b(k) = 1;
    b(1) = 0;
    a(1) = 1;
    num = 1;
    den = [b(k),a(k),1];
    B(k) = tf(num,den);
end
b5 = B(1)*B(2)*B(3);
n = 6;     %VI rząd
for k = 1:1:(n/2)
    a(k) = 2*cos((2*k-1)/(2*n)*pi);
    b(k) = 1;
    num = 1;
    den = [b(k),a(k),1]
    B(k) = tf(num,den);
end
b6 = B(1)*B(2)*B(3);
n = 7;     %VII rząd
for k = 1:1:((n+1)/2)
    a(k) = 2*cos(((k-1)/n)*pi);
    b(k) = 1;
    b(1) = 0;
    a(1) = 1;
    num = 1;
    den = [b(k),a(k),1];
    B(k) = tf(num,den);
end
b7 = B(1)*B(2)*B(3)*B(4);
n = 8;     %VIII rząd
for k = 1:1:(n/2)
    a(k) = 2*cos((2*k-1)/(2*n)*pi);
    b(k) = 1;
    num = 1;
    den = [b(k),a(k),1]
    B(k) = tf(num,den);
end
b8 = B(1)*B(2)*B(3)*B(4);
clear all;
n = 2;     %II rząd
for k = 1:1:(n/2)
    a(k) = 2*cos((2*k-1)/(2*n)*pi);
    b(k) = 1;
    num = 1;
    den = [b(k),a(k),1]
    B(k) = tf(num,den);
end
b2 = B(1);
n = 3;     %III rząd
for k = 1:1:((n+1)/2)
    a(k) = 2*cos(((k-1)/n)*pi);
    b(k) = 1;
    b(1) = 0;
    a(1) = 1;
    num = 1;
    den = [b(k),a(k),1];
    B(k) = tf(num,den);
end
b3 = B(1)*B(2);
n = 4;     %IV rząd
for k = 1:1:(n/2)
    a(k) = 2*cos((2*k-1)/(2*n)*pi);
    b(k) = 1;
    num = 1;
    den = [b(k),a(k),1];
    B(k) = tf(num,den);
end
b4 = B(1)*B(2);
n = 5;     %V rząd
for k = 1:1:((n+1)/2)
    a(k) = 2*cos(((k-1)/n)*pi);
    b(k) = 1;
    b(1) = 0;
    a(1) = 1;
    num = 1;
    den = [b(k),a(k),1];
    B(k) = tf(num,den);
end
b5 = B(1)*B(2)*B(3);
n = 6;     %VI rząd
for k = 1:1:(n/2)
    a(k) = 2*cos((2*k-1)/(2*n)*pi);
    b(k) = 1;
    num = 1;
    den = [b(k),a(k),1]
    B(k) = tf(num,den);
end
b6 = B(1)*B(2)*B(3);
n = 7;     %VII rząd
for k = 1:1:((n+1)/2)
    a(k) = 2*cos(((k-1)/n)*pi);
    b(k) = 1;
    b(1) = 0;
    a(1) = 1;
    num = 1;
    den = [b(k),a(k),1];
    B(k) = tf(num,den);
end
b7 = B(1)*B(2)*B(3)*B(4);
n = 8;     %VIII rząd
for k = 1:1:(n/2)
    a(k) = 2*cos((2*k-1)/(2*n)*pi);
    b(k) = 1;
    num = 1;
    den = [b(k),a(k),1]
    B(k) = tf(num,den);
end
b8 = B(1)*B(2)*B(3)*B(4);
bode(b2,b3,b4,b5,b6,b7,b8)
grid on;