``` r
library(forcedMigration)
library(microbenchmark)

microbenchmark({res1 = test_me(delta = 0.2, tau_ell = 0.6, tau_n = 0.6,
                               r = 0.1, a0 = 5.0, a1 = 0, p = 0.68,
                               q = 3859, H_bar = 331293, omega_n = 0.67,
                               gamma = 0.25, alpha = 0.005)}, times = 10, unit = "s")
#> Unit: seconds
#>                                                                                                                                                                                         expr
#>  {     res1 = test_me(delta = 0.2, tau_ell = 0.6, tau_n = 0.6, r = 0.1,          a0 = 5, a1 = 0, p = 0.68, q = 3859, H_bar = 331293, omega_n = 0.67,          gamma = 0.25, alpha = 0.005) }
#>          min        lq      mean     median         uq        max neval
#>   0.08146182 0.0816337 0.0825702 0.08219371 0.08320416 0.08527273    10


microbenchmark({res2 = test_me(delta = 0.2, tau_ell = 0.6, tau_n = 0.6,
                               r = 0.1, a0 = 5.0, a1 = 0, p = 0.68,
                               q = 3859, H_bar = 331293, omega_n = 0.67,
                               gamma = 0.25, alpha = 0.005)}, times = 10, unit = "s")
#> Unit: seconds
#>                                                                                                                                                                                         expr
#>  {     res2 = test_me(delta = 0.2, tau_ell = 0.6, tau_n = 0.6, r = 0.1,          a0 = 5, a1 = 0, p = 0.68, q = 3859, H_bar = 331293, omega_n = 0.67,          gamma = 0.25, alpha = 0.005) }
#>         min         lq       mean     median         uq        max neval
#>  0.04567314 0.04602822 0.04679136 0.04647838 0.04759033 0.04894398    10
```