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 / rem.m < prev    next >
Text File  |  1996-09-28  |  1KB  |  47 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 = rem (x, y)
  16.  
  17. # usage: rem (x, y)
  18. #
  19. # Return remainder (x, y).
  20.  
  21.   if (nargin != 2)
  22.     usage ("rem (x, y)");
  23.   endif
  24.  
  25.   if (any (size (x) != size (y)) && ! (is_scalar (x) || is_scalar (y)))
  26.     error ("rem: argument sizes must agree");
  27.   endif
  28.  
  29. # Matlab allows complex arguments, but as far as I can tell, that's a
  30. # bunch of hooey.
  31.  
  32.   if (any (any (imag (x))) || any (any (imag (y))))
  33.     error ("rem: complex arguments are not allowed");
  34.   endif
  35.  
  36.   if (nargin == 2)
  37.     retval = x - y .* fix (x ./ y);
  38.   else
  39.     error ("usage: rem (x, y)");
  40.   endif
  41.  
  42. endfunction
  43.