The steady state cannot be computed

dear all:
Can anyone help me? I have stuck on running my model with dynare. here are the mistakes:

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

DYNARE_SOLVE (solve_algo=2|4): number of blocks = 17
DYNARE_SOLVE (solve_algo=2|4): solving block 16, of size 1
DYNARE_SOLVE (solve_algo=2|4): solving block 15, of size 1
DYNARE_SOLVE (solve_algo=2|4): solving block 14, of size 1
DYNARE_SOLVE (solve_algo=2|4): solving block 13, of size 1
DYNARE_SOLVE (solve_algo=2|4): solving block 12, of size 1
DYNARE_SOLVE (solve_algo=2|4): solving block 11, of size 1
DYNARE_SOLVE (solve_algo=2|4): solving block 10, of size 24
DYNARE_SOLVE (solve_algo=2|4): number of blocks = 17
DYNARE_SOLVE (solve_algo=2|4): solving block 16, of size 1
DYNARE_SOLVE (solve_algo=2|4): solving block 15, of size 1
DYNARE_SOLVE (solve_algo=2|4): solving block 14, of size 1
DYNARE_SOLVE (solve_algo=2|4): solving block 13, of size 1
DYNARE_SOLVE (solve_algo=2|4): solving block 12, of size 1
DYNARE_SOLVE (solve_algo=2|4): solving block 11, of size 1
DYNARE_SOLVE (solve_algo=2|4): solving block 10, of size 24
MODEL_DIAGNOSTICS: The steady state cannot be computed

Residuals of the static equations:

Equation number 1 : -0.065677 : 1
Equation number 2 : -0.17273 : lamda
Equation number 3 : 0 : 3
Equation number 4 : -0.082236 : k
Equation number 5 : 0.0091524 : 5
Equation number 6 : -0.0077002 : piwstar
Equation number 7 : -7.3052e-05 : f
Equation number 8 : 0.0018893 : 8
Equation number 9 : -0.15167 : 9
Equation number 10 : 0.058681 : g1
Equation number 11 : 0.0022349 : g2
Equation number 12 : 0.00020686 : 12
Equation number 13 : -0.0029377 : mc0
Equation number 14 : -0.013605 : mc
Equation number 15 : -0.094847 : 15
Equation number 16 : 0.0013818 : w
Equation number 17 : 0.13568 : r
Equation number 18 : 0.040313 : 18
Equation number 19 : 0.011237 : em
Equation number 20 : 0.058198 : emcost
Equation number 21 : 0.082243 : ce
Equation number 22 : 0.00081204 : y1
Equation number 23 : 0.0018955 : y2
Equation number 24 : 0.67179 : omega
Equation number 25 : -0.0036403 : Xat
Equation number 26 : -0.00070777 : Xup
Equation number 27 : -4.1356e-05 : Xlo
Equation number 28 : 0 : F
Equation number 29 : 3.2496 : Tat
Equation number 30 : 0 : Tlo
Equation number 31 : 0 : en
Equation number 32 : 0.063642 : y
Equation number 33 : 0.002838 : 33
Equation number 34 : 0.72026 : l
Equation number 35 : 0.0073632 : vp
Equation number 36 : 0.0020464 : vw
Equation number 37 : 0 : 37
Equation number 38 : 0 : 38
Equation number 39 : 0 : 39
Equation number 40 : 0 : 40
Equation number 41 : 0 : 41

错误使用 print_info (第 32 行)
Impossible to find the steady state (the sum of square residuals of the static equations is
11.6413). 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

If I may ask, what type of model are you running…?

nk-dsge model

Oh ok. If you wrote the model from scratch, maybe you can start from some baseline NK model that works and extend it little by little…

yes,dear!I did follow your advice to do it!but i’m really confused about the codes,could you please solve the problems? many thanks!

Oh ok. If you add features little by little, I guess you could know which stage the problem occurred, and then try to find out what caused it.

For a model that size, I would recommend either starting with a simpler model that works or computing the steady state analytically. In any case, without the file it is impossible to tell what is going on.

sorry,professor.I put the file there.
oceansink.mod (8.2 KB)
myfunction.m (120 Bytes)

I am getting

ERROR: There are 41 equations but 40 endogenous variables!

oh!sorry professor! I forget it.thanks so much!