params <- function(n = 2,
                   K=20,
                   k0_AI=3.5,
                   k1_AI=0.033,
                   eta=0.7,
                   k_deg=0.525,
                   k_diff=2.1,
                   alpha_tetR=24,
                   alpha_cI=24,
                   alpha_lacI=24,
                   beta=1.05){
  paramlist <- c(n,
            K,
            k0_AI,
            k1_AI,
            eta,
            k_deg,
            k_diff,
            alpha_tetR,
            alpha_cI,
            alpha_lacI,
            beta)
  return(paramlist)

}