R<-cumprod(1+r[(timeR-2):(timeR+30)])

sum((C_yun[(timeR-2):(timeR+30)]-C_yan[(timeR-2):(timeR+30)])/R)+sum((C_oun[(timeR-2):(timeR+30)]-C_oan[(timeR-2):(timeR+30)])/R)