num1 = [1 2]; num2 = [1 2]; den1 = [1 20 0 0]; den2 = [1 6 25]; num = conv(num1,num2); den = conv(den1,den2); sys = tf(num, den); bode(sys); nyquist(sys);