diff --git a/matlab/homotopy1.m b/matlab/homotopy1.m
index 8b8305d827b61ad793498b537b6eac1540ee7c36..f0913b86bb3bd0ff9f219c257f245f9b2d29bf23 100644
--- a/matlab/homotopy1.m
+++ b/matlab/homotopy1.m
@@ -77,5 +77,5 @@ for i=1:step_nbr+1
     oo_.exo_steady_state(values(ix,2)) = points(ix,i);
     oo_.exo_det_steady_state(values(ixd,2)) = points(ixd,i);
 
-    steady_(M_,options_,oo_);
+    oo_.steady_state = steady_(M_,options_,oo_);
 end
diff --git a/matlab/homotopy2.m b/matlab/homotopy2.m
index 3c972a0ee4d9c523510528a5f6cc322cf839342f..726a8aac09e6491d8309fe91f450316ba143b97b 100644
--- a/matlab/homotopy2.m
+++ b/matlab/homotopy2.m
@@ -100,6 +100,6 @@ for i = 1:nv
 
         disp([ 'HOMOTOPY mode 2: lauching solver with ' deblank(varname) ' = ' num2str(v) ' ...'])
         
-        steady_(M_,options_,oo_);
+        oo_.steady_state = steady_(M_,options_,oo_);
     end
 end
diff --git a/matlab/homotopy3.m b/matlab/homotopy3.m
index 9075a4a1062e280ecc8ca5c264828adacc99febd..5359680bfbb92c07f9e68374f3739bfce68fddf6 100644
--- a/matlab/homotopy3.m
+++ b/matlab/homotopy3.m
@@ -91,7 +91,7 @@ while iter < step_nbr
     old_ss = oo_.steady_state;
 
     try
-        steady_(M_,options_,oo_);
+        oo_.steady_state = steady_(M_,options_,oo_);
         
         if length([kplus; kminus]) == nv
             return