model;
(1/c) = beta*(1/c(+1))(1 - delta + theta(k^(theta-1)(ki^phi)(h(+1)^(1-theta-phi))));
(1/c)q = beta(1/c(+1))(1 - deltai + phi(k^theta)(ki^(phi-1))(h(+1)^(1-theta-phi)));
psic/(1-h) = ((1-theta-phi)(k(-1)^theta)(ki(-1)^phi))/(h^(theta+phi));
c = y-i;
i = k+kiq+(1-delta+theta*(k^(theta-1)(ki^phi)(h^(1-theta-phi)))(k(-1))+(1-deltai+(phi(k^theta)(ki^(phi-1))(h^(1-theta-phi))))(ki(-1)));
y = (k(-1)^theta)(ki(-1)^phi)*(h^(1-theta-phi));

end;

initval;
ki = ((deltai + (1/beta) - 1)/((phi*((delta + (1/beta)- 1)/theta)^(theta/(theta-1))(h^(1-theta-phi)))))^(1/(theta+phi-1));
k = ((delta + (1/beta)- 1)/(theta(ki^phi)(h^(1-theta-phi))))^(1/(theta-1));
i = (delta(((delta + (1/beta)- 1)/(theta*(ki^phi)(h^(1-theta-phi))))^(1/(theta-1))))+(deltai((delta + (1/beta)- 1)/(theta*(ki^phi)(h^(1-theta-phi))))^(1/(theta-1)));
y = (k^theta)(ki^phi)(h^(1-theta-phi));
c = y - i;
psi = (c(1-theta-phi)(k^theta)(ki^phi))/((h^(theta+phi))*(1-h));
end;

resid;
steady;[/code]
As you can see, your computed values do not match the model equations. So either your steady state computation or the equations are still wrong. focus on the equations with a large residual

I checked my equations and there were a few errors. Now, I have managed to get all the residuals down to zero except for the second equation. I checked my model equation and there seems to be nothing wrong with it. So I changed one of my parameter values (deltai) for convergence to steady state.

When I run the command steady, there seems to be an error: “Not enough input arguments”.

Starting Dynare (version 4.4.3).
Starting preprocessing of the model file …
Found 6 equation(s).
Evaluating expressions…done
Computing static model derivatives:

order 1
Computing dynamic model derivatives:

order 1

order 2
Processing outputs …done
Preprocessing completed.
Starting MATLAB/Octave computing.

Residuals of the static equations:

Equation number 1 : 0
Equation number 2 : 3.0713e-05
Equation number 3 : 0
Equation number 4 : 0
Equation number 5 : 0
Equation number 6 : 0

EIGENVALUES:
Modulus Real Imaginary

3.33e-16 3.33e-16 0
0.9506 0.9506 0
1.075 1.075 0

There are 1 eigenvalue(s) larger than 1 in modulus
for 1 forward-looking variable(s)

The rank condition is verified.

Error using steady (line 4)
Not enough input arguments.