diff --git a/examples/var/clean b/examples/var/clean
index 1ac7e3d6a..e0d1cd114 100755
--- a/examples/var/clean
+++ b/examples/var/clean
@@ -1,6 +1,7 @@
 #!/bin/sh
 
-rm -rf m1 m2 m3
+rm -rf m1 m2 m3 m4
 rm -f m1_dynamic.m  m1.log  m1.m  m1_results.mat  m1_set_auxiliary_variables.m  m1_static.m
 rm -f m2_dynamic.m  m2.log  m2.m  m2_results.mat  m2_set_auxiliary_variables.m  m2_static.m
 rm -f m3_dynamic.m  m3.log  m3.m  m3_results.mat  m3_set_auxiliary_variables.m  m3_static.m
+rm -f m4_dynamic.m  m4.log  m4.m  m4_results.mat  m4_set_auxiliary_variables.m  m4_static.m
diff --git a/examples/var/m4.mod b/examples/var/m4.mod
new file mode 100644
index 000000000..f92fc75ed
--- /dev/null
+++ b/examples/var/m4.mod
@@ -0,0 +1,29 @@
+var y x z;
+
+varexo ex ey ez;
+
+parameters a_x_0 a_x_1 a_x_2 b_z_0 b_z_1 b_z_2;
+
+a_x_0 =  .2;
+a_x_1 =  .9;
+a_x_2 = -.2;
+b_z_0 =  .3;
+b_z_1 =  .7;
+b_z_2 = -.4;
+
+var_model(model_name=toto, eqtags=['eq:x', 'eq:z', 'eq:y']);
+
+model;
+
+[name='eq:y']
+y = y(-1) + ey ;
+
+[name='eq:x']
+diff(x) = -a_x_0*(x(-1)-y(-1)) + a_x_1*diff(x(-1)) + a_x_2*diff(x(-2)) + ex ;
+
+[name='eq:z']
+diff(z) = -b_z_0*(z(-1)-x(-1)) + b_z_1*diff(z(-1)) + b_z_2*diff(z(-2)) + ez ;
+
+end;
+
+get_ar_matrices('toto');
\ No newline at end of file