diff --git a/matlab/stochastic_solvers.m b/matlab/stochastic_solvers.m
index 9967319b31c328c49f1bd997126f80e07c918cdc..edecb1f860caad56dbeb83bbe1bbc24121559a35 100644
--- a/matlab/stochastic_solvers.m
+++ b/matlab/stochastic_solvers.m
@@ -48,6 +48,10 @@ function [dr,info] = stochastic_solvers(dr,task,M_,options_,oo_)
 
 info = 0;
 
+if options_.linear
+    options_.order = 1;
+end
+
 if (options_.aim_solver == 1) && (options_.order > 1)
         error('Option "aim_solver" is incompatible with order >= 2')
 end