steady state BGG financial accelerator

Hello everyone!
I am trying to replicate the steady state risk premium and capital/net worth ratio in BGG(1999). In particular, I am writing a Matlab m-file to compute the steady state.
I would like to program equations (A.1) and (A.2) in the appendix to obtain the BGG steady state values for the risk premium (200bp), the business failure rate (F=0.03 annualized) and the capital/net worth ratio (K/N = 2).
I programmed equations (A.1) and (A.2) as explained in the text, and used the BGG calibration: sigma = 0.28, mu = 0.12, gamma = 0.9728. The problem is that I do not get the steady state risk premium (s) obtained by BGG.
My code is in attachment: can anyone spot the mistake in it?
Thank you very much for your help!
mary
I am trying to replicate the steady state risk premium and capital/net worth ratio in BGG(1999). In particular, I am writing a Matlab m-file to compute the steady state.
I would like to program equations (A.1) and (A.2) in the appendix to obtain the BGG steady state values for the risk premium (200bp), the business failure rate (F=0.03 annualized) and the capital/net worth ratio (K/N = 2).
I programmed equations (A.1) and (A.2) as explained in the text, and used the BGG calibration: sigma = 0.28, mu = 0.12, gamma = 0.9728. The problem is that I do not get the steady state risk premium (s) obtained by BGG.
My code is in attachment: can anyone spot the mistake in it?
Thank you very much for your help!
mary