Estimate the Laffer curve with a DSGE on Dynare

General DSGE Modeling
13

Dear Professor @jpfeifer ,
I need to do the previous looping for two parameters (varpi & psi).

tau_parameters = 0.1:0.3:0.8;
temp_mat=NaN(size(oo_.steady_state,1),length(tau_parameters));
for i=1:length(tau_parameters)
for j=1:length(tau_parameters)
set_param_value(‘varpi’,tau_parameters(i));
set_param_value(‘psi’,tau_parameters(j));
[oo_.steady_state,M_.params,info] = steady_(M_,options_,oo_);
if ~info
temp_mat(:,i,j)=oo_.steady_state;
end
end
end;

Then plot the isocurve of Y. On matlab, we usually do
“contour(varpi,psi,Y)”

But on Dynare, which I need, I don’t know what to do and thanks in advance.