Robert_Hoo ifelse (u < prob/(1+prob), z[j]<-1, z[j]<-1)应该出错了；后面两个选项一样的 而且ifelse也不是这么赋值的，应该是 z[j] <- ifelse(u < prob/(1+prob), 1, 0) 再次看一个程序可不可以用parallel运算，首先要确定第n个运算不依赖于n-1的运算结果......也就是说每次运算都是独立的
