diff --git a/eo/src/es/eoReal.h b/eo/src/es/eoReal.h
new file mode 100644
index 00000000..1ca1247f
--- /dev/null
+++ b/eo/src/es/eoReal.h
@@ -0,0 +1,58 @@
+/*
+   eoReal.h
+// (c) Marc Schoenauer, Maarten Keijzer and GeNeura Team, 2000
+
+    This library is free software; you can redistribute it and/or
+    modify it under the terms of the GNU Lesser General Public
+    License as published by the Free Software Foundation; either
+    version 2 of the License, or (at your option) any later version.
+
+    This library 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
+    Lesser General Public License for more details.
+
+    You should have received a copy of the GNU Lesser General Public
+    License along with this library; if not, write to the Free Software
+    Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
+
+    Contact: Marc.Schoenauer@polytechnique.fr
+             todos@geneura.ugr.es, http://geneura.ugr.es 
+             mkeijzer@dhi.dk
+*/
+
+#ifndef eoReal_h
+#define eoReal_h
+
+//-----------------------------------------------------------------------------
+
+#include <iostream>    // ostream, istream
+#include <string>      // string
+
+#include <eoFixedLength.h>
+
+/** eoReal: implementation of simple real-valued chromosome.
+  * based on eoFixedLength class
+*/
+template <class F> class eoReal: public eoFixedLength<F, double>
+{
+ public:
+
+  /**
+   * (Default) Constructor.
+   * @param size Size of the vector
+   */
+  eoReal(unsigned size = 0, double value = 0.0): 
+    eoFixedLength<F, double>(size, value) {}
+      
+  /// My class name.
+  string className() const 
+    { 
+      return "eoReal"; 
+    }
+  
+};
+
+//-----------------------------------------------------------------------------
+
+#endif //eoReal_h
diff --git a/eo/src/es/eoRealOp.h b/eo/src/es/eoRealOp.h
new file mode 100644
index 00000000..a7b66173
--- /dev/null
+++ b/eo/src/es/eoRealOp.h
@@ -0,0 +1,272 @@
+// -*- mode: c++; c-indent-level: 4; c++-member-init-indent: 8; comment-column: 35; -*-
+
+//-----------------------------------------------------------------------------
+// eoRealOp.h
+// (c) EEAAX 2000 - Maarten Keijzer 2000
+/* 
+    This library is free software; you can redistribute it and/or
+    modify it under the terms of the GNU Lesser General Public
+    License as published by the Free Software Foundation; either
+    version 2 of the License, or (at your option) any later version.
+
+    This library 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
+    Lesser General Public License for more details.
+
+    You should have received a copy of the GNU Lesser General Public
+    License along with this library; if not, write to the Free Software
+    Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
+
+    Contact: Marc.Schoenauer@polytechnique.fr
+             mak@dhi.dk
+ */
+//-----------------------------------------------------------------------------
+
+#ifndef eoRealOp_h
+#define eoRealOp_h
+
+//-----------------------------------------------------------------------------
+
+#include <algorithm>    // swap_ranges
+#include <utils/eoRNG.h>
+#include <es/eoReal.h>
+
+//-----------------------------------------------------------------------------
+
+/** eoUniformMutation --> changes all values of the vector 
+                          by uniform choice with range epsilon  
+                          with probability p_change per variable
+\class eoUniformMutation eoRealOp.h Tutorial/eoRealOp.h
+\ingroup parameteric
+*/
+
+template<class Chrom> class eoUniformMutation: public eoMonOp<Chrom>
+{
+ public:
+  /**
+   * (Default) Constructor.
+   * @param _epsilon the range for uniform nutation
+   * @param _p_change the probability to change a given coordinate
+   */
+  eoUniformMutation(const double& _epsilon, const double& _p_change = 1.0): 
+    epsilon(_epsilon), p_change(_p_change) {}
+
+  /// The class name.
+  string className() const { return "eoUniformMutation"; }
+  
+  /**
+   * Do it!
+   * @param chrom The cromosome undergoing the mutation
+   */
+  void operator()(Chrom& chrom) 
+    {
+      chrom.invalidate();
+      for (unsigned lieu=0; lieu<chrom.size(); lieu++) 
+	{
+	  if (rng.flip(p_change)) 
+	    {
+	      chrom[lieu] += 2*epsilon*rng.uniform()-epsilon;
+	    }
+	}
+    }
+  
+private:
+  double epsilon;
+  double p_change;
+};
+
+/** eoDetUniformMutation --> changes exactly k values of the vector 
+                          by uniform choice with range epsilon  
+\class eoDetUniformMutation eoRealOp.h Tutorial/eoRealOp.h
+\ingroup parameteric
+*/
+
+template<class Chrom> class eoDetUniformMutation: public eoMonOp<Chrom>
+{
+ public:
+  /**
+   * (Default) Constructor.
+   * @param _epsilon the range for uniform nutation
+   * @param number of coordinate to modify
+   */
+  eoDetUniformMutation(const double& _epsilon, const unsigned& _no = 1): 
+    epsilon(_epsilon), no(_no) {}
+
+  /// The class name.
+  string className() const { return "eoDetUniformMutation"; }
+  
+  /**
+   * Do it!
+   * @param chrom The cromosome undergoing the mutation
+   */
+  void operator()(Chrom& chrom) 
+    {
+      chrom.invalidate();
+      for (unsigned i=0; i<no; i++)
+	{
+	  unsigned lieu = rng.random(chrom.size());
+	  // actually, we should test that we don't re-modify same variable!
+	  chrom[lieu] += 2*epsilon*rng.uniform()-epsilon;
+	}
+    }
+
+private:
+  double epsilon;
+  unsigned no;
+};
+
+
+// two arithmetical crossovers
+
+/** eoSegmentCrossover --> uniform choice in segment 
+                 == arithmetical with same value along all coordinates
+\class eoSegmentCrossover eoRealOp.h Tutorial/eoRealOp.h
+\ingroup parameteric
+*/
+
+template<class Chrom> class eoSegmentCrossover: public eoQuadraticOp<Chrom>
+{
+ public:
+  /**
+   * (Default) Constructor.
+   * @param _alpha the amount of exploration OUTSIDE the parents 
+   *               as in BLX-alpha notation (Eshelman and Schaffer)
+   *               0 == contractive application
+   */
+  eoSegmentCrossover(const double& _alpha = 0.0) : 
+    alpha(_alpha), range(1+2*alpha) {}
+
+  /// The class name.
+  string className() const { return "eoSegmentCrossover"; }
+
+  /**
+   * segment crossover - modifies both parents
+   * @param chrom1 The first parent
+   * @param chrom2 The first parent
+   */
+  void operator()(Chrom& chrom1, Chrom& chrom2) 
+    {
+      unsigned i;
+      double r1, r2, fact;
+      fact = rng.uniform(range);	   // in [0,range)
+      for (i=0; i<chrom1.size(); i++)
+	{
+	  r1=chrom1[i];
+	  r2=chrom2[i];
+	  chrom1[i] = fact * r1 + (1-fact) * r2;
+	  chrom2[i] = (1-fact) * r1 + fact * r2;
+	}
+      chrom1.invalidate();	   // shoudl test if fact was 0 or 1 :-)))
+      chrom2.invalidate();
+    }
+
+protected:
+  double alpha;
+  double range;			   // == 1+2*alpha
+};
+  
+/** eoArithmeticCrossover --> uniform choice in hypercube  
+                 == arithmetical with different values for each coordinate
+\class eoArithmeticCrossover eoRealOp.h Tutorial/eoRealOp.h
+\ingroup parameteric
+*/
+
+template<class Chrom> class eoArithmeticCrossover: public eoQuadraticOp<Chrom>
+{
+ public:
+  /**
+   * (Default) Constructor.
+   * @param _alpha the amount of exploration OUTSIDE the parents 
+   *               as in BLX-alpha notation (Eshelman and Schaffer)
+   *               0 == contractive application
+   */
+  eoArithmeticCrossover(const double& _alpha = 0.0): 
+    alpha(_alpha), range(1+2*alpha) {}
+
+  /// The class name.
+  string className() const { return "eoArithmeticCrossover"; }
+
+  /**
+   * arithmetical crossover - modifies both parents
+   * @param chrom1 The first parent
+   * @param chrom2 The first parent
+   */
+  void operator()(Chrom& chrom1, Chrom& chrom2) 
+    {
+      unsigned i;
+      double r1, r2, fact;
+      for (i=0; i<chrom1.size(); i++)
+	{
+	  r1=chrom1[i];
+	  r2=chrom2[i];
+	  fact = rng.uniform(range);	   // in [0,range)
+	  chrom1[i] = fact * r1 + (1-fact) * r2;
+	  chrom2[i] = (1-fact) * r1 + fact * r2;
+	}
+    }
+
+protected:
+  double alpha;
+  double range;			   // == 1+2*alpha
+};
+  
+
+/** eoRealUxOver --> Uniform crossover, also termed intermediate crossover
+\class eoRealUxOver eoRealOp.h Tutorial/eoRealOp.h
+\ingroup parameteric
+*/
+
+template<class Chrom> class eoRealUxOver: public eoQuadraticOp<Chrom>
+{
+ public:
+  /**
+   * (Default) Constructor.
+   * @param _preference bias in the choice (usually, no bias == 0.5)
+   */
+  eoRealUxOver(const float& _preference = 0.5): preference(_preference)
+    { 
+      if ( (_preference <= 0.0) || (_preference >= 1.0) )
+	runtime_error("UxOver --> invalid preference");
+    }
+
+  /// The class name.
+  string className() const { return "eoRealUxOver"; }
+
+  /**
+   * Uniform crossover for real vectors
+   * @param chrom1 The first parent
+   * @param chrom2 The second parent
+   *    @runtime_error if sizes don't match
+   */
+  void operator()(Chrom& chrom1, Chrom& chrom2) 
+    {
+      if ( chrom1.size() != chrom2.size()) 
+	    runtime_error("UxOver --> chromosomes sizes don't match" ); 
+      bool changed = false;
+      for (unsigned int i=0; i<chrom1.size(); i++)
+	{
+	  if (rng.flip(preference))
+	    if (chrom1[i] == chrom2[i])
+	      {
+		double tmp = chrom1[i];
+	      chrom1[i]=chrom2[i];
+	      chrom2[i] = tmp;
+	      changed = true;
+	    }
+	}
+      if (changed)
+	  {
+	    chrom1.invalidate();
+	    chrom2.invalidate();
+	  }
+    }
+    private:
+      float preference;
+};
+  
+
+//-----------------------------------------------------------------------------
+//@}
+#endif eoRealOp_h
+