Thank you jpfeifer. I don't know what I have missed in my steadystate file, but I just copied the steady states of all my variables into my .mod and it worked.
I am having another problem now, I am trying to get the irf but I get the following error message:
Error using chol
Matrix must be square.
Error in irf (line 80)
VC_epschol = chol(VC_eps)';
I found on the forum that if I use two different but correlated shocks, Dynare performs a Cholesky decomposition. In my .mod file, I have several shocks but they are not correlated.
I looked into the error message, it seems that the matrix that is the source of problems is the covariance matrix of reduced-form residuals. Is my problem in my shocks declaration or could it be something else?