dzeta(1:20, p = 2)
## Not run: 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") ## End(Not run)
Run the code above in your browser using DataLab