@Frederic_Gilt wrote:
Dear All,
I have a system of six SDEs I solved withprob_sde = SDEProblem(qdelta,qsigma,u0,(0.0,600.0),noise=W) ensembleprob = EnsembleProblem(prob_sde,prob_func=prob_func) sol = solve(ensembleprob,SRIW1(),EnsembleDistributed(),trajectories=50) plotly() plot(sol,vars=(0,6),linealpha=0.4)
which seems to work well, unfortunately I need the histogram of u[6], but I merely don’t know how to extract from ‘sol’ all the values from the 50 trajectories so that I can simply call the histogram ploting tool… I tried
convert(Array,sol)
which failed, thenv=hcat(sol.u[6,:])
which failed too…sorry again for this naive question
thank you !
Frederic
Posts: 2
Participants: 2