ucp

 UCP(C) tests if the polynomial C is a Unit-Circle-Polynomial.
    It tests if all roots lie inside the unit circle like
       B=ucp(C) does the same as
    B=all(abs(roots(C))<1) but much faster.
    The Algorithm is based on the Jury-Scheme.
    C are the elements of the Polynomial
    C(1)*X^N + ... + C(N)*X + C(N+1).
 
 REFERENCES:
  O. Foellinger "Lineare Abtastsysteme", Oldenburg Verlag, Muenchen, 1986.
  F. Gausch "Systemtechnik", Textbook, University of Technology Graz, 1993.