diff --git a/tests/Makefile.am b/tests/Makefile.am
index ecf041ea48c51194949c9743c3db4b260b0859fe..d2d5b09d5fabf5a9b22001f5899a2fb850712902 100644
--- a/tests/Makefile.am
+++ b/tests/Makefile.am
@@ -8,6 +8,7 @@ OCTAVE_MODS = \
 	example2.mod \
 	example1_use_dll.mod \
 	t_sgu_ex1.mod \
+	osr_example.mod \
 	ramsey.mod \
 	ramst_initval_file.mod \
 	example1_varexo_det.mod \
diff --git a/tests/osr_example.mod b/tests/osr_example.mod
new file mode 100644
index 0000000000000000000000000000000000000000..538d142f9de8779820b416a16912663cd587b941
--- /dev/null
+++ b/tests/osr_example.mod
@@ -0,0 +1,41 @@
+// Example of optimal simple rule
+
+var y inflation r;
+varexo y_ inf_;
+
+parameters delta sigma alpha kappa gammarr gammax0 gammac0 gamma_y_ gamma_inf_;
+
+delta =  0.44;
+kappa =  0.18;
+alpha =  0.48;
+sigma = -0.06;
+
+
+model(linear);
+y  = delta * y(-1)  + (1-delta)*y(+1)+sigma *(r - inflation(+1)) + y_; 
+inflation  =   alpha * inflation(-1) + (1-alpha) * inflation(+1) + kappa*y + inf_;
+r = gammax0*y(-1)+gammac0*inflation(-1)+gamma_y_*y_+gamma_inf_*inf_;
+end;
+
+shocks;
+var y_;
+stderr 0.63;
+var inf_;
+stderr 0.4;
+end;
+
+
+optim_weights;
+inflation 1;
+y 1;
+end;
+
+osr_params gammax0 gammac0 gamma_y_ gamma_inf_;
+
+gammarr = 0;
+gammax0 = 0.2;
+gammac0 = 1.5;
+gamma_y_ = 8;
+gamma_inf_ = 3;
+
+osr;