s<-numeric()
for(i in 1:100)
{
x<-runif(4,min=0,max=10)
dim(x)<-c(1,4)
Y<-x[,1]
Z<-x[,2]
R<-x[,3]
D<-x[,4]
y<-min(20,Y)
z<-min(20-y,Z)
r<-min(20-y-z,R)
s[i]<-20-y-z-r
}
mean(s)

感谢！我以后会好好学R的。。