fish <- function(n) choose(50,3)*choose(n-50,47)/choose(n,50) n <- 50:2000 plot(n, fish(n), type="l")