I have a model that I am trying to love. The following shows the error for each of the equation being close to zero.

Equation number 1 : -0.00019977 : Lambda

Equation number 2 : 0 : M

Equation number 3 : 0 : p

Equation number 4 : 0 : 4

Equation number 5 : -0.0010725 : 5

Equation number 6 : -0.0027062 : 6

Equation number 7 : 0 : 7

Equation number 8 : -8.6782e-05 : z_r

Equation number 9 : -7.9136e-05 : z_lr

Equation number 10 : -6.6192e-06 : y

Equation number 11 : -6.6173e-06 : pi

Equation number 12 : -6.9881e-06 : d

Equation number 13 : 0 : D

Equation number 14 : 0 : 14

Equation number 15 : 0 : 15

Equation number 16 : 0.001224 : 16

Equation number 17 : 0.0013341 : 17

Equation number 18 : -0.003505 : 18

Equation number 19 : -0.00094701 : 19

Equation number 20 : 0 : lev

Equation number 21 : 0 : log_p

Equation number 22 : 0 : log_e

Equation number 23 : 0 : log_md

Equation number 24 : 0 : log_mc

Equation number 25 : 0 : log_r

Equation number 26 : 0 : log_lev

Equation number 27 : 0 : log_ar

Equation number 28 : 0 : log_alr

Equation number 29 : 0 : log_as

However, there is still the “impossible to find steady state” error and it shows that “sum of square residuals of the static equations” is 0.0000. The initial guess is solved directly using an external .m file.

Why is this happening? Can you give me some help to make this go through?

I followed suggestion from another post by adding

options_.dynatol.f=5e-4;

before the steady command, but I this does not seem to do anything.

Can you please give me some help? Thank you!