dzeta(1:20, p = 2)
plot(1:6, dzeta(1:6, p = 4), type = "h", las = 1, ylab = "Probability",
main = "zeta probability function; black: p = 4; blue: p = 2")
points(0.10 + 1:6, dzeta(1:6, p = 2), type = "h", col = "blue")
Run the code above in your browser using DataLab