#!/usr/bin/Rscript n <- 30 bacteria <- n viruses <- 1 steps <- 0 # b and v are only for plots: b <- bacteria v <- viruses while(bacteria>0) { bacteria <- bacteria-viruses viruses <- 2*viruses bacteria <- 2*bacteria steps <- steps+1 # only for plots: b <- c(b, bacteria) v <- c(v, viruses) } cat("n =", n, "\nsteps =", steps, "\n") # plot: plot(b,col=c("orange", "green"),xlab="steps",ylab="count") points(v,col=c("orange", "green")) lines(b,col="blue") lines(v,col="red") legend(1,max(b), c("viruses","bacteria"), fill=c("red", "blue"))