set.seed(1)
p <- numeric(5200)
x <- (ifelse (rbinom (5200,1,0.5) == 1,1.005,0.995))  
for (i in 1:5200) {
    if (i==1) {
    p[i] <- 10*x[i]
}
    else {
    p[i] = x[i]*p[i-1]
}
}
plot(p, type="l")


s <- numeric()
n <- numeric()
share.v <- numeric()
j <- exp(0.07)/(260*20)-1
for (i in 1:40) {
   
    if (i==1) {
    s[i] <- 0
    n[i] <- 100
    share.v[i] <- p[i] * n[i-1]
}
   
    else {
         if (p[i] >= 10.1 && p[i] < 10.1 && share.v > 0) {
         s[i]=p[1]*n[i-1]*(1-0.001)
}
         else {
              if (p[i] <= 9.9 && p[i-1] < 9.9 && s[i-1] > 0） {
              n[i]=s[i-1]/p[i]
}
               else {
                    p[i] = p[i-1] * x[i]
                    s[i] = s[i-1] * (1+j)
}
}
}
}


到底哪错了呀？跪求正解！！！