Routh Hurwitz Criteria using user defined Function

Routh matrix can be derived using the user defined function'routh'

RA=ROUTH(R,EPSILON) returns the symbolic Routh array RA for polynomial. The following special cases are considered:
1) If the first element of a row becomes zero
OR
2) If one encounters a row full of zeros.

>>syms ep
>>a=routh([1 1 2 2 3 5],ep)

The above given case is for encountering a zero in the first column.