1.0 1.0 1.0 1.0 0.0 9.0 (0.2 * (((((-5.0 + x) ^ 2.0) + (y ^ 2.0)) < 1.0) || ((((-5.0 + x) ^ 2.0) + ((-10.0 + y) ^ 2.0)) < 1.0))) ((S_N * I_N * b) - (g * I_N)) (g * I_N) - (S_N * I_N * b) J_Infpop I_N_init J_Recovpop diffusionRate R_N_init J_Suspop diffusionRate S_N_init 0.0, 1.0, 2.0, 3.0