home *** CD-ROM | disk | FTP | other *** search
/ Geek Gadgets 1 / ADE-1.bin / ade-dist / octave-1.1.1p1-src.tgz / tar.out / fsf / octave / scripts / general / is_symmetric.m < prev    next >
Text File  |  1996-09-28  |  1KB  |  44 lines

  1. # Copyright (C) 1993, 1994, 1995 John W. Eaton
  2. # This file is part of Octave.
  3. # Octave is free software; you can redistribute it and/or modify it
  4. # under the terms of the GNU General Public License as published by the
  5. # Free Software Foundation; either version 2, or (at your option) any
  6. # later version.
  7. # Octave is distributed in the hope that it will be useful, but WITHOUT
  8. # ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
  9. # FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
  10. # for more details.
  11. # You should have received a copy of the GNU General Public License
  12. # along with Octave; see the file COPYING.  If not, write to the Free
  13. # Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA.
  14.  
  15. function retval = is_symmetric (x,tol)
  16.  
  17. # Usage: is_symmetric (x {,tol})
  18. #
  19. # If x is symmetric, return the dimension of x, otherwise, return 0.
  20. #
  21. # See also: size, rows, columns, length, is_matrix, is_scalar, 
  22. # is_square, is_vector
  23.  
  24. # Written by A. S. Hodel (scotte@eng.auburn.edu) August 1993.
  25.  
  26.   if (nargin == 1 || nargin == 2)
  27.     if ((retval = is_square (x)))
  28.       if (nargin == 1)
  29.     tol = eps;
  30.       endif
  31.       if (norm (x - x') / norm(x) > tol)
  32.         retval = 0;
  33.       endif
  34.     endif
  35.   else
  36.     usage ("is_symmetric (x {,tol})");
  37.   endif
  38.  
  39. endfunction
  40.