Unable to find out steady state

[mhsjnu]

Residuals of the static equations:

Equation number 1 : -10.4446
Equation number 2 : 3.8
Equation number 3 : 1.6955
Equation number 4 : Inf
Equation number 5 : 0.86647
Equation number 6 : -0.31565
Equation number 7 : 0
Equation number 8 : Inf
Equation number 9 : -0.95
Equation number 10 : 2.8183



STEADY: The Jacobian contains Inf or NaN. The problem arises from:

STEADY: Derivative of Equation 4 with respect to Variable y (initial value of y: 5000)
STEADY: Derivative of Equation 8 with respect to Variable y (initial value of y: 5000)

STEADY: The problem most often occurs, because a variable with
STEADY: exponent smaller than 1 has been initialized to 0. Taking the derivative
STEADY: and evaluating it at the steady state then results in a division by 0.

[jpfeifer]

You cannot use 5000 as a starting value for log output as you will be using

Code: Select all

exp(5000)=Inf

in the model. Also

Code: Select all

exp(g)=rhog*exp(g(+1)/g)+eg;

does not look like a correct equation. Why is there a division by g on the RHS?

[mhsjnu]

Dear Professor Jpfeifer, Thanks for valuable suggestion and guide line. But I am unable to solve after correction. I try to solve for steady state value by using MATLAB under .m file. its looks very difficult. Is there any helpful reference for getting initial value? Please help for solving this model.

[jpfeifer]

Your equations are still wrong.

Code: Select all

exp(z)=rhoz*exp(z(-1))+e;

has steady state

Code: Select all

eps(z)=0

implying that

Code: Select all

z=-Inf

Please make sure the exp() is consistently used in a correct way.

[mhsjnu]

Please help, Regarding this issues. Thanks in advanced.
"An infinite element was encountered when trying to solve equation(s) 3
with respect to the variable(s): pi."

[jpfeifer]

Again, make sure the exp()-substitution is consistent. In the first equation,

Code: Select all

rn

is not in exp(). In

Code: Select all

exp(Ln((1+rn(+1))/(1+rn)))

it is. For the start, I would make the model run without this type of substitution and only then add the exp()

[mhsjnu]

Now I try to find out Steady state value by using "steady_state_model" command for the instant of the initial value. Plz help and thanks in Advanced.

[jpfeifer]

I can only repeat myself. You did not head any of my previous advice. Note that a

Code: Select all

steady_state_model

block requires the full analytical steady state to be computed. It does not seem you did this. For example, z has steady state 0 in the entered AR-process. You have a value for c hard-coded, which also does not seem plausible given that you used an exp()-substitution that implies c is in logs.

[mhsjnu]

Dear Professor Johannes Pfeifer, I follow your suggestions from the first. But I can not find out the solution. So I try to another way. Now I can understand that it's very difficult. Plz see the attachment, I have corrected by your last suggestions.

[jpfeifer]

See my very first comment on

Code: Select all

(ln((1+exp(rn))/(1+rn)))=(((1-thetar)*((thetapi*(ln(exp(pi)/pi))+thetay*(ln(exp(y)/y))+thetau*(ln(exp(u))/u))))+thetar*ln((1+exp(rn(-1)))/(1+rn)));

The way you do the exp()-substitution is still not consistent

[mhsjnu]

I have done. Is it ok?

[jpfeifer]

This is leading nowhere. For the last time, please sit down and check every single equation on whether the exp()-substitution has been done consistently. The last equation for example has an

Code: Select all

n(-1)

not in exp() while it is in exp() everywhere else. As I suggested in the beginning, you should try to solve the model first without any substitution and make sure it runs. That makes debugging much easier.

[mhsjnu]

I rewrite this equation as follows

(1+exp(rn(+1))/(1+rn))=((1-thetar)*((thetapi*(exp(pi(+1))/pi))+thetay*(exp(y(+1))/y)+thetau*(exp(u(+1))/u)))+thetar*(1+exp(rn(-1))/(1+rn));

Please see the attachment. techT.mod

And Alternative file without substitution: Attachment: techW.mod