I need a matlab code that uses the Raible method to prove that the roots of a discrete trasnfer function are inside the unit circle. This is a simplification of Jurys stability method. The matlab code should ask the user to input the coeficients of $a_n$ and $b_n,i$t must show the Raible table that shows the values of $c_n$ and determine if the trasnfer function is stable $or$ not.

A picture below shows the example $1$ transfer function and how its Raible table should look like:

Ejemplo $1$: $H(z)=\frac{z^3+4z^2+z-1}{2z^4+5\mathrm{.}41z^3+10\mathrm{.}24z^2+4\mathrm{.}83z+1\mathrm{.}414}$

For the example shown above since $C_0=0,we$ can conclude that the transfer function $is$ unstable.

Criteria:

$-$All roots of C$($z$)$ have a magnitud of less than $1($inside the unit circle$)$ if the n coefficients $\{a_0,b_0,c_0,...h_0,f_0\}$ are all positive and greater than zero.

$-$If n coefficients from $a_0{f}_0$ are less than zero, then these $n$ roots are outside $of$ the unit circle.

$-a_0$ must always $be$ positive from the start, $if$ its not, then multiply above and below $by(-1).$