clear; close all; s = tf('s'); Ns = 1:9; for idx = 1:length(Ns) n = Ns(idx); D = 1; for jdx=1:n j = jdx; sigma_j = -sin((2*j-1)/(2*n)*pi); omega_j = cos((2*j-1)/(2*n)*pi); sj = sigma_j + i*omega_j; D = D*(s-sj); end coef(D) end function n …