Routh's array in symbolic way

Routh's array of symbolic polynomials.

[M,L,sol] = routh_hurwitz(P,par)

This function gives the Routh's array from a numerical or symbolic polynomial and includes two special cases: (1) the first element of the row is zero; (b) a row of zeros.

P - Numerical or symbolic array of coeficients. In the case of symbolic variables it is necesarry to define them as: syms a b c ...

par - Parameter in relation to which stability is determined (optional)

M - Routh's array without any simplification (e.g., with epsilon notation)

L - First column of Routh's array with simplification (e.g., using the limit when epsilon tends to zero)

sol - Stability condition