Index: lemon/random.h
===================================================================
--- lemon/random.h	(revision 92)
+++ lemon/random.h	(revision 102)
@@ -576,4 +576,24 @@
     }
 
+    /// \brief Seeding random sequence
+    ///
+    /// Seeding the random sequence. The current number type will be
+    /// converted to the architecture word type.
+    template <typename Number>
+    void seed(Number seed) { 
+      _random_bits::Initializer<Number, Word>::init(core, seed);
+    }
+
+    /// \brief Seeding random sequence
+    ///
+    /// Seeding the random sequence. The given range should contain
+    /// any number type and the numbers will be converted to the
+    /// architecture word type.
+    template <typename Iterator>
+    void seed(Iterator begin, Iterator end) { 
+      typedef typename std::iterator_traits<Iterator>::value_type Number;
+      _random_bits::Initializer<Number, Word>::init(core, begin, end);
+    }
+
     /// \brief Returns a random real number from the range [0, 1)
     ///
Index: test/random_test.cc
===================================================================
--- test/random_test.cc	(revision 92)
+++ test/random_test.cc	(revision 102)
@@ -25,4 +25,6 @@
 ///
 
+int seed_array[] = {1, 2};
+
 int main()
 {
@@ -35,3 +37,9 @@
   a=lemon::rnd.gamma(4.0,0);
   a=lemon::rnd.poisson(.5);
+
+  lemon::rnd.seed(100);
+  lemon::rnd.seed(seed_array, seed_array + 
+		  (sizeof(seed_array) / sizeof(seed_array[0])));
+
+  return 0;
 }