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