<p>I am new to Dynare++ and I have really quick question I cannot seem to find the answer too.</p>

<p>What is the difference between these two commands and why is the output different?</p>

<pre><code>!dynare++ --per 20 --sim 3 file_name.mod

dynare file_name.mod
</code></pre>

<p>In the first command its unable to find steady state values based on my initial values and in the second it can. Why?</p>

<p>Here is the mod file of the basic RBC model I want to simulate with the first command. I can do so with the second, but not the first. </p>

<pre><code>var c, l, k, w, r, y, i, z;
varexo eps;

parameters beta, delta, alpha, a, rho, rbar, wbar;

beta=.99;
delta=.025;
alpha=.33;
a=1;
rho=.9;
rbar=1/beta+delta-1; 
wbar=(1-alpha)*(alpha/rbar)^(alpha/(1-alpha)); 


model;
//Inter-Euler
1/c=beta*(1/c(+1))*(1+r(+1)-delta);
//Intra-Euler
a/(1-l)=w/c;
//H.H. Budget Constraint
k=w(-1)*l(-1)+(1-r(-1)-delta)*k(-1)-c(-1);
//Firm FOC Capital
r=exp(z)*alpha*exp(z)*(l/k)^(1-alpha);
//Firm FOC Labor
w=exp(z)*(1-alpha)*exp(z)*(k/l)^(alpha);
// Production Function
y=exp(z)*(k^alpha)*(l^(1-alpha));
i=y-c;
z=rho*z(-1)+eps;
end;

initval;
k=((wbar/a)/(wbar*(1+1/a)-(r+delta)*(alpha/rbar)^(1/(1-alpha))))*(alpha/rbar)^(1/(1-alpha));  
c=wbar*(1-((wbar/a)/(wbar*(1+1/a)-(r+delta)*(alpha/rbar)^(1/(1-alpha)))))/a;  
l=(wbar/a)/(wbar*(1+1/a)-(r+delta)*(alpha/rbar)^(1/(1-alpha))); 
w=wbar; 
r=rbar; 
y=(k^alpha)*(l^(1-alpha));
i=(k^alpha)*(l^(1-alpha))-(wbar*(1-((wbar/a)/(wbar*(1+1/a)-(r+delta)*(alpha/rbar)^(1/(1-alpha)))))/a);
z=0;
end;

shocks;
var eps;
stderr .02;
end;
</code></pre>

CC BY-SA 3.0
2017-09-24T03:43:04.190Z
Hi. Would you mind sharing the code in the mod file? It might help to determine what the issue is. You might consider producing a MWE ("minimum working example" that is as short as possible while still producing the problem in question).
jmbejara
CC BY-SA 3.0
2017-09-24T04:50:03.557Z
I just updated my post with the mod file
jessica
<dynare>
Running Dynare from Matlab
jessica
# Running Dynare from Matlab

(Asked by: username_0 on 2017-09-24)

I am new to Dynare++ and I have really quick question I cannot seem to find the answer too.

What is the difference between these two commands and why is the output different?

```
!dynare++ --per 20 --sim 3 file_name.mod

dynare file_name.mod

```

In the first command its unable to find steady state values based on my initial values and in the second it can. Why?

Here is the mod file of the basic RBC model I want to simulate with the first command. I can do so with the second, but not the first.

```
var c, l, k, w, r, y, i, z;
varexo eps;

parameters beta, delta, alpha, a, rho, rbar, wbar;

beta=.99;
delta=.025;
alpha=.33;
a=1;
rho=.9;
rbar=1/beta+delta-1; 
wbar=(1-alpha)*(alpha/rbar)^(alpha/(1-alpha)); 

model;
//Inter-Euler
1/c=beta*(1/c(+1))*(1+r(+1)-delta);
//Intra-Euler
a/(1-l)=w/c;
//H.H. Budget Constraint
k=w(-1)*l(-1)+(1-r(-1)-delta)*k(-1)-c(-1);
//Firm FOC Capital
r=exp(z)*alpha*exp(z)*(l/k)^(1-alpha);
//Firm FOC Labor
w=exp(z)*(1-alpha)*exp(z)*(k/l)^(alpha);
// Production Function
y=exp(z)*(k^alpha)*(l^(1-alpha));
i=y-c;
z=rho*z(-1)+eps;
end;

initval;
k=((wbar/a)/(wbar*(1+1/a)-(r+delta)*(alpha/rbar)^(1/(1-alpha))))*(alpha/rbar)^(1/(1-alpha));  
c=wbar*(1-((wbar/a)/(wbar*(1+1/a)-(r+delta)*(alpha/rbar)^(1/(1-alpha)))))/a;  
l=(wbar/a)/(wbar*(1+1/a)-(r+delta)*(alpha/rbar)^(1/(1-alpha))); 
w=wbar; 
r=rbar; 
y=(k^alpha)*(l^(1-alpha));
i=(k^alpha)*(l^(1-alpha))-(wbar*(1-((wbar/a)/(wbar*(1+1/a)-(r+delta)*(alpha/rbar)^(1/(1-alpha)))))/a);
z=0;
end;

shocks;
var eps;
stderr .02;
end;

```

**Tags:** dynare

> **Comments:**
> **username_1**: Hi. Would you mind sharing the code in the mod file? It might help to determine what the issue is. You might consider producing a MWE ("minimum working example" that is as short as possible while still producing the problem in question).  
> **username_0**: I just updated my post with the mod file  

---

## 0 answers

