The global oo_.steady_state should also be updated in the homotopy loop

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

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

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