Attached is a very simple mod file in Dynare. and we can calculate the steady state value of capital as

K=((1/beta-(1-delta))/alpha)^(1/(alpha-1));

Then compare this solution with the one dynare gives out, please notice that dynare gives out log(K) since I use a log-linearized version in dynare

It looks the solution from dynare is as big as 4 times of the right solution, or I made a mistake in the mod file.

dyn_jh_ch_ex1.mod (662 Bytes)