Jury's array in symbolic way

[M, L] = jury(P,N)This function gives the Jury's array from a numerical or SYMBOLIC polynomial and includes the two special cases: (1) the first element of the second row of a block is zero; (b) a row of zeros (when there are roots on unit circle...