diff --git a/examples/integrate_times.cpp b/examples/integrate_times.cpp
new file mode 100644
index 0000000..9e6b581
--- /dev/null
+++ b/examples/integrate_times.cpp
@@ -0,0 +1,54 @@
+/* Boost libs/numeric/odeint/examples/integrate_times.cpp
+
+ Copyright 2009-2014 Karsten Ahnert
+ Copyright 2009-2014 Mario Mulansky
+
+ example for the use of integrate_times
+
+ Distributed under the Boost Software License, Version 1.0.
+(See accompanying file LICENSE_1_0.txt or
+ copy at http://www.boost.org/LICENSE_1_0.txt)
+ */
+
+
+#include <iostream>
+#include <boost/numeric/odeint.hpp>
+
+using namespace std;
+using namespace boost::numeric::odeint;
+
+
+/*
+ * simple 1D ODE
+ */
+
+void rhs( const double x , double &dxdt , const double t )
+{
+    dxdt = 3.0/(2.0*t*t) + x/(2.0*t);
+}
+
+void write_cout( const double &x , const double t )
+{
+    cout << t << '\t' << x << endl;
+}
+
+// state_type = double
+typedef runge_kutta_dopri5< double > stepper_type;
+
+const double dt = 0.1;
+
+int main()
+{
+    // create a vector with observation time points
+    std::vector<double> times( 91 );
+    for( size_t i=0 ; i<times.size() ; ++i )
+        times[i] = 1.0 + dt*i;
+
+    double x = 0.0; //initial value x(1) = 0
+    // we can provide the observation time as a boost range (i.e. the vector)
+    integrate_times( make_controlled( 1E-12 , 1E-12 , stepper_type() ) , rhs ,
+                     x , times , dt , write_cout );
+    // or as two iterators
+    //integrate_times( make_controlled( 1E-12 , 1E-12 , stepper_type() ) , rhs ,
+    //                   x , times.begin() , times.end() , dt , write_cout );
+}