Quantcast
Channel: First steps - JuliaLang
Viewing all articles
Browse latest Browse all 2795

Histogram from a set of values from EnsembleProblem

$
0
0

@Frederic_Gilt wrote:

Dear All,
I have a system of six SDEs I solved with

prob_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, then v=hcat(sol.u[6,:]) which failed too…

sorry again for this naive question
thank you !
Frederic

Posts: 2

Participants: 2

Read full topic


Viewing all articles
Browse latest Browse all 2795