#First, define variables

T<-200
timeR<-101

eps<-10^(-10)


alpha<-0.4
rho<-(1+0.01)^35-1
tau<-rep(0.2,T)

n<-rep(0,T)
L<-cumprod((1+n))

k<-rep(1,T)
K<-rep(1,T)
L<-rep(1,T)

Y<-rep(1,T)
w<-rep(0,T)
r<-rep(0,T)

s<-rep(0.5,T)
C_y<-rep(0.5,T)
C_o<-rep(0.5,T)

U<-rep(0.5,T)