basselFilter

a = 1
s = tf('s');
for n = 2:1:8
    i = n-1;
    a = a * (i*2+1);
    b = (i*2+1)/s
    for k = 1:1:n-1
        K_TEMP = (n-k)*2 - 1
        b = 1/b + (K_TEMP/s)
    end

    num = cell2mat(b.numerator)
    den = cell2mat(b.denominator)
    bassel = den + num

    figure(1)
    bode(a, bassel)
    grid on;
    hold on;

    figure(2)
    step(a,bassel)
    grid on;
    hold on;
end