full.binomial <- function(n,p){
x <- matrix(0:n,n+1,1)
y <- dbinom(x,n,p) 
w <- pbinom(x,n,p)
result<-cbind(x,y,w,1-w)
rownames(result) <- rep("",n+1)
colnames(result)<-c('x','Pr(X = x)','Pr(X <= x)','Pr(X > x)')
return(result)
}