Steady state exists but not for Dynare

Dear Jpfeifer,

Thanks in advance for your help.

We are trying to simulate a medium scale model. When we compute the steady state via fsolve in a simple .m file, the model converges immediately with no issue. When computing it via Dynare we get the typical message of steady state having NANs or Infs, despite we use as initial values the solution of the .m file. We thought reading your answer in another topic that our residual were too large but after setting
options_.dynatol.f=5e-3;
or even
options_.dynatol.f=5e-2;

nothing changes and we still get a message like:

Using 64-bit preprocessor
Starting Dynare (version 4.6.2).
Calling Dynare with arguments: none
Starting preprocessing of the model file …
Found 42 equation(s).
Evaluating expressions…done
Computing static model derivatives (order 1).
Computing dynamic model derivatives (order 1).
Processing outputs …
done
Preprocessing completed.

Residuals of the static equations:

Equation number 1 : 5.6712e-05 : 1
Equation number 2 : 1.0213e-05 : 2
Equation number 3 : -4.0923e-05 : 3
Equation number 4 : 7.9075e-05 : 4
Equation number 5 : -1.7714e-06 : 5
Equation number 6 : -0.00044968 : 6
Equation number 7 : 1e-05 : k_h
Equation number 8 : 2.5e-06 : k_c
Equation number 9 : -3.4076e-06 : u_lc
Equation number 10 : -3.4076e-06 : 10
Equation number 11 : 0.00011049 : 11
Equation number 12 : -0.00014858 : 12
Equation number 13 : -1.2992e-05 : 13
Equation number 14 : 0 : 14
Equation number 15 : 0 : 15
Equation number 16 : -3.4076e-06 : omega_h
Equation number 17 : -3.4076e-06 : omega_c
Equation number 18 : -0.0012 : rk_h
Equation number 19 : -0.0012 : rk_c
Equation number 20 : 0 : 20
Equation number 21 : 0 : d_b
Equation number 22 : 7.2e-07 : r
Equation number 23 : 3.6477e-05 : y_c
Equation number 24 : -8.2943e-05 : y_h
Equation number 25 : 5.0012e-05 : 25
Equation number 26 : -0.0001879 : 26
Equation number 27 : 0.00010764 : 27
Equation number 28 : -0.00028484 : 28
Equation number 29 : -7.8856e-06 : 29
Equation number 30 : -1.4658e-05 : 30
Equation number 31 : 0 : 31
Equation number 32 : -4.12e-05 : 32
Equation number 33 : 0 : c
Equation number 34 : 3.0624e-05 : u_bc
Equation number 35 : -3.834e-05 : u_bh
Equation number 36 : 0 : u_bcn
Equation number 37 : 0.0002 : u_bhn
Equation number 38 : 9.0914e-06 : 38
Equation number 39 : 1.2905e-05 : u_lh
Equation number 40 : 0 : u_lcn
Equation number 41 : 0 : u_lhn
Equation number 42 : 0 : chi

Randomize initial guess…

Residuals of the static equations:

Equation number 1 : NaN : 1
Equation number 2 : NaN : 2
Equation number 3 : NaN : 3
Equation number 4 : NaN : 4
Equation number 5 : NaN : 5
Equation number 6 : NaN : 6
Equation number 7 : NaN : k_h
Equation number 8 : NaN : k_c
Equation number 9 : NaN : u_lc
Equation number 10 : NaN : 10
Equation number 11 : NaN : 11
Equation number 12 : NaN : 12
Equation number 13 : NaN : 13
Equation number 14 : NaN : 14
Equation number 15 : NaN : 15
Equation number 16 : NaN : omega_h
Equation number 17 : NaN : omega_c
Equation number 18 : NaN : rk_h
Equation number 19 : NaN : rk_c
Equation number 20 : NaN : 20
Equation number 21 : NaN : d_b
Equation number 22 : NaN : r
Equation number 23 : NaN : y_c
Equation number 24 : NaN : y_h
Equation number 25 : NaN : 25
Equation number 26 : NaN : 26
Equation number 27 : NaN : 27
Equation number 28 : NaN : 28
Equation number 29 : NaN : 29
Equation number 30 : NaN : 30
Equation number 31 : NaN : 31
Equation number 32 : NaN : 32
Equation number 33 : NaN : c
Equation number 34 : NaN : u_bc
Equation number 35 : NaN : u_bh
Equation number 36 : NaN : u_bcn
Equation number 37 : NaN : u_bhn
Equation number 38 : NaN : 38
Equation number 39 : NaN : u_lh
Equation number 40 : NaN : u_lcn
Equation number 41 : NaN : u_lhn
Equation number 42 : NaN : chi

Error using print_info (line 32)
The steady state has NaNs or Inf.
Error in steady (line 102)
print_info(info,options_.noprint, options_);
Error in ModelCapital_v0_2_1.driver (line 551)
steady;
Error in dynare (line 293)
evalin(‘base’,[fname ‘.driver’]) ;

Do you have a suggestion?

Thanks again and best regards

Giorgio

Could you provide me with the files?

Hi Johannes,

In the meantime, reading the forum I discovered the model_diagnostics .

I solved the problem, a derivative was implying dividing by 0.

Thanks for the prompt response and for your amazing job in this forum.

Best

Giorgio