#source("c:/Work07/ShortCourse07/Programs/BetaAR.txt",print.eval=TRUE)#
#This gives an accept-reject algorithm for a beta rv#
a<-2.7; b<-6.3; c<-2.669;nsim<-2500;
#Generate u and v#
u<-runif(nsim);v<-runif(nsim);
#Generate Y, the beta random variable#
test<-dbeta(v, a, b, ncp=0, log = FALSE)/c; ind<-(u<test);
plot(function(x)(dbeta(x, a, b, ncp=0, log =FALSE)/c),xlim=c(0,1),ylim=c(0,1),ylab="Density");
par(new=T);
plot(v,u,xlim=c(0,1),ylim=c(0,1),xlab="",ylab="");par(new=T);
plot(v*(u<test),u*(u<test),xlim=c(0,1),,ylim=c(0,1),xlab="",ylab="",col="red",pch=19)