b1  <- seq(0, 1, by = 0.02)    
by1 <- dbeta(b1, shape1 = 5, shape2 = 20) 
plot(by1)