Mercurial > octave-nkf
view scripts/elfun/lcm.m @ 4839:c1cdf2a32cdb ss-2-1-57
[project @ 2004-03-12 19:38:50 by jwe]
author | jwe |
---|---|
date | Fri, 12 Mar 2004 19:38:50 +0000 |
parents | e0b7a493e5a8 |
children | 2cbc6f37b0c2 |
line wrap: on
line source
## Copyright (C) 1996, 1997 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, 59 Temple Place - Suite 330, Boston, MA ## 02111-1307, USA. ## -*- texinfo -*- ## @deftypefn {Mapping Function} {} lcm (@var{x}, @code{...}) ## Compute the least common multiple of the elements elements of @var{x}, or ## the list of all the arguments. For example, ## ## @example ## lcm (a1, ..., ak) ## @end example ## ## @noindent ## is the same as ## ## @example ## lcm ([a1, ..., ak]). ## @end example ## @end deftypefn ## @seealso{gcd, min, max, ceil, and floor} ## Author: KH <Kurt.Hornik@ci.tuwien.ac.at> ## Created: 16 September 1994 ## Adapted-By: jwe function l = lcm (a, varargin) if (nargin == 0) usage ("lcm (a, ...)"); endif if (nargin > 1) k = 1; for i = 2:nargin a = [a, varargin{k++}]; endfor endif if (round (a) != a) error ("lcm: all arguments must be integer"); endif if (any (a) == 0) l = 0; else a = abs (a); l = a (1); for k = 1:(length (a) - 1) l = l * a(k+1) / gcd (l, a(k+1)); endfor endif endfunction