n<-100
y<-20

xVal <- seq(1000)/1000

par(mfrow=c(2,2))


a<-1
b<-1

plot(xVal,dbeta(xVal,a,b),xlim=c(0,1),ylim=c(0,10),main="",xlab="p",col=2,type="l",lwd=3)
par(new=TRUE)
plot(xVal,dbeta(xVal,a+y,b+n-y),xlim=c(0,1),ylim=c(0,10),main="1,1",xlab="p",type="l",lwd=3)

a<-10
b<-10

plot(xVal,dbeta(xVal,a,b),xlim=c(0,1),ylim=c(0,10),main="",xlab="p",col=2,type="l",lwd=3)
par(new=TRUE)
plot(xVal,dbeta(xVal,a+y,b+n-y),xlim=c(0,1),ylim=c(0,10),main="10,10",xlab="p",type="l",lwd=3)

a<-1
b<-5

plot(xVal,dbeta(xVal,a,b),xlim=c(0,1),ylim=c(0,10),main="",xlab="p",col=2,type="l",lwd=3)
par(new=TRUE)
plot(xVal,dbeta(xVal,a+y,b+n-y),xlim=c(0,1),ylim=c(0,10),main="1,5",xlab="p",type="l",lwd=3)

a<-5
b<-1

plot(xVal,dbeta(xVal,a,b),xlim=c(0,1),ylim=c(0,10),main="",xlab="p",col=2,type="l",lwd=3)
par(new=TRUE)
plot(xVal,dbeta(xVal,a+y,b+n-y),xlim=c(0,1),ylim=c(0,10),main="5,1",xlab="p",type="l",lwd=3)