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)];