@Albert_Zevelev: I appreciate very much your help; extremely useful.
Just three comments:
- This example, I made it up, so I am sorry that I can
t provide the source because there isnt one. - The steady state solution (Y_rec = Ym, in your notation) is:
Y = (1/1-alpha) * (a + G). This means that the value of beta should not affect the steady state solution. I thus wonder why, with the code, changing the value of beta changes the results. - In the simulation on the first graph, Y, C, and I seem to move down to zero, while the closed solution indicates that, in steady state, Y=600, C=500 and I should be zero from the investment equation.
Thanks again for your help, and I apologize for asking so many questions.