Scilab Function
Last update : 14/2/2006

roots - roots of polynomials

Calling Sequence

[x]=roots(p)
[x]=roots(p,'e')

Parameters

Description

x=roots(p) returns in the complex vector x the roots of the polynomial p . For real polynomials of degree <=100 the fast RPOLY algorithm is used. In the other cases the roots are computed as the eigenvalues of the associated companion matrix. Use x=roots(p,'e') to force this algorithm in any cases.

Examples


p=poly([0,10,1+%i,1-%i],'x');
roots(p)
A=rand(3,3);roots(poly(A,'x'))    // Evals by characteristic polynomial
spec(A) 
 
  

See Also

poly ,   spec ,   companion ,  

Authors

Serge Steer (INRIA)

Bibliography

The RPOLY algorithm is described in ACM TOMS 1 (1975) 178-189

Used Function

The rpoly.f source codes can be found in the directory routines/control of a Scilab source distribution. The eigenvalue computation is perfomed using DGEEV and ZGEEV LAPACK codes.