Table 2 -

Approximate average run times for a single simulation for the completely deterministic ODE system (Euler method with step-sizeh = 10−5), the hybrid method with different threshold values (), and the completely stochastic Gillespie algorithm (rows). Each system is run for the single mutation, double mutation, and triple mutation models (columns). In each system we assume all strains are neutral (Fi = 1 for all i). The other parameters are N = 3, μ = 3 × 10−5, λ = 1.59 × 107, β = 3.60 × 10−9, γ = 0, d = 0.016, and a = 0.45.