clear all; alpha=0.35; delta=0.06; sigma=0.20; kstar = (sigma/delta)^(1/(1-alpha)); k(1) = 0.5*kstar; for t = 1:50 y(t) = k(t)^alpha; i(t) = sigma*y(t); k(t+1) = (1-delta)*k(t) + i(t); end disp(k(51)/kstar)