The simplified RBC model below seems to satisfy BK when worked out by hand (at great personal expense :)) but Dynare balks. Since I am new to both Dynare and RBC, maybe I am missing something obvious, so any help would be appreciated:

```
var k, c, epsilon;
varexo u;
parameters alpha, rho, delta, rhoeps, sigma;
alpha=1/3.0;
delta = 0.025;
rho = 0.02;
rhoeps = 0.95;
sigma = 0.05;
model;
log(c) = (1-delta + epsilon(+1)*k(+1)^alpha)*log(c(+1))/(1+rho);
c+k(+1)-(1-delta)*k = epsilon*k^alpha;
log(epsilon(+1)) = rhoeps * log(epsilon) +u(+1);
end;
initval;
k = (rho+delta)^(1/alpha);
c = rho + delta*(1-(rho+delta)^(1/alpha));
epsilon=1;
u=0;
end;
steady_state_model;
epsilon=1;
k = (rho+delta)^(1/alpha);
c = rho + delta - delta*(rho+delta)^(1/alpha);
end;
shocks;
var u; stderr 0.05;
end;
```