<br />
Num<-2500<br />
M<-2.67<br />
u<-runif(Num,max=M)<br />
y<-runif(Num)<br />
plot(y,u)<br />
col.y <- ifelse(u > dbeta(y,2.7,6.3),'red','blue')<br />
points(y,u,col = col.y,pch = 20)<br />
c=seq(0,1,by=0.01)<br />
lines(c,dbeta(c,2.7,6.3))<br />
</p>