qsize <- 1; a <- 0.5; x <- qsize:(qsize+10)
plot(x, dbort(x, qsize, a), type = "h", las = 1, col = "blue",
ylab = paste("fbort(qsize=", qsize, ", a=", a, ")"),
main = "Borel-Tanner density function")
Run the code above in your browser using DataLab