--- a/scripts/polynomial/roots.m
+++ b/scripts/polynomial/roots.m
@@ -1,12 +1,31 @@
+# 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 r = roots (v)
 
+# usage: roots (v)
+#
 # For a vector v with n components, return the roots of the
 # polynomial v(1) * z^(n-1) + ... + v(n-1) * z + v(n).
   
-# Written by KH (Kurt.Hornik@neuro.tuwien.ac.at) on Dec 24, 1993
-# Copyright Dept of Probability Theory and Statistics TU Wien
-  
-  [nr, nc] = size(v);
+# Written by KH (Kurt.Hornik@neuro.tuwien.ac.at) on Dec 24, 1993.
+
+  [nr, nc] = size (v);
   if (nr <= 1 && nc <= 1)
     r = [];
     return;
@@ -23,14 +42,14 @@
   f = find (v);
   m = max (size (f));
   if (m > 0)
-    v = v(f(1):f(m));
+    v = v (f(1):f(m));
     l = max (size (v));
     if (l > 1)
       A = diag (ones (1, l-2), -1);
-      A(1,:) = -v(2:l) ./ v(1);
+      A (1, :) = -v (2:l) ./ v (1);
       r = eig (A);    
       if (f(m) < n)
-	r = [r; zeros (n - f(m), 1)];
+	r = [r; zeros (n - f (m), 1)];
       endif
     else
       r = zeros (n - f(m), 1);