N = 10; s=0.5; y = 1:N
proby = dzipf(y, N=N, s=s)
plot(y, proby, type="h", col="blue", ylab="P[Y=y]", ylim=c(0,0.2),
main=paste("Zipf(N=",N,", s=",s,")", sep=""))
sum(proby) # Should be 1
max(abs(cumsum(proby) - pzipf(y, N=N, s=s))) # Should be 0
Run the code above in your browser using DataLab