--- /dev/null
+## Copyright (C) 2007-2012 David Bateman
+##
+## 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 3 of the License, 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, see
+## <http://www.gnu.org/licenses/>.
+
+## -*- texinfo -*-
+## @deftypefn {Function File} {} scatter (@var{x}, @var{y})
+## @deftypefnx {Function File} {} scatter (@var{x}, @var{y}, @var{s})
+## @deftypefnx {Function File} {} scatter (@var{x}, @var{y}, @var{c})
+## @deftypefnx {Function File} {} scatter (@var{x}, @var{y}, @var{s}, @var{c})
+## @deftypefnx {Function File} {} scatter (@var{x}, @var{y}, @var{s}, @var{c}, @var{style})
+## @deftypefnx {Function File} {} scatter (@var{x}, @var{y}, @var{s}, @var{c}, @var{prop}, @var{val})
+## @deftypefnx {Function File} {} scatter (@dots{}, "filled")
+## @deftypefnx {Function File} {} scatter (@var{h}, @dots{})
+## @deftypefnx {Function File} {@var{h} =} scatter (@dots{})
+##
+## Plot a scatter plot of the data. A marker is plotted at each point
+## defined by the points in the vectors @var{x} and @var{y}. The size of
+## the markers used is determined by the @var{s}, which can be a scalar,
+## a vector of the same length of @var{x} and @var{y}. If @var{s} is not
+## given or is an empty matrix, then the default value of 8 points is used.
+##
+## The color of the markers is determined by @var{c}, which can be a string
+## defining a fixed color; a 3-element vector giving the red, green,and blue
+## components of the color; a vector of the same length as @var{x} that gives
+## a scaled index into the current colormap; or a @var{n}-by-3 matrix defining
+## the colors of each of the markers individually.
+##
+## The marker to use can be changed with the @var{style} argument, that is a
+## string defining a marker in the same manner as the @code{plot} command.
+## If the argument @code{"filled"} is given then the markers as filled. All
+## additional arguments are passed to the underlying patch command.
+##
+## The optional return value @var{h} provides a handle to the patch object
+##
+## @example
+## @group
+## x = randn (100, 1);
+## y = randn (100, 1);
+## scatter (x, y, [], sqrt(x.^2 + y.^2));
+## @end group
+## @end example
+##
+## @seealso{plot, patch, scatter3}
+## @end deftypefn
+
+function retval = scatter (varargin)
+
+ [h, varargin, nargin] = __plt_get_axis_arg__ ("scatter", varargin{:});
+
+ if (nargin < 2)
+ print_usage ();
+ else
+ oldh = gca ();
+ unwind_protect
+ axes (h);
+ newplot ();
+ tmp = __scatter__ (h, 2, "scatter", varargin{:});
+ unwind_protect_cleanup
+ axes (oldh);
+ end_unwind_protect
+ endif
+
+ if (nargout > 0)
+ retval = tmp;
+ endif
+
+endfunction
+
+
+%!demo
+%! clf
+%! x = randn (100, 1);
+%! y = randn (100, 1);
+%! scatter (x, y, "r");
+
+%!demo
+%! clf
+%! x = randn (100, 1);
+%! y = randn (100, 1);
+%! scatter (x, y, [], sqrt (x.^2 + y.^2));
+
+%!demo
+%! clf
+%! rand_10x1_data1 = [0.171577, 0.404796, 0.025469, 0.335309, 0.047814, 0.898480, 0.639599, 0.700247, 0.497798, 0.737940];
+%! rand_10x1_data2 = [0.75495, 0.83991, 0.80850, 0.73603, 0.19360, 0.72573, 0.69371, 0.74388, 0.13837, 0.54143];
+%! x = rand_10x1_data1;
+%! y = rand_10x1_data2;
+%! s = 10 - 10*log (x.^2 + y.^2);
+%! h = scatter (x, y, s, s, "s", "filled");
+
+%!demo
+%! clf
+%! rand_10x1_data3 = [0.42262, 0.51623, 0.65992, 0.14999, 0.68385, 0.55929, 0.52251, 0.92204, 0.19762, 0.93726];
+%! rand_10x1_data4 = [0.020207, 0.527193, 0.443472, 0.061683, 0.370277, 0.947349, 0.249591, 0.666304, 0.134247, 0.920356];
+%! x = rand_10x1_data3;
+%! y = rand_10x1_data4;
+%! s = 10 - 10*log (x.^2 + y.^2);
+%! h = scatter (x, y, [], "r", "s", "filled");
+
+%!demo
+%! clf
+%! rand_10x1_data5 = [0.777753, 0.093848, 0.183162, 0.399499, 0.337997, 0.686724, 0.073906, 0.651808, 0.869273, 0.137949];
+%! rand_10x1_data6 = [0.37460, 0.25027, 0.19510, 0.51182, 0.54704, 0.56087, 0.24853, 0.75443, 0.42712, 0.44273];
+%! x = rand_10x1_data5;
+%! y = rand_10x1_data6;
+%! s = 10 - 10*log (x.^2 + y.^2);
+%! h = scatter (x, y, [], "r", "s");
+
+%!demo
+%! k = 1;
+%! clf
+%! for m = [1, 3]
+%! for n = [101, 50, 1]
+%! x = rand (n, 1);
+%! y = rand (n, 1);
+%! if (m > 1)
+%! str = "Three Colors";
+%! idx = ceil (rand (n, 1) * 3);
+%! colors = eye(3);
+%! colors = colors(idx, :);
+%! else
+%! str = "Random Colors";
+%! colors = rand (n, m);
+%! endif
+%! if (n == 1)
+%! str = sprintf ("%s: 1 point", str);
+%! elseif (n < 100)
+%! str = sprintf ("%s: < 100 points", str);
+%! else
+%! str = sprintf ("%s: > 100 points", str);
+%! endif
+%! subplot (2, 3, k)
+%! k = k + 1;
+%! scatter (x, y, 15, colors, "filled")
+%! axis ([0 1 0 1])
+%! title (str)
+%! endfor
+%! endfor
+
+%!demo
+%! k = 1;
+%! clf
+%! for m = [1, 3]
+%! for n = [101, 50, 1]
+%! x = rand (n, 1);
+%! y = rand (n, 1);
+%! if (m > 1)
+%! str = "Three Colors";
+%! idx = ceil (rand (n, 1) * 3);
+%! colors = eye(3);
+%! colors = colors(idx, :);
+%! else
+%! str = "Random Colors";
+%! colors = rand (n, m);
+%! endif
+%! if (n == 1)
+%! str = sprintf ("%s: 1 point", str);
+%! elseif (n < 100)
+%! str = sprintf ("%s: < 100 points", str);
+%! else
+%! str = sprintf ("%s: > 100 points", str);
+%! endif
+%! subplot (2, 3, k)
+%! k = k + 1;
+%! scatter (x, y, 15, colors)
+%! axis ([0 1 0 1])
+%! title (str)
+%! endfor
+%! endfor