new file mode 100644
--- /dev/null
+++ b/liboctave/oct-rand.h
@@ -0,0 +1,67 @@
+/*
+
+Copyright (C) 2003 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.
+
+*/
+
+#if !defined (octave_rand_h)
+#define octave_rand_h 1
+
+#include <string>
+
+#include "dMatrix.h"
+
+struct
+octave_rand
+{
+  // Return the current seed.
+  static double seed (void);
+
+  // Set the seed.
+  static void seed (double s);
+
+  // Return the current distribution.
+  static std::string distribution (void);
+
+  // Set the current distribution.  May be either "uniform" (the
+  // default) or "normal".
+  static void distribution (const std::string& d);
+
+  static void uniform_distribution (void);
+
+  static void normal_distribution (void);
+
+  // Return the next number from the sequence.
+  static double scalar (void);
+
+  // Return a matrix of numbers from the sequence, filled in column
+  // major order.
+  static Matrix matrix (int r, int c);
+
+  // Return an array of numbers from the sequence.
+  static Array<double> vector (int n);
+};
+
+#endif
+
+/*
+;;; Local Variables: ***
+;;; mode: C++ ***
+;;; End: ***
+*/