Mercurial > hg > octave-lyh
diff scripts/polynomial/poly.m @ 1025:f558749713f1
[project @ 1995-01-11 20:52:10 by jwe]
author | jwe |
---|---|
date | Wed, 11 Jan 1995 20:52:10 +0000 |
parents | 3470f1e25a79 |
children | 611d403c7f3d |
line wrap: on
line diff
--- a/scripts/polynomial/poly.m +++ b/scripts/polynomial/poly.m @@ -1,13 +1,37 @@ +# Copyright (C) 1995 John W. Eaton +# +# This file is part of Octave. +# +# Octave is free software; you can redistribute it and/or modify it +# under the terms of the GNU General Public License as published by the +# Free Software Foundation; either version 2, or (at your option) any +# later version. +# +# Octave is distributed in the hope that it will be useful, but WITHOUT +# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or +# FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License +# for more details. +# +# You should have received a copy of the GNU General Public License +# along with Octave; see the file COPYING. If not, write to the Free +# Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA. + function y = poly (x) +# usage: poly (x) +# # If A is a square n-by-n matrix, poly (A) is the row vector of # the coefficients of det (z * eye(n) - A), the characteristic # polynomial of A. +# # If x is a vector, poly (x) is a vector of coefficients of the # polynomial whose roots are the elements of x. -# Written by KH (Kurt.Hornik@neuro.tuwien.ac.at) on Dec 24, 1993 -# Copyright Dept of Probability Theory and Statistics TU Wien +# Written by KH (Kurt.Hornik@neuro.tuwien.ac.at) Dec 24, 1993. + + if (nargin != 1) + usage ("poly (x)"); + endif m = min (size (x)); n = max (size (x)); @@ -18,7 +42,7 @@ elseif (m == n) v = eig (x); else - usage ("poly(x), where x is a vector or a square matrix"); + usage ("poly (x), where x is a vector or a square matrix"); endif y = [1, zeros (1, n)];