Figure 4

A spatial predator prey system with a safe haven. First moments (top row) and second moments (bottom row) of the predator–prey system described in (9). Voxels (1,1), (3,3) and (5,5) are shown in subplots a, b and c, respectively. Parameters used are k ₁  = 1[A]⁻¹  s ⁻¹, k ₂  = 0.02[A]⁻¹[B]⁻¹  s ⁻¹  k ₃  = 0.25[A]⁻¹  s ⁻¹, k ₄  = 1[B]⁻¹  s ⁻¹ k _diff   = 0.5[A]⁻¹  s ⁻¹ for A and k _diff   = 0.5[B]⁻¹  s ⁻¹. Initial conditions are 100 individuals of A and B in every voxel. Time steps for the SCLE and deterministic approaches were 0.001 seconds. RDME results are shown as blue lines. SCLE and deterministic results are shown as red and black dots, respectively.