check if initial states are empty, and set them to zero.

matlab/+occbin/mkdata.m

@@ -30,7 +30,7 @@ function [data_mat]=mkdata(n_periods,dr_A,dr_B,endo_names,exo_names,wish_list,sh
 % given decision rule
 neqs = size(endo_names,1);
 
-if nargin<9
+if nargin<9 || isempty(var_init)
     var_init = zeros(neqs,1);
 end
 

matlab/+occbin/mkdatap_anticipated_dyn.m

@@ -48,7 +48,7 @@ T = DM.decrulea;
 CONST = zeros(n_vars,1);
 R = DM.decruleb;
 
-if nargin<7
+if nargin<7 || isempty(init)
     init=zeros(n_vars,1);
 end