# 2005 年 10 月 25 日作成
rep.no<-10000
no.sample<-4
x<-rep(0,rep.no)
for (i in 1:rep.no){
r<-runif(no.sample,0,1)
x[i]<-max(r)
}
plot(density(no.sample*(1-x)),xlim=c(-1,10),ylim=c(0,0.8),col="red")
par(new=T)
no.sample<-16
x<-rep(0,rep.no)
for (i in 1:rep.no){
r<-runif(no.sample,0,1)
x[i]<-max(r)
}
plot(density(no.sample*(1-x)),xlim=c(-1,10),ylim=c(0,0.8),col="blue",xlab="",ylab="")
par(new=T)
no.sample<-64
x<-rep(0,rep.no)
for (i in 1:rep.no){
r<-runif(no.sample,0,1)
x[i]<-max(r)
}
plot(density(no.sample*(1-x)),xlim=c(-1,10),ylim=c(0,0.8),col="green",xlab="",ylab="")
par(new=T)
z<-seq(0,10,0.01)
plot(z,dexp(z,1),xlim=c(-1,10),ylim=c(0,0.8),type="l",xlab="",ylab="")