Isolate single species when plotting ensemble problem

Specific Domains Visualization
4

Here’s an example if you are using Catalyst:

using Catalyst, DifferentialEquations, Plots
rn = @reaction_network begin
       1.0, A --> 0
       2.0, B --> 0
       4.0, C --> 0
       end
dprob = DiscreteProblem(rn, [:A => 100, :B => 10, :C => 5], (0.0,20.0))
jprob = JumpProblem(rn, dprob, Direct())
eprob = EnsembleProblem(jprob)
sims = solve(eprob, SSAStepper(), trajectories=100)
@unpack C = rn
plot(sims, idxs=C)