# Drawing an isocurve on Dynare

**URL:** <https://forum.dynare.org/t/drawing-an-isocurve-on-dynare/19612>\
**Category:** Stochastic simulations\
**Created:** [10 January 2022 13:15 UTC](https://forum.dynare.org/t/drawing-an-isocurve-on-dynare/19612 "2022-01-10T13:15:10Z")\
**Posts on this page:** 1\
**Showing post:** 6

<div class="post-metadata">

**Author:** ![jpfeifer](https://forum.dynare.org/user_avatar/forum.dynare.org/jpfeifer/32/5044_2.png) [@jpfeifer](https://forum.dynare.org/u/jpfeifer)\
**Post date:** [11 January 2022 16:18 UTC](https://forum.dynare.org/t/drawing-an-isocurve-on-dynare/19612/6 "2022-01-11T16:18:36Z")

</div>

The basic code should be something along the lines of

```auto
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.

---

_[View the full topic](https://forum.dynare.org/t/drawing-an-isocurve-on-dynare/19612)._