Hi everyone,

I have the following system of equations, which is my model block:

C^(-sigma)=(beta/Q)*C(+1)^(-sigma);
C^(-sigma) = N^phi/W;
C^(-sigma) =beta*(C(+1))^(-sigma)*(1-delta+R(+1));

K(+1) = (1-delta)

*K+I;*

Y = AN^(1-alpha)

Y = A

*K^alpha;*

W = (1-alpha)N^(1-alpha)

W = (1-alpha)

*A*N^(-alpha)*K^alpha;*

R = alphaAR = alpha

*K^(alpha-1);*

Y = C + I;

B = 0;

A = A(-1)^rhoexp(epsilon);

Y = C + I;

B = 0;

A = A(-1)^rho

I need to calculate the steady state by hand, such that the steady state values depend only on parameters like phi, epsilon, alpha etc.

As this is kind of tedious i was wondering if there is a way to let dynare/matlab solve this system, to compare with my own results?

Thanks in advance!

Regards,

kartoffel