# EXAMPLE 5: Introduction to Statistics for Astrophysicists
# SS 2012
# JOCHEN WELLER

# BINOMIAL PROBABILITY

N=20
r=seq(0,N,by=1)

p=0.1

plot(r,dbinom(r,N,p),type='h',xlim=c(0,20),xlab='r',ylab='f(r;N,p)',main=paste("N=",N,"; p=",p))