--- a/scripts/control/series.m
+++ b/scripts/control/series.m
@@ -1,43 +1,43 @@
-# Copyright (C) 1996 Auburn University.  All Rights Reserved.
-#
-# 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, 59 Temple Place, Suite 330, Boston, MA 02111 USA. 
+## Copyright (C) 1996 Auburn University.  All Rights Reserved.
+##
+## 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, 59 Temple Place, Suite 330, Boston, MA 02111 USA. 
  
 function [a,b,c,d] = series(a1,b1,c1,d1,a2,b2,c2,d2)
-# Forms the series connection of two systems.
-#
-# Superseded by sysmult.  Do not use this routine!
-# used internally in zp2ss
-#
-# Type of input: Transfer functions
-# Command:       [num,den]=series(num1,den1,num2,den2)
-# Forms the series representation of the two transfer functions.
-#
-# Type of input: State space systems
-# Command:       [a,b,c,d]=series(a1,b1,c1,d1,a2,b2,c2,d2)
-# Forms the series representation of the two state space system arguments.
-# The series connected system will have the inputs of system 1 and the 
-# outputs of system 2.
-#
-# Type of input: system data structure
-# Command:       syst=series(syst1,syst2)
-# Forms the series representation of the two mu system arguments.
-# Written by David Clem August 15, 1994
+## Forms the series connection of two systems.
+##
+## Superseded by sysmult.  Do not use this routine!
+## used internally in zp2ss
+##
+## Type of input: Transfer functions
+## Command:       [num,den]=series(num1,den1,num2,den2)
+## Forms the series representation of the two transfer functions.
+##
+## Type of input: State space systems
+## Command:       [a,b,c,d]=series(a1,b1,c1,d1,a2,b2,c2,d2)
+## Forms the series representation of the two state space system arguments.
+## The series connected system will have the inputs of system 1 and the 
+## outputs of system 2.
+##
+## Type of input: system data structure
+## Command:       syst=series(syst1,syst2)
+## Forms the series representation of the two mu system arguments.
+## Written by David Clem August 15, 1994
 
-# If two arguments input, take care of mu system case
+## If two arguments input, take care of mu system case
 
   warning("series is superseded by sysmult; use sysmult instead.")
 
@@ -49,7 +49,7 @@
     muflag = 1;
   endif
 
-# If four arguments input, put two transfer functions in series
+## If four arguments input, put two transfer functions in series
 
   if(nargin == 4)
     a = conv(a1,c1);	% was conv1
@@ -57,11 +57,11 @@
     c = 0;
     d = 0;
 
-# Find series combination of 2 state space systems
+## Find series combination of 2 state space systems
 
   elseif((nargin == 8)||(muflag == 1))
 
-# check matrix dimensions
+## check matrix dimensions
   
     [n1,m1,p1] = abcddim(a1,b1,c1,d1);
     [n2,m2,p2] = abcddim(a2,b2,c2,d2);
@@ -70,21 +70,21 @@
       error("Incorrect matrix dimensions");
     endif
 
-# check to make sure the number of outputs of system1 equals the number
-# of inputs of system2
+## check to make sure the number of outputs of system1 equals the number
+## of inputs of system2
 
    if(p1 ~= m2)
      error("System 1 output / System 2 input connection sizes do not match");
    endif
 
-# put the two state space systems in series
+## put the two state space systems in series
 
     a = [a1, zeros(rows(a1),columns(a2));b2*c1, a2];
     b = [b1;b2*d1];
     c = [d2*c1, c2];
     d = [d2*d1];
 
-# take care of mu output
+## take care of mu output
 
     if(muflag == 1)
       a=ss2sys(a,b,c,d);