When you want to get the steady state of ld in the NK_baseline_steadystate file, after the function fzero(@(ld)(1-betta*thetaw…, you set the initial value to 0.25, why?
How was this number obtained?
Just a guess. Usually people work around 1/5 to 1/3 of their time