The basic code should be something along the lines of
y_pos=strmatch('Y',M_.endo_names,'exact');
alpha_parameters = 0.1:0.05:1;
delta_parameters = 0.1:0.3:0.8;
[alpha_mat, delta_mat]=meshgrid(alpha_parameters,delta_parameters);
temp_mat=NaN(size(alpha_mat));
for alpha_iter=1:length(alpha_parameters)
for delta_iter=1:length(delta_parameters)
set_param_value('alpha',alpha_mat(delta_iter,alpha_iter));
set_param_value('delta',delta_mat(delta_iter,alpha_iter));
[oo_.steady_state,M_.params,info] = steady_(M_,options_,oo_);
if ~info
temp_mat(delta_iter,alpha_iter)=oo_.steady_state(y_pos);
end
end
end;
contour(alpha_mat,delta_mat,temp_mat)
But you most probably need a steady_state_model
-block as the steady state is not found for many parameter combinations.