Request for help: Warning: the matrix is close to singular value, or scaling error. The results may not be accurate

Hello! I’m new. I can’t find the steady state of my model. Could you help me?
model.mod (2.2 KB)

Using 64-bit preprocessor
Starting Dynare (version 4.6.1).
Calling Dynare with arguments: none
Starting preprocessing of the model file …
Found 17 equation(s).
Evaluating expressions…done
Computing static model derivatives (order 1).
Computing dynamic model derivatives (order 1).
Processing outputs …
done
Preprocessing completed.
Warning: the matrix is close to singular value, or scaling error. The results may not be accurate.RCOND = 2.831100e-21.
In trust_region>dogleg at 198
In trust_region at 115
In dynare_solve at 253
In evaluate_steady_state at 221
In steady_ at 55
In steady at 80
In driver at 365
In dynare at 293
In trust_region>dogleg at 198
In trust_region at 115
In dynare_solve at 253
In evaluate_steady_state at 221
In steady_ at 55
In steady at 80
In driver at 365
In dynare at 293
Residuals of the static equations:

Equation number 1 : 0.0058541 : 1
Equation number 2 : 0.014134 : k
Equation number 3 : 0.0010479 : 3
Equation number 4 : 0.025977 : 4
Equation number 5 : -2.4743e-05 : mc
Equation number 6 : 0.0017528 : 6
Equation number 7 : -0.0015593 : la
Equation number 8 : 0.0029416 : p_star
Equation number 9 : -0.065476 : x1
Equation number 10 : -0.016099 : x2
Equation number 11 : 0 : pi
Equation number 12 : 0 : 12
Equation number 13 : 0 : s
Equation number 14 : 0.0018893 : y
Equation number 15 : -0.014115 : g
Equation number 16 : -0.014463 : 16
Equation number 17 : 0 : 17

Error print_info (line 32)
Impossible to find the steady state (the sum of square residuals of the static equations is 0.0059). Either the model doesn’t have a steady
state, there are an infinity of steady states, or the guess values are too far from the solution

Maybe this video by Willi Mutschler might help.

Thanks for your help, I’ll learn from this video.

It needs to be

p_star=(theta_j/(theta_j-1))*(x1/x2);

instead of

p_star=(theta_j/(1-theta_j))*(x1/x2);
1 Like

I have modified my code and find the steady state. Thanks a lot!