See
In your steady state file, in
y_k = (aux1 - aux2 + aux3) * 1/ ((1-tau)) * 1/alpha_p;
the variable aux2 is bigger than aux1+aux3, which results in complex values.
See
In your steady state file, in
y_k = (aux1 - aux2 + aux3) * 1/ ((1-tau)) * 1/alpha_p;
the variable aux2 is bigger than aux1+aux3, which results in complex values.