Plotting only certain trajectories of ensemble solution

Specific Domains Modelling & Simulations
, ,
1

If I soleve an ensemble problem:

using DifferentialEquations
using Plots

prob = ODEProblem((u,p,t)->1.01u,0.5,(0.0,1.0))

function prob_func(prob,i,repeat)
  remake(prob,u0=rand()*prob.u0)
end

ensemble_prob = EnsembleProblem(prob,prob_func=prob_func)
sim = solve(ensemble_prob,Tsit5(),EnsembleDistributed(),trajectories=10)

plot(sim,linealpha=0.4)

I will plot 10 trajectories. How do I do if I only want to plot the first, or every second, trajectory? I belive it used to be idxs=1, but that is now used to designate which variable you want to plot.

1 Like
2

Just poking around:

plot(sim.u[1], linealpha=0.4)
plot(sim.u[1:2:end], linealpha=0.4)
4 Likes