Dear Michel,

when moving from 4.2 to 4.3

we changed line 122 of solve1 from

p=-(fjac2+sqrt(nn*eps)*max(sum(abs(fjac2)))*eye(nn))\(fjac'*fvec);

p=-(fjac2+1e6*sqrt(nn*eps)*max(sum(abs(fjac2)))*eye(nn))\(fjac'*fvec);

That is, we added a constant of 1e6. Is there a good reason for this? I tried to use the replication code of Eggertsson, Ferrero und Raffo (2014), written for 4.2.5 and it is unable to find the steady state in 4.3 because this constant was added. As this is not the first time people report that 4.2 succeeds in finding steady states that 4.3 did not find, I would suggest moving back.

Best,

Johannes

Johannes Pfeifer

Friesenwall 104

50672 Köln

Tel.: +49-221-29873852

Mobil.: +49-170-6936820

Germany