Hi, guys:

Can I use fsolve inside a mod file, since the computation of the steady state involves some nonlinear roots finding process.

To be exact, can I just put the following line in the “model” block?

#para1=1.0039;

#para2=0.0001;

#ha_ja=.038;

#xx=fsolve(@(x) para2*x^(sigma-1)-1/ha_ja*x^sigma-para1,.4);

#ha_oz=xx(1);

If I cannot do that, how can I solve the steady state if it cannot be done with paper and pen?