# Copyright (C) 2000 by Peter J. Gawthrop

###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.59  2004/06/20 13:53:06  gawthrop
## Initialise for numpar as well
##
## Revision 1.58  2004/06/20 13:34:11  gawthrop
## Create vectors using zeros(N,1) - avoids new octave
## default to row vector - more reps now included
##
## Revision 1.57  2004/06/20 10:18:00  gawthrop
## Explicitly size vectors - avoid default row vector bug
##

    Npar=`wc -l $system\_sympar.txt | gawk '{print $1}'`
    }

zero_matrices()
{
## Set matrices to zero
echo
 ## echo '## Set matrices to zero'
for name in $zeromatrices; do
    case $name in
	a)
	    N=$Nx; M=$Nx
	    ;;
	b)
	    N=$Nx; M=$Nu