I have the following system of equations, which is my model block:
C^(-sigma) = N^phi/W;
K(+1) = (1-delta)K+I;
Y = AN^(1-alpha)K^alpha;
W = (1-alpha)AN^(-alpha)K^alpha;
R = alphaAN^(1-alpha)K^(alpha-1);
Y = C + I;
B = 0;
A = A(-1)^rhoexp(epsilon);
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!