Mercurial > hg > octave-lyh
diff scripts/polynomial/compan.m @ 904:3470f1e25a79
[project @ 1994-11-09 21:22:15 by jwe]
| author | jwe |
|---|---|
| date | Wed, 09 Nov 1994 21:22:15 +0000 |
| parents | e79ff1f4df3c |
| children | f558749713f1 |
line wrap: on
line diff
--- a/scripts/polynomial/compan.m +++ b/scripts/polynomial/compan.m @@ -1,27 +1,28 @@ function A = compan(c) -#compan (c) -#Compute the companion matrix corresponding to polynomial vector c. + +# compan (c) +# Compute the companion matrix corresponding to polynomial vector c. # -#In octave a polynomial is represented by it's coefficients (arranged -#in descending order). For example a vector c of length n+1 corresponds -#to the following nth order polynomial +# In octave a polynomial is represented by it's coefficients (arranged +# in descending order). For example a vector c of length n+1 corresponds +# to the following nth order polynomial # -# p(x) = c(1) x^n + ... + c(n) x + c(n+1). +# p(x) = c(1) x^n + ... + c(n) x + c(n+1). # -#The corresponding companion matrix is -# _ _ -# | -c(2)/c(1) -c(3)/c(1) ... -c(n)/c(1) -c(n+1)/c(1) | -# | 1 0 ... 0 0 | -# | 0 1 ... 0 0 | -# A = | . . . . . | -# | . . . . . | -# | . . . . . | -# |_ 0 0 ... 1 0 _| +# The corresponding companion matrix is +# _ _ +# | -c(2)/c(1) -c(3)/c(1) ... -c(n)/c(1) -c(n+1)/c(1) | +# | 1 0 ... 0 0 | +# | 0 1 ... 0 0 | +# A = | . . . . . | +# | . . . . . | +# | . . . . . | +# |_ 0 0 ... 1 0 _| # -#The eigenvalues of the companion matrix are equal to the roots of the -#polynomial. +# The eigenvalues of the companion matrix are equal to the roots of the +# polynomial. # -#SEE ALSO: poly, roots, residue, conv, deconv, polyval, polyderiv, polyinteg +# SEE ALSO: poly, roots, residue, conv, deconv, polyval, polyderiv, polyinteg # Author: # Tony Richardson @@ -29,7 +30,7 @@ # June 1994 if(nargin != 1) - error("usage: compan(vector)"); + usage ("compan(vector)"); endif if(!is_vector(c))