Scilab Function
Last update : 23/10/2007
roots - racines d'un polynôme
Calling Sequence
-
[x]=roots(p)
Parameters
-
p
: polynôme à coefficients réels ou complexes ou vecteur des coefficients du
polynŽôme ordonnés par puissance décroissante (compatibilité avec
Matlab).
Description
x=roots(p)
renvoie dans le vecteur complexe
x
les
racines du polynôme
p
. Pour les polynŽômes à coéfficients
réels et de degré <=100, l'algorithme rapide RPOLY est utilisé. Dans
les autres cas les racines sont calculées comme valeurs propres de la
matrice compagnion du polynŽôme. Pour forcer ce dernier algorithme dans
tous les cas utilisez
x=roots(p,'e')
.
Examples
p=poly([0,10,1+%i,1-%i],'x');
roots(p)
A=rand(3,3);roots(poly(A,'x')) // comparaison via le polynôme caractéristique
spec(A)
See Also
poly
,
spec
,
companion
,
Authors
-
Serge Steer (INRIA)
Bibliography
La routine RPOLY est decrite dans ACM TOMS 1 (1975) 178-189
Used Function
Le code source de rpoly.f peut être trouvé dans le repertoire
routines/control de la distribution source de Scilab. Le calcul des
valeurs propores est effectué en utilisant les routines DGEEV et ZGEEV de
LAPACK.