Collatz in R

gistfile1.txt

#Collatz Conjecture
#Will produce graph and number of steps to converge to 1


collatz <- function(start_val, show_me) {
  if(missing(show_me)){
    show_me <- start_val
  }
  n <- start_val
  comp <- c()
  j <- 0
  for(i in 1:n){
    steps <- 0
    j <- j + 1
    while(i!=1) {
      #print(n)
      if(i%%2 == 0){
        i <- i/2
      } else i <- (3*i) + 1
      steps <- steps + 1
    }
    comp[j] <- steps
  }
  numbers <- c(1:n)
  plot(comp~numbers, xlab="Starting Number", ylab="Steps", pch=19, col="blue", cex=0.2)
  paste(show_me, "took", comp[show_me], "steps to converge to 1")
}